# Production Asymmetry at the LHC from Heavy-Quark Recombination

###### Abstract

The asymmetry in the forward region production cross section of is calculated using the heavy-quark recombination mechanism for collisions at TeV. By suitable choices of four non-perturbative parameters, our calculated results can reproduce those obtained at LHCb. We find when integrated over and , which agrees with as measured by LHCb. Furthermore, the calculated distributions in and agree reasonably well with those obtained at LHCb.

Observation and proper interpretation of CP-violation in charm system could provide an outstanding opportunity for indirect searches for physics beyond the standard model. Even already available bounds on CP-violating interactions provide rather stringent constraints on the models of new physics because of availability of large statistical samples of charm data from LHCb, Belle, and BaBar experiments. Larger samples will be available soon from both and machines Butler:2013kdw ().

One of the simplest signals for CP-violation in charm is obtained by comparing partial decay widths of charm mesons to those of anti-charm mesons. While CPT-symmetry requires the total widths of and to be the same, the partial decay widths and are different in the presence of CP-violation, which is signaled by a non-zero value of the asymmetry

(1) |

This signal is reasonably robust for mesons, provided that the number of decaying particles and anti-particles is the same. However, at the Large Hadron Collider (LHC), the number of produced and mesons might not be the same due to the fact that the initial state contains two protons. With CP-violating asymmetries expected to be at the per mille levels Artuso:2008vf (), it is important to examine production asymmetry of -mesons both experimentally and theoretically.

Indeed, fixed-target experiments have already observed large asymmetries of charmed mesons and baryons in the forward region. In hadroproduction, the charmed hadrons are preferentially produced with a light valence quark of the same type as what appears in the hadronic beam, for example Aitala:1996hf (). This has been termed the leading particle effect. More recently, a similar asymmetry in production, defined as

(2) |

has been measured in the forward region to be by the LHCb collaboration LHCb:2012fb (). What are the theoretical expectations for this asymmetry?

Factorization theorems of perturbative QCD Collins:1985gm () state that heavy hadron production cross section can be written in a factorized form. At the LHC, the cross section for producing a () meson in a collision, at leading order in a expansion, is given by

(3) |

where is the partion distribution function for parton in the proton, is the partonic cross section and is the fragmentation function describing hadronization of a quark into a meson. The corresponding equation for is obtained by replacing by . Charge conjugation is expected to be a good symmetry in QCD, so . Thus, perturbative QCD predicts that , which is at least true at leading order in expansion.

This conclusion led theorists to examine other mechanisms for generating production asymmetry of Eq. (2), including attempts to describe the effect phenomenologically. The main idea of those approaches is to identify phenomenological mechanisms that can lead to enhanced production asymmetries, such as “meson cloud” effects. The results of these model-dependent calculations can be found in Refs. Norrbin:2000zc (); Cazaroto:2013wy (). We note that it might be challenging to interpret some of those mechanisms in QCD.

To reconcile the experimental observations with QCD, we note that there are corrections to Eq. (3) that scale as powers of
and . In principle, one can expect non-vanishing power-suppressed contributions to
at low . A QCD-based model for these power corrections is the heavy-quark recombination
mechanism Braaten:2001bf (); Braaten:2001uu (); Braaten:2002yt (); Braaten:2003vy (). In this scenario, a light quark involved in the hard
scattering process combines with the heavy quark produced in that interaction to form the final state meson, leading to corrections of order
. In what follows, after a quick review of the heavy-quark recombination mechanism,^{4}^{4}4For a full
review, please see Refs. Braaten:2001bf (); Braaten:2001uu (); Braaten:2002yt (); Braaten:2003vy (). we calculate
due to heavy quark recombination.

Imagine production of a heavy meson with the light quark of the same flavor as that appears in the beam. For instance, for a proton beam we could have or states, which we shall generically call . The recombination process, shown in Fig. 1 (a), comes in as a power-suppressed correction to Eq. (3). As mentioned, the light quark in the production of comes from the incident proton. The contribution to the cross section is given by

(4) |

where indicates that the light quark of flavor with momentum in the rest frame is produced in the states , where labels the color and angular momentum quantum numbers of the quark pair. The cross section is factored into a perturbatively calculable piece and the nonperturbative factor encoding the probability for the quark pair with quantum number to hadronize into the final state including the . The perturbative piece was calculated to lowest order in Braaten:2001bf (). Equation (4) must then be convoluted with the proton parton distribution functions to get the final hadronic cross section.

Besides the process, there are also contributions from , as shown in Fig. 1 (b). Using the method introduced in Braaten:2001bf (), the partonic cross sections from initial state charm are calculated to be

(5) | |||||

where we have defined and .

The quark in Eq. (4) could fragment into a meson, this time of opposite flavor, i.e., a or , generically labeled . Thus, to get the full rate due to recombination for producing mesons, we also need to account for the contribution where a light antiquark comes from the proton, while the fragments into a . We thus have three contributions,

(6a) | ||||

(6b) | ||||

(6c) |

where can be any hadron. The recombination cross section for producing a is obtained by taking the charge conjugate of the above equations. Below, we will neglect -violation and take . For simplicity, in process c) we will restrict to be only and sum over and with flavor symmetry assumed.

As discussed in Braaten:2002yt (), the nonperturbative parameters with the same flavor and angular momentum quantum numbers as the scale as . However, the amplitudes for production with are suppressed relative to the S-wave states. On the other hand, transition is achieved via emission of magnetic-type gluons, which, contrary to the heavy quarkonia case, is not suppressed for -mesons. Thus, the leading contributions to productions of mesons by heavy-quark recombination consists of four possible options of :

(7) |

These nonperturbative parameters must be extracted from data. Neglecting and , the combination was determined to be by fitting to the E687 and E691 fixed-target photoproduction data Braaten:2001uu (). Neglecting , and , the parameter was determined to be by fitting to data from the E791 experiment Braaten:2002yt (). In this paper, we take and . It turns out that these two contributions only account for of the measured asymmetry at LHCb in Ref. LHCb:2012fb (). Therefore, we include and and choose values of similar size as the spin-matched parameters. We use MSTW 2008 LO PDFs with GeV and the Peterson parametrization for the fragmentation function Peterson:1982ak () is used for :

(8) |

is chosen to be . are determined by the averages of the measured fragmentation probabilities listed in Abramowicz:2013eja (). For the perturbative QCD rate, Eq. (3), which has no asymmetry if we ignore violation but enters into the denominator of Eq. (2), we use the LO cross section and include feed down from . The factorization scale is set to be .

When integrated over and , excluding the region with , , the asymmetry for is found to be with , , and . Figure 2 shows as a function of pseudorapidity and transverse momentum of the mesons as predicted by the heavy-quark recombination mechanism. Data from Ref. LHCb:2012fb () are shown as well. The calculated distributions are reasonably consistent with the data.

In summary, we have calculated the asymmetry using the heavy-quark recombination mechanism for production at the LHCb experiment. The measured asymmetry of in the kinematic range and LHCb:2012fb () can be reproduced using reasonably sized non-perturbative parameters . Further, the and distributions are simultaneously reproduced by the heavy-quark recombination mechanism.

## Acknowledgements

We thank Tao Han and Thomas Mehen for useful discussions. AKL and WKL are supported in part by the National Science Foundation under Grant No. PHY-1212635. AAP is supported in part by the U.S. Department of Energy under contract DE-FG02-12ER41825.

## References

- (1) J. N. Butler et al. [Quark Flavor Physics Working Group Collaboration], “Working Group Report: Quark Flavor Physics,” arXiv:1311.1076 [hep-ex];
- (2) M. Artuso, B. Meadows and A. A. Petrov, Ann. Rev. Nucl. Part. Sci. 58, 249 (2008) [arXiv:0802.2934 [hep-ph]].
- (3) E. M. Aitala et al. [E791 Collaboration], Phys. Lett. B 371, 157 (1996) ; ibid. 411, 230 (1997); G. A. Alves et al. [E769 Collaboration], Phys. Rev. Lett. 72, 812 (1994); ibid. 77, 2392 (1996); M. Adamovich et al. [BEATRICE Collaboration], Nucl. Phys. B 495, 3 (1997); M. Adamovich et al. [WA82 Collaboration], Phys. Lett. B 305, 402 (1993); F. G. Garcia et al. [SELEX Collaboration], Phys. Lett. B 528, 49 (2002); E. M. Aitala et al. [E791 Collaboration], Phys. Lett. B 495, 42 (2000); M. I. Adamovich et al. [WA89 Collaboration], Eur. Phys. J. C 8, 593 (1999).
- (4) R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 718, 902 (2013) [arXiv:1210.4112 [hep-ex]].
- (5) J. C. Collins, D. E. Soper and G. F. Sterman, Nucl. Phys. B 263, 37 (1986).
- (6) E. Norrbin and T. Sjostrand, Eur. Phys. J. C 17, 137 (2000) [hep-ph/0005110].
- (7) E. R. Cazaroto, V. P. Goncalves, F. S. Navarra and M. Nielsen, Phys. Lett. B 724, 108 (2013) [arXiv:1302.0035 [hep-ph]].
- (8) E. Braaten, Y. Jia and T. Mehen, Phys. Rev. D 66, 034003 (2002) [hep-ph/0108201].
- (9) E. Braaten, Y. Jia and T. Mehen, Phys. Rev. D 66, 014003 (2002) [hep-ph/0111296].
- (10) E. Braaten, Y. Jia and T. Mehen, Phys. Rev. Lett. 89, 122002 (2002) [hep-ph/0205149].
- (11) E. Braaten, M. Kusunoki, Y. Jia and T. Mehen, Phys. Rev. D 70, 054021 (2004) [hep-ph/0304280].
- (12) C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas, Phys. Rev. D 27, 105 (1983).
- (13) H. Abramowicz et al. [ZEUS Collaboration], JHEP 1309, 058 (2013) [arXiv:1306.4862 [hep-ex]].