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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

Measurements of anisotropic flow coefficients with two- and multi-particle cumulants for inclusive charged particles in Pb-Pb collisions at $\sqrt{{\textit s}_\text{NN}} = 5.02$ and 2.76 TeV are reported in the pseudorapidity range $|\eta| < 0.8$ and transverse momentum $0.2 < p_\text{T} < 50$ GeV/$c$. The full data sample collected by the ALICE detector in 2015 (2010), corresponding to an integrated luminosity of 12.7 (2.0) $\mu$b$^{-1}$ in the centrality range 0-80%, is analysed. Flow coefficients up to the sixth flow harmonic ($v_6$) are reported and a detailed comparison among results at the two energies is carried out. The $p_\text{T}$ dependence of anisotropic flow coefficients and its evolution with respect to centrality and harmonic number $n$ are investigated. An approximate power-law scaling of the form $v_n(p_\text{T}) \sim p_\text{T}^{n/3}$ is observed for all flow harmonics at low $p_\text{T}$ ($0.2 < p_\text{T} < 3$ GeV/$c$). At the same time, the ratios $v_n/v_m^{n/m}$ are observed to be essentially independent of $p_\text{T}$ for most centralities up to about $p_\text{T} = 10$ GeV/$c$. Analysing the differences among higher-order cumulants of elliptic flow ($v_2$), which have different sensitivities to flow fluctuations, a measurement of the standardised skewness of the event-by-event $v_2$ distribution $P(v_2)$ is reported and constraints on its higher moments are provided. The Elliptic Power distribution is used to parametrise $P(v_2)$, extracting its parameters from fits to cumulants. The measurements are compared to different model predictions in order to discriminate among initial-state models and to constrain the temperature dependence of the shear viscosity to entropy-density ratio.

The B$^0_s$ and B$^+$ production yields are measured in PbPb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data sample, collected with the CMS detector at the LHC, corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The mesons are reconstructed in the exclusive decay channels B$^0_s$$\to$ J/$\psi(\mu^+\mu^-)\phi($K$^+$K$^-)$ and B$^+$$\to$ J/$\psi(\mu^+\mu^-)$K$^+$, in the transverse momentum range 7-50 GeV/c and absolute rapidity 0-2.4. The B$^0_s$ meson is observed with a statistical significance in excess of five standard deviations for the first time in nucleus-nucleus collisions. The measurements are performed as functions of the transverse momentum of the B mesons and of the PbPb collision centrality. The ratio of production yields of B$^0_s$ and B$^+$ is measured and compared to theoretical models that include quark recombination effects.

The production of prompt $\Lambda^+_\mathrm{c}$ baryons is measured via the exclusive decay channel $\Lambda^+_\mathrm{c}$$\to$ pK$^-\pi^+$ at a center-of-mass energy per nucleon pair of 5.02 TeV, using proton-proton (pp) and lead-lead (PbPb) collision data collected by the CMS experiment at the CERN LHC. The pp and PbPb data were obtained in 2017 and 2018 with integrated luminosities of 252 and 0.607 nb$^{-1}$, respectively. The measurements are performed within the $\Lambda^+_\mathrm{c}$ rapidity interval $\vert y \vert$$\lt$ 1 with transverse momentum ($p_\mathrm{T}$) ranges of 3-30 and 6-40 GeV/$c$ for pp and PbPb collisions, respectively. Compared to the yields in pp collisions scaled by the expected number of nucleon-nucleon interactions, the observed yields of $\Lambda^+_\mathrm{c}$ with $p_\mathrm{T}$$\gt$ 10 GeV/$c$ are strongly suppressed in PbPb collisions. The level of suppression depends significantly on the collision centrality. The $\Lambda^+_\mathrm{c}$ / D$^0$ production ratio is similar in PbPb and pp collisions at $p_\mathrm{T}$$\gt$ 10 GeV/$c$, suggesting that the coalescence process does not play a dominant role in prompt $\Lambda^+_\mathrm{c}$ baryon production at higher $p_\mathrm{T}$.

The first measurement of the ${\rm e}^{+}{\rm e}^{-}$ pair production at low lepton pair transverse momentum ($p_{\rm T,ee}$) and low invariant mass ($m_{\rm ee}$) in non-central Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV at the LHC is presented. The dielectron production is studied with the ALICE detector at midrapidity ($|\eta_{\rm e}| < 0.8$) as a function of invariant mass ($0.4 \leq m_{\rm ee} < 2.7$ GeV/$c^2$) in the 50$-$70% and 70$-$90% centrality classes for $p_{\rm T,ee} < 0.1$ GeV/$c$, and as a function of $p_{\rm T,ee}$ in three $m_{\rm ee}$ intervals in the most peripheral Pb$-$Pb collisions. Below a $p_{\rm T,ee}$ of 0.1 GeV/$c$, a clear excess of ${\rm e}^{+}{\rm e}^{-}$ pairs is found compared to the expectations from known hadronic sources and predictions of thermal radiation from the medium. The $m_{\rm ee}$ excess spectra are reproduced, within uncertainties, by different predictions of the photon$-$photon production of dielectrons, where the photons originate from the extremely strong electromagnetic fields generated by the highly Lorentz-contracted Pb nuclei. Lowest-order quantum electrodynamic (QED) calculations, as well as a model that takes into account the impact-parameter dependence of the average transverse momentum of the photons, also provide a good description of the $p_{\rm T,ee}$ spectra. The measured $\sqrt{\langle p_{\rm T,ee}^{2} \rangle}$ of the excess $p_{\rm T,ee}$ spectrum in peripheral Pb$-$Pb collisions is found to be comparable to the values observed previously at RHIC in a similar phase-space region.

Quasireal photons exchanged in relativistic heavy ion interactions are powerful probes of the gluonic structure of nuclei. The coherent J/$\psi$ photoproduction cross section in ultraperipheral lead-lead collisions is measured as a function of photon-nucleus center-of-mass energies per nucleon (W$^\text{Pb}_{\gamma\text{N}}$), over a wide range of 40 $\lt$ W$^\text{Pb}_{\gamma\text{N}}$$\lt$ 400 GeV. Results are obtained using data at the nucleon-nucleon center-of-mass energy of 5.02 TeV collected by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 1.52 nb$^{-1}$. The cross section is observed to rise rapidly at low W$^\text{Pb}_{\gamma\text{N}}$, and plateau above W$^\text{Pb}_{\gamma\text{N}}$$\approx$ 40 GeV, up to 400 GeV, a new regime of small Bjorken-$x$ ($\approx$ 6 $\times$ 10$^{-5}$) gluons being probed in a heavy nucleus. The observed energy dependence is not predicted by current quantum chromodynamic models.

The second-order Fourier coefficients ($v_2$) characterizing the azimuthal distribution of $\Upsilon$(1S) and $\Upsilon$(2S) mesons arising from PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV are studied. The $\Upsilon$ mesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The data set corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The scalar product method is used to extract the $v_2$ coefficients of the azimuthal distribution. Results are reported for the rapidity range $|y|$$\lt$ 2.4, with the transverse momentum 0 $\lt$$p_\mathrm{T}$$\lt$ 50 GeV/$c$, and in three centrality ranges of 10-30%, 30-50% and 50-90%. In contrast to the J/$\psi$ mesons, the measured $v_2$ values for the $\Upsilon$ mesons are found to be consistent with zero.

We report the measured transverse momentum ($p_{\rm T}$) spectra of primary charged particles from pp, p-Pb and Pb-Pb collisions at a center-of-mass energy $\sqrt{s_{\rm NN}} = 5.02$ TeV in the kinematic range of $0.15<p_{\rm T}<50$ GeV/$c$ and $|\eta|< 0.8$. A significant improvement of systematic uncertainties motivated the reanalysis of data in pp and Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV, as well as in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, which is also presented. Spectra from Pb-Pb collisions are presented in nine centrality intervals and are compared to a reference spectrum from pp collisions scaled by the number of binary nucleon-nucleon collisions. For central collisions, the $p_{\rm T}$ spectra are suppressed by more than a factor of 7 around 6-7 GeV/$c$ with a significant reduction in suppression towards higher momenta up to 30 GeV/$c$. The nuclear modification factor $R_{\rm pPb}$, constructed from the pp and p-Pb spectra measured at the same collision energy, is consistent with unity above 8 GeV/$c$. While the spectra in both pp and Pb-Pb collisions are substantially harder at $\sqrt{s_{\rm NN}} = 5.02$ TeV compared to 2.76 TeV, the nuclear modification factors show no significant collision energy dependence. The obtained results should provide further constraints on the parton energy loss calculations to determine the transport properties of the hot and dense QCD matter.

The production of $\Upsilon$(2S) and $\Upsilon$(3S) mesons in lead-lead (PbPb) and proton-proton (pp) collisions is studied in their dimuon decay channel using the CMS detector at the LHC. The $\Upsilon$(3S) meson is observed for the first time in PbPb collisions, with a significance above five standard deviations. The ratios of yields measured in PbPb and pp collisions are reported for both the $\Upsilon$(2S) and $\Upsilon$(3S) mesons, as functions of transverse momentum and PbPb collision centrality. These ratios, when appropriately scaled, are significantly less than unity, indicating a suppression of $\Upsilon$ yields in PbPb collisions. This suppression increases from peripheral to central PbPb collisions. Furthermore, the suppression is stronger for $\Upsilon$(3S) mesons compared to $\Upsilon$(2S) mesons, extending the pattern of sequential suppression of quarkonium states in nuclear collisions previously seen for the $\psi$/J, $\psi$(2S), $\Upsilon$(1S), and $\Upsilon$(2S) mesons.

Jets containing a prompt J$/\psi$ meson are studied in lead-lead collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV, using the CMS detector at the LHC. Jets are selected to be in the transverse momentum range of 30 $\lt$$p_\mathrm{T}$$\lt$ 40 GeV. The J$/\psi$ yield in these jets is evaluated as a function of the jet fragmentation variable $z$, the ratio of the J$/\psi$$p_\mathrm{T} $ to the jet $p_\mathrm{T}$. The nuclear modification factor, $R_\mathrm{AA}$, is then derived by comparing the yield in lead-lead collisions to the corresponding expectation based on proton-proton data, at the same nucleon-nucleon center-of-mass energy. The suppression of the J$/\psi$ yield shows a dependence on $z$, indicating that the interaction of the J$/\psi$ with the quark-gluon plasma formed in heavy ion collisions depends on the fragmentation that gives rise to the J$/\psi$ meson.

The production of Z$^0$ bosons at large rapidities in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV is reported. Z$^0$ candidates are reconstructed in the dimuon decay channel (${\rm Z}^0 \rightarrow \mu^+\mu^-$), based on muons selected with pseudo-rapidity $-4.0<\eta<-2.5$ and $p_{\rm T}>20$ GeV/$c$. The invariant yield and the nuclear modification factor, $R_{\rm AA}$, are presented as a function of rapidity and collision centrality. The value of $R_{\rm AA}$ for the 0-20% central Pb-Pb collisions is $0.67 \pm 0.11 \, \mbox{(stat.)} \, \pm 0.03 \, \mbox{(syst.)} \, \pm 0.06 \, \mbox{(corr. syst.)}$, exhibiting a deviation of $2.6 \sigma$ from unity. The results are well-described by calculations that include nuclear modifications of the parton distribution functions, while the predictions using vacuum PDFs deviate from data by $2.3\sigma$ in the 0-90% centrality class and by $3\sigma$ in the 0-20% central collisions.

The second-order ($v_2$) and third-order ($v_3$) Fourier coefficients describing the azimuthal anisotropy of prompt and nonprompt (from b-hadron decays) J/$\psi$, as well as prompt $\psi$(2S) mesons are measured in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The analysis uses a data set corresponding to an integrated luminosity of 1.61 nb$^{-1}$ recorded with the CMS detector. The J/$\psi$ and $\psi$(2S) mesons are reconstructed using their dimuon decay channel. The $v_2$ and $v_3$ coefficients are extracted using the scalar product method and studied as functions of meson transverse momentum and collision centrality. The measured $v_2$ values for prompt J/$\psi$ mesons are found to be larger than those for nonprompt J/$\psi$ mesons. The prompt J/$\psi$$v_2$ values at high $p_\mathrm{T}$ are found to be underpredicted by a model incorporating only parton energy loss effects in a quark-gluon plasma medium. Prompt and nonprompt J/$\psi$ meson $v_3$ and prompt $\psi$(2S) $v_2$ and $v_3$ values are also reported for the first time, providing new information about heavy quark interactions in the hot and dense medium created in heavy ion collisions.

The differential cross sections of $\Lambda_\mathrm{c}^+$ baryon production are measured via the exclusive decay channel $\Lambda_\mathrm{c}^+ \to $pK$^-\pi^+$, as a function of transverse momentum ($p_\mathrm{T}$) in proton-proton (pp) and lead-lead (PbPb) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV with the CMS detector at the LHC. The measurement is performed within the $\Lambda_\mathrm{c}^+$ rapidity interval $|y|<$1.0 in the $p_\mathrm{T}$ range of 5-20 GeV/$c$ in pp and 10-20 GeV/$c$ in PbPb collisions. The observed yields of $\Lambda_\mathrm{c}^+$ for $p_\mathrm{T}$ of 10-20 GeV/$c$ suggest a possible suppression in central PbPb collisions compared to pp collisions. The $\Lambda_\mathrm{c}^+/$D$^0$ production ratio in pp collisions is compared to theoretical models. In PbPb collisions, this ratio is consistent with the result from pp collisions in their common $p_\mathrm{T}$ range.

The first evidence for X(3872) production in relativistic heavy ion collisions is reported. The X(3872) production is studied in lead-lead (PbPb) collisions at a center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV per nucleon pair, using the decay chain X(3872) $\to$ J$/\psi\, \pi^+\pi^- \to$ $\mu^+\mu^-\pi^+\pi^-$. The data were recorded with the CMS detector in 2018 and correspond to an integrated luminosity of 1.7 nb$^{-1}$. The measurement is performed in the rapidity and transverse momentum ranges $|y|$ $\lt$ 1.6 and 15 $\lt$ $p_\mathrm{T}$ $\lt$ 50 GeV$/c$. The significance of the inclusive X(3872) signal is 4.2 standard deviations. The prompt X(3872) to $\psi$(2S) yield ratio is found to be $\rho^\mathrm{PbPb} = $ 1.08 $\pm$ 0.49 (stat) $\pm$ 0.52 (syst), to be compared with typical values of 0.1 for pp collisions. This result provides a unique experimental input to theoretical models of the X(3872) production mechanism, and of the nature of this exotic state.

The first measurement of the dependence of $\gamma\gamma$$\to$$\mu^{+}\mu^{-}$ production on the multiplicity of neutrons emitted very close to the beam direction in ultraperipheral heavy ion collisions is reported. Data for lead-lead interactions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV, with an integrated luminosity of approximately 1.5 nb$^{-1}$, were collected using the CMS detector at the LHC. The azimuthal correlations between the two muons in the invariant mass region 8 $\lt$$m_{\mu\mu}$$\lt$ 60 GeV are extracted for events including 0, 1, or at least 2 neutrons detected in the forward pseudorapidity range $|\eta|$$\gt$ 8.3. The back-to-back correlation structure from leading-order photon-photon scattering is found to be significantly broader for events with a larger number of emitted neutrons from each nucleus, corresponding to interactions with a smaller impact parameter. This observation provides a data-driven demonstration that the average transverse momentum of photons emitted from relativistic heavy ions has an impact parameter dependence. These results provide new constraints on models of photon-induced interactions in ultraperipheral collisions. They also provide a baseline to search for possible final-state effects on lepton pairs caused by traversing a quark-gluon plasma produced in hadronic heavy ion collisions.

The $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

Two-particle Bose-Einstein momentum correlation functions are studied for charged-hadron pairs in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data sample, containing 4.27 $\times$$10^{9}$ minimum bias events corresponding to an integrated luminosity of 0.607 nb$^{-1}$, was collected by the CMS experiment in 2018. The experimental results are discussed in terms of a Lévy-type source distribution. The parameters of this distribution are extracted as functions of particle pair average transverse mass and collision centrality. These parameters include the Lévy index or shape parameter ($\alpha$), the Lévy scale parameter ($R$), and the correlation strength parameter ($\lambda$). The source shape, characterized by $\alpha$, is found to be neither Cauchy nor Gaussian, implying the need for a full Lévy analysis. Similarly to what was previously found for systems characterized by Gaussian source radii, a hydrodynamical scaling is observed for the Lévy $R$ parameter. The $\lambda$ parameter is studied in terms of the core-halo model.

Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark-gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we extracted the speed of sound in this medium created using lead-lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 $\pm$ 0.002 (stat) $\pm$ 0.016 (syst) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 $\pm$ 8 (syst) MeV. The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

The momentum-weighted sum of the electric charges of particles inside a jet, known as jet charge, is sensitive to the electric charge of the particle initiating the parton shower. This paper presents jet charge distributions in $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV lead-lead (PbPb) and proton-proton (pp) collisions recorded with the CMS detector at the LHC. These data correspond to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$ for PbPb and pp collisions, respectively. Leveraging the sensitivity of the jet charge to fundamental differences in the electric charges of quarks and gluons, the jet charge distributions from simulated events are used as templates to extract the quark- and gluon-like jet fractions from data. The modification of these jet fractions is examined by comparing pp and PbPb data as a function of the overlap of the colliding Pb nuclei (centrality). This measurement tests the color charge dependence of jet energy loss due to interactions with the quark-gluon plasma. No significant modification between different centrality classes and with respect to pp results is observed in the extracted fractions of quark- and gluon-like jet fractions.

Central exclusive and semiexclusive production of $\pi^+\pi^-$ pairs is measured with the CMS detector in proton-proton collisions at the LHC at center-of-mass energies of 5.02 and 13 TeV. The theoretical description of these nonperturbative processes, which have not yet been measured in detail at the LHC, poses a significant challenge to models. The two pions are measured and identified in the CMS silicon tracker based on specific energy loss, whereas the absence of other particles is ensured by calorimeter information. The total and differential cross sections of exclusive and semiexclusive central $\pi^+\pi^-$ production are measured as functions of invariant mass, transverse momentum, and rapidity of the $\pi^+\pi^-$ system in the fiducial region defined as transverse momentum $p_\mathrm{T}(\pi)$ $>$ 0.2 GeV and pseudorapidity $|\eta(\pi)|$ $<$ 2.4. The production cross sections for the four resonant channels f$_0(500)$, $\rho^0(770)$, f$_0(980)$, and f$_2(1270)$ are extracted using a simple model. These results represent the first measurement of this process at the LHC collision energies of 5.02 and 13 TeV.

The elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, ${\rm p+\overline{p}}$, ${\Lambda+\overline{\Lambda}}$, ${\rm K}^{\rm 0}_{\rm S}$, and the $\phi$-meson are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. Results obtained with the scalar product method are reported for the rapidity range $\vert y \vert <$ 0.5 as a function of transverse momentum, $p_\text{T}$, at different collision centrality intervals between 0-70%, including ultra-central (0-1%) collisions for $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$. For $p_\text{T} < 3$ GeV$\kern-0.05em/\kern-0.02em c$, the flow coefficients exhibit a particle mass dependence. At intermediate transverse momenta ($3<p_\text{T}<$~8-10 GeV$\kern-0.05em/\kern-0.02em c$), particles show an approximate grouping according to their type (i.e., mesons and baryons). The $\phi$-meson $v_2$, which tests both particle mass dependence and type scaling, follows ${\rm p+\overline{p}}$ $v_2$ at low $p_\text{T}$ and $\pi^{\pm}$ $v_2$ at intermediate $p_\text{T}$. The evolution of the shape of $v_{\rm n}(p_{\mathrm{T}})$ as a function of centrality and harmonic number $n$ is studied for the various particle species. Flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$ for $p_\text{T}<3$ GeV$\kern-0.05em/\kern-0.02em c$ are compared to iEBE-VISHNU and MUSIC hydrodynamical calculations coupled to a hadronic cascade model (UrQMD). The iEBE-VISHNU calculations describe the results fairly well for $p_\text{T} < 2.5$ GeV$\kern-0.05em/\kern-0.02em c$, while MUSIC calculations reproduce the measurements for $p_\text{T} < 1$ GeV$\kern-0.05em/\kern-0.02em c$. A comparison to $v_{\rm n}$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is also provided.