#### Abstract

We have estimated centrality variation of chemical freeze-out parameters from yield data at midrapidity of , and , for collision energies of RHIC (Relativistic Heavy Ion Collider), Beam Energy Scan (RHIC-BES) program, and LHC (Large Hadron Collider). We have considered a simple hadron resonance gas model and employed a formalism involving conserved charges () of QCD for parameterization. Along with temperature and three chemical potentials (), a strangeness undersaturation factor () has been used to incorporate the partial equilibration in the strange sector. Our obtained freeze-out temperature does not vary much with centrality, whereas chemical potentials and seem to have a significant dependence. The strange hadrons are found to deviate from a complete chemical equilibrium at freeze-out at the peripheral collisions. This deviation appears to be more prominent as the collision energy decreases at lower RHIC-BES energies. We have also shown that this departure from equilibrium reduces towards central collisions, and strange particle equilibration may happen after a threshold number of participants in - collision.

#### 1. Introduction

In the last few decades, several ion collider experiment collaborations have been developed to explore the phase diagram of quantum chromodynamics (QCD). Relativistic Heavy Ion Collider (RHIC) started to investigate the signatures of deconfined quark-gluon plasma, whereas the RHIC Beam Energy Scan (BES) program became motivated in searching the QCD critical point [1]. The Large Hadron Collider (LHC) is trying to investigate the medium created in zero baryon density, where a crossover transition from a hadronic state to a deconfined state of quarks and gluons may happen [2].

The hot and dense fireball created in these collisions experiences a fast expansion due to the initial pressure gradient. Assuming that the system starts from a strongly interacting quark-gluon state, a fast thermalization can occur. The thermodynamic parameters like temperature and chemical potentials can describe this thermalized medium. The matter and energy density dilute with the expansion and the temperature declines. As the energy density(temperature) drops below the hadronization threshold, the matter evolves as a state of hadrons and their resonances. The mean free path increases with further expansion and various collisions among particles abate. In this context, one can define freeze-out as the boundary, onwards which no interaction is supposed to happen among hadrons. In the standard description, two freeze-out surfaces are described, depending on the interaction type. The chemical freeze-out (CFO) happens when inelastic scattering stops and the particle abundances become fixed. The kinetic freeze-out (KFO) is the point where elastic collisions cease. In this free noninteracting limit, the ideal hadron resonance gas model may give a reasonable description of the hadrons at freeze-out.

Yields of strange hadrons help to understand the extent of chemical equilibrium achieved in these collisions. The strange quark equilibrates later than the , quarks due to its larger mass [3]. Equilibrated strangeness spectra may be a suitable signature to understand the existence of a deconfined partonic phase [4]. Recently, a charged particle multiplicity () dependent production of strangeness has been observed in LHC [5]. This strangeness production is related to the collision centrality and number of participants () of the collision system [6, 7]. Here, denotes the average number of participating nucleons of a particular collision system.

Following the success of the Statistical Hadronization Model (SHM) [8], studies have determined the chemical freeze-out parameters, considering the Hadron Resonance Gas (HRG) model. In this context, a fitting with the available yield data is well practiced [9–20]. Generally, one extracts the chemical freeze-out temperature and the baryon chemical potential by the minimization procedure, whereas the charge chemical potential and the strange chemical potential get fixed from the constraints of the colliding nuclei. To scale the possible nonequilibration of strange hadrons, a strangeness under saturation factor can be introduced [13, 21–27]. This parameter scales the deviation of strange hadrons from a complete equilibrium in the Grand Canonical Ensemble (GCE). In a recent work [28], we have shown that in analysis, a larger systematic variation can arise depending on the chosen set of ratios. A conserved-charge-dependent extraction of thermal parameters has been proposed [29], which seems to suitably estimate thermal parameters and predict equilibration in the most central collision. It will be interesting to check the centrality variation of thermal parameters and the equilibration of strange particles in this framework.

In this manuscript, we have tried to study the centrality variation of freeze-out parameters, with an emphasis on the saturation of strangeness equilibration [5] in heavy nuclei collision. We have observed a similar saturation with centrality () for collision energies ranging from 7.7 GeV of RHIC-BES to LHC (2.76 TeV). We have employed a recently developed parameter extraction process [29, 30] with a strangeness suppression factor to measure the possible deviation of strange hadrons from respective equilibrium yield. We have found the kaons to deviate from equilibrium at chemical freeze-out for the peripheral and semiperipheral collision, though the temperature does not change much with centrality. We have further studied the scaling behavior of all freeze-out parameters. The parameters attain a saturation onwards . This flattening indicates that the system created in the heavy-ion collision reaches a grand canonical limit corresponding to the most central value, in which thermodynamic description becomes independent of the system size. Finally, we have verified the efficacy of our parametrization by comparing our estimated hadron yield ratios with experimental data.

We have organized the manuscript as follows. A short description of the Hadron resonance gas model (HRG) is given in Section 2. In Section 3, we have briefly discussed the parameter extraction method in our approach. 4 describes our results followed by discussion. We summarise our results in Section 5.

#### 2. Hadron Resonance Gas Model

The hadron resonance gas (HRG) model describes the system as a mixture of hadrons and their resonances. It is a standard exercise to incorporate all available hadron yields for obtaining a good description of the medium. In recent years various studies have been performed using HRG model [11, 12, 14, 31–57]. This model has successfully described hadron yields from AGS to LHC energies [11, 12, 34, 35, 37–39, 43]. Bulk properties of hadronic matter have also been studied in this model [14, 41, 42].

In the present work, we have considered the ideal HRG model in which hadrons are treated as point-like particles. A grand canonical ensemble can describe the partition function of hadron resonance gas as [14]

The sum runs over all hadrons and resonances. In the idealistic scenario of a chemical freeze-out, we can neglect all dissipative interactions and finite volume corrections. The thermodynamic potential for th species is given as where the upper sign is for baryons and lower for mesons. Here, is the volume and is the temperature of the system. For the species of hadron, , , and are, respectively, the degeneracy factor, energy, and mass, while is the chemical potential, with , , and denoting the baryon number, strangeness, and the electric charge, respectively. For a thermalized system, the number density can be calculated from partition function as

#### 3. Application to Freeze-Out

We first outline the usual application of the HRG model for characterizing the freeze-out temperature and chemical potentials in the context of heavy-ion collision experiments. The rapidity density for th hadron may be related to the corresponding number density as [13] where the subscript denotes the detected hadrons. Here, the total number density of any hadron is

The summation is over the heavier resonances () that decay to the hadron. This number density is calculated using Eq. (3).

In this context, it is also important to consider the constraints regarding the conserved charges. Following the assumption of an isentropic evolution, one can employ conservation conditions like strangeness neutrality (Eq. (6)) and baryon density to charge density (Eq. (7)) to restrict the values of chemical potentials [16]. Here, is the net-charge to net-baryon number ratio of the colliding nuclei. For example, in Au + Au collisions, , with and denoting the number of protons and neutrons in the colliding nuclei. In a proton-proton collision, this ratio is .

The usual approach should be solving Eq. (4) to extract thermal parameters. The freeze-out description will be more reasonable if we include data for a larger number of detected particles in our solving mechanism. So a minimization is performed with all the available yields. One may avoid the volume systematics by taking ratios of two hadrons. In this approach, the effects of hydrodynamical flow also disappear [38]. Further, performing a minimization procedure with available yield ratios, one can parameterize the chemical freeze-out surface. In a recent work Ref. [28], we have shown that significant systematic uncertainty may arise in analysis due to variation in the chosen set of ratios.

##### 3.1. Our Approach

Following the complication regarding the chosen set of ratios, we have introduced an alternative method in Ref. [29]. The individual hadrons are not a conserved quantity in the strong interaction. So we opted to introduce ratios regarding conserved net charge densities like . Along with the constraints Eqs. (6) and (7), we have proposed two new independent equations, the net baryon number normalized to the total baryon number and the net baryon number normalized to the total hadron yield [58], as given below.

We want to mention that the left-hand side consists of the particle yields data from the heavy-ion collision, and those on the right are the number densities calculated from the thermal HRG model. runs only over detected () hadrons with available experimental yield data.

##### 3.2. Application to Centrality

The geometric information of the collision system is crucial to understand, as different final observables like eccentricity, elliptical flow, and charge particle multiplicity are directly dependent on the initial conditions like impact parameter (), the number of participating nucleons () [59]. We employ the centrality bins to differentiate collision events according to their impact parameters. As there is no direct method to measure , the centrality bins can be calibrated from charge particle multiplicity with the Glauber model [6, 59, 60]. Each centrality bin is represented by a corresponding . Most central ( centrality) collisions correspond to the events with the lowest value of impact parameters (highest value of ), whereas the most peripheral () are with the largest impact parameter and smallest . The degree of equilibration of the created medium should strongly depend on centralities as the system’s initial volume and initial energy, baryon deposition depends on these initial specifications.

##### 3.3. Strangeness with Centrality

Introducing a strangeness suppression factor is optional for the most central collisions [28, 29], whereas this appears essential when we deal with a peripheral or semicentral collision. Complete chemical equilibrium may not be achieved in the strangeness sector due to the higher mass threshold of strange particles and their hadronic counterpart [3]. Initially, this suppression factor was introduced considering the phase space undersaturation [22]. Ref. [23, 24] has discussed this undersaturation as an effect of the canonical ensemble consideration of strangeness, where exact strangeness conservation should be considered for a smaller collision system. In Ref. [61], a core corona-dependent model has also tried to discuss this suppression of strangeness. Irrespective of the reason for this undersaturation, considering this factor , gives rise to a better agreement to the thermal description of heavy-ion data. It seems that the strange sector may have a deviation from the respective grand canonical picture, and this parameter is a measure of that departure [13, 21–27]. In the presence of this factor, the number density is modified in the following manner [13]:

Here, denotes the number of valence strange quarks or antiquarks in the th hadron. In this work, we have calculated the number densities following 10. for all nonstrange particles. A smaller value of denotes a larger deviation from the grand canonical limit of equilibrium.

As we have introduced one added parameter , an extra equation is needed to close our system of equations. This parameter is not related to any conserved quantity, rather it is used to describe the possible nonequilibrium of the strange sector. Keeping in mind that we have used only yields of kaons among the strange particles, we have utilized kaon to pion ratio to evaluate the value of in Eq. (11).

Here, stands for two possible charges, i.e., and . Here, we want to reiterate that, for smaller system size (peripheral collisions), the exact strangeness conservation demands the canonical treatment. To study the systematic variation with centrality, we have approached within a Grand Canonical Ensemble (GCE) with the to scale the possible deviation from equilibrium picture. This exercise is well practiced in the context of freeze-out parameter extraction for various centrality [25–27, 62–65]. Finally, we have solved all these five equations Eq. (6)–(9) and Eq. (11) to extract the five parameters (T, , , , and ).

These three quantities are independent so as the constructed total charges. So the correlated uncertainties, arising from repeated entries of a single yield (addressed in Ref. [66]), are reduced in this formalism.

In this analysis, we have used yield data of , , and . In this context, can be constructed out of net-proton, whereas net charge is the sum of net-pion, net-kaon, and net-proton. This consideration is in line with the general approximation of taking the net proton as a proxy for net baryon number [67]. For these set of particles, the abovementioned equations will be

Here, the denotes the total number density of th particles, considering the relevant decay channels. We have considered all the strong decay channels from higher mass resonances, whereas weak decay corrections have been performed depending on the experimental specification [6, 66, 68]. In LHC, we have not included weak decay contribution into protons, whereas in RHIC energies, they are present in the total density.

#### 4. Results and Discussions

In this analysis, we have used the yields of (139.57 MeV), (493.68 MeV), and , (938.27 MeV) for collision energies () ranging from RHIC-BES (7.7 GeV) to LHC (2.76 TeV). For convenience, we have represented the centrality bins by their corresponding number of participants (). The collision system is Au-Au at higher RHIC and RHIC-BES energies and Pb-Pb at LHC. Data have been used following RHIC [6], RHIC-BES [66, 69], and LHC [68]. Data for - collision is available in RHIC for GeV and included in our analysis for completeness. In the present analysis, we have only taken midrapidity data. The details of the experimental yields used in the analysis are listed in Ref. [6, 66, 68, 69].

In our HRG spectrum, we have used all confirmed hadrons up to 2 GeV, with masses and branching ratios following the Particle Data Group [70] and THERMUS [71], which is a numerical thermal model package for the root framework. Finally, we solve Eqs. (6)–(9) and (11) numerically, using Broyden’s method with a minimum convergence criterion of [72]. We have estimated the variances of thermal parameters by repeating the analysis at the given extremum value of hadrons yields.

##### 4.1. Freeze-Out Parameters

We have described the variation of our extracted freeze-out parameters with the number of participants() for various collision energies in Figures 1(a), 1(b), 2(a), 2(b), and 3. In plots, the horizontal axis is the number of participants. Results for collision energy LHC-2.76 TeV to RHIC-BES-7.7 GeV have been shown in different columns in descending order from left to right. For the clarity of discussion, we shall discuss variation concerning first and then try to understand the changes with centrality. For completeness, we have also presented available results for , , and from other studies alongside our findings. We have included results from Ref. [6, 66] for RHIC and RHIC-BES and Ref. [73] for LHC.

**(a)**

**(b)**

**(a)**

**(b)**

The variation of the chemical freeze-out temperatures () in Figure 1(a) has good agreement with general understanding [29]. At most central collisions, the freeze-out temperature increases with collision energy, and near GeV saturates around the value of 160 MeV as it reaches the Hagedorn limit [74]. At LHC, is lower than the expected value as the proton yield is lower than their preceding RHIC energies [68]. The freeze-out temperatures seem to have a weaker dependence on the centrality and appear to be independent of and maintains a flat pattern at all collision energies. For GeV, the variation of is a little different at peripheral bins. The horn-like behavior is arising from a relatively lower yield of the proton (evident from the ratio and ), which also reflects in the extracted values of . There is good agreement with other results (black points) for . In Ref. [6, 66], the parametrization has been performed utilizing only pion, kaon, and proton. Results from these analyses match with our findings, whereas the little differences in LHC energy may arise from the fitting procedure and particle species used for fitting.

We have plotted the baryon chemical potential as a function of for all in Figure 1(b). The general expectation is that at lower collision energies, a larger amount of nucleons deposit in the collision region due to the baryon stopping [75, 76]. But at very high collision energies, the nuclei are transparent to each other [77, 78]. So at higher RHIC and LHC energies, the medium is created having almost zero net baryon number. Therefore, the net baryon density and hence the estimated chemical potential would decrease with increasing due to baryon transparency [78]. In the same manner, one should expect a rise of for higher . In central collisions, the value of the deposited net baryon number increases due to baryon stopping among a large number of the participating nucleons. Contrarily for a peripheral collision, a lesser number of nucleons get deposited in the collision zone, creating a dilute system of net baryon, which results in a smaller value of . We have observed this trend in all . Our resulted agrees with previous findings from Ref. [6, 66].

Strangeness chemical potential shows a similar trend as in Figure 2(a). It decreases as collision energy increases and becomes zero at LHC energy. On the other hand, escalates as one goes from peripheral to the central collision. The correlation between and can be described in the following manner. A higher baryon density demands hyperons to be produced more than antihyperons. To maintain the strange neutrality, this excess amount of strangeness from the baryonic sector has to be nullified from the mesonic sector. So in the mesonic sector, is more abundant than . Being the lightest strange particles, this difference between charged kaons determines the sign and trend of .

The general trend of charge chemical potential is the same as other s except for the sign. As increases, it becomes more negative and the magnitude decreases with in Figure 2(b). However, is more negative for larger baryon densities. We can understand this as following. The neutrons are more abundant than protons in the colliding heavy nuclei. This abundance generates a net negative isospin value in the collision system and produces more than , to conserve the isospin. As the lightest charged particle, these pions determine the negative . This reasoning will be more clear if we look into the value of for GeV at . In this case of - collision, the isospin dominance should not act in favor of , as there is no neutron in the colliding particles. So one should expect the to be positive for this case. Indeed, we have observed a positive value of for the - collision of 200 GeV RHIC energy. The net value of isospin increases with the , thus increases asymmetry between the yield of charged pions. So the magnitude of rises following the .

If strange particles achieve chemical equilibrium, then the thermal abundance of kaons should be described by equilibrium thermal parameters (T, s) of a grand canonical ensemble. This is not observed in cases of a small collision system like -, -, and even in - with a smaller [13]. Several models [23–25, 62, 64] have tried to describe this source of strangeness undersaturation in the smaller system and advocated the use of . The common perception from all this work is that scales the deviation of strange particles from their respective equilibrium thermal yield of a grand canonical ensemble, while denotes the equilibration in the strange sector.

We have shown the variation of in Figure 3. It is interesting to notice that even in LHC and high RHIC energies, the has an increasing trend from lower peripheral to a central collision, though the temperature and other chemical potentials do not change much. Initially, it starts from a lower value in the case of peripheral collisions and increases with participants. Around , tends to saturate to the most central values. It appears that the strangeness tends to be closer to the grand canonical limit as the system size increases. The saturation of with the colliding system size for central - collisions suggests that the strangeness suppression may be independent of the hadronic scatterings, which happens in the later time of the evolution [62]. From the pattern, one can also conclude that the strangeness equilibration has a prominent dependence on the number of participants and system volume. The differences between the values between peripheral and central are larger for lower RHIC-BES energies. The general understanding from the above study is that the strangeness sector may be closer to equilibrium in a peripheral collision of higher , whereas the deviation from equilibrium is larger for peripheral cases in lower collision energy. We have found that strangeness equilibration tends to happen in collisions with a more than , which may be the threshold for the creation of a deconfined phase in a - collision, which drives the system close to strangeness equilibration [3].

We want to mention that even at high RHIC and LHC energies, the central value of lies below (around 0.9). This finding is in agreement with previous analyses from RHIC collaboration [6, 66]. In Ref. [66], the is shown to increase and saturate near as more hyperon species are included in chemical freeze-out parametrization. In this context, our method would be similar to the standard chi-square analysis, where the parameterization depends on the chosen hadronic ratios. Future studies with other heavy ions and data of hyperons may help to understand this.

##### 4.2. Scaling Nature of CFO Parameters

The scaling behavior with is important to calibrate the chemical composition at freeze-out with system size. To simplify the discussion, we have normalized the obtained parameters by their corresponding value for the most central collision of individual . As an example, to understand the scaling of temperature at 200 GeV, we have divided the extracted of each centrality bin () with that of the most central collision (maximum ). Figure 4 shows the variation for all five freeze-out parameters. For simplicity, we have plotted for two collision energies from both RHIC (200, 62.4) and RHIC-BES (27, 11.5). These scaled quantities should lie around if parameters do not vary much with centrality (). The scaled freeze-out temperature shows this pattern for all collision energies. It seems that for given incident energy, the freeze-out temperature does not vary much with the system size, whereas the scaled baryon chemical potential () has an increasing trend as it becomes maximum at most central collisions. In the case of equilibrium among all the charges, all the chemical potentials should commensurate with each other. Scaled and should follow the pattern of with both the number of participants and collision energy, which we have already discussed in Section 4.1. We have indeed observed a similar trend for all three s. There is a trend of saturation near around . Onward this point, the system may achieve a thermodynamical state which is independent of the system size. Future analysis with other colliding ions at these c.m energies may shed light on this issue. Nontriviality could have appeared in the case of as it is a nonequilibrium parameter. But the observed trend is similar to the chemical potentials. It starts from a smaller magnitude and saturates onward . The system may have enough energy and number density for strangeness equilibration onward this centrality bin [79], and we can employ a grand canonical description to describe the yield at freeze-out. The deviation of scaled from central value is larger for peripheral collisions in lower RHIC-BES energy, which indicates that colliding energy has a crucial contribution to decide the equilibration of strangeness.

##### 4.3. Particle Yield Ratios

In this section, we shall discuss ratios regarding detected particles to check the efficiency of our parameterization. We have estimated particle ratios from our extracted freeze-out parameters and plotted them alongside their experimental values. Variances in the detected yield ratios have been obtained using the standard error propagation method [80], considering both the systematic and statistical uncertainties of data. We have calculated the variance of thermally estimated ratios by evaluating them at the extrema of the obtained freeze-out parameters.

In Figure 5(a), we have plotted the particle to antiparticle ratios for pions and kaons. There is good agreement between model estimation and experimental data for both ratios. No notable variation has been observed for with and . Freeze-out temperature and determine the chemical abundance of pions. For the ratio of negatively charged to positively charged pions, the variation should depend on only. Here, we want to mention that there is no prominent variation of with and , as the value of is much smaller (around 5 MeV) than the value of (about 150 MeV). On the other hand, is much lower than the mass of the pion itself. So it does not differentiate between the thermal yield of and , and the ratio lies near unity for all collision energies and centrality classes.

**(a)**

**(b)**

The asymmetry between and depends on net baryon density and net strangeness from the hyperon sector. One should expect a larger yield of than at higher baryon density (). We have observed this pattern with both centrality and . At lower RHIC-BES energies, the ratio is far from unity due to higher net baryon density and approaches as the collision energy increases. The yields of particle and antiparticle become equal at LHC as baryon transparency takes over, and the medium starts with zero net baryon density. With , a commensurable trend has been observed following the value of . At peripheral bins of lower collision energies, our model has overestimated the ratio. This overestimation has occurred as an interplay between the ratio and the constraint . The decreases towards peripheral and noncentral collision following the ratio, which results in a lower density of and other hyperons. On the other hand, to maintain the strangeness neutrality, thermal parameters adjust to produce a larger thermal density of , which results in this overestimation.

The charge independent ratios , are important observables for understanding the strangeness production in high-energy collisions. As the lightest mass hadrons, pions may act as the proxy of entropy, whereas the kaons carry the signature of strangeness. Strange particles are important in studying chemical equilibrium in heavy-ion collisions due to their late production [3]. A charged particle multiplicity () dependent saturation of strangeness normalized to pions has already been observed in LHC [5], which can be utilized to investigate the system size-dependent strangeness production. Ref. [79] has related this saturation to equilibration with a threshold . In heavy nuclei collisions, the overlap region and both are related to [6]. We have observed the same saturation trend with here for all in Figure 5(b) and have suitably reproduced it with our parametrization. This saturation starts around in higher RHIC and LHC. Here, we want to mention that there is no variation of with centrality and collision energy, but has a strict dependence on both. At lower collision energy, is much higher than due to the excess yield of . The difference between the ratios decreases with increasing collision energies, and they become equal at LHC, as the particle-antiparticle yields become the same. The pattern of has a close resemblance to both the ratios. It seems that as decreases, the kaon yields deviate far from their equilibrium yield. So a nonequilibrium parameter had to be introduced in our thermal model. Lower the value of , larger is the deviation from equilibrium for kaons.

Here, we want to reemphasize that both and have no significant variation with centrality bins at LHC. This symmetry between particle and anti-particle demands and to be almost zero. A centrality variation in the ratio cannot be reproduced with zero without introducing a like parameters. This centrality variation of the ratio indicates that the strange particles are out of equilibrium at peripheral collisions of LHC.

From the discussion of freeze-out parameters, it appears that the variation of antiproton (anti baryon) to proton (baryon) is a guideline to understand the variation of . This ratio becomes at upper RHIC and LHC energies as the colliding nuclei pass through each other, and the hadrons are created out of a medium having zero net baryon density. In lower collision energy, baryon stopping motivates a larger net baryon density. As a result, protons are more abundant than antiproton and advocate a smaller at lower collision energy. Our thermal model estimations have good agreement with the experimental data in Figure 6(a). Towards peripheral collisions, this ratio tends to increase and symbolizes the decrease of baryon dominance over antibaryons. Initial net baryon number density decreases as one goes from central to peripheral collision due to nuclear distribution [81] of the colliding nuclei and induces a smaller baryon antibaryon asymmetry in their yield.

**(a)**

**(b)**

At this point, we also want to mention that for the ratio , there is a chance of overfitting as our constructed ratio net-baryon to total baryon Eq. (8) reduces to , as we have utilized centrality dependent data for (anti-)proton from the baryon sector. This deficiency of our formalism will reduce when centrality data for other (anti-)baryons are considered (as discussed in Ref. [29]). On that occasion, will be an independent prediction.

We have discussed proton to positively charged pion ratio in Figure 6(b). As we have already discussed, pions may act as the measure of entropy. So the ratio will describe the variation of baryon production with entropy. If the particles are produced only from deposited energy, then pions will be highly abundant than massive protons. But if the medium starts to evolve from a finite baryon density, then per pion, proton production will be larger to conserve the net baryon density. That is why a clear increasing trend for takes place as the collision energy decreases. This same increment is expected with , as more baryon deposition happens in the more central collisions. This variation is prominent in lower due to the higher efficacy of baryon stopping.

#### 5. Summary and Conclusion

The equilibration of the system created in a heavy-ion collision should strongly depend on the system size and number of participants. A comparison among the chemical freeze-out conditions of high and low multiplicity - collisions may shed light in that direction. Instead of the general minimization, a conserved charge-dependent parametrization process has been adopted, utilizing the midrapidity yield of the pion, kaon, and proton to explore the freeze-out parameters of various centrality bins of - collision for LHC (2.76 TeV), RHIC (200 GeV, 130 GeV, 62.4 GeV), and RHIC-BES (39 GeV, 27 GeV, 19.6 GeV, 14.5 GeV, 11.5 GeV, and 7.7 GeV). We have incorporated a strangeness suppression factor () to estimate the possible nonequilibrium in strange hadrons in the peripheral collisions. We have discussed the variation of these chemical freeze-out parameters with both centrality and collision energy.

The variation of parameters with collision energies has good agreement with general understanding, whereas there are significant variations with the number of participants. The extracted freeze-out temperature has no strong dependence on the number of participants (centrality), whereas the chemical potentials show a wide variation. We have presented the behavior of scaled parameters to have a better understanding of the centrality variation. These parameters have been normalized with the value obtained in most central collisions and compared along with the other collision energies. Scaled and appear to follow the scaling behavior of , which may be a signal of equilibration among three conserved charges. The strangeness suppression factor () deviates from the equilibrium value at peripheral collisions and tends to saturate near unity in central collisions. In the peripheral collisions, starts around .7 and increases towards most central bins. This variation indicates that the strange hadrons are deviated from equilibrium at low multiplicity collision, whereas there is a sign of equilibration as the increases. The flattening of the scaled parameter and appears around a threshold of . So we can apply a grand canonical description for the systems created out of - collisions with more participants than . We have found the to lie below (around ) even at most central collisions. In this study, we have only used kaons, so the variation of is an artifact of the kaon to pion ratio. Future analysis with yields data of other strange hyperons may help to understand this further.

Further, we have estimated different particle ratios to cross-check the effectiveness of our parameterization. Our estimated hadron ratios seem to have good agreement with experimental data. We have only reproduced ratios regarding pions, kaons, and protons as they are present in our analysis. A saturating trend with has been observed for the kaon to pion ratio and explained with the .

We want to mention that the centrality variation has previously been investigated in RHIC-BES energy [66] and in LHC [56, 57, 73] with the approach. Rather than the conventional practice, we have followed a fitting procedure that relies on the conserved quantities and produces similar parameter sets. The agreement with other studies will act as a benchmark for the future application of this parameterization. Further, we have found a threshold , which is significant to study the bulk properties in a thermodynamic picture.

#### Data Availability

All data used in the analysis are properly cited.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work is funded by UGC and DST. The author thanks Sumana Bhattacharyya, Sanjay K. Ghosh, and Rajarshi Ray for the helpful discussions, and Pratik Goshal and Pracheta Singha for the critical reading of the manuscript. The publication and article processing charges are funded by SCOAP3.