The ALICE experiment has measured the inclusive J/$\psi$ production in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}} } = 2.76$ TeV down to zero transverse momentum in the rapidity range $2.5 < y < 4$. A suppression of the inclusive J/$\psi$ yield in Pb-Pb is observed with respect to the one measured in pp collisions scaled by the number of binary nucleon-nucleon collisions. The nuclear modification factor, integrated over the 0-80% most central collisions, is $0.545 \pm 0.032 \rm{(stat.)} \pm 0.083 \rm{(syst.)}$ and does not exhibit a significant dependence on the collision centrality. These features appear significantly different from measurements at lower collision energies. Models including J/$\psi$ production from charm quarks in a deconfined partonic phase can describe our data.

The PHENIX experiment has measured $\phi$ meson production in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV using the dimuon and dielectron decay channels. The $\phi$ meson is measured in the forward (backward) $d$-going (Au-going) direction, $1.2<y<2.2$ ($-2.2<y<-1.2$) in the transverse-momentum ($p_T$) range from 1--7 GeV/$c$, and at midrapidity $|y|<0.35$ in the $p_T$ range below 7 GeV/$c$. The $\phi$ meson invariant yields and nuclear-modification factors as a function of $p_T$, rapidity, and centrality are reported. An enhancement of $\phi$ meson production is observed in the Au-going direction, while suppression is seen in the $d$-going direction, and no modification is observed at midrapidity relative to the yield in $p$$+$$p$ collisions scaled by the number of binary collisions. Similar behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating similar cold-nuclear-matter effects.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has measured $\phi$ meson production and its nuclear modification in asymmetric Cu$+$Au heavy-ion collisions at $\sqrt{s_{NN}}=200$ GeV at both forward Cu-going direction ($1.2<y<2.2$) and backward Au-going direction ($-2.2<y<-1.2$), rapidities. The measurements are performed via the dimuon decay channel and reported as a function of the number of participating nucleons, rapidity, and transverse momentum. In the most central events, 0\%--20\% centrality, the $\phi$ meson yield integrated over $1<p_T<5$ GeV/$c$ prefers a smaller value, which means a larger nuclear modification, in the Cu-going direction compared to the Au-going direction. Additionally, the nuclear-modification factor in Cu$+$Au collisions averaged over all centrality is measured to be similar to the previous PHENIX result in $d$$+$Au collisions for these rapidities.

The elliptic, $v_2$, triangular, $v_3$, and quadrangular, $v_4$, azimuthal anisotropic flow coefficients are measured for unidentified charged particles, pions and (anti-)protons in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the Large Hadron Collider. Results obtained with the event plane and four-particle cumulant methods are reported for the pseudo-rapidity range $|\eta|<0.8$ at different collision centralities and as a function of transverse momentum, $p_{\rm T}$, out to $p_{\rm T}=20$ GeV/$c$. The observed non-zero elliptic and triangular flow depends only weakly on transverse momentum for $p_{\rm T}>8$ GeV/$c$. The small $p_{\rm T}$ dependence of the difference between elliptic flow results obtained from the event plane and four-particle cumulant methods suggests a common origin of flow fluctuations up to $p_{\rm T}=8$ GeV/$c$. The magnitude of the (anti-)proton elliptic and triangular flow is larger than that of pions out to at least $p_{\rm T}=8$ GeV/$c$ indicating that the particle type dependence persists out to high $p_{\rm T}$.

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 production of the prompt charmed mesons $D^0$, $D^+$ and $D^{*+}$ relative to the reaction plane was measured in Pb-Pb collisions at a centre-of-mass energy per nucleon-nucleon collision of $\sqrt{s_{\rm NN}} = 2.76$ TeV with the ALICE detector at the LHC. D mesons were reconstructed via their hadronic decays at central rapidity in the transverse momentum ($p_{\rm T}$) interval of 2-16 GeV/$c$. The azimuthal anisotropy is quantified in terms of the second coefficient $v_2$ in a Fourier expansion of the D meson azimuthal distribution, and in terms of the nuclear modification factor $R_{\rm AA}$, measured in the direction of the reaction plane and orthogonal to it. The $v_2$ coefficient was measured with three different methods and in three centrality classes in the interval 0-50%. A positive $v_2$ is observed in mid-central collisions (30-50% centrality class), with an mean value of $0.204_{-0.036}^{+0.099}$ (tot.unc.) in the interval $2 < p_{\rm T} < 6$ GeV/$c$, which decreases towards more central collisions (10-30% and 0-10% classes). The positive $v_2$ is also reflected in the nuclear modification factor, which shows a stronger suppression in the direction orthogonal to the reaction plane for mid-central collisions. The measurements are compared to theoretical calculations of charm quark transport and energy loss in high-density strongly-interacting matter at high temperature. The models that include substantial elastic interactions with an expanding medium provide a good description of the observed anisotropy. However, they are challenged to simultaneously describe the strong suppression of high-$p_{\rm T}$ yield of D mesons in central collisions and their azimuthal anisotropy in non-central collisions.

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.

Charged-pion-interferometry measurements were made with respect to the 2$^{\rm nd}$- and 3$^{\rm rd}$-order event plane for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. A strong azimuthal-angle dependence of the extracted Gaussian-source radii was observed with respect to both the 2$^{\rm nd}$- and 3$^{\rm rd}$-order event planes. The results for the 2$^{\rm nd}$-order dependence indicate that the initial eccentricity is reduced during the medium evolution, but not reversed in the final state, which is consistent with previous results. In contrast, the results for the 3$^{\rm rd}$-order dependence indicate that the initial triangular shape is significantly reduced and potentially reversed by the end of the medium evolution, and that the 3$^{\rm rd}$-order oscillations are largely dominated by the dynamical effects from triangular flow.