According to first-principle lattice QCD calculations, the transition from quark-gluon plasma to hadronic matter is a smooth crossover in the region μB ≤ T c. In this range the ratio, C6=C2, of net-baryon distributions are predicted to be negative. In this Letter, we report the first measurement of the midrapidity net-proton C6=C2 from 27, 54.4, and 200 GeV Au þ Au collisions at the Relativistic Heavy Ion Collider (RHIC). The dependence on collision centrality and kinematic acceptance in (p T , y) are analyzed. While for 27 and 54.4 GeV collisions the C6=C2 values are close to zero within uncertainties, it is observed that for 200 GeV collisions, the C6=C2 ratio becomes progressively negative from peripheral to central collisions. Transport model calculations without critical dynamics predict mostly positive values except for the most central collisions within uncertainties. These observations seem to favor a smooth crossover in the high-energy nuclear collisions at top RHIC energy.

We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/d\eta$) and follows a scaling behavior. The $dN_{ch}/d\eta$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$\sigma$ and 3.4$\sigma$, respectively, giving a combined significance of 4.1$\sigma$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

We report the first measurement of rapidity-odd directed flow ($v_{1}$) for $D^{0}$ and $\overline{D^{0}}$ mesons at mid-rapidity ($|y| < 0.8$) in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200\,GeV using the STAR detector at the Relativistic Heavy Ion Collider. In 10--80\% Au+Au collisions, the slope of the $v_{1}$ rapidity dependence ($dv_{1}/dy$), averaged over $D^{0}$ and $\overline{D^{0}}$ mesons, is -0.080 $\pm$ 0.017 (stat.) $\pm$ 0.016 (syst.) for transverse momentum $p_{\rm T}$ above 1.5~GeV/$c$. The absolute value of $D^0$-meson $dv_1/dy$ is about 25 times larger than that for charged kaons, with 3.4$\sigma$ significance. These data give a unique insight into the initial tilt of the produced matter, and offer constraints on the geometric and transport parameters of the hot QCD medium created in relativistic heavy-ion collisions.

A decisive experimental test of the Chiral Magnetic Effect (CME) is considered one of the major scientific goals at the Relativistic Heavy-Ion Collider (RHIC) towards understanding the nontrivial topological fluctuations of the Quantum Chromodynamics vacuum. In heavy-ion collisions, the CME is expected to result in a charge separation phenomenon across the reaction plane, whose strength could be strongly energy dependent. The previous CME searches have been focused on top RHIC energy collisions. In this Letter, we present a low energy search for the CME in Au+Au collisions at $\sqrt{s_{_{\rm{NN}}}}=27$ GeV. We measure elliptic flow scaled charge-dependent correlators relative to the event planes that are defined at both mid-rapidity $|\eta|<1.0$ and at forward rapidity $2.1 < |\eta|<5.1$. We compare the results based on the directed flow plane ($\Psi_1$) at forward rapidity and the elliptic flow plane ($\Psi_2$) at both central and forward rapidity. The CME scenario is expected to result in a larger correlation relative to $\Psi_1$ than to $\Psi_2$, while a flow driven background scenario would lead to a consistent result for both event planes. In 10-50% centrality, results using three different event planes are found to be consistent within experimental uncertainties, suggesting a flow driven background scenario dominating the measurement. We obtain an upper limit on the deviation from a flow driven background scenario at the 95% confidence level. This work opens up a possible road map towards future CME search with the high statistics data from the RHIC Beam Energy Scan Phase-II.

The elliptic, $v_2$, triangular, $v_3$, and quadrangular, $v_4$, azimuthal anisotropic flow coefficients are measured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event plane and four-particle cumulant methods are reported for the pseudo-rapidity range $|\eta|<0.8$ at different collision centralities and as a function of transverse momentum, $p_{\rm T}$, out to $p_{\rm T}=20$ GeV/$c$. The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for $p_{\rm T}>8$ GeV/$c$. The small $p_{\rm T}$ dependence of the difference between elliptic flow results obtained from the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations up to $p_{\rm T}=8$ GeV/$c$. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that of pions out to at least $p_{\rm T}=8$ GeV/$c$ indicating that the particle type dependence persists out to high $p_{\rm T}$.

The chiral magnetic effect (CME) is predicted to occur as a consequence of a local violation of $\cal P$ and $\cal CP$ symmetries of the strong interaction amidst a strong electro-magnetic field generated in relativistic heavy-ion collisions. Experimental manifestation of the CME involves a separation of positively and negatively charged hadrons along the direction of the magnetic field. Previous measurements of the CME-sensitive charge-separation observables remain inconclusive because of large background contributions. In order to better control the influence of signal and backgrounds, the STAR Collaboration performed a blind analysis of a large data sample of approximately 3.8 billion isobar collisions of $^{96}_{44}$Ru+$^{96}_{44}$Ru and $^{96}_{40}$Zr+$^{96}_{40}$Zr at $\sqrt{s_{\rm NN}}=200$ GeV. Prior to the blind analysis, the CME signatures are predefined as a significant excess of the CME-sensitive observables in Ru+Ru collisions over those in Zr+Zr collisions, owing to a larger magnetic field in the former. A precision down to 0.4% is achieved, as anticipated, in the relative magnitudes of the pertinent observables between the two isobar systems. Observed differences in the multiplicity and flow harmonics at the matching centrality indicate that the magnitude of the CME background is different between the two species. No CME signature that satisfies the predefined criteria has been observed in isobar collisions in this blind analysis.

The chiral magnetic effect (CME) refers to charge separation along a strong magnetic field due to imbalanced chirality of quarks in local parity and charge-parity violating domains in quantum chromodynamics. The experimental measurement of the charge separation is made difficult by the presence of a major background from elliptic azimuthal anisotropy. This background and the CME signal have different sensitivities to the spectator and participant planes, and could thus be determined by measurements with respect to these planes. We report such measurements in Au+Au collisions at a nucleon-nucleon center-of-mass energy of 200 GeV at the Relativistic Heavy-Ion Collider. It is found that the charge separation, with the flow background removed, is consistent with zero in peripheral (large impact parameter) collisions. Some indication of finite CME signals is seen in mid-central (intermediate impact parameter) collisions. Significant residual background effects may, however, still be present.

We report a systematic measurement of cumulants, $C_{n}$, for net-proton, proton and antiproton multiplicity distributions, and correlation functions, $\kappa_n$, for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{\mathrm {NN}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4 and 200 GeV. The $C_{n}$ and $\kappa_n$ are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, $y$, and transverse momentum, $p_{T}$. The data were taken during the first phase of the Beam Energy Scan (BES) program (2010 -- 2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity ($|y| <$ 0.5) and transverse momentum 0.4 $<$$p_{\rm T}$$<$ 2.0 GeV/$c$, using the STAR detector at RHIC. We observe a non-monotonic energy dependence ($\sqrt{s_{\mathrm {NN}}}$ = 7.7 -- 62.4 GeV) of the net-proton $C_{4}$/$C_{2}$ with the significance of 3.1$\sigma$ for the 0-5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{\mathrm {NN}}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, $\kappa_2$, of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, $\kappa_4$, of protons plays a role in determining the energy dependence of proton $C_4/C_1$ below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.

Understanding gluon density distributions and how they are modified in nuclei are among the most important goals in nuclear physics. In recent years, diffractive vector meson production measured in ultra-peripheral collisions (UPCs) at heavy-ion colliders has provided a new tool for probing the gluon density. In this Letter, we report the first measurement of $J/\psi$ photoproduction off the deuteron in UPCs at the center-of-mass energy $\sqrt{s_{_{\rm NN}}}=200~\rm GeV$ in d$+$Au collisions. The differential cross section as a function of momentum transfer $-t$ is measured. In addition, data with a neutron tagged in the deuteron-going Zero-Degree Calorimeter is investigated for the first time, which is found to be consistent with the expectation of incoherent diffractive scattering at low momentum transfer. Theoretical predictions based on the Color Glass Condensate saturation model and the gluon shadowing model are compared with the data quantitatively. A better agreement with the saturation model has been observed. With the current measurement, the results are found to be directly sensitive to the gluon density distribution of the deuteron and the deuteron breakup, which provides insights into the nuclear gluonic structure.

In high-energy heavy-ion collisions, the correlations between the emitted particles can be used as a probe to gain insight into the charge creation mechanisms. In this Letter, we report the first results of such studies using the electric charge balance function in the relative pseudorapidity ($\Delta\eta$) and azimuthal angle ($\Delta\varphi$) in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV with the ALICE detector at the Large Hadron Collider. The width of the balance function decreases with growing centrality (i.e. for more central collisions) in both projections. This centrality dependence is not reproduced by HIJING, while AMPT, a model which incorporates strings and parton rescattering, exhibits qualitative agreement with the measured correlations in $\Delta\varphi$ but fails to describe the correlations in $\Delta\eta$. A thermal blast-wave model incorporating local charge conservation and tuned to describe the $p_{\rm T}$ spectra and v$_2$ measurements reported by ALICE, is used to fit the centrality dependence of the width of the balance function and to extract the average separation of balancing charges at freeze-out. The comparison of our results with measurements at lower energies reveals an ordering with $\sqrt{s_{\rm NN}}$: the balance functions become narrower with increasing energy for all centralities. This is consistent with the effect of larger radial flow at the LHC energies but also with the late stage creation scenario of balancing charges. However, the relative decrease of the balance function widths in $\Delta\eta$ and $\Delta\varphi$ with centrality from the highest SPS to the LHC energy exhibits only small differences. This observation cannot be interpreted solely within the framework where the majority of the charge is produced at a later stage in the evolution of the heavy--ion collision.