Yields for J/psi production in Cu+Cu collisions at sqrt (s_NN)= 200 GeV have been measured by the PHENIX experiment over the rapidity range |y| < 2.2 at transverse momenta from 0 to beyond 5 GeV/c. The invariant yield is obtained as a function of rapidity, transverse momentum and collision centrality, and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data provide greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, providing a key constraint that is needed for disentangling cold and hot nuclear matter effects.

The production of $\Upsilon (nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{-1}$. The $\Upsilon (nS)$ mesons are reconstructed through their decays into two opposite-sign muons. The measurements comprise the differential production cross-sections of the $\Upsilon (1S)$ and $\Upsilon (2S)$ states, their forward-to-backward ratios and nuclear modification factors, performed as a function of the transverse momentum $p_{\mathrm{T}}$ and rapidity in the nucleon-nucleon centre-of-mass frame $y^*$ of the $\Upsilon (nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($-5.0<y^*<-2.5$) for $p$Pb (Pb$p$) collisions. In addition, production cross-sections for $\Upsilon (3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon (nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.

The cross-section for prompt antiproton production in collisions of protons with an energy of $6.5$ TeV incident on helium nuclei at rest is measured with the LHCb experiment from a data set corresponding to an integrated luminosity of $0.5\,nb^{-1}$. The target is provided by injecting helium gas into the LHC beam line at the LHCb interaction point. The reported results, covering antiproton momenta between $12$ and $110\,\mathrm{GeV/}c$, represent the first direct determination of the antiproton production cross-section in ${\rm p He}$ collisions, and impact the interpretation of recent results on antiproton cosmic rays from space-borne experiments.

The first measurement of the $Z$ boson production cross-section at centre-of-mass energy $\sqrt{s} = 5.02\,$TeV in the forward region is reported, using $pp$ collision data collected by the LHCb experiment in year 2017, corresponding to an integrated luminosity of $100 \pm 2\,\rm{pb^{-1}}$. The production cross-section is measured for final-state muons in the pseudorapidity range $2.0<\eta<4.5$ with transverse momentum $p_{\rm{T}}> 20\,\rm{GeV/}\it{c}$. The integrated cross-section is determined to be \[ \sigma_{Z \rightarrow \mu^{+}\mu^{-}} = 39.6 \pm 0.7\,(\rm{stat}) \pm 0.6\,(\rm{syst}) \pm 0.8\,(\rm{lumi}) \ \rm{pb} \] for the di-muon invariant mass in the range $60<M_{\mu\mu}<120\,\rm{GeV/}\it{c^{2}}$. This result and the differential cross-section results are in good agreement with theoretical predictions at next-to-next-to-leading order in the strong coupling. Based on a previous LHCb measurement of the $Z$ boson production cross-section in $p$Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV, the nuclear modification factor $R_{p\rm{Pb}}$ is measured for the first time at this energy. The measured values are $1.2^{+0.5}_{-0.3}\,(\rm{stat}) \pm 0.1\,(\rm{syst})$ in the forward region ($1.53<y^*_{\mu}<4.03$) and $3.6^{+1.6}_{-0.9}\,(\rm{stat}) \pm 0.2\,(\rm{syst})$ in the backward region ($-4.97<y^*_{\mu}<-2.47$), where $y^*_{\mu}$ represents the muon rapidity in the centre-of-mass frame.

Using a data sample corresponding to an integrated luminosity of $2.0\,fb^{-1}$, collected by the LHCb experiment, the production of the $\eta_c(1S)$ state in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13 \text{ TeV}$ is studied in the rapidity range ${2.0 < y < 4.5}$ and in the transverse momentum range ${6.5 < p_{T} < 14.0\text{ GeV}}$. The cross-section for prompt production of $\eta_c(1S)$ mesons relative to that of the $J/\psi$ meson is measured using the ${p\bar{p}}$ decay mode and is found to be ${\sigma_{\eta_c(1S)}/\sigma_{J/\psi} = 1.69 \pm 0.15 \pm 0.10 \pm 0.18}$. The quoted uncertainties are, in order, statistical, systematic and due to uncertainties on the branching fractions of the ${J/\psi\to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The prompt $\eta_c(1S)$ production cross-section is determined to be ${\sigma_{\eta_c(1S)} = 1.26 \pm 0.11\pm 0.08 \pm 0.14 \,\mu b}$, where the last uncertainty includes that on the ${J/\psi}$ meson cross-section. The ratio of the branching fractions of $b$-hadron decays to the $\eta_c(1S)$ and ${J/\psi}$ states is measured to be ${\mathcal{B}_{b\to\eta_c X}/\mathcal{B}_{b\to J/\psi X} = 0.48 \pm 0.03 \pm 0.03 \pm 0.05}$, where the last uncertainty is due to those on the branching fractions of the ${J/\psi \to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The difference between the ${J/\psi}$ and $\eta_c(1S)$ masses is also determined to be ${113.0 \pm 0.7 \pm 0.1\text{ MeV}}$, which is the most precise single measurement of this quantity to date.

A three-dimensional (3D) correlation function obtained from mid-rapidity, low pT pion pairs in central Au+Au collisions at sqrt(s_NN)=200 GeV is studied. The extracted model-independent source function indicates a long range tail in the directions of the pion pair transverse momentum (out) and the beam (long). Model comparisons to these distensions indicate a proper breakup time \tau_0 ~ 9 fm/c and a mean proper emission duration \Delta\tau ~ 2 fm/c, leading to sizable emission time differences (<|\Delta \tau_LCM |> ~ 12 fm/c), partly due to resonance decays. They also suggest an outside-in 'burning' of the emission source reminiscent of many hydrodynamical models.

The differential cross-section of prompt inclusive production of long-lived charged particles in proton-proton collisions is measured using a data sample recorded by the LHCb experiment at a centre-of-mass energy of ${\sqrt{s} = 13\,\mathrm{TeV}}$. The data sample, collected with an unbiased trigger, corresponds to an integrated luminosity of ${5.4\,\mathrm{nb}^{-1}}$. The differential cross-section is measured as a function of transverse momentum and pseudorapidity in the ranges ${p_\mathrm{T} \in [0.08, 10)\,\mathrm{GeV}\,c^{-1}}$ and ${\eta \in [2.0, 4.8)}$ and is determined separately for positively and negatively charged particles. The results are compared with predictions from various hadronic-interaction models.

The inclusive $b \bar{b}$- and $c \bar{c}$-dijet production cross-sections in the forward region of $pp$ collisions are measured using a data sample collected with the LHCb detector at a centre-of-mass energy of 13 TeV in 2016. The data sample corresponds to an integrated luminosity of 1.6 fb$^{-1}$. Differential cross-sections are measured as a function of the transverse momentum and of the pseudorapidity of the leading jet, of the rapidity difference between the jets, and of the dijet invariant mass. A fiducial region for the measurement is defined by requiring that the two jets originating from the two $b$ or $c$ quarks are emitted with transverse momentum greater than 20 GeV$/c$, pseudorapidity in the range $2.2 < \eta < 4.2$, and with a difference in the azimuthal angle between the two jets greater than 1.5. The integrated $b \bar{b}$-dijet cross-section is measured to be $53.0 \pm 9.7$ nb, and the total $c \bar{c}$-dijet cross-section is measured to be $73 \pm 16$ nb. The ratio between $c \bar{c}$- and $b \bar{b}$-dijet cross-sections is also measured and found to be $1.37 \pm 0.27$. The results are in agreement with theoretical predictions at next-to-leading order.

The PHENIX experiment has measured the suppression of semi-inclusive single high transverse momentum pi^0's in Au+Au collisions at sqrt(s_NN) = 200 GeV. The present understanding of this suppression is in terms of energy-loss of the parent (fragmenting) parton in a dense color-charge medium. We have performed a quantitative comparison between various parton energy-loss models and our experimental data. The statistical point-to-point uncorrelated as well as correlated systematic uncertainties are taken into account in the comparison. We detail this methodology and the resulting constraint on the model parameters, such as the initial color-charge density dN^g/dy, the medium transport coefficient <q^hat>, or the initial energy-loss parameter epsilon_0. We find that high transverse momentum pi^0 suppression in Au+Au collisions has sufficient precision to constrain these model dependent parameters at the +/1 20%-25% (one standard deviation) level. These constraints include only the experimental uncertainties, and further studies are needed to compute the corresponding theoretical uncertainties.

The production of prompt charged particles in proton-lead collisions and in proton-proton collisions at the nucleon-nucleon centre-of-mass energy ${\sqrt{s_{\scriptscriptstyle\mathrm{NN}}}=5\,\mathrm{TeV}}$ is studied at LHCb as a function of pseudorapidity ($\eta$) and transverse momentum ($p_{\mathrm{T}}$) with respect to the proton beam direction. The nuclear modification factor for charged particles is determined as a function of $\eta$ between ${-4.8<\eta<-2.5}$ (backward region) and ${2.0<\eta<4.8}$ (forward region), and $p_{\mathrm{T}}$ between ${0.2<p_{\mathrm{T}}<8.0\,\mathrm{GeV}/c}$. The results show a suppression of charged particle production in proton-lead collisions relative to proton-proton collisions in the forward region and an enhancement in the backward region for $p_{\mathrm{T}}$ larger than $1.5\,\mathrm{GeV}/c$. This measurement constrains nuclear PDFs and saturation models at previously unexplored values of the parton momentum fraction down to $10^{-6}$.

Azimuthal angle (Delta phi) correlations are presented for a broad range of transverse momentum (0.4 < pT < 10 GeV/c) and centrality (0-92%) selections for charged hadrons from di-jets in Au+Au collisions at sqrt(s_NN) = 200 GeV. With increasing pT, the away-side Delta phi distribution evolves from a broad and relatively flat shape to a concave shape, then to a convex shape. Comparisons to p+p data suggest that the away-side distribution can be divided into a partially suppressed head region centered at Delta phi ~ \pi, and an enhanced shoulder region centered at Delta phi ~ \pi \pm 1:1. The pT spectrum for the associated hadrons in the head region softens toward central collisions. The spectral slope for the shoulder region is independent of centrality and trigger pT . The properties of the near-side distributions are also modified relative to those in p + p collisions, reflected by the broadening of the jet shape in Delta phi and Delta eta, and an enhancement of the per-trigger yield. However, these modifications seem to be limited to pT < 4 GeV/c, above which both the dihadron pair shape and per-trigger yield become similar to p + p collisions. These observations suggest that both the away- and near-side distributions contain a jet fragmentation component which dominates for pT \ge 5GeV and a medium-induced component which is important for pT \le 4 GeV/c. We also quantify the role of jets at intermediate and low pT through the yield of jet-induced pairs in comparison to binary scaled p + p pair yield. The yield of jet-induced pairs is suppressed at high pair proxy energy (sum of the pT magnitudes of the two hadrons) and is enhanced at low pair proxy energy. The former is consistent with jet quenching/ the latter is consistent with the enhancement of soft hadron pairs due to transport of lost energy to lower pT.

High-energy proton- and deuteron-nucleus collisions provide an excellent tool for studying a wide array of physics effects, including modifications of parton distribution functions in nuclei, gluon saturation, and color neutralization and hadronization in a nuclear environment, among others. All of these effects are expected to have a significant dependence on the size of the nuclear target and the impact parameter of the collision, also known as the collision centrality. In this article, we detail a method for determining centrality classes in p(d)+A collisions via cuts on the multiplicity at backward rapidity (i.e., the nucleus-going direction) and for determining systematic uncertainties in this procedure. For d+Au collisions at sqrt(s_NN) = 200 GeV we find that the connection to geometry is confirmed by measuring the fraction of events in which a neutron from the deuteron does not interact with the nucleus. As an application, we consider the nuclear modification factors R_{p(d)+A}, for which there is a potential bias in the measured centrality dependent yields due to auto-correlations between the process of interest and the backward rapidity multiplicity. We determine the bias correction factor within this framework. This method is further tested using the HIJING Monte Carlo generator. We find that for d+Au collisions at sqrt(s_NN)=200 GeV, these bias corrections are small and vary by less than 5% (10%) up to p_T = 10 (20) GeV. In contrast, for p+Pb collisions at sqrt(s_NN) = 5.02 TeV we find these bias factors are an order of magnitude larger and strongly p_T dependent, likely due to the larger effect of multi-parton interactions.

We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.

Direct photons have been measured in sqrt(s_NN)=200 GeV d+Au collisions at midrapidity. A wide p_T range is covered by measurements of nearly-real virtual photons (1<p_T<6 GeV/c) and real photons (5<p_T<16 GeV/c). The invariant yield of the direct photons in d+Au collisions over the scaled p+p cross section is consistent with unity. Theoretical calculations assuming standard cold nuclear matter effects describe the data well for the entire p_T range. This indicates that the large enhancement of direct photons observed in Au+Au collisions for 1.0<p_T<2.5 GeV/c is due to a source other than the initial-state nuclear effects.

J/psi production has been measured in proton-proton collisions at sqrt(s)= 200 GeV over a wide rapidity and transverse momentum range by the PHENIX experiment at RHIC. Distributions of the rapidity and transverse momentum, along with measurements of the mean transverse momentum and total production cross section are presented and compared to available theoretical calculations. The total J/psi cross section is 3.99 +/- 0.61(stat) +/- 0.58(sys) +/- 0.40(abs) micro barns. The mean transverse momentum is 1.80 +/- 0.23(stat) +/- 0.16(sys) GeV/c.

We report the measurement of direct photons at midrapidity in Au+Au collisions at sqrt{s_NN} = 200 GeV. The direct photon signal was extracted for the transverse-momentum range of 4 GeV/c < p_T < 22 GeV/c, using a statistical method to subtract decay photons from the inclusive-photon sample. The direct-photon nuclear-modification factor R_AA was calculated as a function of p_T for different Au+Au collision centralities using the measured p+p direct-photon spectrum and compared to theoretical predictions. R_AA was found to be consistent with unity for all centralities over the entire measured p_T range. Theoretical models that account for modifications of initial-direct-photon production due to modified-parton-distribution functions in Au and the different isospin composition of the nuclei, predict a modest change of R_AA from unity and are consistent with the data. Models with compensating effects of the quark-gluon plasma on high-energy photons, such as suppression of jet-fragmentation photons and induced-photon bremsstrahlung from partons traversing the medium, are also consistent with this measurement.

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.

The differential cross section for the production of direct photons in p+p collisions at sqrt(s)=200 GeV at midrapidity was measured in the PHENIX detector at the Relativistic Heavy Ion Collider. Inclusive-direct photons were measured in the transverse-momentum range from 5.5--25 GeV/c, extending the range beyond previous measurements. Event structure was studied with an isolation criterion. Next-to-leading-order perturbative-quantum-chromodynamics calculations give a good description of the spectrum. When the cross section is expressed versus x_T, the PHENIX data are seen to be in agreement with measurements from other experiments at different center-of-mass energies.

The second Fourier component v_2 of the azimuthal anisotropy with respect to the reaction plane was measured for direct photons at midrapidity and transverse momentum (p_T) of 1--13 GeV/c in Au+Au collisions at sqr(s_NN)=200 GeV. Previous measurements of this quantity for hadrons with p_T < 6 GeV/c indicate that the medium behaves like a nearly perfect fluid, while for p_T > 6 GeV/c a reduced anisotropy is interpreted in terms of a path-length dependence for parton energy loss. In this measurement with the PHENIX detector at the Relativistic Heavy Ion Collider we find that for p_T > 4 GeV/c the anisotropy for direct photons is consistent with zero, as expected if the dominant source of direct photons is initial hard scattering. However, in the p_T < 4 GeV/c region dominated by thermal photons, we find a substantial direct photon v_2 comparable to that of hadrons, whereas model calculations for thermal photons in this kinematic region significantly underpredict the observed v_2.

The invariant yields for $J/\psi$ production at forward rapidity $(1.2<|y|<2.2)$ in U$+$U collisions at $\sqrt{s_{_{NN}}}$=193 GeV have been measured as a function of collision centrality. The invariant yields and nuclear-modification factor $R_{AA}$ are presented and compared with those from Au$+$Au collisions in the same rapidity range. Additionally, the direct ratio of the invariant yields from U$+$U and Au$+$Au collisions within the same centrality class is presented, and used to investigate the role of $c\bar{c}$ coalescence. Two different parameterizations of the deformed Woods-Saxon distribution were used in Glauber calculations to determine the values of the number of nucleon-nucleon collisions in each centrality class, $N_{\rm coll}$, and these were found to give significantly different $N_{\rm coll}$ values. Results using $N_{\rm coll}$ values from both deformed Woods-Saxon distributions are presented. The measured ratios show that the $J/\psi$ suppression, relative to binary collision scaling, is similar in U$+$U and Au$+$Au for peripheral and midcentral collisions, but that $J/\psi$ show less suppression for the most central U$+$U collisions. The results are consistent with a picture in which, for central collisions, increase in the $J/\psi$ yield due to $c\bar{c}$ coalescence becomes more important than the decrease in yield due to increased energy density. For midcentral collisions, the conclusions about the balance between $c\bar{c}$ coalescence and suppression depend on which deformed Woods-Saxon distribution is used to determine $N_{\rm coll}$.

The measurement of the direct-photon spectrum from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV is presented by the PHENIX collaboration using the external-photon-conversion technique for 0%--93% central collisions in a transverse-momentum ($p_T$) range of 0.8--10 GeV/$c$. An excess of direct photons, above prompt-photon production from hard-scattering processes, is observed for $p_T<6$ GeV/$c$. Nonprompt direct photons are measured by subtracting the prompt component, which is estimated as $N_{\rm coll}$-scaled direct photons from $p$ $+$ $p$ collisions at 200 GeV, from the direct-photon spectrum. Results are obtained for $0.8<p_T<6.0$ GeV/$c$ and suggest that the spectrum has an increasing inverse slope from ${\approx}0.2$ to 0.4 GeV/$c$ with increasing $p_T$, which indicates a possible sensitivity of the measurement to photons from earlier stages of the evolution of the collision. In addition, like the direct-photon production, the $p_T$-integrated nonprompt direct-photon yields also follow a power-law scaling behavior as a function of collision-system size. The exponent, $\alpha$, for the nonprompt component is found to be consistent with 1.1 with no apparent $p_T$ dependence.

The invariant differential cross section for inclusive neutral pion production in p+p collisions at sqrt(s_NN) = 200 GeV has been measured at mid-rapidity |eta| < 0.35 over the range 1 < p_T <~ 14 GeV/c by the PHENIX experiment at RHIC. Predictions of next-to-leading order perturbative QCD calculations are consistent with these measurements. The precision of our result is sufficient to differentiate between prevailing gluon-to-pion fragmentation functions.

Transverse momentum distributions and yields for $\pi^{\pm}$, $K^{\pm}$, $p$ and $\bar{p}$ in $p+p$ collisions at $\sqrt{s}$=200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter $T_{\rm inv}$, mean transverse momentum $<p_T>$ and yield per unit rapidity $dN/dy$ at each energy, and compare them to other measurements at different $\sqrt{s}$ in $p+p$ and $p+\bar{p}$ collisions. We also present the scaling properties such as $m_T$ scaling, $x_T$ scaling on the $p_T$ spectra between different energies. To discuss the mechanism of the particle production in $p+p$ collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.

New measurements by the PHENIX experiment at RHIC for eta production at midrapidity as a function of transverse momentum (p_T) and collision centrality in sqrt(s_NN) = 200 GeV Au+Au and p+p collisions are presented. They indicate nuclear modification factors (R_AA) that are similar both in magnitude and trend to those found in earlier pi^0 measurements. Linear fits to R_AA in the 5--20 GeV/c p_T region show that the slope is consistent with zero within two standard deviations at all centralities although a slow rise cannot be excluded. Having different statistical and systematic uncertainties the pi^0 and eta measurements are complementary at high p_T/ thus, along with the extended p_T range of these data they can provide additional constraints for theoretical modeling and the extraction of transport properties.

The production cross-sections of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=5$ TeV are measured using a data sample corresponding to an integrated luminosity of $9.13\pm0.18~\text{pb}^{-1}$, collected by the LHCb experiment. The cross-sections are measured differentially as a function of transverse momentum, $p_{\text{T}}$, and rapidity, $y$, and separately for $J/\psi$ mesons produced promptly and from beauty hadron decays (nonprompt). With the assumption of unpolarised $J/\psi$ mesons, the production cross-sections integrated over the kinematic range $0<p_{\text{T}}<20~\text{GeV}/c$ and $2.0<y<4.5$ are $8.154\pm0.010\pm0.283~\mu\text{b}$ for prompt $J/\psi$ mesons and $0.820\pm0.003\pm0.034~\mu\text{b}$ for nonprompt $J/\psi$ mesons, where the first uncertainties are statistical and the second systematic. These cross-sections are compared with those at $\sqrt{s}=8$ TeV and $13$ TeV, and are used to update the measurement of the nuclear modification factor in proton-lead collisions for $J/\psi$ mesons at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\text{NN}}}=5$ TeV. The results are compared with theoretical predictions.

We report on charmonium measurements [J/psi(1S), psi'(2S), and chi_c(1P)] in p+p collisions at sqrt(s)=200 GeV. We find that the fraction of J/psi coming from the feed-down decay of psi' and chi_c in the midrapidity region ($|\eta|<0.35$) is 9.6+/-2.4% and 32+/-9%, respectively. We also report new, higher statistics p_T and rapidity dependencies of the J/psi yield via dielectron decay in the same midrapidity range and at forward rapidity (1.2<|eta|<2.4) via dimuon decay. These results are compared with measurements from other experiments and discussed in the context of current charmonium production models.

The production of deuterons and antideuterons in the transverse momentum range 1.1 < p_T < 4.3 GeV/c at mid-rapidity in Au + Au collisions at sqrt(s_NN)=200 GeV has been studied by the PHENIX experiment at RHIC. A coalescence analysis comparing the deuteron and antideuteron spectra with those of protons and antiprotons, has been performed. The coalescence probability is equal for both deuterons and antideuterons and increases as a function of p_T, which is consistent with an expanding collision zone. Comparing (anti)proton yields p_bar/p = 0.73 +/- 0.01, with (anti)deuteron yields: d_bar/d = 0.47 +/- 0.03, we estimate that n_bar/n = 0.64 +/- 0.04.

Flow coefficients v_n for n = 2, 3, 4, characterizing the anisotropic collective flow in Au+Au collisions at sqrt(s_NN) = 200 GeV, are measured relative to event planes Ψ_n determined at large rapidity. We report v_n as a function of transverse momentum and collision centrality, and study the correlations among the event planes of different order n. The v_n are well described by hydrodynamic models which employ a Glauber Monte Carlo initial state geometry with fluctuations, providing additional constraining power on the interplay between initial conditions and the effects of viscosity as the system evolves. This new constraint improves precision of the extracted viscosity to entropy density ratio eta/s.

Heavy quarkonia are observed to be suppressed in relativistic heavy ion collisions relative to their production in p+p collisions scaled by the number of binary collisions. In order to determine if this suppression is related to color screening of these states in the produced medium, one needs to account for other nuclear modifications including those in cold nuclear matter. In this paper, we present new measurements from the PHENIX 2007 data set of J/psi yields at forward rapidity (1.2<|y|<2.2) in Au+Au collisions at sqrt(s_NN)=200 GeV. The data confirm the earlier finding that the suppression of J/psi at forward rapidity is stronger than at midrapidity, while also extending the measurement to finer bins in collision centrality and higher transverse momentum (pT). We compare the experimental data to the most recent theoretical calculations that incorporate a variety of physics mechanisms including gluon saturation, gluon shadowing, initial-state parton energy loss, cold nuclear matter breakup, color screening, and charm recombination. We find J/psi suppression beyond cold-nuclear-matter effects. However, the current level of disagreement between models and d+Au data precludes using these models to quantify the hot-nuclear-matter suppression.

The PHENIX collaboration presents first measurements of low-momentum ($0.4<p_T<3$ GeV/$c$) direct-photon yields from Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=39 and 62.4 GeV. For both beam energies the direct-photon yields are substantially enhanced with respect to expectations from prompt processes, similar to the yields observed in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200. Analyzing the photon yield as a function of the experimental observable $dN_{\rm ch}/d\eta$ reveals that the low-momentum ($>$1\,GeV/$c$) direct-photon yield $dN_{\gamma}^{\rm dir}/d\eta$ is a smooth function of $dN_{\rm ch}/d\eta$ and can be well described as proportional to $(dN_{\rm ch}/d\eta)^\alpha$ with $\alpha{\sim}$1.25. This new scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and Large Hadron Collider, for centrality selected samples, as well as for different, $A$$+$$A$ collision systems. At a given beam energy the scaling also holds for high $p_T$ ($>5$\,GeV/$c$) but when results from different collision energies are compared, an additional $\sqrt{s_{_{NN}}}$-dependent multiplicative factor is needed to describe the integrated-direct-photon yield.

Cross sections for mid-rapidity production of direct photons in p+p collisions at the Relativistic Heavy Ion Collider (RHIC) are reported for 3 < p_T < 16 GeV/c. Next-to-leading order (NLO) perturbative QCD (pQCD) describes the data well for p_T > 5 GeV/c, where the uncertainties of the measurement and theory are comparable. We also report on the effect of requiring the photons to be isolated from parton jet energy. The observed fraction of isolated photons is well described by pQCD for p_T > 7 GeV/c.

We have measured direct photons for $p_T<5~$GeV/$c$ in minimum bias and 0\%--40\% most central events at midrapidity for Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The $e^{+}e^{-}$ contribution from quasi-real direct virtual photons has been determined as an excess over the known hadronic contributions in the $e^{+}e^{-}$ mass distribution. A clear enhancement of photons over the binary scaled $p$$+$$p$ fit is observed for $p_T<4$ GeV/$c$ in Cu$+$Cu data. The $p_T$ spectra are consistent with the Au$+$Au data covering a similar number of participants. The inverse slopes of the exponential fits to the excess after subtraction of the $p$$+$$p$ baseline are 285$\pm$53(stat)$\pm$57(syst)~MeV/$c$ and 333$\pm$72(stat)$\pm$45(syst)~MeV/$c$ for minimum bias and 0\%--40\% most central events, respectively. The rapidity density, $dN/dy$, of photons demonstrates the same power law as a function of $dN_{\rm ch}/d\eta$ observed in Au$+$Au at the same collision energy.

We present inclusive charged hadron elliptic flow v_2 measured over the pseudorapidity range |\eta| < 0.35 in Au+Au collisions at sqrt(s_NN) = 200 GeV. Results for v_2 are presented over a broad range of transverse momentum (p_T = 0.2-8.0 GeV/c) and centrality (0-60%). In order to study non-flow effects that are not correlated with the reaction plane, as well as the fluctuations of v_2, we compare two different analysis methods: (1) event plane method from two independent sub-detectors at forward (|\eta| = 3.1-3.9) and beam (|\eta| > 6.5) pseudorapidities and (2) two-particle cumulant method extracted using correlations between particles detected at midrapidity. The two event-plane results are consistent within systematic uncertainties over the measured p_T and in centrality 0-40%. There is at most 20% difference of the v_2 between the two event plane methods in peripheral (40-60%) collisions. The comparisons between the two-particle cumulant results and the standard event plane measurements are discussed.

A measurement of direct photons in p+p collisions at sqrt(s)=200 GeV is presented. A photon excess above background from pi^0 --> gamma+gamma, eta --> gamma+gamma, and other decays is observed in the transverse momentum range 5.5 < p_T < 7 GeV/c. The result is compared to a next-to-leading-order perturbative QCD calculation. Within errors, good agreement is found between the QCD calculation and the measured result.

The properties of jets produced in p+p collisions at sqrt(s)=200 GeV are measured using the method of two particle correlations. The trigger particle is a leading particle from a large transverse momentum jet while the associated particle comes from either the same jet or the away-side jet. Analysis of the angular width of the near-side peak in the correlation function determines the jet fragmentation transverse momentum j_T . The extracted value, sqrt(<j_T^2>)= 585 +/- 6(stat) +/- 15(sys) MeV/c, is constant with respect to the trigger particle transverse momentum, and comparable to the previous lower sqrt(s) measurements. The width of the away-side peak is shown to be a convolution of j_T with the fragmentation variable, z, and the partonic transverse momentum, k_T . The <z> is determined through a combined analysis of the measured pi^0 inclusive and associated spectra using jet fragmentation functions measured in e^+e^-. collisions. The final extracted values of k_T are then determined to also be independent of the trigger particle transverse momentum, over the range measured, with value of sqrt(<k_T^2>) = 2.68 +/- 0.07(stat) +/- 0.34(sys) GeV/c.

Fast parton probes produced by hard scattering and embedded within collisions of large nuclei have shown that partons suffer large energy loss and that the produced medium may respond collectively to the lost energy. We present measurements of neutral pion trigger particles at transverse momenta p^t_T = 4-12 GeV/c and associated charged hadrons (p^a_T = 0.5-7 GeV/c) as a function of relative azimuthal angle Delta Phi at midrapidity in Au+Au and p+p collisions at sqrt(s_NN) = 200 GeV. These data lead to two major observations. First, the relative angular distribution of low momentum hadrons, whose shape modification has been interpreted as a medium response to parton energy loss, is found to be modified only for p^t_T &lt; 7 GeV/c. At higher p^t_T, the data are consistent with unmodified or very weakly modified shapes, even for the lowest measured p^a_T. This observation presents a quantitative challenge to medium response scenarios. Second, the associated yield of hadrons opposite to the trigger particle in Au+Au relative to that in p+p (I_AA) is found to be suppressed at large momentum (IAA ~ 0.35-0.5), but less than the single particle nuclear modification factor (R_AA ~0.2).

The PHENIX collaboration at the Relativistic Heavy Ion Collider (RHIC) reports measurements of azimuthal dihadron correlations near midrapidity in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. These measurements complement recent analyses by experiments at the Large Hadron Collider (LHC) involving central $p$$+$Pb collisions at $\sqrt{s_{_{NN}}}$=5.02 TeV, which have indicated strong anisotropic long-range correlations in angular distributions of hadron pairs. The origin of these anisotropies is currently unknown. Various competing explanations include parton saturation and hydrodynamic flow. We observe qualitatively similar, but larger, anisotropies in $d$$+$Au collisions compared to those seen in $p$$+$Pb collisions at the LHC. The larger extracted $v_2$ values in $d$$+$Au collisions at RHIC are consistent with expectations from hydrodynamic calculations owing to the larger expected initial-state eccentricity compared with that from $p$$+$Pb collisions. When both are divided by an estimate of the initial-state eccentricity the scaled anisotropies follow a common trend with multiplicity that may extend to heavy ion data at RHIC and the LHC, where the anisotropies are widely thought to arise from hydrodynamic flow.

We report the observation at the Relativistic Heavy Ion Collider (RHIC) of suppression of back-to-back correlations in the direct photon+jet channel in Au+Au relative to p+p collisions. Two-particle correlations of direct photon triggers with associated hadrons are obtained by statistical subtraction of the decay photon-hadron background. The initial momentum of the away-side parton is tightly constrained, because the parton-photon pair exactly balance in momentum at leading order in perturbative quantum chromodynamics (pQCD), making such correlations a powerful probe of the in-medium parton energy loss. The away-side nuclear suppression factor, I_AA, in central Au+Au collisions, is 0.32 +/- 0.12(stat) +/- 0.09(syst) for hadrons of 3 < p_T < 5 in coincidence with photons of 5 < p_T < 15 GeV/c. The suppression is comparable to that observed for high-p_T single hadrons and dihadrons. The direct photon associated yields in p+p collisions scale approximately with the momentum balance, z_T = p_T^hadron/p_T^photon, as expected for a measure of the away-side parton fragmentation function. We compare to Au+Au collisions for which the momentum balance dependence of the nuclear modification should be sensitive to the path-length dependence of parton energy loss.

The interpretation of cosmic antiproton flux measurements from space-borne experiments is currently limited by the knowledge of the antiproton production cross-section in collisions between primary cosmic rays and the interstellar medium. Using collisions of protons with an energy of 6.5 TeV incident on helium nuclei at rest in the proximity of the interaction region of the LHCb experiment, the ratio of antiprotons originating from antihyperon decays to prompt production is measured for antiproton momenta between 12 and 110 GeV/c. The dominant antihyperon contribution, namely $\bar{\Lambda} \to \bar{p} \pi^+$ decays from promptly produced $\bar{\Lambda}$ particles, is also exclusively measured. The results complement the measurement of prompt antiproton production obtained from the same data sample. At the energy scale of this measurement, the antihyperon contributions to antiproton production are observed to be significantly larger than predictions of commonly used hadronic production models.

Neutral-pion, pi^0, spectra were measured at midrapidity (|y|<0.35) in Au+Au collisions at sqrt(s_NN) = 39 and 62.4 GeV and compared to earlier measurements at 200 GeV in the 1<p_T<10 GeV/c transverse-momentum (p_T) range. The high-p_T tail is well described by a power law in all cases and the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very similar to that observed in the corresponding p+p-collision spectra. The nuclear-modification factors (R_AA) show significant suppression and a distinct energy dependence at moderate p_T in central collisions. At high p_T, R_AA is similar for 62.4 and 200 GeV at all centralities. Perturbative-quantum-chromodynamics calculations that describe R_AA well at 200 GeV, fail to describe the 39 GeV data, raising the possibility that the relative importance of initial-state effects and soft processes increases at lower energies. A conclusion that the region where hard processes are dominant is reached only at higher p_T, is also supported by the x_T dependence of the x_T-scaling power-law exponent.

The first measurement of direct photons in Au+Au collisions at sqrt(s_NN) = 200 GeV is presented. The direct photon signal is extracted as a function of the Au+Au collision centrality and compared to NLO pQCD calculations. The direct photon yield is shown to scale with the number of nucleon-nucleon collisions for all centralities.

A precision measurement of the $Z$ boson production cross-section at $\sqrt{s} = 13$ TeV in the forward region is presented, using $pp$ collision data collected by the LHCb detector, corresponding to an integrated luminosity of 5.1 fb$^{-1}$. The production cross-section is measured using $Z\rightarrow\mu^+\mu^-$ events within the fiducial region defined as pseudorapidity $2.0<\eta<4.5$ and transverse momentum $p_{T}>20$ GeV/$c$ for both muons and dimuon invariant mass $60<M_{\mu\mu}<120$ GeV/$c^2$. The integrated cross-section is determined to be $\sigma (Z \rightarrow \mu^+ \mu^-)$ = 196.4 $\pm$ 0.2 $\pm$ 1.6 $\pm$ 3.9~pb, where the first uncertainty is statistical, the second is systematic, and the third is due to the luminosity determination. The measured results are in agreement with theoretical predictions within uncertainties.

The invariant differential cross section for inclusive electron production in $p + p$ collisions at $\sqrt{s} = 200$~GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range $0.4 \le p_T \le 5.0$~GeV/$c$ in the central rapidity region ($|\eta| \le 0.35$). The contribution to the inclusive electron spectrum from semileptonic decays of hadrons carrying heavy flavor, {\it i.e.} charm quarks or, at high $p_T$, bottom quarks, is determined via three independent methods. The resulting electron spectrum from heavy flavor decays is compared to recent leading and next-to-leading order perturbative QCD calculations. The total cross section of charm quark-antiquark pair production is determined to be $\sigma_{c\bar{c}} = 0.92 \pm 0.15 {\rm (stat.)} \pm 0.54 {\rm (sys.)}$~mb.

Bose-Einstein correlations of charged kaons are measured for Au+Au collisions at sqrt(s_NN) = 200 GeV and are compared to charged pion probes, which have a larger hadronic scattering cross section. Three dimensional Gaussian source radii are extracted, along with a one-dimensional kaon emission source function. The centrality dependences of the three Gaussian radii are well described by a single linear function if N_part^1/3 with zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r >~ 10 fm, although the bulk emission at lower radius is well-described by a Gaussian. The presence of a non-Gaussian tail in the kaon source reaffirms that the particle emission region in a heavy ion collision is extended, and that similar measurements with pions are not solely due to the decay of long-lived resonances.

The anisotropy parameter v_2, the second harmonic of the azimuthal particles distribution, has been measured with the PHENIX detector in Au+Au collisions at sqrt(s_NN) = 200 GeV for identified and inclusive charged particles at central rapidities (|eta| < 0.35) with respect to the reaction plane defined at high rapidities (|eta| = 3-4). The v_2 for all particles reaches a maximum at mid-centrality, and increases with p_T up to 2 GeV/c and then saturates or decreases slightly. Our results depart from hydrodynamically predicted behavior above 2 GeV/c. A quark coalescence model is also investigated.

We present measurements of electrons and positrons from the semileptonic decays of heavy-flavor hadrons at midrapidity ($|y|<$ 0.35) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. The data were collected in 2010 by the PHENIX experiment that included the new hadron-blind detector. The invariant yield of electrons from heavy-flavor decays is measured as a function of transverse momentum in the range $1<p_T^e<5$ GeV/$c$. The invariant yield per binary collision is slightly enhanced above the $p$$+$$p$ reference in Au$+$Au 0%--20%, 20%--40% and 40%--60% centralities at a comparable level. This may be a result of the interplay between initial-state Cronin effects, final-state flow, and energy loss for heavy-quark production at this low beam energy. The $v_2$ of electrons from heavy-flavor decays is nonzero when averaged between $1.3<p_T^e<2.5$ GeV/$c$ from $0<{\rm centrality}<40$% collisions at $\sqrt{s_{_{NN}}}=62.4$ GeV. For 20%--40% centrality collisions, the $v_2$ at $\sqrt{s_{_{NN}}}=62.4$ GeV is smaller than that for heavy flavor decays at $\sqrt{s_{_{NN}}}=200$ GeV. The $v_2$ of the electrons from heavy-flavor decay at the lower beam energy is also smaller than $v_2$ for pions. Both results indicate that the heavy-quarks interact with the medium formed in these collisions, but they may not be at the same level of thermalization with the medium as observed at $\sqrt{s_{_{NN}}}=200$ GeV.

Differential measurements of the elliptic (v_2) and hexadecapole (v_4) Fourier flow coefficients are reported for charged hadrons as a function of transverse momentum (p_T) and collision centrality or the number of participant nucleons (N_part) for Au+Au collisions at sqrt(s_NN)=200 GeV. The v_{2,4} measurements at pseudorapidity |\eta|<=0.35 obtained with four separate reaction plane detectors positioned in the range 1.0<|\eta|<3.9 show good agreement, indicating the absence of significant \eta-dependent nonflow perturbations. Sizable values for v_4(p_T) are observed with a ratio v_4(p_T,N_part)/v_2^2(p_T,N_part)~0.8 for 50<N_part<200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N_part>200 this ratio increases up to 1.7 in the most central collisions.

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.

Two particle correlations between identified meson and baryon trigger particles with 2.5 < p_T < 4.0 GeV/c and lower p_T charged hadrons have been measured at midrapidity by the PHENIX experiment at RHIC in p+p, d+Au and Au+Au collisions at sqrt(s_NN) = 200 GeV. The probability of finding a hadron near in azimuthal angle to the trigger particle is almost identical for leading mesons and baryons for non-central Au+Au. The yield for both trigger baryons and mesons is significantly higher in Au+Au than in p+p and d+Au, except for trigger baryons in central collisions. The baryon excess is likely to arise predominantly from hard scattering processes.

Dihadron correlations at high transverse momentum in d+Au collisions at sqrt(s_NN) = 200 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). From these correlations we extract several structural characteristics of jets; the root-mean-squared (RMS) transverse momentum of fragmenting hadrons with respect to the jet sqrt(<j_T^2>), the mean sine-squared angle between the scattered partons <sin^2(phi_jj)>, and the number of particles produced within the dijet that are associated with a high-p_T particle (dN/dx_E distributions). We observe that the fragmentation characteristics of jets in d+Au collisions are very similar to those in p+p collisions and that there is also little dependence on the centrality of the d+Au collision. This is consistent with the nuclear medium having little influence on the fragmentation process. Furthermore, there is no statistically significant increase in the value of <sin^2(phi_jj)> from p+p to d+Au collisions. This constrains the amount of multiple scattering that partons undergo in the cold nuclear medium before and after a hard-collision.