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.

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.

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 measurements of two-particle angular correlations between high-transverse-momentum ($2<p_T<11$ GeV/$c$) $\pi^0$ observed at midrapidity ($|\eta|<0.35$) and particles produced either at forward ($3.1<\eta<3.9$) or backward ($-3.7<\eta<-3.1$) rapidity in $d$$+$Au and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The azimuthal angle correlations for particle pairs with this large rapidity gap in the Au-going direction exhibit a ridge-like structure that persists up to $p_T{\approx}6$ GeV/$c$ and which strongly depends on collision centrality, which is a similar characteristic to the hydrodynamical particle flow in A+A collisions. The ridge-like structure is absent in the $d$-going direction as well as in $p$$+$$p$ collisions, in the transverse-momentum range studied. The results indicate that the ridge-like structure is shifted in the Au-going direction toward more central collisions, similar to the charged-particle pseudorapidity distributions.

The PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC) has performed systematic measurements of phi meson production in the K+K- decay channel at midrapidity in p+p, d+Au, Cu+Cu and Au+Au collisions at sqrt(S_NN)=200 GeV. Results are presented on the phi invariant yield and the nuclear modification factor R_AA for Au+Au and Cu+Cu, and R_dA for d+Au collisions, studied as a function of transverse momentum (1&lt;p_T&lt;7 GeV/c) and centrality. In central and mid-central Au+Au collisions, the R_AA of phi exhibits a suppression relative to expectations from binary scaled p+p results. The amount of suppression is smaller than that of the neutral pion and the eta meson in the intermediate p_T range (2--5 GeV/c); whereas at higher p_T the phi, pi^0, and eta show similar suppression. The baryon (protons and anti-protons) excess observed in central Au+Au collisions at intermediate p_T is not observed for the phi meson despite the similar mass of the proton and the phi. This suggests that the excess is linked to the number of constituent quarks rather than the hadron mass. The difference gradually disappears with decreasing centrality and for peripheral collisions the R_AA values for both particles are consistent with binary scaling. Cu+Cu collisions show the same yield and suppression as Au+Au collisions for the same number of N_part. The R_dA of phi shows no evidence for cold nuclear effects within uncertainties.

For Au + Au collisions at 200 GeV we measure neutral pion production with good statistics for transverse momentum, p_T, up to 20 GeV/c. A fivefold suppression is found, which is essentially constant for 5 < p_T < 20 GeV/c. Experimental uncertainties are small enough to constrain any model-dependent parameterization for the transport coefficient of the medium, e.g. \mean(q^hat) in the parton quenching model. The spectral shape is similar for all collision classes, and the suppression does not saturate in Au+Au collisions/ instead, it increases proportional to the number of participating nucleons, as N_part^2/3.

The production of low mass e+e- pairs for m_{e+e-} < 300 MeV/c^2 and 1 < p_T <5 GeV/c is measured in p+p and Au+Au collisions at sqrt(s_NN)=200 GeV. Enhanced yield above hadronic sources is observed. Treating the excess as internal conversions, the invariant yield of direct photons is deduced. In central Au+Au collisions, the excess of direct photon yield over p+p is exponential in transverse momentum, with inverse slope T = 221 +/- 19 (stat) +/- 19 (syst) MeV. Hydrodynamical models with initial temperatures ranging from 300--600 MeV at times of ~ 0.6 - 0.15 fm/c after the collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to quark gluon plasma at ~ 170 MeV.

We present measurements of the transverse-momentum dependence of elliptic flow $v_2$ for identified pions and (anti)protons at midrapidity ($|\eta|<0.35$), in 0%--5% central $p$$+$Au and $^3$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. When taken together with previously published measurements in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, the results cover a broad range of small-collision-system multiplicities and intrinsic initial geometries. We observe a clear mass-dependent splitting of $v_2(p_{T})$ in $d$$+$Au and $^3$He$+$Au collisions, just as in large nucleus-nucleus ($A$$+$$A$) collisions, and a smaller splitting in $p$$+$Au collisions. Both hydrodynamic and transport model calculations successfully describe the data at low $p_T$ ($< 1.5$ GeV/$c$), but fail to describe various features at higher $p_T$. In all systems, the $v_2$ values follow an approximate quark-number scaling as a function of the hadron transverse kinetic energy per constituent quark($KE_T/n_q$), which was also seen previously in $A$$+$$A$ collisions.

We report the first measurement of the fraction of $J/\psi$ mesons coming from $B$-meson decay ($F_{B{\rightarrow}J/\psi}$) in $p$+$p$ collisions at $\sqrt{s}=$ 510 GeV. The measurement is performed using the forward silicon vertex detector and central vertex detector at PHENIX, which provide precise tracking and distance-of-closest-approach determinations, enabling the statistical separation of $J/\psi$ due to $B$-meson decays from prompt $J/\psi$. The measured value of $F_{B{\rightarrow}J/\psi}$ is 8.1\%$\pm$2.3\% (stat)$\pm$1.9\% (syst) for $J/\psi$ with transverse momenta $0<p_T<5$ GeV/$c$ and rapidity $1.2<|y|<2.2$. The measured fraction $F_{B{\rightarrow}J/\psi}$ at PHENIX is compared to values measured by other experiments at higher center of mass energies and to fixed-order-next-to-leading-logarithm and color-evaporation-model predictions. The $b\bar{b}$ cross section per unit rapidity ($d\sigma/dy(pp{\rightarrow}b\bar{b})$) extracted from the obtained $F_{B{\rightarrow}J/\psi}$ and the PHENIX inclusive $J/\psi$ cross section measured at 200 GeV scaled with color-evaporation-model calculations, at the mean $B$ hadron rapidity $y={\pm}1.7$ in 510 GeV $p$$+$$p$ collisions, is $3.63^{+1.92}_{-1.70}\mu$b, and it is consistent with the fixed-order-next-to-leading-logarithm calculations.

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section, mean transverse momentum, mean transverse momentum squared of inclusive $J/\psi$ and cross-section ratio of $\psi(2S)$ to $J/\psi$ at forward rapidity in \pp collisions at \sqrts = 510 GeV via the dimuon decay channel. Comparison is made to inclusive $J/\psi$ cross sections measured at \sqrts = 200 GeV and 2.76--13 TeV. The result is also compared to leading-order nonrelativistic QCD calculations coupled to a color-glass-condensate description of the low-$x$ gluons in the proton at low transverse momentum ($p_T$) and to next-to-leading order nonrelativistic QCD calculations for the rest of the $p_T$ range. These calculations overestimate the data at low $p_T$. While consistent with the data within uncertainties above $\approx3$ GeV/$c$, the calculations are systematically below the data. The total cross section times the branching ratio is BR $d\sigma^{J/\psi}_{pp}/dy (1.2<|y|<2.2, 0<p_T<10~\mbox{GeV/$c$}) =$ 54.3 $\pm$ 0.5 (stat) $\pm$ 5.5 (syst) nb.