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.

We present the first measurements of transverse momentum spectra of $\pi^{\pm}$, $K^{\pm}$, $p(\bar{p})$ at midrapidity ($|y| < 0.1$) in U+U collisions at $\sqrt{s_{NN}}$ = 193 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC). The centrality dependence of particle yields, average transverse momenta, particle ratios and kinetic freeze-out parameters are discussed. The results are compared with the published results from Au+Au collisions at $\sqrt{s_{NN}} =$ 200 GeV in STAR. The results are also compared to those from A Multi Phase Transport (AMPT) model.

We present measurements of $\pi^-$ and $\pi^+$ elliptic flow, $v_2$, at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}} =$ 200, 62.4, 39, 27, 19.6, 11.5 and 7.7 GeV, as a function of event-by-event charge asymmetry, $A_{ch}$, based on data from the STAR experiment at RHIC. We find that $\pi^-$ ($\pi^+$) elliptic flow linearly increases (decreases) with charge asymmetry for most centrality bins at $\sqrt{s_{_{\rm NN}}} = \text{27 GeV}$ and higher. At $\sqrt{s_{_{\rm NN}}} = \text{200 GeV}$, the slope of the difference of $v_2$ between $\pi^-$ and $\pi^+$ as a function of $A_{ch}$ exhibits a centrality dependence, which is qualitatively similar to calculations that incorporate a chiral magnetic wave effect. Similar centrality dependence is also observed at lower energies.

Dihadron azimuthal correlations containing a high transverse momentum ($p_T$) trigger particle are sensitive to the properties of the nuclear medium created at RHIC through the strong interactions occurring between the traversing parton and the medium, i.e. jet-quenching. Previous measurements revealed a strong modification to dihadron azimuthal correlations in Au+Au collisions with respect to p+p and d+Au collisions. The modification increases with the collision centrality, suggesting a path-length or energy density dependence to the jet-quenching effect. This paper reports STAR measurements of dihadron azimuthal correlations in mid-central (20-60%) Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=200$ GeV as a function of the trigger particle's azimuthal angle relative to the event plane, $\phi_s=|\phi_t-\psi_{\rm EP}|$. The azimuthal correlation is studied as a function of both the trigger and associated particle $p_T$. The subtractions of the combinatorial background and anisotropic flow, assuming Zero Yield At Minimum (ZYAM), are described. The correlation results are first discussed with subtraction of the even harmonic (elliptic and quadrangular) flow backgrounds. The away-side correlation is strongly modified, and the modification varies with $\phi_s$, with a double-peak structure for out-of-plane trigger particles. The near-side ridge (long range pseudo-rapidity $\Delta\eta$ correlation) appears to drop with increasing $\phi_s$ while the jet-like component remains approximately constant. The correlation functions are further studied with subtraction of odd harmonic triangular flow background arising from fluctuations. It is found that the triangular flow, while responsible for the majority of the amplitudes, is not sufficient to explain the $\phi_s$-dependence of the ridge or the away-side double-peak structure. ...

STAR measurements of dihadron azimuthal correlations ($\Delta\phi$) are reported in mid-central (20-60\%) Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=200$ GeV as a function of the trigger particle's azimuthal angle relative to the event plane, $\phi_{s}=|\phi_{t}-\psi_{\rm EP}|$. The elliptic ($v_2$), triangular ($v_3$), and quadratic ($v_4$) flow harmonic backgrounds are subtracted using the Zero Yield At Minimum (ZYAM) method. The results are compared to minimum-bias d+Au collisions. It is found that a finite near-side ($|\Delta\phi|<\pi/2$) long-range pseudorapidity correlation (ridge) is present in the in-plane direction ($\phi_{s}\sim 0$). The away-side ($|\Delta\phi|>\pi/2$) correlation shows a modification from d+Au data, varying with $\phi_{s}$. The modification may be a consequence of pathlength-dependent jet-quenching and may lead to a better understanding of high-density QCD.

We report on measurements of sequential $\Upsilon$ suppression in Au+Au collisions at $\sqrt{s_{_\mathrm{NN}}}$ = 200 GeV with the STAR detector at the Relativistic Heavy Ion Collider (RHIC) through both the dielectron and dimuon decay channels. In the 0-60% centrality class, the nuclear modification factors ($R_{\mathrm{AA}}$), which quantify the level of yield suppression in heavy-ion collisions compared to $p$+$p$ collisions, for $\Upsilon$(1S) and $\Upsilon$(2S) are $0.40 \pm 0.03~\textrm{(stat.)} \pm 0.03~\textrm{(sys.)} \pm 0.09~\textrm{(norm.)}$ and $0.26 \pm 0.08~\textrm{(stat.)} \pm 0.02~\textrm{(sys.)} \pm 0.06~\textrm{(norm.)}$, respectively, while the upper limit of the $\Upsilon$(3S) $R_{\mathrm{AA}}$ is 0.17 at a 95% confidence level. This provides experimental evidence that the $\Upsilon$(3S) is significantly more suppressed than the $\Upsilon$(1S) at RHIC. The level of suppression for $\Upsilon$(1S) is comparable to that observed at the much higher collision energy at the Large Hadron Collider. These results point to the creation of a medium at RHIC whose temperature is sufficiently high to strongly suppress excited $\Upsilon$ states.

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.

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.

We report a measurement of cumulants and correlation functions of event-by-event proton multiplicity distributions from fixed-target Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV measured by the STAR experiment. Protons are identified within the rapidity ($y$) and transverse momentum ($p_{\rm T}$) region $-0.9 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$ in the center-of-mass frame. A systematic analysis of the proton cumulants and correlation functions up to sixth-order as well as the corresponding ratios as a function of the collision centrality, $p_{\rm T}$, and $y$ are presented. The effect of pileup and initial volume fluctuations on these observables and the respective corrections are discussed in detail. The results are compared to calculations from the hadronic transport UrQMD model as well as a hydrodynamic model. In the most central 5% collisions, the value of proton cumulant ratio $C_4/C_2$ is negative, drastically different from the values observed in Au+Au collisions at higher energies. Compared to model calculations including Lattice QCD, a hadronic transport model, and a hydrodynamic model, the strong suppression in the ratio of $C_4/C_2$ at 3 GeV Au+Au collisions indicates an energy regime dominated by hadronic interactions.

Quark interactions with topological gluon configurations can induce local chirality imbalance and parity violation in quantum chromodynamics, which can lead to the chiral magnetic effect (CME) -- an electric charge separation along the strong magnetic field in relativistic heavy-ion collisions. The CME-sensitive azimuthal correlator observable ($\Delta\gamma$) is contaminated by background arising, in part, from resonance decays coupled with elliptic anisotropy ($v_{2}$). We report here differential measurements of the correlator as a function of the pair invariant mass ($m_{\rm inv}$) in 20-50% centrality Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$= 200 GeV by the STAR experiment at RHIC. Strong resonance background contributions to $\Delta\gamma$ are observed. At large $m_{\rm inv}$ where this background is significantly reduced, the $\Delta\gamma$ value is found to be significantly smaller. An event-shape-engineering technique is deployed to determine the $v_{2}$ background shape as a function of $m_{\rm inv}$. We extract a $v_2$-independent and $m_{\rm inv}$-averaged signal $\Delta\gamma_{\rm sig}$ = (0.03 $\pm$ 0.06 $\pm$ 0.08) $\times10^{-4}$, or $(2\pm4\pm5)\%$ of the inclusive $\Delta\gamma(m_{\rm inv}>0.4$ GeV/$c^2$)$ =(1.58 \pm 0.02 \pm 0.02) \times10^{-4}$, within pion $p_{T}$ = 0.2 - 0.8~\gevc and averaged over pseudorapidity ranges of $-1 < \eta < -0.05$ and $0.05 < \eta < 1$. This represents an upper limit of $0.23\times10^{-4}$, or $15\%$ of the inclusive result, at $95\%$ confidence level for the $m_{\rm inv}$-integrated CME contribution.