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.

A measurement of $W^\pm$ boson production in lead-lead collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV is reported using data recorded by the ATLAS experiment at the LHC in 2015, corresponding to a total integrated luminosity of $0.49\;\mathrm{nb^{-1}}$. The $W^\pm$ bosons are reconstructed in the electron or muon leptonic decay channels. Production yields of leptonically decaying $W^\pm$ bosons, normalised by the total number of minimum-bias events and the nuclear thickness function, are measured within a fiducial region defined by the detector acceptance and the main kinematic requirements. These normalised yields are measured separately for $W^+$ and $W^-$ bosons, and are presented as a function of the absolute value of pseudorapidity of the charged lepton and of the collision centrality. The lepton charge asymmetry is also measured as a function of the absolute value of lepton pseudorapidity. In addition, nuclear modification factors are calculated using the $W^\pm$ boson production cross-sections measured in $pp$ collisions. The results are compared with predictions based on next-to-leading-order calculations with CT14 parton distribution functions as well as with predictions obtained with the EPPS16 and nCTEQ15 nuclear parton distribution functions. No dependence of normalised production yields on centrality and a good agreement with predictions are observed for mid-central and central collisions. For peripheral collisions, the data agree with predictions within 1.7 (0.9) standard deviations for $W^-$ ($W^+$) bosons.

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.

Jets created in association with a photon can be used as a calibrated probe to study energy loss in the medium created in nuclear collisions. Measurements of the transverse momentum balance between isolated photons and inclusive jets are presented using integrated luminosities of 0.49 nb$^{-1}$ of Pb+Pb collision data at $\sqrt{s_\mathrm{NN}}=5.02$ TeV and 25 pb$^{-1}$ of $pp$ collision data at $\sqrt{s}=5.02$ TeV recorded with the ATLAS detector at the LHC. Photons with transverse momentum $63.1 < p_\mathrm{T}^{\gamma} < 200$ GeV and $\left|\eta^{\gamma}\right| < 2.37$ are paired inclusively with all jets in the event that have $p_\mathrm{T}^\mathrm{jet} > 31.6$ GeV and pseudorapidity $\left|\eta^\mathrm{jet}\right| < 2.8$. The transverse momentum balance given by the jet-to-photon $p_\mathrm{T}$ ratio, $x_\mathrm{J\gamma}$, is measured for pairs with azimuthal opening angle $\Delta\phi > 7\pi/8$. Distributions of the per-photon jet yield as a function of $x_\mathrm{J\gamma}$, $(1/N_\gamma)(\mathrm{d}N/\mathrm{d}x_\mathrm{J\gamma})$, are corrected for detector effects via a two-dimensional unfolding procedure and reported at the particle level. In $pp$ collisions, the distributions are well described by Monte Carlo event generators. In Pb+Pb collisions, the $x_\mathrm{J\gamma}$ distribution is modified from that observed in $pp$ collisions with increasing centrality, consistent with the picture of parton energy loss in the hot nuclear medium. The data are compared with a suite of energy-loss models and calculations.

The modification of the production of $J/\psi$, $\psi(\mathrm{2S})$, and $\mit{\Upsilon}(n\mathrm{S})$ ($n = 1, 2, 3$) in $p$+Pb collisions with respect to their production in $pp$ collisions has been studied. The $p$+Pb and $pp$ datasets used in this paper correspond to integrated luminosities of $28$ $\mathrm{nb}^{-1}$ and $25$ $\mathrm{pb}^{-1}$ respectively, collected in 2013 and 2015 by the ATLAS detector at the LHC, both at a centre-of-mass energy per nucleon pair of 5.02 TeV. The quarkonium states are reconstructed in the dimuon decay channel. The yields of $J/\psi$ and $\psi(\mathrm{2S})$ are separated into prompt and non-prompt sources. The measured quarkonium differential cross sections are presented as a function of rapidity and transverse momentum, as is the nuclear modification factor, $R_{p\mathrm{Pb}}$ for $J/\psi$ and $\mit{\Upsilon}(\mathrm{1S})$. No significant modification of the $J/\psi$ production is observed while $\mit{\Upsilon}(\mathrm{1S})$ production is found to be suppressed at low transverse momentum in $p$+Pb collisions relative to $pp$ collisions. The production of excited charmonium and bottomonium states is found to be suppressed relative to that of the ground states in central $p$+Pb collisions.

To assess the properties of the quark-gluon plasma formed in heavy-ion collisions, the ATLAS experiment at the LHC measures a correlation between the mean transverse momentum and the magnitudes of the flow harmonics. The analysis uses data samples of lead-lead and proton-lead collisions obtained at the centre-of-mass energy per nucleon pair of 5.02 TeV, corresponding to total integrated luminosities of $22 ~\mu b^{-1}$ and $28~nb^{-1}$, respectively. The measurement is performed using a modified Pearson correlation coefficient with the charged-particle tracks on an event-by-event basis. The modified Pearson correlation coefficients for the $2^{nd}$-, 3$^{rd}$-, and 4$^{th}$-order harmonics are measured as a function of event centrality quantified as the number of charged particles or the number of nucleons participating in the collision. The measurements are performed for several intervals of the charged-particle transverse momentum. The correlation coefficients for all studied harmonics exhibit a strong centrality evolution in the lead-lead collisions, which only weakly depends on the charged-particle momentum range. In the proton-lead collisions, the modified Pearson correlation coefficient measured for the second harmonics shows only weak centrality dependence. The data is qualitatively described by the predictions based on the hydrodynamical model.

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.

Multi-particle cumulants and corresponding Fourier harmonics are measured for azimuthal angle distributions of charged particles in $pp$ collisions at $\sqrt{s}$ = 5.02 and 13 TeV and in $p$+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, and compared to the results obtained for low-multiplicity Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. These measurements aim to assess the collective nature of particle production. The measurements of multi-particle cumulants confirm the evidence for collective phenomena in $p$+Pb and low-multiplicity Pb+Pb collisions. On the other hand, the $pp$ results for four-particle cumulants do not demonstrate collective behaviour, indicating that they may be biased by contributions from non-flow correlations. A comparison of multi-particle cumulants and derived Fourier harmonics across different collision systems is presented as a function of the charged-particle multiplicity. For a given multiplicity, the measured Fourier harmonics are largest in Pb+Pb, smaller in $p$+Pb and smallest in $pp$ collisions. The $pp$ results show no dependence on the collision energy, nor on the multiplicity.

This paper presents a measurement of jet fragmentation functions in 0.49 nb$^{-1}$ of Pb+Pb collisions and 25 pb$^{-1}$ of $pp$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV collected in 2015 with the ATLAS detector at the LHC. These measurements provide insight into the jet quenching process in the quark-gluon plasma created in the aftermath of ultra-relativistic collisions between two nuclei. The modifications to the jet fragmentation functions are quantified by dividing the measurements in Pb+Pb collisions by baseline measurements in $pp$ collisions. This ratio is studied as a function of the transverse momentum of the jet, the jet rapidity, and the centrality of the collision. In both collision systems, the jet fragmentation functions are measured for jets with transverse momentum between 126 GeV and 398 GeV and with an absolute value of jet rapidity less than 2.1. An enhancement of particles carrying a small fraction of the jet momentum is observed, which increases with centrality and with increasing jet transverse momentum. Yields of particles carrying a very large fraction of the jet momentum are also observed to be enhanced. Between these two enhancements of the fragmentation functions a suppression of particles carrying an intermediate fraction of the jet momentum is observed in Pb+Pb collisions. A small dependence of the modifications on jet rapidity is observed.

Correlations of two flow harmonics $v_n$ and $v_m$ via three- and four-particle cumulants are measured in 13 TeV $pp$, 5.02 TeV $p$+Pb, and 2.76 TeV peripheral Pb+Pb collisions with the ATLAS detector at the LHC. The goal is to understand the multi-particle nature of the long-range collective phenomenon in these collision systems. The large non-flow background from dijet production present in the standard cumulant method is suppressed using a method of subevent cumulants involving two, three and four subevents separated in pseudorapidity. The results show a negative correlation between $v_2$ and $v_3$ and a positive correlation between $v_2$ and $v_4$ for all collision systems and over the full multiplicity range. However, the magnitudes of the correlations are found to depend strongly on the event multiplicity, the choice of transverse momentum range and collision system. The relative correlation strength, obtained by normalisation of the cumulants with the $\langle v_n^2\rangle$ from a two-particle correlation analysis, is similar in the three collision systems and depends weakly on the event multiplicity and transverse momentum. These results based on the subevent methods provide strong evidence of a similar long-range multi-particle collectivity in $pp$, $p$+Pb and peripheral Pb+Pb collisions.