Cold Nuclear Matter Effects on J/psi Yields as a Function of Rapidity and Nuclear Geometry in d plus A Collisions at root S(NN)=200 GeV

We present measurements of J/psi yields in d + Au collisions at root 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/psi 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.

Elliptic and Hexadecapole Flow of Charged Hadrons in Au plus Au Collisions at root s(NN)=200 GeV

Differential measurements of the elliptic (upsilon(2)) and hexadecapole (upsilon(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 root s(NN) = 200 GeV/ The upsilon(2,4) measurements at pseudorapidity vertical bar eta vertical bar <= 0.35, obtained with four separate reaction-plane detectors positioned in the range 1.0 < vertical bar eta vertical bar < 3.9, show good agreement, indicating the absence of significant Delta eta-dependent nonflow correlations. Sizable values for upsilon(4)(p(T)) are observed with a ratio upsilon(4)(p(T), N(part))/upsilon(2)(2)(p(T), N(part)) approximate to 0.8 for 50 less than or similar to N(part) less than or similar to 200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N(part) greater than or similar to 200 this ratio increases up to 1.7 in the most central collisions.

Transverse momentum dependence of eta meson suppression in Au plus Au collisions at root s(NN)=200 GeV

New measurements by the PHENIX experiment at the Relativistic Heavy Ion Collider for. production at midrapidity as a function of transverse momentum ((PT)) and collision centrality in root s(NN) = 200 GeV Au + Au and p + p collisions are presented. They indicate nuclear modification factors (R(AA)) which are similar in both magnitude and trend to those found in earlier pi(0) measurements. Linear fits to R(AA) as a function of (PT) in 5-20 GeV/c 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 (PT); thus, along with the extended (PT) range of these data they can provide additional constraints for theoretical modeling and the extraction of transport properties.

Transverse momentum dependence of J/psi polarization at midrapidity in p plus p collisions at root s=200 GeV

We report the measurement of the transverse momentum dependence of inclusive J/psi polarization in p + p collisions at root s = 200 GeV performed by the PHENIX Experiment at the Relativistic Heavy Ion Collider. The J/psi polarization is studied in the helicity, Gottfried-Jackson, and Collins-Soper frames for p(T) < 5 GeV/c and vertical bar y vertical bar < 0.35. The polarization in the helicity and Gottfried-Jackson frames is consistent with zero for all transverse momenta, with a slight (1.8 sigma) trend towards longitudinal polarization for transverse momenta above 2 GeV/c. No conclusion is allowed due to the limited acceptance in the Collins-Soper frame and the uncertainties of the current data. The results are compared to observations for other collision systems and center of mass energies and to different quarkonia production models.

Enhanced Production of Direct Photons in Au plus Au Collisions at root s(NN)=200 GeV and Implications for the Initial Temperature

The production of e(+)e(-) pairs for m(e+e-) < 0.3 GeV/c(2) and 1< p(T) < 5 GeV/c is measured in p + p and Au + Au collisions at 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) similar to 300-600 MeV at times of similar to 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 similar to 170 MeV.

Detailed measurement of the e(+)e(-) pair continuum in p plus p and Au plus Au collisions at root s(NN)=200 GeV and implications for direct photon production

PHENIX has measured the e(+)e(-) pair continuum in root s(NN) = 200 GeV Au+Au and p+p collisions over a wide range of mass and transverse momenta. The e(+)e(-) yield is compared to the expectations from hadronic sources, based on PHENIX measurements. In the intermediate-mass region, between the masses of the phi and the J/psi meson, the yield is consistent with expectations from correlated c (c) over bar production, although other mechanisms are not ruled out. In the low-mass region, below the phi, the p+p inclusive mass spectrum is well described by known contributions from light meson decays. In contrast, the Au+Au minimum bias inclusive mass spectrum in this region shows an enhancement by a factor of 4.7 +/- 0.4(stat) +/- 1.5(syst) +/- 0.9(model). At low mass (m(ee) < 0.3 GeV/c(2)) and high p(T) (1 < p(T) < 5 GeV/c) an enhanced e(+)e(-) pair yield is observed that is consistent with production of virtual direct photons. This excess is used to infer the yield of real direct photons. In central Au+Au collisions, the excess of the direct photon yield over the p+p is exponential in p(T), with inverse slope T = 221 +/- 19(stat) +/- 19(syst) MeV. Hydrodynamical models with initial temperatures ranging from T(init) similar or equal to 300-600 MeV at times of 0.6-0.15 fm/c after the collision are in qualitative agreement with the direct photon data in Au+Au. For low p(T) < 1 GeV/c the low-mass region shows a further significant enhancement that increases with centrality and has an inverse slope of T similar or equal to 100 MeV. Theoretical models underpredict the low-mass, low-p(T) enhancement.

High-p(T) pi(0) production with respect to the reaction plane in Au plus Au collisions at s(NN)=200 GeV

Measurements of the azimuthal anisotropy of high-p(T) neutral pion (pi(0)) production in Au+Au collisions at s(NN)=200 GeV by the PHENIX experiment are presented. The data included in this article were collected during the 2004 Relativistic Heavy Ion Collider running period and represent approximately an order of magnitude increase in the number of analyzed events relative to previously published results. Azimuthal angle distributions of pi(0) mesons detected in the PHENIX electromagnetic calorimeters are measured relative to the reaction plane determined event-by-event using the forward and backward beam-beam counters. Amplitudes of the second Fourier component (v(2)) of the angular distributions are presented as a function of pi(0) transverse momentum (p(T)) for different bins in collision centrality. Measured reaction plane dependent pi(0) yields are used to determine the azimuthal dependence of the pi(0) suppression as a function of p(T), R(AA)(Delta phi,p(T)). A jet-quenching motivated geometric analysis is presented that attempts to simultaneously describe the centrality dependence and reaction plane angle dependence of the pi(0) suppression in terms of the path lengths of hypothetical parent partons in the medium. This set of results allows for a detailed examination of the influence of geometry in the collision region and of the interplay between collective flow and jet-quenching effects along the azimuthal axis.

Charged Kaon Interferometric Probes of Space-Time Evolution in Au plus Au Collisions at s(NN)=200 GeV

Bose-Einstein correlations of charged kaons are used to probe Au+Au collisions at 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 of N(part)(1/3) with a zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r greater than or similar to 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.

Suppression Pattern of Neutral Pions at High Transverse Momentum in Au plus Au Collisions at root S(NN)=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, p(T), up to 20 GeV/c. A fivefold suppression is found, which is essentially constant for 5 < p(T) < 20 GeV/c. Experimental uncertainties are small enough to constrain any model-dependent parametrization for the transport coefficient of the medium, e. g., <(q) over cap > in the parton quenching model. The spectral shape is similar for all collision classes, and the suppression does not saturate in Au + Au collisions.

Event structure and double helicity asymmetry in jet production from polarized p plus p collisions at root s=200 GeV

We report on the event structure and double helicity asymmetry (A(LL)) of jet production in longitudinally polarized p + p collisions at root s = 200 GeV. Photons and charged particles were measured by the PHENIX experiment at midrapidity vertical bar eta vertical bar < 0.35 with the requirement of a high-momentum (> 2 GeV/c) photon in the event. Event structure, such as multiplicity, p(T) density and thrust in the PHENIX acceptance, were measured and compared with the results from the PYTHIA event generator and the GEANT detector simulation. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet A(LL), photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster pT sum (p(T)(reco)). The effect of detector response and the underlying events on p(T)(reco) was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the next-to-leading-order and perturbative quantum chromodynamics jet production cross section. For 4< p(T)(reco) < 12 GeV/c with an average beam polarization of < P > = 49% we measured Lambda(LL) = -0.0014 +/- 0.0037(stat) at the lowest p(T)(reco) bin (4-5 GeV= c) and -0.0181 +/- 0.0282(stat) at the highest p(T)(reco) bin (10-12 GeV= c) with a beam polarization scale error of 9.4% and a pT scale error of 10%. Jets in the measured p(T)(reco) range arise primarily from hard-scattered gluons with momentum fraction 0: 02 < x < 0: 3 according to PYTHIA. The measured A(LL) is compared with predictions that assume various Delta G(x) distributions based on the Gluck-Reya-Stratmann-Vogelsang parameterization. The present result imposes the limit -a.1 < integral(0.3)(0.02) dx Delta G(x, mu(2) = GeV2) < 0.4 at 95% confidence level or integral(0.3)(0.002) dx Delta G(x, mu(2) = 1 GeV2) < 0.5 at 99% confidence level.

Charged hadron multiplicity fluctuations in Au+Au and Cu+Cu collisions from root S(NN) = 22.5 to 200 GeV

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 s(NN)=62.4 and 200 GeV, and Cu+Cu collisions at 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 nondynamical 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. A comparison of the data with a model where hadrons are independently emitted from a number of hadron clusters suggests that the mean number of hadrons per cluster is small in heavy ion collisions.

Cold nuclear matter effects on J/psi production as constrained by deuteron-gold measurements at root S(NN)=200 GeV

We present a new analysis of J/psi production yields in deuteron-gold collisions at root s(NN) =200 GeV using data taken from the PHENIX experiment in 2003 and previously published in S. S. Adler [Phys. Rev. Lett 96, 012304 (2006)]. The high statistics proton-proton J/psi data taken in 2005 are used to improve the baseline measurement and thus construct updated cold nuclear matter modification factors (R(dAu)). A suppression of J/psi in cold nuclear matter is observed as one goes forward in rapidity (in the deuteron-going direction), corresponding to a region more sensitive to initial-state low-x gluons in the gold nucleus. The measured nuclear modification factors are compared to theoretical calculations of nuclear shadowing to which a J/psi (or precursor) breakup cross section is added. Breakup cross sections of sigma(breakup)=2.8(-1.4)(+1.7) (2.2(-1.5)(+1.6)) mb are obtained by fitting these calculations to the data using two different models of nuclear shadowing. These breakup cross-section values are consistent within large uncertainties with the 4.2 +/- 0.5 mb determined at lower collision energies. Projecting this range of cold nuclear matter effects to copper-copper and gold-gold collisions reveals that the current constraints are not sufficient to firmly quantify the additional hot nuclear matter effect.

Centrality dependence of charged hadron production in deuteron plus gold and nucleon plus gold collisions at root S(NN)=200 GeV

We present transverse momentum (p(T)) spectra of charged hadrons measured in deuteron-gold and nucleon-gold collisions at root s(NN)=200 GeV for four centrality classes. Nucleon-gold collisions were selected by tagging events in which a spectator nucleon was observed in one of two forward rapidity detectors. The spectra and yields were investigated as a function of the number of binary nucleon-nucleon collisions, nu, suffered by deuteron nucleons. A comparison of charged particle yields to those in p+p collisions show that yield per nucleon-nucleon collision saturates with nu for high momentum particles. We also present the charged hadron to neutral pion ratios as a function of p(T).

KINETICS VERSUS HYDRODYNAMICS: GENERALIZATION OF LANDAU/COOPER-FRYE PRESCRIPTION FOR FREEZE-OUT

The problem of spectra formation in hydrodynamic approach to A + A collisions is considered within the Boltzmann equations. It is shown analytically and illustrated by numerical calculations that the particle momentum spectra can be presented in the Cooper-R-ye form despite freeze-out is not sharp and has the finite temporal width. The latter is equal to the inverse of the particle collision rate at points (t(sigma) (r, p), r) of the maximal emission at a fixed momentum p. The set of these points forms the hypersurfaces t(sigma)(r,p) which strongly depend on the values of p and typically do not enclose completely the initially dense matter. This is an important difference from the standard Cooper-Frye prescription (CFp), with a common freeze-out hypersurface for all p, that affects significantly the predicted spectra. Also, the well known problem of CFp as for negative contributions to the spectra from non-space-like parts of the freeze-out hypersurface is naturally eliminated in this improved prescription.

Directed flow at midrapidity in event-by-event hydrodynamics

Fluctuations in the initial geometry of a nucleus-nucleus collision have been recently shown to result in a new type of directed flow (v(1)) that, unlike the usual directed flow, is also present at midrapidity. We compute this new v(1) versus transverse momentum and centrality for Au-Au collisions at RHIC using the hydrodynamic code NeXSPheRIO. We find that the event plane of v(1) is correlated with the angle of the initial dipole of the distribution, as predicted, though with a large dispersion. It is uncorrelated with the reaction plane. Our results are in excellent agreement with results inferred from STAR correlation data.