#### Abstract

We present a subset of experimental results on charge fluctuation from the heavy-ion collisions to search for phase transition and location of critical point in the QCD phase diagram. Measurements from the heavy-ion experiments at the SPS and RHIC energies observe that total charge fluctuations increase from central to peripheral collisions. The net-charge fluctuations in terms of dynamical fluctuation measure are studied as a function of collision energy () and centrality of the collisions. The product of and shows a monotonic decrease with collision energies, which indicates that at LHC energy the fluctuations have their origin in the QGP phase. The fluctuations in terms of higher moments of net-proton, net-electric charge, and net-kaon have been measured for various . Deviations are observed in both and for net-proton multiplicity distributions from the Skellam and hadron resonance gas model for GeV. Higher moment results of the net-electric charge and net-kaon do not observe any significant nonmonotonic behavior as a function of collision energy. We also discuss the extraction of the freeze-out parameters using particle ratios and experimentally measured higher moments of net-charge fluctuations. The extracted freeze-out parameters from experimentally measured moments and lattice calculations are found to be in agreement with the results obtained from the fit of particle ratios to the thermal model calculations.

#### 1. Introduction

The main goal of the high energy heavy-ion collisions is to study the phase structure of the quantum chromodynamic (QCD) phase diagram at finite temperature () and baryon chemical potential () [1–5]. Several theoretical models suggest that the QCD phase diagram may contain a first-order phase transition line between the hadron gas (HG) phase and Quark-Gluon phase which ends at the critical point towards high and lower [6–10]. Experimental programs have been performed at SPS and beam-energy scan (BES-I) program at RHIC to search for critical point and QGP-HG phase transition. In future, the upcoming program at RHIC (BES-II) [11], FAIR [12], NICA [13], and J-PARC [14] will also contribute to the physics at large . The location of the QCD critical point can be explored by systematically varying and . Experimentally, by changing the center of mass energy one can control and of the system and hence enables us to scan different sectors of the phase diagram.

One of the most striking signatures of such a QGP-HG phase transition could be a strong modification in the fluctuations of specific observables measured on an event-by-event basis [15, 16]. In principle, any observable that is not globally conserved fluctuates. Although most of these fluctuations are trivial and are of statistical origin, it is important to find out the dynamically relevant event-to-event fluctuation that enables the search for a possible critical point and a first-order coexistence region in the QCD phase diagram [17]. Over the past two decades quite a number of such observables have been suggested for clarifying the evolution of the system formed in heavy-ion collisions. These either refer to the signals from the plasma that are supposed to survive the phase transition or to the observables that experience strong fluctuation during the phase transition or close to the critical point. Most commonly measured event-by-event fluctuations in heavy-ion collision experiments are particle ratios (, , etc.), transverse energy , mean transverse momentum , and particle multiplicity fluctuations [18–21]. Predictions suggest that enhanced multiplicity fluctuations are connected to the production of QGP droplets, and suppression of fluctuation is connected to the nucleation process in a first-order QGP-HG phase transition. This may happen because of the rapid freeze-out just after the phase transition [15, 22]. An isothermal compressibility of the system can be considered to understand the sensitivity of the measured particle multiplicity to the phase transition [23]. The isothermal compressibility is defined as , where , and are the pressure, volume, and temperature of the system, respectively. In order to relate the compressibility to the measurements of multiplicity fluctuations, we assume that relativistic heavy-ion collisions can be described as a thermal system in the Grand Canonical Ensemble (GCE) [24]. The GCE is the most appropriate description as only part of the particles from the system around mid-rapidity are measured by the experiments. The energy and conserved quantum numbers in this region can be exchanged with the rest of the system that serves as a heat bath [25]. Several other studies have applied Canonical and Microcanonical ensembles to the multiplicity fluctuations too [26–28]. In the GCE, the isothermal compressibility is directly related to the variance of the particle multiplicity as follows:where is the particle multiplicity, is the mean of the multiplicity distribution, and is the Boltzmann constant [29]. Here, multiplicity fluctuation measurements are presented in terms of the scaled variance, , as [30]

In a continuous or second-order phase transition, the compressibility diverges at the critical point. Near the critical point, this divergence is described by a power law in the variable , where is the critical temperature. Hence, the relationship between multiplicity fluctuations and the compressibility can be exploited to search for a signature of critical behavior by looking for the expected power law scaling of the compressibility:where is the critical exponent for isothermal compressibility [29]. Recent studies [31, 32] show the behavior of the quark number susceptibility, , which is related to the value of the isothermal compressibility of the system. They predict that its value will increase by at least an order of magnitude close to the QCD critical point. As discussed above, the scaled variance is proportional to ; hence, the measurements of charged particle multiplicity are expected to be a sensitive probe for critical behavior.

If the system approaches close enough to the critical line for a long enough time period, then critical phenomena could be observed through the measurement of multiplicity fluctuations [8]. Subsequently, it may also be possible to determine the critical exponents of the system. Observation of critical behavior in heavy-ion collisions and the subsequent measurement of the critical exponents could determine the universality class in which QCD is grouped, providing essential constraints for the models [30].

The fluctuations of conserved quantities are predicted to be one of the most sensitive signals of the QGP formation and phase transition, which may provide complementary understanding of strong interactions, apart from other QGP signatures [15, 33]. It has been argued that entropy conserving hadronization of plasma of quarks and gluons should produce a final state characterized by a dramatic reduction of the net-charge fluctuations in QGP phase as compared to that of a hadron gas. Further, prediction relies on the notion that quark-quark correlations can be neglected, and hadronization of gluons produces pairs of positive and negative particles not contributing to the net-charge fluctuations. It has also been suggested that the excitation function of conserved numbers like net-baryon, net-charge, and net-strangeness fluctuations should show a nonmonotonic behavior, as a possible signature of QCD critical end point (CEP) [1, 34, 35]. In the thermodynamic limit, the correlation length () diverges at CEP [1]. The experimentally measured moments of the net-baryon, net-charge, and net-strangeness distributions are related to the higher power of of the system and hence these moments can be used to look for signals of a phase transition and critical point [36, 37]. Also, the comparison of experimentally measured cumulants with the lattice calculations enables us to extract the freeze-out parameters, that is, freeze-out temperature () and of the system produced in heavy-ion collisions [38, 39]. In recent years, lots of efforts have been put on both theoretical and experimental fronts to study the fluctuation of conserved quantities.

This review is organized as follows. In the following section, we discuss the total charge fluctuations from various experiments. In Section 3, the results on net-charge fluctuation are presented, which include dynamical fluctuation measure , higher moments of net-proton, net-electric charge, and net-kaon fluctuations. Towards end of Section 3, extraction of freeze-out parameters using higher moments is discussed. Finally, in Section 4, we summarize our observations.

#### 2. Total Charge Fluctuation

In a thermodynamical system of strongly interacting matter formed in the heavy-ion collisions, the fluctuations of particle multiplicities, mean transverse momentum (), transverse energy (), and other global observables are related to the fundamental properties of the system, such as specific heat, chemical potential, and compressibility. These observables either refer to signals from the plasma that are supposed to survive the phase transition or to observables that experience strong fluctuations during the phase transition or close to the critical point. The existence of a critical point at the QCD phase transition has been associated with the large event-by-event fluctuations of above observables. As far as observables are concerned, electric charge fluctuations have been measured over a wide range of collision energies, from CERN SPS to RHIC and LHC energies. Enhanced fluctuations in neutral to charged pions have been predicted as a signature of the formation of Disoriented Chiral Condensates (DCC) [40, 41]. The relative fluctuation which can be extracted from experimental data has contributions from both statistical and dynamical sources. In order to extract the dynamical part associated with new physics from the observed fluctuations, one has to understand the contributions from statistical and other known sources. Some of the known sources of fluctuations contributing to the observed experimental value of scaled variance () include finite particle multiplicity, effect of limited acceptance of the detectors, impact parameter fluctuations, fluctuations in the number of primary collisions, effects of rescattering of secondaries, resonance decays, and Bose-Einstein correlations [42].

The relative fluctuation is defined as where is the charged particle multiplicity and is also known as scaled variance. If the multiplicity distribution is Poissonian, the scaled variance is 1.0. Figure 1 shows the comparison of the relative fluctuation of the charge particle multiplicity as a function of collision centrality which is related to number of participants () in Pb + Pb collisions at 158 AGeV [40]. The error on calculated in the model is mainly due to the error on the mean number of charged particles in nucleon-nucleon interactions, the error in the number of participants calculated, and the uncertainty in the calculated transverse energy [40]. The experimental data are compared with the model calculations. It is observed that the relative fluctuations increase from central to peripheral collisions. The observed charge particles multiplicity fluctuations have been found to well agree with the results obtained from a simple participant model [45]. In the participant model, the particle multiplicity may be expressed as where is the number of participants in the collision and is the number of particles produced by the th participant within the detector acceptance. The mean value of is the ratio of the average multiplicity measured in the detector acceptance to the average number of participants, . Hence, fluctuation in the particle multiplicity will have contributions due to fluctuations in , () and also due to the fluctuations in the number of particles produced per participant (). The multiplicity fluctuation in the participant model can be expressed as [40] Another experiment at SPS also performed similar study shown in Figure 2 [43]. The scaled variance (), where is the variance of the distribution and is the multiplicity of the particles. The scaled variance of positive, negative, and total charged particles as a function of centrality is shown in Figure 2. The experimental data is compared with the model calculations. The results from different models (HIJING [46], HSD [47], UrQMD [48], and VENUS [49]) are almost independent of centrality and behave like a Poisson expectation. However, the experimental data points indicate strong dependence on centrality. The scaled variance increases from central to peripheral collisions. The measured centrality dependence can be reproduced in superposition models with the assumption of contributions from target participants to the particle production in the forward hemisphere [43, 44]. Figure 3 shows the centrality () dependence of scaled variance at = 62.4 and 200 GeV in Au + Au collisions at RHIC [30]. The shaded regions represent the systematic uncertainties from the reference range. The statistical uncertainties are shown along with the data points. Here, represents the estimate of the remaining dynamical multiplicity fluctuations. For all centralities, the scaled variance lies above the Poisson expectation of 1.0. At these energies, also the scaled variance increases from central to peripheral collisions. Hence, similar centrality dependence has been observed by the experiments at the SPS and RHIC energies. The absence of large dynamical fluctuation in excess of the participant superposition model indicates that there is no evidence of critical behavior related to the compressibility observable.

**(a)**

**(b)**

**(c)**

#### 3. Net-Charge Fluctuations

One of the proposed signatures to search for the phase transition from hadronic to partonic medium is to study the net-charge fluctuations in heavy-ion collisions. The fluctuation in the net-charge depends on the square of the charges and hence strongly depend on which phase it originates from. The charge carriers in the QGP phase are quarks having fractional charges, while in hadronic phase the constituents have unit charge; hence, the measure of the fluctuations in the net-charge is expected to be different in these two cases [22]. In this section, we discuss the net-charge fluctuations using different fluctuation measures such as and higher moments.

##### 3.1. Fluctuation Study Using

The net-charge fluctuations are expected to be smaller in the QGP phase as compared to the hadron gas (HG) phase [22]. The net-charge fluctuations may get affected by uncertainties due to volume fluctuation, exact local charge conservation, or repulsive forces among hadrons [39]. However, it is important to know whether these fluctuations may or may not survive the evolution of the system in the heavy-ion collisions. The collision volume is not directly measured in the experiment which may lead to additional geometrical fluctuations. One can get rid of volume fluctuation by considering the ratios of the number of positive () to negative () particles, . The ratio fluctuation is defined by . Since the fluctuation of the number of charged particles is a measure of the entropy of the system, another fluctuation observable -measure of the net-charge provides a measure of charge fluctuations per unit entropy and is related to the ratio as where is the variance of the net-charge with being the difference between +ve and −ve particles and is the average number of charged particles measured within the experimental acceptance. Assuming the quark-quark interactions to be negligible, was found to be approximately 4 times smaller in QGP phase as compared to HG phase. It has been shown that for uncorrelated pion gas and reduces to 3 after taking the resonance decay into account [25]. For a QGP phase, is estimated between 1.0 and 1.5. Hence, can be used as a probe to distinguish between the QGP and HG phases. Unfortunately, the quantity depends on the experimental efficiency. In the experiments, the net-charge fluctuations are studied in terms of dynamical fluctuation measure , which is found to be independent of detection efficiency. The quantity is a measure of the relative correlation strength of , , and particle pairs. A positive value of signifies the correlation of same charge pair, where a negative value indicates the dominant contribution from correlations between opposite charges. is defined by where and are the average number of positively and negatively charged particles within the detector acceptance. The -measure and are related as [50].

Figure 4 shows the centrality () dependence of dynamical net-charge fluctuations () in Au + Au and Cu + Cu collisions at different . values exhibit a monotonic dependence on and have small dependence on collision energy. For all the studied energies, the values of are negative, which indicates the dominance of correlation of positive and negative charged particle term in (8). The observed monotonic reduction of the magnitude of with increasing number of participants arises mainly due to the progressive dilution of two-particle correlation when the number of particle sources is increased.

In view of prediction for critical point in the QCD phase diagram in the range GeV [8, 51], it can be argued that the reduction of fluctuation might be larger at lower beam energies. On the other hand, one may also argue that the volume of the QGP formed in Au + Au collisions might increase at higher beam energies leading to reduced fluctuations. Hence, it is important to understand whether the fluctuations may be found to vary with beam energy thereby indicating the production of QGP above a critical threshold or with progressively increasing probability at higher energies. Figure 5 shows the product of and (average number of charged particles) as a function of collision energies for 0–5% central collisions using the combined data from SPS, RHIC, and LHC energies [51–53]. Also, the collision energy dependence of -measure is shown in the same Figure 5. It is observed that the fluctuation observable shows monotonic decrease in magnitude with increasing and approaches expectation for a simple QGP-like scenario [22] as we move from RHIC to LHC energies. It has been argued that measurements in lower (below 10 GeV) are dominated by baryons while at higher energies the meson and resonance production becomes increasingly dominant. This suggests that the change in dynamical net-charge fluctuations below GeV might be partially due to this shift in particle production dominance [1, 8, 51, 54–56]. It is also argued that the differences between the fluctuation values below and above 19.6 GeV may result from changes in the collision dynamics and final state interaction effects [8, 54–57]. For the highest RHIC energy, the measured value of fluctuation observable is close to the HG prediction, whereas at lower energies the results are higher than HG value. This may be due to the fact that at highest RHIC energy ( GeV) the fluctuation may not be strong enough to be measured or because of the dilution of fluctuation during the evolution process. At LHC energy TeV, the fluctuation observable value is significantly lower as compared to lower energies results. The fluctuations at the LHC energy might also have been diluted because of various effects; still these fluctuations are smaller than the theoretical expectations. In [22], it is shown that the -measure value for hadron gas with resonance decay is 3 and for QGP phase is 1.0–1.5. The measured -measure value is for at TeV. There is a clear decreasing tendency of the -measure value in the HG phase and approaches toward QGP expectation. This may indicate that the fluctuations have their origin in the QGP phase [53].

Several other observables already indicate that a hot and dense medium of color charges has been formed at RHIC and LHC energies. The net-charge fluctuation results suggest that either the observable is not sensitive enough to QGP physics or the process of hadronization washed out the QGP signal for this observable. It may be also noted that the theoretical results do not incorporate the acceptance effects and dynamical evolution of the system like, for example, the dilution of the signals in the hadronization process [58].

##### 3.2. Fluctuation Study Using Higher Moments

In recent years, the beam-energy scan (BES) programs at SPS and RHIC have drawn much attention to map the quantum chromodynamics (QCD) phase diagram in terms of temperature () and baryon chemical potential (). The location of the critical point can be explored by systematically varying and , which can be experimentally achieved by varying the of the colliding ions. Several theoretical models suggest that excitation function of conserved numbers such as net-baryon, net-charge, and net-strangeness fluctuations should show a nonmonotonic behavior as a possible signature of QCD critical end point (CEP) [6–8]. In the vicinity of the QCD critical point, these variances are proportional to the square of the correlation length (), which is expected to diverge at the critical point [34, 59]. However, the magnitude of the correlation length is limited by the system size and by finite time effects (critical slowing down), which could be as small as 2 to 3 fm. Hence, the contribution to the fluctuations from the critical point might be too weak as to be discovered experimentally, if only the second moments are measured. Therefore, it has been proposed to measure higher moments of the fluctuations which are expected to be much more sensitive to the critical point [35, 60, 61]. The moments of the net-baryon, net-charge, and net-strangeness distributions are related to of the system and hence these moments can be used to look for signals of a phase transition and critical point [1, 36, 37]. The variance () of these distributions is related to as , the skewness () goes as , and the kurtosis () is related as . Hence, it is proposed to study the higher moments of conserved number distribution due to their stronger dependence on .

Experimentally, the net-baryon number fluctuations are not directly measured, as all neutral baryons are not detected by most of the experiments. Hence, net-baryon fluctuations are accessible via measuring the net-proton distributions [62, 63]. The net-charge fluctuations are accessible by measuring the stable charged particles such as pions, kaons, and protons along with their antiparticles [64, 65]. Similarly, the measurement of net-kaon fluctuations acts as a proxy for net-strangeness fluctuations, because higher mass strange particles are not directly measured [66, 67]. In the following subsections we discuss each of the conserved number fluctuations separately.

The moments of the event-by-event experimentally measured net distributions are related to the different order of the cumulants of the distribution as mean () = ; ; , and , where is the multiplicity of the distribution and = . Hence, the ratios of the cumulants are related to the moments as , , , and . Further, the ratios of moments can be related to the susceptibilities of th order () obtained from the lattice QCD or from the HRG model calculations as , , , and [68]. One advantage of measuring the ratios is that the volume dependence (which is not directly measured by the experiment) on the experimentally measured individual cumulants cancels out to first order in the ratios. Hence, experimentally measured quantities can be directly compared with the theoretical calculations [68, 69].

###### 3.2.1. Net-Proton Fluctuation

Theoretical calculations have shown that net-proton fluctuations reflect the singularity of the charge and baryon number susceptibility as expected at the critical point [61]. Figure 6 shows the typical uncorrected net-proton () distributions in Au + Au collisions for 0–5% centrality measured on event-by-event basis. The protons and antiprotons are measured within = 0.4 to 0.8 GeV/ and 0.5. At lower energies, the net-proton distributions are dominated by the contributions from the proton distributions. The mean value of the distribution increases with decrease in collision energy. Different moments which describe the shape of the distribution are extracted from the distribution.

Figure 7 shows the collision energy dependence of and of net-proton distribution for 0–5% and 70–80% centralities in Au + Au collisions. The statistical uncertainties are calculated using Delta theorem approach [70]. In the hot and dense medium, the baryon chemical potential decreases with increasing collision energies; hence, Figure 7 can be interpreted as dependence of the moments over the large range of (20–450 MeV) [63]. Deviations are observed for both and from the Skellam and hadron resonance gas model for GeV. Maximum deviations from Skellam expectation are observed for and 27 GeV. The experimental results are reasonably described by assuming independent production (IP) of and indicating there is no apparent correlations between the proton and antiprotons for presented observables [63]. Therefore, one may ask, in spite of significantly correlated production due to baryon number, electric charge conservation, and kinematical correlations of proton and antiprotons, why do the measured cumulants follow the independent production model. This has been studied by introducing the correlation between two independently produced distributions. It is observed that experimentally measured cumulants will follow the IP model calculations even if the correlation coefficient is less than 20% [71]. However, = and = values will follow the IP model for all the correlation coefficient values. The observation that the experimental data can be explained by the independent production of particles does not rule out the existence of the critical endpoint.

**(a)**

**(b)**

**(c)**

###### 3.2.2. Net-Charge Fluctuation

The net-electric charge fluctuations are accessible through measuring fluctuations of stable charged particles (, , and ). Net-electric charge fluctuations are more straightforward to measure experimentally than net-baryon number fluctuations. As discussed before, experimentally, net-baryon number fluctuations are accessible only through net-proton number fluctuations, but in lattice QCD calculations net-baryon fluctuations are calculated instead of net-proton fluctuations. While net-charge fluctuations are not as sensitive as net-baryon fluctuations to the theoretical parameters, both measurements are necessary for a full understanding of the theory [65]. Several studies suggest that the net-charge multiplicity distributions are better suited to extract the freeze-out parameters ( and ) and the location of the CEP as they directly probe a conserved quantum number (electric charge) [37, 73, 74].

Figure 8 shows the collision energy dependence of efficiency corrected and of the net-charge distributions for central 0–5% Au + Au collisions at RHIC. and measurements from the PHENIX experiment are shown in (a) [65] and measurements from the STAR experiment are shown in (b) [64]. The statistical errors are calculated using Delta theorem method [70]. Here, we would like to briefly discuss about the statistical error calculations in both the STAR and PHENIX experiments. The errors on and are correlated and are dependent on both variance () of the distribution and the number of events. As mentioned in [70], the statistical errors are more dependent on than the number of events. Hence, experiments having larger multiplicity distribution will have larger errors on and . Since STAR experiment has larger acceptance compared to PHENIX, it contributes to larger statistical error. Further, STAR net-charge distributions have larger compared to net-proton distributions, which contributes to the larger statistical error in the net-charge results than the net-proton results, although the number of analyzed events are similar. values from both the experiments decrease with increasing . values from PHENIX ((a), bottom panel) remain constant and positive, between at all the collision energies within the statistical and systematic uncertainties. However, there is 25% increase of values at lower energies below GeV compared to higher energies. Further, the values from the STAR experiment are constant at all energies within uncertainties, except for GeV which shows a negative value. The STAR experiment reported higher weighted mean () of values as compared to PHENIX values. It is to be noted that results from the PHENIX experiment are measured within and with in azimuth, whereas results from the STAR experiment are within and with full azimuth (2). Different lower cut may be responsible for about 30% of the difference between two data sets [75]. More discussion on acceptance effect on higher moments can be found in [72, 75, 76]. However, net-charge results from both the experiments do not observe any significant nonmonotonic behavior in the products of moments as a function of collision energy.

**(a)**

**(b)**

###### 3.2.3. Net-Kaon Fluctuation

Experimentally, the net-strangeness fluctuations are accessible through measuring the net-kaon fluctuations. Figure 9 shows the preliminary results on collision energy dependence of efficiency corrected and of the net-kaon distributions for 0–5% and 70–80% centralities in Au + Au collisions measured by STAR experiment at RHIC [67]. With large uncertainties in the measurement, no significant deviation of the product of higher moments for net-kaon distributions as compared to the Poisson expectation has been observed at the measured energies. However, the upcoming RHIC BES-II with the upgraded detector system will help to reduce the uncertainties on the measurements and may find the location of the critical point.

##### 3.3. Extraction of Freeze-Out Parameters Using Higher Moments

Product of higher moments can be used to extract the freeze-out parameters ( and ) of the QCD phase diagram [37–39]. From the lattice calculations, it has been observed that the ratio of 1st to 2nd cumulants (=) shows a strong dependence on but varies little with . On the other hand, the ratio of 3rd to 1st cumulants (=) shows strong dependence on and has little dependence on [37, 38]. Hence, lattice calculations in combination with experimentally measured values at different collision energies can be used to extract , whereas the measured can be used to extract . The collision energy dependence of freeze-out parameters and is shown in Figure 10. Freeze-out parameters extracted from different methods such as using particle ratios [77] and using cumulants with lattice [65] agree very well. Further, extracted and values using lattice calculations and experimental data [65] are in agreement with the thermal model parameterization [77]. The freeze-out parameters extracted using combination of experimental data and HRG model [39] are also shown in Figure 10. The HRG calculations are performed in the same acceptance as the experiment. The extracted values using HRG and experimental measured cumulants are in agreement with other measurements, but the extracted are about 7–10 MeV lower as compared to the values extracted using lattice calculations [39, 65]. values are consistent across different extraction procedure but there is some discrepancy in the value of . One of the reasons may be that are extracted using ratio of 1st and 2nd cumulants which are measured more precisely compared to higher cumulants. In [39], and are extracted using HRG for each by simultaneously reproducing the experimentally measured ratios of the lowest-order susceptibilities () for net-protons and net-electric charge. In [38], the freeze-out parameters are extracted using STAR experimental data [64]. has been extracted using experimentally measured and lattice calculations, whereas due to the uncertainties on the lattice results in the low-temperature region, it is only possible to extract an upper value for the freeze-out temperature (MeV). In [65], the freeze-out parameters are extracted using and of net-charge from the PHENIX experiment [65] in combination with lattice calculations. Figure 10 shows the direct combination of experimental data and theoretical model calculations to extract physical quantities. The consistency of the results is of fundamental importance to validate that the experimentally created system is close to thermal equilibrium at the freeze-out and can be described by lattice QCD simulations, at least in the light quark sector [38].

#### 4. Summary

In summary, we have presented a review on the experimental measurements on charge fluctuations over various collision energies to search for phase transition and location of the critical point in the QCD phase diagram. We reviewed the results on total charge fluctuations as a function of collision centrality for GeV in Pb + Pb collisions at SPS and for and 200 GeV in Au + Au collisions at RHIC. It is observed that the relative fluctuations increase from central to peripheral collisions. The fluctuation observables from various models are compared with the experiment data. The results from models show centrality independent behavior. Similar centrality dependence has been observed for total charge fluctuation at the SPS and RHIC experiments. The absence of large dynamical fluctuation in excess of the participant superposition model indicates that there is no evidence of critical behavior related to the compressibility observable.

We have also presented the results from net-charge fluctuations from SPS, RHIC, and LHC energies. The net-charge fluctuations in terms of dynamical fluctuation measure as a function of collision centrality and are studied. increases monotonically from peripheral to central collisions and remains negative indicating the dominance of correlation of positive and negative charged particles. The product of and shows a monotonic decrease with increasing and approaches to the expectation from a QGP like scenario. For the highest RHIC energy, the measured value of fluctuation is close to the HG prediction, whereas at lower energies the results are higher than HG expectations. At LHC energy TeV, the value of fluctuation observable is significantly lower compared to lower energy results indicating that at LHC energy the fluctuations have their origin in the QGP phase.

The fluctuations of net-charge (net-proton, net-electric charge, and net-kaon) have been measured for various . Deviations are observed for both and of net-proton distributions from the Skellam and hadron resonance gas model for GeV. Maximum deviations from Skellam expectation are observed for and 27 GeV. The net-electric charge results from both PHENIX and STAR experiments are presented. values from both the experiments decrease with increasing . values from PHENIX remain constant and positive, between at all the collision energies within the statistical and systematic uncertainties. However, there is 25% increase of values at lower energies below GeV compared to higher energies. values from STAR experiment are constant at all energies within uncertainties, except for GeV which shows a negative value. The net-electric charge and net-kaon results do not observe any significant nonmonotonic behavior as a function of collision energy.

To quantify the excess of net-charge fluctuations due to the critical point, different baseline studies have been done by different authors after considering various physical phenomena. On phenomenological side, effects of resonance decay [78, 79], proton-antiproton correlations [71, 80], kinematic acceptance [72, 75, 76], nonextensive systems [81], regeneration [39, 82], excluded volume models [83], and so forth are studied. Besides, various studies have been carried out using different Monte Carlo based models, for example, UrQMD [84] and HIJING [85], to compare the experimental results. But it is to be noted that all the above-mentioned physical processes contribute differently at various center of mass energies; therefore, to find a unique ideal baseline estimate to compare the experimental results is still an unsettled issue. The freeze-out parameters are extracted using higher moments of net-charge distributions in combination with lattice QCD and HRG models. The extracted freeze-out parameters are in agreement with the thermal model parameterization.

First results from the RHIC BES-I program are intriguing and the statistics collected during BES-I at RHIC are not sufficient to locate the CEP in the QCD phase diagram. Hence, the RHIC BES-II program is proposed to have the precise measurements of these observables to map the QCD phase diagram. The upcoming RHIC BES-II program proposed to be in 2019-20 will cover the range from 5 to 20 GeV, which is the region of interest in the search for a critical point and first-order phase transition, identified by the results from BES-I and by model calculations. In addition, the CBM experiment at FAIR will perform a high-precision study of higher-order fluctuations at various beam energies in order to search for the elusive QCD critical point in the high net-baryon density region: MeV corresponding to GeV at SIS100. There are other programs at NICA to study the above physics in heavy-ion collisions in the range of 4–11 GeV. Hence, in the upcoming experimental programs at RHIC BES-II, FAIR, and NICA with upgraded detector system will help to measure the fluctuation variable with reduced uncertainties on the measurements and may find the exact location of the critical point and investigate whole region of phase diagram in more detail.

#### Disclosure

The present address of A. K. Mohanty is Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata 700064, India.

#### Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this manuscript.