Predicting Quadrupole deformation via anisotropic flow and transverse momentum spectra in isotopic $\mathbf{\prescript{128-135}{54}{\mathrm{Xe}}}$ collisions at LHC

In the hydrodynamical description of heavy-ion collisions, the elliptic flow $\mathrm{v_{2}}$ and triangular flow $\mathrm{v_{3}}$ are sensitive to the quadrupole deformation $\mathrm{\beta_{2}}$ of the colliding nuclei. We produce $\mathrm{v_{2}}$ and $\mathrm{v_{3}}$ ratios qualitatively and quantitatively in most-central Xe-Xe collisions at 5.44 TeV. By employing HYDJET++ model, we study the sensitivity of anisotropic flow coefficients and mean transverse momentum to the quadrupole deformation and system-size in isotopic Xe-Xe collisions. Flow observables strongly depend on the strength of nucleon-nucleon scattering occuring in even-A and odd-A nuclei. Flow for odd-A nuclei is suppressed in comparison to flow in even-A collisions. There exists a linear inter-dependence between $\mathrm{p_{T}}$ integrated anisotropic flow and nuclear deformation. Mean transverse momentum signifies the fireball temperature in body-body and tip-tip collisions. There exists a negative linear correlation of $\mathrm{\langle p_{T} \rangle}$ with collision system-size and a positive correlation with nuclear deformation. Flow measurements in high-energy, heavy-ion collisions using isotopic collision systems, offer a new precision tool to study nuclear structure physics. Observation of nuclear structure properties like nuclear deformation in a heavy-ion collision such as this would be very interesting.


I. INTRODUCTION
Atomic nuclei are characterized by several intrinsic properties, one of them being deformed shapes arising as a result of high precision hyperfine studies in the field of atomic physics, thereby marking a non-spherical nuclear charge distribution [1][2][3].Such atomic nuclei in their ground state are deformed from a well defined spherical shape.The deformation has a significant dependence on the number of protons and neutrons, particularly in the proximity of a full shell or subshell, reflecting collective motion induced by interaction between valence nucleons and shell structure [4].This nuclear deformation emerges from the short-range strong nuclear interaction between the nucleons [2,3,5,6].This lead to a plethora of research on the topics such as shape evolution, quadrupole deformation, triaxiality or shape coexistence, octupole deformation, hexadecapole deformation, etc [2,[6][7][8].
In recent decades, Relativistic Heavy Ion Collider (RHIC) at the BNL and Large Hadron Collider (LHC) at the CERN have conducted collider experiments, developing an understanding that the deformation of the colliding nuclei strongly affect the QGP observables studied in (ultra-) relativistic heavy-ion collisions [9][10][11].One of the main motivations in this context has been elliptic flow v 2 [12].In (ultra-)relativistic nucleus-nucleus collisions, elliptic flow arises as a consequence of the quadrupole asymmetry (ellipticity) of the system created upon collision, in the plane transverse to the beam direction [13].When two nucleus collide at a certain impact parameter, the overlap area has an elliptical deformation and we observe that v 2 grows sharply with the collision impact parameter.However, for most of the isotopes, even in the case of vanishing impact parameter, it is expected that ϵ 2 ̸ = 0, and therefore v 2 ̸ = 0 from nuclear structure domain becomes argumentative as most of the nuclei exhibit a non-vanishing intrinsic quadrupole moment, or to say, an ellipsoidal deformation [14].
In ultra-relativistic heavy-ion collisions, the dynamical response of the QGP to its initial spatial anisotropy produces anisotropic flow.This initial spatial anisotropy is affected by the geometric shape of the colliding nuclei, thereby, forming an intrinsic relationship between the structure of atomic nuclei and the phenomenology of heavy-ion collisions.Thus, if we match high-energy data to the data at low-energy expectations, we can assess whether our knowledge of nuclear physics across energy scales yields consistent results.Anisotropic flow (v n ) plays a major role in elucidating the final-state effect as it constrains the viscous hydrodynamic response to the eccentricity (ϵ n ) of the energy-density distribution produced in the initial stages of the collision [15][16][17][18][19][20][21][22][23][24][25].The initial-state profile for each of the colliding nuclei is characterized by the woods-Saxon distribution for the nuclear density as: where, R = R 0 A 1/3 is the nuclear radius (R 0 = 1.15fm), ρ 0 is the center of the nucleus density, a is the surface thickness (i.e., the effective diffusivity), R ′ is the nuclear radius, Y 0 n are the spherical harmonics corresponding to each harmonic order, and R ′ (θ, ϕ) is the arXiv:2312.05853v1[hep-ph] 10 Dec 2023 nuclear surface with relevant axial symmetric quadrupole (β 2 ), octupole (β 3 ) and hexadecapole (β 4 ) deformations [1,4,[26][27][28][29].In this study, we have ignored all higher moment orders of nuclear deformation and only focussed on quadrupole moment.The other parameters such as R and a are fixed in case of every nuclei.Positive β 2 describes the overall quadrupole deformation.Nuclear deformation is a fundamental property of the atomic nucleus.It explains the correlated behaviour of the associated dynamics of nucleons.Many atomic nuclei have quadrupole or octupole deformation, which may influence the anisotropic flow coefficients ′ magnitude, its fluctuations, and correlations.Recent measurements presented detailed comparisons between Au+Au and U+U collisions [30] as well as Pb+Pb and Xe+Xe collisions [31,33].The studied observables were found compatible with nuclear deformation.However, the analysis performed so far on flow harmonics is not so successful in providing constraints useful for the detailed characterization of the deformation in the colliding nuclei.
An important aspect of the above study is that the overprediction or underestimation of experimental anisotropic flow data from theoretical calculations in larger collision systems signifies the missing physics.For instance, significant enhancement was observed in v 2 {2} in deformed central 129 Xe-129 Xe collisions compared to that in case of spherical 129 Xe nucleus. 129Xe nucleus is proposed to have a small quadrupole deformation, and the enhancement in flow can be seen in the articles [10,[32][33][34][35].In another study, 238 U-238 U collisions with the nuclei having a substantial quadrupole deformation, also were observed to have an increase in the elliptic flow in ultra-central events [30,[36][37][38].Thus, flow studies are heplful in mapping out the quadrupole deformations evoluting in the series of stable 54 Xe isotopes.Now, of course, studying xenon isotopic chain over another nuclei such as samarium isotopes can be argumentative.The reason is atrributed to the fact that, firstly, in 2017, LHC carried out an eight hour run of deformed Xe-Xe collisions at 5.44 TeV.The mass number of Xe nucleus lies in mid-between p and Pb 208 .Thus, collisions of Xe 129 nuclei would bridge the multiplicity gap between the larger Pb-ion systems and smaller systems like p+p and p+Pb.Back then, Xe was considered somewhat prolate.This deformation in xenon allows us to probe a different initial condition.However, further research confirmed Xe 129 to be triaxial.So, an odd isotope of intermediate-mass nuclei has already made a significant research impact in the field of heavy-ion physics.Secondly, Xe 129 is triaxial while Xe 128 is prolate.Differing in shape, xenon with A (mass) changing to even and odd, it becomes obligatory to study a series of even and odd isotopes and be investigated in heavy-ion collisions.Also, β 2 values for odd-and even-A xenon nuclei showed a lot of variation in low energy structure models.Thus, one is motivated to cover a range of both the kinds of nuclei with some significant value of quadrupole moment.Thirdly, in case of a recent study on samarium (Sm) iso-topes, it was seen that as the mass number increased, the quadrupole moment increased [39].However, in case of xenon, the situation was quite the opposite for the range of isotopes chosen for study, thereby making this set of isotopic chain important to be studied.Fourthly, a huge amount of research is going on to study isobaric collisions at RHIC and LHC energy regime where Ru-Ru and Zr-Zr collisions from STAR experiment have become quite interesting area of research.A similar study to perform at LHC energies could be with Nd-Sm pair as suggested in article [39].However, these are even-A nuclei.But, odd-A Xe 129 has a strong background of research in heavy-ion collisions.So, why not pick the odd-A, Xe-Sm pair.We preserve this for future work, but as the full isotopic chain of eight isotopes of Sm have been studied in heavy-ion collisions, it makes Xe worthy enough to be also studied through heavy-ion collisions.The hydrodynamic response is essentially expected to be unchanged over the isotopic chain, however, these collision systems offer a strong lever-arm to probe in detail how the initial condition of QGP responds to varying nuclear shapes.Such a study in high-energy heavy-ion collisions would give new experimental insights onto the role and effect of quadrupole deformations in such nuclei.
Most-central collisions offer the best scenario to test the sensitivity of nuclear deformation [40,41].This is because at zero impact parameter the anisotropy happens not only due to fluctuations but also due to the non-spherical shape, while in other collision centralities the impact region is dominated by the almond shape of the overlap region [42].However, in most-central collisions, the shape contribution is quite smaller in comparison to the contribution from fluctuations.Moreover, other possibilities causing fluctuations, which scale with the collision system size, affect minimally in mostcentral collisions where the overlap region upon collision is largest.Also, the large amount of energy deposited in these most-central collisions form a hot and dense quarkgluon plasma (QGP) [43] in the overlap region, whose shape and size are strongly correlated with nuclear deformation.
In this paper, we aim to address the question of whether the values of β found in low-energy literature are consistent with v 2 and v 3 data at high energy.We investigate this using a simple Monte Carlo HYDJET++ model.Further, we also aim to demonstrate that flow measurements in ultra-relativistic heavy-ion collisions using isotopic Xe systems, offer a new precision tool to study nuclear structure physics.We apply the idea to run isotopes 128−135 54 Xe, colliding symmetrically and make predictions on ratios of v n from the known values of quadrupole deformation β 2 from nuclear structure measurements, as given in table I. Here, we have ignored higher harmonic orders of nuclear deformations.Further, we will also study the transverse momentum as a function of the nuclear deformation.We also aim to study the observables in two extreme geometries of deformed collision systems.In the next section II, we discuss the essentials of Monte Carlo HYDJET++ model followed by anisotropic flow presented in the section II A. Further, we will present our results and discussions in section III followed by the summary of our work in section IV.

II. HYDJET++ MODEL
To understand the hydrodynamic response to nuclear deformations and make predictions, we employ a Monte Carlo HYDJET++ (hydrodynamics plus jets) event generator [44] developed to study heavy-ion collisions at RHIC as well as LHC energies.It performs simulation by superimposing soft (hydro-type) state and the hard state (resulting from multi-parton fragmentation) simultaneously, treating both independently.It meticulously treats soft hadroproduction (collective flow phenomenon and the resonance decays) as well as hard parton production, also examining the known medium effects (jet quenching and nuclear shadowing).The in-depth details of the model and the procedure of simulation can be found in the corresponding articles [44,45] and the references there within.The model parameters have been tuned to reproduce the experimental LHC data on various physical observables measured in Xe-Xe collisions at 5.44 TeV of center-of-mass energy per nucleon.In brief, the model is described as follows:  [46], finite-range droplet macroscopic (FRDM) and the folded-Yukawa single-particle microscopic nuclear-structure models mentioned in article [49], and some calculated quadrupole deformation parameters from [47,48].
The hard state of an event in HYDJET++ is treated using Pythia Quenched (PYQUEN) model [50].PYQUEN model modifies a jet event produced by PYTHIA by producing nucleonic collision vertices in accordance with the Glauber model at a certain impact parameter.This is followed by a rescattering-byrescattering simulation of the parton path in the dense zone and the associated radiative and collision energy losses [51][52][53][54][55]. Then the final hadronization is carried out according to the Lund String Model [56] for hard partons and in medium emitted gluons.An impact parameterdependent parameterization obtained under the frame-work of Glauber-Gribov theory is used to incorporate the medium effect called nuclear shadowing [57,58].The hard state in HYDJET++ is separated from the soft state by a free parameter named p min T .Events for which generated total transverse momentum is greater than p min T are considered as the hard part whereas events for which p min T are considered as the soft part.The soft state of an event in HYDJET++ is a thermal hadronic state generated on the chemical and thermal freeze-out hypersurfaces obtained from a parameterization of relativistic hydrodynamics with preset freeze-out conditions [59,60].Here, it is assumed that the hadronic matter produced in a nuclear collision reaches a local equilibrium after a short period of time (<1fm/c) and then expands hydrodynamically.Basically, the hadronic matter created in heavy-ion collisions is a hydrodynamically expanding fireball with the equation of state of an ideal hadron gas.The total yield of particular particle species is determined by the freeze-out temperature, chemical potential and by the total co-moving volume or effective volume of particle production which is a functional of the field of collective velocities on the hypersurface.The effective volume utilizes the collective velocity profile and the form of hypersurface canceling out in all particle number ratios.This effective volume can be used to calculate the hadronic composition at both chemical and thermal freeze-outs [61].At the former one, which happens soon after hadronization, the chemically equilibrated hadronic composition is assumed to be established and frozen in further evolution.This implies that the chemical concentrations of hadron gas do not change during the evolution.This is justified as the rate of expansion in hadron gas is larger than the rate of inelastic reactions and lesser than that of elastic collisions [62,63].In HYDJET++, the presumption of single freeze-out is eliminated as the particle densities at the chemical freeze-out stage are too high to consider the particles as free streaming [61].Therefore, a more complex scenario of different thermal and chemical freezeouts(T ch ≥ T th ) is opted.The system expands hydrodynamically with frozen chemical composition in between these two freeze-outs, then cools down and the hadrons stream freely as soon as the thermal freeze-out temperature is reached.In non-central collisions, the shape of the emission region in the transverse (x-y) plane is approximately elliptical with the (z-x) plane coinciding with the reaction plane.The two radii of the ellipse are: ) and the scale factor R s = (R 2 x + R 2 y )/2.Then, the transverse radius R(b, ϕ) of the fireball on the azimuthal angle is: where, This means that the dimensionless ratio R s (b)/R s (0) at the freeze-out moment depends on the collision energy, the radius R A of the colliding (identical) nuclei and the impact parameter b through a dimensionless ϵ s (b) only.

A. Anisotropic flow in HYDJET++
Quark-Gluon Plasma, the deconfined state of color charges is believed to be created in relativistic heavyion collisions [64].The pressure gradients generated in the QGP medium via most probably non-central collisions at high relativistic energies convert the initial anisotropies to momentum anisotropies of the produced particles via multiple interactions.This phenomenon is called anisotropic flow.This anisotropic flow is defined by the coefficients from the fourier expansion of the azimuthal distribution of the produced particles [65,66] as- where, ψ= azimuthal angle of the produced particle, n = harmonic value, and For n=2, we have the second order coefficient called as elliptic flow, which is sensitive to the early evolution of the system.In HYDJET++ model, the reaction plane for even order flow is zero for each event.The elliptic flow coefficient v 2 is determined in the hadron distribution over the azimuthal angle ψ relative to the reaction plane ψ R , such that, The above equation in the simplified form for v 2 in terms of particle momenta is expressed by- The soft particle emission which contributes most to the elliptic flow, arises from the freeze-out hypersurface.To simulate higher azimuthal anisotropy harmonics, parameterization of the fluid velocity corresponding to each harmonic on the freeze-out surface is performed in HY-DJET++.The model does not describe the evolution of the fireball from the initial to the final freeze-out stage.It utilizes simple parameterization of the freeze-out hypersurface.The anisotropic elliptic shape of the initial overlap of the colliding nuclei results in a corresponding anisotropy of the outgoing momentum distribution.The second harmonic v 2 is described using the coefficients ϵ 2 (b) and δ 2 (b).These are the second order spatial anisotropy and momentum anisotropy parameters, respectively.ϵ 2 (b) exemplifies the elliptic modulation of the final freeze-out hypersurface at a given impact parameter b, whereas δ 2 (b) deals with the alteration of flow velocity profile.These two parameters can be treated independently for each centrality or can be made interdependent via the dependence on the initial ellipticity ϵ 0 (b) = b/2R A where R A is the nucleus radius.Here, we are treating them independent of each other.
The spatial anisotropy gets transformed into momentum anisotropy at freeze-out, because each of the fluid cells is carrying some momentum.In case where δ 2 (b) ̸ = 0 even the spherically symmetric source can mirror the spatially contracted one.The elliptic flow coefficient v 2 (ϵ 2 , δ 2 ) in the hydrodynamical approach (incorporated in HYDJET++ framework) [44] is given as: For triangular flow v 3 in HYDJET++, the model has another parameter ϵ 3 (b), for spatial triangularity of the fireball.Thus, the modified radius of the freeze-out hypersurface in azimuthal plane reads: where, ϕ = spatial azimuthal angle of the fluid element relative to the direction of the impact parameter.The phase ψ RP   3   gives us the advantage to introduce a third harmonic having its own reaction plane, distributed randomly with respect to the direction of the impact parameter (ψ RP 2 = 0).Such modifications do not affect the elliptic flow (controlled by ϵ 2 (b) and δ 2 (b)).Hence, the triangular dynamical anisotropy can be incorporated by the parameterization of the maximal transverse flow rapidity [67], Hence, we can calculate higher harmonics with respect to the direction of the impact parameter b.
HYDJET++ model provides reasonable description of anisotropic flow in heavy-ion (small, intermediate, and large) collision systems at both RHIC and LHC energies [9][10][11].The results show proper qualitative as well as quantitative explanation for the flow in spherical and deformed collision systems.We simulate 10 5 minimum bias, symmetric, and isotopic (0-5)% most-central 128−135 54 Xe collisions under HYDJET++ framework with R 0 = 1.15A 1/3 fm and surface thickness (i.e., the effective diffusivity) a = 0.59 ± 0.07fm.The observables are calculated under |η| < 0.8 and 0 < p T < 5.0 GeV/c kinematic ranges.Under these circumstances, the number of participants N part and the number of binary collisions N coll as a function of atomic mass have been shown in figure 1, thereby certifying our model.It is observed that N part and N coll linearly increase with atomic mass.This is because, HYDJET++ model uses Glauber-type initial condition for soft (hydro-type) part.This atomic mass, if we see table I, is almost inversely proportional to nuclear deformation.Spherical collisions have higher N part and N coll compared to deformed collision systems.The difference in the number of participants of deformed and spherical isotopic collisions is significantly visible at lower atomic mass but eventually decreases as one moves to higher A. The model parameters which control the execution of our Monte Carlo generator show suitable match with experimental data at LHC as presented in the recent articles [10,11] These input model parameters which control the execution of the Monte Carlo generator are obtained by simulating (0-5)% most-central 129 Xe-129 Xe collisions and then matching it with the experimental data at LHC.

III. RESULTS AND DISCUSSION
The upper panel of figure 2 shows elliptic flow as a function of transverse momentum (p T ) for different isotopic minimum bias Xe-Xe collisions with different prolate-type deformations.The lower panel presents triangular flow for the same scenario.Here, the bin values have been scaled by some factor shown in the figure, to visualize the difference between the flow of different isotopic collisions.Using the same factors, the scaling in the other results also have been introduced.
It is observed that, for even isotopes, anisotropic flow v n (n=2,3) decreases as atomic mass increases and quadrupole deformation decreases.However, for odd isotopes, situation appears to be quite opposite; anisotropic flow increases as atomic mass increases and quadrupole deformation decreases.Thus, for even isotopes, anisotropic flow is directly proportional to nuclear deformation whereas for odd isotopes, this flow is inversely proportional to nuclear quadrupole moment.Thus, for an odd-A isotopic chain, the relationship between nuclear deformation and anisotropic flow (as in magnitude) is opposite to that of the even-A isotopic chain.This hints us at the strength of nucleon-nucleon scattering taking place in the two types of isotopic Xe-Xe collisions.The flow for odd-A nuclei is suppressed compared to even-A nuclei.Thus, v n -p T correlation here directly connects us to the initial collision geometry and thereby back to the shape of the colliding nuclei.Elliptic flow v 2 is more than triangular flow, the difference being very less.Flow increases as we move from low p T to high p T region.The peak value at which flow is maximum is around 3.0 GeV/c for elliptic flow.It shifts towards higher p T as the harmonic value, n increases.In figure 3, we present ratios of v n (β 2 , β 3 , ...) for given values of β 2 (with rest orders as zero) to that for spherical nuclei.In this scenario, three types of ratios are studied namely: (i) v 2 (β 2 , 0)/v 2 (0, 0), (ii) v 3 (β 2 , 0)/v 2 (0, 0), and (iii) v 3 (β 2 , 0)/v 3 (0, 0).Here, we observe that the ratios of elliptic flow and triangular flow of deformed collisions with respect to their respective spherical collisions, show similar qualitative behaviour whereas the ratio of deformed v 3 to spherical v 2 has a different inference.Ratios (i) and (iii) present a somewhat linear relationship with increasing transverse momentum and gradually decrease towards higher p T values.Higher quadrupole moment gives rise to higher flow for even-A isotopic collisions, hence, higher ratios.In odd-A isotopic collisions, the situation is reversed, higher deformation results in lower anisotropic flow, eventually giving smaller ratios.On the other hand, ratio (ii) of deformed v 3 to spherical v 2 presents a gaussian-type of behaviour.It increases at low p T and reaches its zenith at around 3.0 GeV/c, then decreases (very slowly) towards higher transverse momentum.Thus, the figures bring into picture the effect of deformation and collision system-size for the various isotopic Xe-Xe collisions.The second ratio v 3 (β 2 , 0)/v 2 (0, 0) presents a significant viscosity dependence.Viscous damping is smaller in v 2 than in v 3 .Therefore, we observe that the ratio increases.Such an inspection with nuclear quadrupole moment suggests that for even-A nuclei, viscous damping increases with increasing nuclear quadrupole moment; while for odd-A nuclei, viscous damping decreases as nuclear deformation increases.Another important aspect is that the ratio is an important parameter for the v 2 − v 3 puzzle.So, it is necessary to measure it as a function of nuclear deformation.In order to visualize the effect of nuclear deformation separately, we normalize our results with the total number of nucleon-pairs in each collision system (please see figure 4) and find the quantitative role of deformation in the three kinds of isotopic Xe-Xe collision systems, where we can easily ascertain the above inferences.Elliptic flow is larger than triangular flow is very evident from the above figures.The quantitative effect of quadrupole deformation irrespective of the atomic mass show the role of quadrupole moment on the initial nuclear density profile of nucleus, and hence on the initial geometry of the collision overlap region.
To study the system-size dependence of anisotropic flow, figure 5, presents the ratio of flow observables in (deformed-) Xe-Xe collisions at 5.44 TeV to spherical Pb-Pb collisions at 5.02 TeV.A baseline expectation is that the flow observables should be equal in the larger (Pb) and smaller (Xe) collision systems.However, several effects break this scaling, and can be taken as corrections to the scale-invariant baseline.Here, we present this test for the various isotopic collisions as a function of transverse momentum p T .We observe that the flow strongly increases in central collisions for the Xe system compared to the Pb system, which distinctly indicates about the nuclear deformation in Xenon.Based upon the strength of nuclear deformation and isotopic collision system-size, the ratio of each flow harmonic corresponding to each isotope exists quantitatively.The scaling shows that the effect of nuclear deformation is visible maximum in low p T region and decreases as we move towards higher p T values.An important fact to be noted here is that elliptic flow v 2 arises not only from fluctu- ations in the initial stages of a collision but also from the initial-state geometry of the nuclear overlap region.In central collisions, fluctuations become quite important while the geometry plays its dominant role in peripheral collisions.Elliptic flow fluctuations due to the smaller size and non-spherical nuclear shape give quantitative increment which evidently can be seen from our results of isotopic collisions.3 ) Xe /(v 2 3 ) Pb as a function of (β 2 2 ) Xe where we find a strict inter-dependence.The ratio for elliptic flow varies linearly with β 2 2 whereas for triangular flow, the ratio appears to be quite independent of β 2 2 , revolving around a constant value as quadrupole moment increases and atomic mass decreases.HYDJET++ model does not include nuclear density profiles for nuclei with some intrinsic deformation.Therefore, the nuclear density function needs to be modified for deformed nucleus like U, Xe, etc.In our present work, we have used the Modified Woods-Saxon (MWS) given by equation 1 to calculate the initial distributions of partons, etc. for deformed nuclei [10,11,35].Within HYDJET++ framework, ρ 0 = ρ const 0 + correction, The correction term is calculated as=ρ const 0 where f = 0.54fm, The values of different parameters have been taken from the reference [49].However, HYDJET++ model works in cylindrical coordinate system unlike the AMPT model that works in spherical coordinate sysytem.Thus, the modified nuclear density profile function in equation 1 is converted from spherical polar coordinates (r, θ, ϕ) to cylindrical polar coordinates (ρ, z, ψ) using the relations θ = tan −1 (z/r) and θ = tan −1 (r/z) for body-body and tip-tip geometrical configurations, respectively.Here, r in the relations is basically ρ of cylindrical polar coordinate system and not spherical polar coordinate r.This representation helps in eliminating the confusion of ρ in cylindrical polar coordinate system with the ρ in the nuclear density function.The Woods Saxon nuclear density profile in tip-type and body-type geometrical configurations is defined in the way two nuclei collide when approaching towards each other along the direction of propagation as shown in figure 7.In figure 8, we present the p T -integrated anisotropic flow in deformed Xe-Xe collisions for the various isotopes at 5.44 TeV centre-of-mass energy.It is observed that, flow for body-body collisions is larger than flow for tip-tip collisions.Minimum bias collisions lie almost in mid-between these two extreme geometries.As nuclear quadrupole moment β 2 increases or atomic mass decreases, flow increases.As a result, a linear interdependence is observed between the observables.The intersection point of the various types of collisions suggests that the geometrical configurations overlap or have same flow.At this point, it is difficult to disentangle the collisions on the basis of nuclear orientation just before collision.The intersection point of the various types of collisions is observed to be at a lower value of β 2 for elliptic flow (intercept around 0.0333-0.0334)than for triangular flow (intercept around 0.0204-0.0206).
In order to study a proper shape-size correlation, we must also look into mean transverse momentum.⟨p T ⟩ carries essential information about the high-density deconfined state of strongly interacting matter, quark-gluon plasma.The slope of p T distribution measures the inverse of the source temperature of the fireball, increases as one moves from most-central to most-peripheral collisions.This testifies the fact that the fireball temperature is maximum in most-central collisions [10,11].However, in deformed collisions it is quite exciting how does the deformation affects this fireball temperature.In figure 9, we study the mean transverse momentum as a function of deformation parameter.Here, we observe that, with increasing deformation, mean transverse momentum increases, but not so strongly.In article [69], it is explained that in ultracentral collisions, the temperature increases as a function of the multiplicity, which in turn implies a rise of the mean transverse momentum of charged hadrons, ⟨p T ⟩, observed in the final state, due to tight correlation with the temperature.As a result, the source temperature increases as nuclear deformation increases.Viewing from the door of geometry, we see that ⟨p T ⟩ is higher for tip-tip collisions than body-body collisions.This implies that fireball temperature is higher in tip-tip collisions than in body-body collisions.Expressing in terms of system-size, we can say that, as systemsize (atomic mass A) increases, mean ⟨p T ⟩ decreases.This is more evidently visualized figure 10 where the ratio of mean transverse momentum in (deformed-) Xe-Xe collisions at 5.44 TeV to mean transverse momentum in spherical Pb-Pb collisions at 5.02 TeV is presented.It is observed that the ratio increases with increasing deformation and decreasing system-size.The mean transverse momentum in Xe-Xe collisions decreases by (1-2)% compared to larger spherical Pb-Pb collisions at 5.02 TeV.Another important point to be noted is that the increasing ⟨p T ⟩ with nuclear deformation is attributed to the increasing transverse radial flow with increasing deformation; and as system-size in terms of atomic mass increases, this transverse radial flow decreases.

IV. SUMMARY AND OUTLOOK
In summary, we performed a study on some isotopes of xenon using Monte-carlo technique under HYD-JET++ framework.Here, we studied two observables, anisotropic flow and transverse momentum, as a function of the two important properties of an isotope, namely the nuclear (quadrupole) deformation and system-size as in atomic mass, (A).Anisotropic flow for odd-A isotopes is suppressed than flow for even-A isotopes.As a function of transverse momentum, for even isotopes, anisotropic flow is directly proportional to nuclear deformation whereas for odd isotopes, this flow is inversely proportional to nuclear quadrupole moment.The ratios of elliptic flow and triangular flow of deformed collisions with respect to their respective spherical collisions present a linear relationship with increasing transverse momentum and gradually decreasing towards higher p T region.Large quadrupole moment gives rise to higher flow for even-A isotopic collisions, hence, higher ratios.Opposite happens for odd-A isotopic collisions.The effect of nuclear deformation is visible maximum in low p T region and decreases moving towards higher p T values.Xenon, being deformed gave us an opportunity to present our study in two extreme geometries of collisions, body-body and tip-tip type of Xe-Xe collisions at 5.44 TeV.We found a linear inter-dependence between p T integrated anisotropic flow and nuclear deforma-tion.Mean transverse momentum marks the fireball temperature in body-body and tip-tip collisions.⟨p T ⟩ is higher in tip-tip collisions than body-body collisions.There exists a negative linear correlation of ⟨p T ⟩ with collision system-size appears and a positive correlation with nuclear deformation.
The modeling of heavy-ion collisions helps in pulling up the hydrodynamic response of v n to the contribution of nuclear deformation and thereby to the eccentricity, which is quite different from the case of spherical nuclei collision systems.Thus, we see that heavy-ion isotopic collisions offer immense potential to constrain the nuclear deformation parameters and help in testing the state-of-the-art nuclear structure models.Our measurements of v 2 and v 3 show strong differences between deformed isotopic 128−135 54 Xe-128−135 54 Xe collisions and spherical 208 Pb-208 Pb collisions.These differences are well explained from the quadrupole deformation β 2 in 128−135 54 Xe isotopes.However, this deformation parameter can be a more useful tool to constrain the initial-state geometry as well as the shape of the nucleus, for example the study of initial state fluctuations in anisotropic flow might curb the geometry of xenon nucleus being either prolate or oblate or purely triaxial such as in [68].This we preserve for future work.Hence, isotopic heavy-ion collisions act as a precision tool to describe the role of nuclear quadrupole deformation in the field of heavy-ion physics and will provide strong support to its role which is quite debative in nuclear structure physics.Such analysis may be hopefully done at the relativistic high-energy collider experiments at LHC in future but for that individual study of either isotope in the field heavy-ion collisions is necessary as a background study.

FIG. 1 .
FIG. 1. Variation of number participants Npart and number of binary collisions N coll with atomic weight (A) in mostcentral (0-5)% Xe-Xe collisions at various values of nuclear deformation parameter (β2).

FIG. 2 .
FIG. 2. Transverse momentum dependence of anisotropic flow vn (n=2,3) in deformed Xe-Xe collisions at 5.44 TeV center of mass energy using HYDJET++ model.The figure shows flow for eight isotopes (A=128-135) of Xenon with different values of quadrupole moment.

FIG. 3 .
FIG. 3. Transverse momentum dependence of the ratios of anisotropic flow vn (n=2,3) from deformed Xe-Xe collisions to their respective spherical Xe-Xe collisions at 5.44 TeV center of mass energy using HYDJET++ model.The figure shows flow for eight isotopes (A=128-135) of Xenon with different values of quadrupole moment.

FIG. 4 .
FIG. 4. Scaling of the ratios of anisotropic flow vn (n=2,3) from deformed Xe-Xe collisions to their respective spherical Xe-Xe collisions at 5.44 TeV center of mass energy using HYDJET++ model by 2A.The figure shows flow for eight isotopes (A=128-135) of Xenon with different values of quadrupole moment.

FIG. 8 .
FIG. 8. Variation of pT-integrated eliptic and triangular flow with respect to nuclear quadrupole moment β2 in deformed Xe-Xe collisions at 5.44 TeV center of mass energy using HY-DJET++ model.The results have been presented in bodybody and tip-tip type of geometric nucleus-nucleus collisions along with the minimum-bias results.

FIG. 9 .
FIG. 9. Variation of mean transverse momentum ⟨pT⟩ as a function of nuclear quadrupole deformation parameter in deformed Xe-Xe collisions at 5.44 TeV center of mass energy using HYDJET++ model.The figure shows flow for isotopes (A=128-135) of Xenon with different values of quadrupole moment.The results have been presented in body-body and tip-tip type of geometric nucleus-nucleus collisions along with the minimum-bias results.

TABLE II .
. The model input parameters for 129 Xe-129 Xe collisions at 5.44 TeV are provided in table II.Model parameters for 129 Xe-129 Xe collisions at 5.44 TeV in most central class of collisions.