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.

The measurement of direct photons from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=39$ and 62.4 GeV in the transverse-momentum range $0.4<p_T<3$ Gev/$c$ is presented by the PHENIX collaboration at the Relativistic Heavy Ion Collider. A significant direct-photon yield is observed in both collision systems. A universal scaling is observed when the direct-photon $p_T$ spectra for different center-of-mass energies and for different centrality selections at $\sqrt{s_{_{NN}}}=62.4$ GeV is scaled with $(dN_{\rm ch}/d\eta)^{\alpha}$ for $\alpha=1.21{\pm}0.04$. This scaling also holds true for direct-photon spectra from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV measured earlier by PHENIX, as well as the spectra from Pb$+$Pb at $\sqrt{s_{_{NN}}}=2760$ GeV published by ALICE. The scaling power $\alpha$ seems to be independent of $p_T$, center of mass energy, and collision centrality. The spectra from different collision energies have a similar shape up to $p_T$ of 2 GeV/$c$. The spectra have a local inverse slope $T_{\rm eff}$ increasing with $p_T$ of $0.174\pm0.018$ GeV/$c$ in the range $0.4<p_T<1.3$ GeV/$c$ and increasing to $0.289\pm0.024$ GeV/$c$ for $0.9<p_T<2.1$ GeV/$c$. The observed similarity of low-$p_T$ direct-photon production from $\sqrt{s_{_{NN}}}= 39$ to 2760 GeV suggests a common source of direct photons for the different collision energies and event centrality selections, and suggests a comparable space-time evolution of direct-photon emission.

The PHENIX collaboration presents first measurements of low-momentum ($0.4<p_T<3$ GeV/$c$) direct-photon yields from Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=39 and 62.4 GeV. For both beam energies the direct-photon yields are substantially enhanced with respect to expectations from prompt processes, similar to the yields observed in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200. Analyzing the photon yield as a function of the experimental observable $dN_{\rm ch}/d\eta$ reveals that the low-momentum ($>$1\,GeV/$c$) direct-photon yield $dN_{\gamma}^{\rm dir}/d\eta$ is a smooth function of $dN_{\rm ch}/d\eta$ and can be well described as proportional to $(dN_{\rm ch}/d\eta)^\alpha$ with $\alpha{\sim}$1.25. This new scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and Large Hadron Collider, for centrality selected samples, as well as for different, $A$$+$$A$ collision systems. At a given beam energy the scaling also holds for high $p_T$ ($>5$\,GeV/$c$) but when results from different collision energies are compared, an additional $\sqrt{s_{_{NN}}}$-dependent multiplicative factor is needed to describe the integrated-direct-photon yield.

We report measurements of the nuclear modification factor, $R_{ \mathrm{CP}}$, for charged hadrons as well as identified $\pi^{+(-)}$, $K^{+(-)}$, and $p(\overline{p})$ for Au+Au collision energies of $\sqrt{s_{_{ \mathrm{NN}}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, and 62.4 GeV. We observe a clear high-$p_{\mathrm{T}}$ net suppression in central collisions at 62.4 GeV for charged hadrons which evolves smoothly to a large net enhancement at lower energies. This trend is driven by the evolution of the pion spectra, but is also very similar for the kaon spectra. While the magnitude of the proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ does depend on collision energy, neither the proton nor the anti-proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ exhibit net suppression at any energy. A study of how the binary collision scaled high-$p_{\mathrm{T}}$ yield evolves with centrality reveals a non-monotonic shape that is consistent with the idea that jet-quenching is increasing faster than the combined phenomena that lead to enhancement.

The extreme temperatures and energy densities generated by ultra-relativistic collisions between heavy nuclei produce a state of matter with surprising fluid properties. Non-central collisions have angular momentum on the order of 1000$\hbar$, and the resulting fluid may have a strong vortical structure that must be understood to properly describe the fluid. It is also of particular interest because the restoration of fundamental symmetries of quantum chromodynamics is expected to produce novel physical effects in the presence of strong vorticity. However, no experimental indications of fluid vorticity in heavy ion collisions have so far been found. Here we present the first measurement of an alignment between the angular momentum of a non-central collision and the spin of emitted particles, revealing that the fluid produced in heavy ion collisions is by far the most vortical system ever observed. We find that $\Lambda$ and $\overline{\Lambda}$ hyperons show a positive polarization of the order of a few percent, consistent with some hydrodynamic predictions. A previous measurement that reported a null result at higher collision energies is seen to be consistent with the trend of our new observations, though with larger statistical uncertainties. These data provide the first experimental access to the vortical structure of the "perfect fluid" created in a heavy ion collision. They should prove valuable in the development of hydrodynamic models that quantitatively connect observations to the theory of the Strong Force. Our results extend the recent discovery of hydrodynamic spin alignment to the subatomic realm.

The inclusive $J/\psi$ transverse momentum ($p_{T}$) spectra and nuclear modification factors are reported at midrapidity ($|y|<1.0$) in Au+Au collisions at $\sqrt{s_{NN}}=$ 39, 62.4 and 200 GeV taken by the STAR experiment. A suppression of $J/\psi$ production, with respect to {\color{black}the production in $p+p$ scaled by the number of binary nucleon-nucleon collisions}, is observed in central Au+Au collisions at these three energies. No significant energy dependence of nuclear modification factors is found within uncertainties. The measured nuclear modification factors can be described by model calculations that take into account both suppression of direct $J/\psi$ production due to the color screening effect and $J/\psi$ regeneration from recombination of uncorrelated charm-anticharm quark pairs.

We present results from a harmonic decomposition of two-particle azimuthal correlations measured with the STAR detector in Au+Au collisions for energies ranging from $\sqrt{s_{NN}}=7.7$ GeV to 200 GeV. The third harmonic $v_3^2\{2\}=\langle \cos3(\phi_1-\phi_2)\rangle$, where $\phi_1-\phi_2$ is the angular difference in azimuth, is studied as a function of the pseudorapidity difference between particle pairs $\Delta\eta = \eta_1-\eta_2$. Non-zero {\vthree} is directly related to the previously observed large-$\Delta\eta$ narrow-$\Delta\phi$ ridge correlations and has been shown in models to be sensitive to the existence of a low viscosity Quark Gluon Plasma (QGP) phase. For sufficiently central collisions, $v_3^2\{2\}$ persist down to an energy of 7.7 GeV suggesting that QGP may be created even in these low energy collisions. In peripheral collisions at these low energies however, $v_3^2\{2\}$ is consistent with zero. When scaled by pseudorapidity density of charged particle multiplicity per participating nucleon pair, $v_3^2\{2\}$ for central collisions shows a minimum near {\snn}$=20$ GeV.

Elliptic flow (v_2) values for identified particles at midrapidity in Au + Au collisions measured by the STAR experiment in the Beam Energy Scan at the Relativistic Heavy Ion Collider at sqrt{s_{NN}}= 7.7--62.4 GeV are presented for three centrality classes. The centrality dependence and the data at sqrt{s_{NN}}= 14.5 GeV are new. Except at the lowest beam energies we observe a similar relative v_2 baryon-meson splitting for all centrality classes which is in agreement within 15% with the number-of-constituent quark scaling. The larger v_2 for most particles relative to antiparticles, already observed for minimum bias collisions, shows a clear centrality dependence, with the largest difference for the most central collisions. Also, the results are compared with A Multiphase Transport Model and fit with a Blast Wave model.

Measurements of the fractional momentum loss ($S_{\rm loss}\equiv{\delta}p_T/p_T$) of high-transverse-momentum-identified hadrons in heavy ion collisions are presented. Using $\pi^0$ in Au$+$Au and Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=62.4$ and 200 GeV measured by the PHENIX experiment at the Relativistic Heavy Ion Collider and and charged hadrons in Pb$+$Pb collisions measured by the ALICE experiment at the Large Hadron Collider, we studied the scaling properties of $S_{\rm loss}$ as a function of a number of variables: the number of participants, $N_{\rm part}$, the number of quark participants, $N_{\rm qp}$, the charged-particle density, $dN_{\rm ch}/d\eta$, and the Bjorken energy density times the equilibration time, $\varepsilon_{\rm Bj}\tau_{0}$. We find that the $p_T$ where $S_{\rm loss}$ has its maximum, varies both with centrality and collision energy. Above the maximum, $S_{\rm loss}$ tends to follow a power-law function with all four scaling variables. The data at $\sqrt{s_{_{NN}}}$=200 GeV and 2.76 TeV, for sufficiently high particle densities, have a common scaling of $S_{\rm loss}$ with $dN_{\rm ch}/d\eta$ and $\varepsilon_{\rm Bj}\tau_{0}$, lending insight on the physics of parton energy loss.

Balance functions have been measured in terms of relative pseudorapidity ($\Delta \eta$) for charged particle pairs at the Relativistic Heavy-Ion Collider (RHIC) from Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV to 200 GeV using the STAR detector. These results are compared with balance functions measured at the Large Hadron Collider (LHC) from Pb+Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV by the ALICE Collaboration. The width of the balance function decreases as the collisions become more central and as the beam energy is increased. In contrast, the widths of the balance functions calculated using shuffled events show little dependence on centrality or beam energy and are larger than the observed widths. Balance function widths calculated using events generated by UrQMD are wider than the measured widths in central collisions and show little centrality dependence. The measured widths of the balance functions in central collisions are consistent with the delayed hadronization of a deconfined quark gluon plasma (QGP). The narrowing of the balance function in central collisions at $\sqrt{s_{\rm NN}}$ = 7.7 GeV implies that a QGP is still being created at this relatively low energy.