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.

PHENIX has measured the azimuthal correlations of muon pairs from charm and bottom semi-leptonic decays in $p$+$p$ collisions at $\sqrt{s}=200$ GeV, using a novel analysis technique utilizing both unlike- and like-sign muon pairs to separate charm, bottom and Drell-Yan contributions. The dimuon measurements combined with the previous electron-muon and dielectron measurements span a wide range in rapidity, and are well described by PYTHIA Tune A. Through a Bayesian analysis based on PYTHIA Tune A, we show that leading order pair creation is the dominant $(76\%\pm^{14}_{19}\%)$ contribution for $b\bar{b}$ production, whereas the data favor the scenario in which next-to-leading-order processes dominate $c\bar{c}$ production. The small contribution of next-to-leading-order processes in $b\bar{b}$ production at the collision energies of the Relativistic Heavy Ion Collider contrasts with the case at Large-Hadron-Collider energies, where next-to-leading-order processes are expected to dominate.

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.

The experimental study of the collisions of heavy nuclei at relativistic energies has established the properties of the quark-gluon plasma (QGP), a state of hot, dense nuclear matter in which quarks and gluons are not bound into hadrons. In this state, matter behaves as a nearly inviscid fluid that efficiently translates initial spatial anisotropies into correlated momentum anisotropies among the produced particles, producing a common velocity field pattern known as collective flow. In recent years, comparable momentum anisotropies have been measured in small-system proton-proton ($p$$+$$p$) and proton-nucleus ($p$$+$$A$) collisions, despite expectations that the volume and lifetime of the medium produced would be too small to form a QGP. Here, we report on the observation of elliptic and triangular flow patterns of charged particles produced in proton-gold ($p$$+$Au), deuteron-gold ($d$$+$Au), and helium-gold ($^3$He$+$Au) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_{NN}}}$~=~200 GeV. The unique combination of three distinct initial geometries and two flow patterns provides unprecedented model discrimination. Hydrodynamical models, which include the formation of a short-lived QGP droplet, provide a simultaneous description of these measurements.

PHENIX reports differential cross sections of $\mu\mu$ pairs from semileptonic heavy-flavor decays and the Drell-Yan production mechanism measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV at forward and backward rapidity ($1.2<|\eta|<2.2$). The $\mu\mu$ pairs from $c\bar{c}$, $b\bar{b}$, and Drell-Yan are separated using a simultaneous fit to unlike- and like-sign muon pair spectra in mass and $p_T$. The azimuthal opening angle correlation between the muons from $c\bar{c}$ and $b\bar{b}$ decays and the pair-$p_T$ distributions are compared to distributions generated using {\sc pythia} and {\sc powheg} models, which both include next-to-leading order processes. The measured distributions for pairs from $c\bar{c}$ are consistent with {\sc pythia} calculations. The $c\bar{c}$ data presents narrower azimuthal correlations and softer $p_T$ distributions compared to distributions generated from {\sc powheg}. The $b\bar{b}$ data are well described by both models. The extrapolated total cross section for bottom production is $3.75{\pm}0.24({\rm stat}){\pm}^{0.35}_{0.50}({\rm syst}){\pm}0.45({\rm global})$[$\mu$b], which is consistent with previous measurements at the Relativistic Heavy Ion Collider in the same system at the same collision energy, and is approximately a factor of two higher than the central value calculated with theoretical models. The measured Drell-Yan cross section is in good agreement with next-to-leading-order quantum-chromodynamics calculations.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The correlations are sensitive to nonperturbative initial-state and final-state transverse momentum $k_T$ and $j_T$ in the azimuthal nearly back-to-back region $\Delta\phi\sim\pi$. In this region, transverse-momentum-dependent evolution can be studied when several different hard scales are measured. To have sensitivity to small transverse momentum scales, nonperturbative momentum widths of $p_{\rm out}$, the out-of-plane transverse momentum component perpendicular to the trigger particle, are measured. These widths are used to investigate possible effects from transverse-momentum-dependent factorization breaking. When accounting for the longitudinal momentum fraction of the away-side hadron with respect to the near-side trigger particle, the widths are found to increase with the hard scale; this is qualitatively similar to the observed behavior in Drell-Yan and semi-inclusive deep-inelastic scattering interactions. The momentum widths are also studied as a function of center-of-mass energy by comparing to previous measurements at $\sqrt{s}=510$ GeV. The nonperturbative jet widths also appear to increase with $\sqrt{s}$ at a similar $x_T$, which is qualitatively consistent to similar measurements in Drell-Yan interactions. To quantify the magnitude of any transverse-momentum-dependent factorization breaking effects, calculations will need to be performed to compare to these measurements.

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 report the transverse single-spin asymmetries of $J/\psi$ production at forward and backward rapidity, $1.2<|y|<2.2$, as a function of $J/\psi$ transverse momentum ($p_T$) and Feynman-$x$ ($x_F$). The data analyzed were recorded by the PHENIX experiment at the Relativistic Heavy Ion Collider in 2015 from $p$$+$$p$, $p$$+$Al, and $p$$+$Au collisions with transversely polarized proton beams at $\sqrt{s_{_{NN}}}=200$ GeV. At this collision energy, single-spin asymmetries for heavy-flavor particle production of $p$$+$$p$ collisions provide access to the spin-dependent gluon distribution and higher-twist correlation functions inside the nucleon, such as the gluon Qiu-Sterman and trigluon correlation functions. Proton+nucleus collisions offer an excellent opportunity to study nuclear effects on the correlation functions. The data indicate negative asymmetries at the two-standard-deviation level in the $p$$+$Au data for $p_T<2$ GeV/$c$ at both forward and backward rapidity, while in $p$$+$$p$ and $p$$+$Al collisions the asymmetries are consistent with zero within the range of experimental uncertainties.

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.

Recently, multiparticle-correlation measurements of relativistic $p/d/^3$He$+$Au, $p$$+$Pb, and even $p$$+$$p$ collisions have shown surprising collective signatures. Here we present beam-energy-scan measurements of 2-, 4-, and 6-particle angular correlations in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200, 62.4, 39, and 19.6 GeV. We also present measurements of 2- and 4-particle angular correlations in $p$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. We find the 4-particle cumulant to be real-valued for $d$$+$Au collisions at all four energies. We also find that the 4-particle cumulant in $p$$+$Au has the opposite sign as that in $d$$+$Au, indicating that the correlations are geometrical and therefore collective in origin. Further we find that the 6-particle cumulant agrees with the 4-particle cumulant in $d$$+$Au collisions at 200 GeV, which indicates that nonflow effects are subdominant and provides strong evidence of collectivity.

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 the elliptic flow, $v_2$, as a function of transverse momentum ($p_T$), pseudorapidity ($\eta$), and centrality in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200, 62.4, 39, and 19.6 GeV. The beam-energy scan of $d$$+$Au collisions provides a testing ground for the onset of flow signatures in small collision systems. We measure a nonzero $v_2$ signal at all four collision energies, which, at midrapidity and low $p_T$, is consistent with predictions from viscous hydrodynamic models. Comparisons with calculations from parton transport models (based on the {\sc ampt} Monte Carlo generator) show good agreement with the data at midrapidity to forward ($d$-going) rapidities and low $p_T$. At backward (Au-going) rapidities and $p_T>1.5$ GeV/$c$, the data diverges from {\sc ampt} calculations of $v_2$ relative to the initial geometry, indicating the possible dominance of nongeometry related correlations, referred to as nonflow. We also present measurements of the charged-particle multiplicity ($dN_{\rm ch}/d\eta$) as a function of $\eta$ in central $d$$+$Au collisions at the same energies. We find that in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=$200 GeV the $v_2$ scales with $dN_{\rm ch}/d\eta$ over all $\eta$ in the PHENIX acceptance. At $\sqrt{s_{_{NN}}}=$62.4, and 39 GeV, $v_2$ scales with $dN_{\rm ch}/d\eta$ at midrapidity and forward rapidity, but falls off at backward rapidity. This departure from the $dN_{\rm ch}/d\eta$ scaling may be a further indication of nonflow effects dominating at backward rapidity.

During 2015 the Relativistic Heavy Ion Collider (RHIC) provided collisions of transversely polarized protons with Au and Al nuclei for the first time, enabling the exploration of transverse-single-spin asymmetries with heavy nuclei. Large single-spin asymmetries in very forward neutron production have been previously observed in transversely polarized $p$$+$$p$ collisions at RHIC, and the existing theoretical framework that was successful in describing the single-spin asymmetry in $p$$+$$p$ collisions predicts only a moderate atomic-mass-number ($A$) dependence. In contrast, the asymmetries observed at RHIC in $p$$+$$A$ collisions showed a surprisingly strong $A$ dependence in inclusive forward neutron production. The observed asymmetry in $p$$+$Al collisions is much smaller, while the asymmetry in $p$$+$Au collisions is a factor of three larger in absolute value and of opposite sign. The interplay of different neutron production mechanisms is discussed as a possible explanation of the observed $A$ dependence.

We report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $p$+$p$ collisions at $\sqrt{s_{NN}}=200$~GeV. The $e^+e^-$ pair yield from $b\bar{b}$ and $c\bar{c}$ is separated by exploiting a double differential fit done simultaneously in dielectron invariant mass and $p_T$. We used three different event generators, {\sc pythia}, {\sc mc@nlo}, and {\sc powheg}, to simulate the $e^+e^-$ spectra from $c\bar{c}$ and $b\bar{b}$ production. The data can be well described by all three generators within the detector acceptance. However, when using the generators to extrapolate to $4\pi$, significant differences are observed for the total cross section. These difference are less pronounced for $b\bar{b}$ than for $c\bar{c}$. The same model dependence was observed in already published $d$+$A$ data. The $p$+$p$ data are also directly compared with $d$+$A$ data in mass and $p_T$, and within the statistical accuracy no nuclear modification is seen.

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

The cross section and transverse single-spin asymmetries of $\mu^{-}$ and $\mu^{+}$ from open heavy-flavor decays in polarized $p$+$p$ collisions at $\sqrt{s}=200$ GeV were measured by the PHENIX experiment during 2012 at the Relativistic Heavy Ion Collider. Because heavy-flavor production is dominated by gluon-gluon interactions at $\sqrt{s}=200$ GeV, these measurements offer a unique opportunity to obtain information on the trigluon correlation functions. The measurements are performed at forward and backward rapidity ($1.4<|y|<2.0$) over the transverse momentum range of $1.25<p_T<7$ GeV/$c$ for the cross section and $1.25<p_T<5$ GeV/$c$ for the asymmetry measurements. The obtained cross section is compared to a fixed-order-plus-next-to-leading-log perturbative-quantum-chromodynamics calculation. The asymmetry results are consistent with zero within uncertainties, and a model calculation based on twist-3 three-gluon correlations agrees with the data.

We report the first measurement of the full angular distribution for inclusive $J/\psi\rightarrow\mu^{+}\mu^{-}$ decays in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. The measurements are made for $J/\psi$ transverse momentum $2<p_{T}<10$ GeV/$c$ and rapidity $1.2<y<2.2$ in the Helicity, Collins-Soper, and Gottfried-Jackson reference frames. In all frames the polar coefficient $\lambda_{\theta}$ is strongly negative at low $p_{T}$ and becomes close to zero at high $p_{T}$, while the azimuthal coefficient $\lambda_{\phi}$ is close to zero at low $p_{T}$, and becomes slightly negative at higher $p_{T}$. The frame-independent coefficient $\tilde{\lambda}$ is strongly negative at all $p_{T}$ in all frames. The data are compared to the theoretical predictions provided by nonrelativistic quantum chromodynamics models.

We report the first measurement of the fraction of $J/\psi$ mesons coming from $B$-meson decay ($F_{B{\rightarrow}J/\psi}$) in $p$+$p$ collisions at $\sqrt{s}=$ 510 GeV. The measurement is performed using the forward silicon vertex detector and central vertex detector at PHENIX, which provide precise tracking and distance-of-closest-approach determinations, enabling the statistical separation of $J/\psi$ due to $B$-meson decays from prompt $J/\psi$. The measured value of $F_{B{\rightarrow}J/\psi}$ is 8.1\%$\pm$2.3\% (stat)$\pm$1.9\% (syst) for $J/\psi$ with transverse momenta $0<p_T<5$ GeV/$c$ and rapidity $1.2<|y|<2.2$. The measured fraction $F_{B{\rightarrow}J/\psi}$ at PHENIX is compared to values measured by other experiments at higher center of mass energies and to fixed-order-next-to-leading-logarithm and color-evaporation-model predictions. The $b\bar{b}$ cross section per unit rapidity ($d\sigma/dy(pp{\rightarrow}b\bar{b})$) extracted from the obtained $F_{B{\rightarrow}J/\psi}$ and the PHENIX inclusive $J/\psi$ cross section measured at 200 GeV scaled with color-evaporation-model calculations, at the mean $B$ hadron rapidity $y={\pm}1.7$ in 510 GeV $p$$+$$p$ collisions, is $3.63^{+1.92}_{-1.70}\mu$b, and it is consistent with the fixed-order-next-to-leading-logarithm calculations.