Collisions between prolate uranium nuclei are used to study how particle production and azimuthal anisotropies depend on initial geometry in heavy-ion collisions. We report the two- and four-particle cumulants, $v_2\{2\}$ and $v_2\{4\}$, for charged hadrons from U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV and Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. Nearly fully overlapping collisions are selected based on the amount of energy deposited by spectators in the STAR Zero Degree Calorimeters (ZDCs). Within this sample, the observed dependence of $v_2\{2\}$ on multiplicity demonstrates that ZDC information combined with multiplicity can preferentially select different overlap configurations in U+U collisions. An initial-state model with gluon saturation describes the slope of $v_2\{2\}$ as a function of multiplicity in central collisions better than one based on Glauber with a two-component multiplicity model.

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.

We present a measurement of the inclusive production of Upsilon mesons in U+U collisions at 193 GeV at mid-rapidity (|y| < 1). Previous studies in central Au+Au collisions at 200 GeV show a suppression of Upsilon(1S+2S+3S) production relative to expectations from the Upsilon yield in p+p collisions scaled by the number of binary nucleon-nucleon collisions (Ncoll), with an indication that the Upsilon(1S) state is also suppressed. The present measurement extends the number of participant nucleons in the collision (Npart) by 20% compared to Au+Au collisions, and allows us to study a system with higher energy density. We observe a suppression in both the Upsilon(1S+2S+3S) and Upsilon(1S) yields in central U+U data, which consolidates and extends the previously observed suppression trend in Au+Au collisions.

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.

We present a measurement of the first-order azimuthal anisotropy, $v_1(y)$, of deuterons from Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV recorded with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The energy dependence of the $v_1(y)$ slope, $dv_{1}/dy|_{y=0}$, for deuterons, where $y$ is the rapidity, is extracted for semi-central collisions (10-40\% centrality) and compared to that of protons. While the $v_1(y)$ slopes of protons are generally negative for $\sqrt{s_{NN}} >$ 10 GeV, those for deuterons are consistent with zero, a strong enhancement of the $v_1(y)$ slope of deuterons is seen at the lowest collision energy (the largest baryon density) at $\sqrt{s_{NN}} =$ 7.7 GeV. In addition, we report the transverse momentum dependence of $v_1$ for protons and deuterons. The experimental results are compared with transport and coalescence models.

We present the first measurement of charge-dependent directed flow in Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The results are presented as a function of the particle transverse momentum and pseudorapidity for different centralities. A finite difference between the directed flow of positive and negative charged particles is observed that qualitatively agrees with the expectations from the effects of the initial strong electric field between two colliding ions with different nuclear charges. The measured difference in directed flow is much smaller than that obtained from the parton-hadron-string-dynamics (PHSD) model, which suggests that most of the electric charges, i.e. quarks and antiquarks, have not yet been created during the lifetime of the strong electric field, which is of the order of, or less than, 1fm/$c$.

We report the STAR measurements of dielectron ($e^+e^-$) production at midrapidity ($|y_{ee}|<$1) in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200\,GeV. The measurements are evaluated in different invariant mass regions with a focus on 0.30-0.76 ($\rho$-like), 0.76-0.80 ($\omega$-like), and 0.98-1.05 ($\phi$-like) GeV/$c^{2}$. The spectrum in the $\omega$-like and $\phi$-like regions can be well described by the hadronic cocktail simulation. In the $\rho$-like region, however, the vacuum $\rho$ spectral function cannot describe the shape of the dielectron excess. In this range, an enhancement of 1.77$\pm$0.11(stat.)$\pm$0.24(sys.)$\pm$0.33(cocktail) is determined with respect to the hadronic cocktail simulation that excludes the $\rho$ meson. The excess yield in the $\rho$-like region increases with the number of collision participants faster than the $\omega$ and $\phi$ yields. Theoretical models with broadened $\rho$ contributions through interactions with constituents in the hot QCD medium provide a consistent description of the dilepton mass spectra for the measurement presented here and the earlier data at the Super Proton Synchrotron energies.

We report on a polarization measurement of inclusive $J/\psi$ mesons in the di-electron decay channel at mid-rapidity at 2 $<p_{T}<$ 6 GeV/$c$ in $p+p$ collisions at $\sqrt{s}$ = 200 GeV. Data were taken with the STAR detector at RHIC. The $J/\psi$ polarization measurement should help to distinguish between different models of the $J/\psi$ production mechanism since they predict different $p_{T}$ dependences of the $J/\psi$ polarization. In this analysis, $J/\psi$ polarization is studied in the helicity frame. The polarization parameter $\lambda_{\theta}$ measured at RHIC becomes smaller towards high $p_{T}$, indicating more longitudinal $J/\psi$ polarization as $p_{T}$ increases. The result is compared with predictions of presently available models.

We report measurements of the third harmonic coefficient of the azimuthal anisotropy, v_3, known as triangular flow. The analysis is for charged particles in Au+Au collisions at sqrtsNN = 200 GeV, based on data from the STAR experiment at the Relativistic Heavy Ion Collider. Two-particle correlations as a function of their pseudorapidity separation are fit with narrow and wide Gaussians. Measurements of triangular flow are extracted from the wide Gaussian, from two-particle cumulants with a pseudorapidity gap, and also from event plane analysis methods with a large pseudorapidity gap between the particles and the event plane. These results are reported as a function of transverse momentum and centrality. A large dependence on the pseudorapidity gap is found. Results are compared with other experiments and model calculations.

We report the direct virtual photon invariant yields in the transverse momentum ranges $1\!<\!p_{T}\!<\!3$ GeV/$c$ and $5\!<\!p_T\!<\!10$ GeV/$c$ at mid-rapidity derived from the dielectron invariant mass continuum region $0.10<M_{ee}<0.28$ GeV/$c^{2}$ for 0-80\% minimum-bias Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV. A clear excess in the invariant yield compared to the number-of-binary-collisions ($N_{bin}$) scaled $p+p$ reference is observed in the $p_T$ range $1\!<\!p_{T}\!<\!3$ GeV/$c$. For $p_T\!>6$ GeV/$c$ the production follows $N_{bin}$ scaling. Model calculations with contributions from thermal radiation and initial hard parton scattering are consistent within uncertainties with the direct virtual photon invariant yield.

We present $\Lambda\Lambda$ correlation measurements in heavy-ion collisions for Au+Au collisions at $\sqrt{s_{NN}}= 200$ GeV using the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). The Lednick\'{y}-Lyuboshitz analytical model has been used to fit the data to obtain a source size, a scattering length and an effective range. Implications of the measurement of the $\Lambda\Lambda$ correlation function and interaction parameters for di-hyperon searches are discussed.

Azimuthal angular correlations of charged hadrons with respect to the axis of a reconstructed (trigger) jet in Au+Au and p+p collisions at $\sqrt{s_{\text{NN}}} = 200 \text{GeV}$ in STAR are presented. The trigger jet population in Au+Au collisions is biased towards jets that have not interacted with the medium, allowing easier matching of jet energies between Au+Au and p+p collisions while enhancing medium effects on the recoil jet. The associated hadron yield of the recoil jet is significantly suppressed at high transverse momentum ($p_{\text{T}}^{\text{assoc}}$) and enhanced at low $p_{\text{T}}^{\text{assoc}}$ in 0-20% central Au+Au collisions compared to p+p collisions, which is indicative of medium-induced parton energy loss in ultrarelativistic heavy-ion collisions.

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.

Rapidity-odd directed flow($v_1$) measurements for charged pions, protons and antiprotons near mid-rapidity ($y=0$) are reported in $\sqrt{s_{NN}} =$ 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV Au + Au collisions as recorded by the STAR detector at the Relativistic Heavy Ion Collider (RHIC). At intermediate impact parameters, the proton and net-proton slope parameter $dv_1/dy|_{y=0}$ shows a minimum between 11.5 and 19.6 GeV. In addition, the net-proton $dv_1/dy|_{y=0}$ changes sign twice between 7.7 and 39 GeV. The proton and net-proton results qualitatively resemble predictions of a hydrodynamic model with a first-order phase transition from hadronic matter to deconfined matter, and differ from hadronic transport calculations.

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.

A study is reported of the same- and opposite-sign charge-dependent azimuthal correlations with respect to the event plane in Au+Au collisions at 200 GeV. The charge multiplicity asymmetries between the up/down and left/right hemispheres relative to the event plane are utilized. The contributions from statistical fluctuations and detector effects were subtracted from the (co-)variance of the observed charge multiplicity asymmetries. In the mid- to most-central collisions, the same- (opposite-) sign pairs are preferentially emitted in back-to-back (aligned on the same-side) directions. The charge separation across the event plane, measured by the difference, $\Delta$, between the like- and unlike-sign up/down $-$ left/right correlations, is largest near the event plane. The difference is found to be proportional to the event-by-event final-state particle ellipticity (via the observed second-order harmonic $v^{\rm obs}_{2}$), where $\Delta=(1.3\pm1.4({\rm stat})^{+4.0}_{-1.0}({\rm syst}))\times10^{-5}+(3.2\pm0.2({\rm stat})^{+0.4}_{-0.3}({\rm syst}))\times10^{-3}v^{\rm obs}_{2}$ for 20-40% Au+Au collisions. The implications for the proposed chiral magnetic effect are discussed.

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.

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.

We report on the first measurements of J/$\psi$ production at very low transverse momentum ($p_{T} <$ 0.2 GeV/c) in hadronic Au+Au collisions at $\sqrt{s_{\rm{NN}}} =$ 200 GeV and U+U collisions at $\sqrt{s_{\rm{NN}}} =$ 193 GeV. Remarkably, the inferred nuclear modification factor of J/$\psi$ at mid-rapidity in Au+Au (U+U) collisions reaches about 24 (52) for $p_{T} <$ 0.05 GeV/c in the 60-80$\%$ collision centrality class. This noteworthy enhancement cannot be explained by hadronic production accompanied by cold and hot medium effects. In addition, the $dN/dt$ distribution of J/$\psi$ for the very low $p_{T}$ range is presented for the first time. The distribution is consistent with that expected from the Au nucleus and shows a hint of interference. Comparison of the measurements to theoretical calculations of coherent production shows that the excess yield can be described reasonably well and reveals a partial disruption of coherent production in semi-central collisions, perhaps due to the violent hadronic interactions. Incorporating theoretical calculations, the results strongly suggest that the dramatic enhancement of J/$\psi$ yield observed at extremely low $p_{T}$ originates from coherent photon-nucleus interactions. In particular, coherently produced J/$\psi$'s in violent hadronic collisions may provide a novel probe of the quark-gluon-plasma.

Dihadron angular correlations in $d$+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV are reported as a function of the measured zero-degree calorimeter neutral energy and the forward charged hadron multiplicity in the Au-beam direction. A finite correlated yield is observed at large relative pseudorapidity ($\Delta\eta$) on the near side (i.e. relative azimuth $\Delta\phi\sim0$). This correlated yield as a function of $\Delta\eta$ appears to scale with the dominant, primarily jet-related, away-side ($\Delta\phi\sim\pi$) yield. The Fourier coefficients of the $\Delta\phi$ correlation, $V_{n}=\langle\cos n\Delta\phi\rangle$, have a strong $\Delta\eta$ dependence. In addition, it is found that $V_{1}$ is approximately inversely proportional to the mid-rapidity event multiplicity, while $V_{2}$ is independent of it with similar magnitude in the forward ($d$-going) and backward (Au-going) directions.

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.

Local parity-odd domains are theorized to form inside a Quark-Gluon-Plasma (QGP) which has been produced in high-energy heavy-ion collisions. The local parity-odd domains manifest themselves as charge separation along the magnetic field axis via the chiral magnetic effect (CME). The experimental observation of charge separation has previously been reported for heavy-ion collisions at the top RHIC energies. In this paper, we present the results of the beam-energy dependence of the charge correlations in Au+Au collisions at midrapidity for center-of-mass energies of 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV from the STAR experiment. After background subtraction, the signal gradually reduces with decreased beam energy, and tends to vanish by 7.7 GeV. The implications of these results for the CME will be discussed.

We report measurements of Upsilon meson production in p+p, d+Au, and Au+Au collisions using the STAR detector at RHIC. We compare the Upsilon yield to the measured cross section in p+p collisions in order to quantify any modifications of the yield in cold nuclear matter using d+Au data and in hot nuclear matter using Au+Au data separated into three centrality classes. Our p+p measurement is based on three times the statistics of our previous result. We obtain a nuclear modification factor for Upsilon(1S+2S+3S) in the rapidity range |y|<1 in d+Au collisions of R_dAu = 0.79 +/- 0.24 (stat.) +/- 0.03 (sys.) +/- 0.10 (pp sys.). A comparison with models including shadowing and initial state parton energy loss indicates the presence of additional cold-nuclear matter suppression. Similarly, in the top 10% most-central Au+Au collisions, we measure a nuclear modification factor of R_AA=0.49 +/- 0.1 (stat.) +/- 0.02 (sys.) +/- 0.06 (pp sys.), which is a larger suppression factor than that seen in cold nuclear matter. Our results are consistent with complete suppression of excited-state Upsilon mesons in Au+Au collisions. The additional suppression in Au+Au is consistent with the level expected in model calculations that include the presence of a hot, deconfined Quark-Gluon Plasma. However, understanding the suppression seen in d+Au is still needed before any definitive statements about the nature of the suppression in Au+Au can be made.

We report on the first measurement of the azimuthal anisotropy ($v_2$) of dielectrons ($e^{+}e^{-}$ pairs) at mid-rapidity from $\sqrt{s_{_{NN}}} = 200$ GeV Au+Au collisions with the STAR detector at RHIC, presented as a function of transverse momentum ($p_T$) for different invariant-mass regions. In the mass region $M_{ee}\!<1.1$ GeV/$c^2$ the dielectron $v_2$ measurements are found to be consistent with expectations from $\pi^{0}$, $\eta$, $\omega$ and $\phi$ decay contributions. In the mass region $1.1\!<M_{ee}\!<2.9$ GeV/$c^2$, the measured dielectron $v_2$ is consistent, within experimental uncertainties, with that from the $c\bar{c}$ contributions.