We present measurements of elliptic and triangular azimuthal anisotropy of charged particles detected at forward rapidity $1<|\eta|<3$ in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, as a function of centrality. The multiparticle cumulant technique is used to obtain the elliptic flow coefficients $v_2\{2\}$, $v_2\{4\}$, $v_2\{6\}$, and $v_2\{8\}$, and triangular flow coefficients $v_3\{2\}$ and $v_3\{4\}$. Using the small-variance limit, we estimate the mean and variance of the event-by-event $v_2$ distribution from $v_2\{2\}$ and $v_2\{4\}$. In a complementary analysis, we also use a folding procedure to study the distributions of $v_2$ and $v_3$ directly, extracting both the mean and variance. Implications for initial geometrical fluctuations and their translation into the final state momentum distributions are discussed.

Transverse momentum spectra and yields of hadrons are measured by the PHENIX collaboration in Au + Au collisions at sqrt(s_NN) = 130 GeV at the Relativistic Heavy Ion Collider (RHIC). The time-of-flight resolution allows identification of pions to transverse momenta of 2 GeV/c and protons and antiprotons to 4 GeV/c. The yield of pions rises approximately linearly with the number of nucleons participating in the collision, while the number of kaons, protons, and antiprotons increases more rapidly. The shape of the momentum distribution changes between peripheral and central collisions. Simultaneous analysis of all the p_T spectra indicates radial collective expansion, consistent with predictions of hydrodynamic models. Hydrodynamic analysis of the spectra shows that the expansion velocity increases with collision centrality and collision energy. This expansion boosts the particle momenta, causing the yield from soft processes to exceed that for hard to large transverse momentum, perhaps as large as 3 GeV/c.

Measurements of anisotropic flow Fourier coefficients ($v_n$) for inclusive charged particles and identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ produced at midrapidity in Cu+Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV are presented. The data were collected in 2012 by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). The particle azimuthal distributions with respect to different order symmetry planes $\Psi_n$, for $n$~=~1, 2, and 3 are studied as a function of transverse momentum $p_T$ over a broad range of collisions centralities. Mass ordering, as expected from hydrodynamic flow, is observed for all three harmonics. The charged-particle results are compared to hydrodynamical and transport model calculations. We also compare these Cu$+$Au results with those in Cu$+$Cu and Au$+$Au collisions at the same $\sqrt{s_{_{NN}}}$, and find that the $v_2$ and $v_3$, as a function of transverse momentum, follow a common scaling with $1/(\varepsilon_n N_{\rm part}^{1/3})$.

The PHENIX Collaboration has measured the ratio of the yields of $\psi(2S)$ to $\psi(1S)$ mesons produced in $p$$+$$p$, $p$$+$Al, $p$$+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV over the forward and backward rapidity intervals $1.2<|y|<2.2$. We find that the ratio in $p$$+$$p$ collisions is consistent with measurements at other collision energies. In collisions with nuclei, we find that in the forward ($p$-going or $^{3}$He-going) direction, the relative yield of $\psi(2S)$ mesons to $\psi(1S)$ mesons is consistent with the value measured in \pp collisions. However, in the backward (nucleus-going) direction, the $\psi(2S)$ is preferentially suppressed by a factor of $\sim$2. This suppression is attributed in some models to breakup of the weakly-bound $\psi(2S)$ through final state interactions with comoving particles, which have a higher density in the nucleus-going direction. These breakup effects may compete with color screening in a deconfined quark-gluon plasma to produce sequential suppression of excited quarkonia states.

We report the measurement of cumulants ($C_n, n=1\ldots4$) of the net-charge distributions measured within pseudorapidity ($|\eta|<0.35$) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=7.7-200$ GeV with the PHENIX experiment at the Relativistic Heavy Ion Collider. The ratios of cumulants (e.g. $C_1/C_2$, $C_3/C_1$) of the net-charge distributions, which can be related to volume independent susceptibility ratios, are studied as a function of centrality and energy. These quantities are important to understand the quantum-chromodynamics phase diagram and possible existence of a critical end point. The measured values are very well described by expectation from negative binomial distributions. We do not observe any nonmonotonic behavior in the ratios of the cumulants as a function of collision energy. The measured values of $C_1/C_2 = \mu/\sigma^2$ and $C_3/C_1 = S\sigma^3/\mu$ can be directly compared to lattice quantum-chromodynamics calculations and thus allow extraction of both the chemical freeze-out temperature and the baryon chemical potential at each center-of-mass energy.

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in \pp and \AuAu collisions at \sqrtsn $= 200$ GeV. The substructure observables include SoftDrop groomed momentum fraction (\zg), groomed jet radius (\rg), and subjet momentum fraction (\zsj) and opening angle (\tsj). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for \pp collisions and are compared to leading order Monte Carlo models. The subjet \tsj~distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in \AuAu collisions. We compare data from \AuAu collisions to those from \pp which are embedded in minimum-bias \AuAu events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than \zg~and \rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuum-like. We also report measurements of the differential di-jet momentum imbalance ($A_{\rm{J}}$) for jets of varying \tsj. We find no qualitative differences in energy loss signatures for varying angular scales in the range $0.1 < $\tsj $ < 0.3$, leading to the possible interpretation that energy loss in this population of high momentum di-jet pairs, is due to soft medium-induced gluon radiation from a single color-charge as it traverses the medium.

We present the first measurements of long-range angular correlations and the transverse momentum dependence of elliptic flow $v_2$ in high-multiplicity $p$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. A comparison of these results with previous measurements in high-multiplicity $d$$+$Au and $^3{\rm He}$$+$Au collisions demonstrates a relation between $v_2$ and the initial collision eccentricity $\varepsilon_2$, suggesting that the observed momentum-space azimuthal anisotropies in these small systems have a collective origin and reflect the initial geometry. Good agreement is observed between the measured $v_2$ and hydrodynamic calculations for all systems, and an argument disfavoring theoretical explanations based on momentum-space domain correlations is presented. The set of measurements presented here allows us to leverage the distinct intrinsic geometry of each of these systems to distinguish between different theoretical descriptions of the long-range correlations observed in small collision systems.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section of $\phi$(1020) meson production at forward rapidity in $p$$+$$p$ collisions at $\sqrt{s}=$510 GeV via the dimuon decay channel. The integrated cross section in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $2<p_T<7$ GeV/$c$ is $\sigma_\phi=2.79 \pm 0.20\,{\rm (stat)} \pm 0.17\,{\rm (syst)} \pm 0.34\, {\rm (norm)} \times 10^{-2}$~mb. The energy dependence of $\sigma_\phi$ ($1.2<|y|<2.2$; $2<p_T<5$ GeV/$c$) is studied using the PHENIX measurements at $\sqrt{s}=$200 and 510 GeV and the Large-Hadron-Collider measurements at $\sqrt{s}=$2.76 and 7 TeV. The experimental results are compared to various event generator predictions (pythia6, pythia8, phojet, ampt, epos3, and epos-lhc).

PHENIX measurements are presented for the cross section and double-helicity asymmetry ($A_{LL}$) in inclusive $\pi^0$ production at midrapidity from $p$$+$$p$ collisions at $\sqrt{s}=510$~GeV from data taken in 2012 and 2013 at the Relativistic Heavy Ion Collider. The next-to-leading-order perturbative-quantum-chromodynamics theory calculation is in excellent agreement with the presented cross section results. The calculation utilized parton-to-pion fragmentation functions from the recent DSS14 global analysis, which prefer a smaller gluon-to-pion fragmentation function. The $\pi^{0}A_{LL}$ results follow an increasingly positive asymmetry trend with $p_T$ and $\sqrt{s}$ with respect to the predictions and are in excellent agreement with the latest global analysis results. This analysis incorporated earlier results on $\pi^0$ and jet $A_{LL}$, and suggested a positive contribution of gluon polarization to the spin of the proton $\Delta G$ for the gluon momentum fraction range $x>0.05$. The data presented here extend to a currently unexplored region, down to $x\sim0.01$, and thus provide additional constraints on the value of $\Delta G$. The results confirm the evidence for nonzero $\Delta G$ using a different production channel in a complementary kinematic region.

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.