We report the beam energy and collision centrality dependence of fifth and sixth order cumulants ($C_{5}$, $C_{6}$) and factorial cumulants ($\kappa_{5}$, $\kappa_{6}$) of net-proton and proton distributions, from $\sqrt{s_{NN}} = 3 - 200$ GeV Au+Au collisions at RHIC. The net-proton cumulant ratios generally follow the hierarchy expected from QCD thermodynamics, except for the case of collisions at $\sqrt{s_{NN}}$ = 3 GeV. $C_{6}/C_{2}$ for 0-40% centrality collisions is increasingly negative with decreasing $\sqrt{s_{NN}}$, while it is positive for the lowest $\sqrt{s_{NN}}$ studied. These observed negative signs are consistent with QCD calculations (at baryon chemical potential, $\mu_{B} \leq$ 110 MeV) that include a crossover quark-hadron transition. In addition, for $\sqrt{s_{NN}} \geq$ 11.5 GeV, the measured proton $\kappa_{n}$, within uncertainties, does not support the two-component shape of proton distributions that would be expected from a first-order phase transition. Taken in combination, the hyper-order proton number fluctuations suggest that the structure of QCD matter at high baryon density, $\mu_{B}\sim 750$ MeV ($\sqrt{s_{NN}}$ = 3 GeV) is starkly different from those at vanishing $\mu_{B}\sim 20$MeV ($\sqrt{s_{NN}}$ = 200 GeV and higher).

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in \pp and \AuAu collisions at \sqrtsn $= 200$ GeV. The substructure observables include SoftDrop groomed momentum fraction (\zg), groomed jet radius (\rg), and subjet momentum fraction (\zsj) and opening angle (\tsj). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for \pp collisions and are compared to leading order Monte Carlo models. The subjet \tsj~distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in \AuAu collisions. We compare data from \AuAu collisions to those from \pp which are embedded in minimum-bias \AuAu events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than \zg~and \rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuum-like. We also report measurements of the differential di-jet momentum imbalance ($A_{\rm{J}}$) for jets of varying \tsj. We find no qualitative differences in energy loss signatures for varying angular scales in the range $0.1 < $\tsj $ < 0.3$, leading to the possible interpretation that energy loss in this population of high momentum di-jet pairs, is due to soft medium-induced gluon radiation from a single color-charge as it traverses the medium.

The deconfined quark-gluon plasma (QGP) created in relativistic heavy-ion collisions enables the exploration of the fundamental properties of matter under extreme conditions. Non-central collisions can produce strong magnetic fields on the order of $10^{18}$ Gauss, which offers a probe into the electrical conductivity of the QGP. In particular, quarks and anti-quarks carry opposite charges and receive contrary electromagnetic forces that alter their momenta. This phenomenon can be manifested in the collective motion of final-state particles, specifically in the rapidity-odd directed flow, denoted as $v_1(\mathsf{y})$. Here we present the charge-dependent measurements of $dv_1/d\mathsf{y}$ near midrapidities for $\pi^{\pm}$, $K^{\pm}$, and $p(\bar{p})$ in Au+Au and isobar ($_{44}^{96}$Ru+$_{44}^{96}$Ru and $_{40}^{96}$Zr+$_{40}^{96}$Zr) collisions at $\sqrt{s_{\rm NN}}=$ 200 GeV, and in Au+Au collisions at 27 GeV, recorded by the STAR detector at the Relativistic Heavy Ion Collider. The combined dependence of the $v_1$ signal on collision system, particle species, and collision centrality can be qualitatively and semi-quantitatively understood as several effects on constituent quarks. While the results in central events can be explained by the $u$ and $d$ quarks transported from initial-state nuclei, those in peripheral events reveal the impacts of the electromagnetic field on the QGP. Our data put valuable constraints on the electrical conductivity of the QGP in theoretical calculations.

The longitudinal and transverse spin transfers to $\Lambda$ ($\overline{\Lambda}$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $\Lambda$ and $\overline{\Lambda}$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|\eta_{\Lambda(\overline{\Lambda})}|$$<$ 1.2 and transverse momentum $p_{T,{\Lambda(\overline{\Lambda})}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions.

The polarization of $\Lambda$ and $\bar{\Lambda}$ hyperons along the beam direction has been measured relative to the second and third harmonic event planes in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. This is the first experimental evidence of the hyperon polarization by the triangular flow originating from the initial density fluctuations. The amplitudes of the sine modulation for the second and third harmonic results are comparable in magnitude, increase from central to peripheral collisions, and show a mild $p_T$ dependence. The azimuthal angle dependence of the polarization follows the vorticity pattern expected due to elliptic and triangular anisotropic flow, and qualitatively disagree with most hydrodynamic model calculations based on thermal vorticity and shear induced contributions. The model results based on one of existing implementations of the shear contribution lead to a correct azimuthal angle dependence, but predict centrality and $p_T$ dependence that still disagree with experimental measurements. Thus, our results provide stringent constraints on the thermal vorticity and shear-induced contributions to hyperon polarization. Comparison to previous measurements at RHIC and the LHC for the second-order harmonic results shows little dependence on the collision system size and collision energy.

Global polarizations ($P$) of $\Lambda$ ($\bar{\Lambda}$) hyperons have been observed in non-central heavy-ion collisions. The strong magnetic field primarily created by the spectator protons in such collisions would split the $\Lambda$ and $\bar{\Lambda}$ global polarizations ($\Delta P = P_{\Lambda} - P_{\bar{\Lambda}} < 0$). Additionally, quantum chromodynamics (QCD) predicts topological charge fluctuations in vacuum, resulting in a chirality imbalance or parity violation in a local domain. This would give rise to an imbalance ($\Delta n = \frac{N_{\text{L}} - N_{\text{R}}}{\langle N_{\text{L}} + N_{\text{R}} \rangle} \neq 0$) between left- and right-handed $\Lambda$ ($\bar{\Lambda}$) as well as a charge separation along the magnetic field, referred to as the chiral magnetic effect (CME). This charge separation can be characterized by the parity-even azimuthal correlator ($\Delta\gamma$) and parity-odd azimuthal harmonic observable ($\Delta a_{1}$). Measurements of $\Delta P$, $\Delta\gamma$, and $\Delta a_{1}$ have not led to definitive conclusions concerning the CME or the magnetic field, and $\Delta n$ has not been measured previously. Correlations among these observables may reveal new insights. This paper reports measurements of correlation between $\Delta n$ and $\Delta a_{1}$, which is sensitive to chirality fluctuations, and correlation between $\Delta P$ and $\Delta\gamma$ sensitive to magnetic field in Au+Au collisions at 27 GeV. For both measurements, no correlations have been observed beyond statistical fluctuations.

In relativistic heavy-ion collisions, a global spin polarization, $P_\mathrm{H}$, of $\Lambda$ and $\bar{\Lambda}$ hyperons along the direction of the system angular momentum was discovered and measured across a broad range of collision energies and demonstrated a trend of increasing $P_\mathrm{H}$ with decreasing $\sqrt{s_{NN}}$. A splitting between $\Lambda$ and $\bar{\Lambda}$ polarization may be possible due to their different magnetic moments in a late-stage magnetic field sustained by the quark-gluon plasma which is formed in the collision. The results presented in this study find no significant splitting at the collision energies of $\sqrt{s_{NN}}=19.6$ and $27$ GeV in the RHIC Beam Energy Scan Phase II using the STAR detector, with an upper limit of $P_{\bar{\Lambda}}-P_{\Lambda}<0.24$% and $P_{\bar{\Lambda}}-P_{\Lambda}<0.35$%, respectively, at a 95% confidence level. We derive an upper limit on the na\"ive extraction of the late-stage magnetic field of $B<9.4\cdot10^{12}$ T and $B<1.4\cdot10^{13}$ T at $\sqrt{s_{NN}}=19.6$ and $27$ GeV, respectively, although more thorough derivations are needed. Differential measurements of $P_\mathrm{H}$ were performed with respect to collision centrality, transverse momentum, and rapidity. With our current acceptance of $|y|<1$ and uncertainties, we observe no dependence on transverse momentum and rapidity in this analysis. These results challenge multiple existing model calculations following a variety of different assumptions which have each predicted a strong dependence on rapidity in this collision-energy range.

We report a new measurement of the production of electrons from open heavy-flavor hadron decays (HFEs) at mid-rapidity ($|y|<$ 0.7) in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. Invariant yields of HFEs are measured for the transverse momentum range of $3.5 < p_{\rm T} < 9$ GeV/$c$ in various configurations of the collision geometry. The HFE yields in head-on Au+Au collisions are suppressed by approximately a factor of 2 compared to that in $p$+$p$ collisions scaled by the average number of binary collisions, indicating strong interactions between heavy quarks and the hot and dense medium created in heavy-ion collisions. Comparison of these results with models provides additional tests of theoretical calculations of heavy quark energy loss in the quark-gluon plasma.

In this erratum we report changes on the $D^0$ $p_T$ spectra and nuclear modification factor ($R_{AA}$) in Au+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 200 GeV by fixing the errors in the efficiency and selection criteria that affected the Au+Au results. The p+p reference spectrum has changed as well and is updated with new fragmentation parameters.

The elliptic ($v_2$) and triangular ($v_3$) azimuthal anisotropy coefficients in central $^{3}$He+Au, $d$+Au, and $p$+Au collisions at $\mbox{$\sqrt{s_{\mathrm{NN}}}$}$ = 200 GeV are measured as a function of transverse momentum ($p_{\mathrm{T}}$) at mid-rapidity ($|\eta|<$0.9), via the azimuthal angular correlation between two particles both at $|\eta|<$0.9. While the $v_2(p_{\mathrm{T}})$ values depend on the colliding systems, the $v_3(p_{\mathrm{T}})$ values are system-independent within the uncertainties, suggesting an influence on eccentricity from sub-nucleonic fluctuations in these small-sized systems. These results also provide stringent constraints for the hydrodynamic modeling of these systems.