Jet quenching is the process of color-charged partons losing energy via interactions with quark-gluon plasma droplets created in heavy-ion collisions. The collective expansion of such droplets is well described by viscous hydrodynamics. Similar evidence of collectivity is consistently observed in smaller collision systems, including $pp$ and $p$+Pb collisions. In contrast, while jet quenching is observed in Pb+Pb collisions, no evidence has been found in these small systems to date, raising fundamental questions about the nature of the system created in these collisions. The ATLAS experiment at the Large Hadron Collider has measured the yield of charged hadrons correlated with reconstructed jets in 0.36 nb$^{-1}$ of $p$+Pb and 3.6 pb$^{-1}$ of $pp$ collisions at 5.02 TeV. The yields of charged hadrons with $p_\mathrm{T}^\mathrm{ch} >0.5$ GeV near and opposite in azimuth to jets with $p_\mathrm{T}^\mathrm{jet} > 30$ or $60$ GeV, and the ratios of these yields between $p$+Pb and $pp$ collisions, $I_{p\mathrm{Pb}}$, are reported. The collision centrality of $p$+Pb events is categorized by the energy deposited by forward neutrons from the struck nucleus. The $I_{p\mathrm{Pb}}$ values are consistent with unity within a few percent for hadrons with $p_\mathrm{T}^\mathrm{ch} >4$ GeV at all centralities. These data provide new, strong constraints which preclude almost any parton energy loss in central $p$+Pb collisions.

Measurements of the suppression and correlations of dijets is performed using 3 $\mu$b$^{-1}$ of Xe+Xe data at $\sqrt{s_{\mathrm{NN}}} = 5.44$ TeV collected with the ATLAS detector at the LHC. Dijets with jets reconstructed using the $R=0.4$ anti-$k_t$ algorithm are measured differentially in jet $p_{\text{T}}$ over the range of 32 GeV to 398 GeV and the centrality of the collisions. Significant dijet momentum imbalance is found in the most central Xe+Xe collisions, which decreases in more peripheral collisions. Results from the measurement of per-pair normalized and absolutely normalized dijet $p_{\text{T}}$ balance are compared with previous Pb+Pb measurements at $\sqrt{s_{\mathrm{NN}}} =5.02$ TeV. The differences between the dijet suppression in Xe+Xe and Pb+Pb are further quantified by the ratio of pair nuclear-modification factors. The results are found to be consistent with those measured in Pb+Pb data when compared in classes of the same event activity and when taking into account the difference between the center-of-mass energies of the initial parton scattering process in Xe+Xe and Pb+Pb collisions. These results should provide input for a better understanding of the role of energy density, system size, path length, and fluctuations in the parton energy loss.

Studies of the correlations of the two highest transverse momentum (leading) jets in individual Pb+Pb collision events can provide information about the mechanism of jet quenching by the hot and dense matter created in such collisions. In Pb+Pb and pp collisions at $\sqrt{s_{_\text{NN}}}$ = 5.02 TeV, measurements of the leading dijet transverse momentum ($p_{\mathrm{T}}$) correlations are presented. Additionally, measurements in Pb+Pb collisions of the dijet pair nuclear modification factors projected along leading and subleading jet $p_{\mathrm{T}}$ are made. The measurements are performed using the ATLAS detector at the LHC with 260 pb$^{-1}$ of pp data collected in 2017 and 2.2 nb$^{-1}$ of Pb+Pb data collected in 2015 and 2018. An unfolding procedure is applied to the two-dimensional leading and subleading jet $p_{\mathrm{T}}$ distributions to account for experimental effects in the measurement of both jets. Results are provided for dijets with leading jet $p_{\mathrm{T}}$ greater than 100 GeV. Measurements of the dijet-yield-normalized $x_{\mathrm{J}}$ distributions in Pb+Pb collisions show an increased fraction of imbalanced jets compared to pp collisions; these measurements are in agreement with previous measurements of the same quantity at 2.76 TeV in the overlapping kinematic range. Measurements of the absolutely-normalized dijet rate in Pb+Pb and pp collisions are also presented, and show that balanced dijets are significantly more suppressed than imbalanced dijets in Pb+Pb collisions. It is observed in the measurements of the pair nuclear modification factors that the subleading jets are significantly suppressed relative to leading jets with $p_{\mathrm{T}}$ between 100 and 316 GeV for all centralities in Pb+Pb collisions.

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.

Two-particle long-range azimuthal correlations are measured in photonuclear collisions using 1.7 nb$^{-1}$ of 5.02 TeV Pb+Pb collision data collected by the ATLAS experiment at the LHC. Candidate events are selected using a dedicated high-multiplicity photonuclear event trigger, a combination of information from the zero-degree calorimeters and forward calorimeters, and from pseudorapidity gaps constructed using calorimeter energy clusters and charged-particle tracks. Distributions of event properties are compared between data and Monte Carlo simulations of photonuclear processes. Two-particle correlation functions are formed using charged-particle tracks in the selected events, and a template-fitting method is employed to subtract the non-flow contribution to the correlation. Significant nonzero values of the second- and third-order flow coefficients are observed and presented as a function of charged-particle multiplicity and transverse momentum. The results are compared with flow coefficients obtained in proton-proton and proton-lead collisions in similar multiplicity ranges, and with theoretical expectations. The unique initial conditions present in this measurement provide a new way to probe the origin of the collective signatures previously observed only in hadronic collisions.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

Exclusive dimuon production in ultraperipheral collisions (UPC), resulting from photon-photon interactions in the strong electromagnetic fields of colliding high-energy lead nuclei, $\mathrm{PbPb}(\gamma\gamma) \rightarrow \mu^+\mu^- (\mathrm{Pb}^{(\star)}\mathrm{Pb}^{(\star)} )$, is studied using $\mathcal{L}_{\mathrm{int}} = 0.48$ nb$^{-1}$ of $\sqrt{s_\mathrm{NN}}=5.02$ TeV lead-lead collision data at the LHC with the ATLAS detector. Dimuon pairs are measured in the fiducial region $p_{\mathrm{T}\mu} > 4$ GeV, $|\eta_{\mu}| < 2.4$, invariant mass $m_{\mu\mu} > 10$ GeV, and $p_{\mathrm{T,\mu\mu}} < 2$ GeV. The primary background from single-dissociative processes is extracted from the data using a template fitting technique. Differential cross sections are presented as a function of $m_{\mu\mu}$, absolute pair rapidity ($|y_{\mu\mu}|$), scattering angle in the dimuon rest frame ($|\cos \vartheta^{\star}_{\mu\mu}|$) and the colliding photon energies. The total cross section of the UPC $\gamma \gamma \rightarrow \mu^{+}\mu^{-}$ process in the fiducial volume is measured to be $\sigma_{\mathrm{fid}}^{\mu\mu} = 34.1 \! \pm \! 0.3 \mathrm{(stat.)} \! \pm \! 0.7 \mathrm{(syst.)}$ $\mu\mathrm{b}$. Generally good agreement is found with calculations from STARlight, which incorporate the leading-order Breit-Wheeler process with no final-state effects, albeit differences between the measurements and theoretical expectations are observed. In particular, the measured cross sections at larger $|y_{\mu\mu}|$ are found to be about 10-20% larger in data than in the calculations, suggesting the presence of larger fluxes of photons in the initial state. Modification of the dimuon cross sections in the presence of forward and/or backward neutron production is also studied and is found to be associated with a harder incoming photon spectrum, consistent with expectations.

Two-particle pseudorapidity correlations are measured in $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV Pb+Pb, $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV $p$+Pb, and $\sqrt{s}$ = 13 TeV $pp$ collisions at the LHC, with total integrated luminosities of approximately 7 $\mu\mathrm{b}^{-1}$, 28 $\mathrm{nb}^{-1}$, and 65 $\mathrm{nb}^{-1}$, respectively. The correlation function $C_{\rm N}(\eta_1,\eta_2)$ is measured as a function of event multiplicity using charged particles in the pseudorapidity range $|\eta|<2.4$. The correlation function contains a significant short-range component, which is estimated and subtracted. After removal of the short-range component, the shape of the correlation function is described approximately by $1+\langle{a_1^2}\rangle \eta_1\eta_2$ in all collision systems over the full multiplicity range. The values of $\sqrt{\langle{a_1^2}\rangle}$ are consistent between the opposite-charge pairs and same-charge pairs, and for the three collision systems at similar multiplicity. The values of $\sqrt{\langle{a_1^2}\rangle}$ and the magnitude of the short-range component both follow a power-law dependence on the event multiplicity. The $\eta$ distribution of the short-range component, after symmetrizing the proton and lead directions in $p$+Pb collisions, is found to be smaller than that in $pp$ collisions with comparable multiplicity.

Light-by-light scattering ($\gamma\gamma\rightarrow\gamma\gamma$) is a quantum-mechanical process that is forbidden in the classical theory of electrodynamics. This reaction is accessible at the Large Hadron Collider thanks to the large electromagnetic field strengths generated by ultra-relativistic colliding lead (Pb) ions. Using 480 $\mu$b$^{-1}$ of Pb+Pb collision data recorded at a centre-of-mass energy per nucleon pair of 5.02 TeV by the ATLAS detector, the ATLAS Collaboration reports evidence for the $\gamma\gamma\rightarrow\gamma\gamma$ reaction. A total of 13 candidate events are observed with an expected background of 2.6$\pm$0.7 events. After background subtraction and analysis corrections, the fiducial cross section of the process $\textrm{Pb+Pb}\,(\gamma\gamma)\rightarrow \textrm{Pb}^{(\ast)}\textrm{+}\textrm{Pb}^{(\ast)}\,\gamma\gamma$, for photon transverse energy $E_{\mathrm{T}}>$3 GeV, photon absolute pseudorapidity $|\eta|<$2.4, diphoton invariant mass greater than 6 GeV, diphoton transverse momentum lower than 2 GeV and diphoton acoplanarity below 0.01, is measured to be 70 $\pm$ 24 (stat.) $\pm$ 17 (syst.) nb, which is in agreement with Standard Model predictions.

The ATLAS Collaboration has measured the inclusive production of $Z$ bosons via their decays into electron and muon pairs in $p+$Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV at the Large Hadron Collider. The measurements are made using data corresponding to integrated luminosities of 29.4 nb$^{-1}$ and 28.1 nb$^{-1}$ for $Z \rightarrow ee$ and $Z \rightarrow \mu\mu$, respectively. The results from the two channels are consistent and combined to obtain a cross section times the $Z \rightarrow \ell\ell$ branching ratio, integrated over the rapidity region $|y^{*}_{Z}|<3.5$, of 139.8 $\pm$ 4.8 (stat.) $\pm$ 6.2 (syst.) $\pm$ 3.8 (lumi.) nb. Differential cross sections are presented as functions of the $Z$ boson rapidity and transverse momentum, and compared with models based on parton distributions both with and without nuclear corrections. The centrality dependence of $Z$ boson production in $p+$Pb collisions is measured and analyzed within the framework of a standard Glauber model and the model's extension for fluctuations of the underlying nucleon-nucleon scattering cross section.