Atomic nuclei are self-organized, many-body quantum systems bound by strong nuclear forces within femtometer-scale space. These complex systems manifest a variety of shapes, traditionally explored using non-invasive spectroscopic techniques at low energies. However, at these energies, their instantaneous shapes are obscured by long-timescale quantum fluctuations, making direct observation challenging. Here we introduce the ``collective flow assisted nuclear shape imaging'' method, which images the nuclear global shape by colliding them at ultrarelativistic speeds and analyzing the collective response of outgoing debris. This technique captures a collision-specific snapshot of the spatial matter distribution within the nuclei, which, through the hydrodynamic expansion, imprints patterns on the particle momentum distribution observed in detectors. We benchmark this method in collisions of ground state Uranium-238 nuclei, known for their elongated, axial-symmetric shape. Our findings show a large deformation with a slight deviation from axial symmetry in the nuclear ground state, aligning broadly with previous low-energy experiments. This approach offers a new method for imaging nuclear shapes, enhances our understanding of the initial conditions in high-energy collisions and addresses the important issue of nuclear structure evolution across energy scales.

We report the measurements of proton-deuteron ($p$-$d$) and deuteron-deuteron ($d$-$d$) correlation functions in Au+Au collisions at $\sqrt{s_\mathrm{NN}}$ = 3 GeV using fixed-target mode with the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). For the first time, the source size ($R_{G}$), scattering length ($f_{0}$), and effective range ($d_{0}$) are extracted from the measured correlation functions with a simultaneous fit. The spin-averaged $f_0$ for $p$-$d$ and $d$-$d$ interactions are determined to be -5.28 $\pm$ 0.11(stat.) $\pm$ 0.82(syst.) fm and -2.62 $\pm$ 0.02(stat.) $\pm$ 0.24(syst.) fm, respectively. The measured $p$-$d$ interaction is consistent with theoretical calculations and low-energy scattering experiment results, demonstrating the feasibility of extracting interaction parameters using the femtoscopy technique. The reasonable agreement between the experimental data and the calculations from the transport model indicates that deuteron production in these collisions is primarily governed by nucleon coalescence.

We report multi-differential measurements of strange hadron production ranging from mid- to target-rapidity in Au+Au collisions at a center-of-momentum energy per nucleon pair of $\sqrt{s_{\rm NN}}=3$ GeV with the STAR experiment at RHIC. $K^0_S$ meson and $\Lambda$ hyperon yields are measured via their weak decay channels. Collision centrality and rapidity dependences of the transverse momentum spectra and particle ratios are presented. Particle mass and centrality dependence of the average transverse momenta of $\Lambda$ and $K^0_S$ are compared with other strange particles, providing evidence of the development of hadronic rescattering in such collisions. The 4$\pi$ yields of each of these strange hadrons show a consistent centrality dependence. Discussions on radial flow, the strange hadron production mechanism, and properties of the medium created in such collisions are presented together with results from hadronic transport and thermal model calculations.

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, $v_1$ and $v_2$, of light nuclei ($d$, $t$, $^{3}$He, $^{4}$He) produced in $\sqrt{s_{\rm NN}}$ = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured $v_1$ slopes of light nuclei at mid-rapidity. For the measured $v_2$ magnitude, a strong rapidity dependence is observed. Unlike $v_2$ at higher collision energies, the $v_2$ values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

Global hyperon polarization, $\overline{P}_\mathrm{H}$, in Au+Au collisions over a large range of collision energy, $\sqrt{s_\mathrm{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While naïve extrapolation of data trends suggests a large $\overline{P}_\mathrm{H}$ as the collision energy is reduced, the behavior of $\overline{P}_\mathrm{H}$ at small $\sqrt{s_\mathrm{NN}}<7.7$ GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of $\Lambda$ hyperons along the direction of global angular momentum in Au+Au collisions at $\sqrt{s_\mathrm{NN}}=3$ GeV. The observation of substantial polarization of $4.91\pm0.81(\rm stat.)\pm0.15(\rm syst.)$% in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

We report on the measurements of directed flow $v_1$ and elliptic flow $v_2$ for hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{S}^0$, $p$, $\phi$, $\Lambda$ and $\Xi^{-}$) from Au+Au collisions at $\sqrt{s_{NN}}$ = 3 GeV and $v_{2}$ for ($\pi^{\pm}$, $K^{\pm}$, $p$ and $\overline{p}$) at 27 and 54.4 GeV with the STAR experiment. While at the two higher energy midcentral collisions the number-of-constituent-quark (NCQ) scaling holds, at 3 GeV the $v_{2}$ at midrapidity is negative for all hadrons and the NCQ scaling is absent. In addition, the $v_1$ slopes at midrapidity for almost all observed hadrons are found to be positive, implying dominant repulsive baryonic interactions. The features of negative $v_2$ and positive $v_1$ slope at 3 GeV can be reproduced with a baryonic mean-field in transport model calculations. These results imply that the medium in such collisions is likely characterized by baryonic interactions.

Notwithstanding decades of progress since Yukawa first developed a description of the force between nucleons in terms of meson exchange, a full understanding of the strong interaction remains a major challenge in modern science. One remaining difficulty arises from the non-perturbative nature of the strong force, which leads to the phenomenon of quark confinement at distances on the order of the size of the proton. Here we show that in relativistic heavy-ion collisions, where quarks and gluons are set free over an extended volume, two species of produced vector (spin-1) mesons, namely $\phi$ and $K^{*0}$, emerge with a surprising pattern of global spin alignment. In particular, the global spin alignment for $\phi$ is unexpectedly large, while that for $K^{*0}$ is consistent with zero. The observed spin-alignment pattern and magnitude for the $\phi$ cannot be explained by conventional mechanisms, while a model with a connection to strong force fields, i.e. an effective proxy description within the Standard Model and Quantum Chromodynamics, accommodates the current data. This connection, if fully established, will open a potential new avenue for studying the behaviour of strong force fields.

Partonic collectivity is one of the necessary signatures for the formation of Quark-Gluon-Plasma in high-energy nuclear collisions. Number of Constituent Quarks (NCQ) scaling has been observed for light hadron elliptic flow $v_2$ in top energy nuclear collisions at RHIC and the LHC, and this has been theoretically suggested as strong evidence for partonic collectivity. In this letter, a systematic analysis of $v_2$ of $\pi^{\pm}$, $K^{\pm}$, $K^{0}_{S}$, $p$ and $\Lambda$ in Au+Au collisions at ${\sqrt{s_{_{\rm{NN}}}}}$ = 3.2, 3.5, 3.9, and 4.5 GeV, with the STAR experiment at RHIC, is presented. NCQ scaling is markedly violated at 3.2 GeV, consistent with a hadronic-interaction dominated equation of state. However, as the collision energy increases, a gradual evolution to NCQ scaling is observed. This beam-energy dependence of $v_2$ for all hadrons studied provides evidence for the onset of dominant partonic interactions by ${\sqrt{s_{_{\rm{NN}}}}}$ = 4.5 GeV.

The chiral magnetic effect (CME) is a phenomenon that arises from the QCD anomaly in the presence of an external magnetic field. The experimental search for its evidence has been one of the key goals of the physics program of the Relativistic Heavy-Ion Collider. The STAR collaboration has previously presented the results of a blind analysis of isobar collisions (${^{96}_{44}\text{Ru}}+{^{96}_{44}\text{Ru}}$, ${^{96}_{40}\text{Zr}}+{^{96}_{40}\text{Zr}}$) in the search for the CME. The isobar ratio ($Y$) of CME-sensitive observable, charge separation scaled by elliptic anisotropy, is close to but systematically larger than the inverse multiplicity ratio, the naive background baseline. This indicates the potential existence of a CME signal and the presence of remaining nonflow background due to two- and three-particle correlations, which are different between the isobars. In this post-blind analysis, we estimate the contributions from those nonflow correlations as a background baseline to $Y$, utilizing the isobar data as well as Heavy Ion Jet Interaction Generator simulations. This baseline is found consistent with the isobar ratio measurement, and an upper limit of 10% at 95% confidence level is extracted for the CME fraction in the charge separation measurement in isobar collisions at $\sqrt{s_{\rm NN}}=200$ GeV.

Partons traversing the strongly interacting medium produced in heavy-ion collisions are expected to lose energy depending on their color charge and mass. We measure the nuclear modification factors for charm- and bottom-decay electrons, defined as the ratio of yields, scaled by the number of binary nucleon-nucleon collisions, in $\sqrt{s_{\rm NN}}$ = 200 GeV Au+Au collisions to $p$+$p$ collisions ($R_{\rm AA}$), or in central to peripheral Au+Au collisions ($R_{\rm CP}$). We find the bottom-decay electron $R_{\rm AA}$ and $R_{\rm CP}$ to be significantly higher than that of charm-decay electrons. Model calculations including mass-dependent parton energy loss in a strongly coupled medium are consistent with the measured data. These observations provide clear evidence of mass ordering of charm and bottom quark energy loss when traversing through the strongly coupled medium created in heavy-ion collisions.