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 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.

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.

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.

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.

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.

We report the differential yields at mid-rapidity of the Breit-Wheeler process ($\gamma\gamma\rightarrow e^{+}e^{-}$) in peripheral Au+Au collisions at $\sqrt{s_{_{\rm{NN}}}} = $ 54.4 GeV and 200 GeV with the STAR experiment at RHIC, as a function of energy $\sqrt{s_{_{\rm{NN}}}}$, $e^{+}e^{-}$ transverse momentum $p_{\rm T}$, $p_{\rm T}^{2}$, invariant mass $M_{ee}$ and azimuthal angle. In the invariant mass range of 0.4 $<$$M_{ee}$$<$ 2.6 GeV/$c^{2}$ at low transverse momentum ($p_{\rm T}$$ < $0.15 GeV/$c$), the yields increase while the pair $\sqrt{\langle p_{\rm T}^{2} \rangle}$ decreases with increasing $\sqrt{s_{_{\rm{NN}}}}$, a feature is correctly predicted by the QED calculation. The energy dependencies of the measured quantities are sensitive to the nuclear form factor, infrared divergence and photon polarization. The data are compiled and used to extract the charge radius of the Au nucleus.

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.