The system created in non-central relativistic nucleus-nucleus collisions possesses large orbital angular momentum. Due to spin-orbit coupling, particles produced in such a system could become globally polarized along the direction of the system angular momentum. We present the results of Lambda and anti-Lambda hyperon global polarization measurements in Au+Au collisions at sqrt{s_NN}=62.4 GeV and 200 GeV performed with the STAR detector at RHIC. The observed global polarization of Lambda and anti-Lambda hyperons in the STAR acceptance is consistent with zero within the precision of the measurements. The obtained upper limit, |P_{Lambda,anti-Lambda}| <= 0.02, is compared to the theoretical values discussed recently in the literature.

We present a systematic analysis of two-pion interferometry in Au+Au collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 GeV and Cu+Cu collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 and 200 GeV using the STAR detector at RHIC. The multiplicity and transverse momentum dependences of the extracted femtoscopic radii are studied. The scaling of the apparent freeze-out volume with charged particle multiplicity is studied for the RHIC energy domain. The multiplicity scaling of the measured radii is found to be independent of colliding system and collision energy.

We present results for the measurement of $\phi$ meson production via its charged kaon decay channel $\phi \to K^+K^-$ in Au+Au collisions at $\sqrt{s_{_{NN}}}=62.4$, 130, 200 GeV, and in p+p and d+Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV from the STAR experiment at RHIC. The mid-rapidity ($|y|<0.5$) $\phi$ meson spectra in central Au+Au collisions are found to be well described by a single exponential distribution. On the other hand, the spectra from p+p, d+Au and peripheral Au+Au collisions show power-law tails at intermediate and high transverse momenta ($p_{T}$) and are described better by Levy distributions. The constant $\phi/K^-$ yield ratio vs. beam species, collision centrality and colliding energy is in contradiction with expectations from models having kaon coalescence as the dominant mechanism for $\phi$ production at RHIC. The $\Omega/\phi$ yield ratio as a function of $p_{T}$ is consistent with a model based on the recombination of thermal $s$ quarks up to $p_{T}\sim 4$ GeV/c, but disagrees at higher transverse momenta. The measured nuclear modification factor, $R_{dAu}$, for the $\phi$ meson increases above unity at intermediate $p_{T}$, similar to that for pions and protons, while $R_{AA}$ is suppressed due to jet quenching in central Au+Au collisions. Number of constituent quark scaling of both $R_{cp}$ and $v_{2}$ for the $\phi$ meson with respect to other hadrons in Au+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV at intermediate $p_{T}$ is observed. These observations support quark coalescence as being the dominant mechanism of hadronization in the intermediate $p_{T}$ region at RHIC.

We study the beam-energy and system-size dependence of \phi meson production (using the hadronic decay mode \phi -- K+K-) by comparing the new results from Cu+Cu collisions and previously reported Au+Au collisions at \sqrt{s_NN} = 62.4 and 200 GeV measured in the STAR experiment at RHIC. Data presented are from mid-rapidity (|y|&lt;0.5) for 0.4 &lt; pT &lt; 5 GeV/c. At a given beam energy, the transverse momentum distributions for \phi mesons are observed to be similar in yield and shape for Cu+Cu and Au+Au colliding systems with similar average numbers of participating nucleons. The \phi meson yields in nucleus-nucleus collisions, normalised by the average number of participating nucleons, are found to be enhanced relative to those from p+p collisions with a different trend compared to strange baryons. The enhancement for \phi mesons is observed to be higher at \sqrt{s_NN} = 200 GeV compared to 62.4 GeV. These observations for the produced \phi(s\bar{s}) mesons clearly suggest that, at these collision energies, the source of enhancement of strange hadrons is related to the formation of a dense partonic medium in high energy nucleus-nucleus collisions and cannot be alone due to canonical suppression of their production in smaller systems.

Parity-odd domains, corresponding to non-trivial topological solutions of the QCD vacuum, might be created during relativistic heavy-ion collisions. These domains are predicted to lead to charge separation of quarks along the orbital momentum of the system created in non-central collisions. To study this effect, we investigate a three particle mixed harmonics azimuthal correlator which is a \P-even observable, but directly sensitive to the charge separation effect. We report measurements of this observable using the STAR detector in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$=200 and 62~GeV. The results are presented as a function of collision centrality, particle separation in rapidity, and particle transverse momentum. A signal consistent with several of the theoretical expectations is detected in all four data sets. We compare our results to the predictions of existing event generators, and discuss in detail possible contributions from other effects that are not related to parity violation.

We present the results of an elliptic flow analysis of Cu+Cu collisions recorded with the STAR detector at 62.4 and 200GeV. Elliptic flow as a function of transverse momentum is reported for different collision centralities for charged hadrons and strangeness containing hadrons $K_{S}^{0}$, $\Lambda$, $\Xi$, $\phi$ in the midrapidity region $|eta|<1.0$. Significant reduction in systematic uncertainty of the measurement due to non-flow effects has been achieved by correlating particles at midrapidity, $|\eta|<1.0$, with those at forward rapidity, $2.5<|\eta|<4.0$. We also present azimuthal correlations in p+p collisions at 200 GeV to help estimating non-flow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au+Au collisions at 200 GeV. We observe that $v_{2}$($p_{T}$) of strange hadrons has similar scaling properties as were first observed in Au+Au collisions, i.e.: (i) at low transverse momenta, $p_T<2GeV/c$, $v_{2}$ scales with transverse kinetic energy, $m_{T}-m$, and (ii) at intermediate $p_T$, $2<p_T<4GeV/c$, it scales with the number of constituent quarks, $n_q$. We have found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of $v_{2}$($p_{T}$) for $K_{S}^{0}$ and $\Lambda$. Eccentricity scaled $v_2$ values, $v_{2}/\epsilon$, are larger in more central collisions, suggesting stronger collective flow develops in more central collisions. The comparison with Au+Au collisions which go further in density shows $v_{2}/\epsilon$ depend on the system size, number of participants $N_{part}$. This indicates that the ideal hydrodynamic limit is not reached in Cu+Cu collisions, presumably because the assumption of thermalization is not attained.

We present STAR measurements of the azimuthal anisotropy parameter $v_2$ for pions, kaons, protons, $\Lambda$, $\bar{\Lambda}$, $\Xi+\bar{\Xi}$, and $\Omega + \bar{\Omega}$, along with $v_4$ for pions, kaons, protons, and $\Lambda + \bar{\Lambda}$ at mid-rapidity for Au+Au collisions at $\sqrt{s_{_{NN}}}=62.4$ and 200 GeV. The $v_2(p_T)$ values for all hadron species at 62.4 GeV are similar to those observed in 130 and 200 GeV collisions. For observed kinematic ranges, $v_2$ values at 62.4, 130, and 200 GeV are as little as 10%--15% larger than those in Pb+Pb collisions at $\sqrt{s_{_{NN}}}=17.3$ GeV. At intermediate transverse momentum ($p_T$ from 1.5--5 GeV/c), the 62.4 GeV $v_2(p_T)$ and $v_4(p_T)$ values are consistent with the quark-number scaling first observed at 200 GeV. A four-particle cumulant analysis is used to assess the non-flow contributions to pions and protons and some indications are found for a smaller non-flow contribution to protons than pions. Baryon $v_2$ is larger than anti-baryon $v_2$ at 62.4 and 200 GeV perhaps indicating either that the initial spatial net-baryon distribution is anisotropic, that the mechanism leading to transport of baryon number from beam- to mid-rapidity enhances $v_2$, or that anti-baryon and baryon annihilation is larger in the in-plane direction.

We study the energy dependence of the transverse momentum (pT) spectra for charged pions, protons and anti-protons for Au+Au collisions at \sqrt{s_NN} = 62.4 and 200 GeV. Data are presented at mid-rapidity (|y| &lt; 0.5) for 0.2 &lt; pT &lt; 12 GeV/c. In the intermediate pT region (2 &lt; pT &lt; 6 GeV/c), the nuclear modification factor is higher at 62.4 GeV than at 200 GeV, while at higher pT (pT >7 GeV/c) the modification is similar for both energies. The p/pi+ and pbar/pi- ratios for central collisions at \sqrt{s_NN} = 62.4 GeV peak at pT ~ 2 GeV/c. In the pT range where recombination is expected to dominate, the p/pi+ ratios at 62.4 GeV are larger than at 200 GeV, while the pbar/pi- ratios are smaller. For pT > 2 GeV/c, the pbar/pi- ratios at the two beam energies are independent of pT and centrality indicating that the dependence of the pbar/pi- ratio on pT does not change between 62.4 and 200 GeV. These findings challenge various models incorporating jet quenching and/or constituent quark coalescence.

We present the directed flow ($v_1$) measured in Au+Au collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV in the mid-pseudorapidity region $|\eta|<1.3$ and in the forward pseudorapidity region $2.5 < |\eta| < 4.0$. The results are obtained using the three-particle cumulant method, the event plane method with mixed harmonics, and for the first time at the Relativistic Heavy Ion Collider (RHIC), the standard method with the event plane reconstructed from spectator neutrons. Results from all three methods are in good agreement. Over the pseudorapidity range studied, charged particle directed flow is in the direction opposite to that of fragmentation neutrons.

We present the scaling properties of Lambda, Xi, Omega and their anti-particles produced at mid-rapidity in Au+Au collisions at RHIC at sqrt(s_NN) = 200 GeV. The yield of multi-strange baryons per participant nucleon increases from peripheral to central collisions more rapidly than the Lambda yield, which appears to correspond to an increasing strange quark density of matter produced. The value of the strange phase space occupancy factor gamma_s, obtained from a thermal model fit to the data, approaches unity for the most central collisions. We also show that the nuclear modification factors, R_CP, of Lambda and Xi are consistent with each other and with that of protons in the transverse momentum range 2.0 < p_T < 5.0 GeV/c. This scaling behaviour is consistent with a scenario of hadron formation from constituent quark degrees of freedom through quark recombination or coalescence.

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.

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

Neutral-pion, pi^0, spectra were measured at midrapidity (|y|<0.35) in Au+Au collisions at sqrt(s_NN) = 39 and 62.4 GeV and compared to earlier measurements at 200 GeV in the 1<p_T<10 GeV/c transverse-momentum (p_T) range. The high-p_T tail is well described by a power law in all cases and the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very similar to that observed in the corresponding p+p-collision spectra. The nuclear-modification factors (R_AA) show significant suppression and a distinct energy dependence at moderate p_T in central collisions. At high p_T, R_AA is similar for 62.4 and 200 GeV at all centralities. Perturbative-quantum-chromodynamics calculations that describe R_AA well at 200 GeV, fail to describe the 39 GeV data, raising the possibility that the relative importance of initial-state effects and soft processes increases at lower energies. A conclusion that the region where hard processes are dominant is reached only at higher p_T, is also supported by the x_T dependence of the x_T-scaling power-law exponent.

We present the first measurement of multiplicity and pseudorapidity distributions of photons in the pseudorapidity region 2.3 $\le$ $\eta$ $\le$ 3.7 for different centralities in Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV. We find that the photon yield in this pseudorapidity range scales with the number of participating nucleons at all collision centralities studied. The pseudorapidity distribution of photons, dominated by neutral pion decays, has been compared to those of identified charged pions, photons, and inclusive charged particles from heavy ion and nucleon-nucleon collisions at various energies. The photon production in the measured pseudorapidity region has been shown to be consistent with the energy and centrality independent limiting fragmentation scenario.

The measurement of direct photons from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=39$ and 62.4 GeV in the transverse-momentum range $0.4<p_T<3$ Gev/$c$ is presented by the PHENIX collaboration at the Relativistic Heavy Ion Collider. A significant direct-photon yield is observed in both collision systems. A universal scaling is observed when the direct-photon $p_T$ spectra for different center-of-mass energies and for different centrality selections at $\sqrt{s_{_{NN}}}=62.4$ GeV is scaled with $(dN_{\rm ch}/d\eta)^{\alpha}$ for $\alpha=1.21{\pm}0.04$. This scaling also holds true for direct-photon spectra from Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV measured earlier by PHENIX, as well as the spectra from Pb$+$Pb at $\sqrt{s_{_{NN}}}=2760$ GeV published by ALICE. The scaling power $\alpha$ seems to be independent of $p_T$, center of mass energy, and collision centrality. The spectra from different collision energies have a similar shape up to $p_T$ of 2 GeV/$c$. The spectra have a local inverse slope $T_{\rm eff}$ increasing with $p_T$ of $0.174\pm0.018$ GeV/$c$ in the range $0.4<p_T<1.3$ GeV/$c$ and increasing to $0.289\pm0.024$ GeV/$c$ for $0.9<p_T<2.1$ GeV/$c$. The observed similarity of low-$p_T$ direct-photon production from $\sqrt{s_{_{NN}}}= 39$ to 2760 GeV suggests a common source of direct photons for the different collision energies and event centrality selections, and suggests a comparable space-time evolution of direct-photon emission.

We present measurements of the J/psi invariant yields in sqrt(s_NN)=39 and 62.4 GeV Au+Au collisions at forward rapidity (1.2<|y|<2.2). Invariant yields are presented as a function of both collision centrality and transverse momentum. Nuclear modifications are obtained for central relative to peripheral Au+Au collisions (R_CP) and for various centrality selections in Au+Au relative to scaled p+p cross sections obtained from other measurements (R_AA). The observed suppression patterns at 39 and 62.4 GeV are quite similar to those previously measured at 200 GeV. This similar suppression presents a challenge to theoretical models that contain various competing mechanisms with different energy dependencies, some of which cause suppression and others enhancement.

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.

The jet fragmentation function is measured with direct photon-hadron correlations in p+p and Au+Au collisions at sqrt(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^\gamma 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 jet fragment yield 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 fragment 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.

The PHENIX experiment presents results from the RHIC 2006 run with polarized proton collisions at sqrt(s) = 62.4 GeV for inclusive pi^0 production at mid-rapidity. Unpolarized cross section results are measured for transverse momenta p_T = 0.5 to 7 GeV/c. Next-to-leading order perturbative quantum chromodynamics calculations are compared with the data, and while the calculations are consistent with the measurements, next-to-leading logarithmic corrections improve the agreement. Double helicity asymmetries A_LL are presented for p_T = 1 to 4 GeV/c and probe the higher range of Bjorken_x of the gluon (x_g) with better statistical precision than our previous measurements at sqrt(s)=200 GeV. These measurements are sensitive to the gluon polarization in the proton for 0.06 < x_g < 0.4.

Unpolarized cross sections and double-helicity asymmetries of single-inclusive positive and negative charged hadrons at midrapidity from p+p collisions at sqrt(s)=62.4 GeV are presented. The PHENIX measurements for 1.0 < p_T < 4.5 GeV/c are consistent with perturbative QCD calculations at next-to-leading order in the strong coupling constant, alpha_s. Resummed pQCD calculations including terms with next-to-leading-log accuracy, yielding reduced theoretical uncertainties, also agree with the data. The double-helicity asymmetry, sensitive at leading order to the gluon polarization in a momentum-fraction range of 0.05 ~< x_gluon ~< 0.2, is consistent with recent global parameterizations disfavoring large gluon polarization.

We report the measurement of cumulants ($C_n, n=1\ldots4$) of the net-charge distributions measured within pseudorapidity ($|\eta|<0.35$) in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=7.7-200$ GeV with the PHENIX experiment at the Relativistic Heavy Ion Collider. The ratios of cumulants (e.g. $C_1/C_2$, $C_3/C_1$) of the net-charge distributions, which can be related to volume independent susceptibility ratios, are studied as a function of centrality and energy. These quantities are important to understand the quantum-chromodynamics phase diagram and possible existence of a critical end point. The measured values are very well described by expectation from negative binomial distributions. We do not observe any nonmonotonic behavior in the ratios of the cumulants as a function of collision energy. The measured values of $C_1/C_2 = \mu/\sigma^2$ and $C_3/C_1 = S\sigma^3/\mu$ can be directly compared to lattice quantum-chromodynamics calculations and thus allow extraction of both the chemical freeze-out temperature and the baryon chemical potential at each center-of-mass energy.

Pair correlations between large transverse momentum neutral pion triggers (p_T=4--7 GeV/c) and charged hadron partners (p_T=3--7 GeV/c) in central (0--20%) and midcentral (20--60%) Au+Au collisions are presented as a function of trigger orientation with respect to the reaction plane. The particles are at larger momentum than where jet shape modifications have been observed, and the correlations are sensitive to the energy loss of partons traveling through hot dense matter. An out-of-plane trigger particle produces only 26+/-20% of the away-side pairs that are observed opposite of an in-plane trigger particle. In contrast, near-side jet fragments are consistent with no suppression or dependence on trigger orientation with respect to the reaction plane. These observations are qualitatively consistent with a picture of little near-side parton energy loss either due to surface bias or fluctuations and increased away-side parton energy loss due to a long path through the medium. The away-side suppression as a function of reaction-plane angle is shown to be sensitive to both the energy loss mechanism in and the space-time evolution of heavy-ion collisions.

We present the centrality dependent measurement of multiplicity and pseudorapidity distributions of charged particles and photons in Au + Au collisions at sqrt{s_NN} = 62.4 GeV. The charged particles and photons are measured in the pseudorapidity region 2.9 < eta < 3.9 and 2.3 < eta < 3.7, respectively. We have studied the scaling of particle production with the number of participating nucleons and the number of binary collisions. The photon and charged particle production in the measured pseudorapidity range has been shown to be consistent with energy independent limiting fragmentation behavior. The photons are observed to follow a centrality independent limiting fragmentation behavior while for the charged particles it is centrality dependent. We have carried out a comparative study of the pseudorapidity distributions of positively charged hadrons, negatively charged hadrons, photons, pions, net protons in nucleus--nucleus collisions and pseudorapidity distributions from p+p collisions. From these comparisons we conclude that baryons in the inclusive charged particle distribution are responsible for the observed centrality dependence of limiting fragmentation. The mesons are found to follow an energy independent behavior of limiting fragmentation while the behavior of baryons seems to be energy dependent.

Measurements of the fractional momentum loss ($S_{\rm loss}\equiv{\delta}p_T/p_T$) of high-transverse-momentum-identified hadrons in heavy ion collisions are presented. Using $\pi^0$ in Au$+$Au and Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=62.4$ and 200 GeV measured by the PHENIX experiment at the Relativistic Heavy Ion Collider and and charged hadrons in Pb$+$Pb collisions measured by the ALICE experiment at the Large Hadron Collider, we studied the scaling properties of $S_{\rm loss}$ as a function of a number of variables: the number of participants, $N_{\rm part}$, the number of quark participants, $N_{\rm qp}$, the charged-particle density, $dN_{\rm ch}/d\eta$, and the Bjorken energy density times the equilibration time, $\varepsilon_{\rm Bj}\tau_{0}$. We find that the $p_T$ where $S_{\rm loss}$ has its maximum, varies both with centrality and collision energy. Above the maximum, $S_{\rm loss}$ tends to follow a power-law function with all four scaling variables. The data at $\sqrt{s_{_{NN}}}$=200 GeV and 2.76 TeV, for sufficiently high particle densities, have a common scaling of $S_{\rm loss}$ with $dN_{\rm ch}/d\eta$ and $\varepsilon_{\rm Bj}\tau_{0}$, lending insight on the physics of parton energy loss.

We present the first study of the energy dependence of $p_t$ angular correlations inferred from event-wise mean transverse momentum $<p_{t} >$ fluctuations in heavy ion collisions. We compare our large-acceptance measurements at CM energies $\sqrt{s_{NN}} =$ 19.6, 62.4, 130 and 200 GeV to SPS measurements at 12.3 and 17.3 GeV. $p_t$ angular correlation structure suggests that the principal source of $p_t$ correlations and fluctuations is minijets (minimum-bias parton fragments). We observe a dramatic increase in correlations and fluctuations from SPS to RHIC energies, increasing linearly with $\ln \sqrt{s_{NN}}$ from the onset of observable jet-related $<p_{t}>$ fluctuations near 10 GeV.

Measurements of transverse-single-spin asymmetries ($A_{N}$) in $p$$+$$p$ collisions at $\sqrt{s}=$62.4 and 200 GeV with the PHENIX detector at RHIC are presented. At midrapidity, $A_{N}$ is measured for neutral pion and eta mesons reconstructed from diphoton decay, and at forward rapidities, neutral pions are measured using both diphotons and electromagnetic clusters. The neutral-pion measurement of $A_{N}$ at midrapidity is consistent with zero with uncertainties a factor of 20 smaller than previous publications, which will lead to improved constraints on the gluon Sivers function. At higher rapidities, where the valence quark distributions are probed, the data exhibit sizable asymmetries. In comparison with previous measurements in this kinematic region, the new data extend the kinematic coverage in $\sqrt{s}$ and $p_T$, and it is found that the asymmetries depend only weakly on $\sqrt{s}$. The origin of the forward $A_{N}$ is presently not understood quantitatively. The extended reach to higher $p_T$ probes the transition between transverse momentum dependent effects at low $p_T$ and multi-parton dynamics at high $p_T$.