# Quarkonium Suppression from Coherent Energy Loss in Fixed-Target Experiments Using LHC Beams

Authors: François Arleo; Stéphane Peigné
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

Quarkonium production in proton-nucleus collisions is a powerful tool to disentangle cold nuclear matter effects. A model based on coherent energy loss is able to explain the available quarkonium suppression data in a broad range of rapidities, from fixed-target to collider energies, suggesting  coherent energy loss in cold nuclear matter to be the dominant effect in quarkonium suppression in p-A collisions. This could be further tested in a high-energy fixed-target experiment using a proton or nucleus beam. The nuclear modification factors ofJ/ψ and Υ as a function of rapidity are computed in p-A collisions at s=114.6 GeV, and in p-Pb and Pb-Pb collisions at s=72 GeV. These center-of-mass energies correspond to the collision on fixed-target nuclei of 7 TeV protons and 2.76 TeV (per nucleon) lead nuclei available at the LHC.


---

## Body

## 1. Introduction

Understanding the physical origin of quarkonium (J/ψ, Υ) suppression in proton-nucleus (p-A) collisions has been a challenge for the past thirty years. This would of course be a prerequisite in order to perform reliablebaseline predictions in heavy-ion collisions, where quarkonia are expected to be dissociated due to Debye screening of the heavy-quark potential at finite temperature [1]. Perhaps even more importantly, the wealth of data (especially for J/ψ and Υ) available in p-A collisions could help to understand generic features of hard QCD processes in a nuclear environment.In everyday language, we often make the distinction between “fixed-target” and “collider” experiments when it comes to quarkonium production. This separation might look a bit artificial but not entirely:(i)
In fixed-target experiments, luminosities are often high, leading to abundant yields and consequently reduced statistical uncertainties. Moreover, thanks to the boost of the center-of-mass frame of the collision, the rapidity coverage of such experiments can extend up to largernegative values of rapidity (we follow the usual convention where positive (negative) rapidities correspond to the proton (nucleus) fragmentation region) or Feynman-x, xF≃2M⊥/s×sinhy, where M⊥ is the quarkonium transverse mass. The highest energies ever reached are of course rather modest, s=38.7 GeV and s=41.6 GeV, using, respectively, 800 GeV and 920 GeV proton beams at the Tevatron and at HERA. This, however, allows for probing also largerpositive values of xF than in collider experiments.(ii)
At collider energies—RHIC and LHC, to quote only the facilities accelerating heavy ions—unprecedented energies can be reached, respectively,s=200 GeV and s=5.02 TeV, for instance, making the production of Υ states, marginally measured in fixed-target experiments, easier. In terms of acceptance, quarkonia are detected in a narrow window in xF, centered around xF=0.Let us illustrate this with a few examples, starting with one of the first experiments which measuredJ/ψ suppression in p-A collisions. The NA3 spectrometer at the CERN SPS collected 1.5 million J/ψ events (!) in hadron-nucleus collisions [2], allowing for precise measurements close to the kinematic edge of phase-space, xF≲0.75 (on the contrary, the coverage at RHIC and LHC is, resp., xF≲0.2 and xF≲0.02 for J/ψ production). It is remarkable that these data, taken in the early 1980s, prove to be as competitive as the most recent LHC results when it comes to understanding J/ψ suppression in nuclei. More than a decade after NA3, the FNAL E866 experiment reported on high-statistics measurements of J/ψ (3×106 events) and ψ′ (105) production on several nuclear targets in the range -0.2≲xF≲0.9 [3]. Other important measurements were reported by experiments at the SPS (NA38 [4], NA50 [5], and NA60 [6]) and HERA (HERA-B [7]), yet on a more restricted xF range.These data are nicely supplemented by those carried out in d-Au collisions at RHIC (PHENIX [8, 9] and STAR [10]) and in p-Pb collisions at LHC (ALICE [11] and LHCb [12]). At LHC, the relative suppression of Υ excited states (2S, 3S) with respect to 1S states has been performed by CMS [13], not to mention open heavy-flavor data (D mesons in ALICE [14], B mesons in CMS [15], and nonprompt J/ψ coming from B decays in ALICE [16] and LHCb [12]).Several cold nuclear matter (CNM) effects could in principle affect quarkonium yields in proton-nucleus collisions. Without being comprehensive, let us mention the following ones:(i)
Quarkonia may interact inelastically with the surrounding nucleons they may encounter while propagating through the nucleus. Suchnuclear absorption may happen when the quarkonium formation time (in the rest frame of the nucleus) is comparable or less than the medium length L, τf×coshylab≲L, where τf is theproper formation time (τf≃0.3 fm for both J/ψ and Υ) and ylab is the quarkonium rapidity in the nucleus rest frame. Note that ylab is directly related to the momentum fraction x2 carried by the target parton, coshylab=M⊥/(2mpx2), where mp is the proton mass.(ii)
Parton distribution functions (PDF) are known to be different in a proton and in a nucleus at all values ofx. Such nuclear PDF (nPDF) effects could either suppress or enhance quarkonium yields in p-A collisions (with respect to p-p collisions) depending on the value of the momentum fraction x2. When x2 is small (typically when the time for the hard process to occur is large in the nucleus rest frame, th≃(1/M⊥)×coshylab=1/(2mpx2)≳L, using L=10 fm, this would correspond to x2≲10-2), the nucleons in the nucleus act coherently leading to a reduction of the quarkonium yield in a nucleus—calledshadowing [17], orsaturation [18] to use a more modern language—as compared to the incoherent sum over A independent nucleons.(iii)
Nuclear transverse momentum broadening of the heavy-quark pair induces coherent gluon radiation, arising from the interference between emission amplitudes off the initial projectile parton and the final color octet quark pair. This coherent medium-induced radiation leads to an average induced energy loss proportional to the quarkonium energy [19]. The consequences of coherent energy loss are quarkonium suppression (resp., enhancement) at large positive (resp., large negative) values of the rapidity and at all center-of-mass energies of the p-A collision.Although each of these CNM effects is plausible, it does not necessarily mean that all play a role in the nuclear dependence of quarkonium production; in particular, the strength of each CNM effect is usually unknowna priori. A sound strategy is to investigate each of these effects separately, through a systematic and quantitative comparison to all available data, while keeping the smallest number of assumptions and free parameters.Quarkonium suppression reported at forward rapidities cannot be reproduced by either nuclear absorption or nPDF effects, nor by a mixture of both. Although the comparison to RHIC and LHC data only may still give the impression that strong nPDF effects could explainJ/ψ data [11] (for examples of nPDF effects on quarkonium production in p-Pb collisions at LHC, see [20, 21]), the significant suppression measured by the fixed-target experiments (NA3 and E866) on a wider xF range is clearly incompatible with the predictions of these two effects. An elegant way to be persuaded is to plot J/ψ suppression data as a function of x2=M⊥/s×exp(-y)[22]. The suppression from either nuclear absorption or nPDF effects is expected to be a function of x2 only, independent of s. This is in strong disagreement with the accumulated data from fixed-target and RHIC experiments (see [3, 23]). Without a doubt, the world data indicate that at least another cold nuclear matter effect is at play.Contrary to nuclear absorption or nPDF effects, the sole effect of coherent energy loss is able to reproduce the data on quarkonium suppression, from fixed-target to collider energies [24–26]. Detailed comparisons were published elsewhere, so let us only highlight the phenomenological successes of this approach:(i)
ThexF (or y) dependence of J/ψ suppression is well reproduced on a very large domain (up to large values of xF≲0.8, when data are available) and at all center-of-mass energies, from s≃20 GeV to s=5 TeV.(ii)
Thep⊥ dependence is well-reproduced too, either at a fixed-target experiment (E866) or at colliders (RHIC, LHC), although the p⊥ dependence seems slightly more abrupt in the model than in collider data. The centrality dependence measured by PHENIX at RHIC is also nicely described.(iii)
Υ measurements in p-A collisions are compatible with the expected mass dependence of coherent energy loss, although the present experimental uncertainties are still fairly large.(iv)
Finally, an original prediction of coherent energy loss is a different magnitude of quarkonium suppression in p-A andπ-A collisions (in contrast with nuclear absorption effects, which should be independent of the projectile hadron), in agreement with the measurements of NA3. The strength of coherent energy loss depends on a single free parameter, namely, the magnitude of the cold nuclear matter transport coefficient, q^0=0.075 GeV2/fm at x=10-2, obtained from a fit of the precise E866 measurements in p-W collisions.We find it appealing that the variety of quarkonium measurements in p-A collisions can be described using asingle CNM effect. Of course, by no means does this imply that no other CNM effects could play a role too, yet these clearly appear to be subleading when the quarkonium is produced at “large enough” rapidity. Both nPDF and coherent energy loss effects could in principle be incorporated consistently in the picture. As a matter of fact, attempts have been made in [25]. However, given the large theoretical uncertainties on nuclear parton distributions, due to the lack of small-x measurements in nuclear collisions (in this respect a high-energy electron-ion collider would be highly beneficial for the community. Let us mention in passing that no coherent energy loss effects are expected in deep-inelastic scattering experiments as the incoming projectile particle is color neutral [19]), especially for gluon densities at small x, we prefer to focus on the single (but in our opinion, leading) effect of coherent energy loss for which rather precise calculation can be performed.An exciting possibility to further constrain cold nuclear matter effects (on quarkonium production, but not only) would be to smash the LHC proton and lead beams on a collection of fixed nuclear targets [27]. We believe that this proposal would combine the above-discussed advantages of fixed-targetand collider experiments.In this paper, the predictions for quarkonium suppression due to coherent energy loss in p-A collisions ats=114.6 GeV (corresponding to the nominal 7 TeV proton beam energy at the LHC) and Pb-A collisions at s=72 GeV (corresponding to the 2.76 TeV lead beam) are given. Before this, we recall in the next section the main ingredients of our approach.


## 2. Coherent Energy Loss Model in a Nutshell

### 2.1. Formulation

We briefly detail in this section the basics of the model based on coherent energy loss used to describeψ (denoting J/ψ or Υ) suppression measured in proton-nucleus collisions. (The model can also be formulated in heavy-ion (A-B) collision; see [28] for details.) The single differential p-A production cross section as a function of the ψ energy reads [25] (1)1AdσpAψdEE=∫0ɛmaxdɛPɛ,E,lA2dσppψdEE+ɛ,where E (resp., ɛ) is the energy (resp., energy loss) of the QQ¯ pair in the rest frame of the nucleus A. The upper limit on the energy loss is ɛmax=minE,Ep-E, where Ep is the beam energy in that frame. P denotes the energy loss probability distribution, orquenching weight.The quenching weight is related to the medium-induced, coherent radiation spectrumdI/dɛ given in [25] (and earlier in [19]), which is a very good approximation to the exact spectrum computed to all orders in the opacity expansion [29]. It depends on the accumulated transverse momentum transfer lA=q^L (assumed to satisfy lA≪M⊥) due to soft rescatterings in the nucleus, where L is the medium path-length and q^ is the transport coefficient in cold nuclear matter (in our picture, the relevant transverse broadening q^L is that of the compact color octet crossing the nucleus, and L is thus the path-length from the “front” surface to the “back” surface of the nucleus for a given impact parameter. Note that, in nuclear absorption models, the relevant length is that from the hard production vertex to the “back” surface of the nucleus, that is, a factor ~2 smaller than the medium length L used here and quoted in Table 1). More precisely [25], (2)q^≡q^010-2min⁡x0,x20.3;x0≡12mpL;x2≡M⊥se-y,where y is the quarkonium rapidity in the center-of-mass frame of the proton-nucleon collision.Table 1
Values ofL used in p, Ca, Cu, and Pb targets.
Nucleus
p
Ca
Cu
Pb
Atomic mass
1
40
63
208
L (fm)
1.5
5.69
6.67
10.11Using the fact that the quenching weight is a scaling function of the variableɛ/E, namely, EP(ɛ,E,l2)=P^(ɛ/E,l2), we can rewrite (1) as [28] (3)1AdσpAψdyy,s=∫0δymaxydδyP^eδy-1,q^yL·dσppψdyy+δy,s.Here, δymax(y)=minln⁡2,ymax-y, with ymax=ln⁡(s/M⊥) the maximal ψ rapidity (in the proton-nucleon c.m. frame) allowed by kinematics. Using (3), we can determine the nuclear modification factor in p-A collision: (4)RpAψy,s=1/AdσpAψ/dyy,sdσppψ/dyy,s.As mentioned in the introduction, quarkonium may suffer inelastic interaction with the surrounding nucleons while escaping the nucleus. Although we do not aim to include such an effect in the present predictions, we nevertheless indicate the critical rapidityycrit: (5)ycrits,L≡ln⁡Lτf·2mps,below which nuclear absorption might come into play.


### 2.2. Ingredients

The medium lengthL is obtained from a Glauber model calculation using realistic nuclear densities. The values are given in [25] and reproduced in Table 1 for the nuclei of interest in the present paper. In addition, (3) requires the knowledge of the p-p cross section. It is given by a simple parameterization dσppψ/dy∝1-(2M⊥/s)coshyn(s), where the exponent n is obtained from a fit to p-p measurements. Lacking p-p data at the energies of interest (s=114.6 GeV in p-A and s=72 GeV in Pb-A collisions), an interpolation between the values obtained at FNAL (s=38.7 GeV) and RHIC (s=200 GeV) energies has been performed. The exponents used in the present paper are given in Table 2. Note that the normalization of the p-p cross section is irrelevant here as it cancels out when computing (4).Table 2
Values ofn used at s=72 GeV and s=114.6 GeV for J/ψ and ϒ.
Mode
Pb-A
p-A
s (GeV)
72
114.6
n
J
/
ψ
5.1
±
0.2
6.0
±
0.3
n
ϒ
4.1
±
0.3
5.0
±
0.4The transport coefficientq^0 is the only free parameter of the model. It is determined by fitting the J/ψ suppression measured by E866 [3] in p-W over p-Be collisions (s=38.7 GeV); see [25]. The obtained value is q^0=0.075-0.005+0.015 GeV2/fm.


## 2.1. Formulation

We briefly detail in this section the basics of the model based on coherent energy loss used to describeψ (denoting J/ψ or Υ) suppression measured in proton-nucleus collisions. (The model can also be formulated in heavy-ion (A-B) collision; see [28] for details.) The single differential p-A production cross section as a function of the ψ energy reads [25] (1)1AdσpAψdEE=∫0ɛmaxdɛPɛ,E,lA2dσppψdEE+ɛ,where E (resp., ɛ) is the energy (resp., energy loss) of the QQ¯ pair in the rest frame of the nucleus A. The upper limit on the energy loss is ɛmax=minE,Ep-E, where Ep is the beam energy in that frame. P denotes the energy loss probability distribution, orquenching weight.The quenching weight is related to the medium-induced, coherent radiation spectrumdI/dɛ given in [25] (and earlier in [19]), which is a very good approximation to the exact spectrum computed to all orders in the opacity expansion [29]. It depends on the accumulated transverse momentum transfer lA=q^L (assumed to satisfy lA≪M⊥) due to soft rescatterings in the nucleus, where L is the medium path-length and q^ is the transport coefficient in cold nuclear matter (in our picture, the relevant transverse broadening q^L is that of the compact color octet crossing the nucleus, and L is thus the path-length from the “front” surface to the “back” surface of the nucleus for a given impact parameter. Note that, in nuclear absorption models, the relevant length is that from the hard production vertex to the “back” surface of the nucleus, that is, a factor ~2 smaller than the medium length L used here and quoted in Table 1). More precisely [25], (2)q^≡q^010-2min⁡x0,x20.3;x0≡12mpL;x2≡M⊥se-y,where y is the quarkonium rapidity in the center-of-mass frame of the proton-nucleon collision.Table 1
Values ofL used in p, Ca, Cu, and Pb targets.
Nucleus
p
Ca
Cu
Pb
Atomic mass
1
40
63
208
L (fm)
1.5
5.69
6.67
10.11Using the fact that the quenching weight is a scaling function of the variableɛ/E, namely, EP(ɛ,E,l2)=P^(ɛ/E,l2), we can rewrite (1) as [28] (3)1AdσpAψdyy,s=∫0δymaxydδyP^eδy-1,q^yL·dσppψdyy+δy,s.Here, δymax(y)=minln⁡2,ymax-y, with ymax=ln⁡(s/M⊥) the maximal ψ rapidity (in the proton-nucleon c.m. frame) allowed by kinematics. Using (3), we can determine the nuclear modification factor in p-A collision: (4)RpAψy,s=1/AdσpAψ/dyy,sdσppψ/dyy,s.As mentioned in the introduction, quarkonium may suffer inelastic interaction with the surrounding nucleons while escaping the nucleus. Although we do not aim to include such an effect in the present predictions, we nevertheless indicate the critical rapidityycrit: (5)ycrits,L≡ln⁡Lτf·2mps,below which nuclear absorption might come into play.


## 2.2. Ingredients

The medium lengthL is obtained from a Glauber model calculation using realistic nuclear densities. The values are given in [25] and reproduced in Table 1 for the nuclei of interest in the present paper. In addition, (3) requires the knowledge of the p-p cross section. It is given by a simple parameterization dσppψ/dy∝1-(2M⊥/s)coshyn(s), where the exponent n is obtained from a fit to p-p measurements. Lacking p-p data at the energies of interest (s=114.6 GeV in p-A and s=72 GeV in Pb-A collisions), an interpolation between the values obtained at FNAL (s=38.7 GeV) and RHIC (s=200 GeV) energies has been performed. The exponents used in the present paper are given in Table 2. Note that the normalization of the p-p cross section is irrelevant here as it cancels out when computing (4).Table 2
Values ofn used at s=72 GeV and s=114.6 GeV for J/ψ and ϒ.
Mode
Pb-A
p-A
s (GeV)
72
114.6
n
J
/
ψ
5.1
±
0.2
6.0
±
0.3
n
ϒ
4.1
±
0.3
5.0
±
0.4The transport coefficientq^0 is the only free parameter of the model. It is determined by fitting the J/ψ suppression measured by E866 [3] in p-W over p-Be collisions (s=38.7 GeV); see [25]. The obtained value is q^0=0.075-0.005+0.015 GeV2/fm.


## 3. Results

### 3.1. p-A Mode

The predictions forJ/ψ and Υ suppression in p-Ca, p-Cu, and p-Pb collisions at s=114.6 GeV are shown in Figure 1. The rapidity range is chosen to match the acceptance of detectors like LHCb [30]. In terms of Feynman-x, the range -3<y<1 (resp., -2<y<1) corresponds to -0.54<xF<0.06 (resp., -0.60<xF<0.19) for J/ψ (resp., Υ (because of the more restricted phase-space from its larger mass, the rapidity acceptance is smaller for Υ than for J/ψ)).Figure 1
J
/
ψ (a) and Υ (b) suppression in p-Ca, p-Cu, and p-Pb collisions at s=114.6 GeV.
(a)
(b)TheJ/ψ suppression is rather moderate, less than 20%, and does not vary too strongly with rapidity except at very negative rapidity values, y<y0≃-2, where J/ψ enhancement (RpA>1) can be seen. In this rapidity region, however, nuclear absorption may come into play as can be seen from the vertical arrows indicating the values of ycrit (ycrit≃-1) for each target. The suppression depends on the atomic mass number of the target (Ca, Cu, and Pb), yet this nuclear dependence is tempered somehow by the slow dependence of the coherent energy loss scale, ω^, on the medium length, ω^∝L [25].The shape ofΥ suppression is similar. The value of the rapidity at which RpA(y0)=1 is y0≃-1, that is, one more unit than in the J/ψ channel. This can be understood from the approximate xF scaling present in the model [25] (at a given rapidity y, the corresponding value of xF is larger for Υ than for J/ψ due to the larger transverse mass, xF∝M⊥) which would predict the difference between these two “crossing points” to be y0Υ-y0J/ψ~ln⁡(MΥ/MJ/ψ)≃1.1. Once more, nuclear absorption may affect Υ suppression, though maybe not as much as the J/ψ because of its smaller radius.


### 3.2. A-p Mode

Let us move now to calculations corresponding to an incoming 2.76 TeV Pb beam on a proton and a Pb target, shown in Figure2. This configuration allows for probing more easily quarkonium suppression in the proton fragmentation, that is, at positive rapidities. (Although the Pb nucleus collides with a proton, we will keep the convention that positive values of y correspond to the proton fragmentation region.) The chosen rapidity range -2<y<2 (resp., -1.5<y<1.5) corresponds to -0.31<xF<0.31 (-0.56<xF<0.56) for J/ψ (Υ). The lower center-of-mass energy however shifts ycrit in p-Pb collisions towards larger values, possibly leading to more pronounced nuclear absorption.Figure 2
J
/
ψ (a) and Υ (b) suppression in p-Pb and Pb-Pb collisions at s=72 GeV.
(a)
(b)In Pb-Pb collisions, the suppression is naturally an even function ofy. In such collisions, one expects a hot medium to be formed leading to extra quarkonium suppression. Therefore, the results should rather be seen as baseline calculations than genuine predictions. Moreover, in A-A collisions, the condition for hadronization taking place outsideboth nuclei reads ycrit<y<-ycrit. This condition is only met in Pb-Pb collisions at s=72 GeV around midrapidity, |y|≲0.1. At larger |y|, the quarkonium state will be fully formed in one of the two nuclei and thus possibly sensitive to nuclear absorption.


## 3.1. p-A Mode

The predictions forJ/ψ and Υ suppression in p-Ca, p-Cu, and p-Pb collisions at s=114.6 GeV are shown in Figure 1. The rapidity range is chosen to match the acceptance of detectors like LHCb [30]. In terms of Feynman-x, the range -3<y<1 (resp., -2<y<1) corresponds to -0.54<xF<0.06 (resp., -0.60<xF<0.19) for J/ψ (resp., Υ (because of the more restricted phase-space from its larger mass, the rapidity acceptance is smaller for Υ than for J/ψ)).Figure 1
J
/
ψ (a) and Υ (b) suppression in p-Ca, p-Cu, and p-Pb collisions at s=114.6 GeV.
(a)
(b)TheJ/ψ suppression is rather moderate, less than 20%, and does not vary too strongly with rapidity except at very negative rapidity values, y<y0≃-2, where J/ψ enhancement (RpA>1) can be seen. In this rapidity region, however, nuclear absorption may come into play as can be seen from the vertical arrows indicating the values of ycrit (ycrit≃-1) for each target. The suppression depends on the atomic mass number of the target (Ca, Cu, and Pb), yet this nuclear dependence is tempered somehow by the slow dependence of the coherent energy loss scale, ω^, on the medium length, ω^∝L [25].The shape ofΥ suppression is similar. The value of the rapidity at which RpA(y0)=1 is y0≃-1, that is, one more unit than in the J/ψ channel. This can be understood from the approximate xF scaling present in the model [25] (at a given rapidity y, the corresponding value of xF is larger for Υ than for J/ψ due to the larger transverse mass, xF∝M⊥) which would predict the difference between these two “crossing points” to be y0Υ-y0J/ψ~ln⁡(MΥ/MJ/ψ)≃1.1. Once more, nuclear absorption may affect Υ suppression, though maybe not as much as the J/ψ because of its smaller radius.


## 3.2. A-p Mode

Let us move now to calculations corresponding to an incoming 2.76 TeV Pb beam on a proton and a Pb target, shown in Figure2. This configuration allows for probing more easily quarkonium suppression in the proton fragmentation, that is, at positive rapidities. (Although the Pb nucleus collides with a proton, we will keep the convention that positive values of y correspond to the proton fragmentation region.) The chosen rapidity range -2<y<2 (resp., -1.5<y<1.5) corresponds to -0.31<xF<0.31 (-0.56<xF<0.56) for J/ψ (Υ). The lower center-of-mass energy however shifts ycrit in p-Pb collisions towards larger values, possibly leading to more pronounced nuclear absorption.Figure 2
J
/
ψ (a) and Υ (b) suppression in p-Pb and Pb-Pb collisions at s=72 GeV.
(a)
(b)In Pb-Pb collisions, the suppression is naturally an even function ofy. In such collisions, one expects a hot medium to be formed leading to extra quarkonium suppression. Therefore, the results should rather be seen as baseline calculations than genuine predictions. Moreover, in A-A collisions, the condition for hadronization taking place outsideboth nuclei reads ycrit<y<-ycrit. This condition is only met in Pb-Pb collisions at s=72 GeV around midrapidity, |y|≲0.1. At larger |y|, the quarkonium state will be fully formed in one of the two nuclei and thus possibly sensitive to nuclear absorption.



---

*Source: 961951-2015-10-04.xml*