Measurements of the differential production cross-sections of prompt and non-prompt $J/\psi$ and $\psi(2$S$)$ mesons with transverse momenta between 8 and 360 GeV and rapidity in the range $|y|<2$ are reported. Furthermore, measurements of the non-prompt fractions of $J/\psi$ and $\psi(2$S$)$, and the prompt and non-prompt $\psi(2$S$)$-to-$J/\psi$ production ratios, are presented. The analysis is performed using 140 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data recorded by the ATLAS detector at the LHC during the years 2015-2018.

This paper presents measurements of charged-hadron spectra obtained in $pp$, $p$+Pb, and Pb+Pb collisions at $\sqrt{s}$ or $\sqrt{s_{_\text{NN}}}=5.02$ TeV, and in Xe+Xe collisions at $\sqrt{s_{_\text{NN}}}=5.44$ TeV. The data recorded by the ATLAS detector at the LHC have total integrated luminosities of 25 pb${}^{-1}$, 28 nb${}^{-1}$, 0.50 nb${}^{-1}$, and 3 $\mu$b${}^{-1}$, respectively. The nuclear modification factors $R_{p\text{Pb}}$ and $R_\text{AA}$ are obtained by comparing the spectra in heavy-ion and $pp$ collisions in a wide range of charged-particle transverse momenta and pseudorapidity. The nuclear modification factor $R_{p\text{Pb}}$ shows a moderate enhancement above unity with a maximum at $p_{\mathrm{T}} \approx 3$ GeV; the enhancement is stronger in the Pb-going direction. The nuclear modification factors in both Pb+Pb and Xe+Xe collisions feature a significant, centrality-dependent suppression. They show a similar distinct $p_{\mathrm{T}}$-dependence with a local maximum at $p_{\mathrm{T}} \approx 2$ GeV and a local minimum at $p_{\mathrm{T}} \approx 7$ GeV. This dependence is more distinguishable in more central collisions. No significant $|\eta|$-dependence is found. A comprehensive comparison with several theoretical predictions is also provided. They typically describe $R_\text{AA}$ better in central collisions and in the $p_{\mathrm{T}}$ range from about 10 to 100 GeV.

The coherent photoproduction of $\rm{J/\psi}$ and $\rm{\psi'}$ mesons was measured in ultra-peripheral Pb-Pb collisions at a center-of-mass energy $\sqrt{s_{\mathrm{NN}}}~=~5.02$ TeV with the ALICE detector. Charmonia are detected in the central rapidity region for events where the hadronic interactions are strongly suppressed. The $\rm{J/\psi}$ is reconstructed using the dilepton ($l^{+} l^{-}$) and proton-antiproton decay channels, while for the $\rm{\psi'}$, the dilepton and the $l^{+} l^{-} \pi^{+} \pi^{-}$ decay channels are studied. The analysis is based on an event sample corresponding to an integrated luminosity of about 233 ${\mu b}^{-1}$. The results are compared with theoretical models for coherent $\rm{J/\psi}$ and $\rm{\psi'}$ photoproduction. The coherent cross section is found to be in a good agreement with models incorporating moderate nuclear gluon shadowing of about 0.65 at a Bjorken-$x$ of around $6\times 10^{-4}$, such as the EPS09 parametrization, however none of the models is able to fully describe the rapidity dependence of the coherent $\rm{J/\psi}$ cross section including ALICE measurements at forward rapidity. The ratio of $\rm{\psi'}$ to $\rm{J/\psi}$ coherent photoproduction cross sections was also measured and found to be consistent with the one for photoproduction off protons.

The production of $\Upsilon$ mesons in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5 TeV is measured with the muon spectrometer of the ALICE detector at the LHC. The yields as well as the nuclear modification factors are determined in the forward rapidity region $2.5<y<4.0$, as a function of rapidity, transverse momentum and collision centrality. The results show that the production of the $\Upsilon$(1S) meson is suppressed by a factor of about three with respect to the production in proton-proton collisions. For the first time, a significant signal for the $\Upsilon$(2S) meson is observed at forward rapidity, indicating a suppression stronger by about a factor 2-3 with respect to the ground state. The measurements are compared with transport, hydrodynamic, comover and statistical hadronisation model calculations.

The transverse momentum ($p_{\rm{T}}$) differential yields of (anti-)$^3$He and (anti-)$^3$H measured in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV with ALICE at the Large Hadron Collider (LHC) are presented. The ratios of the $p_{\rm{T}}$-integrated yields of (anti-)$^3$He and (anti-)$^3$H to the proton yields are reported, as well as the $p_{\rm{T}}$ dependence of the coalescence parameters $B_3$ for (anti-)$^3$He and (anti-)$^3$H. For (anti-)$^3$He, the results obtained in four classes of the mean charged-particle multiplicity density are also discussed. These results are compared to predictions from a canonical statistical hadronization model and coalescence approaches. An upper limit on the total yield of $^4\bar{\mathrm{He}}$ is determined.

The ALICE Collaboration has measured the energy dependence of exclusive photoproduction of $J/\psi$ vector mesons off proton targets in ultra-peripheral p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV. The e$^+$e$^-$ and $\mu^+\mu^-$ decay channels are used to measure the cross section as a function of the rapidity of the $J/\psi$ in the range $-2.5 < y < 2.7$, corresponding to an energy in the $\gamma$p centre-of-mass in the interval $40 < W_{\gamma\mathrm{p}}<550$ GeV. The measurements, which are consistent with a power law dependence of the exclusive $J/\psi$ photoproduction cross section, are compared to previous results from HERA and the LHC and to several theoretical models. They are found to be compatible with previous measurements.

Measurements of the yield and nuclear modification factor, $R_\mathrm{ AA}$, for inclusive jet production are performed using 0.49 nb$^{-1}$ of Pb+Pb data at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV and 25 pb$^{-1}$ of $pp$ data at $\sqrt{s}=5.02$ TeV with the ATLAS detector at the LHC. Jets are reconstructed with the anti-$k_t$ algorithm with radius parameter $R=0.4$ and are measured over the transverse momentum range of 40-1000 GeV in six rapidity intervals covering $|y|<2.8$. The magnitude of $R_\mathrm{ AA}$ increases with increasing jet transverse momentum, reaching a value of approximately 0.6 at 1 TeV in the most central collisions. The magnitude of $R_\mathrm{ AA}$ also increases towards peripheral collisions. The value of $R_\mathrm{ AA}$ is independent of rapidity at low jet transverse momenta, but it is observed to decrease with increasing rapidity at high transverse momenta.

Inclusive four-jet events produced in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 7$ TeV are analysed for the presence of hard double-parton scattering using data corresponding to an integrated luminosity of 37.3 pb$^{-1}$, collected with the ATLAS detector at the LHC. The contribution of hard double-parton scattering to the production of four-jet events is extracted using an artificial neural network, assuming that hard double-parton scattering can be approximated by an uncorrelated overlaying of dijet events. For events containing at least four jets with transverse momentum $p_{\mathrm{T}} \geq 20$ GeV and pseudorapidity $\eta \leq 4.4$, and at least one having $p_{\mathrm{T}} \geq 42.5$ GeV, the contribution of hard double-parton scattering is estimated to be $f_{\mathrm{DPS}} = 0.092 ^{+0.005}_{-0.011} (\mathrm{stat.}) ^{+0.033}_{-0.037} (\mathrm{syst.})$. After combining this measurement with those of the inclusive dijet and four-jet cross-sections in the appropriate phase space regions, the effective overlap area between the interacting protons, $\sigma_{\mathrm{eff}}$, was determined to be $\sigma_{\mathrm{eff}} = 14.9 ^{+1.2}_{-1.0} (\mathrm{stat.}) ^{+5.1}_{-3.8} (\mathrm{syst.})$ mb. This result is consistent within the quoted uncertainties with previous measurements of $\sigma_{\mathrm{eff}}$, performed at centre-of-mass energies between 63 GeV and 8 TeV using various final states, and it corresponds to $21^{+7}_{-6}$% of the total inelastic cross-section measured at $\sqrt{s} = 7$ TeV. The distributions of the observables sensitive to the contribution of hard double-parton scattering, corrected for detector effects, are also provided.

The production cross sections of the prompt charmed mesons D$^0$, D$^+$, D$^{*+}$ and D$_s$ were measured at mid-rapidity in p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE detector at the LHC. D mesons were reconstructed from their decays D$^0\rightarrow{\rm K}^-\pi^+$, D$^+\rightarrow{\rm K}^-\pi^+\pi^+$, D$^{*+}\rightarrow D^0\pi^+$, D$_s^+\rightarrow\phi\pi^+\rightarrow{\rm K}^-{\rm K}^+\pi^+$, and their charge conjugates. The $p_{\rm T}$-differential production cross sections were measured at mid-rapidity in the interval $1<p_{\rm T}<24$ GeV/$c$ for D$^0$, D$^+$ and D$^{*+}$ mesons and in $2<p_{\rm T}<12$ GeV/$c$ for D$_s$ mesons, using an analysis method based on the selection of decay topologies displaced from the interaction vertex. The production cross sections of the D$^0$, D$^+$ and D$^{*+}$ mesons were also measured in three $p_{\rm T}$ intervals as a function of the rapidity $y_{\rm cms}$ in the centre-of-mass system in $-1.26<y_{\rm cms}<0.34$. In addition, the prompt D$^0$ cross section was measured in pp collisions at $\sqrt{s}=7$ TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV down to $p_{\rm T}=0$ using an analysis technique that is based on the estimation and subtraction of the combinatorial background, without reconstruction of the D$^0$ decay vertex. The nuclear modification factor $R_{\rm pPb}(p_{\rm T})$, defined as the ratio of the $p_{\rm T}$-differential D-meson cross section in p-Pb collisions and that in pp collisions scaled by the mass number of the Pb nucleus, was calculated for the four D-meson species and found to be compatible with unity within experimental uncertainties. The results are compared to theoretical calculations that include cold-nuclear-matter effects and to transport model calculations incorporating the interactions of charm quarks with an expanding deconfined medium.

Measurements of the per-event charged-particle yield as a function of the charged-particle transverse momentum and rapidity are performed using $p+$Pb collision data collected by the ATLAS experiment at the LHC at a centre-of-mass energy of $\sqrt{s_{NN}}=5.02$ TeV. Charged particles are reconstructed over pseudorapidity $|\eta|<2.3$ and transverse momentum between $0.1$ GeV and $22$ GeV in a dataset corresponding to an integrated luminosity of $1$ $\mu b^{-1}$. The results are presented in the form of charged-particle nuclear modification factors, where the $p+$Pb charged-particle multiplicities are compared between central and peripheral $p+$Pb collisions as well as to charged-particle cross sections measured in pp collisions. The $p+$Pb collision centrality is characterized by the total transverse energy measured in $-4.9<\eta<-3.1$, which is in the direction of the outgoing lead beam. Three different estimations of the number of nucleons participating in the $p+$Pb collision are carried out using the Glauber model and two Glauber-Gribov colour-fluctuation extensions to the Glauber model. The values of the nuclear modification factors are found to vary significantly as a function of rapidity and transverse momentum. A broad peak is observed for all centralities and rapidities in the nuclear modification factors for charged-particle transverse momentum values around $3$ GeV. The magnitude of the peak increases for more central collisions as well as rapidity ranges closer to the direction of the outgoing lead nucleus.