We present the first systematic comparison of the charged-particle pseudorapidity densities for three widely different collision systems, pp, p-Pb, and Pb-Pb, at the top energy of the Large Hadron Collider ($\sqrt{s_{\rm NN}} = 5.02$ TeV) measured over a wide pseudorapidity range (${-3.5 <\eta <5}$), the widest possible among the four experiments at that facility. The systematic uncertainties are minimised since the measurements are recorded by the same experimental apparatus (ALICE). The distributions for p-Pb and Pb-Pb collisions are determined as a function of the centrality of the collisions, while results from pp collisions are reported for inelastic events with at least one charged particle at midrapidity. The charged-particle pseudorapidity densities are, under simple and robust assumptions, transformed to charged-particle rapidity densities. This allows for the calculation and the presentation of the evolution of the width of the rapidity distributions and of a lower bound on the Bjorken energy density, as a function of the number of participants in all three collision systems. We find a decreasing width of the particle production, and roughly a smooth ten fold increase in the energy density, as the system size grows, which is consistent with a gradually higher dense phase of matter.

The study of the azimuthal anisotropy of inclusive muons produced in p-Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV, using the ALICE detector at the LHC is reported. The measurement of the second-order Fourier coefficient of the particle azimuthal distribution, $v_2$, is performed as a function of transverse momentum $p_{\rm T}$ in the 0-20% high-multiplicity interval at both forward ($2.03 < y_{\rm CMS} < 3.53$) and backward ($-4.46 < y_{\rm CMS} < -2.96$) rapidities over a wide $p_{\rm T}$ range, $0.5 < p_{\rm T} < 10$ GeV/$c$, in which a dominant contribution of muons from heavy-flavour hadron decays is expected at $p_{\rm T} > 2$ GeV/$c$. The $v_2$ coefficient of inclusive muons is extracted using two different techniques, namely two-particle cumulants, used for the first time for heavy-flavour measurements, and forward-central two-particle correlations. Both techniques give compatible results. A positive $v_2$ is measured at both forward and backward rapidities with a significance larger than $4.7\sigma$ and $7.6\sigma$, respectively, in the interval $2 < p_{\rm T} < 6$ GeV/$c$. Comparisons with previous measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, and with AMPT and CGC-based theoretical calculations are discussed. The findings impose new constraints on the theoretical interpretations of the origin of the collective behaviour in small collision systems.

The azimuthal ($\Delta\varphi$) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p$-$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV. Results are reported for electrons with transverse momentum $4<p_{\rm T}<16$ GeV/$c$ and pseudorapidity $|\eta|<0.6$. The associated charged particles are selected with transverse momentum $1<p_{\rm T}<7$ GeV/$c$, and relative pseudorapidity separation with the leading electron $|\Delta\eta| < 1$. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p$-$Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The $\Delta\varphi$ distribution and the peak observables in pp and p$-$Pb collisions are compared with calculations from various Monte Carlo event generators.

The transverse-momentum ($p_{\rm T}$) spectra of ${\rm K}^*(892)^0$ and $\phi(1020)$ measured with the ALICE detector up to $p_{\rm T}$ = 16 GeV/$c$ in the rapidity range $-1.2 < y < 0.3$, in p-Pb collisions at the center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured $p_{\rm T}$ distributions show a dependence on both multiplicity and rapidity at low $p_{\rm T}$ whereas no significant dependence is observed at high $p_{\rm T}$. A rapidity dependence is observed in the $p_{\rm T}$-integrated yield (d$N$/d$y$), whereas the mean transverse momentum ($\left< p_{\rm T} \right>$) shows a flat behavior as a function of rapidity. The rapidity asymmetry ($Y_{\rm asym}$) at low $p_{\rm T}$ ( < 5 GeV/$c$) is more significant for higher multiplicity classes. At high $p_{\rm T}$, no significant rapidity asymmetry is observed in any of the multiplicity classes. Both ${\rm K}^*(892)^0$ and $\phi(1020)$ show similar $Y_{\rm asym}$. The nuclear modification factor ($Q_{\rm CP}$) as a function of $p_{\rm T}$ shows a Cronin-like enhancement at intermediate $p_{\rm T}$, which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high $p_{\rm T}$ (> 5 GeV/$c$), the $Q_{\rm CP}$ values are greater than unity and no significant rapidity dependence is observed.

Neutral pion ($\pi^{0}$) and $\eta$ meson production cross sections were measured up to unprecedentedly high transverse momenta ($p_{\rm T}$) in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The mesons were reconstructed via their two-photon decay channel in the rapidity interval $-1.3< y <0.3$ in the ranges of $0.4<p_{\rm T}<200$ GeV/$c$ and $1.0<p_{\rm T}<50$ GeV/$c$, respectively. The respective nuclear modification factor ($R_{\rm pPb}$) is presented for $p_{\rm T}$ up to of 200 and 30 GeV/$c$, where the former was achieved by extending the $\pi^{0}$ measurement in pp collisions at $\sqrt{s}$ = 8 TeV using the merged cluster technique. The values of $R_{\rm pPb}$ are below unity for $p_{\rm T}<10$ GeV/$c$, while they are consistent with unity for $p_{\rm T}>10$ GeV/$c$, leaving essentially no room for final state energy loss. The new data provide strong constraints for nuclear parton distribution and fragmentation functions over a broad kinematic range and are compared to model predictions as well as previous results at $\sqrt{s_{\rm NN}}$ = 5.02 TeV.

The measurement of the production of charm jets, identified by the presence of a ${\rm D^0}$ meson in the jet constituents, is presented in proton-proton collisions at centre-of-mass energies of $\sqrt{s}$ = 5.02 and 13 TeV with the ALICE detector at the CERN LHC. The ${\rm D^0}$ mesons were reconstructed from their hadronic decay ${\rm D^0} \rightarrow {\rm K^-}\pi^+$ and the respective charge conjugate. Jets were reconstructed from ${\rm D^0}$-meson candidates and charged particles using the anti-$k_{\rm T}$ algorithm, in the jet transverse momentum range $5<p_{\rm T;chjet}<50$ GeV/$c$, pseudorapidity $|\eta_{\rm jet}| <0.9-R$, and with the jet resolution parameters $R$ = 0.2, 0.4, 0.6. The distribution of the jet momentum fraction carried by a ${\rm D^0}$ meson along the jet axis ($z^{\rm ch}_{||}$) was measured in the range $0.4 < z^{\rm ch}_{||} < 1.0$ in four ranges of the jet transverse momentum. Comparisons of results for different collision energies and jet resolution parameters are also presented. The measurements are compared to predictions from Monte Carlo event generators based on leading-order and next-to-leading-order perturbative quantum chromodynamics calculations. A generally good description of the main features of the data is obtained in spite of a few discrepancies at low $p_{\rm T;chjet}$. Measurements were also done for $R = 0.3$ at $\sqrt{s}$ = 5.02 TeV and are shown along with their comparisons to theoretical predictions in an appendix to this paper.

The production of prompt $\rm \Lambda_{\rm c}^{+}$ baryons at midrapidity ($|y|<0.5$) was measured in central (0-10%) and mid-central (30-50%) Pb-Pb collisions at the center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector. The results are more precise, more differential in centrality, and reach much lower transverse momentum ($p_{\rm T}=1$ GeV/$c$) with respect to previous measurements performed by the ALICE, STAR, and CMS Collaborations in nucleus-nucleus collisions, allowing for an extrapolation down to $p_{\rm T}=0$. The $p_{\rm T}$-differential $\rm \Lambda_{\rm c}^{+}$/D$^0$ production ratio is enhanced with respect to the pp measurement for $4<p_{\rm T}<8$ GeV/$c$ by 3.7 standard deviations ($\sigma$), while the $p_{\rm T}$-integrated ratios are compatible within 1$\sigma$. The observed trend is similar to that observed in the strange sector for the $\Lambda/$K$^0_{\rm S}$ ratio. Model calculations including coalescence or statistical hadronization for charm-hadron formation are compared with the data.

The production yields of non-prompt $\mathrm{D_s}^{+}$ mesons, namely $\mathrm{D_s}^{+}$ mesons from beauty-hadron decays, were measured for the first time as a function of the transverse momentum ($p_{\rm T}$) at midrapidity ($|y|<0.5$) in central and semi-central Pb$-$Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE experiment at the LHC. The $\mathrm{D_s}^{+}$ mesons and their charge conjugates were reconstructed from the hadronic decay channel $\mathrm{D_s}^{+} \to \phi\pi^+$, with $\phi \to {\rm K}^-{\rm K}^+$, in the $4 < p_{\rm T}<36$ GeV/$c$ and $2 < p_{\rm T}<24$ GeV/$c$ intervals for the 0$-$10% and 30$-$50% centrality classes, respectively. The measured yields of non-prompt $\mathrm{D_s}^{+}$ mesons are compared to those of prompt $\mathrm{D_s}^{+}$ and non-prompt $\mathrm{D^0}$ mesons by calculating the ratios of the production yields in Pb$-$Pb collisions and the nuclear modification factor $R_\mathrm{AA}$. The ratio between the $R_\mathrm{AA}$ of non-prompt $\mathrm{D_s}^{+}$ and prompt $\mathrm{D_s}^{+}$ mesons, and that between the $R_\mathrm{AA}$ of non-prompt $\mathrm{D_s}^{+}$ and non-prompt $\mathrm{D^0}$ mesons in central Pb$-$Pb collisions are found to be on average higher than unity in the $4< p_{\rm T}<12$ GeV/$c$ interval with a statistical significance of about $1.6\,\sigma$ and $1.7\,\sigma$, respectively. The measured $R_\mathrm{AA}$ ratios are compared with the predictions of theoretical models of heavy-quark transport in a hydrodynamically expanding QGP that incorporate hadronisation via quark recombination.

The cross section for coherent photonuclear production of J/$\psi$ is presented as a function of the electromagnetic dissociation (EMD) of Pb. The measurement is performed with the ALICE detector in ultra-peripheral Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV. Cross sections are presented in five different J/$\psi$ rapidity ranges within $|y|<4$, with the J/$\psi$ reconstructed via its dilepton decay channels. In some events the J/$\psi$ is not accompanied by EMD, while other events do produce neutrons from EMD at beam rapidities either in one or the other beam direction, or in both. The cross sections in a given rapidity range and for different configurations of neutrons from EMD allow for the extraction of the energy dependence of this process in the range $17 < W_{\gamma\, \mathrm{Pb, n}} <920$ GeV, where $W_{\gamma\, \mathrm{Pb, n}}$ is the centre-of-mass energy per nucleon of the $\gamma\,\mathrm{Pb}$ system. This range corresponds to a Bjorken-$x$ interval spanning about three orders of magnitude: $ 1.1\times10^{-5}<x<3.3\times 10^{-2}$. In addition to the ultra-peripheral and photonuclear cross sections, the nuclear suppression factor is obtained. These measurements point to a strong depletion of the gluon distribution in Pb nuclei over a broad, previously unexplored, energy range. These results, together with previous ALICE measurements, provide unprecedented information to probe quantum chromodynamics at high energies.

Understanding the role of parton mass and Casimir colour factors in the quantum chromodynamics parton shower represents an important step in characterising the emission properties of heavy quarks. Recent experimental advances in jet substructure techniques have provided the opportunity to isolate and characterise gluon emissions from heavy quarks. In this work, the first direct experimental constraint on the charm-quark splitting function is presented, obtained via the measurement of the groomed shared momentum fraction of the first splitting in charm jets, tagged by a reconstructed ${\rm D}^0$ meson. The measurement is made in proton-proton collisions at $\sqrt{s}$ = 13 TeV, in the low jet transverse-momentum interval of $15 \leq p_{\rm T}^{\rm jet~ch} < 30$ GeV/$c$ where the emission properties are sensitive to parton mass effects. In addition, the opening angle of the first perturbative emission of the charm quark, as well as the number of perturbative emissions it undergoes, is reported. Comparisons to measurements of an inclusive-jet sample show a steeper splitting function for charm quarks compared with gluons and light quarks. Charm quarks also undergo fewer perturbative emissions in the parton shower, with a reduced probability of large-angle emissions.