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 can serve to improve the precision of the extracted shear viscosity to entropy density ratio η/s.

Back-to-back hadron pair yields in d+Au and p+p collisions at √s(NN)=200 GeV were measured with the PHENIX detector at the Relativistic Heavy Ion Collider. Rapidity separated hadron pairs were detected with the trigger hadron at pseudorapidity |η|<0.35 and the associated hadron at forward rapidity (deuteron direction, 3.0<η<3.8). Pairs were also detected with both hadrons measured at forward rapidity; in this case, the yield of back-to-back hadron pairs in d+Au collisions with small impact parameters is observed to be suppressed by a factor of 10 relative to p+p collisions. The kinematics of these pairs is expected to probe partons in the Au nucleus with a low fraction x of the nucleon momenta, where the gluon densities rise sharply. The observed suppression as a function of nuclear thickness, p(T), and η points to cold nuclear matter effects arising at high parton densities.

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.

Large parity-violating longitudinal single-spin asymmetries A(L)(e+)=-0.86(-0.14)(+0.30) and A(L)(e-)= 0.88(-0.71)(+0.12) are observed for inclusive high transverse momentum electrons and positrons in polarized p+p collisions at a center-of-mass energy of sqrt[s]= 500 GeV with the PHENIX detector at RHIC. These e± come mainly from the decay of W± and Z0 bosons, and their asymmetries directly demonstrate parity violation in the couplings of the W± to the light quarks. The observed electron and positron yields were used to estimate W± boson production cross sections for the e± channels of σ(pp → W+ X) × BR(W+ → e+ ν(e))= 144.1 ± 21.2(stat)(-10.3)(+3.4)(syst) ± 21.6(norm)  pb, and σ(pp → W- X) × BR(W- → e- ν[over ¯](e))= 31.7 ± 12.1(stat)(-8.2)(+10.1)(syst) ± 4.8(norm)  pb.

This Letter presents the first measurement of event-by-event fluctuations of the elliptic flow parameter v_{2} in Au+Au collisions at sqrt[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.

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 sqrt[s_{NN}]=200 GeV over the uniquely broad longitudinal acceptance of the PHOBOS detector (-4<Deltaeta<2). A broadening of the away-side azimuthal correlation compared to elementary collisions is observed at all Deltaeta. 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 Deltaeta and known as the "ridge" is found to reach at least |Deltaeta| approximately 4. The ridge yield is largely independent of Deltaeta 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.

Bose-Einstein correlations of charged kaons are used to probe 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 of N_{part};{1/3} with a zero intercept. Imaging analysis shows a deviation from a Gaussian tail at r greater, similar10 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 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.

The double helicity asymmetry in neutral pion production for p_{T}=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 mu;{2}=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 Deltachi;{2}=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.

We present the first measurements of the pseudorapidity distribution of primary charged particles in Cu+Cu collisions as a function of collision centrality and energy, sqrt[s_{NN}]=22.4, 62.4, and 200 GeV, over a wide range of pseudorapidity, using the PHOBOS detector. A comparison of Cu+Cu and Au+Au results shows that the total number of produced charged particles and the rough shape (height and width) of the pseudorapidity distributions are determined by the number of nucleon participants. More detailed studies reveal that a more precise matching of the shape of the Cu+Cu and Au+Au pseudorapidity distributions over the full range of pseudorapidity occurs for the same N_{part}/2A rather than the same N_{part}. In other words, it is the collision geometry rather than just the number of nucleon participants that drives the detailed shape of the pseudorapidity distribution and its centrality dependence at RHIC energies.