A systematic study is presented for centrality, transverse momentum ($p_T$) and pseudorapidity ($\eta$) dependence of the inclusive charged hadron elliptic flow ($v_2$) at midrapidity($|\eta| < 1.0$) in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 19.6, 27 and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and 4-particle cumulants ($v_2{4}$), are presented in order to investigate non-flow correlations and $v_2$ fluctuations. We observe that the difference between $v_2{2}$ and $v_2{4}$ is smaller at the lower collision energies. Values of $v_2$, scaled by the initial coordinate space eccentricity, $v_{2}/\varepsilon$, as a function of $p_T$ are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider ($\sqrt{s_{NN}}$ = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV). The $v_2(p_T)$ values for fixed $p_T$ rise with increasing collision energy within the $p_T$ range studied ($< 2 {\rm GeV}/c$). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of $v_{2}(p_{T})$. We also compare the $v_2$ results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at Beam Energy Scan (BES) energies are discussed.

We report first results from an analysis based on a new multi-hadron correlation technique, exploring jet-medium interactions and di-jet surface emission bias at RHIC. Pairs of back-to-back high transverse momentum hadrons are used for triggers to study associated hadron distributions. In contrast with two- and three-particle correlations with a single trigger with similar kinematic selections, the associated hadron distribution of both trigger sides reveals no modification in either relative pseudo-rapidity or relative azimuthal angle from d+Au to central Au+Au collisions. We determine associated hadron yields and spectra as well as production rates for such correlated back-to-back triggers to gain additional insights on medium properties.

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.

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

Mid-rapidity open charm spectra from direct reconstruction of $D^{0}$($\bar{D^0}$)$\to K^{\mp}\pi^{\pm}$ in d+Au collisions and indirect electron/positron measurements via charm semileptonic decays in p+p and d+Au collisions at \srt = 200 GeV are reported. The $D^{0}$($\bar{D^0}$) spectrum covers a transverse momentum ($p_T$) range of 0.1 $<p_T<$ 3 \GeVc whereas the electron spectra cover a range of 1 $<p_T<$ 4 GeV/$c$. The electron spectra show approximate binary collision scaling between p+p and d+Au collisions. From these two independent analyses, the differential cross section per nucleon-nucleon binary interaction at mid-rapidity for open charm production from d+Au collisions at RHIC is $d\sigma^{NN}_{c\bar{c}}/dy$=0.30$\pm$0.04 (stat.)$\pm$0.09(syst.) mb. The results are compared to theoretical calculations. Implications for charmoniumm results in A+A collisions are discussed.

We present $\Lambda\Lambda$ correlation measurements in heavy-ion collisions for Au+Au collisions at $\sqrt{s_{NN}}= 200$ GeV using the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). The Lednick\'{y}-Lyuboshitz analytical model has been used to fit the data to obtain a source size, a scattering length and an effective range. Implications of the measurement of the $\Lambda\Lambda$ correlation function and interaction parameters for di-hyperon searches are discussed.

We report the measurement of charged $D^*$ mesons in inclusive jets produced in proton-proton collisions at a center of mass energy $\sqrt{s}$ = 200 GeV with the STAR experiment at RHIC. For $D^{*}$ mesons with fractional momenta $0.2 < z < 0.5$ in inclusive jets with 11.5 GeV mean transverse energy, the production rate is found to be $N(D^{*+}+D^{*-})/N(\mathrm{jet}) = 0.015 \pm 0.008 (\mathrm{stat}) \pm 0.007 (\mathrm{sys})$. This rate is consistent with perturbative QCD evaluation of gluon splitting into a pair of charm quarks and subsequent hadronization.

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.

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.