We measure triangular flow relative to the reaction plane at 3 GeV center-of-mass energy in Au+Au collisions at RHIC. A significant $v_3$ signal is observed for protons, whose magnitude increases for higher rapidity, higher transverse momentum, and more peripheral collisions. The triangular flow is essentially rapidity-odd with a rapidity slope at mid-rapidity, $dv_3/dy|_{(y=0)}$, opposite in sign compared to the slope for directed flow. No significant $v_3$ signal is observed for charged pions and kaons. Comparisons with models suggest that a mean field potential is required to describe these results, and that the triangular shape of the participant nucleons is the result of stopping and nuclear geometry.

A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV, and Cu+Cu collisions sqrt(s_NN) = 22.5, 62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for non-dynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. An analysis of Negative Binomial Distribution fits to the multiplicity distributions demonstrates that the heavy ion data exhibit weak clustering properties.

The PHENIX experiment at RHIC has measured transverse energy and charged particle multiplicity at mid-rapidity in Au+Au collisions at sqrt(s_NN) = 19.6, 130 and 200 GeV as a function of centrality. The presented results are compared to measurements from other RHIC experiments, and experiments at lower energies. The sqrt(s_NN) dependence of dE_T/deta and dN_ch/deta per pair of participants is consistent with logarithmic scaling for the most central events. The centrality dependence of dE_T/deta and dN_ch/deta is similar at all measured incident energies. At RHIC energies the ratio of transverse energy per charged particle was found independent of centrality and growing slowly with sqrt(s_NN). A survey of comparisons between the data and available theoretical models is also presented.

We report here the first observation of directed flow ($v_1$) of the hypernuclei $^3_{\Lambda}$H and $^4_{\Lambda}$H in mid-central Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 3 GeV at RHIC. These data are taken as part of the beam energy scan program carried out by the STAR experiment. From 165 $\times$ 10$^{6}$ events in 5%-40% centrality, about 8400 $^3_{\Lambda}$H and 5200 $^4_{\Lambda}$H candidates are reconstructed through two- and three-body decay channels. We observe that these hypernuclei exhibit significant directed flow. Comparing to that of light nuclei, it is found that the midrapidity $v_1$ slopes of $^3_{\Lambda}$H and $^4_{\Lambda}$H follow baryon number scaling, implying that the coalescence is the dominant mechanism for these hypernuclei production in such collisions.

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.

The Beam Energy Scan (BES) program at the Relativistic Heavy Ion Collider (RHIC) was extended to energies below $\sqrt{\textit{s}_{NN}}$ = 7.7 GeV in 2015 by successful implementation of the fixed-target mode of operation in the STAR (Solenoidal Track At RHIC) experiment. In the fixed-target mode, ions circulate in one ring of the collider and interact with a stationary target at the entrance of the STAR Time Projection Chamber. The first results for Au+Au collisions at $\sqrt{\textit{s}_{NN}}$ = 4.5 GeV are presented, including directed and elliptic flow of identified hadrons, and radii from pion femtoscopy. The proton flow and pion femtoscopy results agree quantitatively with earlier measurements by Alternating Gradient Synchrotron experiments at similar energies. This validates running the STAR experiment in the fixed-target configuration. Pion directed and elliptic flow are presented for the first time at this beam energy. Pion and proton elliptic flow show behavior which hints at constituent quark scaling, but large error bars preclude reliable conclusions. The ongoing second phase of BES (BES-II) will provide fixed-target data sets with 100 times more events at each of several energies down to $\sqrt{\textit{s}_{NN}}$ = 3.0 GeV.

Measurements of mass and $\Lambda$ binding energy of $\rm ^4_{\Lambda}H$ and $\rm ^4_{\Lambda}He$ in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}=3$ GeV are presented, with an aim to address the charge symmetry breaking (CSB) problem in hypernuclei systems with atomic number A = 4. The $\Lambda$ binding energies are measured to be $\rm 2.22\pm0.06(stat.) \pm0.14(syst.)$ MeV and $\rm 2.38\pm0.13(stat.) \pm0.12(syst.)$ MeV for $\rm ^4_{\Lambda}H$ and $\rm ^4_{\Lambda}He$, respectively. The measured $\Lambda$ binding-energy difference is $\rm 0.16\pm0.14(stat.)\pm0.10(syst.)$ MeV for ground states. Combined with the $\gamma$-ray transition energies, the binding-energy difference for excited states is $\rm -0.16\pm0.14(stat.)\pm0.10(syst.)$ MeV, which is negative and comparable to the value of the ground states within uncertainties. These new measurements on the $\Lambda$ binding-energy difference in A = 4 hypernuclei systems are consistent with the theoretical calculations that result in $\rm \Delta B_{\Lambda}^4(1_{exc}^{+})\approx -\Delta B_{\Lambda}^4(0_{g.s.}^{+})<0$ and present a new method for the study of CSB effect using relativistic heavy-ion collisions.

We report the first multi-differential measurements of strange hadrons of $K^{-}$, $\phi$ and $\Xi^{-}$ yields as well as the ratios of $\phi/K^-$ and $\phi/\Xi^-$ in Au+Au collisions at ${\sqrt{s_{\rm NN}} = \rm{3\,GeV}}$ with the STAR experiment fixed target configuration at RHIC. The $\phi$ mesons and $\Xi^{-}$ hyperons are measured through hadronic decay channels, $\phi\rightarrow K^+K^-$ and $\Xi^-\rightarrow \Lambda\pi^-$. Collision centrality and rapidity dependence of the transverse momentum spectra for these strange hadrons are presented. The $4\pi$ yields and ratios are compared to thermal model and hadronic transport model predictions. At this collision energy, thermal model with grand canonical ensemble (GCE) under-predicts the $\phi/K^-$ and $\phi/\Xi^-$ ratios while the result of canonical ensemble (CE) calculations reproduce $\phi/K^-$, with the correlation length $r_c \sim 2.7$ fm, and $\phi/\Xi^-$, $r_c \sim 4.2$ fm, for the 0-10% central collisions. Hadronic transport models including high mass resonance decays could also describe the ratios. While thermal calculations with GCE work well for strangeness production in high energy collisions, the change to CE at $\rm{3\,GeV}$ implies a rather different medium property at high baryon density.

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au+Au collisions at 3.0 GeV, measured by the STAR experiment in the kinematic acceptance of rapidity ($y$) and transverse momentum ($p_{\rm T}$) within $-0.5 < y<0$ and $0.4 < p_{\rm T} <2.0 $ GeV/$c$. In the most central 0--5% collisions, a proton cumulant ratio is measured to be $C_4/C_2=-0.85 \pm 0.09 ~(\rm stat.) \pm 0.82 ~(\rm syst.)$, which is less than unity, the Poisson baseline. The hadronic transport UrQMD model reproduces our $C_4/C_2$ in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in $C_4/C_2$ is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV.

We report precision measurements of hypernuclei ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ lifetimes obtained from Au+Au collisions at \snn = 3.0 GeV and 7.2 GeV collected by the STAR experiment at RHIC, and the first measurement of ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ mid-rapidity yields in Au+Au collisions at \snn = 3.0 GeV. ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$, being the two simplest bound states composed of hyperons and nucleons, are cornerstones in the field of hypernuclear physics. Their lifetimes are measured to be $221\pm15(\rm stat.)\pm19(\rm syst.)$ ps for ${}^3_\Lambda \rm{H}$ and $218\pm6(\rm stat.)\pm13(\rm syst.)$ ps for ${}^4_\Lambda \rm{H}$. The $p_T$-integrated yields of ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ are presented in different centrality and rapidity intervals. It is observed that the shape of the rapidity distribution of ${}^4_\Lambda \rm{H}$ is different for 0--10% and 10--50% centrality collisions. Thermal model calculations, using the canonical ensemble for strangeness, describes the ${}^3_\Lambda \rm{H}$ yield well, while underestimating the ${}^4_\Lambda \rm{H}$ yield. Transport models, combining baryonic mean-field and coalescence (JAM) or utilizing dynamical cluster formation via baryonic interactions (PHQMD) for light nuclei and hypernuclei production, approximately describe the measured ${}^3_\Lambda \rm{H}$ and ${}^4_\Lambda \rm{H}$ yields. Our measurements provide means to precisely assess our understanding of the fundamental baryonic interactions with strange quarks, which can impact our understanding of more complicated systems involving hyperons, such as the interior of neutron stars or exotic hypernuclei.