This paper presents a search for decays of the Higgs boson with a mass of 125 GeV into a pair of new pseudoscalar particles, $H\rightarrow aa$, where one $a$-boson decays into a $b$-quark pair and the other into a muon pair. The search uses 139 fb$^{-1}$ of proton-proton collision data at a center-of-mass energy of $\sqrt{s}=13$ TeV recorded between 2015 and 2018 by the ATLAS experiment at the LHC. A narrow dimuon resonance is searched for in the invariant mass spectrum between 16 GeV and 62 GeV. The largest excess of events above the Standard Model backgrounds is observed at a dimuon invariant mass of 52 GeV and corresponds to a local (global) significance of $3.3 \sigma$ ($1.7 \sigma$). Upper limits at 95% confidence level are placed on the branching ratio of the Higgs boson to the $bb\mu\mu$ final state, $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$, and are in the range $\text{(0.2-4.0)} \times 10^{-4}$, depending on the signal mass hypothesis.

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 resonant Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is presented. The analysis uses 126-139 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector at the Large Hadron Collider. The analysis is divided into two channels, targeting Higgs boson decays which are reconstructed as pairs of small-radius jets or as individual large-radius jets. Spin-0 and spin-2 benchmark signal models are considered, both of which correspond to resonant $HH$ production via gluon$-$gluon fusion. The data are consistent with Standard Model predictions. Upper limits are set on the production cross-section times branching ratio to Higgs boson pairs of a new resonance in the mass range from 251 GeV to 5 TeV.

A summary of the constraints from searches performed by the ATLAS Collaboration for the electroweak production of charginos and neutralinos is presented. Results from eight separate ATLAS searches are considered, each using 140 fb$^{-1}$ of proton-proton data at a centre-of-mass energy of $\sqrt{s}$=13 TeV collected at the Large Hadron Collider during its second data-taking run. The results are interpreted in the context of the 19-parameter phenomenological minimal supersymmetric standard model, where R-parity conservation is assumed and the lightest supersymmetric particle is assumed to be the lightest neutralino. Constraints from previous electroweak, flavour and dark matter related measurements are also considered. The results are presented in terms of constraints on supersymmetric particle masses and are compared with limits from simplified models. Also shown is the impact of ATLAS searches on parameters such as the dark matter relic density and the spin-dependent and spin-independent scattering cross-sections targeted by direct dark matter detection experiments. The Higgs boson and Z boson `funnel regions', where a low-mass neutralino would not oversaturate the dark matter relic abundance, are almost completely excluded by the considered constraints. Example spectra for non-excluded supersymmetric models with light charginos and neutralinos are also presented.

A search for the exotic decay of the Higgs boson ($H$) into a $b\bar{b}$ resonance plus missing transverse momentum is described. The search is performed with the ATLAS detector at the Large Hadron Collider using 139 $\mathrm{fb}^{-1}$ of $pp$ collisions at $\sqrt{s} = 13$ TeV. The search targets events from $ZH$ production in an NMSSM scenario where $H \rightarrow \tilde{\chi}^{0}_{2}\tilde{\chi}^{0}_{1}$, with $\tilde{\chi}^{0}_{2} \rightarrow {a} \tilde{\chi}^{0}_{1}$, where $a$ is a light pseudoscalar Higgs boson and $\tilde{\chi}^{0}_{1,2}$ are the two lightest neutralinos. The decay of the $a$ boson into a pair of $b$-quarks results in a peak in the dijet invariant mass distribution. The final-state signature consists of two leptons, two or more jets, at least one of which is identified as originating from a $b$-quark, and missing transverse momentum. Observations are consistent with Standard Model expectations and upper limits are set on the product of cross section times branching ratio for a three-dimensional scan of the masses of the $\tilde{\chi}^{0}_{2}$, $\tilde{\chi}^{0}_{1}$ and $a$ boson.

Narrow resonances decaying into $WW$, $WZ$ or $ZZ$ boson pairs are searched for in 139 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the Large Hadron Collider from 2015 to 2018. The diboson system is reconstructed using pairs of high transverse momentum, large-radius jets. These jets are built from a combination of calorimeter- and tracker-inputs compatible with the hadronic decay of a boosted $W$ or $Z$ boson, using jet mass and substructure properties. The search is performed for diboson resonances with masses greater than 1.3 TeV. No significant deviations from the background expectations are observed. Exclusion limits at the 95% confidence level are set on the production cross-section times branching ratio into dibosons for resonances in a range of theories beyond the Standard Model, with the highest excluded mass of a new gauge boson at 3.8 TeV in the context of mass-degenerate resonances that couple predominantly to gauge bosons.

A search for dark matter is conducted in final states containing a photon and missing transverse momentum in proton$-$proton collisions at $\sqrt{s}$ = 13 TeV. The data, collected during 2015$-$2018 by the ATLAS experiment at the CERN LHC, correspond to an integrated luminosity of 139 fb$^{-1}$. No deviations from the predictions of the Standard Model are observed and 95% confidence-level upper limits between 2.45 fb and 0.5 fb are set on the visible cross section for contributions from physics beyond the Standard Model, in different ranges of the missing transverse momentum. The results are interpreted as 95% confidence-level limits in models where weakly interacting dark-matter candidates are pair-produced via an s-channel axial-vector or vector mediator. Dark-matter candidates with masses up to 415 (580) GeV are excluded for axial-vector (vector) mediators, while the maximum excluded mass of the mediator is 1460 (1470) GeV. In addition, the results are expressed in terms of 95% confidence-level limits on the parameters of a model with an axion-like particle produced in association with a photon, and are used to constrain the coupling $g_{aZ\gamma}$ of an axion-like particle to the electroweak gauge bosons.

A measurement of the $B^0_s \to J/\psi\phi$ decay parameters using 80.5 $\mathrm{fb}^{-1}$ of integrated luminosity collected with the ATLAS detector from 13 TeV proton-proton collisions at the LHC is presented. The measured parameters include the $CP$-violating phase $\phi_s$, the width difference $\Delta\Gamma_{s}$ between the $B^0_s$ meson mass eigenstates and the average decay width $\Gamma_{s}$. The values measured for the physical parameters are combined with those from 19.2 $\mathrm{fb}^{-1}$ of 7 TeV and 8 TeV data, leading to the following: \begin{eqnarray*} \phi_s & = & -0.087\phantom{0} \pm 0.036\phantom{0} ~\mathrm{(stat.)} \pm 0.021\phantom{0} ~\mathrm{(syst.)~rad} \\ \Delta\Gamma_{s} & = & \phantom{-}0.0657 \pm 0.0043 ~\mathrm{(stat.)} \pm 0.0037 ~\mathrm{(syst.)~ps}^{-1} \\ \Gamma_{s} & = & \phantom{-}0.6703 \pm 0.0014 ~\mathrm{(stat.)} \pm 0.0018 ~\mathrm{(syst.)~ps}^{-1} \\ \end{eqnarray*} Results for $\phi_s$ and $\Delta\Gamma_{s}$ are also presented as 68% confidence level contours in the $\phi_s$-$\Delta\Gamma_{s}$ plane. Furthermore, the transversity amplitudes and corresponding strong phases are measured. $\phi_s$ and $\Delta\Gamma_{s}$ measurements are in agreement with the Standard Model predictions.

The production of dark matter in association with Higgs bosons is predicted in several extensions of the Standard Model. An exploration of such scenarios is presented, considering final states with missing transverse momentum and $b$-tagged jets consistent with a Higgs boson. The analysis uses proton-proton collision data at a centre-of-mass energy of 13 TeV recorded by the ATLAS experiment at the LHC during Run 2, amounting to an integrated luminosity of 139 fb$^{-1}$. The analysis, when compared with previous searches, benefits from a larger dataset, but also has further improvements providing sensitivity to a wider spectrum of signal scenarios. These improvements include both an optimised event selection and advances in the object identification, such as the use of the likelihood-based significance of the missing transverse momentum and variable-radius track-jets. No significant deviation from Standard Model expectations is observed. Limits are set, at 95% confidence level, in two benchmark models with two Higgs doublets extended by either a heavy vector boson $Z'$ or a pseudoscalar singlet $a$ and which both provide a dark matter candidate $\chi$. In the case of the two-Higgs-doublet model with an additional vector boson $Z'$, the observed limits extend up to a $Z'$ mass of 3 TeV for a mass of 100 GeV for the dark matter candidate. The two-Higgs-doublet model with a dark matter particle mass of 10 GeV and an additional pseudoscalar $a$ is excluded for masses of the $a$ up to 520 GeV and 240 GeV for $\tan \beta = 1$ and $\tan \beta = 10$ respectively. Limits on the visible cross-sections are set and range from 0.05 fb to 3.26 fb, depending on the missing transverse momentum and $b$-quark jet multiplicity requirements.

A search for lepton-flavor-violating $Z\to e\tau$ and $Z\to\mu\tau$ decays with $pp$ collision data recorded by the ATLAS detector at the LHC is presented. This analysis uses 139 fb$^{-1}$ of Run 2 $pp$ collisions at $\sqrt{s}=13$ TeV and is combined with the results of a similar ATLAS search in the final state in which the $\tau$-lepton decays hadronically, using the same data set as well as Run 1 data. The addition of leptonically decaying $\tau$-leptons significantly improves the sensitivity reach for $Z\to\ell\tau$ decays. The $Z\to\ell\tau$ branching fractions are constrained in this analysis to $\mathcal{B}(Z\to e\tau)<7.0\times10^{-6}$ and $\mathcal{B}(Z\to \mu\tau)<7.2\times10^{-6}$ at 95% confidence level. The combination with the previously published analyses sets the strongest constraints to date: $\mathcal{B}(Z\to e\tau)<5.0\times10^{-6}$ and $\mathcal{B}(Z\to \mu\tau)<6.5\times10^{-6}$ at 95% confidence level.