# Neutron-skin effect in direct-photon and charged-hadron production in Pb+Pb collisions at the LHC

###### Abstract

A well-established observation in nuclear physics is that in neutron-rich spherical nuclei the distribution of neutrons extends farther than the distribution of protons. In this work, we scrutinize the influence of this so called neutron-skin effect on the centrality dependence of high- direct-photon and charged-hadron production. We find that due to the estimated spatial dependence of the nuclear parton distribution functions, it will be demanding to unambiguously expose the neutron-skin effect with direct photons. However, when taking a ratio between the cross sections for negatively and positively charged high- hadrons, even centrality-dependent nuclear-PDF effects cancel, making this observable a better handle on the neutron skin. Up to 10 % effects can be expected for the most peripheral collisions in the measurable region.

###### Keywords:

Heavy-ion collisions Neutron-skin effect Centrality dependence Direct photons High- hadrons^{†}

^{†}journal: Eur. Phys. J. C\thankstext

e1email: ilkka.helenius@uni-tuebingen.de \thankstexte2email: hannu.t.paukkunen@jyu.fi \thankstexte3email: kari.eskola@jyu.fi

## 1 Introduction

In ultra-relativistic heavy-ion collisions the concept of centrality plays an important role in phenomena such as the jet energy loss Abelev:2013kqa (); Chatrchyan:2012nia (); Aad:2014bxa () or the systematics of azimuthal anisotropies Aad:2014eoa (); Chatrchyan:2012xq (). Experimentally, the centrality of a collision is usually defined according to the amount of energy seen in a specific part of the detector, typically at large pseudorapidities Aad:2014eoa (); Chatrchyan:2011sx (); Abelev:2013qoq (): the more energy observed, the more central the collision. The theoretical centrality categorizations are based on Glauber models Miller:2007ri (), in which the centrality is related to impact parameter (optical Glauber) or to the number of nucleon–nucleon collisions (Monte-Carlo Glauber). While there is no direct, unambiguous relation between the experimental and theoretical prescriptions, it is yet generally accepted that a correspondence exists in collisions of two heavy nuclei. In nucleon–nucleus collisions, however, the same experimental procedure has led to rather unexpected results Chatrchyan:2014hqa (); ATLAS:2014cpa (); Adare:2015gla () and it is now commonly believed that such a centrality classification induces a non-trivial bias on the hard process whose centrality dependence was to be measured Martinez-Garcia:2014ada (); Bzdak:2014rca (); Alvioli:2014eda (); Perepelitsa:2014yta (); Armesto:2015kwa ().

The Glauber models take the nuclear density distribution as an input and it is typically assumed to be identical for protons and neutrons. However, the measurements at lower energies indicate that the tail of the neutron density distribution extends farther than that of the proton density Tarbert:2013jze (); Tsang:2012se (); Zenihiro:2010zz (). While this so-called neutron-skin (NS) effect Horowitz:2000xj () should not have a great importance in the centrality classification itself, it leads to a growth of the relative number of neutrons at high impact parameters and thereby influences the observables sensitive to electroweak effects in peripheral (large impact parameter) collisions of two heavy nuclei. The impact of the NS effect to production in Pb+Pb and p+Pb collisions at the LHC was studied in Ref. Paukkunen:2015bwa ().

In this work, we extend the study of Ref. Paukkunen:2015bwa () to direct-photon and charged-hadron production at high transverse momenta () in Pb+Pb collisions at the LHC. The goal is to study whether the NS effect has a measurable impact on these observables and to quantify at which centralities and kinematics (transverse momentum, rapidity) the effect would be most pronounced. Our hope is that, later on, the NS effect could help to calibrate the centrality classification in collisions involving heavy ions.

## 2 Centrality-dependent hard-process cross section

Centrality classification is done here using the optical Glauber model as in Refs. Paukkunen:2015bwa (); Helenius:2012wd (). For the nuclear density distribution we use the two-parameter Fermi (2pF) distribution,

(1) |

where describes the radius of the nucleus and the thickness of the nuclear surface (skin) in nucleus with a mass number . To account for the NS effect the nuclear density is written as
where now the parameters of the 2pF distribution are different for protons and neutrons. Here we use the parameters from Ref. Tarbert:2013jze (), and , for neutrons and and for protons.^{1}^{1}1In the analysis of Ref. Tarbert:2013jze (), the proton density was taken as fixed when fitting the neutron parameters to the data. Therefore, the proton density has no uncertainty here.

The hard-process cross section in an collision for a given centrality class corresponding to an impact parameter interval (where ) can be calculated from

(2) |

where the nuclear thickness functions are obtained by integrating the density over the longitudinal (i.e. beam) direction,

(3) |

and are defined according to Fig. 20 of Ref. Helenius:2012wd (). The indices and run over combinations

. The impact-parameter intervals required in Eq. (2) correspond to the fractions of the total inelastic cross section , obtained as in Refs. Paukkunen:2015bwa (); Helenius:2012wd () by

(4) |

where

(5) |

We take Antchev:2013iaa (). The spatial dependence of the hard-process cross section arises here from the spatial dependence of nPDFs,

(6) |

where are perturbative coefficient functions and and are parton flavour indices. The nPDFs appearing in Eq. (6) above are defined as

(7) |

where is the free nucleon PDF (here CT10NLO Lai:2010vv ()) and the nuclear modification which depends on the transverse position of the nucleon inside the nucleus.^{2}^{2}2Currently, there is no coherent way to treat the PDF nuclear modifications within the Monte-Carlo Glauber model. This is actually why, in this work, we stick to the optical version of the Glauber model. Here we use EPS09s nuclear modifications from Ref. Helenius:2012wd () in which

(8) |

where the coefficients are obtained by analyzing the -dependence of the EPS09 Eskola:2009uj () nPDFs. The neutron PDFs are obtained from the proton PDFs by the isospin symmetry. By combining all, Eq. (2) factorizes into purely geometric and purely momentum-dependent parts which can be evaluated separately thereby reducing the dimensions of the required numerical integrations. We use the Incnlo program Aurenche:1987fs (); Aversa:1988vb (); Aurenche:1998gv (); Aurenche:1999nz () to calculate the momentum-dependent parts at next-to-leading order in perturbative QCD.

## 3 Results

### 3.1 Direct-photon production

Direct photons are produced either in the hard process or by the fragmentation of high- partons from the hard process. To obtain the latter contribution we convolute the partonic spectra with the BFG (set II) parton-to-photon fragmentation functions (FFs) Bourhis:1997yu (). Since the photon coupling is stronger in the case of up-type quarks than with the down-type quarks, the production rate of direct photons is larger in p+p collisions than in n+n collisions. This leads to a lower per-nucleon rate of direct photons in heavy-ion collisions than in p+p collisions due to the presence of neutrons. This is often referred to as the isospin effect and it becomes important at large values of where the valence quarks dominate. Furthermore, since the relative fraction of neutrons grows towards the edge of nucleus due to the NS effect, an additional suppression of direct photons in peripheral collisions is expected.

A canonical way to quantify the nuclear effects is to compute the nuclear modification factor, defined in a given centrality class of a Pb+Pb collision as

(9) |

where the normalization is related to the amount of interacting nuclear matter,

(10) |

where the impact parameters and define the centrality class as in Eq. (2). For the cross-section calculations we have set the renormalization, factorization and fragmentation scales to photon . The uncertainties related to the scale ambiguities are not considered here in more detail since they largely cancel out in the ratio, especially at large values of relevant here Helenius:2013bya (). The isolation criterion, often used by experiments to suppress secondary photons from hadronic decays, is not applied here since the effect to is negligible at the very high values of considered here.^{3}^{3}3The valence quark-gluon channel dominates irrespectively of the isolation.

Figure 1 shows at mid-rapidity for two centrality classes, 70–80 % and 90–100 % with and without the NS effect, compared also to the minimum-bias (0–100 %, MB) result and to the isospin effect. The uncertainties considered here are the EPS09s uncertainty (light colour band) and the one related to the uncertainty of neutrons 2pF parameters (dark colour band), obtained by evaluating with the quoted parameter variations, and adding the differences to the central prediction in quadrature.

The different regions are sensitive to different nPDF effects. First, comparing the MB result to the result without nPDF effects (only isospin), at some suppression due to shadowing is observed which then turns into an enhancement due to anti-shadowing. At a suppression due to the EMC effect is observed. The spatial dependence of the nPDFs always decreases the nuclear effects towards more peripheral collisions whereas the NS effect generates additional suppression with increasing . Therefore, at high values of where the impact of NS is more pronounced, these two effects pull towards opposite directions thereby “softening” the aggregate centrality dependence.

As the nuclear modifications of the PDFs gradually disappear with increasing peripherality, also the uncertainty becomes smaller for more peripheral events. However, even in the 90–100 % bin the nPDF uncertainty is of the same order as the NS effect which further complicates the separation of different effects. The non-zero nPDF effects even at the most peripheral bin are due to the power series ansatz in EPS09s, see Eq. (8), which, by construction, gives zero nuclear modifications only when . The uncertainty from the neutron 2pF parametrization turns out to be rather small as the contribution of the n+n channel is inferior e.g. to the contribution of the p+p channel and thus the variations in the neutron density are not that important. Since the fraction of neutrons (and therefore the n+n channel contribution) grows towards more peripheral collisions, also the uncertainty grows accordingly.

The centrality dependence can also be studied using the central-to-peripheral ratio defined as

(11) |

The advantage is that there is no need for a separate p+p baseline measurement and also that some uncertainties are expected to cancel out. Here, we have used the bin 0–10 % as the central result and compared it to the 70–80 % and 90–100 % bins. The results are shown in Fig. 2 again with and without the NS effect. Since the peripheral bins are now in the denominator, the NS effect increases the ratio and therefore decreases the centrality dependence at high- region ( closer to unity). The nPDF originating uncertainties are now larger with the most peripheral bin red (90–100 %) because the uncertainties in the central bin (similar to the MB uncertainty in Fig. 1) do not cancel here as effectively as with the less peripheral bin (70–80 %). Even though the interpretation of this observable is easier, the NS effect is still of the same order as the nPDF uncertainties.

At forward/backward rapidities (the lower panels in Fig. 2) the nPDF uncertainties are smaller. This is because here the dominant contribution comes from +g initial state where the gluon is at shadowing region with only mild uncertainty (at high factorization scale), and the quark is also a well-constrained high- valence quark. The modifications, however, are quite small and since there is an additional uncertainty due to modelling of the spatial dependence of the nPDFs, it is difficult to unambiguously study the NS effect with this observable. The most accurate centrality-dependent measurement for photons in Pb+Pb comes from the ATLAS collaboration Aad:2015lcb (). However, their most peripheral bin 40–80 % is still too central, and also the experimental uncertainties are large, to see any effects of NS.

### 3.2 Charged-hadron production

An observable in which the nPDF effects should cancel out very efficiently but yet be sensitive to the NS effect, is the ratio between negatively and positively charged hadrons,

(12) |

Since the relative number of neutron-involving (p+n, n+p, n+n) collisions is higher in peripheral than in central collisions, the increased d-quark contribution produces less positively charged hadrons and more negatively charged hadrons during the fragmentation. Here we do not consider any additional final-state effects that may affect the hadron production even though a significant suppression for the production of high- hadrons has been observed Abelev:2014laa (); Adam:2015kca (); Khachatryan:2016odn () in all centralities. Indeed, the measurements in Refs. Abelev:2014laa (); Adam:2015kca () show that the suppression at high () is very similar for all light charged hadrons (pions, kaons, protons) and, consequently, the particle ratios and are the same in p+p and Pb+Pb collisions. This motivates us to conjecture that final-state effects would have only a relatively small influence on the ratio of Eq. (12). Moreover, at the very high- region () where the current measurements are still statistically limited Khachatryan:2016odn (), the suppression effect in peripheral bins may be even negligible.

The cross section for hadron production is calculated by convoluting the partonic spectra with non-perturbative parton-to-hadron FFs. We consider three options, dss deFlorian:2007ekg (), kretzer Kretzer:2000yf () and akk08 Albino:2008fy (). To better understand the variations seen using different FFs, the ratios in p+p and n+n collisions at are shown in Fig. 3. The first observation is that with akk08 FFs the ratio in p+p actually turns negative at high-, caused by the cross section for becoming negative. This clearly unphysical result implies that the considered kinematic region is out of the validity region of akk08. The results using dss and kretzer are not that different in p+p collisions but for n+n collisions almost a factor of two difference at the very highest values of is observed. These differences between the FF analyses generate some further theoretical uncertainty for the considered observable. Turning this around, a measurement of in p+Pb or Pb+Pb collisions would clearly provide additional constraints for future FF analyses (modulo the possible final-state effects in Pb+Pb).

The ratios in Pb+Pb collisions at for centrality classes 70–80 % and 90–100 % are shown in the upper right panel of Fig. 3 together with the MB result with dss and kretzer FFs. The nPDF effects, including the centrality dependence and the uncertainties, are found to be negligible as expected. The uncertainty in the 2pF parametrization is negligible for the MB case but increases towards more peripheral collisions. The uncertainties are larger than in the case of direct-photon production as the cross section for now gets a large contribution from the n+n channel and thus carries more sensitivity to the parameter uncertainties in the neutron density. More importantly, the centrality dependence from the NS effect is clearly visible in this observable. However, the different FFs still yield rather different results but normalizing the ratio with the MB result, the FF dependence largely cancels out. This is demonstrated in the lower panels of Fig. 3, where the ratios in 70–80 % and 90–100 % classes are normalized with the 0–100 % result for and . Some FF dependence persists with but it is still smaller or of the same order than the uncertainty in 2pF parametrization. Also the nPDFs yield a few-percent uncertainty for the observable.

To estimate the achievable experimental precision for the ratios discussed above, we multiply the cross sections from Eq. (2) by the nominal Pb-Pb nucleon–nucleon luminosity of (also with which would correspond to the luminosity targeted after the LHC Long Shutdown 2 highlumi ())
^{4}^{4}4In these rough estimates we do not consider in detail the uncertainties related e.g. to the use of different FFs, next-to-NLO corrections, choices for the fragmentation/factorization/renormalization scales, or suppression of the hadron yields in Pb+Pb relative to p+p baseline.
.
From the resulting number of events we compute the relative statistical uncertainty by . To have better statistics, we consider here the rapidity bin (excluding the mid-rapidity to have a larger effect from the neutron skin) and suitably wide intervals. The results are shown in Fig. 4 where the total statistical uncertainty follows from combining the statistical uncertainties for and quadratically. As can be seen from Fig. 4, the realistically measurable region is for the 70–80 % bin and for the 90–100 % bin. We would expect that systematic uncertainties are small in these ratios like they are e.g. in the case of W charge asymmetry.

## 4 Summary and Outlook

We have studied the impact of the NS effect to direct-photon and charged-hadron production in Pb+Pb collisions at the LHC. In the case of photon production the NS effect has a 5–10 % impact on though the uncertainties in the nPDFs and their spatial dependence are of the same order or even larger than the expected effect. With some of the nPDF uncertainties cancel out making the NS effect more transparent. Also, going to larger rapidities decreases the nPDF uncertainties, but still the smallness of the NS effect and the ambiguities due to the centrality dependence of the nPDFs makes the direct-photon production a challenging observable to study the NS effect.

A more promising observable is the ratio between negatively and positively charged high- hadrons, for which we find up to 10 % effects in the statistically relevant region. In this case, the spatial dependence of the nPDFs cancel out very efficiently and, in general, the NS effect has a more pronounced impact than in the case of direct photons. The downsides here are the sensitivity to the applied fragmentation functions and, towards smaller , possible final-state modifications due to the produced strongly interacting medium. The first one can be cured by normalizing the ratio with the minimum bias result, but for a more detailed study of the latter, further modelling would be required. However, as discussed, there are indications that the final-state effects may largely disappear when considering particle ratios like the ones we have done here and, after all, the disparity between the amount of initial-state up and down quarks should strongly correlate with the balance of produced negatively and positively charged hadrons, irrespectively of the exact way the produced hard partons hadronize.

We hope that in near future the NS effect could provide an additional handle to control the centrality classification and help to bridge the theoretical and experimental centrality definitions. As a further prospect, we plan to study the NS effect in the future high-luminosity lepton-ion colliders.

###### Acknowledgements.

I. H. has been supported by the MCnetITN FP7 Marie Curie Initial Training Network, Contract PITN-GA-2012-315877 and has received funding from the European Research Council (ERC) under the European Unionâs Horizon 2020 research and innovation programme (Grant Agreement No 668679). This research was supported by the European Research Council Grant HotLHC ERC-2011-StG-279579 and by Xunta de Galicia (Conselleria de Educacion)–H. P. is part of the Strategic Unit AGRUP2015/11.## References

- (1) B. Abelev et al. (ALICE Collaboration), JHEP 1403 (2014) 013 doi:10.1007/JHEP03(2014)013 [arXiv:1311.0633 [nucl-ex]].
- (2) S. Chatrchyan et al. (CMS Collaboration), Phys. Lett. B 712 (2012) 176 doi:10.1016/j.physletb.2012.04.058 [arXiv:1202.5022 [nucl-ex]].
- (3) G. Aad et al. (ATLAS Collaboration), Phys. Rev. Lett. 114 (2015) no.7, 072302 doi:10.1103/PhysRevLett.114.072302 [arXiv:1411.2357 [hep-ex]].
- (4) G. Aad et al. (ATLAS Collaboration), Eur. Phys. J. C 74 (2014) no.8, 2982 doi:10.1140/epjc/s10052-014-2982-4 [arXiv:1405.3936 [hep-ex]].
- (5) S. Chatrchyan et al. (CMS Collaboration), Phys. Rev. Lett. 109 (2012) 022301 doi:10.1103/PhysRevLett.109.022301 [arXiv:1204.1850 [nucl-ex]].
- (6) S. Chatrchyan et al. (CMS Collaboration), Phys. Rev. C 84 (2011) 024906 doi:10.1103/PhysRevC.84.024906 [arXiv:1102.1957 [nucl-ex]].
- (7) B. Abelev et al. (ALICE Collaboration), Phys. Rev. C 88 (2013) no.4, 044909 doi:10.1103/PhysRevC.88.044909 [arXiv:1301.4361 [nucl-ex]].
- (8) M. L. Miller, K. Reygers, S. J. Sanders and P. Steinberg, Ann. Rev. Nucl. Part. Sci. 57 (2007) 205 doi:10.1146/annurev.nucl.57.090506.123020 [nucl-ex/0701025].
- (9) S. Chatrchyan et al. (CMS Collaboration), Eur. Phys. J. C 74 (2014) no.7, 2951 doi:10.1140/epjc/s10052-014-2951-y [arXiv:1401.4433 [nucl-ex]].
- (10) G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 748 (2015) 392 doi:10.1016/j.physletb.2015.07.023 [arXiv:1412.4092 [hep-ex]].
- (11) A. Adare et al. (PHENIX Collaboration), Phys. Rev. Lett. 116 (2016) no.12, 122301 doi:10.1103/PhysRevLett.116.122301 [arXiv:1509.04657 [nucl-ex]].
- (12) G. Martinez-Garcia, arXiv:1408.3108 [hep-ph].
- (13) A. Bzdak, V. Skokov and S. Bathe, Phys. Rev. C 93 (2016) no.4, 044901 doi:10.1103/PhysRevC.93.044901 [arXiv:1408.3156 [hep-ph]].
- (14) M. Alvioli, B. A. Cole, L. Frankfurt, D. V. Perepelitsa and M. Strikman, Phys. Rev. C 93 (2016) no.1, 011902 doi:10.1103/PhysRevC.93.011902 [arXiv:1409.7381 [hep-ph]].
- (15) D. V. Perepelitsa and P. A. Steinberg, arXiv:1412.0976 [nucl-ex].
- (16) N. Armesto, D. C. GÃ¼lhan and J. G. Milhano, Phys. Lett. B 747 (2015) 441 doi:10.1016/j.physletb.2015.06.032 [arXiv:1502.02986 [hep-ph]].
- (17) C. M. Tarbert et al., Phys. Rev. Lett. 112 (2014) no.24, 242502 doi:10.1103/PhysRevLett.112.242502 [arXiv:1311.0168 [nucl-ex]].
- (18) M. B. Tsang et al., Phys. Rev. C 86 (2012) 015803 doi:10.1103/PhysRevC.86.015803 [arXiv:1204.0466 [nucl-ex]].
- (19) J. Zenihiro et al., Phys. Rev. C 82 (2010) 044611. doi:10.1103/PhysRevC.82.044611
- (20) C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86 (2001) 5647 doi:10.1103/PhysRevLett.86.5647 [astro-ph/0010227].
- (21) H. Paukkunen, Phys. Lett. B 745 (2015) 73 doi:10.1016/j.physletb.2015.04.037 [arXiv:1503.02448 [hep-ph]].
- (22) I. Helenius, K. J. Eskola, H. Honkanen and C. A. Salgado, JHEP 1207 (2012) 073 doi:10.1007/JHEP07(2012)073 [arXiv:1205.5359 [hep-ph]].
- (23) G. Antchev et al. (TOTEM Collaboration), Europhys. Lett. 101 (2013) 21004. doi:10.1209/0295-5075/101/21004
- (24) H. L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin and C.-P. Yuan, Phys. Rev. D 82 (2010) 074024 doi:10.1103/PhysRevD.82.074024 [arXiv:1007.2241 [hep-ph]].
- (25) K. J. Eskola, H. Paukkunen and C. A. Salgado, JHEP 0904 (2009) 065 doi:10.1088/1126-6708/2009/04/065 [arXiv:0902.4154 [hep-ph]].
- (26) P. Aurenche, R. Baier, M. Fontannaz and D. Schiff, Nucl. Phys. B 297 (1988) 661. doi:10.1016/0550-3213(88)90553-6
- (27) F. Aversa, P. Chiappetta, M. Greco and J. P. Guillet, Nucl. Phys. B 327 (1989) 105. doi:10.1016/0550-3213(89)90288-5
- (28) P. Aurenche, M. Fontannaz, J. P. Guillet, B. A. Kniehl, E. Pilon and M. Werlen, Eur. Phys. J. C 9 (1999) 107 doi:10.1007/s100529900018 [hep-ph/9811382].
- (29) P. Aurenche, M. Fontannaz, J. P. Guillet, B. A. Kniehl and M. Werlen, Eur. Phys. J. C 13 (2000) 347 doi:10.1007/s100520000309 [hep-ph/9910252].
- (30) L. Bourhis, M. Fontannaz and J. P. Guillet, Eur. Phys. J. C 2 (1998) 529 doi:10.1007/s100520050158 [hep-ph/9704447].
- (31) I. Helenius, K. J. Eskola and H. Paukkunen, JHEP 1305 (2013) 030 doi:10.1007/JHEP05(2013)030 [arXiv:1302.5580 [hep-ph]].
- (32) G. Aad et al. (ATLAS Collaboration), Phys. Rev. C 93 (2016) no.3, 034914 doi:10.1103/PhysRevC.93.034914 [arXiv:1506.08552 [hep-ex]].
- (33) B. B. Abelev et al. (ALICE Collaboration), Phys. Lett. B 736 (2014) 196 doi:10.1016/j.physletb.2014.07.011 [arXiv:1401.1250 [nucl-ex]].
- (34) J. Adam et al. (ALICE Collaboration), Phys. Rev. C 93 (2016) no.3, 034913 doi:10.1103/PhysRevC.93.034913 [arXiv:1506.07287 [nucl-ex]].
- (35) V. Khachatryan et al. (CMS Collaboration), arXiv:1611.01664 [nucl-ex].
- (36) D. de Florian, R. Sassot and M. Stratmann, Phys. Rev. D 76 (2007) 074033 doi:10.1103/PhysRevD.76.074033 [arXiv:0707.1506 [hep-ph]].
- (37) S. Kretzer, Phys. Rev. D 62 (2000) 054001 doi:10.1103/PhysRevD.62.054001 [hep-ph/0003177].
- (38) S. Albino, B. A. Kniehl and G. Kramer, Nucl. Phys. B 803 (2008) 42 doi:10.1016/j.nuclphysb.2008.05.017 [arXiv:0803.2768 [hep-ph]].
- (39) A. Uras, PoS LHCP 2016 (2016) 177