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.

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.

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.

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.

Ultrarelativistic heavy ion collisions recreate in the laboratory the thermodynamical conditions prevailing in the early universe up to 10$^{-6}$ seconds, thereby allowing the study of the quark-gluon plasma (QGP), a state of quantum chromodynamics (QCD) matter with deconfined partons. The top quark, the heaviest elementary particle known, is accessible in nucleus-nucleus collisions at the CERN LHC, and constitutes a novel probe of the QGP. Here, we report the first-ever evidence for the production of top quarks in nucleus-nucleus collisions, using lead-lead collision data at a nucleon-nucleon centre-of-mass energy of 5.02 TeV recorded by the CMS experiment. Two methods are used to measure the cross section for top quark pair production ($\sigma_\mathrm{t\bar{t}}$) via the decay into charged leptons (electrons or muons) and bottom quarks. One method relies on the leptonic information alone, and the second one exploits, in addition, the presence of bottom quarks. The measured cross sections, $\sigma_\mathrm{t\bar{t}} = $ 2.54 $^{+0.84}_{-0.74}$ and 2.03 $^{+0.71}_{-0.64}$ $\mu$b, respectively, are compatible with expectations from scaled proton-proton data and QCD predictions.

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 spectra measurements of charged hadrons produced in the shower of a parton originating in the same hard scattering with a leptonically decaying Z boson, are reported in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. Both PbPb and pp data sets are recorded by the CMS experiment at the LHC, and correspond to an integrated luminosity of 1.7 nb$^{-1}$ and 320 pb$^{-1}$, respectively. Hadronic collision data with one reconstructed Z boson candidate with the transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV/$c$ are analyzed. The Z boson constrains the initial energy and direction of the associated parton. In heavy ion events, azimuthal angular distributions of charged hadrons with respect to the direction of a Z boson are sensitive to modifications of the in-medium parton shower and medium response. Compared to reference data from pp interactions, the results for central PbPb collisions indicate a modification of the angular correlations. The measurements of the fragmentation functions and $p_\mathrm{T}$ spectra of charged particles in Z boson events, which are sensitive to medium modifications of the parton shower longitudinal structure, are also reported. Significant modifications in central PbPb events compared to pp reference data are also found for these observables.

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.

Jet production in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV is studied with the CMS detector at the LHC, using PbPb and pp data samples corresponding to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$, respectively. Jets with different areas are reconstructed using the anti-$k_\mathrm{T}$ algorithm by varying the distance parameter $R$. The measurements are performed using jets with transverse momenta ($p_\mathrm{T}$) greater than 200 GeV and in a pseudorapidity range of $|\eta|$$\lt$ 2. To reveal the medium modification of the jet spectra in PbPb collisions, the properly normalized ratio of spectra from PbPb and pp data is used to extract jet nuclear modification factors as functions of the PbPb collision centrality, $p_\mathrm{T}$ and, for the first time, as a function of $R$ up to 1.0. For the most central collisions, a strong suppression is observed for high-$p_\mathrm{T}$ jets reconstructed with all distance parameters, implying that a significant amount of jet energy is scattered to large angles. The dependence of jet suppression on $R$ is expected to be sensitive to both the jet energy loss mechanism and the medium response, and so the data are compared to several modern event generators and analytic calculations. The models considered do not fully reproduce the data.

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