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.

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 with ${\rm K}^{0}_{\rm{S}}$, $\Lambda$/$\bar{\Lambda}$, and charged hadrons as trigger particles in the transverse momentum range $8 < p_\mathrm{T,trig}<16$ GeV/$c$, and associated charged particles within $1 < p_\mathrm{T,assoc}<8$ GeV/$c$, are studied at mid-rapidity in pp and central Pb-Pb collisions at a centre-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV with the ALICE detector at the LHC. After subtracting the contributions of the flow background, the per-trigger yields are extracted on both the near and away sides, and the ratio in Pb-Pb collisions with respect to pp collisions ($I_{\mathrm {AA}}$) is computed. The per-trigger yield in Pb-Pb collisions on the away side is strongly suppressed to the level of $I_{\mathrm {AA}} \approx 0.6$ for $p_\mathrm{T,assoc}>3$ GeV/$c$ as expected from strong in-medium energy loss, while an enhancement develops at low $p_\mathrm{T,assoc}$ on both the near and away sides, reaching $I_{\mathrm {AA}} \approx 1.8$ and $2.7$ respectively. These findings are in good agreement with previous ALICE measurements from two-particle correlations triggered by neutral pions ($\pi^{0}$-h) and charged hadrons (h-h) in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 2.76$ TeV. Moreover, the correlations with ${\rm K}^{0}_{\rm{S}}$ mesons and $\Lambda$/$\bar{\Lambda}$ baryons as trigger particles are compared to those of inclusive charged hadrons. The results are compared with the predictions of Monte Carlo models.

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 cross sections for $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) production in lead-lead (PbPb) and proton-proton (pp) collisions at $\sqrt{s_{_\mathrm{NN}}} =$ 5.02 TeV have been measured using the CMS detector at the LHC. The nuclear modification factors, R$_\mathrm{AA}$, derived from the PbPb-to-pp ratio of yields for each state, are studied as functions of meson rapidity and transverse momentum, as well as PbPb collision centrality. The yields of all three states are found to be significantly suppressed, and compatible with a sequential ordering of the suppression, R$_\mathrm{AA}$($\Upsilon$(1S)) $>$ R$_\mathrm{AA}$($\Upsilon$(2S)) $>$ R$_\mathrm{AA}$($\Upsilon$(3S)) . The suppression of $\Upsilon$(1S) is larger than that seen at $\sqrt{s_{_\mathrm{NN}}} =$ 2.76 TeV, although the two are compatible within uncertainties. The upper limit on the R$_\mathrm{AA}$ of $\Upsilon$(3S) integrated over $p_\mathrm{T}$ and rapidity is 0.094 at 95% confidence level, which is the strongest suppression observed for any hadron species in heavy ion collisions to date.

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.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

The azimuthal ($\Delta\varphi$) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p$-$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV. Results are reported for electrons with transverse momentum $4<p_{\rm T}<16$ GeV/$c$ and pseudorapidity $|\eta|<0.6$. The associated charged particles are selected with transverse momentum $1<p_{\rm T}<7$ GeV/$c$, and relative pseudorapidity separation with the leading electron $|\Delta\eta| < 1$. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p$-$Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The $\Delta\varphi$ distribution and the peak observables in pp and p$-$Pb collisions are compared with calculations from various Monte Carlo event generators.

The transverse-momentum ($p_{\rm T}$) spectra of ${\rm K}^*(892)^0$ and $\phi(1020)$ measured with the ALICE detector up to $p_{\rm T}$ = 16 GeV/$c$ in the rapidity range $-1.2 < y < 0.3$, in p-Pb collisions at the center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured $p_{\rm T}$ distributions show a dependence on both multiplicity and rapidity at low $p_{\rm T}$ whereas no significant dependence is observed at high $p_{\rm T}$. A rapidity dependence is observed in the $p_{\rm T}$-integrated yield (d$N$/d$y$), whereas the mean transverse momentum ($\left< p_{\rm T} \right>$) shows a flat behavior as a function of rapidity. The rapidity asymmetry ($Y_{\rm asym}$) at low $p_{\rm T}$ ( < 5 GeV/$c$) is more significant for higher multiplicity classes. At high $p_{\rm T}$, no significant rapidity asymmetry is observed in any of the multiplicity classes. Both ${\rm K}^*(892)^0$ and $\phi(1020)$ show similar $Y_{\rm asym}$. The nuclear modification factor ($Q_{\rm CP}$) as a function of $p_{\rm T}$ shows a Cronin-like enhancement at intermediate $p_{\rm T}$, which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high $p_{\rm T}$ (> 5 GeV/$c$), the $Q_{\rm CP}$ values are greater than unity and no significant rapidity dependence is observed.

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.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_\mathrm{NN}} =$ 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

The production of Z bosons is studied in the dimuon and dielectron decay channels in PbPb and pp collisions at sqrt(s[NN]) = 2.76 TeV, using data collected by the CMS experiment at the LHC. The PbPb data sample corresponds to an integrated luminosity of about 150 inverse microbarns, while the pp data sample collected in 2013 at the same nucleon-nucleon centre-of-mass energy has an integrated luminosity of 5.4 inverse picobarns. The Z boson yield is measured as a function of rapidity, transverse momentum, and collision centrality. The ratio of PbPb to pp yields, scaled by the number of inelastic nucleon-nucleon collisions, is found to be 1.06 +/- 0.05 (stat) +/- 0.08 (syst) in the dimuon channel and 1.02 +/- 0.08 (stat) +/- 0.15 (syst) in the dielectron channel, for centrality-integrated Z boson production. This binary collision scaling is seen to hold in the entire kinematic region studied, as expected for a colourless probe that is unaffected by the hot and dense QCD medium produced in heavy ion collisions.

A measurement is presented of the charged hadron multiplicity in hadronic PbPb collisions, as a function of pseudorapidity and centrality, at a collision energy of 2.76 TeV per nucleon pair. The data sample is collected using the CMS detector and a minimum-bias trigger, with the CMS solenoid off. The number of charged hadrons is measured both by counting the number of reconstructed particle hits and by forming hit doublets of pairs of layers in the pixel detector. The two methods give consistent results. The charged hadron multiplicity density dN(ch)/d eta, evaluated at eta=0 for head-on collisions, is found to be 1612 +/- 55, where the uncertainty is dominated by systematic effects. Comparisons of these results to previous measurements and to various models are also presented.

Evidence for the light-by-light scattering process, $\gamma\gamma$ $\to$ $\gamma\gamma$, in ultraperipheral PbPb collisions at a centre-of-mass energy per nucleon pair of 5.02 TeV is reported. The analysis is conducted using a data sample corresponding to an integrated luminosity of 390 $\mu$b$^{-1}$ recorded by the CMS experiment at the LHC. Light-by-light scattering processes are selected in events with two photons exclusively produced, each with transverse energy E$_\mathrm{T}^{\gamma}$ $>$ 2 GeV, pseudorapidity $|\eta^{\gamma}|$ $\lt$ 2.4, diphoton invariant mass $m^{\gamma\gamma}$ $\gt$ 5 GeV, diphoton transverse momentum $p_\mathrm{T}^{\gamma\gamma}$ $\lt$ 1 GeV, and diphoton acoplanarity below 0.01. After all selection criteria are applied, 14 events are observed, compared to expectations of 9.0 $\pm$ 0.9 (theo) events for the signal and 4.0 $\pm$ 1.2 (stat) for the background processes. The excess observed in data relative to the background-only expectation corresponds to a significance of 3.7 standard deviations, and has properties consistent with those expected for the light-by-light scattering signal. The measured fiducial light-by-light scattering cross section, $\sigma_\mathrm{fid} (\gamma\gamma$ $\to$ $\gamma\gamma) =$ 120 $\pm$ 46 (stat) $\pm$ 28 (syst) $\pm$ 12 (theo) nb, is consistent with the standard model prediction. The $m^{\gamma\gamma}$ distribution is used to set new exclusion limits on the production of pseudoscalar axion-like particles, via the $\gamma\gamma$ $\to$ a $\to$ $\gamma\gamma$ process, in the mass range $m_{\mathrm{a}} =$ 5-90 GeV.

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.

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 anisotropy of the azimuthal distributions of charged particles produced in PbPb collisions with a nucleon-nucleon center-of-mass energy of 2.76 TeV is studied with the CMS experiment at the LHC. The elliptic anisotropy parameter defined as the second coefficient in a Fourier expansion of the particle invariant yields, is extracted using the event-plane method, two- and four-particle cumulants, and Lee--Yang zeros. The anisotropy is presented as a function of transverse momentum (pt), pseudorapidity (eta) over a broad kinematic range: 0.3 < pt < 20 GeV, abs(eta) < 2.4, and in 12 classes of collision centrality from 0 to 80%. The results are compared to those obtained at lower center-of-mass energies, and various scaling behaviors are examined. When scaled by the geometric eccentricity of the collision zone, the elliptic anisotropy is found to obey a universal scaling with the transverse particle density for different collision systems and center-of-mass energies.

The relative modification of the prompt psi(2S) and J/psi yields from pp to PbPb collisions, at the center of mass energy of 5.02 TeV per nucleon pair, is presented. The analysis is based on pp and PbPb data samples collected by the CMS experiment at the LHC in 2015, corresponding to integrated luminosities of 28.0 inverse picobarns and 464 inverse microbarns, respectively. The double ratio of measured yields of prompt charmonia reconstructed through their decays into muon pairs, (N[psi(2S)]/N[J/psi])[PbPb] / (N[psi(2S)]/N[J/psi])[pp], is determined as a function of PbPb collision centrality and charmonium transverse momentum pt, in two kinematic intervals: abs(y) < 1.6 covering 6.5 < pt < 30 GeV/c and 1.6 < abs(y) < 2.4 covering 3 < pt < 30 GeV/c. The centrality-integrated double ratios are 0.36 +/- 0.08 (stat) +/-0.05 (syst) in the first interval and 0.24 +/- 0.22 (stat) +/- 0.09 (syst) in the second. The double ratio is lower than unity in all the measured bins, suggesting that the psi(2S) yield is more suppressed than the J/psi yield in the explored phase space.

Measurements of azimuthal angle and transverse momentum ($p_\mathrm{T}$) correlations of isolated photons and associated jets are reported for pp and PbPb collisions at $\sqrt{s_{_{\mathrm{NN}}}} =$ 5.02 TeV. The data were recorded with the CMS detector at the CERN LHC. For events containing a leading isolated photon with $p_\mathrm{T}^\gamma >$ 40 GeV$/c$ and an associated jet with $p_\mathrm{T}^\text{jet} >$ 30 GeV$/c$, the photon+jet azimuthal correlation and $p_\mathrm{T}$ imbalance in PbPb collisions are studied as functions of collision centrality and $p_\mathrm{T}^\gamma$. The results are compared to pp reference data collected at the same collision energy and to predictions from several theoretical models for parton energy loss. No evidence of broadening of the photon+jet azimuthal correlations is observed, while the ratio $p_\mathrm{T}^\text{jet}/p_\mathrm{T}^\gamma$ decreases significantly for PbPb data relative to the pp reference. All models considered agree within uncertainties with the data. The number of associated jets per photon with $p_\mathrm{T}^\gamma >$ 80 GeV$/c$ is observed to be shifted towards lower $p_\mathrm{T}^\text{jet}$ values in central PbPb collisions compared to pp collisions.

The ratio between the prompt psi(2S) and J/psi yields, reconstructed via their decays into muon pairs, is measured in PbPb and pp collisions at sqrt(s[NN]) = 2.76 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC, corresponding to integrated luminosities of 150 inverse microbarns and 5.4 inverse picobarns, respectively. The double ratio of measured yields, (N[psi(2S)]/N[J/psi])[PbPb] / (N[psi(2S)]/ N[J/psi])[pp], is computed in three PbPb collision centrality bins and two kinematic ranges: one at midrapidity, abs(y) < 1.6, covering the transverse momentum range 6.5 < pt < 30 GeV/c, and the other at forward rapidity, 1.6 < abs(y) < 2.4, extending to lower pt values, 3 < pt < 30 GeV/c. The centrality-integrated double ratio changes from 0.45 +/- 0.13 (stat) +/- 0.07 (syst) in the first range to 1.67 +/- 0.34 (stat) +/- 0.27 (syst) in the second. This difference is most pronounced in the most central collisions.

Measurements of prompt $\psi$(2S) meson production cross sections in proton-lead (pPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{_\mathrm{NN}}}=$ 5.02 TeV are reported. The results are based on pPb and pp data collected by the CMS experiment at the LHC, corresponding to integrated luminosities of 34.6 nb$^{-1}$ and 28.0 pb$^{-1}$, respectively. The nuclear modification factor $R_\mathrm{pPb}$ is measured for prompt $\psi$(2S) in the transverse momentum range 4 $<$ p$_\mathrm{T}$ $<$ 30 GeV$/c$ and the center-of-mass rapidity range $-$2.4 $< y_\mathrm{cm} <$ 1.93. The results on $\psi$(2S) $R_\mathrm{pPb}$ are compared to the corresponding modification factor for prompt J$/\psi$ mesons and are found to be more suppressed than the J$/\psi$ states over the entire kinematic range studied.

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.

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.

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