The production of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{K}^{0}_{S}$, $\rm{K}^{*}(892)^{0}$, $\rm{p}$, $\phi(1020)$, $\Lambda$, $\Xi^{-}$, $\Omega^{-}$, and their antiparticles was measured in inelastic proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV at midrapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm{T}}$) using the ALICE detector at the CERN LHC. Furthermore, the single-particle $p_{\rm{T}}$ distributions of $\rm{K}^{0}_{S}$, $\Lambda$, and $\overline{\Lambda}$ in inelastic pp collisions at $\sqrt{s}$ = 7 TeV are reported here for the first time. The $p_{\rm{T}}$ distributions are studied at midrapidity within the transverse momentum range $0\leq p_{\rm{T}}\leq20$ GeV/$c$, depending on the particle species. The $p_{\rm{T}}$ spectra, integrated yields, and particle yield ratios are discussed as a function of collision energy and compared with measurements at lower $\sqrt{s}$ and with results from various general-purpose QCD-inspired Monte Carlo models. A hardening of the spectra at high $p_{\rm{T}}$ with increasing collision energy is observed, which is similar for all particle species under study. The transverse mass and $x_{\rm{T}}\equiv2p_{\rm{T}}/\sqrt{s}$ scaling properties of hadron production are also studied. As the collision energy increases from $\sqrt{s}$ = 7 to 13 TeV, the yields of non- and single-strange hadrons normalized to the pion yields remain approximately constant as a function of $\sqrt{s}$, while ratios for multi-strange hadrons indicate enhancements. The $p_{\rm{T}}$-differential cross sections of $\pi^{\pm}$, $\rm{K}^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) are compared with next-to-leading order perturbative QCD calculations, which are found to overestimate the cross sections for $\pi^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) at high $p_{\rm{T}}$.

The multiplicity and pseudorapidity distributions of inclusive photons have been measured at forward rapidities ($2.3 < \eta < 3.9$) in proton-proton collisions at three center-of-mass energies, $\sqrt{s}=0.9$, 2.76 and 7 TeV using the ALICE detector. It is observed that the increase in the average photon multiplicity as a function of beam energy is compatible with both a logarithmic and a power-law dependence. The relative increase in average photon multiplicity produced in inelastic pp collisions at 2.76 and 7 TeV center-of-mass energies with respect to 0.9 TeV are 37.2% $\pm$ 0.3% (stat) $\pm$ 8.8% (sys) and 61.2% $\pm$ 0.3% (stat) $\pm$ 7.6% (sys), respectively. The photon multiplicity distributions for all center-of-mass energies are well described by negative binomial distributions. The multiplicity distributions are also presented in terms of KNO variables. The results are compared to model predictions, which are found in general to underestimate the data at large photon multiplicities, in particular at the highest center-of-mass energy. Limiting fragmentation behavior of photons has been explored with the data, but is not observed in the measured pseudorapidity range.

Correlations between mean transverse momentum $[p_{\rm T}]$ and anisotropic flow coefficients $v_{\rm 2}$ or $v_{\rm 3}$ are measured as a function of centrality in Pb$-$Pb and Xe$-$Xe collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and 5.44 TeV, respectively, with ALICE. In addition, the recently proposed higher-order correlation between $[p_{\rm T}]$, $v_{\rm 2}$, and $v_{\rm 3}$ is measured for the first time, which shows an anticorrelation for the presented centrality ranges. These measurements are compared with hydrodynamic calculations using IP-Glasma and $\rm T_{R}ENTo$ initial-state shapes, the former based on the Color Glass Condensate effective theory with gluon saturation, and the latter a parameterized model with nucleons as the relevant degrees of freedom. The data are better described by the IP-Glasma rather than the $\rm T_{R}ENTo$ based calculations. In particular, Trajectum and JETSCAPE predictions, both based on the $\rm T_{R}ENTo$ initial state model but with different parameter settings, fail to describe the measurements. As the correlations between $[p_{\rm T}]$ and $v_{\rm n}$ are mainly driven by the correlations of the size and the shape of the system in the initial state, these new studies pave a novel way to characterize the initial state and help pin down the uncertainty of the extracted properties of the quark$-$gluon plasma recreated in relativistic heavy-ion collisions.

The first measurements of the production of muons and electrons from heavy-flavour hadron decays in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV, using the ALICE detector at the LHC, are reported. The measurement of the nuclear modification factor $R_{\rm AA}$ is performed as a function of transverse momentum $p_{\rm T}$ in several centrality classes at forward rapidity ($2.5 < y <4$) and midrapidity ($\vert y \vert < 0.8$) for muons and electrons from heavy-flavour hadron decays, respectively. A suppression by a factor up to about 2.5 compared to the binary-scaled pp reference is observed in central collisions at both central and forward rapidities. The $R_{\rm AA}$ of muons from heavy-flavour hadron decays is compared to previous measurements in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. When the nuclear modification factors are compared in the centrality classes 0-10% for Xe-Xe collisions and 10-20% for Pb-Pb collisions, which have similar charged-particle multiplicity density, a similar suppression, with $R_{\rm AA} \sim 0.4$ in the $p_{\rm T}$ interval $4 < p_{\rm T} < 8$ GeV/$c$, is observed. The comparison of the measured $R_{\rm AA}$ values in the two collision systems brings new insights on the properties of the quark-gluon plasma by investigating the system-size and path-length dependence of medium-induced parton energy loss. The results of muons and electrons from heavy-flavour hadron decays provide new constraints to model calculations.

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons $v_2${$\Psi_{\rm SP}$} in collisions of Pb ions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$=2.76 TeV and Xe ions at $\sqrt{s_{\rm NN}}$=5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum. The ratio between $v_2${$\Psi_{\rm SP}$} and the elliptic flow coefficient relative to the participant plane $v_2$4, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the TRENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models of initial state fluctuations. A significant transverse momentum dependence of the ratio $v_2${$\Psi_{\rm SP}$}/$v_2${4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of $v_2${$\Psi_{\rm SP}$} and $v_2${4} to the corresponding initial state eccentricities for Xe-Xe and Pb-Pb collisions at similar initial entropy density show a difference of $(7.0 \pm 0.9)$% with an additional variation of +1.8% when including RHIC data in the TRENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark-gluon plasma.

Measurements of inclusive and direct photon production at mid-rapidity in pp collisions at $\sqrt{s}=2.76$ and 8 TeV are presented by the ALICE experiment at the LHC. The results are reported in transverse momentum ranges of $0.4<p_{T}<10$ GeV/$c$ and $0.3<p_{T}<16$ GeV/$c$, respectively. Photons are detected with the electromagnetic calorimeter~(EMCal) and via reconstruction of e$^+$e$^-$ pairs from conversions in the ALICE detector material using the central tracking system. For the final measurement of the inclusive photon spectra the results are combined in the overlapping $p_{T}$ interval of both methods. Direct photon spectra, or their upper limits at 90% C.L. are extracted using the direct photon excess ratio $R_{\gamma}$, which quantifies the ratio of inclusive photons over decay photons generated with a decay-photon simulation. An additional hybrid method, combining photons reconstructed from conversions with those identified in the EMCal, is used for the combination of the direct photon excess ratio $R_{\gamma}$, as well as the extraction of direct photon spectra or their upper limits. While no significant signal of direct photons is seen over the full $p_{T}$ range, $R_{\gamma}$ for $p_{T}>7$ GeV/$c$ is at least one $\sigma$ above unity and consistent with expectations from next-to-leading order pQCD calculations.

Multiplicity ($N_{\rm ch}$) distributions and transverse momentum ($p_{\rm T}$) spectra of inclusive primary charged particles in the kinematic range of $|\eta| < 0.8$ and 0.15 GeV/$c$$< p_{T} <$ 10 GeV/$c$ are reported for pp, p-Pb, Xe-Xe and Pb-Pb collisions at centre-of-mass energies per nucleon pair ranging from $\sqrt{s_{\rm NN}} = 2.76$ TeV up to $13$ TeV. A sequential two-dimensional unfolding procedure is used to extract the correlation between the transverse momentum of primary charged particles and the charged-particle multiplicity of the corresponding collision. This correlation sharply characterises important features of the final state of a collision and, therefore, can be used as a stringent test of theoretical models. The multiplicity distributions as well as the mean and standard deviation derived from the $p_{\rm T}$ spectra are compared to state-of-the-art model predictions. Providing these fundamental observables of bulk particle production consistently across a wide range of collision energies and system sizes can serve as an important input for tuning Monte Carlo event generators.

We present measurements of Underlying Event observables in pp collisions at $\sqrt{s}$ = 0.9 and 7 TeV. The analysis is performed as a function of the highest charged-particle transverse momentum $p_{\rm T, LT}$ in the event. Different regions are defined with respect to the azimuthal direction of the leading (highest transverse momentum) track: Toward, Transverse and Away. The Toward and Away regions collect the fragmentation products of the hardest partonic interaction. The Transverse region is expected to be most sensitive to the Underlying Event activity. The study is performed with charged particles above three different $p_{\rm T}$ thresholds: 0.15, 0.5 and 1.0 GeV/$c$. In the Transverse region we observe an increase in the multiplicity of a factor 2-3 between the lower and higher collision energies, depending on the track $p){\rm T}$ threshold considered. Data are compared to Pythia 6.4, Pythia 8.1 and Phojet. On average, all models considered underestimate the multiplicity and summed $p_{\rm T}$ in the Transverse region by about 10-30%.

The first measurement of the production of pions, kaons, (anti-)protons and $\phi$ mesons at midrapidity in Xe-Xe collisions at $\sqrt{s_{\rm NN}} = 5.44$ TeV is presented. Transverse momentum ($p_{\rm T}$) spectra and $p_{\rm T}$-integrated yields are extracted in several centrality intervals bridging from p-Pb to mid-central Pb-Pb collisions in terms of final-state multiplicity. The study of Xe-Xe and Pb-Pb collisions allows systems at similar charged-particle multiplicities but with different initial geometrical eccentricities to be investigated. A detailed comparison of the spectral shapes in the two systems reveals an opposite behaviour for radial and elliptic flow. In particular, this study shows that the radial flow does not depend on the colliding system when compared at similar charged-particle multiplicity. In terms of hadron chemistry, the previously observed smooth evolution of particle ratios with multiplicity from small to large collision systems is also found to hold in Xe-Xe. In addition, our results confirm that two remarkable features of particle production at LHC energies are also valid in the collision of medium-sized nuclei: the lower proton-to-pion ratio with respect to the thermal model expectations and the increase of the $\phi$-to-pion ratio with increasing final-state multiplicity.

The pseudorapidity dependence of elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of charged particles measured in Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}=5.02$ TeV and in Xe-Xe collisions at $\sqrt{s_{\rm NN}}=5.44$ TeV with ALICE at the LHC are presented. The measurements are performed in the pseudorapidity range $-3.5 < \eta < 5$ for various centrality intervals using two- and multi-particle cumulants with the subevent method. The flow probability density function (p.d.f.) is studied with the ratio of flow coefficient $v_2$ calculated with four- and two-particle cumulant, and suggests that the variance of flow p.d.f. is independent of pseudorapidity. The decorrelation of the flow vector in the longitudinal direction is probed using two-particle correlations. The results measured with respect to different reference regions in pseudorapidity exhibit differences, argued to be a result of saturating decorrelation effect above a certain pseudorapidity separation, in contrast to previous publications which assign this observation to non-flow effects. The results are compared to $3+1$ dimensional hydrodynamic and the AMPT transport model calculations. Neither of the models is able to simultaneously describe the pseudorapidity dependence of measurements of anisotropic flow and its fluctuations. The results presented in this work highlight shortcomings in our current understanding of initial conditions and subsequent system expansion in the longitudinal direction. Therefore, they provide input for its improvement.