Probe nuclear structure using the anisotropic flow at the Large Hadron Collider

Recent studies have shown that the shape and radial profile of the colliding nuclei have strong influences on the initial condition of the heavy ion collisions and the subsequent development of the anisotropic flow. Using A Multi-Phase Transport model (AMPT) model, we investigated the impact of nuclear quadrupole deformation β2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _2$$\end{document} and nuclear diffuseness a0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0$$\end{document} of 129\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{129}$$\end{document}Xe on various of flow observables in Xe–Xe collisions at sNN=\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}}}=$$\end{document} 5.44 TeV. We found that β2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _2$$\end{document} has a strong influence on central collisions while a0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0$$\end{document} mostly influences the mid-central collisions. The relative change of flow observables induced by a change in β2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _2$$\end{document} and a0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0$$\end{document} are also found to be insensitive to the values of parameters controlling the strength of the interaction among final state particles. Our study demonstrates the potential for constraining the initial condition of heavy ion collisions using future system scans at the LHC.


Introduction
Ultra-relativistic heavy-ion collisions conducted at both the Relativistic Heavy-Ion Collider (RHIC) and the Large Hadron Collider (LHC) provide a pivotal platform for the comprehensive study of Quantum Chromodynamics (QCD) across both perturbative and nonperturbative regimes.These collisions afford the remarkable opportunity to recreate a novel state of matter, quark-gluon plasma (QGP), characterized by extreme temperatures and densities in the early stages of high-energy heavy-ion interactions.Over the past two decades, a dedicated endeavor has been directed towards extracting precise insights into the properties and the dynamic evolution of QGP [1][2][3][4][5][6][7][8][9][10][11].Anisotropic flow, which quantifies the anisotropic expansion of the produced particles, has been a powerful tool for QGP studies [5][6][7][8][9][10][11].It is characterized by the Fourier coefficients v n of the azimuthal particle distribution [12]: where φ is the azimuthal angle of the final particles, Ψ n is the n th -order flow symmetry plane, and v n is called flow coefficient.From Eq(1), the flow coefficient v n can be a e-mail: you.zhou@cern.charXiv:2309.09663v1[nucl-th] 18 Sep 2023 defined as: Here, the angular bracket ⟨⟩ denotes the average over all particles in one event.The v n and flow angle Ψ n are the magnitude (amplitude) and angle (orientation) of the flow vector, defined as: Systematic measurements on v n , event-by-event fluctuations of v n , and correlations of different flow coefficients v n , v m , v k have been previously reported in Refs.[13][14][15][16][17].
By performing the Bayesian fit on the extensive flow data, critical information on the temperature dependence of shear and bulk viscosity over entropy density ratios of the QGP, η/s(T ) and ζ/s(T ), can be extracted [18][19][20][21].
In addition, the flow measurements give direct access to the event-averaged initialstate shape of the nuclear overlap region and their event-by-event fluctuations.For typical heavy-ion collisions, the nuclear density profile in the initial-state can be described by Woods-Saxon distribution: where a 0 denotes the nuclear diffuseness, while R 0 represents the half-width radius.The nuclear surface, denoted as R(θ, ϕ), is expanded in terms of spherical harmonics Y n,m , where we retain terms up to n = 4 as expressed in the Eq. ( 4).Furthermore, β 2 , β 3 , and β 4 stand for the quadrupole, octupole, and hexadecapole deformation parameters, respectively.The parameter γ characterizes the triaxial shape, depicting any imbalance present within the axes of the spheroid.Analogous to γ, α 3,m and α 4,m describe the inequality of axes and satisfy the normalization condition.
In recent years, various observables have been investigated for their sensitivities to nuclear structure parameters in heavy-ion collisions [22][23][24][25][26][27][28][29].As mentioned above, the anisotropic flow reflects the initial spatial anisotropies in the overlap region of the colliding nucleus.Thus, it serves as an ideal probe of initial conditions and can be utilized for the nuclear structure study.The flow coefficient v n has been found to be sensitive to the deformation, characterized by deformation parameters β n , in 96 Ru-96 Ru, 96 Zr-96 Zr, 238 U-238 U, 197 Au-197 Au collisions [22][23][24][25][26]. Beyond v n , the spotlight extends to multi-particle cumulants of v n and nonlinear flow, underlining their potential for discerning the parameters β n , a 0 , and R 0 [25][26][27].Moreover, the mean transverse momentum of the produced charged hadrons, denoted as [p T ], which reflects the initial overlap region's size, is a valuable probe for exploring neutron skin thickness and nuclear deformation [28].The fluctuations of v n and [p T ], as well as the correlations between them, quantified by Pearson correlation coefficient(PCC) and denoted as ρ(v 2 n , [p T ]), emerge as a nuanced avenue for constraining deformation parameters β n and the triaxial parameter γ [29].
Almost all these studies are conducted within the RHIC energies (GeV energy scale), while a similar study at the LHC energies ranges (TeV energy scale) is still lacking at the moment.In particular, 129 Xe is the nucleus believed to have quadrupole and triaxial deformation [30], determined from low-energy nuclear theory and experiments [31][32][33][34][35][36].Nevertheless, only very few selected flow observables, such as v n {2} and ρ(v 2 n , [p T ]), have been studied so far [37,38].The realm of systematic investigations targeting more intricate flow observables, based on multi-particle correlations

A Multi-Phase Transport Model
A Multi-Phase Transport (AMPT) Model [39] holds widespread application in the realm of ultra-relativistic nuclear collisions for investigating the initial conditions and transport characteristics of Quark-Gluon Plasma (QGP) [39][40][41][42][43][44][45][46][47].This paper employs the AMPT model incorporating the string melting scenario.The model encompasses a sequence of processes, including the initial conditions of parton production, interactions among partons, hadronization through coalescence, and, finally, hadronic rescattering.More specifically, the nucleons are generated using the HIJING model [48] to establish their spatial and momentum distributions, after which they convert into partons.The interactions among these partons are governed by Zhang's Parton Cascade model (ZPC) [49].In this model, the parton-scattering cross-section σ can be characterized by the following: where α s is the QCD coupling constant and µ is the screening mass.This crosssection delineates the dynamic expansion of the QGP phase.Subsequent to the parton cascade, partons combine to form hadrons, termed hadronization, employing a coalescence model [50].Following hadronization, the interactions among the resulting hadrons in the final state are elucidated by the ART model [51].
The AMPT simulations for the Xe-Xe collisions are employed by parameterizing the nucleon density profile with Woods-Saxon distribution shown in Eq (4).To investigate the effect of nuclear deformation and diffuseness, and also the sensitivity to the system's dynamic evolution, several sets of values for β 2 , γ, a 0 , and σ are used for comparisons.Here we use sets 1-4 to study the effect of nuclear deformations, i.e., by changing β 2 from 0 (spherical) to 0.18 (deformed, obtained from Ref. [30]) and changing γ from 0 (prolate) to 27 • (triaxial) to 60 • (oblate).Also, comparing the results from sets 2 (a 0 = 0.57 from Ref. [30]) and 5 (a 0 = 0.492 from Ref. [52]) can provide information for nuclear diffuseness.For the impact of transport properties, we use σ = 3.0 mb (set 6) [53,54] instead of σ = 6.0 mb (set 3) [55].The observables are shown in centrality dependence, where the centrality in this study is determined by the impact parameter in Xe-Xe collisions.More detailed information concerning the input parameters of the AMPT model can be found in table 1: 3 Analysis details

Observables
Experimentally, flow coefficients cannot be obtained directly via Eq.( 2) but from two-and multi-particle correlations/cumulants [14,[56][57][58]: where c n {2} is the two-particle cumulant.In this analysis, v 2 {2},v 3 {2},v 4 {2} are studied.In addition, the four-particle cumulants of v n , denoted as v n {4}, can be obtained via four-particle cumulants c n {4}: The higher order cumulants of v n are defined in a similar way, noted as v n {6}, v n {8}, etc.The two-and multi-particle cumulants of v n have different contributions from flow fluctuations σ vn .It is well known that for the Gaussian type flow fluctuations and in the case σ vn ≪ vn , we have [59]: here σ vn is the standard deviation of v n distribution, which represents the event-byevent fluctuations of v n .Then mean flow coefficient vn (which is also known as v n from the flow symmetry plane) and the flow fluctuation σ vn can be extracted from the combination of v n {2} and v n {4} according to the Eq. ( 8): For central and semi-central collisions, the lower order flow coefficients v n (for n = 2, 3) are linearly correlated with the initial eccentricity coefficients ε n [13,60].While higher harmonic flow v n (for n > 3) not only has the linear response to the corresponding initial ε n but also has contributions from the lower order ε 2 and/or ε 3 [61][62][63].The latter is called the nonlinear flow mode.For example, V 4 and V 5 can be decomposed into the linear and nonlinear components: Here V NL n and V L n are the nonlinear and linear (or called leftover) components, respectively.Their magnitudes are denoted as v n,mk and v L n .Besides, χ n,mk is the nonlinear coefficient representing the strength of nonlinear response from lower order eccentricities [63].The correlation between different order flow symmetry planes can be studied by calculating the ratio between v n,mk and v n {2} [64]: ρ 4,22 and ρ 5,32 can be used to study the correlations between Ψ 2 and Ψ 4 , as well as the correlations between three planes of Ψ 2 , Ψ 3 and Ψ 5 .The study of nonlinear flow modes, i.e., v n,mk , ρ n,mk , χ n,mk have been performed before, they could provide further constraints on the initial conditions [14,61,65].
The correlations between v 2 n and v 2 m can be quantified via normalized symmetric cumulants NSC(m, n), defined as [14]: where the angular bracket ⟨⟩ represents an average over all events.It allows to study if v 2 n and v 2 m are correlated, anti-correlated or uncorrelated, if NSC(m, n) > 0, < 0 and = 0, respectively.In this paper, NSC (3,2) and NSC(4, 2) will be studied with the AMPT model to see if the results could bring extra information into the initial conditions and the structure of 129 Xe.

Multi-particle correlation
All the flow observables introduced in section 3.1 can be obtained via the multiparticle correlation method [14,[56][57][58].To begin with, the flow coefficient v n can be calculated using two-particle correlations: where φ 1 and φ 2 are azimuthal angles from different particles.Double brackets ⟨⟨⟩⟩ denote the average over all particles in an event and then the average over all events.For multi-particle cumulants of v n [58]: where : and The magnitude of nonlinear flow mode v n,mk could be obtained via the multi-particle correlations [63,64].For n = 4, 5, which is studied in this paper, we have: They quantify the magnitude of the nonlinear mode in high-order flow coefficients.By subtracting v 4,22 and v 5,32 from v 4 , and v 5 , respectively, we can easily calculate the magnitude of linear modes: Nonlinear coefficient χ n,mk also describes the nonlinear contributions but is independent of v 2 and v 3 .It can be derived by taking the ratio of v n,mk and corresponding lower-order v n (n = 2, 3): The ρ n,mk , correlation between different order flow symmetry planes can be obtained by taking Eqs. ( 13) and ( 17): Expanding Eq. ( 12) with multi-particle correlations, normalized symmetric cumulants NSC(m, n) can be expressed as: All these observables are now expressed in terms of 2-and multi-particle correlations and then can be calculated using the Generic Framework [14,66] or its latest implementation Generic Algorithm [58].

Study on the nuclear deformation
The centrality dependence of v 2 {2} and v 3 {2} in Xe-Xe collisions at 5.44 TeV are shown in Fig. 1.Here, v 2 {2} increases significantly in the Ultra-Central Collisions (UCC) region when changing β 2 from 0 (red diamonds) to 0.18 (the other markers).In the UCC region, where the two colliding nuclei almost fully overlap, the eccentricity of the overlapping region is determined by the shape or nuclear structure of the colliding nuclei.As it's well known, the β 2 parameter in the Woods-Saxon distribution characterizes an elliptical shape of the Woods-Saxon nucleon density profile.A nonzero β 2 of 129 Xe enhances the initial eccentricity ε 2 compared to the one with a spherical shape where β 2 = 0 and consequently leads to an increase in v 2 because of the linear correlation between ε n and v n , i.e., v n ∝ ε n , for n = 2 or 3 [67,68].In Fig. 1(b), v 3 {2} has negligible differences when changing the β 2 values, because the triangularity ε 3 of the overlapping region does not depend on β 2 [24].Figure 1 also shows consistent results of v n {2}(n = 2, 3) with variations in the triaxial parameter γ across the 0-30% centrality range.This is not a surprise, as probing the 3-D structure usually requests correlation involving more than two particles.In light of preceding investigations [22], a potential sensitivity of ε 2 to γ was indicated for centrality range 0-0.2%, as concluded from simulations using the initial-state model.However, the current study refrains from probing such phenomena, relying on the AMPT model's final state information.The intrinsic computational demands inherent to capturing the entirety of the dynamic evolution process within AMPT serve as a significant impediment.
The value of v 2 {2} receives the contributions not only from the initial eccentricity ε 2 but also from the event-by-event eccentricity fluctuations.To undertake a more comprehensive exploration of the ramifications imposed by the initial event-by-event eccentricity fluctuations and correspondingly the final state elliptic flow fluctuations, the utilization of both v 2 {2} and v 2 {4} is employed.This combined approach is particularly insightful as these two observables carry opposite contributions from flow fluctuations, as shown in Eqs.(8).As a result, the centrality dependence of v 2 {4}, as well as the v2 and σ v2 in Xe-Xe collisions at 5.44 TeV are presented in Fig. 2. As shown in Fig. 2 diamonds) to 0.18 (the other markers), despite the sizable statistical uncertainties in the presented centrality region.It was previously suggested that v 2 {4} is influenced by the nuclear diffuseness a 0 , with a weak sensitivity to β 2 using different nuclei at the RHIC energy [25].This aligns with our AMPT results shown in Figs.2(a) and 9(a).In Fig. 2(b)(c), v2 and σ v2 have significant increases in UCC region when changing β 2 from 0 to 0.18.The sensitivity of v2 to β 2 primarily arises from its linear correlation with ε 2 .In the case of flow fluctuations σ v2 , a deformed nucleus assumes varied orientations compared to a spherical shape.These random orientations result in stronger fluctuations of the elliptic flow.
In Fig. 2, v 2 {4}, v2 , and σ v2 stay unchanged with variations in γ.These observables do not lend themselves to the task of constraining the unique triaxial parameter characterizing the nucleus 129 Xe.
An evident sensitivity of v 2 {6} to β 3 is reported in a previous study [25].The same work also discusses the correlation between the sensitivity of β n and the number of particles used in the multi-particle correlations.Following the same concept, it is expected that observables like v 2 {6} and v 2 {8} may yield distinctive insights into the nuclear structure, compared to v 2 {2} and v 2 {4} discussed above.The centrality dependence of v 2 {6} and v 2 {8} are shown in Fig. 3. Amidst the sizable uncertainties, the current study refrains from drawing a firm conclusion regarding the sensitivity of these two observables to either β 2 or γ across the 0-30% centrality interval.
The v n with multi-particle correlations, i.e., v 2 {4}, v 2 {6}, and v 2 {8}, were reported to exhibit a slight difference at the precision of 1-2% in Pb-Pb collisions [69].This difference is believed to have originated from the deviations from a Bessel-Gaussian shape, particularly a non-zero skewness of the event-by-event v 2 distribution.To determine whether a similar pattern persists in Xe-Xe collisions and explore the discrepancies between higher-order multi-particle cumulants, the ratios v 2 {m}/v 2 {4} for m = 2, 6, and 8 are presented as a function of centrality in Fig. 4. The ratio of v 2 {2}/v 2 {4} increases dramatically toward the central region in Fig. 4 193 GeV [26].Subsequent investigations with increased statistics are anticipated to elucidate the presence of this sensitivity.However, the current study is constrained by the substantial computing demands involved.Concerning higher-order cumulants with m=6,8, the ratios v 2 {6}/v 2 {4} and v 2 {8}/v 2 {4} are close to unity within uncertainties and they stay unchanged when varying the nuclear structure configurations β 2 and γ. .This is consistent with previous measurements in Pb-Pb collisions at the LHC [64].The linear component v L 4 shows an absence of sensitivity to the quadrupole deformation parameter β 2 across the 0-30% centrality range.A similar observation has been reported recently in Ref. [26], where both v L 4 and v 4 {2} show no sensitivity to β 2 .In contrast, the sensitivity of both v 4,22 and v L 4 to the hexadecapole deformation parameter β 4 is seen [26], owing to their shared provenance from ε 4 .The nonlinear component v 4,22 , as illustrated in Fig. 5(b), manifests a degree of 20% reductions with zero β 2 within the UCC region, albeit with a magnitude smaller than that of v 4 {2}.Since the nonlinear flow v n,mk probes a smaller spatial distribution than the flow coefficient v n , it can impose more stringent constraints on initial conditions and nuclear structure parameters.Variations in γ parameter show that v 4,22 , v L 4 , and the v 4 {2} all exhibit diminished sensitivities to γ within the 0-30% centrality range, and thus they are unable to probe the triaxial nuclear structure.
The nonlinear coefficient χ 4,22 , defined in Eq. ( 10), is not solely affected by transport properties but also by initial conditions [70].It previously exhibits a weak dependence on nuclear deformation β 2 in central collisions at the RHIC isobar runs [26,27,71].Verifying whether the same conclusion holds in the context of Xe-Xe collisions remains pertinent.Nonlinear correlation ρ 4,22 , introduced in Eq. (11), demonstrates an insensitivity to the influences of final state interactions [70], rendering it an effective tool for probing initial conditions.In Fig. 6, the centrality dependence of χ sensitivity of v 4,22 to β 2 primarily stems from the sensitivity conveyed by (V 2 ) 2 .A similar scenario can be observed in the case of v 5,32 and χ 5,32 , as shown in Fig. 15 in the appendix, suggesting that the sensitivity arises from v 2 v 3 .Furthermore, the results in Fig. 1 have already indicated that v 3 is insensitive to β 2 , thereby emphasizing that the dominant sensitivity of v 5,32 to β 2 emerges from v 2 .In Fig. 6 When the initial geometry is spherical (β 2 = 0), Ψ 4 and Ψ 2 could be arbitrary orientations and thus have negligible correlations.When the geometry has an anisotropy (β 2 > 0), they will have specific orientations.As a result, their correlation increases as well.Changing triaxial parameter γ, the results of χ 4,22 and ρ 4,22 remain consistent within errors, restraining themselves from an ideal probe of the triaxial structure of 129 Xe.
Normalized symmetric cumulant (NSC) quantifies the correlation between different order flow coefficients and has been systematically studied in Pb-Pb collisions [72][73][74].It was found that NSC(3, 2) is insensitive to the dynamic evolution but carries unique sensitivity to initial conditions [72,73].Thus, it is potentially a good probe of the nuclear structure.In fact, NSC(3, 2) has been previously studied in the U-U collisions, showing its sensitivity to β 2 [26].While for NSC(4, 2), it is sensitive to both initial conditions and transport properties [72,73].In Figure 7(a), a tiny (and insignificant) increase can be observed in NSC (3,2) by decreasing the β 2 value from 0.18 to 0, and no significant variation is observed when changing γ.This indicates that the anti-correlation between ⟨v 2 3 ⟩ and ⟨v 2 2 ⟩ is reduced by quadrupole deformation and is independent of γ.For Fig. 7(b), it is evident that NSC(4, 2) exhibits little to no sensitivity to either β 2 or γ, suggests that the correlation between ⟨v 4 ⟩ and ⟨v 2 ⟩ remains relatively independent of nuclear deformation.

Study on the nuclear diffuseness
The radial profile of the nucleus is determined by either the nuclear diffuseness a 0 or the radius R 0 , both of which impact the initial spatial anisotropy and consequently influence the flow observables.Recent studies have highlighted the significant impact of a 0 on various flow observables in the mid-central region (20-60% centrality) [25], underscoring the importance of verifying this within the context of Xe-Xe collisions using the AMPT model.This section delves into the potential of flow observables in investigating the a 0 introduced in Eq. ( 4).Generally, a larger nuclear diffuseness a 0 results in a significantly increased total hadronic cross-section, leading to a more diffused QGP.
Fig. 8 shows the centrality dependence of v 2 {2} and v 3 {2} in Xe-Xe collisions at 5.44 TeV using two distinct values of a 0 .Calculations with larger a 0 values (presented by black markers) lead to smaller v 2 results in the semi-central collisions, which is agreed with the results in the isobar runs of Ru-Ru/Zr-Zr [25].Meanwhile, no obvious change in v 3 is seen after changing the a 0 values across the 0-60% centrality range.
As discussed in section 4.1, the increase of σ v2 is mainly due to the random orientations of deformed nuclei and is expected to be less influenced by the diffuseness of nuclei.To verify this, the results of v 2 {4}, v2 , and σ v2 are presented in Fig. 9.It is observed that v 2 {4} and v2 are reduced by larger a 0 in the semi-central region, while σ v2 remains the same.This observation aligns with our earlier discussion, affirming that the diffuseness primarily impacts the eccentricity rather than its fluctuations.
In exploring a 0 , we can split v 4 into its linear and nonlinear components to see how they respond differently to changes in a 0 .Figure 10 shows the centrality dependence of v 4 {2} and its linear and nonlinear components.A larger a 0 value reduces both v 4 {2} and v 4,22 in 20-40% centrality.Figure 10(c   smaller than v L 4 in the most central region, so we do not see any effect of β 2 on v 4 {2} even though β 2 increases v 4,22 significantly in that region.For a 0 , more sensitivity is found in the mid-central collisions, where the nonlinear component v 4,22 is about half of the total v 4 {2}, and how it changes with a 0 has a bigger effect on v 4 {2} in this region.As also introduced in section 4.1, χ 4,22 is not only affected by dynamic evolution but also by initial conditions, while ρ 4,22 only depends on the initial state.It is crucial to check if the nonlinear modes can probe the initial conditions, particularly the diffuseness of nuclei.Here χ 4,22 and ρ 4,22 with different a 0 are presented in Fig. 11.Neither of them is affected by a 0 in 0-60% centrality region.As ρ 4,22 = v 4,22 /v 4 {2}, the sensitivity to a 0 is, to a large extent, canceled by taking this ratio.The insensitivity of χ 4,22 to a 0 can lead to similar conclusion as in the discussion of β 2 that the sensitivity of v 4,22 to a 0 is mainly due to the V 2 2 term in V NL 4 ≈ χ 4,22 (V 2 ) 2 .Because of the limited statistics, there are large uncertainties in the results of NSC(3,2) and NSC(4,2), and no firm conclusion about how sensitive they are to the a 0 parameter can be drawn.Thus, the results are not presented and discussed here but added in the appendix.

Influence from the transport properties
The results shown in this study have indicated that v 2 {2}, v2 , σ v2 , v 4,22 and ρ 4,22 can potentially probe the nuclear deformation, meanwhile v 2 {2}, v 2 {4}, v2 , v 4 {2}, v 4,22 have the ability to study nuclear diffuseness a 0 .In order to achieve unbiased constraints on the nuclear structure parameters in experiments, one must ensure that the chosen observables are only sensitive to the initial conditions and not influenced by the dynamic evolution of the created system.In particular, it has been found in the study at RHIC isobar runs that the ratio observables in Zr-Zr and Ru-Ru can largely cancel out the effects of final state interactions and thus reflect mainly the impact from the initial stages.Here we present the mentioned observables (v 2 {2}, v 2 {4}, v2 , σ v2 , v 4 {2}, v 4,22 , ρ 4,22 ) in the ratio of Xe-Xe/Pb-Pb, checking if the ratio of two collision system minimizes the influence of dynamic evolution.
In Fig. 12, the centrality dependence of v 2 {2}, v 2 {4}, v2 and σ v2 in the ratio of Xe-Xe/Pb-Pb collisions are presented, using two different partonic cross sections.The bottom panels of each figure show the ratio of the results using 3 mb and 6 mb ("double ratio").The ratios in Fig. 12 are all consistent with unity, indicating that these observables, including v 2 {2}, v 2 {4}, v2 , σ v2 , in the ratio of Xe-Xe/Pb-Pb are less dependent on the dynamic evolution and potentially are good probes to initial conditions.
In addition, v 4 {2} and its nonlinear modes v 4,22 and ρ 4,22 in the ratio of Xe-Xe/Pb-Pb are also investigated, the centrality dependence is presented in Fig. 13.Despite large uncertainties, the results with different partonic cross sections are consistent with each other, and the "double ratio" is compatible with unity.Such results show that v 4 {2}, v 4,22 , ρ 4,22 in the ratio of Xe-Xe/Pb-Pb might be not affected by dynamic evolution.Considering v 4 {2} and v 4,22 are sensitive to a 0 in the mid-central region, while v 4,22 and ρ 4,22 are enhanced by non-zero β 2 in the most central region, these observables are valuable in the exploration of initial conditions and nuclear structure parameters.

Summary
This paper presents a comprehensive exploration of the influence of nuclear deformation and nuclear diffuseness on various flow observables in Xe-Xe collisions at √ s NN = 5.44 TeV, using AMPT event generator.We observe that the elliptic flow coefficients v 2 {2}, v2 , elliptic flow fluctuation σ v2 , as well as the nonlinear flow mode observables v 4,22 and ρ 4,22 , are enhanced in central collisions in the presence of nuclear quadrupole deformation β 2 .These enhancements come from the increased eccentricities and their fluctuations in the initial geometry in the presence of a deformed shape of 129 Xe.Thus, the aforementioned flow observables can potentially constrain nuclear deformation through data-model comparisons in high-energy heavy-ion collisions in the near future.On the other hand, most flow observables do not exhibit much sensitivity to the triaxiality parameter, underscoring the importance of utilizing combined flow observables like v n − [p T ] correlations in future studies.In addition, larger nuclear diffuseness a 0 for 129 Xe leads to smaller values in v 2 {2}, v2 , v 4 {2}, v 4,22 , and v 2 with multi-particle correlations (v 2 {4}, v 2 {6}, v 2 {8}) in midcentral collisions.Among them, v 4 {2}, v 2 {4}, v 2 {6}, and v 2 {8} are only sensitive to a 0 , while σ v2 and ρ 4,22 are only sensitive to β 2 .The differing sensitivities of these observables allow for separating constraints on the values of a 0 and β 2 , which will help fine-tune other model parameters and enable more precise predictions.
We also investigate the influence of dynamic evolution on flow observables by varying the parton cross-section in the AMPT model.The parton cross-section does not significantly affect most of the observables in the ratio of Xe-Xe/Pb-Pb, indicating that the final state effect is canceled out when comparing the two collision systems.Future collisions of different nuclear species at varying energies will provide more insights into heavy ion collisions and improve our nuclear structure knowledge.

Acknowledgements
M. Zhao and Y.Zhou are funded by the European Union (ERC, InitialConditions), VILLUM FONDEN (grant number 00025462), and Danmarks Frie Forskningsfond (Independent Research Fund Denmark).J. Jia's work is supported by the US Department of Energy (grant number DE-FG02-87ER40331).
) does not show the sensitivity of v L 4 to a 0 due to the large uncertainties.The reduction of the nonlinear part v 4,22 affects v 4 {2}.This is different from what we saw for β 2 .Also, it is seen that v 4,22 is much

Table 1 .
parameter sets used in this study