A search for the production of three massive vector bosons in proton--proton collisions is performed using data at $\sqrt{s}=13\,TeV$ recorded with the ATLAS detector at the Large Hadron Collider in the years 2015--2017, corresponding to an integrated luminosity of $79.8\,\text{fb}^{-1}$. Events with two same-sign leptons $\ell$ (electrons or muons) and at least two reconstructed jets are selected to search for $WWW\to\ell\nu\ell\nu qq$. Events with three leptons without any same-flavour opposite-sign lepton pairs are used to search for $WWW\to\ell\nu\ell\nu\ell\nu$, while events with three leptons and at least one same-flavour opposite-sign lepton pair and one or more reconstructed jets are used to search for $WWZ\to\ell\nu qq \ell\ell$. Finally, events with four leptons are analysed to search for $WWZ\to\ell\nu\ell\nu\ell\ell$ and $WZZ\to qq \ell\ell\ell\ell$. Evidence for the joint production of three massive vector bosons is observed with a significance of 4.0 standard deviations, where the expectation is 3.1 standard deviations.

This letter describes the observation of the light-by-light scattering process, $\gamma\gamma\rightarrow\gamma\gamma$, in Pb+Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The analysis is conducted using a data sample corresponding to an integrated luminosity of 1.73 nb$^{-1}$, collected in November 2018 by the ATLAS experiment at the LHC. Light-by-light scattering candidates are selected in events with two photons produced exclusively, each with transverse energy $E_{\textrm{T}}^{\gamma} > 3$ GeV and pseudorapidity $|\eta_{\gamma}| < 2.37$, diphoton invariant mass above 6 GeV, and small diphoton transverse momentum and acoplanarity. After applying all selection criteria, 59 candidate events are observed for a background expectation of 12 $\pm$ 3 events. The observed excess of events over the expected background has a significance of 8.2 standard deviations. The measured fiducial cross section is 78 $\pm$ 13 (stat.) $\pm$ 7 (syst.) $\pm$ 3 (lumi.) nb.

An inclusive search for a new-physics signature of lepton-jet resonances has been performed by the ATLAS experiment. Scalar leptoquarks, pair-produced in $pp$ collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider, have been considered. An integrated luminosity of 3.2 fb$^{-1}$, corresponding to the full 2015 dataset was used. First (second) generation leptoquarks were sought in events with two electrons (muons) and two or more jets. The observed event yield in each channel is consistent with Standard Model background expectations. The observed (expected) lower limits on the leptoquark mass at 95% confidence level are 1100 GeV and 1050 GeV (1160 GeV and 1040 GeV) for first and second generation leptoquarks, respectively, assuming a branching ratio into a charged lepton and a quark of 100%. Upper limits on the aforementioned branching ratio are also given as a function of leptoquark mass. Compared with the results of earlier ATLAS searches, the sensitivity is increased for leptoquark masses above 860 GeV, and the observed exclusion limits confirm and extend the published results.

A search is presented for photonic signatures, motivated by generalized models of gauge-mediated supersymmetry breaking. This search makes use of proton-proton collision data at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ recorded by the ATLAS detector at the LHC, and it explores models dominated by both strong and electroweak production of supersymmetric partner states. Experimental signatures incorporating an isolated photon and significant missing transverse momentum are explored. These signatures include events with an additional photon or additional jet activity not associated with any specific underlying quark flavor. No significant excess of events is observed above the Standard Model prediction, and 95% confidence-level upper limits of between 0.083 fb and 0.32 fb are set on the visible cross section of contributions from physics beyond the Standard Model. These results are interpreted in terms of lower limits on the masses of gluinos, squarks, and gauginos in the context of generalized models of gauge-mediated supersymmetry, which reach as high as 2.3 TeV for strongly produced and 1.3 TeV for weakly produced supersymmetric partner pairs.

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.

The inclusive production rates of isolated, prompt photons in $p$+Pb collisions at $\sqrt{s_\mathrm{NN}} = 8.16$ TeV are studied with the ATLAS detector at the Large Hadron Collider using a dataset with an integrated luminosity of 165 nb$^{-1}$ recorded in 2016. The cross-section and nuclear modification factor $R_{p\mathrm{Pb}}$ are measured as a function of photon transverse energy from 20 GeV to 550 GeV and in three nucleon-nucleon centre-of-mass pseudorapidity regions, (-2.83,-2.02), (-1.84,0.91), and (1.09,1.90). The cross-section and $R_{p\mathrm{Pb}}$ values are compared with the results of a next-to-leading-order perturbative QCD calculation, with and without nuclear parton distribution function modifications, and with expectations based on a model of the energy loss of partons prior to the hard scattering. The data disfavour a large amount of energy loss and provide new constraints on the parton densities in nuclei.

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.

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.

A measurement of the four-lepton invariant mass spectrum is made with the ATLAS detector, using an integrated luminosity of 36.1 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}$ = 13 TeV delivered by the Large Hadron Collider. The differential cross-section is measured for events containing two same-flavour opposite-sign lepton pairs. It exhibits a rich structure, with different mass regions dominated in the Standard Model by single $Z$ boson production, Higgs boson production, and $Z$ boson pair production, and non-negligible interference effects at high invariant masses. The measurement is compared with state-of-the-art Standard Model calculations, which are found to be consistent with the data. These calculations are used to interpret the data in terms of $gg\rightarrow ZZ \rightarrow 4\ell$ and $Z \rightarrow 4\ell$ subprocesses, and to place constraints on a possible contribution from physics beyond the Standard Model.

The result of a search for pair production of the supersymmetric partner of the Standard Model bottom quark ($\tilde{b}_1$) is reported. The search uses 3.2 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}=$13 TeV collected by the ATLAS experiment at the Large Hadron Collider in 2015. Bottom squarks are searched for in events containing large missing transverse momentum and exactly two jets identified as originating from $b$-quarks. No excess above the expected Standard Model background yield is observed. Exclusion limits at 95% confidence level on the mass of the bottom squark are derived in phenomenological supersymmetric $R$-parity-conserving models in which the $\tilde{b}_1$ is the lightest squark and is assumed to decay exclusively via $\tilde{b}_1 \rightarrow b \tilde{\chi}_1^0$, where $\tilde{\chi}_1^0$ is the lightest neutralino. The limits significantly extend previous results; bottom squark masses up to 800 (840) GeV are excluded for the $\tilde{\chi}_1^0$ mass below 360 (100) GeV whilst differences in mass above 100 GeV between the $\tilde{b}_1$ and the $\tilde{\chi}_1^0$ are excluded up to a $\tilde{b}_1$ mass of 500 GeV.

The number of charged particles inside jets is a widely used discriminant for identifying the quark or gluon nature of the initiating parton and is sensitive to both the perturbative and non-perturbative components of fragmentation. This paper presents a measurement of the average number of charged particles with $p_\text{T}>500$ MeV inside high-momentum jets in dijet events using 20.3 fb$^{-1}$ of data recorded with the ATLAS detector in $pp$ collisions at $\sqrt{s}=8$ TeV collisions at the LHC. The jets considered have transverse momenta from 50 GeV up to and beyond 1.5 TeV. The reconstructed charged-particle track multiplicity distribution is unfolded to remove distortions from detector effects and the resulting charged-particle multiplicity is compared to several models. Furthermore, quark and gluon jet fractions are used to extract the average charged-particle multiplicity for quark and gluon jets separately.

The $ZZ$ production cross section in proton-proton collisions at 13 TeV center-of-mass energy is measured using 3.2 fb$^{-1}$ of data recorded with the ATLAS detector at the Large Hadron Collider. The considered $Z$ boson candidates decay to an electron or muon pair of mass 66-116 GeV. The cross section is measured in a fiducial phase space reflecting the detector acceptance. It is also extrapolated to a total phase space for $Z$ bosons in the same mass range and of all decay modes, giving $16.7^{+2.2}_{-2.0}$(stat.)$^{+0.9}_{-0.7}$(syst.)$^{+1.0}_{-0.7}$(lumi.) pb. The results agree with standard model predictions.

Results from a search for supersymmetry in events with four or more charged leptons (electrons, muons and taus) are presented. The analysis uses a data sample corresponding to 36.1 fb$^{-1}$ of proton-proton collisions delivered by the Large Hadron Collider at $\sqrt{s}=13$ TeV and recorded by the ATLAS detector. Four-lepton signal regions with up to two hadronically decaying taus are designed to target a range of supersymmetric scenarios that can be either enriched in or depleted of events involving the production and decay of a $Z$ boson. Data yields are consistent with Standard Model expectations and results are used to set upper limits on the event yields from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of General Gauge Mediated supersymmetry, where higgsino masses are excluded up to 295 GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of 1.46 TeV, 1.06 TeV, and 2.25 TeV are placed on wino, slepton and gluino masses, respectively.

Two-particle pseudorapidity correlations are measured in $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV Pb+Pb, $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV $p$+Pb, and $\sqrt{s}$ = 13 TeV $pp$ collisions at the LHC, with total integrated luminosities of approximately 7 $\mu\mathrm{b}^{-1}$, 28 $\mathrm{nb}^{-1}$, and 65 $\mathrm{nb}^{-1}$, respectively. The correlation function $C_{\rm N}(\eta_1,\eta_2)$ is measured as a function of event multiplicity using charged particles in the pseudorapidity range $|\eta|<2.4$. The correlation function contains a significant short-range component, which is estimated and subtracted. After removal of the short-range component, the shape of the correlation function is described approximately by $1+\langle{a_1^2}\rangle \eta_1\eta_2$ in all collision systems over the full multiplicity range. The values of $\sqrt{\langle{a_1^2}\rangle}$ are consistent between the opposite-charge pairs and same-charge pairs, and for the three collision systems at similar multiplicity. The values of $\sqrt{\langle{a_1^2}\rangle}$ and the magnitude of the short-range component both follow a power-law dependence on the event multiplicity. The $\eta$ distribution of the short-range component, after symmetrizing the proton and lead directions in $p$+Pb collisions, is found to be smaller than that in $pp$ collisions with comparable multiplicity.

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.

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.

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.

A search for direct pair production of top squarks in final states with two tau leptons, $b$-jets, and missing transverse momentum is presented. The analysis is based on proton-proton collision data at $\sqrt{s} = 13$ TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ recorded with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. Two exclusive channels with either two hadronically decaying tau leptons or one hadronically and one leptonically decaying tau lepton are considered. No significant deviation from the Standard Model predictions is observed in the data. The analysis results are interpreted in terms of model-independent limits and used to derive exclusion limits on the masses of the top squark $\tilde t_1$ and the tau slepton $\tilde \tau_1$ in a simplified model of supersymmetry with a nearly massless gravitino. In this model, masses up to $m(\tilde t_1) = 1.16$ TeV and $m(\tilde \tau_1) = 1.00$ TeV are excluded at 95% confidence level.

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.

Results of a search for the pair production of photon-jets$-$collimated groupings of photons$-$in the ATLAS detector at the Large Hadron Collider are reported. Highly collimated photon-jets can arise from the decay of new, highly boosted particles that can decay to multiple photons collimated enough to be identified in the electromagnetic calorimeter as a single, photonlike energy cluster. Data from proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 36.7 fb$^{-1}$, were collected in 2015 and 2016. Candidate photon-jet pair production events are selected from those containing two reconstructed photons using a set of identification criteria much less stringent than that typically used for the selection of photons, with additional criteria applied to provide improved sensitivity to photon-jets. Narrow excesses in the reconstructed diphoton mass spectra are searched for. The observed mass spectra are consistent with the Standard Model background expectation. The results are interpreted in the context of a model containing a new, high-mass scalar particle with narrow width, $X$, that decays into pairs of photon-jets via new, light particles, $a$. Upper limits are placed on the cross section times the product of branching ratios $\sigma \times \mathcal{B}(X \rightarrow aa) \times \mathcal {B}(a \rightarrow \gamma \gamma)^{2}$ for 200 GeV $< m_{X} <$ 2 TeV and for ranges of $ m_a $ from a lower mass of 100 MeV up to between 2 and 10 GeV, depending upon $ m_X $. Upper limits are also placed on $\sigma \times \mathcal{B}(X \rightarrow aa) \times \mathcal {B}(a \rightarrow 3\pi^{0})^{2}$ for the same range of $ m_X $ and for ranges of $ m_a $ from a lower mass of 500 MeV up to between 2 and 10 GeV.

A search for magnetic monopoles and high-electric-charge objects is presented using 34.4 fb$^{-1}$ of 13 TeV $pp$ collision data collected by the ATLAS detector at the LHC during 2015 and 2016. The considered signature is based upon high ionization in the transition radiation tracker of the inner detector associated with a pencil-shape energy deposit in the electromagnetic calorimeter. The data were collected by a dedicated trigger based on the tracker high-threshold hit capability. The results are interpreted in models of Drell-Yan pair production of stable particles with two spin hypotheses (0 and 1/2) and masses ranging from 200 GeV to 4000 GeV. The search improves by approximately a factor of five the constraints on the direct production of magnetic monopoles carrying one or two Dirac magnetic charges and stable objects with electric charge in the range $20\le|z|\le60$ and extends the charge range to $60<|z|\le100$.

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.

Multi-particle azimuthal cumulants are measured as a function of centrality and transverse momentum using 470 $\mu$b$^{-1}$ of Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ATLAS detector at the LHC. These cumulants provide information on the event-by-event fluctuations of harmonic flow coefficients $v_n$ and correlated fluctuations between two harmonics $v_n$ and $v_m$. For the first time, a non-zero four-particle cumulant is observed for dipolar flow, $v_1$. The four-particle cumulants for elliptic flow, $v_2$, and triangular flow, $v_3$, exhibit a strong centrality dependence and change sign in ultra-central collisions. This sign change is consistent with significant non-Gaussian fluctuations in $v_2$ and $v_3$. The four-particle cumulant for quadrangular flow, $v_4$, is found to change sign in mid-central collisions. Correlations between two harmonics are studied with three- and four-particle mixed-harmonic cumulants, which indicate an anti-correlation between $v_2$ and $v_3$, and a positive correlation between $v_2$ and $v_4$. These correlations decrease in strength towards central collisions and either approach zero or change sign in ultra-central collisions. To investigate the possible flow fluctuations arising from intrinsic centrality or volume fluctuations, the results are compared between two different event classes used for centrality definitions. In peripheral and mid-central collisions where the cumulant signals are large, only small differences are observed. In ultra-central collisions, the differences are much larger and transverse momentum dependent. These results provide new information to disentangle flow fluctuations from the initial and final states, as well as new insights on the influence of centrality fluctuations.