Fluctuations of conserved quantities such as baryon number, charge, and strangeness are sensitive to the correlation length of the hot and dense matter created in relativistic heavy-ion collisions and can be used to search for the QCD critical point. We report the first measurements of the moments of net-kaon multiplicity distributions in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4, and 200 GeV. The collision centrality and energy dependence of the mean ($M$), variance ($\sigma^2$), skewness ($S$), and kurtosis ($\kappa$) for net-kaon multiplicity distributions as well as the ratio $\sigma^2/M$ and the products $S\sigma$ and $\kappa\sigma^2$ are presented. Comparisons are made with Poisson and negative binomial baseline calculations as well as with UrQMD, a transport model (UrQMD) that does not include effects from the QCD critical point. Within current uncertainties, the net-kaon cumulant ratios appear to be monotonic as a function of collision energy.

Global polarization of $\Lambda$ hyperons has been measured to be of the order of a few tenths of a percent in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV, with no significant difference between $\Lambda$ and $\bar{\Lambda}$. These new results reveal the collision energy dependence of the global polarization together with the results previously observed at $\sqrt{s_{_{NN}}}$ = 7.7 -- 62.4 GeV and indicate noticeable vorticity of the medium created in non-central heavy-ion collisions at the highest RHIC collision energy. The signal is in rough quantitative agreement with the theoretical predictions from a hydrodynamic model and from the AMPT (A Multi-Phase Transport) model. The polarization is larger in more peripheral collisions, and depends weakly on the hyperon's transverse momentum and pseudorapidity $\eta^H$ within $|\eta^H|<1$. An indication of the polarization dependence on the event-by-event charge asymmetry is observed at the $2\sigma$ level, suggesting a possible contribution to the polarization from the axial current induced by the initial magnetic field.

We present a measurement of inclusive $J/\psi$ production at mid-rapidity ($|y|<1$) in $p+p$ collisions at a center-of-mass energy of $\sqrt{s}$ = 200 GeV with the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The differential production cross section for $J/\psi$ as a function of transverse momentum ($p_T$) for $0<p_T<14$ GeV/$c$ and the total cross section are reported and compared to calculations from the color evaporation model and the non-relativistic Quantum Chromodynamics model. The dependence of $J/\psi$ relative yields in three $p_T$ intervals on charged-particle multiplicity at mid-rapidity is measured for the first time in $p+p$ collisions at $\sqrt{s}$ = 200 GeV and compared with that measured at $\sqrt{s}$ = 7 TeV, PYTHIA8 and EPOS3 Monte Carlo generators, and the Percolation model prediction.

New measurements of directed flow for charged hadrons, characterized by the Fourier coefficient \vone, are presented for transverse momenta $\mathrm{p_T}$, and centrality intervals in Au+Au collisions recorded by the STAR experiment for the center-of-mass energy range $\mathrm{\sqrt{s_{_{NN}}}} = 7.7 - 200$ GeV. The measurements underscore the importance of momentum conservation and the characteristic dependencies on $\mathrm{\sqrt{s_{_{NN}}}}$, centrality and $\mathrm{p_T}$ are consistent with the expectations of geometric fluctuations generated in the initial stages of the collision, acting in concert with a hydrodynamic-like expansion. The centrality and $\mathrm{p_T}$ dependencies of $\mathrm{v^{even}_{1}}$, as well as an observed similarity between its excitation function and that for $\mathrm{v_3}$, could serve as constraints for initial-state models. The $\mathrm{v^{even}_{1}}$ excitation function could also provide an important supplement to the flow measurements employed for precision extraction of the temperature dependence of the specific shear viscosity.

We present measurements of three-particle correlations for various harmonics in Au+Au collisions at energies ranging from $\sqrt{s_{{\rm NN}}}=7.7$ to 200 GeV using the STAR detector. The quantity $\langle\cos(m\phi_1+n\phi_2-(m+n)\phi_3)\rangle$ is evaluated as a function of $\sqrt{s_{{\rm NN}}}$, collision centrality, transverse momentum, $p_T$, pseudo-rapidity difference, $\Delta\eta$, and harmonics ($m$ and $n$). These data provide detailed information on global event properties like the three-dimensional structure of the initial overlap region, the expansion dynamics of the matter produced in the collisions, and the transport properties of the medium. A strong dependence on $\Delta\eta$ is observed for most harmonic combinations consistent with breaking of longitudinal boost invariance. Data reveal changes with energy in the two-particle correlation functions relative to the second-harmonic event-plane and provide ways to constrain models of heavy-ion collisions over a wide range of collision energies.

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.

Global spin alignment of vector mesons is a sensitive probe of system vorticity and particle production mechanism in relativistic heavy ion collisions. The measurement of global spin alignment is gaining increasing interest and deserves careful considerations. In this paper, we lay out a few practical issues that need to be taken care of when measuring global spin alignment of vector mesons. They are, the correction for event plane resolution, reconciling measurements made with different event planes, the correction for the effect of finite acceptance in pseudorapidity, and the consideration for measuring the azimuthal angle dependence. Insights and methodologies offered in this paper will help experiments to measure the global spin alignment properly and accurately.

The transversity distribution, which describes transversely polarized quarks in transversely polarized nucleons, is a fundamental component of the spin structure of the nucleon, and is only loosely constrained by global fits to existing semi-inclusive deep inelastic scattering (SIDIS) data. In transversely polarized $p^\uparrow+p$ collisions it can be accessed using transverse polarization dependent fragmentation functions which give rise to azimuthal correlations between the polarization of the struck parton and the final state scalar mesons. This letter reports on spin dependent di-hadron correlations measured by the STAR experiment. The new dataset corresponds to 25 pb$^{-1}$ integrated luminosity of $p^\uparrow+p$ collisions at $\sqrt{s}=500$ GeV, an increase of more than a factor of ten compared to our previous measurement at $\sqrt{s}=200$ GeV. Non-zero asymmetries sensitive to transversity are observed at a $Q^2$ of several hundred GeV and are found to be consistent with the former measurement and a model calculation. %we observe consistent with the former measurement are observed.} We expect that these data will enable an extraction of transversity with comparable precision to current SIDIS datasets but at much higher momentum transfers where subleading effects are suppressed.

Rapidity-odd directed flow measurements at midrapidity are presented for $\Lambda$, $\bar{\Lambda}$, $K^\pm$, $K^0_s$ and $\phi$ at $\sqrt{s_{NN}} =$ 7.7, 11.5, 14.5, 19.6, 27, 39, 62.4 and 200 GeV in Au+Au collisions recorded by the STAR detector at the Relativistic Heavy Ion Collider. These measurements greatly expand the scope of data available to constrain models with differing prescriptions for the equation of state of quantum chromodynamics. Results show good sensitivity for testing a picture where flow is assumed to be imposed before hadron formation and the observed particles are assumed to form via coalescence of constituent quarks. The pattern of departure from a coalescence-inspired sum-rule can be a valuable new tool for probing the collision dynamics.

We report the first measurements of transverse single-spin asymmetries for inclusive jet and jet + $\pi^{\pm}$ production at midrapidity from transversely polarized proton-proton collisions at $\sqrt{s} = 500$ GeV. The data were collected in 2011 with the STAR detector sampled from 23 pb$^{-1}$ integrated luminosity with an average beam polarization of 53%. Asymmetries are reported for jets with transverse momenta $6 < p_{T, jet} < 55$ GeV/c and pseudorapidity $|\eta| < 1$. Presented are measurements of the inclusive-jet azimuthal transverse single-spin asymmetry, sensitive to twist-3 initial-state quark-gluon correlators; the Collins asymmetry, sensitive to quark transversity coupled to the polarized Collins fragmentation function; and the first measurement of the "Collins-like" asymmetry, sensitive to linearly polarized gluons. Within the present statistical precision, inclusive-jet and Collins-like asymmetries are small, with the latter allowing the first experimental constraints on gluon linear polarization in a polarized proton. At higher values of jet transverse momenta, we observe the first non-zero Collins asymmetries in polarized-proton collisions, with a statistical significance of greater than $5\sigma$. The results span a range of x similar to results from SIDIS but at much higher $Q^{2}$. The Collins results enable tests of universality and factorization-breaking in the transverse momentum-dependent formulation of perturbative quantum chromodynamics.

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

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.

We present measurements of elliptic flow ($v_2$) of electrons from the decays of heavy-flavor hadrons ($e_{HF}$) by the STAR experiment. For Au+Au collisions at $\sqrt{s_{\rm NN}} = $ 200 GeV we report $v_2$, for transverse momentum ($p_T$) between 0.2 and 7 GeV/c using three methods: the event plane method ($v_{2}${EP}), two-particle correlations ($v_2${2}), and four-particle correlations ($v_2${4}). For Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 62.4 and 39 GeV we report $v_2${2} for $p_T< 2$ GeV/c. $v_2${2} and $v_2${4} are non-zero at low and intermediate $p_T$ at 200 GeV, and $v_2${2} is consistent with zero at low $p_T$ at other energies. The $v_2${2} at the two lower beam energies is systematically lower than at $\sqrt{s_{\rm NN}} = $ 200 GeV for $p_T < 1$ GeV/c. This difference may suggest that charm quarks interact less strongly with the surrounding nuclear matter at those two lower energies compared to $\sqrt{s_{\rm NN}} = 200$ GeV.

We report the first measurement of the longitudinal double-spin asymmetry $A_{LL}$ for mid-rapidity di-jet production in polarized $pp$ collisions at a center-of-mass energy of $\sqrt{s} = 200$ GeV. The di-jet cross section was measured and is shown to be consistent with next-to-leading order (NLO) perturbative QCD predictions. $A_{LL}$ results are presented for two distinct topologies, defined by the jet pseudorapidities, and are compared to predictions from several recent NLO global analyses. The measured asymmetries, the first such correlation measurements, support those analyses that find positive gluon polarization at the level of roughly 0.2 over the region of Bjorken-$x > 0.05$.

The STAR Collaboration reports on the photoproduction of $\pi^+\pi^-$ pairs in gold-gold collisions at a center-of-mass energy of 200 GeV/nucleon-pair. These pion pairs are produced when a nearly-real photon emitted by one ion scatters from the other ion. We fit the $\pi^+\pi^-$ invariant mass spectrum with a combination of $\rho$ and $\omega$ resonances and a direct $\pi^+\pi^-$ continuum. This is the first observation of the $\omega$ in ultra-peripheral collisions, and the first measurement of $\rho-\omega$ interference at energies where photoproduction is dominated by Pomeron exchange. The $\omega$ amplitude is consistent with the measured $\gamma p\rightarrow \omega p$ cross section, a classical Glauber calculation and the $\omega\rightarrow\pi^+\pi^-$ branching ratio. The $\omega$ phase angle is similar to that observed at much lower energies, showing that the $\rho-\omega$ phase difference does not depend significantly on photon energy. The $\rho^0$ differential cross section $d\sigma/dt$ exhibits a clear diffraction pattern, compatible with scattering from a gold nucleus, with 2 minima visible. The positions of the diffractive minima agree better with the predictions of a quantum Glauber calculation that does not include nuclear shadowing than with a calculation that does include shadowing.

The STAR Collaboration reports the measurement of semi-inclusive distributions of charged-particle jets recoiling from a high transverse momentum hadron trigger, in central and peripheral Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Charged jets are reconstructed with the anti-kT algorithm for jet radii R between 0.2 and 0.5 and with low infrared cutoff of track constituents ($p_T>0.2$ GeV/c). A novel mixed-event technique is used to correct the large uncorrelated background present in heavy ion collisions. Corrected recoil jet distributions are reported at mid-rapidity, for charged-jet transverse momentum $p_T^\mathrm{jet,ch}<30$ GeV/c. Comparison is made to similar measurements for Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV, to calculations for p+p collisions at $\sqrt{s}$ = 200 GeV based on the PYTHIA Monte Carlo generator and on a Next-to-Leading Order perturbative QCD approach, and to theoretical calculations incorporating jet quenching. The recoil jet yield is suppressed in central relative to peripheral collisions, with the magnitude of the suppression corresponding to medium-induced charged energy transport out of the jet cone of $2.8\pm0.2\mathrm{(stat)}\pm1.5\mathrm{(sys)}$ GeV/c, for $10<p_T^\mathrm{jet,ch}<20$ GeV/c and R = 0.5. No medium-induced change in jet shape is observed for $R<0.5$. The azimuthal distribution of low-$p_T^\mathrm{jet,ch}$ recoil jets may be enhanced at large azimuthal angles to the trigger axis, due to scattering off quasi-particles in the hot QCD medium. Measurement of this distribution gives a 90% statistical confidence upper limit to the yield enhancement at large deflection angles in central Au+Au collisions of $50\pm30\mathrm{(sys)\%}$ of the large-angle yield in p+pcollisions predicted by PYTHIA.

This Workshop brought top experts, researchers, postdocs, and students from high-energy heavy ion interactions, lattice QCD and hadronic physics communities together. YSTAR2016 discussed the impact of "missing" hyperon resonances on QCD thermodynamics, on freeze-out in heavy ion collisions, on the evolution of early universe, and on the spectroscopy of strange particles. Recent studies that compared lattice QCD predictions of thermodynamic properties of quark-gluon plasma at freeze-out with calculations based on statistical hadron resonance gas models as well as experimentally measured ratios between yields of different hadron species in heavy ion collisions provide indirect evidence for the presence of "missing" resonances in all of these contexts. The aim of the YSTAR2016 Workshop was to sharpen these comparisons and advance our understanding of the formation of strange hadrons from quarks and gluons microseconds after the Big Bang and in todays experiments at LHC and RHIC as well as at future facilities like FAIR, J-PARC and KL at JLab. It was concluded that the new initiative to create a secondary beam of neutral kaons at JLab will make a bridge between the hardron spectroscopy, heavy-ion experiments and lattice QCD studies addressing some major issues related to thermodynamics of the early universe and cosmology in general.

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

We report on the measurement of $\rm{J}/\psi$ production in the dielectron channel at mid-rapidity (|y|<1) in p+p and d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV from the STAR experiment at the Relativistic Heavy Ion Collider. The transverse momentum $p_{T}$ spectra in p+p for $p_{T}$ < 4 GeV/c and d+Au collisions for $p_{T}$ < 3 GeV/c are presented. These measurements extend the STAR coverage for $\rm{J}/\psi$ production in p+p collisions to low $p_{T}$. The $<p_{T}^{2}>$ from the measured $\rm{J}/\psi$ invariant cross section in p+p and d+Au collisions are evaluated and compared to similar measurements at other collision energies. The nuclear modification factor for $\rm{J}/\psi$ is extracted as a function of $p_{T}$ and collision centrality in d+Au and compared to model calculations using the modified nuclear Parton Distribution Function and a final-state $\rm{J}/\psi$ nuclear absorption cross section.