A search for high-mass charged and neutral bosons decaying to $W\gamma$ and $Z\gamma$ final states is presented in this paper. The analysis uses a data sample of $\sqrt{s} = 13$ TeV proton-proton collisions with an integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector during LHC Run 2 operation. The sensitivity of the search is determined using models of the production and decay of spin-1 charged bosons and spin-0/2 neutral bosons. The range of resonance masses explored extends from 1.0 TeV to 6.8 TeV. At these high resonance masses, it is beneficial to target the hadronic decays of the $W$ and $Z$ bosons because of their large branching fractions. The decay products of the high-momentum $W/Z$ bosons are strongly collimated and boosted-boson tagging techniques are employed to improve the sensitivity. No evidence of a signal above the Standard Model backgrounds is observed, and upper limits on the production cross-sections of these bosons times their branching fractions to $W\gamma$ and $Z\gamma$ are derived for various boson production models.

A search for forward proton scattering in association with light-by-light scattering mediated by an axion-like particle is presented, using the ATLAS Forward Proton spectrometer to detect scattered protons and the central ATLAS detector to detect pairs of outgoing photons. Proton-proton collision data recorded in 2017 at a centre-of-mass energy of $\sqrt{s} = 13$ TeV were analysed, corresponding to an integrated luminosity of 14.6 fb$^{-1}$. A total of 441 candidate signal events were selected. A search was made for a narrow resonance in the diphoton mass distribution, corresponding to an axion-like particle (ALP) with mass in the range 150-1600 GeV. No excess is observed above a smooth background. Upper limits on the production cross section of a narrow resonance are set as a function of the mass, and are interpreted as upper limits on the ALP production coupling constant, assuming 100% decay branching ratio into a photon pair. The inferred upper limit on the coupling constant is in the range 0.04-0.09 TeV$^{-1}$ at 95%confidence level.

A search for heavy right-handed Majorana or Dirac neutrinos $N_{\mathrm{R}}$ and heavy right-handed gauge bosons $W_{\mathrm{R}}$ is performed in events with energetic electrons or muons, with the same or opposite electric charge, and energetic jets. The search is carried out separately for topologies of clearly separated final-state products (``resolved'' channel) and topologies with boosted final states with hadronic and/or leptonic products partially overlapping and reconstructed as a large-radius jet (``boosted'' channel). The events are selected from $pp$ collision data at the LHC with an integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector at $\sqrt{s}$ = 13 TeV. No significant deviations from the Standard Model predictions are observed. The results are interpreted within the theoretical framework of a left-right symmetric model, and lower limits are set on masses in the heavy right-handed $W_{\mathrm{R}}$ boson and $N_{\mathrm{R}}$ plane. The excluded region extends to about $m(W_{\mathrm{R}}) = 6.4$ TeV for both Majorana and Dirac $N_{\mathrm{R}}$ neutrinos at $m(N_{\mathrm{R}})<1$ TeV. $N_{\mathrm{R}}$ with masses of less than 3.5 (3.6) TeV are excluded in the electron (muon) channel at $m(W_{\mathrm{R}})=4.8$ TeV for the Majorana neutrinos, and limits of $m(N_{\mathrm{R}})$ up to 3.6 TeV for $m(W_{\mathrm{R}}) = 5.2$ (5.0) TeV in the electron (muon) channel are set for the Dirac neutrinos. These constitute the most stringent exclusion limits to date for the model considered.

Quasireal photons exchanged in relativistic heavy ion interactions are powerful probes of the gluonic structure of nuclei. The coherent J/$\psi$ photoproduction cross section in ultraperipheral lead-lead collisions is measured as a function of photon-nucleus center-of-mass energies per nucleon (W$^\text{Pb}_{\gamma\text{N}}$), over a wide range of 40 $\lt$ W$^\text{Pb}_{\gamma\text{N}}$$\lt$ 400 GeV. Results are obtained using data at the nucleon-nucleon center-of-mass energy of 5.02 TeV collected by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 1.52 nb$^{-1}$. The cross section is observed to rise rapidly at low W$^\text{Pb}_{\gamma\text{N}}$, and plateau above W$^\text{Pb}_{\gamma\text{N}}$$\approx$ 40 GeV, up to 400 GeV, a new regime of small Bjorken-$x$ ($\approx$ 6 $\times$ 10$^{-5}$) gluons being probed in a heavy nucleus. The observed energy dependence is not predicted by current quantum chromodynamic models.

Differential and double-differential distributions of kinematic variables of leptons from decays of top-quark pairs ($t\bar{t}$) are measured using the full LHC Run 2 data sample collected with the ATLAS detector. The data were collected at a $pp$ collision energy of $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 140 fb$^{-1}$. The measurements use events containing an oppositely charged $e\mu$ pair and $b$-tagged jets. The results are compared with predictions from several Monte Carlo generators. While no prediction is found to be consistent with all distributions, a better agreement with measurements of the lepton $p_{\text{T}}$ distributions is obtained by reweighting the $t\bar{t}$ sample so as to reproduce the top-quark $p_{\text{T}}$ distribution from an NNLO calculation. The inclusive top-quark pair production cross-section is measured as well, both in a fiducial region and in the full phase-space. The total inclusive cross-section is found to be \[ \sigma_{t\bar{t}} = 829 \pm 1\;(\textrm{stat}) \pm 13\;(\textrm{syst}) \pm 8\;(\textrm{lumi}) \pm 2\; (\textrm{beam})\ \textrm{pb}, \] where the uncertainties are due to statistics, systematic effects, the integrated luminosity and the beam energy. This is in excellent agreement with the theoretical expectation.

Measurements of the suppression and correlations of dijets is performed using 3 $\mu$b$^{-1}$ of Xe+Xe data at $\sqrt{s_{\mathrm{NN}}} = 5.44$ TeV collected with the ATLAS detector at the LHC. Dijets with jets reconstructed using the $R=0.4$ anti-$k_t$ algorithm are measured differentially in jet $p_{\text{T}}$ over the range of 32 GeV to 398 GeV and the centrality of the collisions. Significant dijet momentum imbalance is found in the most central Xe+Xe collisions, which decreases in more peripheral collisions. Results from the measurement of per-pair normalized and absolutely normalized dijet $p_{\text{T}}$ balance are compared with previous Pb+Pb measurements at $\sqrt{s_{\mathrm{NN}}} =5.02$ TeV. The differences between the dijet suppression in Xe+Xe and Pb+Pb are further quantified by the ratio of pair nuclear-modification factors. The results are found to be consistent with those measured in Pb+Pb data when compared in classes of the same event activity and when taking into account the difference between the center-of-mass energies of the initial parton scattering process in Xe+Xe and Pb+Pb collisions. These results should provide input for a better understanding of the role of energy density, system size, path length, and fluctuations in the parton energy loss.

This Letter reports the observation of single top quarks produced together with a photon, which directly probes the electroweak coupling of the top quark. The analysis uses 139 fb$^{-1}$ of 13 TeV proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider. Requiring a photon with transverse momentum larger than 20 GeV and within the detector acceptance, the fiducial cross section is measured to be 688 $\pm$ 23 (stat.) $^{+75}_{-71}$ (syst.) fb, to be compared with the standard model prediction of 515 $^{+36}_{-42}$ fb at next-to-leading order in QCD.

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.

The production of a $W$ boson in association with a single charm quark is studied using 140 $\mathrm{fb}^{-1}$ of $\sqrt{s} = 13\,\mathrm{TeV}$ proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider. The charm quark is tagged by a charmed hadron, reconstructed with a secondary-vertex fit. The $W$ boson is reconstructed from an electron/muon decay and the missing transverse momentum. The mesons reconstructed are $D^{\pm} \to K^\mp \pi^\pm \pi^\pm$ and $D^{*\pm} \to D^{0} \pi^\pm \to (K^\mp \pi^\pm) \pi^\pm$, where $p_{\text{T}}(e, \mu) > 30\,\mathrm{GeV}$, $|\eta(e, \mu)| < 2.5$, $p_{\text{T}}(D) > 8\,\mathrm{GeV}$, and $|\eta(D)| < 2.2$. The integrated and normalized differential cross-sections as a function of the pseudorapidity of the lepton from the $W$ boson decay, and of the transverse momentum of the meson, are extracted from the data using a profile likelihood fit. The measured fiducial cross-sections are $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{-}{+}D^{+}) = 50.2\pm0.2\,\mathrm{(stat.)}\,^{+2.4}_{-2.3}\,\mathrm{(syst.)}\,\mathrm{pb}$, $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{+}{+}D^{-}) = 48.5\pm0.2\,\mathrm{(stat.)}\,^{+2.3}_{-2.2}\,\mathrm{(syst.)}\,\mathrm{pb}$, $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{-}{+}D^{*+}) = 51.1\pm0.4\,\mathrm{(stat.)}\,^{+1.9}_{-1.8}\,\mathrm{(syst.)}\,\mathrm{pb}$, and $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{+}{+}D^{*-}) = 50.0\pm0.4\,\mathrm{(stat.)}\,^{+1.9}_{-1.8}\,\mathrm{(syst.)}\,\mathrm{pb}$. Results are compared with the predictions of next-to-leading-order quantum chromodynamics calculations performed using state-of-the-art parton distribution functions. The ratio of charm to anti-charm production cross-sections is studied to probe the $s$-$\bar{s}$ quark asymmetry and is found to be $R_c^\pm = 0.971\pm0.006\,\mathrm{(stat.)}\pm0.011\,\mathrm{(syst.)}$.

A search for long-lived particles decaying into hadrons is presented. The analysis uses 139 fb$^{-1}$ of $pp$ collision data collected at $\sqrt{s} = 13$ TeV by the ATLAS detector at the LHC using events that contain multiple energetic jets and a displaced vertex. The search employs dedicated reconstruction techniques that significantly increase the sensitivity to long-lived particles decaying in the ATLAS inner detector. Background estimates for Standard Model processes and instrumental effects are extracted from data. The observed event yields are compatible with those expected from background processes. The results are used to set limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model, and on scenarios with pair-production of supersymmetric particles with long-lived electroweakinos that decay via a small $R$-parity-violating coupling. The pair-production of electroweakinos with masses below 1.5 TeV is excluded for mean proper lifetimes in the range from 0.03 ns to 1 ns. When produced in the decay of $m(\tilde{g})=2.4$ TeV gluinos, electroweakinos with $m(\tilde\chi^0_1)=1.5$ TeV are excluded with lifetimes in the range of 0.02 ns to 4 ns.

A search for flavor-changing neutral-current couplings between a top quark, an up or charm quark and a $Z$ boson is presented, using proton-proton collision data at $\sqrt{s} = 13$ TeV collected by the ATLAS detector at the Large Hadron Collider. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets both single-top-quark events produced as $gq\rightarrow tZ$ (with $q = u, c$) and top-quark-pair events, with one top quark decaying through the $t \rightarrow Zq$ channel. The analysis considers events with three leptons (electrons or muons), a $b$-tagged jet, possible additional jets, and missing transverse momentum. The data are found to be consistent with the background-only hypothesis and 95% confidence-level limits on the $t \rightarrow Zq$ branching ratios are set, assuming only tensor operators of the Standard Model effective field theory framework contribute to the $tZq$ vertices. These are $6.2 \times 10^{-5}$ ($13\times 10^{-5}$) for $t\rightarrow Zu$ ($t\rightarrow Zc$) for a left-handed $tZq$ coupling, and $6.6 \times 10^{-5}$ ($12\times 10^{-5}$) in the case of a right-handed coupling. These results are interpreted as 95% CL upper limits on the strength of corresponding couplings, yielding limits for $|C_{uW}^{(13)*}|$ and $|C_{uB}^{(13)*}|$ ($|C_{uW}^{(31)}|$ and $|C_{uB}^{(31)}|$) of 0.15 (0.16), and limits for $|C_{uW}^{(23)*}|$ and $|C_{uB}^{(23)*}|$ ($|C_{uW}^{(32)}|$ and $|C_{uB}^{(32)}|$) of 0.22 (0.21), assuming a new-physics energy scale $\Lambda_\text{NP}$ of 1 TeV.

Measurements of transverse energy$-$energy correlations and their associated azimuthal asymmetries in multijet events are presented. The analysis is performed using a data sample corresponding to 139 $\mbox{fb\(^{-1}\)}$ of proton$-$proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV, collected with the ATLAS detector at the Large Hadron Collider. The measurements are presented in bins of the scalar sum of the transverse momenta of the two leading jets and unfolded to particle level. They are then compared to next-to-next-to-leading-order perturbative QCD calculations for the first time, which feature a significant reduction in the theoretical uncertainties estimated using variations of the renormalisation and factorisation scales. The agreement between data and theory is good, thus providing a precision test of QCD at large momentum transfers $Q$. The strong coupling constant $\alpha_s$ is extracted differentially as a function of $Q$, showing a good agreement with the renormalisation group equation and with previous analyses. A simultaneous fit to all transverse energy$-$energy correlation distributions across different kinematic regions yields a value of $\alpha_\mathrm{s}(m_Z) = 0.1175 \pm 0.0006 \mbox{ (exp.)} ^{+0.0034}_{-0.0017} \mbox{ (theo.)}$, while the global fit to the asymmetry distributions yields $\alpha_{\mathrm{s}}(m_Z) = 0.1185 \pm 0.0009 \mbox{ (exp.)} ^{+0.0025}_{-0.0012} \mbox{ (theo.)}$.

This paper presents a search for a new Z' vector gauge boson with the ATLAS experiment at the Large Hadron Collider using pp collision data collected at $\sqrt{s} = 13$ TeV, corresponding to an integrated luminosity of 139 fb$^{-1}$. The new gauge boson Z' is predicted by $L_{\mu}-L_{\tau}$ models to address observed phenomena that can not be explained by the Standard Model. The search examines the four-muon (4$\mu$) final state, using a deep learning neural network classifier to separate the Z' signal from the Standard Model background events. The di-muon invariant masses in the $4\mu$ events are used to extract the Z' resonance signature. No significant excess of events is observed over the predicted background. Upper limits at a 95% confidence level on the Z' production cross-section times the decay branching fraction of $pp \rightarrow Z'\mu\mu \rightarrow 4\mu$ are set from 0.31 to 4.3 fb for the Z' mass ranging from 5 to 81 GeV. The corresponding common coupling strengths, $g_{Z'}$, of the Z' boson to the second and third generation leptons above 0.003 - 0.2 have been excluded.

For the first time at LHC energies, the forward rapidity gap spectra from proton-lead collisions for both proton and lead dissociation processes are presented. The analysis is performed over 10.4 units of pseudorapidity at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV, almost 300 times higher than in previous measurements of diffractive production in proton-nucleus collisions. For lead dissociation processes, which correspond to the pomeron-lead event topology, the EPOS-LHC generator predictions are a factor of two below the data, but the model gives a reasonable description of the rapidity gap spectrum shape. For the pomeron-proton topology, the EPOS-LHC, QGSJET II, and HIJING predictions are all at least a factor of five lower than the data. The latter effect might be explained by a significant contribution of ultra-peripheral photoproduction events mimicking the signature of diffractive processes. These data may be of significant help in understanding the high energy limit of quantum chromodynamics and for modeling cosmic ray air showers.

A search for the leptonic charge asymmetry ($A_\text{c}^{\ell}$) of top-quark$-$antiquark pair production in association with a $W$ boson ($t\bar{t}W$) is presented. The search is performed using final states with exactly three charged light leptons (electrons or muons) and is based on $\sqrt{s} = 13$ TeV proton$-$proton collision data collected with the ATLAS detector at the Large Hadron Collider at CERN during the years 2015$-$2018, corresponding to an integrated luminosity of 139 fb$^{-1}$. A profile-likelihood fit to the event yields in multiple regions corresponding to positive and negative differences between the pseudorapidities of the charged leptons from top-quark and top-antiquark decays is used to extract the charge asymmetry. At reconstruction level, the asymmetry is found to be $-0.123 \pm 0.136$ (stat.) $\pm \, 0.051$ (syst.). An unfolding procedure is applied to convert the result at reconstruction level into a charge-asymmetry value in a fiducial volume at particle level with the result of $-0.112 \pm 0.170$ (stat.) $\pm \, 0.054$ (syst.). The Standard Model expectations for these two observables are calculated using Monte Carlo simulations with next-to-leading-order plus parton shower precision in quantum chromodynamics and including next-to-leading-order electroweak corrections. They are $-0.084 \, ^{+0.005}_{-0.003}$ (scale) $\pm\, 0.006$ (MC stat.) and $-0.063 \, ^{+0.007}_{-0.004}$ (scale) $\pm\, 0.004$ (MC stat.) respectively, and in agreement with the measurements.

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 search for long-lived particles decaying in the outer regions of the CMS silicon tracker or in the calorimeters is presented. The search is based on a data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018, corresponding to an integrated luminosity of 138 fb$^{-1}$. A novel technique, using trackless and out-of-time jet information combined in a deep neural network discriminator, is employed to identify decays of long-lived particles. The results are interpreted in a simplified model of chargino-neutralino production, where the neutralino is the next-to-lightest supersymmetric particle, is long-lived, and decays to a gravitino and either a Higgs or Z boson. This search is most sensitive to neutralino proper decay lengths of approximately 0.5 m, for which masses up to 1.18 TeV are excluded at 95% confidence level. The current search is the best result to date in the mass range from the kinematic limit imposed by the Higgs mass up to 1.8 TeV.

This paper presents measurements of charged-hadron spectra obtained in $pp$, $p$+Pb, and Pb+Pb collisions at $\sqrt{s}$ or $\sqrt{s_{_\text{NN}}}=5.02$ TeV, and in Xe+Xe collisions at $\sqrt{s_{_\text{NN}}}=5.44$ TeV. The data recorded by the ATLAS detector at the LHC have total integrated luminosities of 25 pb${}^{-1}$, 28 nb${}^{-1}$, 0.50 nb${}^{-1}$, and 3 $\mu$b${}^{-1}$, respectively. The nuclear modification factors $R_{p\text{Pb}}$ and $R_\text{AA}$ are obtained by comparing the spectra in heavy-ion and $pp$ collisions in a wide range of charged-particle transverse momenta and pseudorapidity. The nuclear modification factor $R_{p\text{Pb}}$ shows a moderate enhancement above unity with a maximum at $p_{\mathrm{T}} \approx 3$ GeV; the enhancement is stronger in the Pb-going direction. The nuclear modification factors in both Pb+Pb and Xe+Xe collisions feature a significant, centrality-dependent suppression. They show a similar distinct $p_{\mathrm{T}}$-dependence with a local maximum at $p_{\mathrm{T}} \approx 2$ GeV and a local minimum at $p_{\mathrm{T}} \approx 7$ GeV. This dependence is more distinguishable in more central collisions. No significant $|\eta|$-dependence is found. A comprehensive comparison with several theoretical predictions is also provided. They typically describe $R_\text{AA}$ better in central collisions and in the $p_{\mathrm{T}}$ range from about 10 to 100 GeV.

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.

This paper presents a statistical combination of searches targeting final states with two top quarks and invisible particles, characterised by the presence of zero, one or two leptons, at least one jet originating from a $b$-quark and missing transverse momentum. The analyses are searches for phenomena beyond the Standard Model consistent with the direct production of dark matter in $pp$ collisions at the LHC, using 139 fb$^{-\text{1}}$ of data collected with the ATLAS detector at a centre-of-mass energy of 13 TeV. The results are interpreted in terms of simplified dark matter models with a spin-0 scalar or pseudoscalar mediator particle. In addition, the results are interpreted in terms of upper limits on the Higgs boson invisible branching ratio, where the Higgs boson is produced according to the Standard Model in association with a pair of top quarks. For scalar (pseudoscalar) dark matter models, with all couplings set to unity, the statistical combination extends the mass range excluded by the best of the individual channels by 50 (25) GeV, excluding mediator masses up to 370 GeV. In addition, the statistical combination improves the expected coupling exclusion reach by 14% (24%), assuming a scalar (pseudoscalar) mediator mass of 10 GeV. An upper limit on the Higgs boson invisible branching ratio of 0.38 (0.30$^{+\text{0.13}}_{-\text{0.09}}$) is observed (expected) at 95% confidence level.

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.

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

A measurement of observables sensitive to effects of colour reconnection in top-quark pair-production events is presented using 139 fb$^{-1}$ of 13$\,$TeV proton-proton collision data collected by the ATLAS detector at the LHC. Events are selected by requiring exactly one isolated electron and one isolated muon with opposite charge and two or three jets, where exactly two jets are required to be $b$-tagged. For the selected events, measurements are presented for the charged-particle multiplicity, the scalar sum of the transverse momenta of the charged particles, and the same scalar sum in bins of charged-particle multiplicity. These observables are unfolded to the stable-particle level, thereby correcting for migration effects due to finite detector resolution, acceptance and efficiency effects. The particle-level measurements are compared with different colour reconnection models in Monte Carlo generators. These measurements disfavour some of the colour reconnection models and provide inputs to future optimisation of the parameters in Monte Carlo generators.

Three searches are presented for signatures of physics beyond the standard model (SM) in $\tau\tau$ final states in proton-proton collisions at the LHC, using a data sample collected with the CMS detector at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. Upper limits at 95% confidence level (CL) are set on the products of the branching fraction for the decay into $\tau$ leptons and the cross sections for the production of a new boson $\phi$, in addition to the H(125) boson, via gluon fusion (gg$\phi$) or in association with b quarks, ranging from $\mathcal{O}$(10 pb) for a mass of 60 GeV to 0.3 fb for a mass of 3.5 TeV each. The data reveal two excesses for gg$\phi$ production with local $p$-values equivalent to about three standard deviations at $m_\phi$ = 0.1 and 1.2 TeV. In a search for $t$-channel exchange of a vector leptoquark U$_1$, 95% CL upper limits are set on the dimensionless U$_1$ leptoquark coupling to quarks and $\tau$ leptons ranging from 1 for a mass of 1 TeV to 6 for a mass of 5 TeV, depending on the scenario. In the interpretations of the $M_\mathrm{h}^{125}$ and $M_\mathrm{h, EFT}^{125}$ minimal supersymmetric SM benchmark scenarios, additional Higgs bosons with masses below 350 GeV are excluded at 95% CL.