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.

This paper presents a measurement of the triple-differential cross section for the Drell--Yan process $Z/\gamma^*\rightarrow \ell^+\ell^-$ where $\ell$ is an electron or a muon. The measurement is performed for invariant masses of the lepton pairs, $m_{\ell\ell}$, between $46$ and $200$ GeV using a sample of $20.2$ fb$^{-1}$ of $pp$ collisions data at a centre-of-mass energy of $\sqrt{s}=8$ TeV collected by the ATLAS detector at the LHC in 2012. The data are presented in bins of invariant mass, absolute dilepton rapidity, $|y_{\ell\ell}|$, and the angular variable $\cos\theta^{*}$ between the outgoing lepton and the incoming quark in the Collins--Soper frame. The measurements are performed in the range $|y_{\ell\ell}|<2.4$ in the muon channel, and extended to $|y_{\ell\ell}|<3.6$ in the electron channel. The cross sections are used to determine the $Z$ boson forward-backward asymmetry as a function of $|y_{\ell\ell}|$ and $m_{\ell\ell}$. The measurements achieve high-precision, below the percent level in the pole region, excluding the uncertainty in the integrated luminosity, and are in agreement with predictions. These precision data are sensitive to the parton distribution functions and the effective weak mixing angle.

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.

A measurement of $b$-hadron pair production is presented, based on a data set corresponding to an integrated luminosity of 11.4 fb$^{-1}$ of proton--proton collisions recorded at $\sqrt{s}=8$ TeV with the ATLAS detector at the LHC. Events are selected in which a $b$-hadron is reconstructed in a decay channel containing $J/\psi \rightarrow \mu\mu$, and a second $b$-hadron is reconstructed in a decay channel containing a muon. Results are presented in a fiducial volume defined by kinematic requirements on three muons based on those used in the analysis. The fiducial cross section is measured to be $17.7 \pm 0.1 ($stat.$) \pm 2.0 ($syst.$)$ nb. A number of normalised differential cross sections are also measured, and compared to predictions from the Pythia8, Herwig++, MadGraph5\_aMC@NLO+Pythia8 and Sherpa event generators, providing new constraints on heavy flavour production.

Inclusive jet production cross-sections are measured in proton--proton collisions at a centre-of-mass energy of $\sqrt{s}=$8 TeV recorded by the ATLAS experiment at the Large Hadron Collider at CERN. The total integrated luminosity of the analysed data set amounts to $20.2$ fb$^{-1}$. Double-differential cross-sections are measured for jets defined by the anti-$k_{t}$ jet clustering algorithm with radius parameters of $R=0.4$ and $R=0.6$ and are presented as a function of the jet transverse momentum, in the range between 70 GeV and 2.5 TeV and in six bins of the absolute jet rapidity, between 0 and 3.0. The measured cross-sections are compared to predictions of quantum chromodynamics, calculated at next-to-leading order in perturbation theory, and corrected for non-perturbative and electroweak effects. The level of agreement with predictions, using a selection of different parton distribution functions for the proton, is quantified. Tensions between the data and the theory predictions are observed.

A search for the direct production of charginos and neutralinos in final states with at least two hadronically decaying tau leptons is presented. The analysis uses a dataset of $pp$ collisions corresponding to an integrated luminosity of 36.1 fb$^{-1}$, recorded with the ATLAS detector at the Large Hadron Collider at a centre-of-mass energy of 13 TeV. No significant deviation from the Standard Model background expectation is observed. Limits are derived in scenarios of $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ pair production and of $\tilde{\chi}_{1}^{\pm}\tilde{\chi}_{2}^{0}$ and $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ production in simplified models where the neutralinos and charginos decay solely via intermediate left-handed staus and tau sneutrinos, and the mass of the $\tilde{\tau}_{\mathrm L}$ state is set to be halfway between the masses of the $\tilde{\chi}_{1}^{\pm}$ and the $\tilde{\chi}_{1}^{0}$. Chargino masses up to 630 GeV are excluded at 95% confidence level in the scenario of direct production of $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ for a massless $\tilde{\chi}_{1}^{0}$. Common $\tilde{\chi}_{1}^{\pm}$, $\tilde{\chi}_{2}^{0}$ masses up to 760 GeV are excluded in the case of production of $\tilde{\chi}_{1}^{\pm}\tilde{\chi}_{2}^{0}$ and $\tilde{\chi}_{1}^{+}\tilde{\chi}_{1}^{-}$ assuming a massless $\tilde{\chi}_{1}^{0}$. Exclusion limits for additional benchmark scenarios with large and small mass-splitting between the $\tilde{\chi}_{1}^{\pm}$ and the $\tilde{\chi}_{1}^{0}$ are also studied by varying the $\tilde{\tau}_{\mathrm L}$ mass between the masses of the $\tilde{\chi}_{1}^{\pm}$ and the $\tilde{\chi}_{1}^{0}$.

Narrow resonances decaying into $WW$, $WZ$ or $ZZ$ boson pairs are searched for in 36.7 fb $^{-1}$ of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. The diboson system is reconstructed using pairs of large-radius jets with high transverse momentum and tagged as compatible with the hadronic decay of high-momentum $W$ or $Z$ bosons, using jet mass and substructure properties. The search is sensitive to diboson resonances with masses in the range 1.2-5.0 TeV. No significant excess is observed in any signal region. Exclusion limits are set at the 95% confidence level on the production cross section times branching ratio to dibosons for a range of theories beyond the Standard Model. Model-dependent lower limits on the mass of new gauge bosons are set, with the highest limit set at 3.5 TeV in the context of mass-degenerate resonances that couple predominantly to bosons.

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.

Jet substructure observables have significantly extended the search program for physics beyond the Standard Model at the Large Hadron Collider. The state-of-the-art tools have been motivated by theoretical calculations, but there has never been a direct comparison between data and calculations of jet substructure observables that are accurate beyond leading-logarithm approximation. Such observables are significant not only for probing the collinear regime of QCD that is largely unexplored at a hadron collider, but also for improving the understanding of jet substructure properties that are used in many studies at the Large Hadron Collider. This Letter documents a measurement of the first jet substructure quantity at a hadron collider to be calculated at next-to-next-to-leading-logarithm accuracy. The normalized, differential cross-section is measured as a function of log$_{10}\rho^2$, where $\rho$ is the ratio of the soft-drop mass to the ungroomed jet transverse momentum. This quantity is measured in dijet events from 32.9 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collisions recorded by the ATLAS detector. The data are unfolded to correct for detector effects and compared to precise QCD calculations and leading-logarithm particle-level Monte Carlo simulations.

A selection of searches by the ATLAS experiment at the LHC for the electroweak production of SUSY particles are used to study their impact on the constraints on dark matter candidates. The searches use $20\,{\rm fb}^{-1}$ of proton-proton collision data at $\sqrt{s}=8$ TeV. A likelihood-driven scan of a five-dimensional effective model focusing on the gaugino--higgsino and Higgs sector of the phenomenological minimal supersymmetric Standard Model is performed. This scan uses data from direct dark matter detection experiments, the relic dark matter density and precision flavour physics results. Further constraints from the ATLAS Higgs mass measurement and SUSY searches at LEP are also applied. A subset of models selected from this scan are used to assess the impact of the selected ATLAS searches in this five-dimensional parameter space. These ATLAS searches substantially impact those models for which the mass $m(\tilde{\chi}^0_1)$ of the lightest neutralino is less than 65 GeV, excluding 86% of such models. The searches have limited impact on models with larger $m(\tilde{\chi}^0_1)$ due to either heavy electroweakinos or compressed mass spectra where the mass splittings between the produced particles and the lightest supersymmetric particle is small.