This Letter reports the observation of $WZ\gamma$ production and a measurement of its cross-section using 140.1 $\pm$ 1.2 fb$^{-1}$ of proton-proton collision data recorded at a center-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The $WZ\gamma$ production cross-section, with both the $W$ and $Z$ bosons decaying leptonically, $pp \rightarrow WZ\gamma \rightarrow {\ell'}^{\pm}\nu\ell^{+}\ell^{-}\gamma$ ($\ell^{(')} = e, \mu$), is measured in a fiducial phase-space region defined such that the leptons and the photon have high transverse momentum and the photon is isolated. The cross-section is found to be 2.01 $\pm$ 0.30 (stat.) $\pm$ 0.16 (syst) fb. The corresponding Standard Model predicted cross-section calculated at next-to-leading order in perturbative quantum chromodynamics and at leading order in the electroweak coupling constant is 1.50 $\pm$ 0.06 fb. The observed significance of the $WZ\gamma$ signal is 6.3$\sigma$, compared with an expected significance of 5.0$\sigma$.

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 Majorana neutrinos in same-sign $WW$ scattering events is presented. The analysis uses $\sqrt{s}= 13$ TeV proton-proton collision data with an integrated luminosity of 140 fb$^{-1}$ recorded during 2015-2018 by the ATLAS detector at the Large Hadron Collider. The analysis targets final states including exactly two same-sign muons and at least two hadronic jets well separated in rapidity. The modelling of the main backgrounds, from Standard Model same-sign $WW$ scattering and $WZ$ production, is constrained with data in dedicated signal-depleted control regions. The distribution of the transverse momentum of the second-hardest muon is used to search for signals originating from a heavy Majorana neutrino with a mass between 50 GeV and 20 TeV. No significant excess is observed over the background expectation. The results are interpreted in a benchmark scenario of the Phenomenological Type-I Seesaw model. In addition, the sensitivity to the Weinberg operator is investigated. Upper limits at the 95% confidence level are placed on the squared muon-neutrino-heavy-neutrino mass-mixing matrix element $\vert V_{\mu N} \vert^{2}$ as a function of the heavy Majorana neutrino's mass $m_N$, and on the effective $\mu\mu$ Majorana neutrino mass $|m_{\mu\mu}|$.

Differential cross-section measurements of $Z\gamma$ production in association with hadronic jets are presented, using the full 139 fb$^{-1}$ dataset of $\sqrt{s}=13$ TeV proton-proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the $Z$ boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLO$_\text{PS}$, as well as next-to-leading-order predictions from MadGraph5_aMC@NLO and Sherpa.

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.

A search for leptoquarks decaying into the $b\tau$ final state is performed using Run 2 proton-proton collision data from the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb$^{-1}$ at $\sqrt{s} = 13$ TeV recorded by the ATLAS detector. The benchmark models considered in this search are vector leptoquarks with electric charge of 2/3e and scalar leptoquarks with an electric charge of 4/3e. No significant excess above the Standard Model prediction is observed, and 95% confidence level upper limits are set on the cross-section times branching fraction of leptoquarks decaying into $b\tau$. For the vector leptoquark production two models are considered: the Yang-Mills and Minimal coupling models. In the Yang-Mills (Minimal coupling) scenario, vector leptoquarks with a mass below 1.58 (1.35) TeV are excluded for a gauge coupling of 1.0 and below 2.05 (1.99) TeV for a gauge coupling of 2.5. In the case of scalar leptoquarks, masses below 1.28 TeV (1.53 TeV) are excluded for a Yukawa coupling of 1.0 (2.5). Finally, an interpretation of the results with minimal model dependence is performed for each of the signal region categories, and limits on the visible cross-section for beyond the Standard Model processes are provided.

We present a search for magnetic monopoles and high-electric-charge objects using LHC Run 2 $\sqrt{s} =$13 TeV proton$-$proton collisions recorded by the ATLAS detector. A total integrated luminosity of 138 fb$^{-1}$ was collected by a specialized trigger. No highly ionizing particle candidate was observed. Considering the Drell-Yan and photon-fusion pair production mechanisms as benchmark models, cross-section upper limits are presented for spin-0 and spin-$\frac{1}{2}$ magnetic monopoles of magnetic charge $1g_\textrm{D}$ and $2g_\textrm{D}$ and for high-electric-charge objects of electric charge $20 \leq |z| \leq 100$, for masses between 200 GeV and 4000 GeV. The search improves by approximately a factor of three the previous cross-section limits on the Drell-Yan production of magnetic monopoles and high-electric charge objects. Also, the first ATLAS limits on the photon-fusion pair production mechanism of magnetic monopoles and high-electric-charge objects have been obtained.

A search for the pair production of heavy leptons as predicted by the type-III seesaw mechanism is presented. The search uses proton-proton collision data at a centre-of-mass energy of 13 TeV, corresponding to 139 fb$^{-1}$ of integrated luminosity recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. The analysis focuses on final states with three or four electrons or muons from the possible decays of new heavy leptons via intermediate electroweak bosons. No significant deviations above the Standard Model expectation are observed; upper and lower limits on the heavy lepton production cross-section and masses are derived respectively. These results are then combined for the first time with the ones already published by ATLAS using the channel with two leptons in the final state. The observed lower limit on the mass of the type-III seesaw heavy leptons combining two, three and four lepton channels together is 910 GeV at the 95% confidence level.

A measurement of single top-quark production in the s-channel is performed in proton$-$proton collisions at a centre-of-mass energy of 13 TeV with the ATLAS detector at the CERN Large Hadron Collider. The dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The analysis is performed on events with an electron or muon, missing transverse momentum and exactly two $b$-tagged jets in the final state. A discriminant based on matrix element calculations is used to separate single-top-quark s-channel events from the main background contributions, which are top-quark pair production and $W$-boson production in association with jets. The observed (expected) signal significance over the background-only hypothesis is 3.3 (3.9) standard deviations, and the measured cross-section is $\sigma=8.2^{+3.5}_{-2.9}$ pb, consistent with the Standard Model prediction of $\sigma^{\mathrm{SM}}=10.32^{+0.40}_{-0.36}$ pb.

Measurements of differential cross sections are presented for inclusive isolated-photon production in $pp$ collisions at a centre-of-mass energy of 13 TeV provided by the LHC and using 139 fb$^{-1}$ of data recorded by the ATLAS experiment. The cross sections are measured as functions of the photon transverse energy in different regions of photon pseudorapidity. The photons are required to be isolated by means of a fixed-cone method with two different cone radii. The dependence of the inclusive-photon production on the photon isolation is investigated by measuring the fiducial cross sections as functions of the isolation-cone radius and the ratios of the differential cross sections with different radii in different regions of photon pseudorapidity. The results presented in this paper constitute an improvement with respect to those published by ATLAS earlier: the measurements are provided for different isolation radii and with a more granular segmentation in photon pseudorapidity that can be exploited in improving the determination of the proton parton distribution functions. These improvements provide a more in-depth test of the theoretical predictions. Next-to-leading-order QCD predictions from JETPHOX and SHERPA and next-to-next-to-leading-order QCD predictions from NNLOJET are compared to the measurements, using several parameterisations of the proton parton distribution functions. The measured cross sections are well described by the fixed-order QCD predictions within the experimental and theoretical uncertainties in most of the investigated phase-space region.