Transverse momentum spectra of $\pi^{\pm}$, $p$ and $\bar{p}$ up to 12 GeV/c at mid-rapidity in centrality selected Au+Au collisions at $\sqrt{s_{_{NN}}} = 200$ GeV are presented. In central Au+Au collisions, both $\pi^{\pm}$ and $p(\bar{p})$ show significant suppression with respect to binary scaling at $p_T > $ 4 GeV/c. Protons and anti-protons are less suppressed than $\pi^{\pm}$, in the range 1.5 $< p_{T} <$6 GeV/c. The $\pi^-/\pi^+$ and $\bar{p}/p$ ratios show at most a weak $p_T$ dependence and no significant centrality dependence. The $p/\pi$ ratios in central Au+Au collisions approach the values in p+p and d+Au collisions at $p_T >$ 5 GeV/c. The results at high $p_T$ indicate that the partonic sources of $\pi^{\pm}$, $p$ and $\bar{p}$ have similar energy loss when traversing the nuclear medium.

The dependence of $\mathrm{f}_{0}$(980) production on the final-state charged-particle multiplicity in p$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV is reported. The production of $\mathrm{f}_{0}$(980) is measured with the ALICE detector via the $\mathrm{f}_0 (980) \rightarrow \pi^{+}\pi^{-}$ decay channel in a midrapidity region of $-0.5<y<0$. Particle yield ratios of $\mathrm{f}_{0}$(980) to $\pi$ and $\mathrm{K}^{*}$(892)$^{0}$ are found to be decreasing with increasing charged-particle multiplicity. The magnitude of the suppression of the $\mathrm{f}_{0}$(980)/$\pi$ and $\mathrm{f}_{0}$(980)/$\mathrm{K}^{*}$(892)$^{0}$ yield ratios is found to be dependent on the transverse momentum $p_{\mathrm{T}}$, suggesting different mechanisms responsible for the measured effects. Furthermore, the nuclear modification factor $Q_{\mathrm{pPb}}$ of $\mathrm{f}_{0}$(980) is measured in various multiplicity ranges. The $Q_{\mathrm{pPb}}$ shows a strong suppression of the $\mathrm{f}_{0}$(980) production in the $p_{\mathrm{T}}$ region up to about 4 GeV/$c$. The results on the particle yield ratios and $Q_{\mathrm{pPb}}$ for $\mathrm{f}_{0}$(980) may help to understand the late hadronic phase in p$-$Pb collisions and the nature of the internal structure of $\mathrm{f}_{0}$(980) particle.

We present measurements of elliptic flow ($v_2$) of electrons from the decays of heavy-flavor hadrons ($e_{HF}$) by the STAR experiment. For Au+Au collisions at $\sqrt{s_{\rm NN}} = $ 200 GeV we report $v_2$, for transverse momentum ($p_T$) between 0.2 and 7 GeV/c using three methods: the event plane method ($v_{2}${EP}), two-particle correlations ($v_2${2}), and four-particle correlations ($v_2${4}). For Au+Au collisions at $\sqrt{s_{\rm NN}}$ = 62.4 and 39 GeV we report $v_2${2} for $p_T< 2$ GeV/c. $v_2${2} and $v_2${4} are non-zero at low and intermediate $p_T$ at 200 GeV, and $v_2${2} is consistent with zero at low $p_T$ at other energies. The $v_2${2} at the two lower beam energies is systematically lower than at $\sqrt{s_{\rm NN}} = $ 200 GeV for $p_T < 1$ GeV/c. This difference may suggest that charm quarks interact less strongly with the surrounding nuclear matter at those two lower energies compared to $\sqrt{s_{\rm NN}} = 200$ GeV.

The energy deposited at very forward rapidities (very forward energy) is a powerful tool for characterising proton fragmentation in pp and p$-$Pb collisions. The correlation of very forward energy with particle production at midrapidity provides direct insights into the initial stages and the subsequent evolution of the collision. Furthermore, the correlation with the production of particles with large transverse momenta at midrapidity provides information complementary to the measurements of the underlying event, which are usually interpreted in the framework of models implementing centrality-dependent multiple parton interactions. Results about very forward energy, measured by the ALICE zero degree calorimeters (ZDCs), and its dependence on the activity measured at midrapidity in pp collisions at $\sqrt{s}=13$ TeV and in p$-$Pb collisions at $\sqrt{s_{\rm{NN}}}=8.16$ TeV are discussed. The measurements performed in pp collisions are compared with the expectations of three hadronic interaction event generators: PYTHIA 6 (Perugia 2011 tune), PYTHIA 8 (Monash tune), and EPOS LHC. These results provide new constraints on the validity of models in describing the beam remnants at very forward rapidities, where perturbative QCD cannot be used.

We present measurements of $e^+e^-$ production at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 200 GeV. The invariant yield is studied within the PHENIX detector acceptance over a wide range of mass ($m_{ee} <$ 5 GeV/$c^2$) and pair transverse momentum ($p_T$ $<$ 5 GeV/$c$), for minimum bias and for five centrality classes. The \ee yield is compared to the expectations from known sources. In the low-mass region ($m_{ee}=0.30$--0.76 GeV/$c^2$) there is an enhancement that increases with centrality and is distributed over the entire pair \pt range measured. It is significantly smaller than previously reported by the PHENIX experiment and amounts to $2.3\pm0.4({\rm stat})\pm0.4({\rm syst})\pm0.2^{\rm model}$ or to $1.7\pm0.3({\rm stat})\pm0.3({\rm syst})\pm0.2^{\rm model}$ for minimum bias collisions when the open-heavy-flavor contribution is calculated with {\sc pythia} or {\sc mc@nlo}, respectively. The inclusive mass and $p_T$ distributions as well as the centrality dependence are well reproduced by model calculations where the enhancement mainly originates from the melting of the $\rho$ meson resonance as the system approaches chiral symmetry restoration. In the intermediate-mass region ($m_{ee}$ = 1.2--2.8 GeV/$c^2$), the data hint at a significant contribution in addition to the yield from the semileptonic decays of heavy-flavor mesons.

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.

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.

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.

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.

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.