ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015--2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons ($N_\mu>4$) and in the zenith angle range $0^{\circ}<\theta<50^{\circ}$. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about $ 4 \times 10^{15}<E_\mathrm{prim}< 6 \times 10^{16}$~eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET-II-04 is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL 2.3d and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than $20\%$ and $30\%$, respectively. The rate of high muon multiplicity events ($N_\mu>100$) obtained with QGSJET-II-04 and SIBYLL 2.3d is compatible with the data, while EPOS-LHC produces a significantly lower rate ($55\%$ of the measured rate). For both QGSJET-II-04 and SIBYLL 2.3d, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy $E_\mathrm{prim} \sim 10^{17}$ eV of these events. This result places significant constraints on more exotic production mechanisms.

The production yields of the orbitally excited charm-strange mesons $\mathrm{D_{s1}(1^{+})(2536)^+}$ and $\mathrm{D_{s2}^{*}(2^{+})(2573)^+}$were measured for the first time in proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}=13$ TeV with the ALICE experiment at the LHC. The $\mathrm{D_{s1}^+}$ and $\mathrm{D_{s2}^{*+}}$ mesons were measured at midrapidity ($|y|<0.5$) in minimum-bias and high-multiplicity pp collisions in the transverse-momentum interval $2<p_{\rm T}<24$ GeV/$c$. Their production yields relative to the $\mathrm{D_{s}^{+}}$ ground-state yield were found to be compatible between minimum-bias and high-multiplicity collisions, as well as with previous measurements in ${\rm e}^{\pm}$p and ${\rm e^+ e^-}$ collisions. The measured $\mathrm{D_{s1}^+/D_{s}^{+}}$ and $\mathrm{D_{s2}^{*+}/D_{s}^{+}}$ yield ratios are described by statistical hadronization models and can be used to tune the parameters governing the production of excited charm-strange hadrons in Monte Carlo generators, such as PYTHIA 8.

We report on the measurement of inclusive, non-prompt, and prompt J/$ψ$-hadron correlations by the ALICE Collaboration at the CERN Large Hadron Collider in pp collisions at a center-of-mass energy of $13~\text{TeV}$. The correlations are studied at midrapidity ($|y| < 0.9$) in the transverse momentum ranges $p_{\rm T} < 40~\text{GeV}/c$ for the J/$ψ$ and $0.15 < p_{\rm T} < 10~\text{GeV}/c$ and $|η|<0.9$ for the associated hadrons. The measurement is based on minimum bias and high multiplicity data samples corresponding to integrated luminosities of $L_{\text{int}} = 34~\text{nb}^{-1}$ and $L_{\text{int}} = 6.9~\text{pb}^{-1}$, respectively. In addition, two more data samples are employed, requiring, on top of the minimum bias condition, a threshold on the tower energy of $E = 4$ and $9~\text{GeV}$ in the ALICE electromagnetic calorimeters, which correspond to integrated luminosities of $L_{\text{int}} = 0.9~\text{pb}^{-1}$ and $L_{\text{int}} = 8.4~\text{pb}^{-1}$, respectively. The azimuthally integrated near and away side yields of associated charged hadrons per J/$ψ$ trigger are presented as a function of the J/$ψ$ and associated hadron transverse momentum. The measurements are discussed in comparison to PYTHIA calculations.

Differential two-particle normalized cumulants ($R_2$) and transverse momentum correlations ($P_2$) are measured as a function of the relative pseudorapidity and azimuthal angle difference $( Δη, Δφ)$ of charged particle pairs in minimum bias pp collisions at $\sqrt{\textit{s}}$ = 13 TeV. The measurements use charged hadrons in the pseudorapidity region of $|η| < 0.8$ and the transverse momentum range \mbox{0.2 $< \textit{p}_{\mathrm T} < $ 2.0 $\mathrm{GeV}/\textit{c}$} in order to focus on soft multiparticle interactions and to complement prior measurements of these correlation functions in p-Pb and Pb-Pb collisions. The correlation functions are reported for both unlike-sign and like-sign pairs and their charge-independent and charge-dependent combinations. Both the $R_2$ and $P_2$ measured in pp collisions exhibit features qualitatively similar to those observed in p--Pb and Pb--Pb collisions. The $Δη$ and $Δφ$ root mean square widths of the near-side peak of the correlation functions are evaluated and compared with those observed in p-Pb and Pb-Pb collisions and show smooth evolution with the multiplicity of charged particles produced in the collision. The comparison of the measured correlation functions with predictions from PYTHIA8 shows that this model qualitatively captures their basic structure and characteristics but feature important differences. In addition, the $R_2^{\rm CD}$ is used to determine the charge balance function of hadrons produced within the detector acceptance of the measurements. The integral of the balance function is found to be compatible with those reported by a previous measurement in Pb--Pb collisions.

Ultrarelativistic heavy-ion collisions create a quark-gluon plasma (QGP), a hot and dense state of strongly interacting QCD matter. In ultracentral collisions, the QGP volume remains nearly constant event-by-event, while its total entropy can fluctuate due to quantum effects, leading to temperature variations. These features allow the correlation between the mean transverse momentum $(\langle p_{\mathrm{T}} \rangle)$ of charged hadrons and their multiplicity to serve as a probe of the QGP's speed of sound, $c_{s}$. This study extracts $c_{s}$ by analyzing the relative increase in $\langle p_{\mathrm{T}} \rangle$ with respect to the charged-particle density $(\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle)$ at midrapidity in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$, using data from the ALICE detector. Centrality is determined with estimators based on multiplicity and transverse energy, applying a pseudorapidity gap to reduce selection biases. The extracted value of $c_{s}^{2}$ is found to strongly depend on the employed centrality estimator and ranges between $0.1146 \pm 0.0028 \,\mathrm{(stat.)} \pm 0.0065 \,\mathrm{(syst.)}$ and $0.4374 \pm 0.0006 \mathrm{(stat.)} \pm 0.0184 \mathrm{(syst.)}$ in natural units. Additionally, the event-by-event $[p_{\mathrm{T}}]$ distribution is studied through its variance, skewness, and kurtosis. A pronounced decrease in the self-normalized variance and a peak followed by a drop in skewness suggest the suppression of impact-parameter fluctuations in ultracentral collisions. These observations provide new insights into the thermodynamic properties and initial-state fluctuations of the QGP.

The transverse momentum spectra and integrated yields of $\overlineΣ^{\pm}$ have been measured in pp and p-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV with the ALICE experiment. Measurements are performed via the newly accessed decay channel $\overlineΣ^{\pm} \rightarrow {\rm\overline{n}}π^{\pm}$. A new method of antineutron reconstruction with the PHOS electromagnetic spectrometer is developed and applied to this analysis. The $p_{\rm T}$ spectra of $\overlineΣ^{\pm}$ are measured in the range $0.5 < p_{\rm T} < 3$ GeV/$c$ and compared to predictions of the PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4 models. The EPOS LHC and EPOS4 models provide the best descriptions of the measured spectra both in pp and p-Pb collisions, while models which do not account for multiparton interactions provide a considerably worse description at high $p_{\rm T}$. The total yields of $\overlineΣ^{\pm}$ in both pp and p-Pb collisions are compared to predictions of the Thermal-FIST model and dynamical models PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4. All models reproduce the total yields in both colliding systems within uncertainties. The nuclear modification factors $R_{\rm pPb}$ for both $\overlineΣ^{+}$ and $\overlineΣ^{-}$ are evaluated and compared to those of protons, $Λ$ and $Ξ$ hyperons, and predictions of EPOS LHC and EPOS4 models. No deviations of $R_{\rm pPb}$ for $\overlineΣ^{\pm}$ from the model predictions or measurements for other hadrons are found within uncertainties.

The coherent J/$\psi$ photoproduction cross section is measured for the first time at midrapidity in peripheral to semicentral Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV. The centrality differential cross section ${\rm d} \sigma/ {\rm d}y$ is reported for the centrality range 40-90%, together with the doubly-differential cross section ${\rm d}^2 \sigma /{\rm d}y {\rm d} p_{\rm T}$, extracted in two peripheral centrality classes. The J/$\psi$ mesons are reconstructed in the dielectron channel, in the rapidity interval $|y| <$ 0.9 using the ALICE central barrel detectors. The J/$\psi$ cross section at midrapidity is statistically compatible to the earlier ALICE measurement at forward rapidity and at the same centre-of-mass energy, and shows only a mild centrality dependence over the covered range. Several sets of theoretical calculations taking into account the hadronic overlap in the collisions but ignoring possible final-state effects from a hot expanding medium are found to give a fairly good description of the current measurements within uncertainties.

The transverse momentum ($p_{\rm T}$) differential production cross section of the promptly produced charm-strange baryon $\Xi_{\rm c}^{0}$ (and its charge conjugate $\overline{\Xi_{\rm c}^{0}}$) is measured at midrapidity via its hadronic decay into ${\rm \pi^{+}}\Xi^{-}$ in p$-$Pb collisions at a centre-of-mass energy per nucleon$-$nucleon collision $\sqrt{s_{\mathrm{NN}}}~=~5.02$ TeV with the ALICE detector at the LHC. The $\Xi_{\rm c}^{0}$ nuclear modification factor ($R_{\rm pPb}$), calculated from the cross sections in pp and p$-$Pb collisions, is presented and compared with the $R_{\rm pPb}$ of $\Lambda_{\rm c}^{+}$ baryons. The ratios between the $p_{\rm T}$-differential production cross section of $\Xi_{\rm c}^{0}$ baryons and those of $\mathrm {D^0}$ mesons and $\Lambda_{\rm c}^{+}$ baryons are also reported and compared with results at forward and backward rapidity from the LHCb Collaboration. The measurements of the production cross section of prompt $\Xi^0_{\rm c}$ baryons are compared with a model based on perturbative QCD calculations of charm-quark production cross sections, which includes only cold nuclear matter effects in p$-$Pb collisions, and underestimates the measurement by a factor of about 50. This discrepancy is reduced when the data is compared with a model that includes string formation beyond leading-colour approximation or in which hadronisation is implemented via quark coalescence. The $p_{\rm T}$-integrated cross section of prompt $\Xi^0_{\rm c}$-baryon production at midrapidity extrapolated down to $p_{\rm T}$ = 0 is also reported. These measurements offer insights and constraints for theoretical calculations of the hadronisation process. Additionally, they provide inputs for the calculation of the charm production cross section in p$-$Pb collisions at midrapidity.

This Letter presents the first measurement of the impact-parameter dependent angular anisotropy in the decay of coherently photoproduced $\rho^0$ mesons. The $\rho^0$ mesons are reconstructed through their decay into pion pairs. The measured anisotropy corresponds to the amplitude of the $\cos(2\phi)$ modulation, where $\phi$ is the angle between the two vectors formed by the sum and the difference of the transverse momenta of the pions, respectively. The measurement was performed by the ALICE Collaboration at the LHC using data from ultraperipheral Pb$-$Pb collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}}~=~5.02$ TeV per nucleon pair. Different impact-parameter regions are selected by classifying the events in nuclear-breakup classes. The amplitude of the $\cos(2\phi)$ modulation is found to increase by about one order of magnitude from large to small impact parameters. Theoretical calculations describe the measured $\cos(2\phi)$ anisotropy and its impact-parameter dependence as the result of a quantum interference effect at the femtometer scale, arising from the ambiguity regarding which of the nuclei is the photon source in the interaction.

The $p_{\rm T}$-differential production cross sections of non-prompt ${\rm D^0}$, ${\rm D^+}$, and ${\rm D_s^+}$ mesons originating from beauty-hadron decays are measured in proton$-$proton collisions at a centre-of-mass energy $\sqrt{s}=13$ TeV. The measurements are performed at midrapidity, $|y| < 0.5$, with the data sample collected by ALICE from 2016 to 2018. The results are in agreement with predictions from several perturbative QCD calculations. The fragmentation fraction of beauty quarks to strange mesons divided by the one to non-strange mesons, $f_{\rm{s}}/(f_{\rm{u}} + f_{\rm{d}})$, is found to be $0.114 \pm 0.016~{\rm (stat.)} \pm 0.006~{\rm (syst.)} \pm 0.003~{\rm (BR)} \pm 0.003~{\rm (extrap.)}$. This value is compatible with previous measurements at lower centre-of-mass energies and in different collision systems in agreement with the assumption of universality of fragmentation functions. In addition, the dependence of the non-prompt D meson production on the centre-of-mass energy is investigated by comparing the results obtained at $\sqrt{s} = 5.02$ and 13 TeV, showing a hardening of the non-prompt D-meson $p_{\rm T}$-differential production cross section at higher $\sqrt{s}$. Finally, the ${\rm b\overline{b}}$ production cross section per unit of rapidity at midrapidity is calculated from the non-prompt ${\rm D^0}$, ${\rm D^+}$, ${\rm D_s^+}$, and $\Lambda_{\rm c}^+$ hadron measurements, obtaining ${\rm d}\sigma/{\rm d}y = 75.2\pm 3.2~(\mathrm{stat.}) \pm 5.2~(\mathrm{syst.})^{+12.3}_{-3.2} ~(\mathrm{extrap.})\text{ } \rm \mu b \;.$