Measurements of elliptic ($v_2$) and triangular ($v_3$) flow coefficients of $\pi^{\pm}$, K$^{\pm}$, p+$\rm \overline{p}$, K$^0_{\rm S}$, and $\Lambda + \overline{\Lambda}$ obtained with the scalar product method in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV are presented. The results are obtained in the rapidity range $\left | y \right |<0.5$ and reported as a function of transverse momentum, $p_{\rm T}$, for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for $p_{\rm T}<3$ GeV/$c$, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3< $p_{\rm T}$ <8 GeV/$c$). The magnitude of the baryon $v_{2}$ is larger than that of mesons up to $p_{\rm T}$ = 6 GeV/$c$. The centrality dependence of the shape evolution of the $p_{\rm T}$-differential $v_2$ is studied for the various hadron species. The $v_2$ coefficients of $\pi^{\pm}$, K$^{\pm}$, and p+$\rm \overline{p}$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for $p_{\rm T} <1$ GeV/$c$. A comparison with $v_{\rm n}$ measurements in the corresponding centrality intervals in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV yields an enhanced $v_2$ in central collisions and diminished value in semicentral collisions.

Inclusive J/$\psi$ production is studied in Xe-Xe interactions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}= 5.44$ TeV, using the ALICE detector at the CERN LHC. The J/$\psi$ meson is reconstructed via its decay into a muon pair, in the centre-of-mass rapidity interval $2.5<y<4$ and down to zero transverse momentum. In this Letter, the nuclear modification factors $R_{\rm AA}$ for inclusive J/$\psi$, measured in the centrality range 0-90% as well as in the centrality intervals 0-20% and 20-90% are presented. The $R_{\rm AA}$ values are compared to previously published results for Pb-Pb collisions at $\sqrt{s_{\rm NN}}= 5.02$ TeV and to the calculation of a transport model. A good agreement is found between Xe-Xe and Pb-Pb results as well as between data and the model.

An invariant differential cross section measurement of inclusive $\pi^{0}$ and $\eta$ meson production at mid-rapidity in pp collisions at $\sqrt{s}=8$ TeV was carried out by the ALICE experiment at the LHC. The spectra of $\pi^{0}$ and $\eta$ mesons were measured in transverse momentum ranges of $0.3<p_{\rm T}<35$ GeV/$c$ and $0.5<p_{\rm T}<35$ GeV/$c$, respectively. Next-to-leading order perturbative QCD calculations using fragmentation functions DSS14 for the $\pi^{0}$ and AESSS for the $\eta$ overestimate the cross sections of both neutral mesons, although such calculations agree with the measured $\eta/\pi^{0}$ ratio within uncertainties. The results were also compared with PYTHIA~8.2 predictions for which the Monash~2013 tune yields the best agreement with the measured neutral meson spectra. The measurements confirm a universal behavior of the $\eta/\pi^{0}$ ratio seen for NA27, PHENIX and ALICE data for pp collisions from $\sqrt{s}=27.5$ GeV to $\sqrt{s}=8$ TeV within experimental uncertainties. A relation between the $\pi^{0}$ and $\eta$ production cross sections for pp collisions at $\sqrt{s}=8$ TeV is given by $m_{\rm T}$ scaling for $p_{\rm T}>3.5$ GeV/$c$. However, a deviation from this empirical scaling rule is observed for transverse momenta below $p_{\rm T}<3.5$ GeV/$c$ in the $\eta/\pi^0$ ratio with a significance of $6.2\sigma$.

The ALICE collaboration at the CERN LHC reports novel measurements of jet substructure in pp collisions at $\sqrt{s}$= 7 TeV and central Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV. Jet substructure of track-based jets is explored via iterative declustering and grooming techniques. We present the measurement of the momentum sharing of two-prong substructure exposed via grooming, the $z_{\rm{g}}$, and its dependence on the opening angle, in both pp and Pb-Pb collisions. We also present the first measurement of the distribution of the number of branches obtained in the iterative declustering of the jet, which is interpreted as the number of its hard splittings. In Pb-Pb collisions, we observe a suppression of symmetric splittings at large opening angles and an enhancement of splittings at small opening angles relative to pp collisions, with no significant modification of the number of splittings. The results are compared to predictions from various Monte Carlo event generators to test the role of important concepts in the evolution of the jet in the medium such as color coherence.

Production cross sections of muons from semi-leptonic decays of charm and beauty hadrons were measured at forward rapidity ($2.5<y<4$) in proton--proton (pp) collisions at a centre-of-mass energy $\sqrt{s}=5.02$ TeV with the ALICE detector at the CERN LHC. The results were obtained in an extended transverse momentum interval, $2 < p_{\rm T} < 20$ GeV/$c$, and with an improved precision compared to previous measurements performed in the same rapidity interval at centre-of-mass energies $\sqrt{s}= 2.76$ and 7 TeV. The $p_{\rm T}$- and $y$-differential production cross sections as well as the $p_{\rm T}$-differential production cross section ratios between different centre-of-mass energies and different rapidity intervals are described, within experimental and theoretical uncertainties, by predictions based on perturbative QCD.

We present results on transverse momentum ($p_{\rm T}$) and rapidity ($y$) differential production cross sections, mean transverse momentum and mean transverse momentum square of inclusive J/$\psi$ and $\psi(2S)$ at forward rapidity ($2.5<y<4$) as well as $\psi(2S)$-to-J/$\psi$ cross section ratios. These quantities are measured in pp collisions at center of mass energies $\sqrt{s}=5.02$ and 13 TeV with the ALICE detector. Both charmonium states are reconstructed in the dimuon decay channel, using the muon spectrometer. A comprehensive comparison to inclusive charmonium cross sections measured at $\sqrt{s}=2.76$, 7 and 8 TeV is performed. A comparison to non-relativistic quantum chromodynamics and fixed-order next-to-leading logarithm calculations, which describe prompt and non-prompt charmonium production respectively, is also presented. A good description of the data is obtained over the full $p_{\rm T}$ range, provided that both contributions are summed. In particular, it is found that for $p_{\rm T}>15$ GeV/$c$ the non-prompt contribution reaches up to 50% of the total charmonium yield.

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

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.