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.

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.

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.

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.

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.

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

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.