#### Abstract

We present a review of the measurements of elliptic flow () of light nuclei (, , , , and ) from the RHIC and LHC experiments. Light (anti)nuclei have been compared with that of (anti)proton. We observed a similar trend in light nuclei to that in identified hadron with respect to the general observations such as dependence, low mass ordering, and centrality dependence. We also compared the difference of nuclei and antinuclei with the corresponding difference of of proton and antiproton at various collision energies. Qualitatively they depict similar behavior. We also compare the data on light nuclei to various theoretical models such as blast-wave and coalescence. We then present a prediction of for and using coalescence and blast-wave models.

#### 1. Introduction

The main goals of high energy heavy-ion collision experiments have primarily been to study the properties of Quark Gluon Plasma (QGP) and the other phase structures in the QCD phase diagram [1–11]. The energy densities created in these high energy collisions are similar to that found in the universe, microseconds after the Big Bang [5–8, 12–14]. Subsequently, the universe cooled down to form nuclei. It is expected that high energy heavy-ion collisions will allow studying the production of light nuclei such as , , , and their corresponding antinuclei. There are two possible production mechanisms for light (anti)nuclei. The first mechanism is thermal production of nucleus-antinucleus pairs in elementary nucleon-nucleon or parton-parton interactions [15–21]. However, due to their small (~few MeV) binding energies, the directly produced nuclei or antinuclei are likely to break up in the medium before escaping. The second mechanism is via final state coalescence of produced (anti)nucleons or from transported nucleons [22–36]. The quark coalescence as a mechanism of hadron production at intermediate transverse momentum has been well established by studying the number of constituent quarks (NCQ) scaling for of identified hadrons measured at RHIC [37–45]. Light nuclei may also be produced via coalescence of quarks similar to the hadrons. But the nuclei formed via quark coalescence are highly unlikely to survive in the high temperature environment due to their small binding energies. In case of hadron formation by quark coalescence, the momentum space distribution of quarks is not directly measurable in experiments. However, in case of nucleon coalescence, momentum space distributions of both the constituents (nucleons) and the products (nuclei) are measurable in heavy-ion collision experiments. Therefore, measurements of of light nuclei provide a tool to understand the production mechanism of light nuclei and freeze-out properties at a later stage of the evolution. It also provides an excellent opportunity to understand the mechanism of coalescence at work in high energy heavy-ion collisions.

The production of light (anti)nuclei has been studied extensively at lower energies in Bevelac at LBNL [24, 46–49], AGS at RHIC [50–53], and SPS at CERN [54–58]. In the AGS experiments, it was found that the coalescence parameter () is of similar magnitude for both and indicating similar freeze-out hypersurface of nucleons and antinucleons. Furthermore, the dependence of on collision energy and indicated that light nuclei production is strongly influenced by the source volume and transverse expansion profile of the system [58, 59]. In this paper, we review the results of elliptic flow of light nuclei measured at RHIC and LHC and discuss the possible mechanisms for the light nuclei production.

The paper is organized as follows. Section 2 briefly describes the definition of , identification of light (anti)nuclei in the experiments and measurement of of light (anti)nuclei. In Section 3, we present the results for minimum bias collisions from various experiments. We also discuss the centrality dependence, difference between nuclei and antinuclei , and the energy dependence of deuteron . In the same section, we present the atomic mass number scaling and also compare the experimental results with various theoretical models. Finally in Section 4, we summarize our observations and discuss the main conclusions of this review.

#### 2. Experimental Method

##### 2.1. Elliptic Flow

The azimuthal distribution of produced particles in heavy-ion collision can be expressed in terms of a Fourier series,where is the azimuthal angle of produced particle, is called the reaction plane angle, and the Fourier coefficients , , and so on are called flow coefficients [60]. is defined as the angle between the impact parameter vector and the -axis of the reference detector in the laboratory frame. Since it is impossible to measure the direction of impact parameter in heavy-ion collisions, a proxy of , namely, the event plane angle , is used for the flow analysis in heavy-ion collisions [61]. is measured with respect to the 2nd-order event plane angle [61]. is calculated using the azimuthal distribution of the produced particles. In an event with particles, the event plane angle is defined as [61] and are defined aswhere are weights given to each particle to optimise the event plane resolution [61, 62]. Usually the magnitude of particle transverse momentum is used as weights as increases with . Special techniques are followed while calculating the event plane angle so that it does not contain the particle of interest whose is to be calculated (self-correlation) and also the nonflow effects (e.g., jets and short range correlations) are removed as much as possible [41, 42, 61, 63]. Heavy-ion experiments use the -subevent plane method to calculate the elliptic flow of identified hadrons as well as for light nuclei. In this method, each event is divided into two subevents in two different -windows (e.g., positive and negative ). Then two subevent plane angles are calculated with the particles in each subevent. Each particle with a particular is then correlated with the subevent plane of the opposite . This ensures that the particle of interest is not included in the calculation of event plane angle. A finite gap is applied between the two subevents to reduce short range correlations which does not originate from flow.

The distribution of the event plane angles should be isotropic in the laboratory frame for a azimuthally isotropic detector. If the distribution of the event plane angles is not flat in the laboratory frame (due to detector acceptance and/or detector inefficiency) then special techniques are applied to make the distribution uniform. The popular methods to make the distribution uniform is the -weight and recentering [64, 65]. In the -weight method, one takes the actual azimuthal distribution of the produced particle, averaged over large sample of events, and then uses inverse of this distribution as weights while calculating the correlation of the particles with the event plane angle [64, 65]. In the recentering method, one subtracts and from the event-by-event and , respectively, where and denote the average of and over a large sample of similar events. The main disadvantage of applying one of these methods is that it does not remove the contribution from higher flow harmonics. Therefore, another correction method known as the shift correction is used to remove the effects coming from higher flow harmonics. In this method, one fits the distribution (after applying -weight and/or recentering method) averaged over all events, with a Fourier function. The Fourier coefficients from this fit (obtained as fit parameters) are used to shift of each event, to make the distribution uniform in the laboratory frame [64, 65].

Since the number of particles produced in heavy-ion collisions are finite, the calculated event plane angle may not coincide with . For this reason, the measured with respect to is corrected with the event plane resolution factor , where

In order to calculate the event plane resolution, one calculates two subevent plane angles and , where and correspond to two independent subevents. If the multiplicities of each subevent are approximately half of the full event plane, then the resolution of each of subevent plane can be calculated as [60, 61],However, the full event plane resolution can be expressed aswhere = and , are modified Bessel functions [60, 61]. The parameter is inversely proportional to the square-root of , the number of particles used to determine the event plane [60, 61]. To calculate the resolution for full event plane (), one has to solve (6) iteratively for the value of using the subevent plane resolution () which is calculated experimentally using (5). The value is then multiplied with as is proportional to and reused in (6) to calculate the resolution of the full event plane. In a case of very low magnitudes, the full event plane resolution can be approximately given as [60, 61]The procedure for calculating full and subevent plane resolutions using subevent plane angles and various approximations is discussed in detail in [60, 61].

##### 2.2. Data on Light Nuclei

For this review, we have collected light nuclei data from the STAR [63] and PHENIX [66] experiments at RHIC and ALICE experiment at LHC [67]. Table 1 summarizes the measurement of light nuclei available till date.

##### 2.3. Extraction of Light Nuclei

In heavy-ion collisions, light nuclei are primarily identified by comparing the mean ionization energy loss per unit length () in the Time Projection Chamber (TPC) with that from the theoretical predictions () [41, 42, 63, 67–71]. Light nuclei are also identified via the time-of-flight measurement techniques using the Time-of-Flight (TOF) detector [66, 67, 71–74].

In the STAR experiment, to identify light nuclei using TPC, a variable [63] is defined asThen the light nuclei yields are extracted from these -distributions in differential and bins either for minimum bias collisions or in selected centrality classes. The distribution is then fitted with a 2nd-order Fourier function; namely,The Fourier coefficient is called elliptic flow and is extracted from the fit. As we discussed in the previous subsection this measured is then corrected with the event plane resolution factor () [41, 42, 63].

In the ALICE experiment, light nuclei in the low region (<1.0 GeV/c for , ) are identified by comparing the variance () of the measured in the TPC with the corresponding theoretical estimate () [67, 71]. Light nuclei are considered identified if the measured lies within of the . On the other hand, the light nuclei yield is extracted from the mass squared () distribution using the TOF detector. The mass of each particle is calculated using the time-of-flight () from the TOF detector and the momentum () information from the TPC [66, 67, 71]. Both the ALICE and PHENIX experiments use the TOF detector to identify light nuclei at high (>1.0 GeV/c). The mass of a particle can be calculated using the TOF detector aswhere the track length and momentum are determined with the tracking detectors placed inside magnetic field [44, 66, 67, 71]. After getting for each particle, a selection cut is implemented to reject tracks which have their several away from the true value of the light nuclei, as done in the STAR experiment [63]. The ALICE experiment, on the other hand, defines a quantity such that , where is the mass of the light nuclei under study. The distribution of is then fitted with a Gaussian + exponential function for signal and an exponential function for the background [67]. Then of light nuclei is calculated by fitting () with the weighted function,where the total measured is the weighted sum of that from the signal () and background (). of the candidate particles are calculated using the scalar product method and corrected for the event plane resolution [67].

The PHENIX experiment calculates charged average of (anti) deuterons asThe quantity can readily be identified as the resolution of the event plane angle [66]. The resolution of full event plane is calculated with subevent planes (, ) estimated using two Beam-Beam Counter (BBC) detectors [44, 66]. The detailed procedure of calculating the full event plane resolution from subevents is already mentioned in the previous subsection. The large gap between the central TOF and the BBCs () reduces the effects of nonflow significantly [44, 66]. The nuclei calculated in PHENIX are also corrected for the contribution coming from backgrounds, mainly consisting of misidentification of other particles (e.g., protons) as nuclei. A dependent correction factor was applied on the total (referred to as ) such thatwhere is the measured for at a given , is the corrected of , and is the ratio of signal and signal + background.

#### 3. Results and Discussion

##### 3.1. General Aspects of Light Nuclei

Figure 1 shows the energy dependence of light (anti)nuclei for , 11.5, 19.6, 27, 39, 62.4, 200, and 2760 GeV. The panels are arranged by increasing energy from left to right and top to bottom. The dependence of of , , , , and is shown for 0–80% centrality in STAR, 20–60% centrality in PHENIX, and 30–40% centrality in ALICE. Since PHENIX and ALICE do not have measurements in the minimum bias collisions, we only show the results for mid-central collisions. The data points of PHENIX and ALICE correspond to inclusive . The general trend of nuclei of all species is the same: it increases with increasing . The slight difference of between STAR and PHENIX is due to the difference in centrality ranges. The centrality range for PHENIX is 20–60% and that for STAR is 0–80%.

From the trend in Figure 1 it seems that light nuclei show mass ordering; that is, heavier particles have smaller value compared to lighter ones, similar to of identified particles [41, 42, 44, 75]. In order to see the mass ordering effect more clearly, we restrict the –axis range to 2.5 GeV/c and compare of with of identified particles such as , ( in Pb + Pb), and as shown in Figure 2. We see that at all collision energies is lower than of the identified hadrons at a fixed value of . Although mass ordering is a theoretical expectation from the hydrodynamical approach to heavy-ion collisions [76], coalescence formalism for light nuclei can also give rise to this effect. Recent studies using AMPT and VISHNU hybrid model suggest that mass ordering is also expected from transport approach to heavy-ion collisions [77–79]. of light nuclei is negative for some collision energies as shown in Figure 1. This negative is expected to be the outcome of strong radial flow in heavy-ion collisions [80].

In order to study the energy dependence of light nuclei , we compare the deuteron from collision energy GeV to 2760 GeV as shown in Figure 3. The deuteron () shows energy dependence prominently for high ( GeV/c), where is highest for top collision energy ( GeV) and gradually decreases with decreasing collision energy. This energy dependent trend of light nuclei is similar to the energy dependence of identified hadron where () also decreases with decreasing collision energy [41, 42].

The STAR experiment has measured the difference of nuclei () and antinuclei () for collision energies , 27, 39, 62.4, and 200 GeV [63]. Figure 4 shows the difference of and as a function of collision energy. For comparison, the difference of proton and antiproton is also shown [41, 42]. We observe that the difference of and remains positive for GeV. However, for GeV the difference of and is almost zero. The difference of and qualitatively follows the same trend as seen for difference of and [41, 42]. It is easy to infer from simple coalescence model that light (anti)nuclei formed via coalescence of (anti)nucleons will also reflect similar difference in as the constituent nucleon and antinucleon. The difference in between particles and their antiparticles has been attributed to the chiral magnetic effect in finite baryon-density matter [81], different of produced and transported particles [82], different rapidity distributions for quarks and antiquarks [83], the conservation of baryon number, strangeness, and isospin [84], and different mean-field potentials acting on particles and their antiparticles [85–88].

The centrality dependence of light nuclei measured by the STAR and ALICE is shown in Figure 5. STAR has measured and in two different centrality ranges, namely, 0–30% and 30–80% for collision energies below GeV. In case of GeV, the light nuclei are measured in three different centrality ranges, namely, 0–10% (central), 10–40% (mid-central), and 40–80% (peripheral) as high statistics data were available. ALICE has measured inclusive in 6 different centrality ranges, namely, 0–5%, 5–10%, 10–20%, 20–30%, 30–40%, and 40–50%. We only present the results from 0–5%, 20–30%, and 40–50% centrality from ALICE as shown in Figure 5. of shows strong centrality dependence for all collision energies studied in the STAR experiment. We observe that more central events have lower compared to peripheral events. shows the same trend as for collision energies down to GeV.

The STAR experiment could not study centrality dependence of below GeV due to limited event statistics [63]. Comparing the centrality dependence of from STAR and ALICE we can see that both experiments show strong centrality dependence of light nuclei . The centrality dependence of light nuclei is analogous to the centrality dependence observed for identified nucleon () [89, 90].

##### 3.2. Mass Number Scaling and Model Comparison

It is expected from the formulations of coalescence model that if light nuclei are formed via the coalescence of nucleons then the elliptic flow of light nuclei, when divided by atomic mass number (), should scale with the elliptic flow of nucleons [91, 92]. Therefore, we expect that the light (anti)nuclei divided by should scale with () . Here, we essentially assume that of (anti)proton and (anti)neutron are the same as expected from the observed NCQ scaling of identified particle [41, 42]. Figure 6 shows the atomic mass number scaling of light nuclei from STAR, PHENIX, and ALICE experiments. Since ALICE does not have results in minimum bias events, we used both and from 30–40% centrality range. We observe that light nuclei from STAR and PHENIX show atomic mass number scaling up to GeV/c. However, deviation of the scaling of the order of 20% is observed for from ALICE. The scaling of light (anti)nuclei with (anti)proton suggests that light (anti)nuclei might have formed via coalescence of (anti)nucleons at a later stage of the evolution at RHIC energies for up to 1.5 GeV/c [22–32]. However, this simple picture of coalescence may not be holding for ALICE experiment at LHC energies. On the contrary, there is another method to produce light nuclei, for example, by thermal production in which it is assumed that light nuclei are produced thermally like any other primary particles [17–21]. Various thermal model studies have successfully reproduced the different ratios of produced particles as well as light nuclei in heavy-ion collisions [17–21].

In order to investigate the success of these models, both STAR and ALICE have compared the elliptic flow of light nuclei with the predictions from blast-wave models [63, 67]. Figure 7 shows the of light nuclei predicted from blast-wave model using the parameters obtained from fits to the identified particles [67, 93]. We observe that blast-wave model reproduces of light nuclei from STAR with moderate success except for low (<1.0 GeV/c), where of () are underpredicted for all collision energies. However, the blast-wave model seems to successfully reproduce the from ALICE. The low relative production of light nuclei compared to identified nucleons at RHIC collisions energies supports the procedure of light nuclei production via coalescence mechanism [22–32]. However, the success of blast-wave model in reproducing the nuclei at LHC and moderate success at RHIC suggest that the light nuclei production is also supported by thermal process [17–21]. The light nuclei production in general might be a more complicated coalescence process, for example, coalescence of nucleons in the local rest frame of the fluid cell. This scenario might give rise to deviations from simple scaling [63].

At RHIC energies the light nuclei have been compared with results from a hybrid AMPT + coalescence model [63]. A Multiphase Transport (AMPT) model is an event generator with Glauber Monte Carlo initial state [94]. The AMPT model includes Zhang’s Partonic Cascade (ZPC) model for initial partonic interactions and A Relativistic Transport (ART) model for later hadronic interactions [94]. The nucleon phase-space information from the AMPT model is fed to the coalescence model to generate light nuclei [63, 95]. Figure 8 shows the light nuclei from the coalescence model and compared to the data. The coalescence model prediction for in Pb + Pb collisions at GeV is taken from [96]. The coalescence model fairly reproduces the measurement from data for all collision energies except for the lowest energy GeV. The AMPT model generates nucleon from both partonic and hadronic interactions for all the collision energies presented. However, increased hadronic interactions compared to partonic, at lowest collision energies, are not implemented in the AMPT + coalescence model. This could be the reason behind the deviation of the data from the model predictions at lowest collision energy [41, 42].

We have performed simultaneous fit to the and spectra of identified hadrons + light nuclei using the same blast-wave model as used in [67, 75]. The simultaneous fit of and spectra for measurements from the PHENIX and the ALICE experiment is shown in Figure 9. We find that the inclusion of light nuclei results to the fit does not change the fit results compared to the blast-wave fit performed only on identified hadron and spectra. This indicates that the light nuclei and spectra are well described by the blast-wave model.

**(a)**

**(b)**

**(c)**

**(d)**

##### 3.3. Model Prediction of and

We have predicted of and using the simple coalescence and blast-wave model. Since protons and neutrons have similar masses and the same number of constituent quarks, they should exhibit similar collective behavior and, hence, similar magnitude of . Therefore, we parametrize the elliptic flow of using the fit formula [97],where , , , and are fit parameters and is the constituent quark number of the particle [97]. The fit to (solid lines) from the PHENIX and ALICE experiment is shown in Figures 10(a) and 10(b), respectively. Assuming similar magnitude of neutron as that of proton, we then predict of and as

**(a)**

**(b)**

This simplified coalescence model prediction of and is shown in Figures 10(a) and 10(b) as blue (dotted and dashed) lines. For comparison, the blast-wave model predicted of and from the fit parameters obtained in Figure 9 are also shown in red (dotted and dashed) lines. We observe characteristic difference in the prediction of and from the coalescence and the blast-wave model. As one expects from the mass ordering effect of blast-wave model, of and are predicted to be almost zero in the intermediate range (1.0 GeV/c). On the other hand, the simple coalescence model predicts orders of magnitude higher compared to blast-wave for both and in the same range. Hence, experimental measurements of and in future would significantly improve our knowledge on the mechanisms of light nuclei formation in heavy-ion collisions [71, 98–100].

#### 4. Summary and Conclusions

We have presented a review of elliptic flow of light nuclei (, , and ) and antinuclei ( and ) from STAR experiment and inclusive from PHENIX at RHIC and ALICE at LHC. Similar to identified hadrons, the light nuclei show a monotonic rise with increasing and mass ordering at low for all measured collision energies. The beam energy dependence of is small at intermediate and only prominent at high , which is similar to the trend as observed for the charged hadron . The of nuclei and antinuclei are of similar magnitude for top collision energies at RHIC but at lower collision energies; the difference in between nuclei and antinuclei qualitatively follows the difference in proton and antiproton . The centrality dependence of light (anti)nuclei is similar to that of identified hadrons .

Light (anti)nuclei is found to follow the atomic mass number () scaling for almost all collision energies at RHIC suggesting coalescence as the underlying process for the light nuclei production in heavy-ion collisions. However, a deviation from mass number scaling at the level of 20% is observed at LHC. This indicates that a simple coalescence process may not be the only underlying mechanism for light nuclei production. Furthermore, a transport-plus-coalescence model study is found to approximately reproduce the light nuclei measured at RHIC and LHC. The agreement of coalescence model with the data from PHENIX and STAR is imperceptibly better than the blast-wave model. However, at the LHC energy, the light nuclei are better described by blast-wave model rather than the simple coalescence model. The coalescence mechanism, intuitively, should be the prominent process of light nuclei production. However, the breaking of mass scaling at LHC energy and success of blast-wave model prevent us from drawing any definitive conclusion on the light nuclei production mechanism.

We observed orders of magnitude difference in and as predicted by blast-wave and coalescence model. The blast-wave model predicts almost zero for and up to GeV/c, whereas the coalescence model predicts significant for and at same range. Hence, the precise measurements of and in the future can significantly improve the knowledge of the light nuclei production mechanism in heavy-ion collisions.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The authors thank STAR collaboration, PHENIX collaboration, and ALICE collaboration for providing the light nuclei data and the model predictions. This work is supported by DAE-BRNS Project (Grant no. 2010/21/15-BRNS/2026) and Dr. C. Jena is supported by 12th plan project (PIC no. 12-R&D-NIS-5.11-0300).