Beam Energy and Centrality Dependence of Direct-Photon Emission from Ultrarelativistic Heavy-Ion Collisions

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_{ch}/dη reveals that the low-momentum (>1  GeV/c) direct-photon yield dN_{γ}^{dir}/dη is a smooth function of dN_{ch}/dη and can be well described as proportional to (dN_{ch}/dη)^{α} with α≈1.25. This scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and the 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.

Measurements of Elliptic and Triangular Flow in High-Multiplicity 3He+Au Collisions at √(s(NN))=200 GeV

We present the first measurement of elliptic (v(2)) and triangular (v(3)) flow in high-multiplicity (3)He+Au collisions at √(s(NN))=200  GeV. Two-particle correlations, where the particles have a large separation in pseudorapidity, are compared in (3)He+Au and in p+p collisions and indicate that collective effects dominate the second and third Fourier components for the correlations observed in the (3)He+Au system. The collective behavior is quantified in terms of elliptic v(2) and triangular v(3) anisotropy coefficients measured with respect to their corresponding event planes. The v(2) values are comparable to those previously measured in d+Au collisions at the same nucleon-nucleon center-of-mass energy. Comparisons with various theoretical predictions are made, including to models where the hot spots created by the impact of the three (3)He nucleons on the Au nucleus expand hydrodynamically to generate the triangular flow. The agreement of these models with data may indicate the formation of low-viscosity quark-gluon plasma even in these small collision systems.

Azimuthal-angle dependence of charged-pion-interferometry measurements with respect to second- and third-order event planes in Au+Au collisions at √[S(NN)]=200 GeV

Charged-pion-interferometry measurements were made with respect to the second- and third-order event plane for Au+Au collisions at sqrt[s_{NN}]=200  GeV. A strong azimuthal-angle dependence of the extracted Gaussian-source radii was observed with respect to both the second- and third-order event planes. The results for the second-order dependence indicate that the initial eccentricity is reduced during the medium evolution, which is consistent with previous results. In contrast, the results for the third-order dependence indicate that the initial triangular shape is significantly reduced and potentially reversed by the end of the medium evolution, and that the third-order oscillations are largely dominated by the dynamical effects from triangular flow.

Medium modification of jet fragmentation in Au+Au collisions at √[s(NN)]=200 GeV measured in direct photon-hadron correlations

The jet fragmentation function is measured with direct photon-hadron correlations in p+p and Au+Au collisions at √[s(NN)]=200 GeV. The p(T) of the photon is an excellent approximation to the initial p(T) of the jet and the ratio z(T)=p(T)(h)/p(T)(γ) is used as a proxy for the jet fragmentation function. A statistical subtraction is used to extract the direct photon-hadron yields in Au+Au collisions while a photon isolation cut is applied in p+p. I(AA), the ratio of hadron yield opposite the photon in Au+Au to that in p+p, indicates modification of the jet fragmentation function. Suppression, most likely due to energy loss in the medium, is seen at high z(T). The associated hadron yield at low z(T) is enhanced at large angles. Such a trend is expected from redistribution of the lost energy into increased production of low-momentum particles.

Observation of direct-photon collective flow in Au + Au collisions at √s(NN)] = 200 GeV

The second Fourier component v(2) of the azimuthal anisotropy with respect to the reaction plane is measured for direct photons at midrapidity and transverse momentum (p(T)) of 1-12 GeV/c in Au + Au collisions at √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, which is 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 underpredict the observed v(2).

Measurements of higher order flow harmonics in Au+Au collisions at √s(NN)=200 GeV

Flow coefficients ν(n) for n=2, 3, 4, characterizing the anisotropic collective flow in Au+Au collisions at √s(NN)=200 GeV, are measured relative to event planes Ψ(n), determined at large rapidity. We report ν(n) as a function of transverse momentum and collision centrality, and study the correlations among the event planes of different order n. The ν(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 can serve to improve the precision of the extracted shear viscosity to entropy density ratio η/s.

Cold nuclear matter effects on J/ψ yields as a function of rapidity and nuclear geometry in d+A collisions at sqrt[s(NN)]=200 GeV

We present measurements of J/ψ yields in d+Au collisions at sqrt[s(NN)]=200  GeV recorded by the PHENIX experiment and compare them with yields in p+p collisions at the same energy per nucleon-nucleon collision. The measurements cover a large kinematic range in J/ψ rapidity (-2.2<y<2.4) with high statistical precision and are compared with two theoretical models: one with nuclear shadowing combined with final state breakup and one with coherent gluon saturation effects. In order to remove model dependent systematic uncertainties we also compare the data to a simple geometric model. The forward rapidity data are inconsistent with nuclear modifications that are linear or exponential in the density weighted longitudinal thickness, such as those from the final state breakup of the bound state.

Azimuthal anisotropy of π⁰ production in Au+Au collisions at sqrt((s)NN)=200 GeV: path-length dependence of jet quenching and the role of initial geometry

We have measured the azimuthal anisotropy of π⁰ production for 1<p(T)<18  GeV/c for Au+Au collisions at sqrt((s)NN)=200  GeV. The observed anisotropy shows a gradual decrease for 3≲p(T)≲7-10  GeV/c, but remains positive beyond 10  GeV/c. The magnitude of this anisotropy is underpredicted, up to at least ∼10  GeV/c, by current perturbative QCD (PQCD) energy-loss model calculations. An estimate of the increase in anisotropy expected from initial-geometry modification due to gluon saturation effects and fluctuations is insufficient to account for this discrepancy. Calculations that implement a path-length dependence steeper than what is implied by current PQCD energy-loss models show reasonable agreement with the data.

Elliptic and hexadecapole flow of charged hadrons in Au+Au collisions at sq.rt(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 number of participant nucleons (N(part)) for Au+Au collisions at sq.rt(s(NN))=200  GeV. The v(2,4) measurements at pseudorapidity |η|≤0.35, obtained with four separate reaction-plane detectors positioned in the range 1.0<|η|<3.9, show good agreement, indicating the absence of significant Δη-dependent nonflow correlations. 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.

Transition in yield and azimuthal shape modification in dihadron correlations in relativistic heavy ion collisions

Hard-scattered parton probes produced in collisions of large nuclei indicate large partonic energy loss, possibly with collective produced-medium response to the lost energy. We present measurements of π^{0} trigger particles at transverse momenta p{T}{t}=4-12  GeV/c and associated charged hadrons (p{T}{a}=0.5-7  GeV/c) vs relative azimuthal angle Δϕ in Au+Au and p+p collisions at sqrt[s{NN}]=200  GeV. The Au+Au distribution at low p{T}{a}, whose shape has been interpreted as a medium effect, is modified for p{T}{t}<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{T}{a}, which quantitatively challenges some medium response models. The associated yield of hadrons opposing the trigger particle in Au+Au relative to p+p (I{AA}) is suppressed at high p{T} (I{AA}≈0.35-0.5), but less than for inclusive suppression (R{AA}≈0.2).

Event-by-event fluctuations of azimuthal particle anisotropy in Au+Au collisions at square root(sNN) = 200 GeV

This Letter presents the first measurement of event-by-event fluctuations of the elliptic flow parameter v(2) in Au+Au collisions at square root(s(NN))=200 GeV as a function of collision centrality. The relative nonstatistical fluctuations of the v(2) parameter are found to be approximately 40%. The results, including contributions from event-by-event elliptic flow fluctuations and from azimuthal correlations that are unrelated to the reaction plane (nonflow correlations), establish an upper limit on the magnitude of underlying elliptic flow fluctuations. This limit is consistent with predictions based on spatial fluctuations of the participating nucleons in the initial nuclear overlap region. These results provide important constraints on models of the initial state and hydrodynamic evolution of relativistic heavy ion collisions.

Enhanced production of direct photons in Au + Au collisions at square root(S(NN)) = 200 GeV and implications for the initial temperature

The production of e+ e- pairs for m(e+ e-)<0.3 GeV/c2 and 1<p(T)<5 GeV/c is measured in p+p and Au+Au collisions at square root(S(NN))=200 GeV. An enhanced yield above hadronic sources is observed. Treating the excess as photon internal conversions, the invariant yield of direct photons is deduced. In central Au+Au collisions, the excess of the direct photon yield over p+p is exponential in transverse momentum, with an inverse slope T=221+/-19(stat)+/-19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) approximately 300-600 MeV at times of approximately 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 approximately 170 MeV.

High transverse momentum triggered correlations over a large pseudorapidity acceptance in Au + Au collisions at square root(s(NN)) = 200 GeV

A measurement of two-particle correlations with a high transverse momentum trigger particle (p(T)(trig) > 2.5 GeV/c) is presented for Au+Au collisions at square root(s(NN)) = 200 GeV over the uniquely broad longitudinal acceptance of the PHOBOS detector (-4 < Delta eta < 2). A broadening of the away-side azimuthal correlation compared to elementary collisions is observed at all Delta eta. As in p+p collisions, the near side is characterized by a peak of correlated partners at small angle relative to the trigger particle. However, in central Au+Au collisions an additional correlation extended in Delta eta and known as the "ridge" is found to reach at least |Delta eta| approximately = 4. The ridge yield is largely independent of Delta eta over the measured range, and it decreases towards more peripheral collisions. For the chosen (p(T)(trig) cut, the ridge yield is consistent with zero for events with less than roughly 100 participating nucleons.

Measurement of bottom versus charm as a function of transverse momentum with electron-hadron correlations in p + p collisions at square root of s = 200 GeV

The momentum distribution of electrons from semileptonic decays of charm and bottom quarks for midrapidity |y|<0.35 in p+p collisions at square root of s=200 GeV is measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range 2<pT<7 GeV/c. The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is determined using partial D/D-->e(+/-)K(-/+)X (K unidentified) reconstruction. It is found that the yield of electrons from bottom becomes significant above 4 GeV/c in pT. A fixed-order-plus-next-to-leading-log perturbative quantum chromodynamics calculation agrees with the data within the theoretical and experimental uncertainties. The extracted total bottom production cross section at this energy is sigma(bb)=3.2(-1.1)(+1.2)(stat)(-1.3)(+1.4)(syst)mub.

Gluon-spin contribution to the proton spin from the double-helicity asymmetry in inclusive pi0 production in polarized p+p collisions at [sqrt]s=200 GeV

The double helicity asymmetry in neutral pion production for pT=1 to 12 GeV/c was measured with the PHENIX experiment to access the gluon-spin contribution, DeltaG, to the proton spin. Measured asymmetries are consistent with zero, and at a theory scale of micro2=4 GeV2 a next to leading order QCD analysis gives DeltaG[0.02,0.3]=0.2, with a constraint of -0.7<DeltaG[0.02,0.3]<0.5 at Deltachi2=9 (approximately 3sigma) for the sampled gluon momentum fraction (x) range, 0.02 to 0.3. The results are obtained using predictions for the measured asymmetries generated from four representative fits to polarized deep inelastic scattering data. We also consider the dependence of the DeltaG constraint on the choice of the theoretical scale, a dominant uncertainty in these predictions.

Suppression pattern of neutral pions at high transverse momentum in Au + Au collisions at sqrt[sNN]=200 GeV and constraints on medium transport coefficients

For Au + Au collisions at 200 GeV, we measure neutral pion production with good statistics for transverse momentum, pT, up to 20 GeV/c. A fivefold suppression is found, which is essentially constant for 5 < pT < 20 GeV/c. Experimental uncertainties are small enough to constrain any model-dependent parametrization for the transport coefficient of the medium, e.g., q in the parton quenching model. The spectral shape is similar for all collision classes, and the suppression does not saturate in Au + Au collisions.

Onset of pi(0) suppression studied in Cu+Cu collisions at sqrt S NN=22.4, 62.4, and 200 GeV

Neutral pion transverse momentum (p(T)) spectra at midrapidity (|y| less than or approximately 0.35) were measured in Cu+Cu collisions at sqrt[s(NN)]=22.4, 62.4, and 200 GeV. Relative to pi_(0) yields in p+p collisions scaled by the number of inelastic nucleon-nucleon collisions (N(coll) the pi_(0) yields for p(T) more than or approximately 2 GeV/c in central Cu+Cu collisions are suppressed at 62.4 and 200 GeV 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 and 200 GeV, while the enhancement at 22.4 GeV is consistent with nuclear modifications in the initial state alone.

J/psi production in sqrt s_NN=200 GeV Cu+Cu collisions

Yields for J/psi production in Cu+Cu collisions at sqrt s_NN=200 GeV have been measured over the rapidity range |y|<2.2 and compared with results in p+p and Au+Au collisions at the same energy. The Cu+Cu data offer greatly improved precision over existing Au+Au data for J/psi production in collisions with small to intermediate numbers of participants, in the range where the quark-gluon plasma transition threshold is predicted to lie. Cold nuclear matter estimates based on ad hoc fits to d+Au data describe the Cu+Cu data up to N_part approximately 50, corresponding to a Bjorken energy density of at least 1.5 GeV/fm(3).

System size, energy, pseudorapidity, and centrality dependence of elliptic flow

This Letter presents measurements of the elliptic flow of charged particles as a function of pseudorapidity and centrality from Cu-Cu collisions at 62.4 and 200 GeV using the PHOBOS detector at the Relativistic Heavy Ion Collider. The elliptic flow in Cu-Cu collisions is found to be significant even for the most central events. For comparison with the Au-Au results, it is found that the detailed way in which the collision geometry (eccentricity) is estimated is of critical importance when scaling out system-size effects. A new form of eccentricity, called the participant eccentricity, is introduced which yields a scaled elliptic flow in the Cu-Cu system that has the same relative magnitude and qualitative features as that in the Au-Au system.

System size and energy dependence of jet-induced hadron pair correlation shapes in Cu+Cu and Au+Au collisions at square root sNN=200 and 62.4 GeV

We present azimuthal angle correlations of intermediate transverse momentum (1-4 GeV/c) hadrons from dijets in Cu+Cu and Au+Au collisions at square root sNN=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 semicentral 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.

J/psi production versus transverse momentum and rapidity in p+p collisions at square root s=200 GeV

J/psi production in p+p collisions at square root s=200 GeV has been measured by the PHENIX experiment at the BNL Relativistic Heavy Ion Collider over a rapidity range of -2.2<y<2.2 and a transverse momentum range of 0<pT<9 GeV/c. The size of the present data set allows a detailed measurement of both the pT and the rapidity distributions and is sufficient to constrain production models. The total cross section times the branching ratio is Bll sigma pp J/psi=178+/-3stat+/-53sys+/-18norm nb.

J/psi production versus centrality, transverse momentum, and rapidity in Au+Au collisions at square root sNN=200 GeV

The PHENIX experiment at the BNL Relativistic Heavy Ion Collider (RHIC) has measured J/psi production for rapidities -2.2<y<2.2 in Au+Au collisions at square root sNN=200 GeV. The J/psi invariant yield and nuclear modification factor RAA as a function of centrality, transverse momentum, and rapidity are reported. A suppression of J/psi relative to binary collision scaling of proton-proton reaction yields is observed. Models which describe the lower energy J/psi data at the CERN Super Proton Synchrotron invoking only J/psi destruction based on the local medium density predict a significantly larger suppression at RHIC and more suppression at midrapidity than at forward rapidity. Both trends are contradicted by our data.

Scaling properties of azimuthal anisotropy in Au+Au and Cu+Cu Collisions at sqrt[s NN]=200 GeV

Differential measurements of elliptic flow (v2) for Au+Au and Cu+Cu collisions at sqrt[sNN]=200 GeV are used to test and validate predictions from perfect fluid hydrodynamics for scaling of v2 with eccentricity, system size, and transverse kinetic energy (KE T). For KE T identical with mT-m up to approximately 1 GeV the scaling is compatible with hydrodynamic expansion of a thermalized fluid. For large values of KE T mesons and baryons scale separately. Quark number scaling reveals a universal scaling of v2 for both mesons and baryons over the full KE T range for Au+Au. For Au+Au and Cu+Cu the scaling is more pronounced in terms of KE T, rather than transverse momentum.

Measurement of high-pT single electrons from heavy-flavor decays in p + p collisions at square root of s = 200 GeV

The momentum distribution of electrons from decays of heavy flavor (charm and bottom) for midrapidity absolute value of y < 0.35 in p + p collisions at square root of s = 200 GeV has been measured by the PHENIX experiment at the BNL Relativistic Heavy Ion Collider over the transverse momentum range 0.3 < pT < 9 GeV/c. Two independent methods have been used to determine the heavy-flavor yields, and the results are in good agreement with each other. A fixed-order-plus-next-to-leading-log perturbative QCD calculation agrees with the data within the theoretical and experimental uncertainties, with the data/theory ratio of 1.71+/-0.02stat+/-0.18sys for 0.3 < pT < 9 GeV/c. The total charm production cross section at this energy has also been deduced to be sigma cc = 567+/-57stat+/-193sys microb.

Energy dependence of directed flow over a wide range of pseudorapidity in Au + Au collisions at the BNL Relativistic Heavy Ion Collider

We report on measurements of directed flow as a function of pseudorapidity in Au + Au collisions at energies of square root of SNN = 19.6, 62.4, 130 and 200 GeV as measured by the PHOBOS detector at the BNL Relativistic Heavy Ion Collider. These results are particularly valuable because of the extensive, continuous pseudorapidity coverage of the PHOBOS detector. There is no significant indication of structure near midrapidity and the data surprisingly exhibit extended longitudinal scaling similar to that seen for elliptic flow and charged particle pseudorapidity density.

Energy dependence of elliptic flow over a large pseudorapidity range in Au+Au collisions at the BNL relativistic heavy ion collider

This Letter describes the measurement of the energy dependence of elliptic flow for charged particles in Au+Au collisions using the PHOBOS detector at the Relativistic Heavy Ion Collider. Data taken at collision energies of square root of s(NN)=19.6, 62.4, 130, and 200 GeV are shown over a wide range in pseudorapidity. These results, when plotted as a function of eta(')=|eta|-y(beam), scale with approximate linearity throughout eta('), implying no sharp changes in the dynamics of particle production as a function of pseudorapidity or increasing beam energy.

Centrality dependence of charged hadron transverse momentum spectra in Au+Au collisions from sqrt[s(NN)]=62.4 to 200 GeV

We have measured transverse momentum distributions of charged hadrons produced in Au+Au collisions at sqrt[s(NN)]=62.4 GeV. The spectra are presented for transverse momenta 0.25<p(T)<4.5 GeV/c, in a pseudorapidity range of 0.2<eta<1.4. The nuclear modification factor R(AA) is calculated relative to p+p data at the same collision energy as a function of collision centrality. For 2<p(T)<4.5 GeV/c, R(AA) is found to be significantly larger than in Au+Au collisions at sqrt[s(NN)]= 130 and 200 GeV. In contrast to the large change in R(AA), we observe a very similar centrality evolution of the p(T) spectra at sqrt[s(NN)]=62.4 and 200 GeV. The dynamical origin of this surprising factorization of energy and centrality dependence of particle production in heavy-ion collisions remains to be understood.

Pseudorapidity distribution of charged particles in d+Au collisions at sqrt[sNN]=200 GeV

The measured pseudorapidity distribution of primary charged particles in minimum-bias d+Au collisions at sqrt[s(NN)]=200 GeV is presented for the first time. This distribution falls off less rapidly in the gold direction as compared to the deuteron direction. The average value of the charged particle pseudorapidity density at midrapidity is <dN(ch)/d eta>|eta|< or =0.6)=9.4+/-0.7(syst) and the integrated primary charged particle multiplicity in the measured region is 82+/-6(syst). Estimates of the total charged particle production, based on extrapolations outside the measured pseudorapidity region, are also presented. The pseudorapidity distribution, normalized to the number of participants in d+Au collisions, is compared to those of Au+Au and p+(-)p systems at the same energy. The d+Au distribution is also compared to the predictions of the parton saturation model, as well as microscopic models.

Centrality dependence of charged-hadron transverse-momentum spectra in d+Au collisions at sqrt[s(NN)]=200 GeV

We have measured transverse momentum distributions of charged hadrons produced in d+Au collisions at sqrt[s(NN)]=200 GeV. The spectra were obtained for transverse momenta 0.25<p(T)<6.0 GeV/c, in a pseudorapidity range of 0.2<eta<1.4 in the deuteron direction. The evolution of the spectra with collision centrality is presented in comparison to p+pmacr; collisions at the same collision energy. With increasing centrality, the yield at high transverse momenta increases more rapidly than the overall particle density, leading to a strong modification of the spectral shape. This change in spectral shape is qualitatively different from observations in Au+Au collisions at the same energy. The results provide important information for discriminating between different models for the suppression of high-p(T) hadrons observed in Au+Au collisions.

Significance of the fragmentation region in ultrarelativistic heavy-ion collisions

We present measurements of the pseudorapidity distribution of primary charged particles produced in Au+Au collisions at three energies, sqrt[s(NN)]=19.6, 130, and 200 GeV, for a range of collision centrali-ties. The distribution narrows for more central collisions and excess particles are produced at high pseudorapidity in peripheral collisions. For a given centrality, however, the distributions are found to scale with energy according to the "limiting fragmentation" hypothesis. The universal fragmentation region described by this scaling grows in pseudorapidity with increasing collision energy, extending well away from the beam rapidity and covering more than half of the pseudorapidity range over which particles are produced. This approach to a universal limiting curve appears to be a dominant feature of the pseudorapidity distribution and therefore of the total particle production in these collisions.