Parton energy loss in the quark-gluon plasma (QGP) is studied with a measurement of photon-tagged jet production in 1.7 nb$^{-1}$ of Pb+Pb data and 260 pb$^{-1}$ of $pp$ data, both at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV, with the ATLAS detector. The process $pp \to \gamma$+jet+$X$ and its analogue in Pb+Pb collisions is measured in events containing an isolated photon with transverse momentum ($p_\mathrm{T}$) above $50$ GeV and reported as a function of jet $p_\mathrm{T}$. This selection results in a sample of jets with a steeply falling $p_\mathrm{T}$ distribution that are mostly initiated by the showering of quarks. The $pp$ and Pb+Pb measurements are used to report the nuclear modification factor, $R_\mathrm{AA}$, and the fractional energy loss, $S_\mathrm{loss}$, for photon-tagged jets. In addition, the results are compared with the analogous ones for inclusive jets, which have a significantly smaller quark-initiated fraction. The $R_\mathrm{AA}$ and $S_\mathrm{loss}$ values are found to be significantly different between those for photon-tagged jets and inclusive jets, demonstrating that energy loss in the QGP is sensitive to the colour-charge of the initiating parton. The results are also compared with a variety of theoretical models of colour-charge-dependent energy loss.

The first measurement of the ${\rm e}^{+}{\rm e}^{-}$ pair production at low lepton pair transverse momentum ($p_{\rm T,ee}$) and low invariant mass ($m_{\rm ee}$) in non-central Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV at the LHC is presented. The dielectron production is studied with the ALICE detector at midrapidity ($|\eta_{\rm e}| < 0.8$) as a function of invariant mass ($0.4 \leq m_{\rm ee} < 2.7$ GeV/$c^2$) in the 50$-$70% and 70$-$90% centrality classes for $p_{\rm T,ee} < 0.1$ GeV/$c$, and as a function of $p_{\rm T,ee}$ in three $m_{\rm ee}$ intervals in the most peripheral Pb$-$Pb collisions. Below a $p_{\rm T,ee}$ of 0.1 GeV/$c$, a clear excess of ${\rm e}^{+}{\rm e}^{-}$ pairs is found compared to the expectations from known hadronic sources and predictions of thermal radiation from the medium. The $m_{\rm ee}$ excess spectra are reproduced, within uncertainties, by different predictions of the photon$-$photon production of dielectrons, where the photons originate from the extremely strong electromagnetic fields generated by the highly Lorentz-contracted Pb nuclei. Lowest-order quantum electrodynamic (QED) calculations, as well as a model that takes into account the impact-parameter dependence of the average transverse momentum of the photons, also provide a good description of the $p_{\rm T,ee}$ spectra. The measured $\sqrt{\langle p_{\rm T,ee}^{2} \rangle}$ of the excess $p_{\rm T,ee}$ spectrum in peripheral Pb$-$Pb collisions is found to be comparable to the values observed previously at RHIC in a similar phase-space region.

We report about the properties of the underlying event measured with ALICE at the LHC in pp and p$-$Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. The event activity, quantified by charged-particle number and summed-$p_{\rm T}$ densities, is measured as a function of the leading-particle transverse momentum ($p_{\rm T}^{\rm trig}$). These quantities are studied in three azimuthal-angle regions relative to the leading particle in the event: toward, away, and transverse. Results are presented for three different $p_{\rm T}$ thresholds (0.15, 0.5, and 1 GeV/$c$) at mid-pseudorapidity ($|\eta|<0.8$). The event activity in the transverse region, which is the most sensitive to the underlying event, exhibits similar behaviour in both pp and p$-$Pb collisions, namely, a steep increase with $p_{\rm T}^{\rm trig}$ for low $p_{\rm T}^{\rm trig}$, followed by a saturation at $p_{\rm T}^{\rm trig} \approx 5$ GeV/$c$. The results from pp collisions are compared with existing measurements at other centre-of-mass energies. The quantities in the toward and away regions are also analyzed after the subtraction of the contribution measured in the transverse region. The remaining jet-like particle densities are consistent in pp and p$-$Pb collisions for $p_{\rm T}^{\rm trig}>10$ GeV/$c$, whereas for lower $p_{\rm T}^{\rm trig}$ values the event activity is slightly higher in p$-$Pb than in pp collisions. The measurements are compared with predictions from the PYTHIA 8 and EPOS LHC Monte Carlo event generators.

A measurement of $J/\psi$ and $\psi(2\mathrm{S})$ production is presented. It is based on a data sample from Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV and $pp$ collisions at $\sqrt{s}$ = 5.02 TeV recorded by the ATLAS detector at the LHC in 2015, corresponding to an integrated luminosity of $0.42\mathrm{nb}^{-1}$ and $25\mathrm{pb}^{-1}$ in Pb+Pb and $pp$, respectively. The measurements of per-event yields, nuclear modification factors, and non-prompt fractions are performed in the dimuon decay channel for $9 < p_{T}^{\mu\mu} < 40$ GeV in dimuon transverse momentum, and $-2.0 < y_{\mu\mu} < 2.0$ in rapidity. Strong suppression is found in Pb+Pb collisions for both prompt and non-prompt $J/\psi$, as well as for prompt and non-prompt $\psi(2\mathrm{S})$, increasing with event centrality. The suppression of prompt $\psi(2\mathrm{S})$ is observed to be stronger than that of $J/\psi$, while the suppression of non-prompt $\psi(2\mathrm{S})$ is equal to that of the non-prompt $J/\psi$ within uncertainties, consistent with the expectation that both arise from \textit{b}-quarks propagating through the medium. Despite prompt and non-prompt $J/\psi$ arising from different mechanisms, the dependence of their nuclear modification factors on centrality is found to be quite similar.

Studies of the fragmentation of jets into charged particles in heavy-ion collisions can provide information about the mechanism of jet-quenching by the hot and dense QCD matter created in such collisions, the quark-gluon plasma. This paper presents a measurement of the angular distribution of charged particles around the jet axis in $\sqrt{s_{\mathrm{NN}}}=$ 5.02 TeV Pb+Pb and $pp$ collisions, using the ATLAS detector at the LHC. The Pb+Pb and $pp$ data sets have integrated luminosities of 0.49 nb$^{-1}$ and 25 pb$^{-1}$, respectively. The measurement is performed for jets reconstructed with the anti-$k_{t}$ algorithm with radius parameter $R = 0.4$ and is extended to an angular distance of $r= 0.8$ from the jet axis. Results are presented as a function of Pb+Pb collision centrality and distance from the jet axis for charged particles with transverse momenta in the 1$-$63 GeV range, matched to jets with transverse momenta in the 126$-$316 GeV range and an absolute value of jet rapidity of less than 1.7. Modifications to the measured distributions are quantified by taking a ratio to the measurements in $pp$ collisions. Yields of charged particles with transverse momenta below 4 GeV are observed to be increasingly enhanced as a function of angular distance from the jet axis, reaching a maximum at $r=0.6$. Charged particles with transverse momenta above 4 GeV have an enhanced yield in Pb+Pb collisions in the jet core for angular distances up to $r = 0.05$ from the jet axis, with a suppression at larger distances.

Multi-particle azimuthal cumulants are measured as a function of centrality and transverse momentum using 470 $\mu$b$^{-1}$ of Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ATLAS detector at the LHC. These cumulants provide information on the event-by-event fluctuations of harmonic flow coefficients $v_n$ and correlated fluctuations between two harmonics $v_n$ and $v_m$. For the first time, a non-zero four-particle cumulant is observed for dipolar flow, $v_1$. The four-particle cumulants for elliptic flow, $v_2$, and triangular flow, $v_3$, exhibit a strong centrality dependence and change sign in ultra-central collisions. This sign change is consistent with significant non-Gaussian fluctuations in $v_2$ and $v_3$. The four-particle cumulant for quadrangular flow, $v_4$, is found to change sign in mid-central collisions. Correlations between two harmonics are studied with three- and four-particle mixed-harmonic cumulants, which indicate an anti-correlation between $v_2$ and $v_3$, and a positive correlation between $v_2$ and $v_4$. These correlations decrease in strength towards central collisions and either approach zero or change sign in ultra-central collisions. To investigate the possible flow fluctuations arising from intrinsic centrality or volume fluctuations, the results are compared between two different event classes used for centrality definitions. In peripheral and mid-central collisions where the cumulant signals are large, only small differences are observed. In ultra-central collisions, the differences are much larger and transverse momentum dependent. These results provide new information to disentangle flow fluctuations from the initial and final states, as well as new insights on the influence of centrality fluctuations.

Azimuthal anisotropies of muons from charm and bottom hadron decays are measured in Pb+Pb collisions at $\sqrt{s_\mathrm{NN}}= 5.02$ TeV. The data were collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2018 with integrated luminosities of $0.5~\mathrm{nb}^{-1}$ and $1.4~\mathrm{nb^{-1}}$, respectively. The kinematic selection for heavy-flavor muons requires transverse momentum $4 < p_\mathrm{T} < 30$ GeV and pseudorapidity $|\eta|<2.0$. The dominant sources of muons in this $p_\mathrm{T}$ range are semi-leptonic decays of charm and bottom hadrons. These heavy-flavor muons are separated from light-hadron decay muons and punch-through hadrons using the momentum imbalance between the measurements in the tracking detector and in the muon spectrometers. Azimuthal anisotropies, quantified by flow coefficients, are measured via the event-plane method for inclusive heavy-flavor muons as a function of the muon $p_\mathrm{T}$ and in intervals of Pb+Pb collision centrality. Heavy-flavor muons are separated into contributions from charm and bottom hadron decays using the muon transverse impact parameter with respect to the event primary vertex. Non-zero elliptic ($v_{2}$) and triangular ($v_{3}$) flow coefficients are extracted for charm and bottom muons, with the charm muon coefficients larger than those for bottom muons for all Pb+Pb collision centralities. The results indicate substantial modification to the charm and bottom quark angular distributions through interactions in the quark-gluon plasma produced in these Pb+Pb collisions, with smaller modifications for the bottom quarks as expected theoretically due to their larger mass.

The elliptic flow of muons from the decay of charm and bottom hadrons is measured in $pp$ collisions at $\sqrt{s}=13$ TeV using a data sample with an integrated luminosity of 150 pb$^{-1}$ recorded by the ATLAS detector at the LHC. The muons from heavy-flavor decay are separated from light-hadron decay muons using momentum imbalance between the tracking and muon spectrometers. The heavy-flavor decay muons are further separated into those from charm decay and those from bottom decay using the distance-of-closest-approach to the collision vertex. The measurement is performed for muons in the transverse momentum range 4-7 GeV and pseudorapidity range $|\eta|<2.4$. A significant non-zero elliptic anisotropy coefficient $v_{2}$ is observed for muons from charm decays, while the $v_{2}$ value for muons from bottom decays is consistent with zero within uncertainties.

This paper presents a measurement of forward-forward and forward-central dijet azimuthal angular correlations and conditional yields in proton-proton ($pp$) and proton-lead ($p$+Pb) collisions as a probe of the nuclear gluon density in regions where the fraction of the average momentum per nucleon carried by the parton entering the hard scattering is low. In these regions, gluon saturation can modify the rapidly increasing parton distribution function of the gluon. The analysis utilizes 25 pb$^{-1}$ of $pp$ data and 360 $\mu \mathrm{b}^{-1}$ of $p$+Pb data, both at $\sqrt{s_{\rm NN}}$ = 5.02 TeV, collected in 2015 and 2016, respectively, with the ATLAS detector at the LHC. The measurement is performed in the center-of-mass frame of the nucleon-nucleon system in the rapidity range between $-$4.0 and 4.0 using the two highest transverse momentum jets in each event, with the highest transverse momentum jet restricted to the forward rapidity range. No significant broadening of azimuthal angular correlations is observed for forward-forward or forward-central dijets in $p$+Pb compared to $pp$ collisions. For forward-forward jet pairs in the proton-going direction, the ratio of conditional yields in $p$+Pb collisions to those in $pp$ collisions is suppressed by approximately 20%, with no significant dependence on the transverse momentum of the dijet system. No modification of conditional yields is observed for forward-central dijets.

The azimuthal anisotropy of charged particles produced in $\sqrt{s_{\mathrm{NN}}}=8.16$ TeV $p$+Pb collisions is measured with the ATLAS detector at the LHC. The data correspond to an integrated luminosity of $165$ $\mathrm{nb}^{-1}$ that was collected in 2016. Azimuthal anisotropy coefficients, elliptic $v_2$ and triangular $v_3$, extracted using two-particle correlations with a non-flow template fit procedure, are presented as a function of particle transverse momentum ($p_\mathrm{T}$) between $0.5$ and $50$ GeV. The $v_2$ results are also reported as a function of centrality in three different particle $p_\mathrm{T}$ intervals. The results are reported from minimum-bias events and jet-triggered events, where two jet $p_\mathrm{T}$ thresholds are used. The anisotropies for particles with $p_\mathrm{T}$ less than about $2$ GeV are consistent with hydrodynamic flow expectations, while the significant non-zero anisotropies for $p_\mathrm{T}$ in the range $9$-$50$ GeV are not explained within current theoretical frameworks. In the $p_\mathrm{T}$ range $2$-$9$ GeV, the anisotropies are larger in minimum-bias than in jet-triggered events. Possible origins of these effects, such as the changing admixture of particles from hard scattering and the underlying event, are discussed.