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.

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.

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.

We present a measurement of two-particle angular correlations in proton-proton collisions at sqrt(s) = 900 GeV and 7 TeV. The collision events were collected during 2009 and 2010 with the ATLAS detector at the Large Hadron Collider using a single-arm minimum bias trigger. Correlations are measured for charged particles produced in the kinematic range of transverse momentum pT > 100 MeV and pseudorapidity |eta| < 2.5. A complex structure in pseudorapidity and azimuth is observed at both collision energies. Results are compared to Pythia 8 and Herwig++ as well as to the AMBT2B, DW and Perugia 2011 tunes of Pythia 6. The data are not satisfactorily described by any of these models.

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.

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.

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.

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.

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.

In order to study further the long-range correlations ("ridge") observed recently in p+Pb collisions at sqrt(s_NN) =5.02 TeV, the second-order azimuthal anisotropy parameter of charged particles, v_2, has been measured with the cumulant method using the ATLAS detector at the LHC. In a data sample corresponding to an integrated luminosity of approximately 1 microb^(-1), the parameter v_2 has been obtained using two- and four-particle cumulants over the pseudorapidity range |eta|<2.5. The results are presented as a function of transverse momentum and the event activity, defined in terms of the transverse energy summed over 3.1<eta<4.9 in the direction of the Pb beam. They show features characteristic of collective anisotropic flow, similar to that observed in Pb+Pb collisions. A comparison is made to results obtained using two-particle correlation methods, and to predictions from hydrodynamic models of p+Pb collisions. Despite the small transverse spatial extent of the p+Pb collision system, the large magnitude of v_2 and its similarity to hydrodynamic predictions provide additional evidence for the importance of final-state effects in p+Pb reactions.

The paper presents studies of Bose-Einstein Correlations (BEC) for pairs of like-sign charged particles measured in the kinematic range $p_{\rm T}>$ 100 MeV and $|\eta|<$ 2.5 in proton--proton collisions at centre-of-mass energies of 0.9 and 7 TeV with the ATLAS detector at the CERN Large Hadron Collider. The integrated luminosities are approximately 7 $\mu$b$^{-1}$, 190 $\mu$b$^{-1}$ and 12.4 nb$^{-1}$ for 0.9 TeV, 7 TeV minimum-bias and 7 TeV high-multiplicity data samples, respectively. The multiplicity dependence of the BEC parameters characterizing the correlation strength and the correlation source size are investigated for charged-particle multiplicities of up to 240. A saturation effect in the multiplicity dependence of the correlation source size is observed using the high-multiplicity 7 TeV data sample. The dependence of the BEC parameters on the average transverse momentum of the particle pair is also investigated.

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.

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.

The azimuthal variation of jet yields in heavy-ion collisions provides information about the path-length dependence of the energy loss experienced by partons passing through the hot, dense nuclear matter known as the quark-gluon plasma. This paper presents the azimuthal anisotropy coefficients $v_2$, $v_3$, and $v_4$ measured for jets in Pb+Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV using the ATLAS detector at the LHC. The measurement uses data collected in 2015 and 2018, corresponding to an integrated luminosity of 2.2 nb$^{-1}$. The $v_n$ values are measured as a function of the transverse momentum of the jets between 71 GeV and 398 GeV and the event centrality. A nonzero value of $v_2$ is observed in all but the most central collisions. The value of $v_2$ is largest for jets with lower transverse momentum, with values up to 0.05 in mid-central collisions. A smaller, nonzero value of $v_3$ of approximately 0.01 is measured with no significant dependence on jet $p_T$ or centrality, suggesting that fluctuations in the initial state play a small but distinct role in jet energy loss. No significant deviation of $v_4$ from zero is observed in the measured kinematic region.

Measurements of the centrality and rapidity dependence of inclusive jet production in $\sqrt{s_\mathrm{NN}} = 5.02$ TeV proton--lead ($p$+Pb) collisions and the jet cross-section in $\sqrt{s} = 2.76$ TeV proton--proton collisions are presented. These quantities are measured in datasets corresponding to an integrated luminosity of 27.8 nb$^{-1}$ and 4.0 pb$^{-1}$, respectively, recorded with the ATLAS detector at the Large Hadron Collider in 2013. The $p$+Pb collision centrality was characterised using the total transverse energy measured in the pseudorapidity interval $-4.9 < \eta < -3.2$ in the direction of the lead beam. Results are presented for the double-differential per-collision yields as a function of jet rapidity and transverse momentum ($p_\mathrm{T}$) for minimum-bias and centrality-selected $p$+Pb collisions, and are compared to the jet rate from the geometric expectation. The total jet yield in minimum-bias events is slightly enhanced above the expectation in a $p_\mathrm{T}$-dependent manner but is consistent with the expectation within uncertainties. The ratios of jet spectra from different centrality selections show a strong modification of jet production at all $p_\mathrm{T}$ at forward rapidities and for large $p_\mathrm{T}$ at mid-rapidity, which manifests as a suppression of the jet yield in central events and an enhancement in peripheral events. These effects imply that the factorisation between hard and soft processes is violated at an unexpected level in proton-nucleus collisions. Furthermore, the modifications at forward rapidities are found to be a function of the total jet energy only, implying that the violations may have a simple dependence on the hard parton-parton kinematics.

Two-particle correlations in relative azimuthal angle ($\Delta-\phi$) and pseudorapidity ($\Delta-\eta$) are measured in $\sqrt{s_{NN}}$ = 5.02 TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1 $\mu b^{-1}$ of data as a function of $p_T$ and the transverse energy ($\sum E_T^{Pb}$) summed over 3.1 < $\eta$ < 4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|$\Delta-\eta$|<5) near-side ($\Delta-\phi$ ~ 0) correlation that grows rapidly with increasing $\sum E_T^{Pb}$. A long-range away-side ($\Delta-\phi$ ~ pi) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small $\sum E_T^{Pb}$, is found to match the near-side correlation in magnitude, shape (in $\Delta-\eta$ and $\Delta-\phi$) and $\sum E_T^{Pb}$ dependence. The resultant $\Delta-\phi$ correlation is approximately symmetric about $\pi$/2, and is consistent with a cos(2$\Delta-\phi$) modulation for all $\sum E_T^{Pb}$ ranges and particle $p_T$. The amplitude of this modulation is comparable in magnitude and $p_T$ dependence to similar modulations observed in heavy-ion collisions, suggestive of final-state collective effects in high multiplicity events.

Measurements of the azimuthal anisotropy in lead-lead collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV are presented using a data sample corresponding to 0.49 $\mathrm{nb}^{-1}$ integrated luminosity collected by the ATLAS experiment at the LHC in 2015. The recorded minimum-bias sample is enhanced by triggers for "ultra-central" collisions, providing an opportunity to perform detailed study of flow harmonics in the regime where the initial state is dominated by fluctuations. The anisotropy of the charged-particle azimuthal angle distributions is characterized by the Fourier coefficients, $v_{2}-v_{7}$, which are measured using the two-particle correlation, scalar-product and event-plane methods. The goal of the paper is to provide measurements of the differential as well as integrated flow harmonics $v_{n}$ over wide ranges of the transverse momentum, 0.5 $

Measurements of two-particle correlation functions and the first five azimuthal harmonics, $v_1$ to $v_5$, are presented, using 28 $\mathrm{nb}^{-1}$ of $p$+Pb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV measured with the ATLAS detector at the LHC. Significant long-range "ridge-like" correlations are observed for pairs with small relative azimuthal angle ($|\Delta\phi|<\pi/3$) and back-to-back pairs ($|\Delta\phi| > 2\pi/3$) over the transverse momentum range $0.4 < p_{\rm T} < 12$ GeV and in different intervals of event activity. The event activity is defined by either the number of reconstructed tracks or the total transverse energy on the Pb-fragmentation side. The azimuthal structure of such long-range correlations is Fourier decomposed to obtain the harmonics $v_n$ as a function of $p_{\rm T}$ and event activity. The extracted $v_n$ values for $n=2$ to 5 decrease with $n$. The $v_2$ and $v_3$ values are found to be positive in the measured $p_{\rm T}$ range. The $v_1$ is also measured as a function of $p_{\rm T}$ and is observed to change sign around $p_{\rm T}\approx 1.5$-2.0 GeV and then increase to about 0.1 for $p_{\rm T}>4$ GeV. The $v_2(p_{\rm T})$, $v_3(p_{\rm T})$ and $v_4(p_{\rm T})$ are compared to the $v_n$ coefficients in Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}} =2.76$ TeV with similar event multiplicities. Reasonable agreement is observed after accounting for the difference in the average $p_{\rm T}$ of particles produced in the two collision systems.

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.

A detailed study of multi-particle azimuthal correlations is presented using $pp$ data at $\sqrt{s}=5.02$ and 13 TeV, and $p$+Pb data at $\sqrt{s_{\rm{NN}}}=5.02$ TeV, recorded with the ATLAS detector at the LHC. The azimuthal correlations are probed using four-particle cumulants $c_{n}\{4\}$ and flow coefficients $v_n\{4\}=(-c_{n}\{4\})^{1/4}$ for $n=2$ and 3, with the goal of extracting long-range multi-particle azimuthal correlation signals and suppressing the short-range correlations. The values of $c_{n}\{4\}$ are obtained as a function of the average number of charged particles per event, $\left\langle N_{\rm{ch}} \right\rangle$, using the recently proposed two-subevent and three-subevent cumulant methods, and compared with results obtained with the standard cumulant method. The three-subevent method is found to be least sensitive to short-range correlations, which originate mostly from jets with a positive contribution to $c_{n}\{4\}$. The three-subevent method gives a negative $c_{2}\{4\}$, and therefore a well-defined $v_2\{4\}$, nearly independent of $\left\langle N_{\rm{ch}} \right\rangle$, which provides direct evidence that the long-range multi-particle azimuthal correlations persist to events with low multiplicity. Furthermore, $v_2\{4\}$ is found to be smaller than the $v_2\{2\}$ measured using the two-particle correlation method, as expected for long-range collective behavior. Finally, the measured values of $v_2\{4\}$ and $v_2\{2\}$ are used to estimate the number of sources relevant for the initial eccentricity in the collision geometry.

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.

Heavy-flavour hadron production provides information about the transport properties and microscopic structure of the quark-gluon plasma created in ultra-relativistic heavy-ion collisions. A measurement of the muons from semileptonic decays of charm and bottom hadrons produced in Pb+Pb and $pp$ collisions at a nucleon-nucleon centre-of-mass energy of 5.02 TeV with the ATLAS detector at the Large Hadron Collider is presented. The Pb+Pb data were collected in 2015 and 2018 with sampled integrated luminosities of $208~\mathrm{\mu b}^{-1}$ and $38~\mathrm{\mu b^{-1}}$, respectively, and $pp$ data with a sampled integrated luminosity of $1.17~\mathrm{pb}^{-1}$ were collected in 2017. Muons from heavy-flavour semileptonic decays are separated from the light-flavour hadronic background using the momentum imbalance between the inner detector and muon spectrometer measurements, and muons originating from charm and bottom decays are further separated via the muon track's transverse impact parameter. Differential yields in Pb+Pb collisions and differential cross sections in $pp$ collisions for such muons are measured as a function of muon transverse momentum from 4 GeV to 30 GeV in the absolute pseudorapidity interval $|\eta| < 2$. Nuclear modification factors for charm and bottom muons are presented as a function of muon transverse momentum in intervals of Pb+Pb collision centrality. The measured nuclear modification factors quantify a significant suppression of the yields of muons from decays of charm and bottom hadrons, with stronger effects for muons from charm hadron decays.

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.

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.

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.

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.

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.

The inclusive jet cross-section has been measured in proton-proton collisions at sqrt(s)=2.76 TeV in a dataset corresponding to an integrated luminosity of 0.20pb-1 collected with the ATLAS detector at the Large Hadron Collider in 2011. Jets are identified using the anti-kt algorithm with two radius parameters of 0.4 and 0.6. The inclusive jet double-differential cross-section is presented as a function of the jet transverse momentum pT and jet rapidity y, covering a range of 20 <= pT < 430 GeV and |y| < 4.4. The ratio of the cross-section to the inclusive jet cross-section measurement at sqrt(s)=7 TeV, published by the ATLAS Collaboration, is calculated as a function of both transverse momentum and the dimensionless quantity xT = 2 pT / sqrt(s), in bins of jet rapidity. The systematic uncertainties on the ratios are significantly reduced due to the cancellation of correlated uncertainties in the two measurements. Results are compared to the prediction from next-to-leading order perturbative QCD calculations corrected for non-perturbative effects, and next-to-leading order Monte Carlo simulation. Furthermore, the ATLAS jet cross-section measurements at sqrt(s)=2.76 TeV and sqrt(s)=7 TeV are analysed within a framework of next-to-leading order perturbative QCD calculations to determine parton distribution functions of the proton, taking into account the correlations between the measurements.