A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV, and Cu+Cu collisions sqrt(s_NN) = 22.5, 62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for non-dynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. An analysis of Negative Binomial Distribution fits to the multiplicity distributions demonstrates that the heavy ion data exhibit weak clustering properties.

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.

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.

The first measurement of two-pion Bose-Einstein correlations in central Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV at the Large Hadron Collider is presented. We observe a growing trend with energy now not only for the longitudinal and the outward but also for the sideward pion source radius. The pion homogeneity volume and the decoupling time are significantly larger than those measured at RHIC.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV. Correlations of charged hadrons of $0.7<p_T<10$ GeV/$c$ with $\pi^0$ mesons of $4<p_T<15$ GeV/$c$ or isolated direct photons of $7<p_T<15$ GeV/$c$ are used to study nonperturbative effects generated by initial-state partonic transverse momentum and final-state transverse momentum from fragmentation. The nonperturbative behavior is characterized by measuring the out-of-plane transverse momentum component $p_{\rm out}$ perpendicular to the axis of the trigger particle, which is the high-$p_T$ direct photon or $\pi^0$. Nonperturbative evolution effects are extracted from Gaussian fits to the away-side inclusive-charged-hadron yields for different trigger-particle transverse momenta ($p_T^{\rm trig}$). The Gaussian widths and root mean square of $p_{\rm out}$ are reported as a function of the interaction hard scale $p_T^{\rm trig}$ to investigate possible transverse-momentum-dependent evolution differences between the $\pi^0$-h$^\pm$ and direct photon-h$^\pm$ correlations and factorization breaking effects. The widths are found to decrease with $p_T^{\rm trig}$, which indicates that the Collins-Soper-Sterman soft factor is not driving the evolution with the hard scale in nearly back-to-back dihadron and direct photon-hadron production in $p$$+$$p$ collisions. This behavior is in contrast to Drell-Yan and semi-inclusive deep-inelastic scattering measurements.

The cross section of bottom quark-antiquark ($b\bar{b}$) production in $p$+$p$ collisions at $\sqrt{s}=510$ GeV is measured with the PHENIX detector at the Relativistic Heavy Ion Collider. The results are based on the yield of high mass, like-sign muon pairs measured within the PHENIX muon arm acceptance ($1.2<|y|<2.2$). The $b\bar{b}$ signal is extracted from like-sign dimuons by utilizing the unique properties of neutral $B$ meson oscillation. We report a differential cross section of $d\sigma_{b\bar{b}\rightarrow \mu^\pm\mu^\pm}/dy = 0.16 \pm 0.01~(\mbox{stat}) \pm 0.02~(\mbox{syst}) \pm 0.02~(\mbox{global})$ nb for like-sign muons in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $p_T>1$ GeV/$c$, and dimuon mass of 5--10 GeV/$c^2$. The extrapolated total cross section at this energy for $b\bar{b}$ production is $13.1 \pm 0.6~(\mbox{stat}) \pm 1.5~(\mbox{syst}) \pm 2.7~(\mbox{global})~\mu$b. The total cross section is compared to a perturbative quantum chromodynamics calculation and is consistent within uncertainties. The azimuthal opening angle between muon pairs from $b\bar{b}$ decays and their $p_T$ distributions are compared to distributions generated using {\sc ps pythia 6}, which includes next-to-leading order processes. The azimuthal correlations and pair $p_T$ distribution are not very well described by {\sc pythia} calculations, but are still consistent within uncertainties. Flavor creation and flavor excitation subprocesses are favored over gluon splitting.

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.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section, mean transverse momentum, mean transverse momentum squared of inclusive $J/\psi$ and cross-section ratio of $\psi(2S)$ to $J/\psi$ at forward rapidity in \pp collisions at \sqrts = 510 GeV via the dimuon decay channel. Comparison is made to inclusive $J/\psi$ cross sections measured at \sqrts = 200 GeV and 2.76--13 TeV. The result is also compared to leading-order nonrelativistic QCD calculations coupled to a color-glass-condensate description of the low-$x$ gluons in the proton at low transverse momentum ($p_T$) and to next-to-leading order nonrelativistic QCD calculations for the rest of the $p_T$ range. These calculations overestimate the data at low $p_T$. While consistent with the data within uncertainties above $\approx3$ GeV/$c$, the calculations are systematically below the data. The total cross section times the branching ratio is BR $d\sigma^{J/\psi}_{pp}/dy (1.2<|y|<2.2, 0<p_T<10~\mbox{GeV/$c$}) =$ 54.3 $\pm$ 0.5 (stat) $\pm$ 5.5 (syst) nb.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the longitudinal double spin asymmetries, $A_{LL}$, for charged pions at midrapidity ($|\eta|<0.35$) in longitudinally polarized $p+p$ collisions at $\sqrt{s}=510$ GeV. These measurements are sensitive to the gluon spin contribution to the total spin of the proton in the parton momentum fraction $x$ range between 0.04 and 0.09. One can infer the sign of the gluon polarization from the ordering of pion asymmetries with charge alone. The asymmetries are found to be consistent with global quantum-chromodynamics fits of deep-inelastic scattering and data at $\sqrt{s}=200$ GeV, which show a nonzero positive contribution of gluon spin to the proton spin.

We report the double helicity asymmetry, $A_{LL}^{J/\psi}$, in inclusive $J/\psi$ production at forward rapidity as a function of transverse momentum $p_T$ and rapidity $|y|$. The data analyzed were taken during $\sqrt{s}=510$ GeV longitudinally polarized $p