Transverse momentum spectra of charged hadrons with p_T < 8 GeV/c and neutral pions with p_T < 10 GeV/c have been measured at mid-rapidity by the PHENIX experiment at RHIC in d+Au collisions at sqrt(s_NN) = 200 GeV. The measured yields are compared to those in p+p collisions at the same sqrt(s_NN) scaled up by the number of underlying nucleon-nucleon collisions in d+Au. The yield ratio does not show the suppression observed in central Au+Au collisions at RHIC. Instead, there is a small enhancement in the yield of high momentum particles.

We present two-dimensional (2D) two-particle angular correlations on relative pseudorapidity $\eta$ and azimuth $\phi$ for charged particles from Au-Au collisions at $\sqrt{s_{\rm NN}} = 62$ and 200 GeV with transverse momentum $p_t \geq 0.15$ GeV/$c$, $|\eta| \leq 1$ and $2\pi$ azimuth. Observed correlations include a {same-side} (relative azimuth $< \pi/2$) 2D peak, a closely-related away-side azimuth dipole, and an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton-proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling in Au-Au collisions until mid-centrality where a transition to a qualitatively different centrality trend occurs within a small centrality interval. Above the transition point the number of same-side and away-side correlated pairs increases rapidly {relative to} binary-collision scaling, the $\eta$ width of the same-side 2D peak also increases rapidly ($\eta$ elongation) and the $\phi$ width actually decreases significantly. Those centrality trends are more remarkable when contrasted with expectations of jet quenching in a dense medium. Observed centrality trends are compared to {\sc hijing} predictions and to the expected trends for semihard parton scattering and fragmentation in a thermalized opaque medium. We are unable to reconcile a semihard parton scattering and fragmentation origin for the observed correlation structure and centrality trends with heavy ion collision scenarios which invoke rapid parton thermalization. On the other hand, if the collision system is effectively opaque to few-GeV partons the observations reported here would be inconsistent with a minijet picture.

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

Results on high transverse momentum charged particle emission with respect to the reaction plane are presented for Au+Au collisions at $\sqrt{s_{_{NN}}}$= 200 GeV. Two- and four-particle correlations results are presented as well as a comparison of azimuthal correlations in Au+Au collisions to those in $p+p$ at the same energy. Elliptic anisotropy, $v_2$, is found to reach its maximum at $p_t \sim 3$ GeV/c, then decrease slowly and remain significant up to $p_t\approx 7$ -- 10 GeV/c. Stronger suppression is found in the back-to-back high-$p_t$ particle correlations for particles emitted out-of-plane compared to those emitted in-plane. The centrality dependence of $v_2$ at intermediate $p_t$ is compared to simple models based on jet quenching.

The results from the STAR Collaboration on directed flow (v_1), elliptic flow (v_2), and the fourth harmonic (v_4) in the anisotropic azimuthal distribution of particles from Au+Au collisions at sqrtsNN = 200 GeV are summarized and compared with results from other experiments and theoretical models. Results for identified particles are presented and fit with a Blast Wave model. Different anisotropic flow analysis methods are compared and nonflow effects are extracted from the data. For v_2, scaling with the number of constituent quarks and parton coalescence is discussed. For v_4, scaling with v_2^2 and quark coalescence is discussed.

The azimuthal anisotropic flow of identified and unidentified charged particles has been systematically studied in Cu+Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV for harmonics $n=$ 1-4 in the pseudorapidity range $|\eta|<1$. The directed flow in Cu+Au collisions is compared with the rapidity-odd and, for the first time, the rapidity-even components of charged particle directed flow in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 200~GeV. The slope of the directed flow pseudorapidity dependence in Cu+Au collisions is found to be similar to that in Au+Au collisions, with the intercept shifted toward positive $\eta$ values, i.e., the Cu-going direction. The mean transverse momentum projected onto the spectator plane, $\langle p_x\rangle$, in Cu+Au collision also exhibits approximately linear dependence on $\eta$ with the intercept at about $\eta\approx-0.4$, closer to the rapidity of the Cu+Au system center-of-mass. The observed dependencies find natural explanation in a picture of the directed flow originating partly due the "tilted source" and partly due to the rapidity dependent asymmetry in the initial density distribution. Charge-dependence of the $\langle p_x\rangle$ was also observed in Cu+Au collisions, indicating an effect of the initial electric field created by charge difference of the spectator protons in two colliding nuclei. The rapidity-even component of directed flow in Au+Au collisions is close to that in Pb+Pb collisions at $\sqrt{s_{_{NN}}}$ = 2.76 TeV, indicating a similar magnitude of dipole-like fluctuations in the initial-state density distribution. Higher harmonic flow in Cu+Au collisions exhibits similar trends to those observed in Au+Au and Pb+Pb collisions and is qualitatively reproduced by a viscous hydrodynamic model and a multi-phase transport model. For all harmonics with $n\ge2$ we observe an approximate scaling of $v_n$ with the number of constituent quarks.

Collisions between prolate uranium nuclei are used to study how particle production and azimuthal anisotropies depend on initial geometry in heavy-ion collisions. We report the two- and four-particle cumulants, $v_2\{2\}$ and $v_2\{4\}$, for charged hadrons from U+U collisions at $\sqrt{s_{\rm NN}}$ = 193 GeV and Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV. Nearly fully overlapping collisions are selected based on the amount of energy deposited by spectators in the STAR Zero Degree Calorimeters (ZDCs). Within this sample, the observed dependence of $v_2\{2\}$ on multiplicity demonstrates that ZDC information combined with multiplicity can preferentially select different overlap configurations in U+U collisions. An initial-state model with gluon saturation describes the slope of $v_2\{2\}$ as a function of multiplicity in central collisions better than one based on Glauber with a two-component multiplicity model.

We have measured the azimuthal anisotropy of pi0's for 1 < pT < 18 GeV/c for Au+Au collisions at sqrt s_NN = 200 GeV. The observed anisotropy shows a gradual decrease in 3 < pT < 7 - 10 GeV/c, but remains positive beyond 10 GeV/c. The magnitude of this anisotropy is under-predicted, 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 initial-geometry fluctuations is insufficient to account for this discrepancy. Calculations which implement a path length dependence steeper than what is implied by current pQCD energy-loss models, show reasonable agreement with the data.

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.

Yields, correlation shapes, and mean transverse momenta \pt{} of charged particles associated with intermediate to high-\pt{} trigger particles ($2.5 < \pt < 10$ \GeVc) in d+Au and Au+Au collisions at $\snn=200$ GeV are presented. For associated particles at higher $\pt \gtrsim 2.5$ \GeVc, narrow correlation peaks are seen in d+Au and Au+Au, indicating that the main production mechanism is jet fragmentation. At lower associated particle $\pt < 2$ \GeVc, a large enhancement of the near- ($\dphi \sim 0$) and away-side ($\dphi \sim \pi$) associated yields is found, together with a strong broadening of the away-side azimuthal distributions in Au+Au collisions compared to d+Au measurements, suggesting that other particle production mechanisms play a role. This is further supported by the observed significant softening of the away-side associated particle yield distribution at $\dphi \sim \pi$ in central Au+Au collisions.