Predictions on global properties in O+O collisions at the Large Hadron Collider using a multi-phase transport model

Oxygen ($^{16}$O) ions are planned to be injected at the Large Hadron Collider (LHC) in its next runs, and a day of physics run is anticipated for O+O collisions at $\sqrt{s_{\rm{NN}}}$ = 7 TeV. As the system size of O+O collisions has the final state multiplicity overlap with those produced in pp, p+Pb and Pb+Pb collisions, the study of global properties in O+O collisions may provide a deeper insight into the heavy-ion-like behavior observed in small collision systems and its similarities/differences with a larger system like Pb+Pb collisions. In the present work, we report the predictions for global properties in O+O collisions at $\sqrt{s_{\rm{NN}}}$ = 7 TeV using a multi-phase transport model (AMPT). We report the mid-rapidity charged-particle multiplicity, transverse mass, Bjorken energy density, pseudo-rapidity distributions, squared speed of sound, transverse momentum ($p_{\rm T}$) spectra, the kinetic freeze-out parameters, and $p_{\rm T}$-differential particle ratio as a function of collision centrality. Further, we have studied the transverse momentum-dependent elliptic flow of charged particles. The results are shown for Woods-Saxon and harmonic oscillator nuclear density profiles. In addition, we have compared the results with an $\alpha$-clustered structure incorporated inside the oxygen nucleus. Average charged-particle multiplicity and the Bjorken energy density show a significant increase in most central collisions for the harmonic oscillator density profile, while other global properties show less dependence on the density profiles considered in this work. The results from the $\alpha$-clustered structure incorporated inside the oxygen nucleus show similar initial energy density and final charged-particle multiplicity as observed for the harmonic oscillator density profile.


I. INTRODUCTION
In order to understand the properties of the hot and dense medium, often referred to as Quark-Gluon Plasma (QGP), formed in ultra-relativistic heavy-ion collisions at the Large Hadron Collider (LHC) and Relativistic Heavy Ion collider (RHIC), several measurements in different collision systems at different center-of-mass energies are performed. Historically in these collider experiments, heavy-ion collisions such as Pb+Pb and Au+Au collisions are the primary focus in the study of QGP while the small collision systems such as proton-proton (pp) collisions act as a baseline. However, recent results from the LHC experiments [1,2] show QGP-like properties in high-multiplicity pp collisions, which raises concerns in the heavy-ion physics community about whether pp collisions can act as a baseline and in addition, if QGPdroplets are produced in the collisions of small systems at the LHC energies. This ambiguity has a serious consequence on the results reported for heavy-ion collisions. Thus, a closer look at the small collision systems is a call of time.
In the upcoming run at the LHC, brief oxygen-oxygen (O+O) collisions are anticipated [3], which has a final individually, the oxygen nucleus is doubly magic and it is expected to be highly compact. Oxygen has also got an isotope ( 24 O), which is doubly magic with eight protons and sixteen neutrons but unstable in nature. Readers can refer to Ref. [16] for the details of the nuclear structure of 24 O, as an outlook.
To get a deeper insight into the impact of nuclear structure inside the oxygen nucleus, we have incorporated simple harmonic oscillator potential and a more realistic Woods-Saxon potential in the oxygen nucleus using a multi-phase transport model (AMPT) [17]. Also, we have compared the results with an α-clustered structure incorporated inside the oxygen nucleus. Global observables like charged-particle multiplicities, transverse energy, particle spectra, and pseudorapidity distributions provide insight into the possible formation of QGP in a system. It is proposed that the equation of state of hot hadronic matter can be probed via the correlation of mean transverse momentum and particle multiplicity [18]. Charged-particle multiplicity provides information about the soft processes in the collision, while the mean transverse mass and momentum give insight into the hard processes. In recent results from the LHC [1], it has been shown that the QGP-like properties seen in high-multiplicity pp collisions are driven by the final state multiplicity in an event. Thus, it would be interesting to confront initial and final state effects in O+O collisions as it has multiplicity overlap with high-multiplicity pp collisions. The initial energy density can be estimated by using Bjorken hydrodynamic model [19], where one uses the transverse energy or charged-particle multiplicity density in rapidity and mean transverse mass for each collision centrality. To explore the final state effects, a study on particle spectra, kinetic freeze-out parameters, and particle ratios can be studied. In this work, the global properties [20] such as Bjorken energy density, squared speed of sound, particle ratios, and kinetic freeze-out parameters are studied for O+O collisions at √ s NN = 7 TeV using AMPT for both harmonic oscillator and Woods-Saxon nuclear density profile as well as with an α-clustered structure incorporated in oxygen nucleus, in order to explore the effect of the nuclear density profile of the final state observables.
The paper is organized as follows. We begin with a brief introduction to O+O collisions in Section I. In Section II, the detailed event generation methodology with AMPT along with different types of nuclear density profiles are discussed. In section III we give a detailed discussion of the obtained results. Finally, the results are summarized in Section VI with possible outlook.

II. EVENT GENERATION AND ANALYSIS METHODOLOGY
In this section, we begin with a brief introduction to the AMPT model and then discuss the charge density profile of oxygen.

A. A Multi-Phase Transport (AMPT) Model
AMPT model primarily consists of four components [17]: 1. initialization of the collisions using HI-JING [21], where the differential cross-section of the produced mini-jets in pp collisions is calculated and converted to heavy-ion collisions with the inbuilt Glauber model, 2. the produced partons are propagated into parton transport part via Zhang's Parton Cascade (ZPC) model [22], 3. hadronization mechanism: in AMPT string melting version, the transported partons are hadronized using spatial coalescence mechanism [23,24], 4. hadron transport: the produced hadrons undergo a final evolution in relativistic transport mechanism [25,26] via baryon-baryon, meson-baryon, and meson-meson interactions. Although there is also a default version available in AMPT, we have used the string melting mode of AMPT (AMPT version 2.26t9b) in the current work due to the fact that the particle flow and spectra at the intermediate-p T regions are well explained by a quark coalescence mechanism available in the string melting version for hadronization [27][28][29]. We have used similar AMPT settings in the current work as reported in Ref. [30] unless specified explicitly. As there is no experimental data available for O+O collisions, we have compared the charged-particle p T -spectra from minimum bias p+Pb collisions at √ s NN = 5.02 TeV in midrapidity with the predictions from AMPT by fixing the partonic scattering cross section, σ gg = 3 mb [4], which can be seen in the Appendix section A of the paper. The collision centrality, the number of participants (N part ), and the number of binary collisions (N coll ) in heavy-ion collisions cannot be determined directly from experiments, instead it is obtained via the impact parameter determination in the Glauber model [31][32][33]. The Glauber model considers a nucleus-nucleus collision as a superposition of several independent nucleon-nucleon collisions and the estimations are dependent on the nuclear density profile. In this work, we have modified the   inbuilt Glauber model in AMPT to incorporate a harmonic oscillator density profile and α-clustered structure in the oxygen nucleus. They are compared with the results with that of the Woods-Saxon nuclear density profile. The values of impact parameters for different centrality classes in O+O collisions are obtained by using the publicly available MC Glauber code (TGlauberMC-3.2) [34][35][36][37]. The values of the impact parameter, the average number of nucleon participants, and average number of binary nucleon-nucleon collisions for different density profiles and in different centrality classes are shown in Table I.

B. Woods-Saxon density profile
In high-energy heavy-ion collisions, the standard method employed for nuclear density profile is in terms of three-parameter Fermi (3pF) distribution, which is often referred as Woods-Saxon distribution. The Woods-Saxon distribution is given by, Here, r is the radial distance from the center of the nucleus, r 0 is the mean radius of the nucleus, a is the skin depth/diffusivity of the nucleus and w is the deformation parameter. For the oxygen nucleus, r 0 is 2.608 fm, a is 0.513 fm, w is -0.051 [37]. ρ 0 is the nuclear density constant at r = 0, which is obtained by the overall normalization condition, Here, Z is the atomic number of the nucleus i.e., 8 for oxygen nucleus. For a hard sphere configuration in r < r 0 , ρ(r) = ρ 0 , and ρ 0 = 3Ze/(4πR 3 ).

C. Harmonic oscillator density profile
The harmonic oscillator charge density distribution is given by, Here, a and α are parameters which are taken as 1.544 fm and 1.833 for oxygen [37], respectively. Similar to the Woods-Saxon case, ρ 0 is the nuclear density constant at r = 0, which should satisfy the normalization condition mentioned in Eq. 2. Figure 1 shows the comparison of different normalized nuclear density profiles of oxygen, namely the hardsphere, Woods-Saxon, and the harmonic oscillator. A significant rise of the normalized density profile at small r is due to the harmonic oscillator potential seen with respect to the Woods-Saxon density profile. As the tetrahedral α-cluster structure of 16 O nucleus is implemented numerically in AMPT, a nuclear density profile is difficult to obtain analytically. However, the probability of the radial position of the nucleons distributed inside the nucleus has been shown in Appendix section B and compared the same for harmonic oscillator potential and the Woods-Saxon density profile. A brief discussion on the details of tetrahedral α-cluster structure is mentioned be-low. Clustering plays a crucial role in studying nuclear structure. Protons and neutrons in many-body nuclear systems tend to form clusters in order to reduce overall energy or boost system stability. Two neutrons and two protons in a nucleus can cluster together to form α particles. The α-clustering is observed in several light nuclei such as 8 Be and 12 C [38]. In these types of nuclei, the mean-field is not strong enough to break the cluster effect. The intrinsic state of the 12 C nucleus is a triangular structure made up of 3 α-particles [39]. Several studies [38,40] have indicated that the structure in 16 O corresponds to a state analogous to the 12 C state with α-clustering. The experimental observation assumes that the 16 O has an α-like cluster at the corners of a tetrahedron (Fig. 2). A signature of α−clustering is seen in recent simulation studies for collisions involving oxygen nuclei [9]. These studies have inspired us to implement the α−clustered structure inside the oxygen nucleus using the AMPT model. Nucleons inside an α−cluster are distributed following the Woods-Saxon distribution for 4 He nucleus with a rms radius 1.676 fm. Such randomized α−clusters are placed on the vertices of a regular tetrahedron with a side length 3.42 fm. The rms radius of such an arrangement gives the rms radius for 16 O to be 2.699 fm [9]. The arrangement of nucleons is randomized event by event by rotating the system along x-y-z directions for both projectile and target nuclei following the tetrahedral structure.

III. RESULTS AND DISCUSSIONS
Here, we discuss the predictions for different global properties such as Bjorken energy density, speed of sound, freeze-out parameters as a function of centrality class in different subsections. We have also studied the transverse momentum-dependent elliptic flow of charged particles. Here onwards for simplicity, we refer π + + π − , K + +K − and p+p as pions, kaons, and protons, respectively. In the study of QGP properties in heavy-ion collisions, one of the key variables is the initial energy density produced in such collisions. The initial density can be estimated via the Bjorken boost-invariant hydrodynamics model [19], where transverse energy density (E T ) in midrapidity gives the quantitative estimation of the initial energy density produced in the interaction. The Bjorken energy density ( Bj ) with the assumption of boost invariance is given as,

A. Transverse energy and Bjorken Energy Density
where, S T is the transverse overlap area of the two colliding nuclei and dE T /dy is the transverse energy density at midrapidity at a formation time τ . As Eq. 4 diverges at τ → 0, a finite formation time (τ ) = 1 fm/c is assumed for the calculation of Bjorken energy density in this work. E T is the total transverse energy produced in an event and S T = πR 2 is the total transverse overlap area of the colliding nuclei. Replacing, R = R 0 A 1/3 and A = N part /2, As most of the transverse energy is carried by the pions, kaons and protons due to their abundance, the total transverse energy (E T ) can be approximated as [41][42][43], The multiplicative factor in each term accounts for corresponding neutral particles. m T = p 2 T + m 2 , is the transverse mass and dN/dy is the integrated yield for π ± , K ± and p +p in mid-rapidity region i.e. |y| < 0.5. Now, Eq. 4 can be written as, Figures 3 and 4 show the integrated yields and mean transverse momenta for pions, kaons and protons as a function of centrality classes for O+O collisions at √ s NN = 7 TeV, respectively. The solid markers represent the Woods-Saxon density profile, the open markers represent the results for the harmonic oscillator density profile, and the markers with a cross represent the results from α-clustered structure. As expected, pions' integrated yield is higher than kaons and protons as they are the most abundant among the identified particles and are understood from a thermalized Boltzmannian production of secondaries in nuclear collisions. The comparison of Woods-Saxon and harmonic oscillator density profiles show that the soft particle production is higher for the harmonic oscillator density profile in central collisions. The results from the α-clustered structure show similar behavior as seen for the harmonic oscillator density profile. However, the mean transverse mass remains nearly the same, indicating a similar spectral shape for all cases. For peripheral collisions, the differences among the density profiles and α-clustered structure diminish.
With the input of integrated yields and mean transverse mass for pions, kaons, and protons, the Bjorken energy density is obtained as a function centrality classes using Eq. 7, which is shown in Fig. 5. The solid markers represent the Woods-Saxon density profile, the open markers represent the harmonic oscillator density profile and the markers with a cross represent the results for the α-clustered structure. The Bjorken energy density is found to be higher for most central collisions and linearly decreases while going from central to peripheral collisions. As evident in Eq. 7, Bjorken energy density strongly depends on the integrated yield, and the difference in integrated yield for different nuclear density profiles is reflected in the Bjorken energy density. The oxygen nucleus with harmonic oscillator density profile show about 15% higher energy density compared to the Woods-Saxon density profile. Thus, for central O+O collisions, the density profile plays a crucial role in studying the initial-state (Bjorken energy density) and final-state (integrated yields) effects. However, going towards the higher impact parameter, i.e., for peripheral collisions the difference is negligible. As observed for integrated yield, α-clustered structure shows similar behavior as seen for harmonic oscillator density profile. The nuclear density profile has a clear effect on the initial energy density, which in fact controls the particle production and the subsequent space-time evolution of the fireball, and hence the equation of state (EoS). The values of the initial energy densities for all collision centralities are observed to be higher than the lattice QCD estimated requirement of 1 GeV/fm 3 energy density for a deconfinement transition [44]. This hints at observing the signals of QGP in oxygen-oxygen collisions at the LHC energies. We shall further strengthen these arguments in the following sections while discussing other global observables in heavyion collisions.  and (70-100)% centrality classes for Woods-Saxon density profiles, harmonic oscillator density profiles and α-clustered structure. Here, the differences in the chargedparticle multiplicity in the mid-pseudorapidity due to modification of nuclear density profiles are clearly visible for central collisions while the difference is relatively smaller for peripheral collisions. However, at forwardpseudorapidity, almost no dependence of charged-particle density on the density profiles and α-clustered structure is seen. In Landau hydrodynamical model [45], the rapidity distributions are expected to follow a Gaussian distribution and in this framework, the speed of sound (c s ) is related to the width of the rapidity distribution via the following expression.
Here, σ y is the width of the rapidity distribution and m p is the mass of a proton. The massless ideal gas limit for the squared speed of sound (c 2 s ) is 0.33. However, the presence of the dip structure in pseudorapidity distributions at around |η| = 0 makes it difficult to fit via a single Gaussian distribution. Usually, in experiments [46,47], the following double Gaussian distribution is used to fit the pseudorapidity distributions to extrapolate the distributions to the unmeasured regions: Here, A 1 and A 2 are the normalization parameters while σ 1 and σ 2 are the widths of each Gaussian distribution. Thus, using Eq. 9 for the fitting of pseudorapidity distributions, we have obtained σ 1 and σ 2 . The values of σ 1 and σ 2 are found to be similar within uncertainties, which can be seen in Table III. Figure 7 shows the squared speed of sound as a function centrality classes for pions, kaons and protons in O+O collisions at √ s NN = 7 TeV for Woods-Saxon, harmonic oscillator density profiles and α-clustered structure obtained using Eq. 8. Here, σ 1 has been used to obtain the absolute value and the maximum deviation of σ 1 and σ 2 is used as uncertainties for c 2 s . Considering the uncertainties, c 2 s is found to be similar as a function of centrality. The negligible dependence of c 2 s on centrality classes could give an indication that the system produced in O+O collisions is significantly less dense compared Pb+Pb collisions and they are more similar to pp collisions. The ideal gas limit is also shown as a red dotted line in Fig. 7 and the observed values of c 2 s is found to be around 27% lower than the massless ideal gas limit.   Woods-Saxon, harmonic oscillator density profiles and α-clustered structure in oxygen nucleus. The fitting ranges for each particle are similar to the range reported by ALICE [48]. The fitting is performed via the χ 2minimisation method and the values of χ 2 per degree of freedom are shown in each of the cases in Fig 8. The expression for invariant yield in the Boltzmann-Gibbs blastwave framework [49] is given as the following: Here C is the normalisation constant. The particle four-momentum is given by, p µ = (m T coshy, p T cos φ, p T sin φ, m T sinhy), (11) and the four-velocity is given by, u µ = cosh ρ (cosh η, tanh ρ cos φ r , tanh ρ sin φ r , sinh η). (12) Finally, the freeze-out surface is parametrised as, here, η is the space-time rapidity. Now, Eq. 10 is expressed as, K 1 m T coshρ T kin and I 0 p T sinhρ T kin are modified Bessel's functions, which are given by, where, ρ is given by ρ = tanh −1 β T and β T (= β s ξ n ) [49][50][51][52] is the radial flow velocity. Here, ξ is given as (r/R 0 ), β s is the maximum surface velocity and r is the radial distance. R 0 is the maximum radius of the fireball at freeze-out. In this model, the particles closer to the center of the fireball are assumed to move slower than the ones at the edges. The average of the transverse velocity is evaluated as [53], For our calculation, we use n as a free parameter. Figure 9 shows the kinetic freeze-out temperature versus transverse radial flow velocity from the simultaneous fit of identified particles' p T -spectra with Boltzmann-Gibbs blastwave distribution (Eq. 14) in O+O collisions at √ s NN = 7 TeV. The solid (open) markers represent the Woods-Saxon (harmonic oscillator) density profile and the markers with a cross represent α-clustered structure nucleus. Within uncertainties, the correlation between kinetic freeze-out temperature (T kin ) and average transverse flow ( β T ) is similar for all cases. For the most central collisions (0-5% class), T kin is the lowest and the transverse flow is the highest. This behavior is expected as the most central collisions have the largest system size due to which the hadronic phase lasts longer, which makes the T kin lowest. Also, due to the largest system size, the radial flow is expected to be the highest. A similar behavior is seen for Pb+Pb collisions at LHC energies [48].   Figure 10 shows p T -differential kaon-to-pion and proton-to-pion ratios in O+O collisions at √ s NN = 7 TeV for (0-5)% and (70-100)% centrality classes. Both the ratios to pions increase as a function of p T . As the kaonto-pion ratio is the measure of strangeness, we see enhancement of strangeness production as a function of p T . This enhancement seems to be the similar for both (0-5)% and (70-100)% centrality classes at low-p T . However, the enhancement is higher for (0-5)% centrality class at intermediate-p T . The proton-to-pion ratio, a ratio between the lightest baryon to lightest meson, acts as a proxy for the baryon to meson ratio. In general, the particle ratios do not show any dependence on the density profiles for (70-100)%. However, for (0-5)% centrality class, a rise in both kaon-to-pion and proton-to-pion ratios is seen at intermediate-p T for harmonic oscillator density profile with respect to Woods-Saxon density profile, which becomes prominent at higher p T . The results for the α-clustred structure are similar to that observed for the harmonic oscillator density profile.

V. ELLIPTIC FLOW
The initial spatial anisotropies of the overlap region in non-central heavy-ion collisions get converted to the momentum space azimuthal anisotropy of the final state particles due to the differential pressure gradients of the produced medium. The azimuthal anisotropy can be expressed as a Fourier series in the azimuthal angle, φ: where, v n is the anisotropic flow of different order n, with n = 2 is elliptic flow. ψ n is the n th harmonic event plane angle [54]. To reduce the non-flow effects, a two-particle correlation method [55,56] can be adopted to estimate the elliptic flow. The correlation function between two particles is obtained in relative pseudorapidity (∆η = η a − η b ) and relative azimuthal angle (∆φ = φ a − φ b ).
Here, a and b denote two separate particles in a pair.
In this study, we have taken charged particles in |η| < 2.5 and p T > 0.5 GeV/c to be consistent with previous studies [55,56]. The 1D correlation function is given as, where, the normalization constant (A) ensures that the number of pairs are the same between signal, S(∆η, ∆φ) and background, B(∆η, ∆φ). The ∆η interval is chosen carefully, as done in previous studies [55,56], by excluding the jet peak region observed in the C(∆η, ∆φ) distribution. The interval is taken to be 2.0 < |∆η| < 4.8. This pseudorapidity cut removes the residual non-flow effects significantly in the estimation of elliptic flow. These non-flow correlations usually arise from jets and shortrange resonance decays, and they are not associated with the anisotropy in the early stage of the collisions. In the current study, the non-flow effects are reduced significantly but they might be still non-zero. Thus, the quantitative interpretation of the results may be taken with caution.
The pair distribution can be expanded into a Fourier series: where, v n,n is the two-particle flow coefficient. Now, Eq. 17 is given as: The definition of harmonics defined in Eq. 16 now enters in Eq. 18, Through this definition, the event plane angle drops out in convolution. Thus, if the elliptic flow is driven by purely collective expansion then v 2,2 should be factorized into the product of two single-particle elliptic flow coefficients.
Using Eq. 21, we calculate the single particle elliptic flow coefficient as, We now proceed for the estimation of the elliptic flow of charged particles as a function of transverse momentum for (0-5)% and (20-30)% centrality classes in O+O collisions at √ s NN = 7 TeV using AMPT for Woods-Saxon and harmonic oscillator density profiles along with α-cluster structure considered in this study. These are shown in Figure 11. v 2 is found to increase at low-p T and saturates at intermediate-p T . The qualitative trend of elliptic flow as a function of p T in O+O collisions is found to be similar to that is observed by ALICE experiment in p+Pb collisions at √ s NN = 5.02 TeV [57].
Within uncertainties, v 2 is found to be similar for both Woods-Saxon and harmonic oscillator type density profiles for (0-5)% and (20-30)% centrality classes. This indicates that, although there is a significant dependence of the initial energy density on the density profiles, the collectivity (both radial and elliptic flow) is less affected by the modification in these two density profiles. However, when compared with results for α-cluster structure in oxygen nucleus, elliptic flow at intermediate-p T is found to have a centrality dependent trend as compared to the other two density profiles considered in this study. It is found that, for the most central case, elliptic flow for the α-cluster nucleus is slightly less compared to the Woods-Saxon type nucleus. However, this trend seems to reverse as one moves to mid-central case, where v 2 is quantitatively more for the α-cluster nucleus. To get a further insight into α-cluster structure in oxygen nucleus and the similarities/differences with different density profiles, a detailed study on the centrality and transverse spherocity dependence of elliptic flow [55] needs to be performed. It would be also interesting to compare the predictions of elliptic flow from AMPT with upcoming experimental results in O+O collisions.

VI. SUMMARY
In summary, we report the predictions for global properties in O + O collisions at √ s NN = 7 TeV using a multi-phase transport model (AMPT). We report the mid-rapidity charged-particle multiplicity, mean transverse mass, Bjorken energy density, pseudorapidity distributions, squared speed of sound, p T spectra, and the kinetic freeze-out parameters as a function of collision centrality. The results are shown for both harmonic oscillator and Woods-Saxon nuclear density profiles along with α-clustered structure incorporated for oxygen nucleus. With the change of the density profile from Woods-Saxon to the harmonic oscillator and with the implementation of α-clustered in the nucleus, a modification of average charged-particle multiplicity is seen, which is also reflected in the initial energy density as one naively expects. However, other global properties show less dependence on the density profile considered in this work. In general, the initial energy density produced in all collision centralities in O+O collisions at √ s NN = 7 TeV stays higher than the lattice QCD predicted value for a deconfinement transition, making oxygen nuclei collisions a potential case of light-nuclei collisions at the LHC energies to create a state of QGP. In addition, a substantial radial flow and a comparable freeze-out temperature with that of Pb-Pb collisions are also observed in the present analysis. Although a significant dependence of the initial energy density on the density profiles is seen, the collectivity (both radial and elliptic flow) is less affected by the modification in the density profiles. When compared to the results from the α-clustered structure incorporated in the oxygen nucleus with different density profiles, the initial energy density, and charged-particle multiplicity are found to be similar for the α-clustered structure and harmonic oscillator density profile. However, the magnitude of elliptic flow at intermediate-p T for α-cluster and Woods-Saxon type nucleus is found to reverse the trends when studied as a function of different centrality classes. It would be interesting to confront these results with the experimental observations, when available, to reveal the density profile of the oxygen nucleus best suitable to describe ultra-relativistic nuclear collisions. As the experimental data will be available for O+O collisions in the future (LHC RUN3), we have compared the charged particle p T -spectra from a collision system having a closer system size to O+O collisions i.e., minimum bias p+Pb collisions at √ s NN = 5.02 TeV with the predictions from AMPT. We have compared the AMPT predictions modifying the partonic scattering cross-section, σ gg to different values such as 3, 5, and 10 mb. From this exercise, we found that fixing σ gg = 3 mb, the AMPT predictions are closer to experimental data in p+Pb collisions at √ s NN = 5.02 TeV. Figure 12 shows the comparison of ALICE data with the predictions from AMPT with σ gg = 3 mb for p-Pb collisions at √ s NN = 5.02 TeV. It is found that the spectral shape from AMPT matches with the experimental data at high-p T while at intermediatep T , 10-30% difference is observed. One should also note here that, to exactly match the experimental data, one can vary the tunes of the AMPT model, which is currently out of the scope of this manuscript.  Figure 13 represents the probability of the radial positions of the nucleons distributed inside the oxygen nucleus for Woods-Saxon, harmonic oscillator density profiles, and α-cluster structure. The α-cluster structure seems to have a compact radial distribution of nucleons as compared to the other two density profiles.