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.

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.

Charge-dependent azimuthal particle correlations with respect to the second-order event plane in pPb and PbPb collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV have been studied with the CMS experiment at the LHC. The measurement is performed with a three-particle correlation technique, using two particles with the same or opposite charge within the pseudorapidity range abs(eta) < 2.4, and a third particle measured in the hadron forward calorimeters (4.4 < abs(eta) < 5). The observed differences between the same and opposite sign correlations, as functions of multiplicity and eta gap between the two charged particles, are of similar magnitude in pPb and PbPb collisions at the same multiplicities. These results pose a challenge for the interpretation of charge-dependent azimuthal correlations in heavy ion collisions in terms of the chiral magnetic effect.

The prompt D$^0$ meson azimuthal anisotropy coefficients, $v_2$ and $v_3$, are measured at midrapidity ($|y| < 1.0$) in PbPb collisions at a center-of-mass energy $\sqrt{s_\mathrm{NN}} = $5.02 TeV per nucleon pair with data collected by the CMS experiment. The measurement is performed in the transverse momentum ($p_\mathrm{T}$) range of 1 to 40 GeV/c, for central and midcentral collisions. The $v_2$ coefficient is found to be positive throughout the $p_\mathrm{T}$ range studied. The first measurement of the prompt D$^0$ meson $v_3$ coefficient is performed, and values up to 0.07 are observed for $p_\mathrm{T}$ around 4 GeV/c. Compared to measurements of charged particles, a similar $p_\mathrm{T}$ dependence, but smaller magnitude for $p_\mathrm{T} < $6 GeV/c, is found for prompt D$^0$ meson $v_2$ and $v_3$ coefficients. The results are consistent with the presence of collective motion of charm quarks at low $p_\mathrm{T}$ and a path length dependence of charm quark energy loss at high $p_\mathrm{T}$, thereby providing new constraints on the theoretical description of the interactions between charm quarks and the quark-gluon plasma.

The nuclear modification factor RAA and the azimuthal anisotropy coefficient v[2] of prompt and nonprompt (i.e. those from decays of b hadrons) J/psi mesons, measured from PbPb and pp collisions at sqrt(s[NN]) = 2.76 TeV at the LHC, are reported. The results are presented in several event centrality intervals and several kinematic regions, for transverse momenta pt > 6.5 GeV/c and rapidity abs(y) < 2.4, extending down to pt = 3 GeV/c in the 1.6 < abs(y) < 2.4 range. The v[2] of prompt J/psi is found to be nonzero, but with no strong dependence on centrality, rapidity, or pt over the full kinematic range studied. The measured v[2] of nonprompt J/psi is consistent with zero. The RAA of prompt J/psi exhibits a suppression that increases from peripheral to central collisions but does not vary as a function of either y or pt in the fiducial range. The nonprompt J/psi RAA shows a suppression which becomes stronger as rapidity or pt increase. The v[2] and nuclear suppression of open and hidden charm, and of open charm and beauty, are compared.

The spectra of charged particles produced within the pseudorapidity window abs(eta) < 1 at sqrt(s[NN]) = 5.02 TeV are measured using 404 inverse microbarns of PbPb and 27.4 inverse picobarns of pp data collected by the CMS detector at the LHC in 2015. The spectra are presented over the transverse momentum ranges spanning 0.5 < pt < 400 GeV in pp and 0.7 < pt < 400 GeV in PbPb collisions. The corresponding nuclear modification factor, R[AA], is measured in bins of collision centrality. The R[AA] in the 5% most central collisions shows a maximal suppression by a factor of 7-8 in the pt region of 6-9 GeV. This dip is followed by an increase, which continues up to the highest pt measured, and approaches unity in the vicinity of pt = 200 GeV. The R[AA] is compared to theoretical predictions and earlier experimental results at lower collision energies. The newly measured pp spectrum is combined with the pPb spectrum previously published by the CMS Collaboration to construct the pPb nuclear modification factor, R[pA], up to 120 GeV. For pt > 20 GeV, R[pA] exhibits weak momentum dependence and shows a moderate enhancement above unity.

The production yields of Upsilon(1S), Upsilon(2S), and Upsilon(3S) quarkonium states are measured through their decays into muon pairs in the CMS detector, in PbPb and pp collisions at the centre-of-mass energy per nucleon pair of 2.76 TeV. The data correspond to integrated luminosities of 166 inverse microbarns and 5.4 inverse picobarns for PbPb and pp collisions, respectively. Differential production cross sections are reported as functions of Upsilon rapidity y up to 2.4, and transverse momentum pT up to 20 GeV/c. A strong centrality-dependent suppression is observed in PbPb relative to pp collisions, by factors of up to approximately 2 and 8, for the Upsilon(1S) and Upsilon(2S) states, respectively. No significant dependence of this suppression is observed as a function of y or pT. The Upsilon(3S) state is not observed in PbPb collisions, which corresponds to a suppression for the centrality-integrated data by at least a factor of approximately 7 at a 95% confidence level. The observed suppression is in agreement with theoretical scenarios modeling the sequential melting of quarkonium states in a quark gluon plasma.

A measurement of the groomed jet mass in PbPb and pp collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV with the CMS detector at the LHC is presented. Jet grooming is a recursive procedure which sequentially removes soft constituents of a jet until a pair of hard subjets is found. The resulting groomed jets can be used to study modifications to the parton shower evolution in the presence of the hot and dense medium created in heavy ion collisions. Predictions of groomed jet properties from the PYTHIA and HERWIG++ event generators agree with the measurements in pp collisions. When comparing the results from the most central PbPb collisions to pp data, a hint of an increase of jets with large jet mass is observed, which could originate from additional medium-induced radiation at a large angle from the jet axis. However, no modification of the groomed mass of the core of the jet is observed for all PbPb centrality classes. The PbPb results are also compared to predictions from the JEWEL and Q-PYTHIA event generators, which predict a large modification of the groomed mass not observed in the data.

The first observation of top quark production in proton-nucleus collisions is reported using proton-lead data collected by the CMS experiment at the CERN LHC at a nucleon-nucleon center-of-mass energy of sqrt(s[NN]) = 8.16 TeV. The measurement is performed using events with exactly one isolated electron or muon and at least four jets. The data sample corresponds to an integrated luminosity of 174 inverse nanobarns. The significance of the tt-bar signal against the background-only hypothesis is above five standard deviations. The measured cross section is sigma[tt-bar] = 45 +/- 8 nb, consistent with predictions from perturbative quantum chromodynamics.

The second- and third-order azimuthal anisotropy Fourier harmonics of charged particles produced in pPb collisions, at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV, are studied over a wide range of event multiplicities. Multiparticle correlations are used to isolate global properties stemming from the collision overlap geometry. The second-order "elliptic" harmonic moment is obtained with high precision through four-, six-, and eight-particle correlations and, for the first time, the third-order "triangular" harmonic moment is studied using four-particle correlations. A sample of peripheral PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV that covers a similar range of event multiplicities as the pPb results is also analyzed. Model calculations of initial-state fluctuations in pPb and PbPb collisions can be directly compared to the high precision experimental results. This work provides new insight into the fluctuation-driven origin of the $v_3$ coefficients in pPb and PbPb collisions, and into the dominating overall collision geometry in PbPb collisions at the earliest stages of heavy ion interactions.