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.

This Letter presents a search for magnetic monopoles with the ATLAS detector at the CERN Large Hadron Collider using an integrated luminosity of 2.0 fb^-1 of pp collisions recorded at a center-of-mass energy of sqrt(s)=7 TeV. No event is found in the signal region, leading to an upper limit on the production cross section at 95% confidence level of 1.6/epsilon fb for Dirac magnetic monopoles with the minimum unit magnetic charge and with mass between 200 GeV and 1500 GeV, where epsilon is the monopole reconstruction efficiency. The efficiency epsilon is high and uniform in the fiducial region given by pseudorapidity |eta|<1.37 and transverse kinetic energy 600-700<E^kin sin(theta)<1400 GeV. The minimum value of 700 GeV is for monopoles of mass 200 GeV, whereas the minimum value of 600 GeV is applicable for higher mass monopoles. Therefore, the upper limit on the production cross section at 95% confidence level is 2 fb in this fiducial region. Assuming the kinematic distributions from Drell-Yan pair production of spin-1/2 Dirac magnetic monopoles, the efficiency is in the range 1%-10%, leading to an upper limit on the cross section at 95% confidence level that varies from 145 fb to 16 fb for monopoles with mass between 200 GeV and 1200 GeV. This limit is weaker than the fiducial limit because most of these monopoles lie outside the fiducial region.

Exclusive dimuon production in ultraperipheral collisions (UPC), resulting from photon-photon interactions in the strong electromagnetic fields of colliding high-energy lead nuclei, $\mathrm{PbPb}(\gamma\gamma) \rightarrow \mu^+\mu^- (\mathrm{Pb}^{(\star)}\mathrm{Pb}^{(\star)} )$, is studied using $\mathcal{L}_{\mathrm{int}} = 0.48$ nb$^{-1}$ of $\sqrt{s_\mathrm{NN}}=5.02$ TeV lead-lead collision data at the LHC with the ATLAS detector. Dimuon pairs are measured in the fiducial region $p_{\mathrm{T}\mu} > 4$ GeV, $|\eta_{\mu}| < 2.4$, invariant mass $m_{\mu\mu} > 10$ GeV, and $p_{\mathrm{T,\mu\mu}} < 2$ GeV. The primary background from single-dissociative processes is extracted from the data using a template fitting technique. Differential cross sections are presented as a function of $m_{\mu\mu}$, absolute pair rapidity ($|y_{\mu\mu}|$), scattering angle in the dimuon rest frame ($|\cos \vartheta^{\star}_{\mu\mu}|$) and the colliding photon energies. The total cross section of the UPC $\gamma \gamma \rightarrow \mu^{+}\mu^{-}$ process in the fiducial volume is measured to be $\sigma_{\mathrm{fid}}^{\mu\mu} = 34.1 \! \pm \! 0.3 \mathrm{(stat.)} \! \pm \! 0.7 \mathrm{(syst.)}$ $\mu\mathrm{b}$. Generally good agreement is found with calculations from STARlight, which incorporate the leading-order Breit-Wheeler process with no final-state effects, albeit differences between the measurements and theoretical expectations are observed. In particular, the measured cross sections at larger $|y_{\mu\mu}|$ are found to be about 10-20% larger in data than in the calculations, suggesting the presence of larger fluxes of photons in the initial state. Modification of the dimuon cross sections in the presence of forward and/or backward neutron production is also studied and is found to be associated with a harder incoming photon spectrum, consistent with expectations.

Cross-sections for the production of a $Z$ boson in association with two photons are measured in proton$-$proton collisions at a centre-of-mass energy of 13 TeV. The data used correspond to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS experiment during Run 2 of the LHC. The measurements use the electron and muon decay channels of the $Z$ boson, and a fiducial phase-space region where the photons are not radiated from the leptons. The integrated $Z(\rightarrow\ell\ell)\gamma\gamma$ cross-section is measured with a precision of 12% and differential cross-sections are measured as a function of six kinematic variables of the $Z\gamma\gamma$ system. The data are compared with predictions from MC event generators which are accurate to up to next-to-leading order in QCD. The cross-section measurements are used to set limits on the coupling strengths of dimension-8 operators in the framework of an effective field theory.

A search for charged leptons with large impact parameters using 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV $pp$ collision data from the ATLAS detector at the LHC is presented, addressing a long-standing gap in coverage of possible new physics signatures. Results are consistent with the background prediction. This search provides unique sensitivity to long-lived scalar supersymmetric lepton-partners (sleptons). For lifetimes of 0.1 ns, selectron, smuon and stau masses up to 720 GeV, 680 GeV, and 340 GeV are respectively excluded at 95% confidence level, drastically improving on the previous best limits from LEP.

Constraints on the Higgs boson self-coupling are set by combining double-Higgs boson analyses in the $b\bar{b}b\bar{b}$, $b\bar{b}\tau^+\tau^-$ and $b\bar{b} \gamma \gamma$ decay channels with single-Higgs boson analyses targeting the $\gamma \gamma$, $ZZ^*$, $WW^*$, $\tau^+ \tau^-$ and $b\bar{b}$ decay channels. The data used in these analyses were recorded by the ATLAS detector at the LHC in proton$-$proton collisions at $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 126$-$139 fb$^{-1}$. The combination of the double-Higgs analyses sets an upper limit of $\mu_{HH} < 2.4$ at 95% confidence level on the double-Higgs production cross-section normalised to its Standard Model prediction. Combining the single-Higgs and double-Higgs analyses, with the assumption that new physics affects only the Higgs boson self-coupling ($\lambda_{HHH}$), values outside the interval $-0.4< \kappa_{\lambda}=(\lambda_{HHH}/\lambda_{HHH}^{\textrm{SM}})< 6.3$ are excluded at 95% confidence level. The combined single-Higgs and double-Higgs analyses provide results with fewer assumptions, by adding in the fit more coupling modifiers introduced to account for the Higgs boson interactions with the other Standard Model particles. In this relaxed scenario, the constraint becomes $-1.4 < \kappa_{\lambda} < 6.1$ at 95% CL.

A search for pair-produced vector-like quarks using events with exactly one lepton ($e$ or $\mu$), at least four jets including at least one $b$-tagged jet, and large missing transverse momentum is presented. Data from proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV, recorded by the ATLAS detector at the LHC from 2015 to 2018 and corresponding to an integrated luminosity of 139 fb$^{-1}$, are analysed. Vector-like partners $T$ and $B$ of the top and bottom quarks are considered, as is a vector-like $X$ with charge +5/3, assuming their decay into a $W$, $Z$, or Higgs boson and a third-generation quark. No significant deviations from the Standard Model expectation are observed. Upper limits on the production cross-section of $T$ and $B$ quark pairs as a function of their mass are derived for various decay branching ratio scenarios. The strongest lower limits on the masses are 1.59 TeV assuming mass-degenerate VLQs and branching ratios corresponding to the weak-isospin doublet model, and 1.47 TeV (1.46 TeV) for exclusive $T \rightarrow Zt$ ($B/X \rightarrow Wt$) decays. In addition, lower limits on the $T$ and $B$ quark masses are derived for all possible branching ratios.

Measurements of joint-polarisation states of $W$ and $Z$ gauge bosons in $W^{\pm}Z$ production are presented. The data set used corresponds to an integrated luminosity of $139$ fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of $13$ TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. The $W^{\pm}Z$ candidate events are reconstructed using leptonic decay modes of the gauge bosons into electrons and muons. The simultaneous pair-production of longitudinally polarised vector bosons is measured for the first time with a significance of $7.1$ standard deviations. The measured joint helicity fractions integrated over the fiducial region are $f_{\mathrm{00}} = 0.067 \pm 0.010$, $f_{\mathrm{0T}} = 0.110 \pm 0.029$, $f_{\mathrm{T0}} = 0.179 \pm 0.023$ and $f_{\mathrm{TT}} = 0.644 \pm 0.032$, in agreement with the next-to-leading-order Standard Model predictions. Individual helicity fractions of the $W$ and $Z$ bosons are also measured and found to be consistent with joint helicity fractions within the expected amount of correlations. Both the joint and individual helicity fractions are also measured separately in $W^+Z$ and $W^-Z$ events. Inclusive and differential cross sections for several kinematic observables sensitive to polarisation are presented.

This paper reports a search for Higgs boson pair ($hh$) production in association with a vector boson ($W$ or $Z$) using 139 $fb^{-1}$ of proton-proton collision data at $\sqrt{s}=$ 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The search is performed in final states in which the vector boson decays leptonically ($W\to\ell\nu, Z\to\ell\ell,\nu\nu$ with $\ell=e, \mu$) and the Higgs bosons each decay into a pair of $b$-quarks. It targets $Vhh$ signals from both non-resonant $hh$ production, present in the Standard Model (SM), and resonant $hh$ production, as predicted in some SM extensions. A 95% confidence-level upper limit of 183 (87) times the SM cross-section is observed (expected) for non-resonant $Vhh$ production when assuming the kinematics are as expected in the SM. Constraints are also placed on Higgs boson coupling modifiers. For the resonant search, upper limits on the production cross-sections are derived for two specific models: one is the production of a vector boson along with a neutral heavy scalar resonance $H$, in the mass range 260-1000 GeV, that decays into $hh$, and the other is the production of a heavier neutral pseudoscalar resonance $A$ that decays into a $Z$ boson and $H$ boson, where the $A$ boson mass is 360-800 GeV and the $H$ boson mass is 260-400 GeV. Constraints are also derived in the parameter space of two-Higgs-doublet models.

A combination of searches for resonant Higgs boson pair production is presented, using up to 139 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the LHC. The combination includes searches performed in three decay channels: $b\bar{b}b\bar{b}$, $bb\tau^+\tau^-$ and $bb\gamma\gamma$. No excess above the expected background is observed and upper limits are set at the 95% confidence level on the production cross section of Higgs boson pairs originating from the decay of a narrow scalar resonance with mass in the range 251 GeV-5 TeV. The observed (expected) limits are in the range 0.96-600 fb (1.2-390 fb). The limits are interpreted in the Type-I Two-Higgs-Doublet Model and the Minimimal Supersymmetric Standard Model, and constrain parameter space not previously excluded by other searches.