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$.

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.

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}|$.

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.

Correlations of two flow harmonics $v_n$ and $v_m$ via three- and four-particle cumulants are measured in 13 TeV $pp$, 5.02 TeV $p$+Pb, and 2.76 TeV peripheral Pb+Pb collisions with the ATLAS detector at the LHC. The goal is to understand the multi-particle nature of the long-range collective phenomenon in these collision systems. The large non-flow background from dijet production present in the standard cumulant method is suppressed using a method of subevent cumulants involving two, three and four subevents separated in pseudorapidity. The results show a negative correlation between $v_2$ and $v_3$ and a positive correlation between $v_2$ and $v_4$ for all collision systems and over the full multiplicity range. However, the magnitudes of the correlations are found to depend strongly on the event multiplicity, the choice of transverse momentum range and collision system. The relative correlation strength, obtained by normalisation of the cumulants with the $\langle v_n^2\rangle$ from a two-particle correlation analysis, is similar in the three collision systems and depends weakly on the event multiplicity and transverse momentum. These results based on the subevent methods provide strong evidence of a similar long-range multi-particle collectivity in $pp$, $p$+Pb and peripheral Pb+Pb collisions.

This paper presents studies of Bose-Einstein correlations (BEC) in proton-proton collisions at a centre-of-mass energy of 13 TeV, using data from the ATLAS detector at the CERN Large Hadron Collider. Data were collected in a special low-luminosity configuration with a minimum-bias trigger and a high-multiplicity track trigger, accumulating integrated luminosities of 151 $\mu$b$^{-1}$ and 8.4 nb$^{-1}$ respectively. The BEC are measured for pairs of like-sign charged particles, each with $|\eta|$ < 2.5, for two kinematic ranges: the first with particle $p_T$ > 100 MeV and the second with particle $p_T$ > 500 MeV. The BEC parameters, characterizing the source radius and particle correlation strength, are investigated as functions of charged-particle multiplicity (up to 300) and average transverse momentum of the pair (up to 1.5 GeV). The double-differential dependence on charged-particle multiplicity and average transverse momentum of the pair is also studied. The BEC radius is found to be independent of the charged-particle multiplicity for high charged-particle multiplicity (above 100), confirming a previous observation at lower energy. This saturation occurs independent of the transverse momentum of the pair.

A search is performed for the electroweak pair production of charginos and associated production of a chargino and neutralino, each of which decays through an $R$-parity-violating coupling into a lepton and a $W$, $Z$, or Higgs boson. The trilepton invariant-mass spectrum is constructed from events with three or more leptons, targeting chargino decays that include an electron or muon and a leptonically decaying $Z$ boson. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collision data produced by the Large Hadron Collider at a center-of-mass energy of $\sqrt{s}$ = 13 TeV and collected by the ATLAS experiment between 2015 and 2018. The data are found to be consistent with predictions from the Standard Model. The results are interpreted as limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model. Limits are also set on the production of charginos and neutralinos for a Minimal Supersymmetric Standard Model with an approximate $B$-$L$ symmetry. Charginos and neutralinos with masses between 100 GeV and 1100 GeV are excluded depending on the assumed decay branching fractions into a lepton (electron, muon, or $\tau$-lepton) plus a boson ($W$, $Z$, or Higgs).

A measurement of prompt photon-pair production in proton-proton collisions at $\sqrt{s}=13$ TeV is presented. The data were recorded by the ATLAS detector at the LHC with an integrated luminosity of 139 fb$^{-1}$. Events with two photons in the well-instrumented region of the detector are selected. The photons are required to be isolated and have a transverse momentum of $p_\mathrm{T,\gamma_{1(2)}} > 40(30)$ GeV for the leading (sub-leading) photon. The differential cross sections as functions of several observables for the diphoton system are measured and compared with theoretical predictions from state-of-the-art Monte Carlo and fixed-order calculations. The QCD predictions from next-to-next-to-leading-order calculations and multi-leg merged calculations are able to describe the measured integrated and differential cross sections within uncertainties, whereas lower-order calculations show significant deviations, demonstrating that higher-order perturbative QCD corrections are crucial for this process. The resummed predictions with parton showers additionally provide an excellent description of the low transverse-momentum regime of the diphoton system.

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.