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 distributions of event-by-event harmonic flow coefficients v_n for n=2-4 are measured in sqrt(s_NN)=2.76 TeV Pb+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using charged particles with transverse momentum pT> 0.5 GeV and in the pseudorapidity range |eta|<2.5 in a dataset of approximately 7 ub^-1 recorded in 2010. The shapes of the v_n distributions are described by a two-dimensional Gaussian function for the underlying flow vector in central collisions for v_2 and over most of the measured centrality range for v_3 and v_4. Significant deviations from this function are observed for v_2 in mid-central and peripheral collisions, and a small deviation is observed for v_3 in mid-central collisions. It is shown that the commonly used multi-particle cumulants are insensitive to the deviations for v_2. The v_n distributions are also measured independently for charged particles with 0.5<pT<1 GeV and pT>1 GeV. When these distributions are rescaled to the same mean values, the adjusted shapes are found to be nearly the same for these two pT ranges. The v_n distributions are compared with the eccentricity distributions from two models for the initial collision geometry: a Glauber model and a model that includes corrections to the initial geometry due to gluon saturation effects. Both models fail to describe the experimental data consistently over most of the measured centrality range.

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.

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.

The integrated elliptic flow of charged particles produced in Pb+Pb collisions at sqrt(s_NN)=2.76 TeV has been measured with the ATLAS detector using data collected at the Large Hadron Collider. The anisotropy parameter, v_2, was measured in the pseudorapidity range |eta| <= 2.5 with the event-plane method. In order to include tracks with very low transverse momentum pT, thus reducing the uncertainty in v_2 integrated over pT, a 1 mu b-1 data sample without a magnetic field in the tracking detectors is used. The centrality dependence of the integrated v_2 is compared to other measurements obtained with higher pT thresholds. A weak pseudorapidity dependence of the integrated elliptic flow is observed for central collisions, and a small decrease when moving away from mid-rapidity is observed only in peripheral collisions. The integrated v2 transformed to the rest frame of one of the colliding nuclei is compared to the lower-energy RHIC data.

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to the physics of jet quenching. In a sample of lead-lead collisions at $\sqrt{s_{NN}}$ = 2.76 TeV corresponding to an integrated luminosity of approximately 7 inverse microbarns, ATLAS has measured jets with a calorimeter over the pseudorapidity interval |$\eta$| < 2.1 and over the transverse momentum range 38 < pT < 210 GeV. Jets were reconstructed using the anti-$k_t$ algorithm with values for the distance parameter that determines the nominal jet radius of R = 0.2, 0.3, 0.4 and 0.5. The centrality dependence of the jet yield is characterized by the jet "central-to-peripheral ratio," $R_{cp}$. Jet production is found to be suppressed by approximately a factor of two in the 10% most central collisions relative to peripheral collisions. $R_{cp}$ varies smoothly with centrality as characterized by the number of participating nucleons. The observed suppression is only weakly dependent on jet radius and transverse momentum. These results provide the first direct measurement of inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse energy imbalance at the LHC.

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.

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.

Measurements of the variation of inclusive jet suppression as a function of relative azimuthal angle, Delta phi, with respect to the elliptic event plane provide insight into the path-length dependence of jet quenching. ATLAS has measured the Delta phi dependence of jet yields in 0.14 nb^-1 of sqrt(s(NN))= 2.76 TeV Pb+Pb collisions at the LHC for jet transverse momenta p_T > 45 GeV in different collision centrality bins using an underlying event subtraction procedure that accounts for elliptic flow. The variation of the jet yield with Delta phi was characterized by the parameter, v_2^jet, and the ratio of out-of-plane (Delta phi ~ pi/2) to in-plane (Delta phi ~ 0) yields. Non-zero v_2^jet values were measured in all centrality bins for p_T < 160 GeV. The jet yields are observed to vary by as much as 20% between in-plane and out-of-plane directions.