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.

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.

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.

We present measurements of three-dimensional correlation functions of like-sign low transverse momentum kaon pairs from sqrt(sNN)=200 GeV Au+Au collisions. A Cartesian surface-spherical harmonic decomposition technique was used to extract the kaon source function. The latter was found to have a three-dimensional Gaussian shape and can be adequately reproduced by Therminator event generator simulations with resonance contributions taken into account. Compared to the pion one, the kaon source function is generally narrower and does not have the long tail along the pair transverse momentum direction. The kaon Gaussian radii display a monotonic decrease with increasing transverse mass m_T over the interval of 0.55<=m_T<=1.15 GeV/c^2. While the kaon radii are adequately described by the m_T-scaling in the outward and sideward directions, in the longitudinal direction the lowest m_T value exceeds the expectations from a pure hydrodynamical model prediction.

We present three-particle mixed-harmonic correlations $\la \cos (m\phi_a + n\phi_b - (m+n) \phi_c)\ra$ for harmonics $m,n=1-3$ for charged particles in $\sqrt{s_{NN}}=$200 GeV Au+Au collisions at RHIC. These measurements provide information on the three-dimensional structure of the initial collision zone and are important for constraining models of a subsequent low-viscosity quark-gluon plasma expansion phase. We investigate correlations between the first, second and third harmonics predicted as a consequence of fluctuations in the initial state. The dependence of the correlations on the pseudorapidity separation between particles show hints of a breaking of longitudinal invariance. We compare our results to a number of state-of-the art hydrodynamic calculations with different initial states and temperature dependent viscosities. These measurements provide important steps towards constraining the temperature dependent transport and the longitudinal structure of the initial state at RHIC.

The measurement of J/{psi} azimuthal anisotropy is presented as a function of transverse momentum for different centralities in Au+Au collisions at {sqrt{s_{NN}}} = 200 GeV. The measured J/{psi} elliptic flow is consistent with zero within errors for transverse momentum between 2 and 10 GeV/c. Our measurement suggests that J/{psi} with relatively large transverse momentum are not dominantly produced by coalescence from thermalized charm quarks, when comparing to model calculations.

Measurements of the elliptic flow, $v_{2}$, of identified hadrons ($\pi^{\pm}$, $K^{\pm}$, $K_{s}^{0}$, $p$, $\bar{p}$, $\phi$, $\Lambda$, $\bar{\Lambda}$, $\Xi^{-}$, $\bar{\Xi}^{+}$, $\Omega^{-}$, $\bar{\Omega}^{+}$) in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, 39 and 62.4 GeV are presented. The measurements were done at mid-rapidity using the Time Projection Chamber and the Time-of-Flight detectors of the STAR experiment during the Beam Energy Scan program at RHIC. A significant difference in the $v_{2}$ values for particles and the corresponding anti-particles was observed at all transverse momenta for the first time. The difference increases with decreasing center-of-mass energy, $\sqrt{s_{NN}}$ (or increasing baryon chemical potential, $\mu_{B}$) and is larger for the baryons as compared to the mesons. This implies that particles and anti-particles are no longer consistent with the universal number-of-constituent quark (NCQ) scaling of $v_{2}$ that was observed at $\sqrt{s_{NN}}=$ 200 GeV. However, for the group of particles NCQ scaling at $(m_{T}-m_{0})/n_{q}>$ 0.4 GeV/$c^{2}$ is not violated within $\pm$10%. The $v_{2}$ values for $\phi$ mesons at 7.7 and 11.5 GeV are approximately two standard deviations from the trend defined by the other hadrons at the highest measured $p_{T}$ values.

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.

We report measurements of single- and double- spin asymmetries for $W^{\pm}$ and $Z/\gamma^*$ boson production in longitudinally polarized $p+p$ collisions at $\sqrt{s} = 510$ GeV by the STAR experiment at RHIC. The asymmetries for $W^{\pm}$ were measured as a function of the decay lepton pseudorapidity, which provides a theoretically clean probe of the proton's polarized quark distributions at the scale of the $W$ mass. The results are compared to theoretical predictions, constrained by recent polarized deep inelastic scattering measurements, and show a preference for a sizable, positive up antiquark polarization in the range $0.05<x<0.2$.

We present measurements of bulk properties of the matter produced in Au+Au collisions at $\sqrt{s_{NN}}=$ 7.7, 11.5, 19.6, 27, and 39 GeV using identified hadrons ($\pi^\pm$, $K^\pm$, $p$ and $\bar{p}$) from the STAR experiment in the Beam Energy Scan (BES) Program at the Relativistic Heavy Ion Collider (RHIC). Midrapidity ($|y|<$0.1) results for multiplicity densities $dN/dy$, average transverse momenta $\langle p_T \rangle$ and particle ratios are presented. The chemical and kinetic freeze-out dynamics at these energies are discussed and presented as a function of collision centrality and energy. These results constitute the systematic measurements of bulk properties of matter formed in heavy-ion collisions over a broad range of energy (or baryon chemical potential) at RHIC.

We report the first measurements of the moments -- mean ($M$), variance ($\sigma^{2}$), skewness ($S$) and kurtosis ($\kappa$) -- of the net-charge multiplicity distributions at mid-rapidity in Au+Au collisions at seven energies, ranging from $\sqrt {{s_{\rm NN}}}$= 7.7 to 200 GeV, as a part of the Beam Energy Scan program at RHIC. The moments are related to the thermodynamic susceptibilities of net-charge, and are sensitive to the proximity of the QCD critical point. We compare the products of the moments, $\sigma^{2}/M$, $S\sigma$ and $\kappa\sigma^{2}$ with the expectations from Poisson and negative binomial distributions (NBD). The $S\sigma$ values deviate from Poisson and are close to NBD baseline, while the $\kappa\sigma^{2}$ values tend to lie between the two. Within the present uncertainties, our data do not show non-monotonic behavior as a function of collision energy. These measurements provide a distinct way of determining the freeze-out parameters in heavy-ion collisions by comparing with theoretical models.

We report the first measurement of the elliptic anisotropy ($v_2$) of the charm meson $D^0$ at mid-rapidity ($|y|$\,$<$\,1) in Au+Au collisions at \sNN = 200\,GeV. The measurement was conducted by the STAR experiment at RHIC utilizing a new high-resolution silicon tracker. The measured $D^0$ $v_2$ in 0--80\% centrality Au+Au collisions can be described by a viscous hydrodynamic calculation for transverse momentum ($p_{\rm T}$) less than 4\,GeV/$c$. The $D^0$ $v_2$ as a function of transverse kinetic energy ($m_{\rm T} - m_0$, where $m_{\rm T} = \sqrt{p_{\rm T}^2 + m_0^2}$) is consistent with that of light mesons in 10--40\% centrality Au+Au collisions. These results suggest that charm quarks have achieved local thermal equilibrium with the medium created in such collisions. Several theoretical models, with the temperature--dependent, dimensionless charm spatial diffusion coefficient ($2{\pi}TD_s$) in the range of $\sim$2--12, are able to simultaneously reproduce our $D^0$ $v_2$ result and our previously published results for the $D^0$ nuclear modification factor.

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.

Flow harmonics ($\textit{v}_{n}$) of the Fourier expansion for the azimuthal distributions of hadrons are commonly employed to quantify the azimuthal anisotropy of particle production relative to the collision symmetry planes. While lower order Fourier coefficients ($\textit{v}_{2}$ and $\textit{v}_{3}$) are more directly related to the corresponding eccentricities of the initial state, the higher-order flow harmonics ($\textit{v}_{n>3}$) can be induced by a mode-coupled response to the lower-order anisotropies, in addition to a linear response to the same-order anisotropies. These higher-order flow harmonics and their linear and mode-coupled contributions can be used to more precisely constrain the initial conditions and the transport properties of the medium in theoretical models. The multiparticle azimuthal cumulant method is used to measure the linear and mode-coupled contributions in the higher-order anisotropic flow, the mode-coupled response coefficients, and the correlations of the event plane angles for charged particles as functions of centrality and transverse momentum in Au+Au collisions at nucleon-nucleon center-of-mass energy \roots = 200 GeV. The results are compared to similar LHC measurements as well as to several viscous hydrodynamic calculations with varying initial conditions.

Global polarization of $\Xi$ and $\Omega$ hyperons has been measured for the first time in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The measurements of the $\Xi^-$ and $\bar{\Xi}^+$ hyperon polarization have been performed by two independent methods, via analysis of the angular distribution of the daughter particles in the parity violating weak decay $\Xi\rightarrow\Lambda+\pi$, as well as by measuring the polarization of the daughter $\Lambda$-hyperon, polarized via polarization transfer from its parent. The polarization, obtained by combining the results from the two methods and averaged over $\Xi^-$ and $\bar{\Xi}^+$, is measured to be $\langle P_\Xi \rangle = 0.47\pm0.10~({\rm stat.})\pm0.23~({\rm syst.})\,\%$ for the collision centrality 20%-80%. The $\langle P_\Xi \rangle$ is found to be slightly larger than the inclusive $\Lambda$ polarization and in reasonable agreement with a multi-phase transport model (AMPT). The $\langle P_\Xi \rangle$ is found to follow the centrality dependence of the vorticity predicted in the model, increasing toward more peripheral collisions. The global polarization of $\Omega$, $\langle P_\Omega \rangle = 1.11\pm0.87~({\rm stat.})\pm1.97~({\rm syst.})\,\%$ was obtained by measuring the polarization of daughter $\Lambda$ in the decay $\Omega \rightarrow \Lambda + K$, assuming the polarization transfer factor $C_{\Omega\Lambda}=1$.

Event-by-event fluctuations of the ratio of inclusive charged to photon multiplicities at forward rapidity in Au+Au collision at $\sqrt{s_{NN}}$=200 GeV have been studied. Dominant contribution to such fluctuations is expected to come from correlated production of charged and neutral pions. We search for evidences of dynamical fluctuations of different physical origins. Observables constructed out of moments of multiplicities are used as measures of fluctuations. Mixed events and model calculations are used as baselines. Results are compared to the dynamical net-charge fluctuations measured in the same acceptance. A non-zero statistically significant signal of dynamical fluctuations is observed in excess to the model prediction when charged particles and photons are measured in the same acceptance. We find that, unlike dynamical net-charge fluctuation, charge-neutral fluctuation is not dominated by correlation due to particle decay. Results are compared to the expectations based on the generic production mechanism of pions due to isospin symmetry, for which no significant (<1%) deviation is observed.

The STAR collaboration presents for the first time two-dimensional di-hadron correlations with identified leading hadrons in 200 GeV central Au+Au and minimum-bias d+Au collisions to explore hadronization mechanisms in the quark gluon plasma. The enhancement of the jet-like yield for leading pions in Au+Au data with respect to the d+Au reference and the absence of such an enhancement for leading non-pions (protons and kaons) are discussed within the context of a quark recombination scenario. The correlated yield at large angles, specifically in the \emph{ridge region}, is found to be significantly higher for leading non-pions than pions. The consistencies of the constituent quark scaling, azimuthal harmonic model and a mini-jet modification model description of the data are tested, providing further constraints on hadronization.

We report on the first measurement of charm-strange meson $D_s^{\pm}$ production at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}}$ = 200 GeV from the STAR experiment. The yield ratio between strange ($D_{s}^{\pm}$) and non-strange ($D^{0}$) open-charm mesons is presented and compared to model calculations. A significant enhancement, relative to a PYTHIA simulation of $p$+$p$ collisions, is observed in the $D_{s}^{\pm}/D^0$ yield ratio in Au+Au collisions over a large range of collision centralities. Model calculations incorporating abundant strange-quark production in the quark-gluon plasma (QGP) and coalescence hadronization qualitatively reproduce the data. The transverse-momentum integrated yield ratio of $D_{s}^{\pm}/D^0$ at midrapidity is consistent with a prediction from a statistical hadronization model with the parameters constrained by the yields of light and strange hadrons measured at the same collision energy. These results suggest that the coalescence of charm quarks with strange quarks in the QGP plays an important role in $D_{s}^{\pm}$ meson production in heavy-ion collisions.

Quark interactions with topological gluon configurations can induce chirality imbalance and local parity violation in quantum chromodynamics. This can lead to electric charge separation along the strong magnetic field in relativistic heavy-ion collisions -- the chiral magnetic effect (CME). We report measurements by the STAR collaboration of a CME-sensitive observable in $p$+Au and $d$+Au collisions at 200 GeV, where the CME is not expected, using charge-dependent pair correlations relative to a third particle. We observe strong charge-dependent correlations similar to those measured in heavy-ion collisions. This bears important implications for the interpretation of the heavy-ion data.

In this erratum we report changes on the $D^0$ $p_T$ spectra and nuclear modification factor ($R_{AA}$) in Au+Au collisions at $\sqrt{s_{_{\mathrm{NN}}}}$ = 200 GeV by fixing the errors in the efficiency and selection criteria that affected the Au+Au results. The p+p reference spectrum has changed as well and is updated with new fragmentation parameters.

We report the first measurement of the inclusive jet and the dijet longitudinal double-spin asymmetries, $A_{LL}$, at midrapidity in polarized $pp$ collisions at a center-of-mass energy $\sqrt{s} = 510$ GeV. The inclusive jet $A_{LL}$ measurement is sensitive to the gluon helicity distribution down to a gluon momentum fraction of $x\approx 0.015$, while the dijet measurements, separated into four jet-pair topologies, provide constraints on the $x$ dependence of the gluon polarization. Both results are consistent with previous measurements made at $\sqrt{s}= 200$ GeV in the overlapping kinematic region, $x > 0.05$, and show good agreement with predictions from recent next-to-leading order global analyses.

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.

We present measurements of $\pi^-$ and $\pi^+$ elliptic flow, $v_2$, at midrapidity in Au+Au collisions at $\sqrt{s_{_{\rm NN}}} =$ 200, 62.4, 39, 27, 19.6, 11.5 and 7.7 GeV, as a function of event-by-event charge asymmetry, $A_{ch}$, based on data from the STAR experiment at RHIC. We find that $\pi^-$ ($\pi^+$) elliptic flow linearly increases (decreases) with charge asymmetry for most centrality bins at $\sqrt{s_{_{\rm NN}}} = \text{27 GeV}$ and higher. At $\sqrt{s_{_{\rm NN}}} = \text{200 GeV}$, the slope of the difference of $v_2$ between $\pi^-$ and $\pi^+$ as a function of $A_{ch}$ exhibits a centrality dependence, which is qualitatively similar to calculations that incorporate a chiral magnetic wave effect. Similar centrality dependence is also observed at lower energies.

Dihadron correlations are analyzed in $\sqrt{s_{_{\rm NN}}} = 200$ GeV $d$+Au collisions classified by forward charged particle multiplicity and zero-degree neutral energy in the Au-beam direction. It is found that the jetlike correlated yield increases with the event multiplicity. After taking into account this dependence, the non-jet contribution on the away side is minimal, leaving little room for a back-to-back ridge in these collisions.

Elliptic flow ($v_{2}$) values for identified particles at mid-rapidity in Au+Au collisions, measured by the STAR experiment in the Beam Energy Scan at RHIC at $\sqrt{s_{NN}}=$ 7.7--62.4 GeV, are presented. A beam-energy dependent difference of the values of $v_{2}$ between particles and corresponding anti-particles was observed. The difference increases with decreasing beam energy and is larger for baryons compared to mesons. This implies that, at lower energies, particles and anti-particles are not consistent with the universal number-of-constituent-quark (NCQ) scaling of $v_{2}$ that was observed at $\sqrt{s_{NN}}=$ 200 GeV.

We present results of analyses of two-pion interferometry in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27, 39, 62.4 and 200 GeV measured in the STAR detector as part of the RHIC Beam Energy Scan program. The extracted correlation lengths (HBT radii) are studied as a function of beam energy, azimuthal angle relative to the reaction plane, centrality, and transverse mass ($m_{T}$) of the particles. The azimuthal analysis allows extraction of the eccentricity of the entire fireball at kinetic freeze-out. The energy dependence of this observable is expected to be sensitive to changes in the equation of state. A new global fit method is studied as an alternate method to directly measure the parameters in the azimuthal analysis. The eccentricity shows a monotonic decrease with beam energy that is qualitatively consistent with the trend from all model predictions and quantitatively consistent with a hadronic transport model.

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.

The differential cross section and spin asymmetries for neutral pions produced within the intermediate pseudorapidity range 0.8 < {\eta} < 2.0 in polarized proton-proton collisions at sqrt{s} = 200 GeV are presented. Neutral pions were detected using the endcap electromagnetic calorimeter in the STAR detector at RHIC. The cross section was measured over a transverse momentum range of 5 < p_T < 16 GeV/c and is found to be within the scale uncertainty of a next-to-leading order perturbative QCD calculation. The longitudinal double-spin asymmetry, A_LL, is measured in the same pseudorapidity range. This quantity is sensitive to the gluonic contribution to the proton spin, {\Delta}g(x), at low Bjorken-x (down to x approx 0.01), where it is less constrained by measurements at central pseudorapidity. The measured A_LL is consistent with model predictions. The parity-violating asymmetry, A_L, is also measured and found to be consistent with zero. The transverse single-spin asymmetry, A_N, is measured within a previously unexplored kinematic range in Feynman-x and p_T. Such measurements may aid our understanding of the on-set and kinematic dependence of the large asymmetries observed at more forward pseudorapidity ({\eta} approx 3) and their underlying mechanisms. The A_N results presented are consistent with a twist-3 model prediction of a small asymmetry within the present kinematic range.