The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.

This article reports measurements characterizing the Underlying Event (UE) associated with hard scatterings at midrapidity in pp collisions at $\sqrt{s}=13$ TeV. The hard scatterings are identified by the leading particle, the charged particle with the highest transverse momentum ($p_{\rm T}^{\rm leading}$) in the event. Charged-particle number and summed transverse-momentum densities are measured in different azimuthal regions defined with respect to the leading particle direction: Toward, Transverse, and Away. The Toward and Away regions contain the fragmentation products of the hard scatterings in addition to the UE contribution, whereas particles in the Transverse region are expected to originate predominantly from the UE. The study is performed as a function of $p_{\rm T}^{\rm leading}$ with three different $p_{\rm T}$ thresholds for the associated particles, $p_{\rm T}^{\rm min} >$ 0.15, 0.5, and 1.0 GeV/$c$. The charged-particle density in the Transverse region rises steeply for low values of $p_{\rm T}^{\rm leading}$ and reaches a plateau. The results confirm the trend that the charged-particle density in the Transverse region shows a stronger increase with $\sqrt{s}$ than the inclusive charged-particle density at midrapidity. The UE activity is increased by approximately 20% when going from 7 to 13 TeV. The plateau in the Transverse region ($5 < p_{\rm T}^{\rm leading} < ~ 40$ GeV/$c$ ) is further characterized by the probability distribution of its charged-particle multiplicity normalized to its average value (relative transverse activity, $R_{T}$) and the mean transverse momentum as a function of $R_{T}$. Experimental results are compared to model calculations using PYTHIA 8 and EPOS LHC. The overall agreement between models and data is within 30%. These measurements provide new insights on the interplay between hard scatterings and the associated UE in pp collisions.

This Letter presents the first direct investigation of the p-$\Sigma^{0}$ interaction, using the femtoscopy technique in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV measured by the ALICE detector. The $\Sigma^{0}$ is reconstructed via the decay channel to $\Lambda \gamma$, and the subsequent decay of $\Lambda$ to p$\pi^-$. The photon is detected via the conversion in material to e$^{+}$e$^{-}$ pairs exploiting the unique capability of the ALICE detector to measure electrons at low transverse momenta. The measured p-$\Sigma^{0}$ correlation indicates a shallow strong interaction. The comparison of the data to several theoretical predictions obtained employing the $Correlation~Analysis~Tool~using~the~Schr\"odinger~Equation$ (CATS) and the Lednick\'y-Lyuboshits approach shows that the current experimental precision does not yet allow to discriminate between different models, as it is the case for the available scattering and hypernuclei data. Nevertheless, the p-$\Sigma^{0}$ correlation function is found to be sensitive to the strong interaction, and driven by the interplay of the different spin and isospin channels. This pioneering study demonstrates the feasibility of a femtoscopic measurement in the p-$\Sigma^{0}$ channel and with the expected larger data samples in LHC Run 3 and Run 4, the p-$\Sigma^{0}$ interaction will be constrained with high precision.

Deuteron production in high-energy collisions is sensitive to the space-time evolution of the collision system, and is typically described by a coalescence mechanism. For the first time, we present results on jet-associated deuteron production in \pp collisions at $\sqrt{s}=13$ TeV, providing an opportunity to test the established picture for deuteron production in events with a hard scattering. Using a trigger particle with high transverse-momentum ($p_{\rm T}>5$ GeV/$c$) as a proxy for the presence of a jet at midrapidity, we observe a measurable population of deuterons being produced around the jet proxy. The associated deuteron yield measured in a narrow angular range around the trigger particle differs by 2.4-4.8 standard deviations from the uncorrelated background. The data are described by PYTHIA model calculations featuring baryon coalescence.

A measurement of dielectron production in proton-proton (pp) collisions at $\sqrt{s} = 13$ TeV, recorded with the ALICE detector at the CERN LHC, is presented in this Letter. The data set was recorded with a reduced magnetic solenoid field. This enables the investigation of a kinematic domain at low dielectron invariant mass $m_{\rm ee}$ and pair transverse momentum $p_{\rm T,ee}$ that was previously inaccessible at the LHC. The cross section for dielectron production is studied as a function of $m_{\rm ee}$, $p_{\rm T,ee}$, and event multiplicity ${\rm d} N_{\rm ch}/{\rm d} \eta$. The expected dielectron rate from hadron decays, called hadronic cocktail, utilizes a parametrization of the measured $\eta/\pi^0$ ratio in pp and proton-nucleus (p-A) collisions, assuming that this ratio shows no strong dependence on collision energy at low transverse momentum. Comparison of the measured dielectron yield to the hadronic cocktail at $0.15<m_{\rm ee}<0.6$ GeV/$c^2$ and for $p_{\rm T,ee} < 0.4$ GeV/$c$ indicates an enhancement of soft dielectrons, reminiscent of the 'anomalous' soft-photon and -dilepton excess in hadron-hadron collisions reported by several experiments under different experimental conditions. The enhancement factor over the hadronic cocktail amounts to $1.61\pm 0.13\,(\rm{stat.})\pm 0.17\,(\rm{syst., data}) \pm 0.34\,(\rm{syst., cocktail})$ in the ALICE acceptance. Acceptance-corrected excess spectra in $m_{\rm ee}$ and $p_{\rm T,ee}$ are extracted and compared with calculations of dielectron production from hadronic bremsstrahlung and thermal radiation within a hadronic many-body approach.

Understanding the production mechanism of light (anti)nuclei is one of the key challenges of nuclear physics and has important consequences for astrophysics, since it provides an input for indirect dark-matter searches in space. In this paper, the latest results about the production of light (anti)nuclei in pp collisions at $\sqrt{s} = 13$ TeV are presented, focusing on the comparison with the predictions of coalescence and thermal models. For the first time, the coalescence parameters $B_2$ for deuterons and $B_3$ for helions are compared with parameter-free theoretical predictions that are directly constrained by the femtoscopic measurement of the source radius in the same event class. A fair description of the data with a Gaussian wave function is observed for both deuteron and helion, supporting the coalescence mechanism for the production of light (anti)nuclei in pp collisions. This method paves the way for future investigations of the internal structure of more complex nuclear clusters, including the hypertriton.

The cross section of jets reconstructed from charged particles is measured in the transverse momentum range of $5<p_\mathrm{T}<100\ \mathrm{GeV}/c$ in pp collisions at the center-of-mass energy of $\sqrt{s} = 5.02\ \mathrm{TeV}$ with the ALICE detector. The jets are reconstructed using the anti-$k_\mathrm{T}$ algorithm with resolution parameters $R=0.2$, $0.3$, $0.4$, and $0.6$ in the pseudorapidity range $|\eta|< 0.9-R$. The charged jet cross sections are compared with the leading order (LO) and to next-to-leading order (NLO) perturbative Quantum ChromoDynamics (pQCD) calculations. It was found that the NLO calculations agree better with the measurements. The cross section ratios for different resolution parameters were also measured. These ratios increase from low $p_\mathrm{T}$ to high $p_\mathrm{T}$ and saturate at high $p_\mathrm{T}$, indicating that jet collimation is larger at high $p_\mathrm{T}$ than at low $p_\mathrm{T}$. These results provide a precision test of pQCD predictions and serve as a baseline for the measurement in Pb$-$Pb collisions at the same energy to quantify the effects of the hot and dense medium created in heavy-ion collisions at the LHC.

The first experimental information on the strong interaction between $\Lambda$ and $\Xi^-$ strange baryons is presented in this Letter. The correlation function of $\Lambda-\Xi^-$ and $\overline{\Lambda}-\overline{\Xi}^{+}$ pairs produced in high-multiplicity proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the LHC is measured as a function of the relative momentum of the pair. The femtoscopy method is used to calculate the correlation function, which is then compared with theoretical expectations obtained using a meson exchange model, chiral effective field theory, and Lattice QCD calculations close to the physical point. Data support predictions of small scattering parameters while discarding versions with large ones, thus suggesting a weak $\Lambda-\Xi^{-}$ interaction. The limited statistical significance of the data does not yet allow one to constrain the effects of coupled channels like $\Sigma-\Xi$ and N$-\Omega$.

We report on the inclusive J/$\psi$ production cross section measured at the CERN Large Hadron Collider in proton-proton collisions at a centre-of-mass energy $\sqrt{s}~=~13$ TeV. The J/$\psi$ mesons are reconstructed in the $\rm e^{+} e^{-}$ decay channel and the measurements are performed at midrapidity ($|y|<0.9$) in the transverse-momentum interval $0< p_{\rm T} <40$ GeV/$c$, using a minimum bias data sample corresponding to an integrated luminosity $L_{\text{int}} = 32.2~\text{nb}^{-1}$ and an Electromagnetic Calorimeter triggered data sample with $L_{\text{int}} = 8.3~\mathrm{pb}^{-1}$. The $p_{\rm T}$-integrated J/$\psi$ production cross section at midrapidity, computed using the minimum bias data sample, is $\text{d}\sigma/\text{d}y|_{y=0} = 8.97\pm0.24~(\text{stat})\pm0.48~(\text{syst})\pm0.15~(\text{lumi})~\mu\text{b}$. An approximate logarithmic dependence with the collision energy is suggested by these results and available world data, in agreement with model predictions. The integrated and $p_{\rm T}$-differential measurements are compared with measurements in pp collisions at lower energies and with several recent phenomenological calculations based on the non-relativistic QCD and Color Evaporation models.

Particle production sensitive to non-factorizable and non-perturbative processes that contribute to the underlying event associated with a high transverse momentum ($p_{T}$) jet in proton+proton collisions at $\sqrt{s}$=200 GeV is studied with the STAR detector. Each event is divided into three regions based on the azimuthal angle with respect to the highest-$p_{T}$ jet direction: in the leading jet direction ("Toward"), opposite to the leading jet ("Away"), and perpendicular to the leading jet ("Transverse"). In the Transverse region, the average charged particle density is found to be between 0.4 and 0.6 and the mean transverse momentum, $\langle p_{T}\rangle$, between 0.5-0.7 GeV/$c$ for particles with $p_{T}$$>$0.2 GeV/$c$ at mid-pseudorapidity ($|\eta|$$<$1) and jet $p_{T}$$>$15 GeV/$c$. Both average particle density and $\langle p_{T}\rangle$ depend weakly on the leading jet $p_{T}$. Closer inspection of the Transverse region hints that contributions to the underlying event from initial- and final-state radiation are significantly smaller in these collisions than at the higher energies, up to 13 TeV, recorded at the LHC. Underlying event measurements associated with a high-$p_{T}$ jet will contribute to our understanding of QCD processes at hard and soft scales at RHIC energies, as well as provide constraints to modeling of underlying event dynamics.

The study of the strength and behaviour of the antikaon-nucleon ($\mathrm{\overline{K}N}$) interaction constitutes one of the key focuses of the strangeness sector in low-energy Quantum Chromodynamics (QCD). In this letter a unique high-precision measurement of the strong interaction between kaons and protons, close and above the kinematic threshold, is presented. The femtoscopic measurements of the correlation function at low pair-frame relative momentum of (K$^+$ p $\oplus$ K$^-$ $\overline{\mathrm{p}}$) and (K$^-$ p $\oplus$ K$^+$ $\overline{\mathrm{p}}$) pairs measured in pp collisions at $\sqrt{s}$ = 5, 7 and 13 TeV are reported. A structure observed around a relative momentum of 58 MeV/$c$ in the measured correlation function of (K$^-$ p $\oplus$ K$^+$ $\overline{\mathrm{p}}$) with a significance of 4.4. $\sigma$ constitutes the first experimental evidence for the opening of the $(\mathrm{\overline{K}^0 n} \oplus \mathrm{K^0 \overline{n}})$ isospin breaking channel due to the mass difference between charged and neutral kaons. The measured correlation functions have been compared to J\"{u}lich and Kyoto models in addition to the Coulomb potential. The high-precision data at low relative momenta presented in this work prove femtoscopy to be a powerful complementary tool to scattering experiments and provide new constraints above the $\mathrm{\overline{K}N}$ threshold for low-energy QCD chiral models.

The inclusive production of the charm-strange baryon $\rm \Omega_c^0$ is measured for the first time via its hadronic decay into $\rm \Omega^-\pi^+$ at midrapidity ($|y|<0.5$) in proton-proton (pp) collisions at the centre-of-mass energy $\sqrt{s}=13$ TeV with the ALICE detector at the LHC. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval $2<p_{\rm T}<12~{\rm GeV}/c$. The $p_{\rm T}$ dependence of the $\rm \Omega_c^0$-baryon production relative to the prompt $\rm D^0$-meson and to the prompt $\rm \Xi_c^0$-baryon production is compared to various models that take different hadronisation mechanisms into consideration. In the measured $p_{\rm T}$ interval, the ratio of the $p_{\rm T}$-integrated cross sections of $\rm \Omega_c^0$ and prompt $\Lambda_{\rm c}^{+}$ baryons multiplied by the $\rm \Omega^-\pi^+$ branching ratio is found to be larger by a factor of about 20 with a significance of about 4$\sigma$ when compared to $\rm e^+e^-$ collisions.

The production of $\phi$ mesons has been studied in pp collisions at LHC energies with the ALICE detector via the dimuon decay channel in the rapidity region $2.5 < y < 4$. Measurements of the differential cross section ${\rm d}^2\sigma/{\rm d}y {\rm d}p_{\rm T}$ are presented as a function of the transverse momentum ($p_{\rm T}$) at the center-of-mass energies $\sqrt{s}=5.02$, 8 and 13 TeV and compared with the ALICE results at midrapidity. The differential cross sections at $\sqrt{s}=5.02$ and 13 TeV are also studied in several rapidity intervals as a function of $p_{\rm T}$, and as a function of rapidity in three $p_{\rm T}$ intervals. A hardening of the $p_{\rm T}$-differential cross section with the collision energy is observed, while, for a given energy, $p_{\rm T}$ spectra soften with increasing rapidity and, conversely, rapidity distributions get slightly narrower at increasing $p_{\rm T}$. The new results, complementing the published measurements at $\sqrt{s}=2.76$ and 7 TeV, allow one to establish the energy dependence of $\phi$ meson production and to compare the measured cross sections with phenomenological models. None of the considered models manages to describe the evolution of the cross section with $p_{\rm T}$ and rapidity at all the energies.

The multiplicity dependence of jet production in pp collisions at the centre-of-mass energy of $\sqrt{s} = 13\ \mathrm{TeV}$ is studied for the first time. Jets are reconstructed from charged particles using the anti-$k_\mathrm{T}$ algorithm with resolution parameters $R$ varying from $0.2$ to $0.7$. The jets are measured in the pseudorapidity range $|\eta_{\rm jet}|< 0.9-R$ and in the transverse momentum range $5<p_\mathrm{T,jet}^{\rm ch}<140\ \mathrm{GeV}/c$. The multiplicity intervals are categorised by the ALICE forward detector V0. The $p_{\mathrm{T}}$ differential cross section of charged-particle jets are compared to leading order (LO) and next-to-leading order (NLO) perturbative quantum chromodynamics (pQCD) calculations. It is found that the data are better described by the NLO calculation, although the NLO prediction overestimates the jet cross section below $20\ \mathrm{GeV}/c$. The cross section ratios for different $R$ are also measured and compared to model calculations. These measurements provide insights into the angular dependence of jet fragmentation. The jet yield increases with increasing self-normalised charged-particle multiplicity. This increase shows only a weak dependence on jet transverse momentum and resolution parameter at the highest multiplicity. While such behaviour is qualitatively described by the present version of PYTHIA, quantitative description may require implementing new mechanisms for multi-particle production in hadronic collisions.

Measurements of the production cross sections of prompt ${\rm D^0}$, ${\rm D^+}$, ${\rm D^{\ast +}}$, ${\rm D_s^+}$, ${\rm \Lambda_{c}^{+}}$, and ${\rm \Xi_{c}^{+}}$ charm hadrons at midrapidity in proton$-$proton collisions at $\sqrt{s}=13$ TeV with the ALICE detector are presented. The D-meson cross sections as a function of transverse momentum ($p_{\rm T}$) are provided with improved precision and granularity. The ratios of $p_{\rm T}$-differential meson production cross sections based on this publication and on measurements at different rapidity and collision energy provide a constraint on gluon parton distribution functions at low values of Bjorken-$x$ ($10^{-5}-10^{-4}$). The measurements of ${\rm \Lambda_{c}^{+}}$ (${\rm \Xi_{c}^{+}}$) baryon production extend the measured $p_{\rm T}$ intervals down to $p_{\rm T}=0(3)$~GeV$/c$. These measurements are used to determine the charm-quark fragmentation fractions and the ${\rm c\overline{c}}$ production cross section at midrapidity ($|y|<0.5$) based on the sum of the cross sections of the weakly-decaying ground-state charm hadrons ${\rm D^0}$, ${\rm D^+}$, ${\rm D_s^+}$, ${\rm \Lambda_{c}^{+}}$, ${\rm \Xi_{c}^{0}}$ and, for the first time, ${\rm \Xi_{c}^{+}}$, and of the strongly-decaying J/$psi$ mesons. The first measurements of ${\rm \Xi_{c}^{+}}$ and ${\rm \Sigma_{c}^{0,++}}$ fragmentation fractions at midrapidity are also reported. A significantly larger fraction of charm quarks hadronising to baryons is found compared to e$^+$e$^-$ and ep collisions. The ${\rm c\overline{c}}$ production cross section at midrapidity is found to be at the upper bound of state-of-the-art perturbative QCD calculations.

The azimuthal anisotropy Fourier coefficients ($v_n$) in 8.16 TeV pPb data are extracted via long-range two-particle correlations as a function of event multiplicity and compared to corresponding results in pp and PbPb collisions. Using a four-particle cumulant technique, $v_n$ correlations are measured for the first time in pp and pPb collisions. The $v_2$ and $v_4$ coefficients are found to be positively correlated in all collision systems. For high multiplicity pPb collisions an anticorrelation of $v_2$ and $v_3$ is observed, with a similar correlation strength as in PbPb data at the same multiplicity. The new correlation results strengthen the case for a common origin of the collectivity seen in pPb and PbPb collisions in the measured multiplicity range.

Measurements of the inclusive J/$\psi$ yield as a function of charged-particle pseudorapidity density ${\rm d}N_{\rm ch}/\rm{d}\eta$ in pp collisions at $\sqrt{s}$ = 13 TeV with ALICE at the LHC are reported. The J/$\psi$ meson yield is measured at midrapidity ($|y|<0.9$) in the dielectron channel, for events selected based on the charged-particle multiplicity at midrapidity ($|\eta|<1$) and at forward rapidity ($-3.7<\eta<-1.7$ and $2.8<\eta<5.1$); both observables are normalized to their corresponding averages in minimum bias events. The increase of the normalized J/$\psi$ yield with normalized ${\rm d}N_{\rm ch}/\rm{d}\eta$ is significantly stronger than linear and dependent on the transverse momentum. The data are compared to theoretical predictions, which describe the observed trends well, albeit not always quantitatively.

The production rates and the transverse momentum distribution of strange hadrons at mid-rapidity ($\ |y\ | < 0.5$) are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV as a function of the charged particle multiplicity, using the ALICE detector at the LHC. The production rates of $\rm{K}^{0}_{S}$, $\Lambda$, $\Xi$, and $\Omega$ increase with the multiplicity faster than what is reported for inclusive charged particles. The increase is found to be more pronounced for hadrons with a larger strangeness content. Possible auto-correlations between the charged particles and the strange hadrons are evaluated by measuring the event-activity with charged particle multiplicity estimators covering different pseudorapidity regions. When comparing to lower energy results, the yields of strange hadrons are found to depend only on the mid-rapidity charged particle multiplicity. Several features of the data are reproduced qualitatively by general purpose QCD Monte Carlo models that take into account the effect of densely-packed QCD strings in high multiplicity collisions. However, none of the tested models reproduce the data quantitatively. This work corroborates and extends the ALICE findings on strangeness production in proton-proton collisions at 7 TeV.

This paper presents the measurements of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{p}$ and $\bar{\rm{p}}$ transverse momentum ($p_{\rm{T}}$) spectra as a function of charged-particle multiplicity density in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector at the LHC. Such study allows us to isolate the center-of-mass energy dependence of light-flavour particle production. The measurements reported here cover a $p_{\rm{T}}$ range from 0.1 GeV/$c$ to 20 GeV/$c$ and are done in the rapidity interval $|y|<0.5$. The $p_{\rm{T}}$-differential particle ratios exhibit an evolution with multiplicity, similar to that observed in pp collisions at $\sqrt{s}$ = 7 TeV, which is qualitatively described by some of the hydrodynamical and pQCD-inspired models discussed in this paper. Furthermore, the $p_{\rm{T}}$-integrated hadron-to-pion yield ratios measured in pp collisions at two different center-of-mass energies are consistent when compared at similar multiplicities. This also extends to strange and multistrange hadrons, suggesting that, at LHC energies, particle hadrochemistry scales with particle multiplicity the same way under different collision energies and colliding systems.

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

The production cross section of inclusive J/$\psi$ pairs in pp collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV is measured with ALICE. The measurement is performed for J/$\psi$ in the rapidity interval $2.5 < y < 4.0$ and for transverse momentum $p_{\rm T} > 0$. The production cross section of inclusive J/$\psi$ pairs is reported to be $10.3 \pm 2.3 {\rm (stat.)} \pm 1.3 {\rm (syst.)}$ nb in this kinematic interval. The contribution from non-prompt J/$\psi$ (i.e. originated from beauty-hadron decays) to the inclusive sample is evaluated. The results are discussed and compared with data.

Two-particle azimuthal correlations are measured with the ALICE apparatus in pp collisions at $\sqrt{s} = 13$ TeV to explore strangeness- and multiplicity-related effects in the fragmentation of jets and the transition regime between bulk and hard production, probed with the condition that a strange meson (K$^{0}_{\rm S}$) or baryon ($\Lambda$) with transverse momentum $p_{\rm T} > 3$ GeV/c is produced. Azimuthal correlations between kaons or $\Lambda$ hyperons with other hadrons are presented at midrapidity for a broad range of the trigger ($3 < p_{\rm T}^{\rm trigg} < 20$ GeV/$c$) and associated particle $p_{\rm T}$ (1 GeV/$c$$< p_{\rm T}^{\rm assoc} < p_{\rm T}^{\rm trigg}$), for minimum-bias events and as a function of the event multiplicity. The near- and away-side peak yields are compared for the case of either K$^{0}_{\rm S}$ or $\Lambda$($\overline{\Lambda}$) being the trigger particle with that of inclusive hadrons (a sample dominated by pions). In addition, the measurements are compared with predictions from PYTHIA 8 and EPOS LHC event generators.

Measurements of the second Fourier harmonic coefficient ($v_2$) of the azimuthal distributions of prompt and nonprompt D$^0$ mesons produced in pp and pPb collisions are presented. Nonprompt D$^0$ mesons come from beauty hadron decays. The data samples are collected by the CMS experiment at nucleon-nucleon center-of-mass energies of 13 and 8.16 TeV, respectively. In high multiplicity pp collisions, $v_2$ signals for prompt charm hadrons are reported for the first time, and are found to be comparable to those for light-flavor hadron species over a transverse momentum ($p_\mathrm{T}$) range of 2-6 GeV. Compared at similar event multiplicities, the prompt D$^0$ meson $v_2$ values in pp and pPb collisions are similar in magnitude. The $v_2$ values for open beauty hadrons are extracted for the first time via nonprompt D$^0$ mesons in pPb collisions. For $p_\mathrm{T}$ in the range of 2-5 GeV, the results suggest that $v_2$ for nonprompt D$^0$ mesons are smaller than those for prompt D$^0$ mesons. These new measurements indicate a positive charm hadron $v_2$ in pp collisions and suggest a mass dependence in $v_2$ between charm and beauty hadrons in the pPb system. These results provide insights into the origin of heavy-flavor quark collectivity in small systems.

Production of inclusive charmonia in pp collisions at center-of-mass energy of $\sqrt{s}$ = 13 TeV and p-Pb collisions at center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}$ = 8.16 TeV is studied as a function of charged-particle pseudorapidity density with ALICE. Ground and excited charmonium states (J/$\psi$, $\psi$(2S)) are measured from their dimuon decays in the interval of rapidity in the center-of-mass frame $2.5 < y_{\rm cms} < 4.0$ for pp collisions, and $2.03 < y_{\rm cms} < 3.53$ and $-4.46 < y_{\rm cms} < -2.96$ for p-Pb collisions. The charged-particle pseudorapidity density is measured around midrapidity ($|\eta|<1.0$). In pp collisions, the measured charged-particle multiplicity extends to about six times the average value, while in p-Pb collisions at forward (backward) rapidity a multiplicity corresponding to about three (four) times the average is reached. The $\psi$(2S) yield increases with the charged-particle pseudorapidity density. The ratio of $\psi$(2S) over J/$\psi$ yield does not show a significant multiplicity dependence in either colliding system, suggesting a similar behavior of J/$\psi$ and $\psi$(2S) yields with respect to charged-particle pseudorapidity density. Results for the $\psi$(2S) yield and its ratio with respect to J/$\psi$ agree with available model calculations.

This work presents new constraints on the existence and the binding energy of a possible $\Lambda$-$\Lambda$ bound state, the H-dibaryon, derived from $\Lambda$-$\Lambda$ femtoscopic measurements by the ALICE collaboration. The results are obtained from a new measurement using the femtoscopy technique in pp collisions at $\sqrt{s}=13$ TeV and p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV, combined with previously published results from p-Pb collisions at $\sqrt{s}=7$ TeV. The $\Lambda$-$\Lambda$ scattering parameter space, spanned by the inverse scattering length $f_0^{-1}$ and the effective range $d_0$, is constrained by comparing the measured $\Lambda$-$\Lambda$ correlation function with calculations obtained within the Lednicky model. The data are compatible with hypernuclei results and lattice computations, both predicting a shallow attractive interaction, and permit to test different theoretical approaches describing the $\Lambda$-$\Lambda$ interaction. The region in the $(f_0^{-1},d_0)$ plane which would accommodate a $\Lambda$-$\Lambda$ bound state is substantially restricted compared to previous studies. The binding energy of the possible $\Lambda$-$\Lambda$ bound state is estimated within an effective-range expansion approach and is found to be $B_{\Lambda\Lambda}=3.2^{+1.6}_{-2.4}\mathrm{(stat)}^{+1.8}_{-1.0}\mathrm{(syst)}$ MeV.

Measurements of the production of electrons from heavy-flavour hadron decays in pp collisions at $\sqrt{s} = 13$ TeV at midrapidity with the ALICE detector are presented down to a transverse momentum ($p_{\rm T}$) of 0.2 GeV$/c$ and up to $p_{\rm T} = 35$ GeV$/c$, which is the largest momentum range probed for inclusive electron measurements in ALICE. In p$-$Pb collisions, the production cross section and the nuclear modification factor of electrons from heavy-flavour hadron decays are measured in the $p_{\rm T}$ range $0.5 < p_{\rm T} < 26$ GeV$/c$ at $\sqrt{s_{\rm NN}} = 8.16$ TeV. The nuclear modification factor is found to be consistent with unity within the statistical and systematic uncertainties. In both collision systems, first measurements of the yields of electrons from heavy-flavour hadron decays in different multiplicity intervals normalised to the multiplicity-integrated yield (self-normalised yield) at midrapidity are reported as a function of the self-normalised charged-particle multiplicity estimated at midrapidity. The self-normalised yields in pp and p$-$Pb collisions grow faster than linear with the self-normalised multiplicity. A strong $p_{\rm T}$ dependence is observed in pp collisions, where the yield of high-$p_{\rm T}$ electrons increases faster as a function of multiplicity than the one of low-$p_{\rm T}$ electrons. The measurement in p$-$Pb collisions shows no $p_{\rm T}$ dependence within uncertainties. The self-normalised yields in pp and p$-$Pb collisions are compared with measurements of other heavy-flavour, light-flavour, and strange particles, and with Monte Carlo simulations.

Measurements of identified hadrons as a function of the charged-particle multiplicity in pp collisions enable a search for the onset of collective effects in small collision systems. With such measurements, it is possible to study the mechanisms that determine the shapes of hadron transverse momentum ($p_{\rm{T}}$) spectra, to search for possible modifications of the yields of short-lived hadronic resonances due to scattering effects in the hadron-gas phase, and to investigate different explanations for the multiplicity evolution of strangeness production provided by phenomenological models. In this paper, these topics are addressed through measurements of the $\rm{K}^{*}(892)^{0}$ and $\phi(1020)$ mesons at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV as a function of the charged-particle multiplicity. The results include the $p_{\rm{T}}$ spectra, $p_{\rm{T}}$-integrated yields, mean transverse momenta, and the ratios of the yields of these resonances to those of longer-lived hadrons. Comparisons with results from other collision systems and energies, as well as predictions from phenomenological models, are also discussed.

Results on two-particle angular correlations for charged particles produced in pp collisions at a center-of-mass energy of 13 TeV are presented. The data were taken with the CMS detector at the LHC and correspond to an integrated luminosity of about 270 inverse nanobarns. The correlations are studied over a broad range of pseudorapidity (abs(eta) < 2.4) and over the full azimuth (phi) as a function of charged particle multiplicity and transverse momentum (pt). In high-multiplicity events, a long-range (abs(Delta eta) > 2.0), near-side (Delta phi approximately 0) structure emerges in the two-particle Delta eta-Delta phi correlation functions. The magnitude of the correlation exhibits a pronounced maximum in the range 1.0 < pt < 2.0 GeV/c and an approximately linear increase with the charged particle multiplicity, with an overall correlation strength similar to that found in earlier pp data at sqrt(s) = 7 TeV. The present measurement extends the study of near-side long-range correlations up to charged particle multiplicities of N[ch] approximately 180, a region so far unexplored in pp collisions. The observed long-range correlations are compared to those seen in pp, pPb, and PbPb collisions at lower collision energies.

The pseudorapidity ($\eta$) and transverse-momentum ($p_{\rm T}$) distributions of charged particles produced in proton-proton collisions are measured at the centre-of-mass energy $\sqrt{s}$ = 13 TeV. The pseudorapidity distribution in $|\eta|<$ 1.8 is reported for inelastic events and for events with at least one charged particle in $|\eta|<$ 1. The pseudorapidity density of charged particles produced in the pseudorapidity region $|\eta|<$ 0.5 is 5.31 $\pm$ 0.18 and 6.46 $\pm$ 0.19 for the two event classes, respectively. The transverse-momentum distribution of charged particles is measured in the range 0.15 $<$ $p_{\rm T}$ $<$ 20 GeV/c and $|\eta|<$ 0.8 for events with at least one charged particle in $|\eta|<$ 1. The correlation between transverse momentum and particle multiplicity is also investigated by studying the evolution of the spectra with event multiplicity. The results are compared with calculations from PYTHIA and EPOS Monte Carlo generators.

The pseudorapidity density of charged particles with minimum transverse momentum ($p_{\rm T}$) thresholds of 0.15, 0.5, 1, and 2 GeV$/c$ is measured in pp collisions at the centre of mass energies of $\sqrt{s} =$ 5.02 and 13 TeV with the ALICE detector. The study is carried out for inelastic collisions with at least one primary charged particle having a pseudorapidity ($\eta$) within $\pm0.8$ and $p_{\rm T}$ larger than the corresponding threshold. In addition, measurements without $p_{\rm T}$-thresholds are performed for inelastic and non-single-diffractive events as well as for inelastic events with at least one charged particle having $|\eta|<1$ in pp collisions at $\sqrt{s} =$ 5.02 TeV for the first time at the LHC. These measurements are compared to the PYTHIA 6, PYTHIA 8, and EPOS-LHC models. In general, the models describe the $\eta$ dependence of particle production well. However, discrepancies are observed for the highest transverse momentum threshold ($p_{\rm T}>2 {\rm\ GeV}/c$), highlighting the importance of such measurements for tuning event generators. The new measurements agree within uncertainties with results from the ATLAS and CMS experiments obtained at $\sqrt{s} = 13$ TeV.

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 production cross sections of ${\rm D^0}$ and $\Lambda^+_{\rm c}$ hadrons originating from beauty-hadron decays (i.e. non-prompt) were measured for the first time at midrapidity ($|y|<0.5$) by the ALICE Collaboration in proton-proton collisions at a center-of-mass energy $\sqrt{s}=13$ TeV. They are described within uncertainties by perturbative QCD calculations employing the fragmentation fractions of beauty quarks to baryons measured at forward rapidity by the LHCb Collaboration. The ${\rm b\overline{b}}$ production cross section per unit of rapidity at midrapidity, estimated from these measurements, is ${\rm d}\sigma_{\rm b\overline{b}}/{\rm d}y|_{|y|<0.5} = 83.1 \pm 3.5 (\mathrm{stat.}) \pm 5.4(\mathrm{syst.}) ^{+12.3}_{-3.2} (\mathrm{extrap.})\,\mu$b. The baryon-to-meson ratios are computed to investigate the hadronization mechanism of beauty quarks. The non-prompt $\Lambda^+_{\rm c}/{\rm D^0}$ production ratio has a similar trend to the one measured for the promptly produced charmed particles and to the p$/\pi^+$ and $\Lambda/{\rm K^0_S}$ ratios, suggesting a similar baryon-formation mechanism among light, strange, charm, and beauty hadrons. The $p_{\rm T}$-integrated non-prompt $\Lambda_{\rm c}/{\rm D^0}$ ratio is found to be significantly higher than the one measured in e$^+$e$^-$ collisions.

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.

The first study of the shapes of jets arising from bottom (b) quarks in heavy ion collisions is presented. Jet shapes are studied using charged hadron constituents as a function of their radial distance from the jet axis. Lead-lead (PbPb) collision data at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV were recorded by the CMS detector at the LHC, with an integrated luminosity of 1.69 nb$^{-1}$. Compared to proton-proton collisions, a redistribution of the energy in b jets to larger distances from the jet axis is observed in PbPb collisions. This medium-induced redistribution is found to be substantially larger for b jets than for inclusive jets.

Measurements of the elliptic ($v_2$) and triangular ($v_3$) azimuthal anisotropy coefficients are presented for D$^0$ mesons produced in b hadron decays (nonprompt D$^0$ mesons) in lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The results are compared with previously published charm meson anisotropies measured using prompt D$^0$ mesons. The data were collected with the CMS detector in 2018 with an integrated luminosity of 0.58 nb$^{-1}$. Azimuthal anisotropy is sensitive to the interactions of quarks with the hot and dense medium created in heavy ion collisions. Comparing results for prompt and nonprompt D$^0$ mesons can assist in understanding the mass dependence of these interactions. The nonprompt results show lower magnitudes of $v_2$ and $v_3$ and weaker dependences on the meson transverse momentum and collision centrality than those found for prompt D$^0$ mesons. The results are in agreement with theoretical predictions that include a mass dependence in the interactions of quarks with the medium.

Multi-particle cumulants and corresponding Fourier harmonics are measured for azimuthal angle distributions of charged particles in $pp$ collisions at $\sqrt{s}$ = 5.02 and 13 TeV and in $p$+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, and compared to the results obtained for low-multiplicity Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. These measurements aim to assess the collective nature of particle production. The measurements of multi-particle cumulants confirm the evidence for collective phenomena in $p$+Pb and low-multiplicity Pb+Pb collisions. On the other hand, the $pp$ results for four-particle cumulants do not demonstrate collective behaviour, indicating that they may be biased by contributions from non-flow correlations. A comparison of multi-particle cumulants and derived Fourier harmonics across different collision systems is presented as a function of the charged-particle multiplicity. For a given multiplicity, the measured Fourier harmonics are largest in Pb+Pb, smaller in $p$+Pb and smallest in $pp$ collisions. The $pp$ results show no dependence on the collision energy, nor on the multiplicity.

The measurement of dielectron production is presented as a function of invariant mass and transverse momentum ($p_{\rm T}$) at midrapidity ($|y_{\rm e}|<0.8$) in proton-proton (pp) collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The contributions from light-hadron decays are calculated from their measured cross sections in pp collisions at $\sqrt{s}=7$ TeV or 13 TeV. The remaining continuum stems from correlated semileptonic decays of heavy-flavour hadrons. Fitting the data with templates from two different MC event generators, PYTHIA and POWHEG, the charm and beauty cross sections at midrapidity are extracted for the first time at this collision energy: ${\rm d}\sigma_{\rm c\bar{c}}/{\rm d}y|_{y=0}=974\pm138(\rm{stat.})\pm140(\rm{syst.})~\mu{\rm b}$ and ${\rm d}\sigma_{\rm b\bar{b}}/{\rm d}y|_{y=0}=79\pm14(\rm{stat.})\pm11(\rm{syst.})~\mu{\rm b}$ using PYTHIA simulations and ${\rm d}\sigma_{\rm c\bar{c}}/{\rm d}y|_{y=0}=1417\pm184(\rm{stat.})\pm204(\rm{syst.})~\mu{\rm b}$ and ${\rm d}\sigma_{\rm b\bar{b}}/{\rm d}y|_{y=0}=48\pm14(\rm{stat.})\pm7(\rm{syst.})~\mu{\rm b}$ for POWHEG. These values, whose uncertainties are fully correlated between the two generators, are consistent with extrapolations from lower energies. The different results obtained with POWHEG and PYTHIA imply different kinematic correlations of the heavy-quark pairs in these two generators. Furthermore, comparisons of dielectron spectra in inelastic events and in events collected with a trigger on high charged-particle multiplicities are presented in various $p_{\rm T}$ intervals. The differences are consistent with the already measured scaling of light-hadron and open-charm production at high charged-particle multiplicity as a function of $p_{\rm T}$. Upper limits for the contribution of virtual direct photons are extracted at 90% confidence level and found to be in agreement with pQCD calculations.

Collective behavior has been observed in high-energy heavy-ion collisions for several decades. Collectivity is driven by the high particle multiplicities that are produced in these collisions. At the Large Hadron Collider (LHC), features of collectivity have also been seen in high-multiplicity proton-proton collisions that can attain particle multiplicities comparable to peripheral Pb-Pb collisions. One of the possible signatures of collective behavior is the decrease of femtoscopic radii extracted from pion and kaon pairs emitted from high-multiplicity collisions with increasing pair transverse momentum. This decrease can be described in terms of an approximate transverse mass scaling. In the present work, femtoscopic analyses are carried out by the ALICE collaboration on charged pion and kaon pairs produced in pp collisions at $\sqrt{s}=13$ TeV from the LHC to study possible collectivity in pp collisions. The event-shape analysis method based on transverse sphericity is used to select for spherical versus jet-like events, and the effects of this selection on the femtoscopic radii for both charged pion and kaon pairs are studied. This is the first time this selection method has been applied to charged kaon pairs. An approximate transverse-mass scaling of the radii is found in all multiplicity ranges studied when the difference in the Lorentz boost for pions and kaons is taken into account. This observation does not support the hypothesis of collective expansion of hot and dense matter that should only occur in high-multiplicity events. A possible alternate explanation of the present results is based on a scenario of common emission conditions for pions and kaons in pp collisions for the multiplicity ranges studied.

This letter reports the first measurement of spin alignment, with respect to the helicity axis, for ${\rm D^{*+}}$ vector mesons and their charge conjugates from charm-quark hadronisation (prompt) and from beauty-meson decays (non-prompt) in hadron collisions. The measurements were performed at midrapidity ($|y|<0.8$) as a function of transverse momentum ($p_{\rm T}$) in proton-proton (pp) collisions collected by ALICE at the centre-of-mass energy $\sqrt{s} = 13$ TeV. The diagonal spin density matrix element $\rho_{00}$ of ${\rm D^{*+}}$ mesons was measured from the angular distribution of the ${\rm D^{*+}} \to {\rm D}^0 (\to {\rm K}^{-}\pi^{+}) \pi^+$ decay products, in the ${\rm D^{*+}}$ rest frame, with respect to the ${\rm D^{*+}}$ momentum direction in the pp centre of mass frame. The $\rho_{00}$ value for prompt ${\rm D^{*+}}$ mesons is consistent with $1/3$, which implies no spin alignment. However, for non-prompt ${\rm D^{*+}}$ mesons an evidence of $\rho_{00}$ larger than $1/3$ is found. The measured value of the spin density element is $\rho_{00}=0.455\pm0.022(\text{stat.})\pm0.035(\text{syst.})$ in the $5 < p_{\rm T} < 20$ GeV/$c$ interval, which is consistent with a PYTHIA 8 Monte Carlo simulation coupled with the EVTGEN package, which implements the helicity conservation in the decay of ${\rm D^{*+}}$ meson from beauty mesons. In non-central heavy-ion collisions, the spin of the ${\rm D^{*+}}$ mesons may be globally aligned with the direction of the initial angular momentum and magnetic field. Based on the results for pp collisions reported in this letter it is shown that alignment of non-prompt ${\rm D^{*+}}$ mesons due to the helicity conservation coupled to the collective anisotropic expansion may mimic the signal of global spin alignment in heavy-ion collisions.

Correlations of two flow harmonics $v_n$ and $v_m$ via three- and four-particle cumulants are measured in 13 TeV $pp$, 5.02 TeV $p$+Pb, and 2.76 TeV peripheral Pb+Pb collisions with the ATLAS detector at the LHC. The goal is to understand the multi-particle nature of the long-range collective phenomenon in these collision systems. The large non-flow background from dijet production present in the standard cumulant method is suppressed using a method of subevent cumulants involving two, three and four subevents separated in pseudorapidity. The results show a negative correlation between $v_2$ and $v_3$ and a positive correlation between $v_2$ and $v_4$ for all collision systems and over the full multiplicity range. However, the magnitudes of the correlations are found to depend strongly on the event multiplicity, the choice of transverse momentum range and collision system. The relative correlation strength, obtained by normalisation of the cumulants with the $\langle v_n^2\rangle$ from a two-particle correlation analysis, is similar in the three collision systems and depends weakly on the event multiplicity and transverse momentum. These results based on the subevent methods provide strong evidence of a similar long-range multi-particle collectivity in $pp$, $p$+Pb and peripheral Pb+Pb collisions.

Measurements of anisotropic flow coefficients ($v_n$) and their cross-correlations using two- and multi-particle cumulant methods are reported in collisions of pp at $\sqrt{s} = 13$ TeV, p-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV, Xe-Xe at $\sqrt{s_{_{\rm NN}}} = 5.44$ TeV, and Pb-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV recorded with the ALICE detector. The multiplicity dependence of $v_n$ is studied in a very wide range from 20 to 3000 particles produced in the mid-rapidity region $|\eta|<0.8$ for the transverse momentum range $0.2 < p_{\rm T} < 3.0$ GeV/$c$. An ordering of the coefficients $v_2 > v_3 > v_4$ is found in pp and p-Pb collisions, similar to that seen in large collision systems, while a weak $v_2$ multiplicity dependence is observed relative to nucleus-nucleus collisions in the same multiplicity range. Using a novel subevent method, $v_{2}$ measured with four-particle cumulants is found to be compatible with that from six-particle cumulants in pp and p-Pb collisions. The magnitude of the correlation between $v_n^2$ and $v_m^2$, evaluated with the symmetric cumulants SC$(m,n)$ is observed to be positive at all multiplicities for $v_2$ and $v_4$, while for $v_2$ and $v_3$ it is negative and changes sign for multiplicities below 100, which may indicate a different $v_n$ fluctuation pattern in this multiplicity range. The observed long-range multi-particle azimuthal correlations in high multiplicity pp and p-Pb collisions can neither be described by PYTHIA 8 nor by IP-Glasma+MUSIC+UrQMD model calculations, and hence provide new insights into the understanding of collective effects in small collision systems.

The production of J/$\psi$ is measured at midrapidity ($|y|<0.9$) in proton-proton collisions at $\sqrt{s}$ = 5.02 and 13 TeV, through the dielectron decay channel, using the ALICE detector at the Large Hadron Collider. The data sets used for the analyses correspond to integrated luminosities of $\mathcal{L}_{\rm int}$ = 19.4 $\pm$ 0.4 nb$^{-1}$ and $\mathcal{L}_{\rm int}$ = 32.2 $\pm$ 0.5 nb$^{-1}$ at $\sqrt{s}$ = 5.02 and 13 TeV, respectively. The fraction of non-prompt J/$\psi$ mesons, i.e. those originating from the decay of beauty hadrons, is measured down to a transverse momentum $p_{\rm T}$ = 2 GeV/$c$ (1 GeV/$c$) at $\sqrt{s}$ = 5.02 TeV (13 TeV). The $p_{\rm T}$ and rapidity ($y$) differential cross sections, as well as the corresponding values integrated over $p_{\rm T}$ and $y$, are carried out separately for prompt and non-prompt J/$\psi$ mesons. The results are compared with measurements from other experiments and theoretical calculations based on quantum chromodynamics (QCD). The shapes of the $p_{\rm T}$ and $y$ distributions of beauty quarks predicted by state-of-the-art perturbative QCD models are used to extrapolate an estimate of the $\rm b\overline{b}$ pair cross section at midrapidity and in the total phase space. The total $\rm b\overline{b}$ cross sections are found to be $\sigma_{\rm b \overline{\rm b}} = 541 \pm 45 (\rm stat.) \pm 69 (\rm syst.)_{-12}^{+10} (\rm extr.)~{\rm \mu b}$ and $\sigma_{\rm b \overline{\rm b}}~=~218 \pm 37 (\rm stat.) \pm 31 (\rm syst.)_{-9.1}^{+8.2} (\rm extr.)~{\rm \mu b}$ at $\sqrt{s}$ = 13 and 5.02 TeV, respectively. The value obtained from the combination of ALICE and LHCb measurements in pp collisions at $\sqrt{s}$ = 13 TeV is also provided.

The production of J/$\psi$ is measured as a function of charged-particle multiplicity at forward rapidity in proton$-$proton (pp) collisions at center-of-mass energies $\sqrt{s} =$ 5.02 and 13 TeV. The J/$\psi$ mesons are reconstructed via their decay into dimuons in the rapidity interval (2.5 $< y <$ 4.0), whereas the charged-particle multiplicity density (${\rm d}N_{\rm{ch}}/{\rm d}\eta$) is measured at midrapidity $(|\eta| < 1)$. The production rate as a function of multiplicity is reported as the ratio of the yield in a given multiplicity interval to the multiplicity-integrated one. This observable shows a linear increase with charged-particle multiplicity normalized to the corresponding average value for inelastic events (${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$), at both the colliding energies. Measurements are compared with available ALICE results at midrapidity and theoretical model calculations. First measurement of the mean transverse momentum ($\langle p_{\mathrm{T}}\rangle$) of J/$\psi$ in pp collisions exhibits an increasing trend as a function of ${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$ showing a saturation towards high charged-particle multiplicities.

Two-particle angular correlations are measured in high-multiplicity proton-proton collisions at $\sqrt{s} =13$ TeV by the ALICE Collaboration. The yields of particle pairs at short-($\Delta\eta$$\sim$ 0) and long-range ($1.6 < |\Delta\eta| < 1.8$) in pseudorapidity are extracted on the near-side ($\Delta\varphi$$\sim$ 0). They are reported as a function of transverse momentum ($p_{\mathrm T}$) in the range $1 <p_{\mathrm T}< 4$ GeV/$c$. Furthermore, the event-scale dependence is studied for the first time by requiring the presence of high-$p_{\rm T}$ leading particles and jets for varying $p_{\rm T}$ thresholds. The results demonstrate that the long-range "ridge" yield, possibly related to the collective behavior of the system, is present in events with high-$p_{\mathrm T}$ processes. The magnitudes of the short- and long-range yields are found to grow with the event scale. The results are compared to EPOS LHC and PYTHIA 8 calculations, with and without string-shoving interactions. It is found that while both models describe the qualitative trends in the data, calculations from EPOS LHC show a better quantitative agreement, in particular for the $p_{\rm T}$ and event-scale dependencies.

The fractions of non-prompt (i.e. originating from beauty-hadron decays) D$^0$ and D$^+$ mesons with respect to the inclusive yield are measured as a function of the charged-particle multiplicity in proton$-$proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV with the ALICE detector at the LHC. The results are reported in intervals of transverse momentum ($p_{\rm T}$) and integrated in the range $1 < p_{\rm T} < 24$ GeV/$c$. The fraction of non-prompt D$^0$ and D$^+$ mesons is found to increase slightly as a function of $p_{\rm T}$ in all the measured multiplicity intervals, while no significant dependence on the charged-particle multiplicity is observed. In order to investigate the production and hadronisation mechanisms of charm and beauty quarks, the results are compared to PYTHIA 8 as well as EPOS 3 and EPOS 4 Monte Carlo simulations, and to calculations based on the colour glass condensate including three-pomeron fusion.

The production of $\pi^\pm$, ${\rm K}^\pm$, and $(\overline{\rm p})$p is measured in pp collisions at $\sqrt{s}=13$ TeV in different topological regions. Particle transverse momentum ($p_{\rm T}$) spectra are measured in the ``toward'', ``transverse'', and ``away'' angular regions defined with respect to the direction of the leading particle in the event. While the toward and away regions contain the fragmentation products of the near-side and away-side jets, respectively, the transverse region is dominated by particles from the Underlying Event (UE). The relative transverse activity classifier, $R_{\rm T}=N_{\rm T}/\langle N_{\rm T}\rangle$, is used to group events according to their UE activity, where $N_{\rm T}$ is the measured charged-particle multiplicity per event in the transverse region and $\langle N_{\rm T}\rangle$ is the mean value over all the analysed events. The first measurements of identified particle $p_{\rm T}$ spectra as a function of $R_{\rm T}$ in the three topological regions are reported. The yield of high transverse momentum particles relative to the $R_{\rm T}$-integrated measurement decreases with increasing $R_{\rm T}$ in both the toward and away regions, indicating that the softer UE dominates particle production as $R_{\rm T}$ increases and validating that $R_{\rm T}$ can be used to control the magnitude of the UE. Conversely, the spectral shapes in the transverse region harden significantly with increasing $R_{\rm T}$. This hardening follows a mass ordering, being more significant for heavier particles. The $p_{\rm T}$-differential particle ratios $({\rm p+\overline{p}})/(\pi^+ +\pi^-)$ and $({\rm K^+ +K^-})/(\pi^+ +\pi^-)$ in the low UE limit $(R_{\rm T}\rightarrow 0)$ approach expectations from Monte Carlo generators such as PYTHIA 8 with Monash 2013 tune and EPOS LHC, where the jet-fragmentation models have been tuned to reproduce ${\rm e^+ e^-}$ results.

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.

Angular correlations of heavy-flavour and charged particles in high-energy proton-proton collisions are sensitive to the production mechanisms of heavy quarks and to their fragmentation as well as hadronisation processes. The measurement of the azimuthal-correlation function of prompt D mesons with charged particles in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV with the ALICE detector is reported, considering D$^{0}$, D$^{+}$, and D$^{*+}$ mesons in the transverse-momentum interval $3 < p_{\rm T} < 36$ GeV/$c$ at midrapidity ($|y| < 0.5$), and charged particles with $p_{\rm T} > 0.3$ GeV/$c$ and pseudorapidity $|\eta| < 0.8$. This measurement has an improved precision and provides an extended transverse-momentum coverage compared to previous ALICE measurements at lower energies. The study is also performed as a function of the charged-particle multiplicity, showing no modifications of the correlation function with multiplicity within uncertainties. The properties and the transverse-momentum evolution of the near- and away-side correlation peaks are studied and compared with predictions from various Monte Carlo event generators. Among those considered, PYTHIA8 and POWHEG+PYTHIA8 provide the best description of the measured observables. The obtained results can provide guidance on tuning the generators.