The inclusive production rates of isolated, prompt photons in $p$+Pb collisions at $\sqrt{s_\mathrm{NN}} = 8.16$ TeV are studied with the ATLAS detector at the Large Hadron Collider using a dataset with an integrated luminosity of 165 nb$^{-1}$ recorded in 2016. The cross-section and nuclear modification factor $R_{p\mathrm{Pb}}$ are measured as a function of photon transverse energy from 20 GeV to 550 GeV and in three nucleon-nucleon centre-of-mass pseudorapidity regions, (-2.83,-2.02), (-1.84,0.91), and (1.09,1.90). The cross-section and $R_{p\mathrm{Pb}}$ values are compared with the results of a next-to-leading-order perturbative QCD calculation, with and without nuclear parton distribution function modifications, and with expectations based on a model of the energy loss of partons prior to the hard scattering. The data disfavour a large amount of energy loss and provide new constraints on the parton densities in nuclei.

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.

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 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 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.

The ATLAS experiment has observed 1995 Z boson candidates in data corresponding to 0.15 inverse nb of integrated luminosity obtained in the 2011 LHC Pb+Pb run at sqrt(s_NN)=2.76 TeV. The Z bosons are reconstructed via di-electron and di-muon decay channels, with a background contamination of less than 3%. Results from the two channels are consistent and are combined. Within the statistical and systematic uncertainties, the per-event Z boson yield is proportional to the number of binary collisions estimated by the Glauber model. The elliptic anisotropy of the azimuthal distribution of the Z boson with respect to the event plane is found to be consistent with zero.

The distributions of transverse momentum and longitudinal momentum fraction of charged particles in jets are measured in Pb+Pb and pp collisions with the ATLAS detector at the LHC. The distributions are measured as a function of jet transverse momentum and rapidity. The analysis utilises an integrated luminosity of 0.14 nb$^{-1}$ of Pb+Pb data and 4.0 pb$^{-1}$ of pp data collected in 2011 and 2013, respectively, at the same centre-of-mass energy of 2.76 TeV per colliding nucleon pair. The distributions measured in pp collisions are used as a reference for those measured in Pb+Pb collisions in order to evaluate the impact on the internal structure of jets from the jet energy loss of fast partons propagating through the hot, dense medium created in heavy-ion collisions. Modest but significant centrality-dependent modifications of fragmentation functions in Pb+Pb collisions with respect to those in pp collisions are seen. No significant dependence of modifications on jet $p_{\mathrm{T}}$ and rapidity selections is observed except for the fragments with the highest transverse momenta for which some reduction of yields is observed for more forward jets.

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.

Differential measurements of charged particle azimuthal anisotropy are presented for lead-lead collisions at sqrt(s_NN) = 2.76 TeV with the ATLAS detector at the LHC, based on an integrated luminosity of approximately 8 mb^-1. This anisotropy is characterized via a Fourier expansion of the distribution of charged particles in azimuthal angle (phi), with the coefficients v_n denoting the magnitude of the anisotropy. Significant v_2-v_6 values are obtained as a function of transverse momentum (0.5<pT<20 GeV), pseudorapidity (|eta|<2.5) and centrality using an event plane method. The v_n values for n>=3 are found to vary weakly with both eta and centrality, and their pT dependencies are found to follow an approximate scaling relation, v_n^{1/n}(pT) \propto v_2^{1/2}(pT). A Fourier analysis of the charged particle pair distribution in relative azimuthal angle (Dphi=phi_a-phi_b) is performed to extract the coefficients v_{n,n}=<cos (n Dphi)>. For pairs of charged particles with a large pseudorapidity gap (|Deta=eta_a-eta_b|>2) and one particle with pT<3 GeV, the v_{2,2}-v_{6,6} values are found to factorize as v_{n,n}(pT^a,pT^b) ~ v_n(pT^a)v_n(pT^b) in central and mid-central events. Such factorization suggests that these values of v_{2,2}-v_{6,6} are primarily due to the response of the created matter to the fluctuations in the geometry of the initial state. A detailed study shows that the v_{1,1}(pT^a,pT^b) data are consistent with the combined contributions from a rapidity-even v_1 and global momentum conservation. A two-component fit is used to extract the v_1 contribution. The extracted v_1 is observed to cross zero at pT\sim1.0 GeV, reaches a maximum at 4-5 GeV with a value comparable to that for v_3, and decreases at higher pT.

Correlations between the elliptic or triangular flow coefficients $v_m$ ($m$=2 or 3) and other flow harmonics $v_n$ ($n$=2 to 5) are measured using $\sqrt{s_{NN}}=2.76$ TeV Pb+Pb collision data collected in 2010 by the ATLAS experiment at the LHC, corresponding to an integrated lumonisity of 7 $\mu$b$^{-1}$. The $v_m$-$v_n$ correlations are measured in midrapidity as a function of centrality, and, for events within the same centrality interval, as a function of event ellipticity or triangularity defined in a forward rapidity region. For events within the same centrality interval, $v_3$ is found to be anticorrelated with $v_2$ and this anticorrelation is consistent with similar anticorrelations between the corresponding eccentricities $\epsilon_2$ and $\epsilon_3$. On the other hand, it is observed that $v_4$ increases strongly with $v_2$, and $v_5$ increases strongly with both $v_2$ and $v_3$. The trend and strength of the $v_m$-$v_n$ correlations for $n$=4 and 5 are found to disagree with $\epsilon_m$-$\epsilon_n$ correlations predicted by initial-geometry models. Instead, these correlations are found to be consistent with the combined effects of a linear contribution to $v_n$ and a nonlinear term that is a function of $v_2^2$ or of $v_2v_3$, as predicted by hydrodynamic models. A simple two-component fit is used to separate these two contributions. The extracted linear and nonlinear contributions to $v_4$ and $v_5$ are found to be consistent with previously measured event-plane correlations.