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.

A search for exotic Higgs boson decays to light pseudoscalars in the final state of two muons and two $\tau$ leptons is performed using proton-proton collision data recorded by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Masses of the pseudoscalar boson between 15.0 and 62.5 GeV are probed, and no significant excess of data is observed above the prediction of the standard model. Upper limits are set on the branching fraction of the Higgs boson to two light pseudoscalar bosons in different types of two-Higgs-doublet models extended with a complex scalar singlet.

A search for an exotic decay of the Higgs boson to a pair of light pseudoscalar bosons is performed for the first time in the final state with two b quarks and two $\tau$ leptons. The search is motivated in the context of models of physics beyond the standard model (SM), such as two Higgs doublet models extended with a complex scalar singlet (2HDM+S), which include the next-to-minimal supersymmetric SM (NMSSM). The results are based on a data set of proton-proton collisions corresponding to an integrated luminosity of 35.9 fb$^{-1}$, accumulated by the CMS experiment at the LHC in 2016 at a center-of-mass energy of 13 TeV. Masses of the pseudoscalar boson between 15 and 60 GeV are probed, and no excess of events above the SM expectation is observed. Upper limits between 3 and 12% are set on the branching fraction $\mathcal{B}$(h $\to$ aa $\to$ 2$\tau$2b) assuming the SM production of the Higgs boson. Upper limits are also set on the branching fraction of the Higgs boson to two light pseudoscalar bosons in different 2HDM+S scenarios. Assuming the SM production cross section for the Higgs boson, the upper limit on this quantity is as low as 20% for a mass of the pseudoscalar of 40 GeV in the NMSSM.

A search for heavy neutral Higgs bosons and $Z^{\prime}$ bosons is performed using a data sample corresponding to an integrated luminosity of 36.1 fb$^{-1}$ from proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS detector at the LHC during 2015 and 2016. The heavy resonance is assumed to decay to $\tau^+\tau^-$ with at least one tau lepton decaying to final states with hadrons and a neutrino. The search is performed in the mass range of 0.2-2.25 TeV for Higgs bosons and 0.2-4.0 TeV for $Z^{\prime}$ bosons. The data are in good agreement with the background predicted by the Standard Model. The results are interpreted in benchmark scenarios. In the context of the hMSSM scenario, the data exclude $\tan\beta > 1.0$ for $m_A$ = 0.25 TeV and $\tan\beta > 42$ for $m_A$ = 1.5 TeV at the 95% confidence level. For the Sequential Standard Model, $Z^{\prime}_\mathrm{SSM}$ with $m_{Z^{\prime}} < 2.42$ TeV is excluded at 95% confidence level, while $Z^{\prime}_\mathrm{NU}$ with $m_{Z^{\prime}} < 2.25$ TeV is excluded for the non-universal $G(221)$ model that exhibits enhanced couplings to third-generation fermions.

A search for heavy resonances decaying into a pair of $Z$ bosons leading to $\ell^+\ell^-\ell^+\ell^-$ and $\ell^+\ell^-\nu\bar\nu$ final states, where $\ell$ stands for either an electron or a muon, is presented. The search uses proton proton collision data at a centre-of-mass energy of 13 TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ collected with the ATLAS detector during 2015 and 2016 at the Large Hadron Collider. Different mass ranges for the hypothetical resonances are considered, depending on the final state and model. The different ranges span between 200 GeV and 2000 GeV. The results are interpreted as upper limits on the production cross section of a spin 0 or spin 2 resonance. The upper limits for the spin 0 resonance are translated to exclusion contours in the context of Type I and Type II two-Higgs-doublet models, while those for the spin 2 resonance are used to constrain the Randall Sundrum model with an extra dimension giving rise to spin 2 graviton excitations.

The associated production of a Higgs boson with a $W$ or $Z$ boson decaying into leptons and where the Higgs boson decays to a $b\bar{b}$ pair is measured in the high vector-boson transverse momentum regime, above 250 GeV, with the ATLAS detector. The analysed data, corresponding to an integrated luminosity of 139 fb$^{-1}$, were collected in proton-proton collisions at the Large Hadron Collider between 2015 and 2018 at a centre-of-mass energy of $\sqrt{s} = 13$ TeV. The measured signal strength, defined as the ratio of the measured signal yield to that predicted by the Standard Model, is $0.72 ^{+0.39}_{-0.36}$ corresponding to an observed (expected) significance of 2.1 (2.7) standard deviations. Cross-sections of associated production of a Higgs boson decaying into $b$ quark pairs with a $W$ or $Z$ gauge boson, decaying into leptons, are measured in two exclusive vector boson transverse momentum regions, 250-400 GeV and above 400 GeV, and interpreted as constraints on anomalous couplings in the framework of a Standard Model effective field theory.

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.

A search for heavy resonances decaying into a pair of $Z$ bosons leading to $\ell^+\ell^-\ell'^+\ell'^-$ and $\ell^+\ell^-\nu\bar\nu$ final states, where $\ell$ stands for either an electron or a muon, is presented. The search uses proton-proton collision data at a centre-of-mass energy of 13 TeV collected from 2015 to 2018 that corresponds to the full integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. Different mass ranges spanning 200 GeV to 2000 GeV for the hypothetical resonances are considered, depending on the final state and model. In the absence of a significant observed excess, the results are interpreted as upper limits on the production cross section of a spin-0 or spin-2 resonance. The upper limits for the spin-0 resonance are translated to exclusion contours in the context of Type-I and Type-II two-Higgs-doublet models, and the limits for the spin-2 resonance are used to constrain the Randall--Sundrum model with an extra dimension giving rise to spin-2 graviton excitations.