Arxiv:1908.10231V2 [Nucl-Th] 25 Oct 2019 Could Disentangle the Background from the Signal
Total Page:16
File Type:pdf, Size:1020Kb
Influence of the neutron-skin effect on nuclear isobar collisions at RHIC Jan Hammelmann,1, 2, ∗ Alba Soto-Ontoso,3, y Massimiliano Alvioli,4, z Hannah Elfner,5, 1, 2, x and Mark Strikman6, { 1Goethe-Universit¨atFrankfurt, Max-von-Laue-Strasse 1, 60438 Frankfurt am Main, Germany. 2Frankfurt Institute for Advanced Studies, Ruth-Moufang-Strasse 1, 60438 Frankfurt am Main, Germany. 3Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA. 4Consiglio Nazionale delle Ricerche, Istituto di Ricerca per la Protezione Idrogeologica, via Madonna Alta 126, I-06128 Perugia, Italy. 5GSI Helmholtzzentrum f¨urSchwerionenforschung GmbH, Planckstr. 1, 64291 Darmstadt, Germany. 6Pennsylvania State University, University Park, PA 16802, USA. The unambiguous observation of a Chiral Magnetic Effect (CME)-driven charge separation is 96 96 96 96 pthe core aim of the isobar program at RHIC consisting of 40Zr+40Zr and 44Ru+44Ru collisions at sNN = 200 GeV. We quantify the role of the spatial distributions of the nucleons in the isobars on both eccentricity and magnetic field strength within a relativistic hadronic transport approach (SMASH, Simulating Many Accelerated Strongly-interacting Hadrons). In particular, we introduce isospin-dependent nucleon-nucleon spatial correlations in the geometric description of both nuclei, 96 96 deformation for 44Ru and the so-called neutron skin effect for the neutron-rich isobar i.e. 40Zr. The 96 96 main result of this study is a reduction of the magnetic field strength difference between 44Ru+44Ru 96 96 and 40Zr+40Zr by a factor of 2, from 10% to 5% in peripheral collisions when the neutron-skin effect is included. Further, we find an increase of the eccentricity ratio between the isobars by up to 10% 96 in ultra-central collisions as due to the deformation of 44Ru while neither the neutron skin effect nor the nucleon-nucleon correlations result into a significant modification of this observable with respect to the traditional Woods-Saxon modeling. Our results suggest a significantly smaller CME signal to background ratio for the experimental charge separation measurement in peripheral collisions with the isobar systems than previously expected. INTRODUCTION efficiently converted into a separation of positive and negative charges along the direction of the magnetic field. This phenomenon, dubbed Chiral Magnetic Effect One of the fundamental properties of Quantum Chro- (CME) [10, 11], manifests itself into charge-dependent modynamics (QCD) is the axial anomaly, which in the azimuthal correlations of the measured hadrons [12]. massless fermion limit, reads as follows 2 A decade after the pioneering analysis of the STAR µ g a a,µν Collaboration [13], the experimental confirmation of the @µj5 = − 2 Fµν Fe (1) 16π CME remains unsettled. Numerous charge separation a a,µν measurements in line with CME expectations were re- where Fµν is the gluon field strength, Fe its dual, g µ ported with different collisions systems and energies the strong coupling constant and j5 the axial current density. The axial anomaly establishes a direct relation- both from RHIC [14–16] and LHC [17, 18]. However, ship between the generation of a net axial charge and these measurements are known to be strongly affected the dynamics of non-Abelian gauge fields. by background contamination arising from flow [19] and local charge conservation [20]. New observables be- Together with condensed matter systems [1], ultra- yond the traditional three particle correlator could help relativistic heavy-ion collisions provide a unique envi- solving the problem [21]. Moreover, new RHIC mea- ronment to experimentally test the chiral anomaly. At surements with different isobars, i.e. 96Zr and 96R[22] least two mechanisms contribute to the right-hand side 40 44 arXiv:1908.10231v2 [nucl-th] 25 Oct 2019 could disentangle the background from the signal. For of Eq.1. On the one hand, in the Color Glass Conden- that purpose, 96Zr+96Zr collisions will provide a precise sate description of the early, non-equilibrium stage of 40 40 characterization of the background contribution to the the collision known as Glasma [2], fluctuations of the experimental charge separation measurement. On the chromo-electric and chromo-magnetic fields give rise to 96 96 a a,µν other hand, the proton-rich isobar system 44Ru+44Ru a non-vanishing Fµν Fe [3,4]. Further, the non-trivial will provide an enhanced sensitivity to the CME com- topological structure of the QCD vacuum results into ponent, due to the formation of larger magnetic fields. another source of net axial charge density known as A correct interpretation of the forthcoming experi- sphaleron transitions, whose rate is enhanced at high mental data requires accurate quantification of back- temperatures such as the ones reached in the Quark- ground and signal from theory. A multiphase trans- Gluon Plasma phase [5–7]. These local fluctuations of port model predicted the magnetic field strength, pro- axial charge density in the transverse plane occur in portional to the CME contribution, to be 10% larger the presence of a strong electromagnetic field in non- for 96Ru+96Ru than for 96Zr+96Zr in peripheral colli- central collisions [8,9]. Then, the chiral imbalance is 44 44 40 40 2 sions [23, 24]. A hydrodynamic framework predicted having either R0;p > R0;n; dn ∼ dp dubbed neutron-skin differences of up to 10% on the elliptic flow of both col- type or R0;p ∼ R0;n; dn > dp refer to as neutron-halo type. lision systems related to deformation [25]. A systematic A remarkable example of the latter category is 208Pb comparison between a Woods-Saxon shape and density with ∆rnp ∼ 0:15 fm [34]. In this case, the implications functional theory calculations shows that the functional of ∆rnp 6=0 in the context of observables relevant for the form of the nuclear density distributions used in the heavy-ion program at the LHC were recently studied simulations also impacts v2 ( ∼ 3%)[26]. All in all, the in [35–38]. Interestingly, one of the nuclei chosen for the 96 results of these studies identify the nuclear structure of isobar run at RHIC, 40Zr, also pertains to the neutron- the two isobar nuclei to be a source of uncertainty for v2 halo category with but not for the magnetic field strength. In this work, we analyze the effect due to an experi- ∆rnp = 0:12 ± 0:03 fm (5) 96Zr mentally measured nuclear phenomenon in the descrip- 40 96 tion of the density distribution of 40Zr i.e. the neutron- as extracted from the experimental analysis performed skin effect [27, 28]. This ingredient leads to an enhance- 96 96 with the Low Energy Antiproton Ring at CERN [27, 28]. ment of the magnetic field in peripheral 40Zr+40Zr col- The goal of this work is to study the consequences of lisions within SMASH [29] consequently undermining considering isospin-dependent Woods-Saxon distribu- the experimental prospects of finding out the Chiral tions fulfilling the upper limit of the constraint given Magnetic Effect with the isobar run. 96 by Eq.5, ∆rnp = 0:15 fm, to describe 40Zr on CME- related observables. Note that we take the upper limit of ∆rnp = 0:15 fm with the purpose of studying the FRAMEWORK neutron-skin impact at its extreme. Neutron-skin effect and nucleon-nucleon correlations 1.2 Traditionally, the spatial distribution of nucleons in- 1.0 ) side nuclei is generated by randomly sampling the r ( 0.8 ) 1 ( Woods-Saxon density distribution [30] n / 0.6 ρ0 ) r ρ(r; θ) = 0 (2) ( (r−R (θ,φ))=d ) e + 1 1 0.4 ( p where 0.2 Zirconium w neutron skin 0 0 R (θ) = R0(1 + β2Y2 (θ)): (3) Zirconium w/o neutron skin 0.0 In Eqs.2-3, ρ =0:168 is the ground state density, d refers 0 2 4 6 8 10 0 r [fm] to the difussiveness, R0 is the nuclear radius and β2 to- 0 gether with the spherical harmonic Y2 control the defor- Figure 1. Ratio of one-body density of protons to neutrons as mation. Two severe simplifications are commonly made 96 when using Eq.2 for the nuclear geometry. First, nu- a function of the radial distance for 40Zr with (solid) or with- out (dashed) considering the neutron skin effect. Error bars cleons are considered to be independent of each other. account for statistical uncertainties. Second, protons and neutrons are treated indistinctly so that they are sampled from the same Woods-Saxon dis- For that purpose, the starting point are the experi- tribution i.e. with identical values for R0 and d. Experi- mental values for (R0; d) of the charge distribution [39], mental measurements and theoretical calculations ruled displayed both for 96Zr and 96Ru in TableI. To extract out both assumptios, as discussed below. 40 44 the values of point distributions (R0;p(n); dp(n)) from the Since the early 80’s the tails of the proton (p) and neu- charged ones, keeping the size of the nucleus fixed, we tron (n) distributions are known to be distinct [31–33] follow the procedure outlined in [37, 40, 41]. It is not the i.e. R0 and d in Eq.2 are isospin dependent. The neu- goal of this paper to repeat the precise derivation de- tron distribution populates the outer region of neutron- tailed in the aforementioned papers. For completeness, rich nuclei. That is, the difference between the neutron in Appendix A we provide the formulae that were used and proton distributions mean square radii, which can to obtain the values of (R0;p(n); dp(n)) shown in TableII be written as follows: where we observe how ∆rnp =0:15 fm is translated into 2 1=2 2 1=2 a larger value of the diffusiveness for the neutron distri- ∆rnp = hr i − hr i ; (4) n p bution while the radius remain the same for both types is positive.