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 search is presented for the pair production of heavy vector-like $B$ quarks, primarily targeting $B$ quark decays into a $W$ boson and a top quark. The search is based on $36.1$ $fb^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 13 TeV recorded in 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider. Data are analysed in the lepton-plus-jets final state, characterised by a high-transverse-momentum isolated electron or muon, large missing transverse momentum, and multiple jets, of which at least one is $b$-tagged. No significant deviation from the Standard Model expectation is observed. The 95% confidence level lower limit on the $B$ mass is 1350 GeV assuming a 100% branching ratio to $Wt$. In the SU(2) singlet scenario, the lower mass limit is 1170 GeV. This search is also sensitive to a heavy vector-like $B$ quark decaying into other final states ($Zb$ and $Hb$) and thus mass limits on $B$ production are set as a function of the decay branching ratios. The 100% branching ratio limits are found to be also applicable to heavy vector-like $X$ production, with charge $+$5/3, that decay into $Wt$.

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 a heavy resonance decaying into $WZ$ in the fully leptonic channel (electrons and muons) is performed. It is based on proton-proton collision data collected by the ATLAS experiment at the Large Hadron Collider at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 36.1 fb$^{-1}$. No significant excess is observed over the Standard Model predictions and limits are set on the production cross section times branching ratio of a heavy vector particle produced either in quark-antiquark fusion or through vector-boson fusion. Constraints are also obtained on the mass and couplings of a singly charged Higgs boson, in the Georgi-Machacek model, produced through vector-boson fusion.

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.