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

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

A search for the Higgs boson decaying into a photon and a pair of electrons or muons with an invariant mass $m_{\ell\ell} < 30$ GeV is presented. The analysis is performed using 139 fb$^{-1}$ of proton-proton collision data, produced by the LHC at a centre-of-mass energy of 13 TeV and collected by the ATLAS experiment. Evidence for the $H \rightarrow \ell \ell \gamma$ process is found with a significance of 3.2$\sigma$ over the background-only hypothesis, compared to an expected significance of 2.1$\sigma$. The best-fit value of the signal strength parameter, defined as the ratio of the observed signal yield to the one expected in the Standard Model, is $\mu = 1.5 \pm 0.5$. The Higgs boson production cross-section times the $H \rightarrow\ell\ell\gamma$ branching ratio for $m_{\ell\ell} <$ 30 GeV is determined to be 8.7 $^{+2.8}_{-2.7}$ fb.

This paper presents direct searches for lepton flavour violation in Higgs boson decays, $H\rightarrow e\tau$ and $H\rightarrow\mu\tau$, performed using data collected 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 138 fb$^{-1}$. Leptonic ($\tau \rightarrow \ell \nu_\ell \nu_\tau$) and hadronic ($\tau \rightarrow $ hadrons $ \nu_\tau$) decays of the $\tau$-lepton are considered. Two background estimation techniques are employed: the MC-template method, based on data-corrected simulation samples, and the Symmetry method, based on exploiting the symmetry between electrons and muons in the Standard Model backgrounds. No significant excess of events is observed and the results are interpreted as upper limits on lepton-flavour-violating branching ratios of the Higgs boson. The observed (expected) upper limits set on the branching ratios at 95% confidence level, $\mathcal{B}(H\rightarrow e\tau)<0.20\%$ (0.12%) and $\mathcal{B}(H\rightarrow \mu\tau)<0.18\%$ (0.09%), are obtained with the MC-template method from a simultaneous measurement of potential $H \rightarrow e\tau$ and $H \rightarrow\mu\tau$ signals. The best-fit branching ratio difference, $\mathcal{B}(H\rightarrow \mu\tau)- \mathcal{B}(H\rightarrow e\tau)$, measured with the Symmetry method in the channel where the $\tau$-lepton decays to leptons, is (0.25 $\pm$ 0.10)%, compatible with a value of zero within 2.5$\sigma$.

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.

Measurements of the production cross-sections of the Standard Model (SM) Higgs boson ($H$) decaying into a pair of $\tau$-leptons are presented. The measurements use data collected with the ATLAS detector from $pp$ collisions produced at the Large Hadron Collider at a centre-of-mass energy of $\sqrt{s}=13\,\text{TeV}$, corresponding to an integrated luminosity of $139\,\text{fb}^{-1}$. Leptonic ($\tau\to\ell\nu_{\ell}\nu_{\tau}$) and hadronic ($\tau\to\text{hadrons}~\nu_{\tau}$) decays of the $\tau$-lepton are considered. All measurements account for the branching ratio of $H\to\tau\tau$ and are performed with a requirement $|y_H|<2.5$, where $y_H$ is the true Higgs boson rapidity. The cross-section of the $pp\to H\to\tau\tau$ process is measured to be $2.94 \pm 0.21 \text{(stat)} ^{+\,0.37}_{-\,0.32} \text{(syst)}$ pb, in agreement with the SM prediction of $3.17\pm0.09~ \mbox{pb}$. Inclusive cross-sections are determined separately for the four dominant production modes: $2.65 \pm 0.41 \text{(stat)} ^{+\,0.91}_{-\,0.67} \text{(syst)}$ pb for gluon$-$gluon fusion, $0.197 \pm 0.028 \text{(stat)} ^{+\,0.032}_{-\,0.026} \text{(syst)}$ pb for vector-boson fusion, $0.115 \pm 0.058 \text{(stat)} ^{+\,0.042}_{-\,0.040} \text{(syst)}$ pb for vector-boson associated production, and $0.033 \pm 0.031 \text{(stat)} ^{+\,0.022}_{-\,0.017} \text{(syst)}$ pb for top-quark pair associated production. Measurements in exclusive regions of the phase space, using the simplified template cross-section framework, are also performed. All results are in agreement with the SM predictions.

A search is presented for flavour-changing neutral-current interactions involving the top quark, the Higgs boson and an up-type quark ($q=u,c$) with the ATLAS detector at the Large Hadron Collider. The analysis considers leptonic decays of the top quark along with Higgs boson decays into two $W$ bosons, two $Z$ bosons or a $\tau^{+}\tau^{-}$ pair. It focuses on final states containing either two leptons (electrons or muons) of the same charge or three leptons. The considered processes are $t\bar{t}$ and $Ht$ production. For the $t\bar{t}$ production, one top quark decays via $t\to Hq$. The proton-proton collision data set analysed amounts to 140 fb$^{-1}$ at $\sqrt{s}=13$ TeV. No significant excess beyond Standard Model expectations is observed and upper limits are set on the $t\to Hq$ branching ratios at 95% confidence level, amounting to observed (expected) limits of $\mathcal{B}(t\to Hu)<2.8\,(3.0) \times 10^{-4}$ and $\mathcal{B}(t\to Hc)<3.3\,(3.8) \times 10^{-4}$. Combining this search with other searches for $tHq$ flavour-changing neutral-current interactions previously conducted by ATLAS, considering $H\to b\bar{b}$ and $H\to\gamma\gamma$ decays, as well as $H\to\tau^{+}\tau^{-}$ decays with one or two hadronically decaying $\tau$-leptons, yields observed (expected) upper limits on the branching ratios of $\mathcal{B}(t\to Hu)<2.6\,(1.8) \times 10^{-4}$ and $\mathcal{B}(t\to Hc)<3.4\,(2.3) \times 10^{-4}$.

This paper describes a search for beyond the Standard Model decays of the Higgs boson into a pair of new spin-0 particles subsequently decaying into $b$-quark pairs, $H \rightarrow aa \rightarrow (b\bar{b})(b\bar{b})$, using proton-proton collision data collected by the ATLAS detector at the Large Hadron Collider at center-of-mass energy $\sqrt{s}=13$ TeV. This search focuses on the regime where the decay products are collimated and in the range $15 \leq m_a \leq 30$ GeV and is complementary to a previous search in the same final state targeting the regime where the decay products are well separated and in the range $20 \leq m_a \leq 60$ GeV. A novel strategy for the identification of the $a \rightarrow b\bar{b}$ decays is deployed to enhance the efficiency for topologies with small separation angles. The search is performed with 36 fb$^{-1}$ of integrated luminosity collected in 2015 and 2016 and sets upper limits on the production cross-section of $H \rightarrow aa \rightarrow (b\bar{b})(b\bar{b})$, where the Higgs boson is produced in association with a $Z$ boson.