Detecting nuclear mass distribution in isobar collisions via charmonium

The collective properties of final state hadrons produced in the high statistics 4496\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{44}^{96}$$\end{document}Ru + 4496\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{44}^{96}$$\end{document}Ru and 4096\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{40}^{96}$$\end{document}Zr + 4096\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{40}^{96}$$\end{document}Zr collisions at sNN=200GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s_\textrm{NN}} = 200~\textrm{GeV}$$\end{document} are found to be significantly different. Such differences were argued to be precise probes of the difference in nucleon distribution in the isobar nuclei. We investigate the J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} production in the isobar collision via a relativistic transport approach. By comparing the isobar systems according to equal centrality bin and equal multiplicity bin, we find that the yield ratio of J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} is sensitive to the differences in both the number of binary collisions and the medium evolution. Besides, the elliptic flow v2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_2$$\end{document} of J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J/\psi $$\end{document} is qualitatively different from the light hadrons, and the ratio between Ru + Ru and Zr + Zr collisions is sensitive to the medium evolution. The charmonium production provides an independent probe to study the nucleon distribution in the isobar system.

Introduction.-At the mean-field level, nuclear properties like spin, parity of ground state, magic number, and β decays can be well described by the nuclear shell model [1].In the shell model, we know most of the nuclei in their ground state are deformed and non-uniform, except the Doubly-Magic Nucleus.Under the polarization effect of valence nucleons, nuclei with an unfilled shell or sub-shell will be deformed and generates collective motion induced by the interaction between valence nucleons and shell structure [2].The deformation depends closely on the number of protons and neutrons in the nucleus and has been studied in low-energy nuclear physics for many years, e.g., see reviews [3,4].In relativistic heavyion collisions, the nuclei are fully destroyed at the beginning of the collisions, but the deformed information will still show its fingerprint in the final observables [5,6].The deformed nuclear-nuclear collisions will lead to an anisotropic overlap region.The spatial anisotropy of the collision overlap region in the initial state will transform into a momentum anisotropy of the produced hadrons in the final state [7][8][9][10][11][12].Therefore, relativistic heavy-ion collisions may supply a new way to probe nuclear deformation.
Recently, a high statistics heavy-ion collision, which collides 96  44 Ru+ 96 44 Ru and 96 40 Zr+ 96 40 Zr with beam energy √ s NN = 200 GeV, is performed by the STAR Collaboration at the Relativistic Heavy Ion Collider (RHIC) [13] 1 .Such a contrast experiment originally proposed to search for the chiral magnetic effect (CME) [14][15][16], under the presumption that the same baryon number would lead to the same non-CME background while the different electric number would induce a sizable difference in the CME signal.Unexpectedly, the bulk properties, namely the charge multiplicity and elliptic and triangular flows, are found to be different in the isobar systems [13].While such a subtlety prevents one to make a conclusive statement on the existence of CME in heavy-ion collisions and calls for more efforts in better quantification of the background [17,18], it brings a new opportunity to study the nucleon distribution in relativistic heavy-ion collisions [19][20][21][22][23][24][25][26][27][28][29][30][31][32].In such experiments, a new phase of matter, called Quark Gluon Plasma (QGP), is created, and the final state particles can be well explained by a hydrodynamic description of the system evolution, (see e.g.Refs.[33][34][35]).The initial state of the evolution is the consequence of multiple nucleon-nucleon collisions, and it provides an unique opportunity of measuring multinucleon correlation within a nucleus, regardless of the fact that the mapping between initial condition and final observables is complicated.In Ru and Zr nuclei, the spatial distribution of nucleons are usually described by a deformed Woods-Saxon distribution [36], where the parameters R 0 and a 0 are respectively called radius and diffusiveness.The maximum nuclear den-arXiv:2211.01971v2[hep-ph] 2 Jun 2023 sity ρ 0 is given by the normalization condition that ρ(r, θ, φ)d 3 r = A. Y ,m 's are spherical harmonics.β 2 and β 3 control the quadrupole and octupole deformations, respectively.Both lower-energy nuclear experiments and high-energy heavy-ion collisions suggest that the Ru nucleus has a sizable quadrupole deformation while the Zr nucleus has an octupole deformation [23,27,37,38].Particularly, it is found in [23] and [32], which respectively describes the medium evolution by a particle transport model and hydrodynamics, that the high-precision measurements in the isobar collisions can be described by the parameter set R Ru Charmonium states like J/ψ, which are boundstates of a charm quark and its anti-quark, are independent probes of the initial condition.Owing to the heavy mass, charm quarks and charmonium states are dominantly produced in the initial hard scattering in heavy-ion collisions.After the initial production, they interact with the QGP and their properties get modified.The spectrum and flow information of charmonium states in heavy-ion collisions reflect both the QGP properties and initial condition [39][40][41][42][43][44][45][46][47][48][49][50][51].In heavy-ion collisions, the initially produced charmonium bound states might dissociate into scattering states due to the static color-screening effect and dynamic dissociation [39,[52][53][54][55].Some of them can survive to the end of the QGP evolution and be detected.Another source of final state charmonium particle comes from the recombination of uncorrelated charm and anticharm quarks in the QGP [40][41][42][43][44][45][46][47].We note that initial production rates of both charm quarks and charmonium states are proportional to the number of binary nucleonnucleon collisions, denoted as N coll .In the hypothetical situation that the bulk background being the same in the isobar system, the survived charmonium production rate is proportional to N coll , whereas the recombination production rate is proportional to the square of charm quark number, hence ∝ N 2 coll .Compared to light flavor observables, charmonium states are sensitive to different aspects of the initial state and bulk evolution, and henceforth serve as an independent probe of the QGP.In this work, we study the differences between the properties of the J/ψ particles produced in Ru+Ru and those in Zr+Zr collisions.We aim to provide an independent probe for the deformations in Ru and Zr nuclei.
Charmonium transport.-Westart by describing the charmonium dissociation and recombination in the QGP.Since charmonium states are heavy and color neutral, they are unlikely to be thermalized with the medium.Their phase space distribution, f ψ (p, x, τ ) with ψ ∈ {J/ψ, χ c , ψ }, is governed by a transport equation which includes both initial production and regeneration [44][45][46][47], where y = (1/2) ln[(E + p z )/(E − p z )] is the momentum rapidity, v T = p T /E T is the transverse velocity, T is the transverse energy, and ∇ T ≡ (∂ x , ∂ y ) is the transverse gradient.The second(third) term on the left hand side arises from the free streaming of ψ which leads to the leakage effect in the longitudinal(transverse) direction.The anomalous suppression and regeneration mechanisms in the QGP medium are respectively reflected in the loss term α and gain term β.Charmonia in hot QGP medium suffer Debye screening [39].With increasing temperature, the interaction between a pair of heavy quarks is more screened, while the averaged size of a charmonium state increases.When the averaged size r diverges, the charmonium is dissociated.We solve the two-body Schrödinger equation with the finite-temperature potential and define the dissociation temperature T d as the temperature when r → ∞.We find T d = (2.3,1.2, 1.1) T c for J/ψ, χ c , and ψ , respectively [56].When the local temperature is higher than the dissociation temperatures T d , the charmonium state disappears, and the regeneration only happens when the temperature is lower than T d .In addition to the Debye screening, charmonium states suffer dynamical dissociation in the QGP, such as gluon dissociation process, g + ψ → c + c [53,54].The charmonium gluodissociation cross-section σ cc gψ can be derived, via the operator-production-expansion (OPE) method, in vacuum [53,54] and extended to the finite temperature [45][46][47].Here we adopt the setup of [45].Taking only the gluo-dissociation as the loss term and its inverse process, c + c → g + ψ, as the gain term in this study, α and β can be explicitly expressed as [44] where E g , E c , and E c are the gluon, charm quark and anti-charm quark energies, respectively, and p g , p c and p c are their momenta.The Mandelstam variable s is the gψ interaction energy.W cc gψ is the dissociation probability while W gψ cc the regeneration probability.Their related due to the detailed balance between the two processes, and the explicit expressions are given in [57].Θ is the Heaviside step function to guarantee the calculation in the QGP phase above the critical temperature T c .The latter is taken to be T c = 160 MeV.The gluon distribution f g follows the Bose-Einstein distribution, whereas charm quarks with low transverse momenta are found to be thermalized according to experimental measurement of open-charm meson [58], despite of their large mass.As a first-order approximation, we take a thermal distribution for charm/anti-charm momentum, , where N c is the normalization factor to ensure Nc e p µ c uµ/T +1 Considering the small transverse velocity and the similarity of the isobar system, we neglect the density change due to transverse expansion, and the density in coordinate space is governed by the nuclear geometry of the colliding system [45], where dσ cc pp /dy is charm quark production cross section per unit rapidity in p+p collisions.We take dσ cc pp /dy = 0.162 mb, which corresponds to √ s NN = 200 GeV [56].
The binary collision density ρ coll (x) is given by the initial condition of the bulk background.In this work, we adopt the boost invariant, event-by-event fluctuating initial conditions generated by the Monte Carlo Glauber model (MCGlauber) [59].We generate millions of initial profiles and bin them into different centrality or multiplicity classes.Events within the same class are then averaged after being overlapped according to the center of mass and aligned according to the second-order participant plane.With such event-averaged initial conditions, we solve the medium evolution taking the data validated music hydrodynamic simulation package [33,34,60,61], with taking the s95p Equation of State [62], a constant shear viscosity η/s = 0.08, and vanishing bulk viscosity [63][64][65].The hydrodynamic evolution provides the space-time profile of temperature (T ) and flow velocity (u µ ) of the bulk background in the calculation of charmonium transport and gluon and charm quark distribution.The centrality and multiplicity dependence of the maximum temperature of the medium produced in Ru+Ru and Zr+Zr collisions are shown in the upper panels of Fig. 2. One can find them to be close to identical in central collisions but visibly different in peripheral collisions.
Besides the hot nuclear matter effect which affects charmonium motion through the above discussed anomalous suppression and regeneration in the QGP, there is also the cold nuclear matter effect which changes the initial condition of the transport equation (2).The cold nuclear matter effect includes mainly the nuclear absorption [66], Cronin effect [67], and nuclear shadowing [68].Considering a finite formation time of charmonia and small nuclear size of Ru and Zr in high energy collisions, the nuclear absorption is neglected in this study.Such an approximation can be further justified by the fact that the absorption effect shall be similar for both nuclei.The initial charmonium distribution in heavy-ion collisions is constructed by a superposition of the charmonium distribution in p+p collisions with considering the cold nuclear matter effects, where the differential cross section [56,69], encodes the Cronin effect that enhance the transverse momentum compared to p+p collisions The parameters a = 2b 2 )/Γ(n − 1), n = 3.93, and p 2 T p = 3.05 GeV 2 are fitted to match the experimental data for p+p collisions [56,69,70].The Cronin parameter, which describes the averaged charmonium transverse momentum square obtained from the gluon scattering with a unit of length of nucleons, is set as a gN = 0.1 GeV 2 /fm [41,56] in order to fit experimental data in p-A and low-energy A-A collisions [71,72].The mean trajectory length of the two gluons in the two nuclei before the cc formation has been estimated as ρ −1 0 ρ coll (x)/σ inel NN , where ρ 0 is the maximum nuclear density in Eq. ( 1).The direct production cross sections dσ {J/ψ,χc,ψ } pp /dy = {0.6,0.3, 0.1} × dσ J/ψ,prompt pp /dy [73], where dσ J/ψ,prompt pp /dy = 716.7 nb is prompt production cross section [70].
R P g and R T g in Eq. ( 5) are spatial-dependent shadowing factors for gluon in projectile and target.Assuming that the inhomogeneous shadowing is proportional to the parton path length through the nucleus, the spatialdependent shadowing factor can be estimated as, where the shadowing modification factor R g can be simulated by the EPS09 package [74].
x P,T g = ( m 2 ψ + p 2 T / √ s NN )e ±y is gluon longitudinal momentum fraction with y being the charmonium rapidity, and the +(−) sign is chose for projectile(target).The momentum transfer squared can be taken as For both Ru and Zr, the nuclear mass number is A = 96.N cen coll is the number of binary collision for the most central collisions.
Results.-Experimental measurement of prompt J/ψ mesons include feed-down contributions from the excited charmonium states, i.e., χ c and ψ .With initial condition given in Eq. ( 5), we solve the transport equations (2) for J/ψ, ψ and χ c on top of event-averaged hydrodynamic background corresponding to different collision systems and different centrality or multiplicity bins.Then we decay χ c and ψ into J/ψ, with branch ratios respectively being 22% and 61% [75], and obtain the properties of prompt J/ψ in the final state.It shall be noted that the non-prompt contribution from B-decays can be neglected at 200 GeV Ru+Ru and Zr+Zr collisions, due to the small production cross-section of the latter at low beam energy [73].Therefore, inclusive production rate of J/ψ is dominated by the prompt one.
We first study the nuclear modification factor of J/ψ in both Ru+Ru and Zr+Zr collisions, which is defined as the ratio of production rate in AA to that in p+p collisions, scaled by the inverse of number of binary collisions, In Fig. 2, we present both the centrality and multiplicity dependence of prompt J/ψ R AA in Ru+Ru and Zr+Zr collisions at √ s NN = 200 GeV, with contributions from initial production and regeneration separated.We observe a clear cusp located the 5 − 10% centrality bin or the ∼ 230 multiplicity bin.For collisions more central than the cusp, the highest temperature is higher than the dissociation temperature of J/ψ, T d = 2.3 T c , which is indicated by the horizontal dashed lines.J/ψ mesons initially produced at the T > T d region suffer a strong suppression.Meanwhile, one can find the regeneration contribution plans an important role in final total yield, especially at central collisions even in small systems like Ru+Ru and Zr+Zr collisions.
Then we move on to compare the J/ψ production in the isobar systems by taking the ratio between observables in Ru+Ru and those in Zr+Zr collisions.In the upper panels of Fig. 3, we show the ratio of total yield as well as those for initial and regeneration productions.In the hypothetical case of identical background in the isobar systems, one would expect the ratio of initial production (thin solid lines) to be the same as that of N coll (filled square), whereas the ratio of regeneration production (thin dash lines) to be identical to that of N2 coll (open square) 2 .We observe that the former is below the N coll ratio, because higher temperature and larger volume in Ru+Ru collisions.Meanwhile, the ratio of regeneration production agree well with the N 2 coll ratio at central collisions, and goes above such a expectation for peripheral collisions owning to the longer evolution time in Ru+Ru collisions, and therefore more J/ψ particles are produced via regeneration.As for the double ratio, by definition, it equals to the yield ratio times the inverse ratio of the N coll .So, the double ratio is smaller than 1 and decreases monotonously at peripheral collisions.In comparing the left and right panels, one can clearly see that all ratios are closer to unity when we compare the isobar system within the same multiplicity bin than the comparison for the same centrality bin.
The nucleon distribution in the isobar nuclei influences the anisotropy of the medium created in heavy-ion collisions.The latter can be measured by the momentum anisotropy of final state hadrons, which is defined as the Fourier coefficient of their azimuthal distribution and referred to as v n [7].The second flow coefficient, called elliptic flow v 2 , is mostly a collective flow response to the ellipsoidal shape of the overlap region in non-central collisions.The difference in the light flavor v 2 has been found to be sensitive to the initial anisotropy contributed by the deformation [22,23], and it would be interesting to measure v J/ψ 2 which probes the medium anisotropy in a different way.Results are shown in Fig. 4. The v 2 of J/ψ is influenced by two origins.One is regeneration in which the anisotropy is similar to that of π in both the centrality trend and order of magnitude.The other is the initial production which is isotropic initially and pick up some small amount of anisotropy via path-dependent dissociation in the anisotropic medium background.In central collisions, the anisotropy of the background is small, whereas in peripheral collisions the temperature is small and the dissociation effect is weak.Both these factors make v 2 of initially produced J/ψ flat and small for all centrality ranges.As is shown in Fig. 2, going from the most central to peripheral collisions, the regeneration fraction decreases and henceforth the overall v 2 of J/ψ decreases, which is opposite to the trend of light hadrons.A small tilde (enhancement) in v 2 for the most central collisions is also observed, because the existence of T > T d region induces extra anisotropy for initially produced J/ψ.With the Ru-to-Zr ratio shown in the lower panel, we expect ∼ 5 − 10% difference in the elliptic flow of J/ψ.It shall be worth noting that such a ratio is greater than that of v 2,π .With details shown in the Appendix, the difference of v 2,J/ψ in the isobar system is contributed dominantly by the fact that there is a greater portion of J/ψ that is created via regeneration production in Ru+Ru collisions than in Zr+Zr.Thus, a precise measurement of v 2,J/ψ difference serves as an accurate probe of the regeneration fraction.
Summary.-Via a relativistic transport approach, we investigate the yield and elliptic flow of J/ψ in isobar collisions.The evolution of the hot medium produced in Ru+Ru and Zr+Zr collisions are described by relativistic viscous hydrodynamics with event averaged initial conditions.Aiming to detect the difference in nucleon distributions between Ru and Zr, we compare the observables in the isobar system according to the same centrality bin and the same multiplicity bin.Due to the charm produced via initial binary collisions, the yield ratio can clearly show the difference between Ru and Zr.The elliptic flow v 2 of J/ψ in isobar collisions is qualitatively different from the light hadrons, and the ratio between Ru+Ru and Zr+Zr collisions is around 5 − 10% in noncentral collisions.These heavy flavor particles serve as an independent probe to characterize the number of binary collisions in the initial state and therefore the nucleon distribution in the isobar nuclei.
The regeneration production of charmonium is proportional to the square of the production cross-section of charm quark, with the latter being proportional to the number of binary collision and also increasing with collision energy.At LHC energy charmonium states are produced dominantly via regeneration and we expect the difference in nucleon distribution to be better reflected there.
In addition, the deformation of the nucleus is better reflected in ultra-central collisions, in which one may select events with preference on the so called body-body collisions.In such events, the deformation axes of the projectile and target are parallel to each other, and both of them are perpendicular to the beam direction.Therefore, the deformation of the nucleon distribution is translated into the eccentricity of the initial condition, and henceforth the final state anisotropy.Comparison of both light and heavy flavor observables will be reported in our future publication.(A1) In the left panel of Fig. 5, we plot these three sources separately, together with the full value of ∆v tot 2 .It turns out that the last term, which is caused by the fact that there is a greater portion of J/ψ that is created via regeneration production in Ru+Ru collisions than in Zr+Zr.As noted in the main text that v 2,reg is much greater than v 2,ini , the difference in f reg makes the greatest contribution to the difference in elliptic flow.

FIG. 2 .
FIG. 2. (Left) The centrality dependence of the maximum temperature of the medium (upper), nuclear modification factor RAA of J/ψ in Ru+Ru (middle) and Zr+Zr (lower) collisions at √ sNN = 200 GeV.Vertical dashed lines indicate the position of the cusp, whereas the horizontal dashed line label the dissociation temperature of J/ψ.(Right) Same as Left but for equal-multiplicity bins.

FIG. 4 .
FIG. 4. (Upper panel) The centrality dependence of the elliptic flow v2 of light hadrons π and J/ψ produced in Ru+Ru and Zr+Zr collisions at √ sNN = 200 GeV.(Lower panel) The ratio of v2 in Ru+Ru and Zr+Zr collisions for J/ψ and π ± .