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

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

We report a high-statistics measurement of differential cross sections for the process gamma gamma -> pi^0 pi^0 in the kinematic range 0.6 GeV <= W <= 4.0 GeV and |cos theta*| <= 0.8, where W and theta* are the energy and pion scattering angle, respectively, in the gamma gamma center-of-mass system. Differential cross sections are fitted to obtain information on S, D_0, D_2, G_0 and G_2 waves. The G waves are important above W ~= 1.6 GeV. For W <= 1.6 GeV the D_2 wave is dominated by the f_2(1270) resonance while the S wave requires at least one additional resonance besides the f_0(980), which may be the f_0(1370) or f_0(1500). The differential cross sections are fitted with a simple parameterization to determine the parameters (the mass, total width and Gamma_{gamma gamma}B(f_0 -> pi^0 pi^0)) of this scalar meson as well as the f_0(980). The helicity 0 fraction of the f_2(1270) meson, taking into account interference for the first time, is also obtained.

Transverse mass spectra of protons emitted in Au+Au collisions at beam energies of 6, 8, and 10.8 GeV/nucleon have been measured as a function of collision centrality over a rapidity range 0.5<ylab<1.5. The spectra are well reproduced by Boltzmann distributions over the measured transverse mass region, which allows for extrapolation of the data to derive the rapidity density and apparent temperature of the emitting source. The shapes of the rapidity distributions suggest significant transparency or substantial longitudinal expansion in even the most central collisions at all three beam energies. The data are analyzed within a simple thermal source plus longitudinal expansion model.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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

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

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured electrons from heavy flavor (charm and bottom) decays for 0.3 < p_T < 9 GeV/c at midrapidity (|y| < 0.35) in Au+Au collisions at sqrt(s_NN) = 200 GeV. The nuclear modification factor R_AA relative to p+p collisions shows a strong suppression in central Au+Au collisions, indicating substantial energy loss of heavy quarks in the medium produced at RHIC. A large azimuthal anisotropy, v_2, with respect to the reaction plane is observed for 0.5 < p_T < 5 GeV/c indicating non-zero heavy flavor elliptic flow. Both R_AA and v_2 show a p_T dependence different from those of neutral pions. A comparison to transport models which simultaneously describe R_AA(p_T) and v_2(p_T) suggests that the viscosity to entropy density ratio is close to the conjectured quantum lower bound, i.e., near a perfect fluid.

The centrality dependence of transverse momentum distributions and yields for pi^+/-, K^+/-, p and p^bar in Au+Au collisions at sqrt(s_NN) = 200 GeV at mid-rapidity are measured by the PHENIX experiment at RHIC. We observe a clear particle mass dependence of the shapes of transverse momentum spectra in central collisions below ~ 2 GeV/c in p_T. Both mean transverse momenta and particle yields per participant pair increase from peripheral to mid-central and saturate at the most central collisions for all particle species. We also measure particle ratios of pi^-/pi^+, K^-/K^+, p^bar/p, K/pi, p/pi and p^bar/pi as a function of p_T and collision centrality. The ratios of equal mass particle yields are independent of p_T and centrality within the experimental uncertainties. In central collisions at intermediate transverse momenta ~ 1.5-4.5 GeV/c, proton and anti-proton yields constitute a significant fraction of the charged hadron production and show a scaling behavior different from that of pions.

The azimuthal anisotropy coefficients v_2 and v_4 of pi^0 and eta mesons are measured in Au+Au collisions at sqrt(s_NN)=200 GeV, as a function of transverse momentum p_T (1-14 GeV/c) and centrality. The extracted v_2 coefficients are found to be consistent between the two meson species over the measured p_T range. The ratio of v_4/v_2^2 for pi^0 mesons is found to be independent of p_T for 1-9 GeV/c, implying a lack of sensitivity of the ratio to the change of underlying physics with p_T. Furthermore, the ratio of v_4/v_2^2 is systematically larger in central collisions, which may reflect the combined effects of fluctuations in the initial collision geometry and finite viscosity in the evolving medium.

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.

The PHENIX experiment at RHIC has measured charged hadron yields at mid-rapidity over a wide range of transverse momentum (0.5 < p_T < 10 GeV/c) in Au+Au collisions at sqrt(s_NN)=200 GeV. The data are compared to pi^zero measurements from the same experiment. For both charged hadrons and neutral pions, the yields per nucleon-nucleon collision are significantly suppressed in central compared to peripheral and nucleon-nucleon collisions. The suppression sets in gradually and increases with increasing centrality of the collisions. Above 4-5 GeV/c in p_T, a constant and almost identical suppression of charged hadrons and pi^zeroes is observed. The p_T spectra are compared to published spectra from Au+Au at sqrt(s_NN)=130 in terms of x_t scaling. Central and peripheral pi^zero as well as peripheral charged spectra exhibit the same x_t scaling as observed in p+p data.

The transverse momentum dependence of the azimuthal anisotropy parameter v_2, the second harmonic of the azimuthal distribution, for electrons at mid-rapidity (|eta| < 0.35) has been measured with the PHENIX detector in Au+Au collisions at sqrt(s_NN) = 200 GeV. The measurement was made with respect to the reaction plane defined at high rapidities (|eta| = 3.1 -- 3.9). From the result we have measured the v_2 of electrons from heavy flavor decay after subtraction of the v_2 of electrons from other sources such as photon conversions and Dalitz decay from light neutral mesons. We observe a non-zero single electron v_2 with a 90% confidence level in the intermediate p_T region.

We report PHENIX measurements of the correlation of a trigger hadron at intermediate transverse momentum (2.5<p_{T,trig}<4 GeV/c), with associated mesons or baryons at lower p_{T,assoc}, in Au+Au collisions at sqrt(s_NN) = 200 GeV. The jet correlations for both baryons and mesons show similar shape alterations as a function of centrality, characteristic of strong modification of the away-side jet. The ratio of jet-associated baryons to mesons for this jet increases with centrality and p_{T,assoc} and, in the most central collisions, reaches a value similar to that for inclusive measurements. This trend is incompatible with in-vacuum fragmentation, but could be due to jet-like contributions from correlated soft partons which recombine upon hadronization.

The PHENIX experiement has measured the electron-positron pair mass spectrum from 0 to 8 GeV/c^2 in p+p collisions at sqrt(s)=200 GeV. The contributions from light meson decays to e^+e^- pairs have been determined based on measurements of hadron production cross sections by PHENIX. They account for nearly all e^+e^- pairs in the mass region below 1 GeV/c^2. The e^+e^- pair yield remaining after subtracting these contributions is dominated by semileptonic decays of charmed hadrons correlated through flavor conservation. Using the spectral shape predicted by PYTHIA, we estimate the charm production cross section to be 544 +/- 39(stat) +/- 142(syst) +/- 200(model) \mu b, which is consistent with QCD calculations and measurements of single leptons by PHENIX.

J/Psi production in p+p collisions at sqrt(s) = 200 GeV has been Measured in the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) over a rapidity range of -2.2 < y < 2.2 and a transverse momentum range of 0 < pT < 9 GeV/c. The statistics available allow a detailed measurement of both the pT and rapidity distributions and are sufficient to constrain production models. The total cross section times branching ratio determined for J/Psi production is B_{ll} sigma_pp^J/psi = 178 +/- 3(stat) +/- 53(syst) +/- 18(norm) nb.

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.

We report the measurement of the transverse momentum dependence of inclusive J/psi polarization in p+p collisions at sqrt(s)=200 GeV performed by the PHENIX Experiment at RHIC. The polarization is studied in the helicity, Gottfried-Jackson, and Collins-Soper frames for p_T < 5 GeV/c and |y| < 0.35. The J/psi 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.

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.

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.

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.

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$.

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.

We present transverse momentum (p_T) spectra of charged hadrons measured in deuteron-gold and nucleon-gold collisions at \sqrts = 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 the 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.

The PHENIX experiment has measured the production of neutral pions in Au+Au collisions at sqrt(s_NN)=200 GeV. The new data offer a fourfold increase in recorded luminosity, providing higher precision and a larger reach in transverse momentum, p_T, to 20 GeV/c. The production ratio of eta/pi^0 is 0.46+/-0.01(stat)+/-0.05(syst), constant with p_T and collision centrality. The observed ratio is consistent with earlier measurements, as well as with the p+p and d+Au values. The production of pi^0 is suppressed by a factor of 5, as in earlier findings. However, with the improved statistical precision a small but significant rise of the nuclear modification factor, R_AA, vs p_T, with a slope of 0.0106+/-^(0.0034)_(0.0029)[GeV/c]^-1, is discernible in central collisions. A phenomenological extraction of the average fractional parton energy loss shows a decrease with increasing p_T. To study the path length dependence of suppression, the pi^0 yield was measured at different angles with respect to the event plane; a strong azimuthal dependence of the pi^0 R_AA is observed. The data are compared to theoretical models of parton energy loss as a function of the path length, L, in the medium. Models based on pQCD are insufficient to describe the data, while a hybrid model utilizing pQCD for the hard interactions and AdS/CFT for the soft interactions is consistent with the data.

Correlations of charged hadrons of 1 < pT < 10 GeV/c with high pT direct photons and pi^ 0 mesons in the range 5 <pT < 15 GeV/c are used to study jet fragmentation in the photon+jet and di-jet channels, respectively. The magnitude of the partonic transverse momentum, kT, is obtained by comparing to a model incorporating a Gaussian kT smearing. The sensitivity of the associated charged hadron spectra to the underlying fragmentation function is tested and the data are compared to calculations using recent global fit results. The shape of the direct photon-associated hadron spectrum as well as its charge asymmetry are found to be consistent with a sample dominated by quark-gluon Compton scattering. No significant evidence of fragmentation photon correlated production is observed within experimental uncertainties.

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

It has been postulated that partonic orbital angular momentum can lead to a significant double-helicity dependence in the net transverse momentum of Drell-Yan dileptons produced in longitudinally polarized p+p collisions. Analogous effects are also expected for dijet production. If confirmed by experiment, this hypothesis, which is based on semi-classical arguments, could lead to a new approach for studying the contributions of orbital angular momentum to the proton spin. We report the first measurement of the double-helicity dependence of the dijet transverse momentum in longitudinally polarized p+p collisions at sqrt(s) = 200 GeV from data taken by the PHENIX experiment in 2005 and 2006. The analysis deduces the transverse momentum of the dijet from the widths of the near- and far-side peaks in the azimuthal correlation of the dihadrons. When averaged over the transverse momentum of the triggered particle, the difference of the root-mean-square of the dijet transverse momentum between like- and unlike-helicity collisions is found to be -37 +/- 88(stat) +/- 14(syst) MeV/c.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of K^0_S , \omega, \eta prime, and \phi mesons in p + p collisions at = 200 GeV. Measurements \omega and \phi production in different decay channels give consistent results. New results for the \phi are in agreement with previously published data and extend the measured pT coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, n and T, determining the high-pT and characterizing the low-pT regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter T extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.

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.