We employ data taken by the JADE and OPAL experiments for an integrated QCD study in hadronic e+e- annihilations at c.m.s. energies ranging from 35 GeV through 189 GeV. The study is based on jet-multiplicity related observables. The observables are obtained to high jet resolution scales with the JADE, Durham, Cambridge and cone jet finders, and compared with the predictions of various QCD and Monte Carlo models. The strong coupling strength, alpha_s, is determined at each energy by fits of O(alpha_s^2) calculations, as well as matched O(alpha_s^2) and NLLA predictions, to the data. Matching schemes are compared, and the dependence of the results on the choice of the renormalization scale is investigated. The combination of the results using matched predictions gives alpha_s(MZ)=0.1187+{0.0034}-{0.0019}. The strong coupling is also obtained, at lower precision, from O(alpha_s^2) fits of the c.m.s. energy evolution of some of the observables. A qualitative comparison is made between the data and a recent MLLA prediction for mean jet multiplicities.

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.

Using data collected with the CLEO II detector at the Cornell Electron Storage Ring, we determine the ratio R(chrg) for the mean charged multiplicity observed in Upsilon(1S)->gggamma events, to the mean charged multiplicity observed in e+e- -> qqbar gamma events. We find R(chrg)=1.04+/-0.02+/-0.05 for jet-jet masses less than 7 GeV.

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.

Measurements of the midrapidity transverse energy distribution, $d\Et/d\eta$, are presented for $p$$+$$p$, $d$$+$Au, and Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV and additionally for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ and 130 GeV. The $d\Et/d\eta$ distributions are first compared with the number of nucleon participants $N_{\rm part}$, number of binary collisions $N_{\rm coll}$, and number of constituent-quark participants $N_{qp}$ calculated from a Glauber model based on the nuclear geometry. For Au$+$Au, $\mean{d\Et/d\eta}/N_{\rm part}$ increases with $N_{\rm part}$, while $\mean{d\Et/d\eta}/N_{qp}$ is approximately constant for all three energies. This indicates that the two component ansatz, $dE_{T}/d\eta \propto (1-x) N_{\rm part}/2 + x N_{\rm coll}$, which has been used to represent $E_T$ distributions, is simply a proxy for $N_{qp}$, and that the $N_{\rm coll}$ term does not represent a hard-scattering component in $E_T$ distributions. The $dE_{T}/d\eta$ distributions of Au$+$Au and $d$$+$Au are then calculated from the measured $p$$+$$p$ $E_T$ distribution using two models that both reproduce the Au$+$Au data. However, while the number-of-constituent-quark-participant model agrees well with the $d$$+$Au data, the additive-quark model does not.

All of the experimental data points presented in the original paper are correct and unchanged (including statistical and systematic uncertainties). However, herein we correct a comparison between the experimental data and a theoretical picture, because we discovered a mistake in the code used. All of the most probable sigma_breakup values differ by less than 0.4 mb from those originally presented. However, the one standard deviation uncertainties (that include contributions from both the statistical and systematic uncertainties on the experimental data points) are approximately 30-60% larger than originally reported. We give a table of the new comparison results and corrected versions of Figs. 8-11 of the original paper and we note that no correction is needed for results from the data-driven method in Fig. 13.

Charmonium is a valuable probe in heavy-ion collisions to study the properties of the quark gluon plasma, and is also an interesting probe in small collision systems to study cold nuclear matter effects, which are also present in large collision systems. With the recent observations of collective behavior of produced particles in small system collisions, measurements of the modification of charmonium in small systems have become increasingly relevant. We present the results of J/ψ measurements at forward and backward rapidity in various small collision systems, p+p, p+Al, p+Au and 3He+Au, at √sNN =200 GeV. The results are presented in the form of the observable RAB, the nuclear modification factor, a measure of the ratio of the J/ψ invariant yield compared to the scaled yield in p+p collisions. We examine the rapidity, transverse momentum, and collision centrality dependence of nuclear effects on J/ψ production with different projectile sizes p and 3He, and different target sizes Al and Au. The modification is found to be strongly dependent on the target size, but to be very similar for p+Au and 3He+Au. However, for 0%–20% central collisions at backward rapidity, the modification for 3He+Au is found to be smaller than that for p+Au, with a mean fit to the ratio of 0.89±0.03(stat)±0.08(syst), possibly indicating final state effects due to the larger projectile size.

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.

We present measurements of azimuthal correlations of charged hadron pairs in $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions after subtracting an underlying event using a model that includes higher-order azimuthal anisotropy $v_2$, $v_3$, and $v_4$. After subtraction, the away-side ($\Delta\phi\sim\pi)$ of the highest transverse-momentum trigger ($p_T>4$ GeV/$c$) correlations is suppressed compared to that of correlations measured in $p$$+$$p$ collisions. At the lowest associated particle $p_T$, the away-side shape and yield are modified. These observations are consistent with the scenario of radiative-jet energy loss. For the lowest-$p_T$ trigger correlations, an away-side yield exists and we explore the dependence of the shape of the away-side within the context of an underlying-event model. Correlations are also studied differentially versus event-plane angle $\Psi_n$. The angular correlations show an asymmetry when selecting the sign of the trigger-particle azimuthal angle with respect to the $\Psi_2$ event plane. This asymmetry and the measured suppression of the pair yield out of plane is consistent with a path-length-dependent energy loss. No $\Psi_3$ dependence can be resolved within experimental uncertainties.

Measurements of midrapidity charged particle multiplicity distributions, $dN_{\rm ch}/d\eta$, and midrapidity transverse-energy distributions, $dE_T/d\eta$, are presented for a variety of collision systems and energies. Included are distributions for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV, Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ and 62.4 GeV, Cu$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, U$+$U collisions at $\sqrt{s_{_{NN}}}=193$ GeV, $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Centrality-dependent distributions at midrapidity are presented in terms of the number of nucleon participants, $N_{\rm part}$, and the number of constituent quark participants, $N_{q{\rm p}}$. For all $A$$+$$A$ collisions down to $\sqrt{s_{_{NN}}}=7.7$ GeV, it is observed that the midrapidity data are better described by scaling with $N_{q{\rm p}}$ than scaling with $N_{\rm part}$. Also presented are estimates of the Bjorken energy density, $\varepsilon_{\rm BJ}$, and the ratio of $dE_T/d\eta$ to $dN_{\rm ch}/d\eta$, the latter of which is seen to be constant as a function of centrality for all systems.

Bose-Einstein correlations of identically charged pion pairs were measured by the PHENIX experiment at mid-rapidity in Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. The Bertsch-Pratt radius parameters were determined as a function of the transverse momentum of the pair and as a function of the centrality of the collision. Using the \it{full} Coulomb correction, the ratio $R_{\rm out}/R_{\rm side}$ is smaller than unity for $<k_{\rm T}>$ from 0.25 to 1.2 GeV/c and for all measured centralities. However, using recently developed partial Coulomb correction methods, we find that $R_{\rm out}/R_{\rm side}$ is 0.8-1.1 for the measured $<k_{\rm T}>$ range, and approximately constant at unity with the number of participants.

The production of low mass e+e- pairs for m_{e+e-} < 300 MeV/c^2 and 1 < p_T <5 GeV/c is measured in p+p and Au+Au collisions at sqrt(s_NN)=200 GeV. Enhanced yield above hadronic sources is observed. Treating the excess as internal conversions, the invariant yield of direct photons is deduced. In central Au+Au collisions, the excess of direct photon yield over p+p is exponential in transverse momentum, with inverse slope T = 221 +/- 19 (stat) +/- 19 (syst) MeV. Hydrodynamical models with initial temperatures ranging from 300--600 MeV at times of ~ 0.6 - 0.15 fm/c after the collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to quark gluon plasma at ~ 170 MeV.

There is strong evidence for the formation of small droplets of quark-gluon plasma in $p/d/^{3}$He+Au collisions at the Relativistic Heavy Ion Collider (RHIC) and in $p$+$p$/Pb collisions at the Large Hadron Collider. In particular, the analysis of data at RHIC for different geometries obtained by varying the projectile size and shape has proven insightful. In the present analysis, we find excellent agreement with the previously published PHENIX at RHIC results on elliptical and triangular flow with an independent analysis via the two-particle correlation method, which has quite different systematic uncertainties and an independent code base. In addition, the results are extended to other detector combinations with different kinematic (pseudorapidity) coverage. These results provide additional constraints on contributions from nonflow and longitudinal decorrelations.

Two-pion correlations in sqrt(s_NN)=130 GeV Au+Au collisions at RHIC have been measured over a broad range of pair transverse momentum k_T by the PHENIX experiment at RHIC. The k_T dependent transverse radii are similar to results from heavy ion collisions at sqrt(s_NN) = 4.1, 4.9, and 17.3 GeV, whereas the longitudinal radius increases monotonically with beam energy. The ratio of the outwards to sidewards transverse radii (R_out/R_side) is consistent with unity and independent of k_T.

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.

PHENIX has measured the centrality dependence of charged hadron p_T spectra from central Au+Au collisions at sqrt(s_NN)=130 GeV. The truncated mean p_T decreases with centrality for p_T > 2 GeV/c, indicating an apparent reduction of the contribution from hard scattering to high p_T hadron production. For central collisions the yield at high p_T is shown to be suppressed compared to binary nucleon-nucleon collision scaling of p+p data. This suppression is monotonically increasing with centrality, but most of the change occurs below 30% centrality, i.e. for collisions with less than about 140 participating nucleons. The observed p_T and centrality dependence is consistent with the particle production predicted by models including hard scattering and subsequent energy loss of the scattered partons in the dense matter created in the collisions.

Emission source functions are extracted from correlation functions constructed from charged pions produced at mid-rapidity in Au+Au collisions at sqrt(s_NN)=200 GeV. The source parameters extracted from these functions at low k_T, give first indications of a long tail for the pion emission source. The source extension cannot be explained solely by simple kinematic considerations. The possible role of a halo of secondary pions from resonance emissions is explored.

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.

New measurements are presented for charged hadron azimuthal correlations at mid-rapidity in Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV. They are compared to earlier measurements obtained at sqrt(s_NN) = 130 GeV and in Pb+Pb collisions at sqrt(s_NN) = 17.2 GeV. Sizeable anisotropies are observed with centrality and transverse momentum (p_T) dependence characteristic of elliptic flow (v_2). For a broad range of centralities, the observed magnitudes and trends of the differential anisotropy, v_2(p_T), change very little over the collision energy range sqrt(s_NN) = 62-200 GeV, indicating saturation of the excitation function for v_2 at these energies. Such a saturation may be indicative of the dominance of a very soft equation of state for sqrt(s_NN) = 62-200 GeV.

Differential elliptic flow (v_2) for phi mesons and (anti)deuterons (d^bar)d is measured for Au+Au collisions at sqrt(s_NN) = 200 GeV. The v_2 for phi mesons follows the trend of lighter pi^+/- and K^+/- mesons, suggesting that ordinary hadrons interacting with standard hadronic cross sections are not the primary driver for elliptic flow development. The v_2 values for (d^bar)d suggest that elliptic flow is additive for composite particles. This further validation of the universal scaling of v_2 per constituent quark for baryons and mesons suggests that partonic collectivity dominates the transverse expansion dynamics.

The PHENIX experiment has measured open heavy-flavor production via semileptonic decay muons over the transverse momentum range 1 < pT < 6 GeV/c at forward and backward rapidity (1.4 < |y| < 2.0) in d+Au and p+p collisions at ?sNN = 200 GeV. In central d+Au collisions an enhancement (suppression) of heavy-flavor muon production is observed at backward (forward) rapidity relative to the yield in p+p collisions scaled by the number of binary collisions. Modification of the gluon density distribution in the Au nucleus contributes in terms of anti-shadowing enhancement and shadowing suppression; however, the enhancement seen at backward rapidity exceeds expectations from this effect alone. These results, implying an important role for additional cold nuclear matter effects, serves as a key baseline for heavy-quark measurements in A+A collisions and in constraining the magnitude of charmonia breakup effects at the Relativistic Heavy Ion Collider and the Large Hadron Collider.

We report on charged hadron production in deuteron-gold reactions at sqrt(s_NN) = 200 GeV. Our measurements in the deuteron-direction cover 1.4 < eta < 2.2, referred to as forward rapidity, and in the gold-direction -2.0 < eta < -1.4, referred to as backward rapidity, and a transverse momentum range p_T = 0.5-4.0 GeV/c. We compare the relative yields for different deuteron-gold collision centrality classes. We observe a suppression relative to binary collision scaling at forward rapidity, sensitive to low momentum fraction (x) partons in the gold nucleus, and an enhancement at backward rapidity, sensitive to high momentum fraction partons in the gold nucleus.

Production of $\pi^0$ and $\eta$ mesons has been measured at midrapidity in Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. Measurements were performed in $\pi^0(\eta)\rightarrow\gamma\gamma$ decay channel in the 1(2)-20 GeV/$c$ transverse momentum range. A strong suppression is observed for $\pi^0$ and $\eta$ meson production at high transverse momentum in central Cu$+$Au collisions relative to the $p$$+$$p$ results scaled by the number of nucleon-nucleon collisions. In central collisions the suppression is similar to Au$+$Au with comparable nuclear overlap. The $\eta/\pi^0$ ratio measured as a function of transverse momentum is consistent with $m_T$-scaling parameterization down to $p_T=$2 GeV/$c$, its asymptotic value is constant and consistent with Au$+$Au and $p$$+$$p$ and does not show any significant dependence on collision centrality. Similar results were obtained in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions as well as in $e^+e^-$ collisions in a range of collision energies $\sqrt{s_{_{NN}}}=$3--1800 GeV. This suggests that the quark-gluon-plasma medium produced in Cu$+$Cu collisions either does not affect the jet fragmentation into light mesons or it affects the $\pi^0$ and $\eta$ the same way.

We present results on the measurement of lambda and lambda^bar production in Au+Au collisions at sqrt(s_NN)=130 GeV with the PHENIX detector at RHIC. The transverse momentum spectra were measured for minimum bias and for the 5% most central events. The lambda^bar/lambda ratios are constant as a function of p_T and the number of participants. The measured net lambda density is significantly larger than predicted by models based on hadronic strings (e.g. HIJING) but in approximate agreement with models which include the gluon junction mechanism.

The PHENIX experiment at the Relativistic Heavy Ion Collider measured $\pi^0$ and $\eta$ mesons at midrapidity in U$+$U collisions at $\sqrt{s_{_{NN}}}=192$ GeV in a wide transverse momentum range. Measurements were performed in the $\pi^0(\eta)\rightarrow\gamma\gamma$ decay modes. A strong suppression of $\pi^0$ and $\eta$ meson production at high transverse momentum was observed in central U$+$U collisions relative to binary scaled $p$$+$$p$ results. Yields of $\pi^0$ and $\eta$ mesons measured in U$+$U collisions show similar suppression pattern to the ones measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for similar numbers of participant nucleons. The $\eta$/$\pi^0$ ratios do not show dependence on centrality or transverse momentum, and are consistent with previously measured values in hadron-hadron, hadron-nucleus, nucleus-nucleus, and $e^+e^-$ collisions.

Asymmetric nuclear collisions of $p$$+$Al, $p$$+$Au, $d$$+$Au, and $^{3}$He$+$Au at $\sqrt{s_{_{NN}}}=200$ GeV provide an excellent laboratory for understanding particle production, as well as exploring interactions among these particles after their initial creation in the collision. We present measurements of charged hadron production $dN_{\rm ch}/d\eta$ in all such collision systems over a broad pseudorapidity range and as a function of collision multiplicity. A simple wounded quark model is remarkably successful at describing the full data set. We also measure the elliptic flow $v_{2}$ over a similarly broad pseudorapidity range. These measurements provide key constraints on models of particle emission and their translation into flow.

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.

Neutral pion transverse momentum (pT) spectra at mid-rapidity (|y| < 0.35) were measured in Cu+Cu collisions at \sqrt s_NN = 22.4, 62.4, and 200 GeV. Relative to pi -zero yields in p+p collisions scaled by the number of inelastic nucleon-nucleon collisions (Ncoll) at the respective energies, the pi-zero yields for pT \ge 2 GeV/c in central Cu+Cu collisions at 62.4 and 200 GeV are suppressed, whereas an enhancement is observed at 22.4 GeV. A comparison with a jet quenching model suggests that final state parton energy loss dominates in central Cu+Cu collisions at 62.4 GeV and 200 GeV, while the enhancement at 22.4 GeV is consistent with nuclear modifications in the initial state alone.

Correlations between p and pbar's at transverse momenta typical of enhanced baryon production in Au+Au collisions are reported. The PHENIX experiment measures same and opposite sign baryon pairs in Au+Au collisions at sqrt(s_NN) = 200 GeV. Correlated production of p and p^bar with the trigger particle from the range 2.5 < p_T < 4.0 GeV/c and the associated particle with 1.8 < p_T < 2.5 GeV/c is observed to be nearly independent of the centrality of the collisions. Same sign pairs show no correlation at any centrality. The conditional yield of mesons triggered by baryons (and anti-baryons) and mesons in the same pT range rises with increasing centrality, except for the most central collisions, where baryons show a significantly smaller number of associated mesons. These data are consistent with a picture in which hard scattered partons produce correlated p and p^bar in the p_T region of the baryon excess.

The standard model (SM) of particle physics is spectacularly successful, yet the measured value of the muon anomalous magnetic moment $(g-2)_\mu$ deviates from SM calculations by 3.6$\sigma$. Several theoretical models attribute this to the existence of a "dark photon," an additional U(1) gauge boson, which is weakly coupled to ordinary photons. The PHENIX experiment at the Relativistic Heavy Ion Collider has searched for a dark photon, $U$, in $\pi^0,\eta \rightarrow \gamma e^+e^-$ decays and obtained upper limits of $\mathcal{O}(2\times10^{-6})$ on $U$-$\gamma$ mixing at 90% CL for the mass range $30<m_U<90$ MeV/$c^2$. Combined with other experimental limits, the remaining region in the $U$-$\gamma$ mixing parameter space that can explain the $(g-2)_\mu$ deviation from its SM value is nearly completely excluded at the 90% confidence level, with only a small region of $29<m_U<32$ MeV/$c^2$ remaining.

Detailed differential measurements of the elliptic flow for particles produced in Au+Au and Cu+Cu collisions at sqrt(s_NN) = 200 GeV are presented. Predictions from perfect fluid hydrodynamics for the scaling of the elliptic flow coefficient v_2 with eccentricity, system size and transverse energy are tested and validated. For transverse kinetic energies KE_T ~ m_T-m up to ~1 GeV, scaling compatible with the hydrodynamic expansion of a thermalized fluid is observed for all produced particles. For large values of KE_T, the mesons and baryons scale separately. A universal scaling for the flow of both mesons and baryons is observed for the full transverse kinetic energy range of the data when quark number scaling is employed. In both cases the scaling is more pronounced in terms of KE_T rather than transverse momentum.

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.

Azimuthal angle \Delta\phi correlations are presented for charged hadrons from dijets for 0.4 < p_T < 10 GeV/c in Au+Au collisions at sqrt(s_NN) = 200 GeV. With increasing p_T, the away-side distribution evolves from a broad to a concave shape, then to a convex shape. Comparisons to p+p data suggest that the away-side can be divided into a partially suppressed 'head' region centered at Delta\phi ~ \pi, and an enhanced 'shoulder' region centered at Delta\phi ~ \pi +/- 1.1. The p_T spectrum for the 'head' region softens toward central collisions, consistent with the onset of jet quenching. The spectral slope for the 'shoulder' region is independent of centrality and trigger p_T, which offers constraints on energy transport mechanisms and suggests that the 'shoulder' region contains the medium response to energetic jets.

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.

The PHENIX experiment at RHIC has measured the invariant cross section for omega-meson production at mid-rapidity in the transverse momentum range 2.5 < p_T < 9.25 GeV/c in p+p and d+Au collisions at sqrt(s_NN) = 200 GeV. Measurements in two decay channels (omega --> pi^0 pi^+ pi^- and omega --> pi^0 gamma) yield consistent results, and the reconstructed omega mass agrees with the accepted value within the p_T range of the measurements. The omega/pi^0 ratio is found to be 0.85 +/- 0.05(stat) +/- 0.09(sys) and 0.94 +/- 0.08(stat) +/- 0.12(sys) in p+p and d+Au collisions respectively, independent of p_T . The nuclear modification factor R_dA is 1.03 +/- 0.12(stat) +/- 0.21(sys) and 0.83 +/- 0.21(stat) +/- 0.17(sys) in minimum bias and central (0-20%) d+Au collisions, respectively.

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.

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.

We present a systematic study of charged pion and kaon interferometry in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. The kaon mean source radii are found to be larger than pion radii in the outward and longitudinal directions for the same transverse mass; this difference increases for more central collisions. The azimuthal-angle dependence of the radii was measured with respect to the second-order event plane and similar oscillations of the source radii were found for pions and kaons. Hydrodynamic models qualitatively describe the similar oscillations of the mean source radii for pions and kaons, but they do not fully describe the transverse-mass dependence of the oscillations.

We report on the yield of protons and anti-protons, as a function of centrality and transverse momentum, in Au+Au collisions at sqrt(s_NN) = 200 GeV measured at mid-rapidity by the PHENIX experiment at RHIC. In central collisions at intermediate transverse momenta (1.5 < p_T < 4.5 GeV/c) a significant fraction of all produced particles are protons and anti-protons. They show a centrality-scaling behavior different from that of pions. The p-bar/pion and p/pion ratios are enhanced compared to peripheral Au+Au, p+p, and electron+positron collisions. This enhancement is limited to p_T < 5 GeV/c as deduced from the ratio of charged hadrons to pi^0 measured in the range 1.5 < p_T < 9 GeV/c.

Transverse momentum spectra of electrons from Au+Au collisions at sqrt(s_NN) = 130 GeV have been measured by the PHENIX experiment at RHIC. The spectra show an excess above the background from photon conversions and light hadron decays. The electron signal is consistent with that expected from semi-leptonic decays of charm. The yield of the electron signal dN_e/dy for p_T > 0.8 GeV/c is 0.025 +/- 0.004 (stat.) +/- 0.010 (sys.) in central collisions, and the corresponding charm cross section is 380 +/- 60 (stat.) +/- 200 (sys.) micro barns per binary nucleon-nucleon collision.

Measurements of transverse-single-spin asymmetries ($A_{N}$) in $p$$+$$p$ collisions at $\sqrt{s}=$62.4 and 200 GeV with the PHENIX detector at RHIC are presented. At midrapidity, $A_{N}$ is measured for neutral pion and eta mesons reconstructed from diphoton decay, and at forward rapidities, neutral pions are measured using both diphotons and electromagnetic clusters. The neutral-pion measurement of $A_{N}$ at midrapidity is consistent with zero with uncertainties a factor of 20 smaller than previous publications, which will lead to improved constraints on the gluon Sivers function. At higher rapidities, where the valence quark distributions are probed, the data exhibit sizable asymmetries. In comparison with previous measurements in this kinematic region, the new data extend the kinematic coverage in $\sqrt{s}$ and $p_T$, and it is found that the asymmetries depend only weakly on $\sqrt{s}$. The origin of the forward $A_{N}$ is presently not understood quantitatively. The extended reach to higher $p_T$ probes the transition between transverse momentum dependent effects at low $p_T$ and multi-parton dynamics at high $p_T$.

The PHENIX Collaboration at the Relativistic Heavy Ion Collider has measured open heavy flavor production in Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}$=200 GeV through the measurement of electrons at midrapidity that originate from semileptonic decays of charm and bottom hadrons. In peripheral Cu$+$Cu collisions an enhanced production of electrons is observed relative to $p$$+$$p$ collisions scaled by the number of binary collisions. In the transverse momentum range from 1 to 5 GeV/$c$ the nuclear modification factor is $R_{AA}$$\sim$1.4. As the system size increases to more central Cu$+$Cu collisions, the enhancement gradually disappears and turns into a suppression. For $p_T>3$ GeV/$c$, the suppression reaches $R_{AA}$$\sim$0.8 in the most central collisions. The $p_T$ and centrality dependence of $R_{AA}$ in Cu$+$Cu collisions agree quantitatively with $R_{AA}$ in $d+$Au and Au$+$Au collisions, if compared at similar number of participating nucleons $\langle N_{\rm part} \rangle$.

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.

Measurements of the fractional momentum loss ($S_{\rm loss}\equiv{\delta}p_T/p_T$) of high-transverse-momentum-identified hadrons in heavy ion collisions are presented. Using $\pi^0$ in Au$+$Au and Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=62.4$ and 200 GeV measured by the PHENIX experiment at the Relativistic Heavy Ion Collider and and charged hadrons in Pb$+$Pb collisions measured by the ALICE experiment at the Large Hadron Collider, we studied the scaling properties of $S_{\rm loss}$ as a function of a number of variables: the number of participants, $N_{\rm part}$, the number of quark participants, $N_{\rm qp}$, the charged-particle density, $dN_{\rm ch}/d\eta$, and the Bjorken energy density times the equilibration time, $\varepsilon_{\rm Bj}\tau_{0}$. We find that the $p_T$ where $S_{\rm loss}$ has its maximum, varies both with centrality and collision energy. Above the maximum, $S_{\rm loss}$ tends to follow a power-law function with all four scaling variables. The data at $\sqrt{s_{_{NN}}}$=200 GeV and 2.76 TeV, for sufficiently high particle densities, have a common scaling of $S_{\rm loss}$ with $dN_{\rm ch}/d\eta$ and $\varepsilon_{\rm Bj}\tau_{0}$, lending insight on the physics of parton energy loss.

We report on two-particle azimuthal angle correlations between charged hadrons at forward/backward (deuteron/gold going direction) rapidity and charged hadrons at mid-rapidity in deuteron-gold (d+Au) and proton-proton (p+p) collisions at sqrt(s_NN) = 200 GeV. Jet structures are observed in the correlations which we quantify in terms of the conditional yield and angular width of away side partners. The kinematic region studied here samples partons in the gold nucleus carrying nucleon momentum fraction x~0.1 to x~0.01. Within this range, we find no x dependence of the jet structure in d+Au collisions.

The azimuthal distribution of identified pi^0 and inclusive photons has been measured in sqrt{s_{NN}} = 200 GeV Au+Au collisions with the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). The second harmonic parameter (v_2) was measured to describe the observed anisotropy of the azimuthal distribution. The measured inclusive photon v_2 is consistent with the value expected for the photons from hadron decay and is also consistent with the lack of direct photon signal over the measured p_T range 1-6 GeV/c. An attempt is made to extract v_2 of direct photons.

Two particle azimuthal correlation functions are presented for charged hadrons produced in Au + Au collisions at RHIC sqrt(s_NN) = 130 GeV. The measurements permit determination of elliptic flow without event-by-event estimation of the reaction plane. The extracted elliptic flow values v_2 show significant sensitivity to both the collision centrality and the transverse momenta of emitted hadrons, suggesting rapid thermalization and relatively strong velocity fields. When scaled by the eccentricity of the collision zone, epsilon, the scaled elliptic flow shows little or no dependence on centrality for charged hadrons with relatively low p_T. A breakdown of this epsilon scaling is observed for charged hadrons with p_T > 1.0 GeV/c for the most central collisions.

We present azimuthal angle correlations of intermediate transverse momentum (1-4 GeV/c) hadrons from {dijets} in Cu+Cu and Au+Au collisions at sqrt(s_NN) = 62.4 and 200 GeV. The away-side dijet induced azimuthal correlation is broadened, non-Gaussian, and peaked away from \Delta\phi=\pi in central and semi-central collisions in all the systems. The broadening and peak location are found to depend upon the number of participants in the collision, but not on the collision energy or beam nuclei. These results are consistent with sound or shock wave models, but pose challenges to Cherenkov gluon radiation models.

Inclusive transverse momentum spectra of eta mesons have been measured within p_T = 2-10 GeV/c at mid-rapidity by the PHENIX experiment in Au+Au collisions at sqrt(s_NN) = 200 GeV. In central Au+Au the eta yields are significantly suppressed compared to peripheral Au+Au, d+Au and p+p yields scaled by the corresponding number of nucleon-nucleon collisions. The magnitude, centrality and p_T dependence of the suppression is common, within errors, for eta and pi^0. The ratio of eta to pi^0 spectra at high p_T amounts to 0.40 < R_eta/pi^0 < 0.48 for the three systems in agreement with the world average measured in hadronic and nuclear reactions and, at large scaled momentum, in e^+e^- collisions.

PHENIX has measured the centrality dependence of mid-rapidity pion, kaon and proton transverse momentum distributions in d+Au and p+p collisions at sqrt(s_NN) = 200 GeV. The p+p data provide a reference for nuclear effects in d+Au and previously measured Au+Au collisions. Hadron production is enhanced in d+Au, relative to independent nucleon-nucleon scattering, as was observed in lower energy collisions. The nuclear modification factor for (anti) protons is larger than that for pions. The difference increases with centrality, but is not sufficient to account for the abundance of baryon production observed in central Au+Au collisions at RHIC. The centrality dependence in d+Au shows that the nuclear modification factor increases gradually with the number of collisions suffered by each participant nucleon. We also present comparisons with lower energy data as well as with parton recombination and other theoretical models of nuclear effects on particle production.