A search for a new heavy scalar or pseudo-scalar Higgs boson ($H/A$) produced in association with a pair of top quarks, with the Higgs boson decaying into a pair of top quarks ($H/A\rightarrow t\bar{t}$) is reported. The search targets a final state with exactly two leptons with same-sign electric charges or at least three leptons. The analysed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collisions collected at a centre-of-mass energy of 13 TeV with the ATLAS detector at the LHC. Two multivariate classifiers are used to separate the signal from the background. No significant excess of events over the Standard Model expectation is observed. The results are interpreted in the context of a type-II two-Higgs-doublet model. The observed (expected) upper limits at 95% confidence level on the $t\bar{t}H/A$ production cross-section times the branching ratio of $H/A\rightarrow t\bar{t}$ range between 14 (10) fb and 6 (5) fb for a heavy Higgs boson with mass between 400 GeV and 1000 GeV, respectively. Assuming that only one particle, either the scalar $H$ or the pseudo-scalar $A$, contributes to the $t\bar{t}t\bar{t}$ final state, values of $\tan\beta$ below 1.2 or 0.5 are excluded for a mass of 400 GeV or 1000 GeV, respectively. These exclusion ranges increase to $\tan\beta$ below 1.6 or 0.6 when both particles are considered.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

A search for nonresonant Higgs boson pair production in the $b\bar{b}\gamma\gamma$ final state is performed using 140 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. This analysis supersedes and expands upon the previous nonresonant ATLAS results in this final state based on the same data sample. The analysis strategy is optimised to probe anomalous values not only of the Higgs ($H$) boson self-coupling modifier $\kappa_\lambda$ but also of the quartic $HHVV$ ($V=W,Z$) coupling modifier $\kappa_{2V}$. No significant excess above the expected background from Standard Model processes is observed. An observed upper limit $\mu_{HH}<4.0$ is set at 95% confidence level on the Higgs boson pair production cross-section normalised to its Standard Model prediction. The 95% confidence intervals for the coupling modifiers are $-1.4<\kappa_\lambda<6.9$ and $-0.5<\kappa_{2V}<2.7$, assuming all other Higgs boson couplings except the one under study are fixed to the Standard Model predictions. The results are interpreted in the Standard Model effective field theory and Higgs effective field theory frameworks in terms of constraints on the couplings of anomalous Higgs boson (self-)interactions.

A search is presented for a heavy resonance $Y$ decaying into a Standard Model Higgs boson $H$ and a new particle $X$ in a fully hadronic final state. The full Large Hadron Collider Run 2 dataset of proton-proton collisions at $\sqrt{s}= 13$ TeV collected by the ATLAS detector from 2015 to 2018 is used, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets the high $Y$-mass region, where the $H$ and $X$ have a significant Lorentz boost in the laboratory frame. A novel signal region is implemented using anomaly detection, where events are selected solely because of their incompatibility with a learned background-only model. It is defined using a jet-level tagger for signal-model-independent selection of the boosted $X$ particle, representing the first application of fully unsupervised machine learning to an ATLAS analysis. Two additional signal regions are implemented to target a benchmark $X$ decay into two quarks, covering topologies where the $X$ is reconstructed as either a single large-radius jet or two small-radius jets. The analysis selects Higgs boson decays into $b\bar{b}$, and a dedicated neural-network-based tagger provides sensitivity to the boosted heavy-flavor topology. No significant excess of data over the expected background is observed, and the results are presented as upper limits on the production cross section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) for signals with $m_Y$ between 1.5 and 6 TeV and $m_X$ between 65 and 3000 GeV.

A search for supersymmetry involving the pair production of gluinos decaying via off-shell third-generation squarks into the lightest neutralino ($\tilde\chi^0_1$) is reported. It exploits LHC proton$-$proton collision data at a centre-of-mass energy $\sqrt{s} = 13$ TeV with an integrated luminosity of 139 fb$^{-1}$ collected with the ATLAS detector from 2015 to 2018. The search uses events containing large missing transverse momentum, up to one electron or muon, and several energetic jets, at least three of which must be identified as containing $b$-hadrons. Both a simple kinematic event selection and an event selection based upon a deep neural-network are used. No significant excess above the predicted background is found. In simplified models involving the pair production of gluinos that decay via off-shell top (bottom) squarks, gluino masses less than 2.44 TeV (2.35 TeV) are excluded at 95% CL for a massless $\tilde\chi^0_1$. Limits are also set on the gluino mass in models with variable branching ratios for gluino decays to $b\bar{b}\tilde\chi^0_1$, $t\bar{t}\tilde\chi^0_1$ and $t\bar{b}\tilde\chi^-_1$ / $\bar{t}b\tilde\chi^+_1$.

This Letter presents the measurement of the fiducial and differential cross-sections of the electroweak production of a $Z \gamma$ pair in association with two jets. The analysis uses 140 fb$^{-1}$ of LHC proton-proton collision data taken at $\sqrt{s}$=13 TeV recorded by the ATLAS detector during the years 2015-2018. Events with a $Z$ boson candidate decaying into either an $e^+e^-$ or $\mu^+ \mu^-$ pair, a photon and two jets are selected. The electroweak component is extracted by requiring a large dijet invariant mass and a large rapidity gap between the two jets and is measured with an observed and expected significance well above five standard deviations. The fiducial $pp \rightarrow Z \gamma jj$ cross-section for the electroweak production is measured to be 3.6 $\pm$ 0.5 fb. The total fiducial cross-section that also includes contributions where the jets arise from strong interactions is measured to be $16.8^{+2.0}_{-1.8}$ fb. The results are consistent with the Standard Model predictions. Differential cross-sections are also measured using the same events and are compared with parton-shower Monte Carlo simulations. Good agreement is observed between data and predictions.

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.

Searches for the exclusive decays of the Higgs boson into $D^*\gamma$ and of the $Z$ boson into $D^0\gamma$ and $K^0_s\gamma$ can probe flavour-violating Higgs and $Z$ boson couplings to light quarks. Searches for these decays are performed with a $pp$ collision data sample corresponding to an integrated luminosity of $136.3$ fb$^{-1}$ collected at $\sqrt{s}=13$ TeV between 2016-2018 with the ATLAS detector at the CERN Large Hadron Collider. In the $D^*\gamma$ and $D^0\gamma$ channels, the observed (expected) 95$\%$ confidence-level upper limits on the respective branching fractions are ${\cal B}(H\rightarrow D^*\gamma)< 1.0 (1.2)\times 10^{-3}$, ${\cal B}(Z\rightarrow D^0\gamma)< 4.0 (3.4)\times 10^{-6}$, while the corresponding results in the $K^0_s\gamma$ channel are ${\cal B}(Z\rightarrow K^0_s\gamma)< 3.1 (3.0)\times 10^{-6}$.

A search for the electroweak production of pairs of charged sleptons or charginos decaying into two-lepton final states with missing transverse momentum is presented. Two simplified models of $R$-parity-conserving supersymmetry are considered: direct pair-production of sleptons ($\tilde{\ell}\tilde{\ell}$), with each decaying into a charged lepton and a $\tilde{\chi}_1^0$ neutralino, and direct pair-production of the lightest charginos $(\tilde{\chi}_1^\pm\tilde{\chi}_1^\mp)$, with each decaying into a $W$-boson and a $\tilde{\chi}_1^0$. The lightest neutralino ($\tilde{\chi}_1^0$) is assumed to be the lightest supersymmetric particle (LSP). The analyses target the experimentally challenging mass regions where $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and $m(\tilde{\chi}_1^\pm)-m(\tilde{\chi}_1^0)$ are close to the $W$-boson mass (`moderately compressed' regions). The search uses 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collisions recorded by the ATLAS detector at the Large Hadron Collider. No significant excesses over the expected background are observed. Exclusion limits on the simplified models under study are reported in the ($\tilde{\ell},\tilde{\chi}_1^0$) and ($\tilde{\chi}_1^\pm,\tilde{\chi}_1^0$) mass planes at 95% confidence level (CL). Sleptons with masses up to 150 GeV are excluded at 95% CL for the case of a mass-splitting between sleptons and the LSP of 50 GeV. Chargino masses up to 140 GeV are excluded at 95% CL for the case of a mass-splitting between the chargino and the LSP down to about 100 GeV.

A search for flavour-changing neutral current (FCNC) $tqH$ interactions involving a top quark, another up-type quark ($q=u$, $c$), and a Standard Model (SM) Higgs boson decaying into a $\tau$-lepton pair ($H\rightarrow \tau^+\tau^-$) is presented. The search is based on a dataset of $pp$ collisions at $\sqrt{s}=13$ TeV that corresponds to an integrated luminosity of 139 fb$^{-1}$ recorded with the ATLAS detector at the Large Hadron Collider. Two processes are considered: single top quark FCNC production in association with a Higgs boson ($pp\rightarrow tH$), and top quark pair production in which one of the top quarks decays into $Wb$ and the other decays into $qH$ through the FCNC interactions. The search selects events with two hadronically decaying $\tau$-lepton candidates ($\tau_{\text{had}}$) or at least one $\tau_{\text{had}}$ with an additional lepton ($e$, $\mu$), as well as multiple jets. Event kinematics is used to separate signal from the background through a multivariate discriminant. A slight excess of data is observed with a significance of 2.3$\sigma$ above the expected SM background, and 95% CL upper limits on the $t\to qH$ branching ratios are derived. The observed (expected) 95% CL upper limits set on the $t\to cH$ and $t\to uH$ branching ratios are $9.4 \times 10^{-4}$ $(4.8^{+2.2}_{-1.4}\times 10^{-4})$ and $6.9\times 10^{-4}$ $(3.5^{+1.5}_{-1.0}\times 10^{-4})$, respectively. The corresponding combined observed (expected) upper limits on the dimension-6 operator Wilson coefficients in the effective $tqH$ couplings are $C_{c\phi} <1.35$ $(0.97)$ and $C_{u\phi} <1.16$ $(0.82)$.