This Letter presents direct searches for lepton flavour violation in Higgs boson decays, $H\rightarrow e\tau$ and $H\rightarrow\mu\tau$, performed with the ATLAS detector at the LHC. The searches are based on a data sample of proton-proton collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV, corresponding to an integrated luminosity of $36.1\,\mathrm{fb}^{-1}$. No significant excess is observed above the expected background from Standard Model processes. The observed (median expected) 95 % confidence-level upper limits on the lepton-flavour-violating branching ratios are $0.47\%$ ($0.34^{+0.13}_{-0.10}\,\%$) and $0.28\%$ ($0.37^{+0.14}_{-0.10}\,\%$) for $H\to e\tau$ and $H\to\mu\tau$, respectively.

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.

This paper presents a search for pair production of higgsinos, the supersymmetric partners of the Higgs bosons, in scenarios with gauge-mediated supersymmetry breaking. Each higgsino is assumed to decay into a Higgs boson and a nearly massless gravitino. The search targets events where each Higgs boson decays into $b\bar{b}$, leading to a reconstructed final state with at least three energetic $b$-jets and missing transverse momentum. Two complementary analysis channels are used, with each channel specifically targeting either low or high values of the higgsino mass. The low-mass (high-mass) channel exploits 126 (139) fb$^{-1}$ of $\sqrt{s}=13$ TeV data collected by the ATLAS detector during Run 2 of the Large Hadron Collider. No significant excess above the Standard Model prediction is found. At 95% confidence level, masses between 130 GeV and 940 GeV are excluded for higgsinos decaying exclusively into Higgs bosons and gravitinos. Exclusion limits as a function of the higgsino decay branching ratio to a Higgs boson are also reported.

A search for the dimuon decay of the Higgs boson was performed using data corresponding to an integrated luminosity of 36.1 fb$^{-1}$ collected with the ATLAS detector in $pp$ collisions at $\sqrt{s}=13$ TeV at the Large Hadron Collider. No significant excess is observed above the expected background. The observed (expected) upper limit on the cross section times branching ratio is 3.0 (3.1) times the Standard Model prediction at the 95% confidence level for a Higgs boson mass of 125 GeV. When combined with the $pp$ collision data at $\sqrt{s}=7$ TeV and $\sqrt{s}=8$ TeV, the observed (expected) upper limit is 2.8 (2.9) times the Standard Model prediction.

A search for heavy neutral Higgs bosons is performed using the LHC Run 2 data, corresponding to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}=13$ TeV recorded with the ATLAS detector. The search for heavy resonances is performed over the mass range 0.2-2.5 TeV for the $\tau^+\tau^-$ decay with at least one $\tau$-lepton decaying into final states with hadrons. The data are in good agreement with the background prediction of the Standard Model. In the $M_{h}^{125}$ scenario of the Minimal Supersymmetric Standard Model, values of $\tan\beta>8$ and $\tan\beta>21$ are excluded at the 95% confidence level for neutral Higgs boson masses of 1.0 TeV and 1.5 TeV, respectively, where $\tan\beta$ is the ratio of the vacuum expectation values of the two Higgs doublets.

A search for pair production of the supersymmetric partners of the Higgs boson (higgsinos $\tilde{H}$) in gauge-mediated scenarios is reported. Each higgsino is assumed to decay to a Higgs boson and a gravitino. Two complementary analyses, targeting high- and low-mass signals, are performed to maximize sensitivity. The two analyses utilize LHC $pp$ collision data at a center-of-mass energy $\sqrt{s} = 13$ TeV, the former with an integrated luminosity of 36.1 fb$^{-1}$ and the latter with 24.3 fb$^{-1}$, collected with the ATLAS detector in 2015 and 2016. The search is performed in events containing missing transverse momentum and several energetic jets, at least three of which must be identified as $b$-quark jets. No significant excess is found above the predicted background. Limits on the cross-section are set as a function of the mass of the $\tilde{H}$ in simplified models assuming production via mass-degenerate higgsinos decaying to a Higgs boson and a gravitino. Higgsinos with masses between 130 and 230 GeV and between 290 and 880 GeV are excluded at the 95% confidence level. Interpretations of the limits in terms of the branching ratio of the higgsino to a $Z$ boson or a Higgs boson are also presented, and a 45% branching ratio to a Higgs boson is excluded for $m_{\tilde{H}} \approx 400$ 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)$.

A search for the $Z\gamma$ decay of the Higgs boson, with $Z$ boson decays into pairs of electrons or muons is presented. The analysis uses proton$-$proton collision data at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS detector at the Large Hadron Collider. The observed data are consistent with the expected background with a $p$-value of 1.3%. An upper limit at 95% confidence level on the production cross-section times the branching ratio for $pp\to H\to Z\gamma$ is set at 3.6 times the Standard Model prediction while 2.6 times is expected in the presence of the Standard Model Higgs boson. The best-fit value for the signal yield normalised to the Standard Model prediction is $2.0^{+1.0}_{-0.9}$ where the statistical component of the uncertainty is dominant.

A search for Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is carried out with up to 36.1 $\mathrm{fb}^{-1}$ of LHC proton--proton collision data collected at $\sqrt{s}$ = 13 TeV with the ATLAS detector in 2015 and 2016. Three benchmark signals are studied: a spin-2 graviton decaying into a Higgs boson pair, a scalar resonance decaying into a Higgs boson pair, and Standard Model non-resonant Higgs boson pair production. Two analyses are carried out, each implementing a particular technique for the event reconstruction that targets Higgs bosons reconstructed as pairs of jets or single boosted jets. The resonance mass range covered is 260--3000 GeV. The analyses are statistically combined and upper limits on the production cross section of Higgs boson pairs times branching ratio to $b\bar{b}b\bar{b}$ are set in each model. No significant excess is observed; the largest deviation of data over prediction is found at a mass of 280 GeV, corresponding to 2.3 standard deviations globally. The observed 95% confidence level upper limit on the non-resonant production is 13 times the Standard Model prediction.

Measurements of the Standard Model Higgs boson decaying into a $b\bar{b}$ pair and produced in association with a $W$ or $Z$ boson decaying into leptons, using proton-proton collision data collected between 2015 and 2018 by the ATLAS detector, are presented. The measurements use collisions produced by the Large Hadron Collider at a centre-of-mass energy of $\sqrt{s} = $13 TeV, corresponding to an integrated luminosity of 139 fb$^{-1}$. The production of a Higgs boson in association with a $W$ or $Z$ boson is established with observed (expected) significances of 4.0 (4.1) and 5.3 (5.1) standard deviations, respectively. Cross-sections of associated production of a Higgs boson decaying into bottom quark pairs with an electroweak gauge boson, $W$ or $Z$, decaying into leptons are measured as a function of the gauge boson transverse momentum in kinematic fiducial volumes. The cross-section measurements are all consistent with the Standard Model expectations, and the total uncertainties vary from 30% in the high gauge boson transverse momentum regions to 85% in the low regions. Limits are subsequently set on the parameters of an effective Lagrangian sensitive to modifications of the $WH$ and $ZH$ processes as well as the Higgs boson decay into $b\bar{b}$.