This paper reports a search for Higgs boson pair ($hh$) production in association with a vector boson ($W$ or $Z$) using 139 $fb^{-1}$ of proton-proton collision data at $\sqrt{s}=$ 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The search is performed in final states in which the vector boson decays leptonically ($W\to\ell\nu, Z\to\ell\ell,\nu\nu$ with $\ell=e, \mu$) and the Higgs bosons each decay into a pair of $b$-quarks. It targets $Vhh$ signals from both non-resonant $hh$ production, present in the Standard Model (SM), and resonant $hh$ production, as predicted in some SM extensions. A 95% confidence-level upper limit of 183 (87) times the SM cross-section is observed (expected) for non-resonant $Vhh$ production when assuming the kinematics are as expected in the SM. Constraints are also placed on Higgs boson coupling modifiers. For the resonant search, upper limits on the production cross-sections are derived for two specific models: one is the production of a vector boson along with a neutral heavy scalar resonance $H$, in the mass range 260-1000 GeV, that decays into $hh$, and the other is the production of a heavier neutral pseudoscalar resonance $A$ that decays into a $Z$ boson and $H$ boson, where the $A$ boson mass is 360-800 GeV and the $H$ boson mass is 260-400 GeV. Constraints are also derived in the parameter space of two-Higgs-doublet models.

A search for decays of pair-produced neutral long-lived particles (LLPs) is presented using 139 fb$^{-1}$ of proton-proton collision data collected by the ATLAS detector at the LHC in 2015-2018 at a centre-of-mass energy of 13 TeV. Dedicated techniques were developed for the reconstruction of displaced jets produced by LLPs decaying hadronically in the ATLAS hadronic calorimeter. Two search regions are defined for different LLP kinematic regimes. The observed numbers of events are consistent with the expected background, and limits for several benchmark signals are determined. For a SM Higgs boson with a mass of 125 GeV, branching ratios above 10% are excluded at 95% confidence level for values of $c$ times LLP mean proper lifetime in the range between 20 mm and 10 m depending on the model. Upper limits are also set on the cross-section times branching ratio for scalars with a mass of 60 GeV and for masses between 200 GeV and 1 TeV.

Constraints on the Higgs boson self-coupling are set by combining double-Higgs boson analyses in the $b\bar{b}b\bar{b}$, $b\bar{b}\tau^+\tau^-$ and $b\bar{b} \gamma \gamma$ decay channels with single-Higgs boson analyses targeting the $\gamma \gamma$, $ZZ^*$, $WW^*$, $\tau^+ \tau^-$ and $b\bar{b}$ decay channels. The data used in these analyses were recorded by the ATLAS detector at the LHC in proton$-$proton collisions at $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 126$-$139 fb$^{-1}$. The combination of the double-Higgs analyses sets an upper limit of $\mu_{HH} < 2.4$ at 95% confidence level on the double-Higgs production cross-section normalised to its Standard Model prediction. Combining the single-Higgs and double-Higgs analyses, with the assumption that new physics affects only the Higgs boson self-coupling ($\lambda_{HHH}$), values outside the interval $-0.4< \kappa_{\lambda}=(\lambda_{HHH}/\lambda_{HHH}^{\textrm{SM}})< 6.3$ are excluded at 95% confidence level. The combined single-Higgs and double-Higgs analyses provide results with fewer assumptions, by adding in the fit more coupling modifiers introduced to account for the Higgs boson interactions with the other Standard Model particles. In this relaxed scenario, the constraint becomes $-1.4 < \kappa_{\lambda} < 6.1$ at 95% CL.

Heavy-flavour hadron production provides information about the transport properties and microscopic structure of the quark-gluon plasma created in ultra-relativistic heavy-ion collisions. A measurement of the muons from semileptonic decays of charm and bottom hadrons produced in Pb+Pb and $pp$ collisions at a nucleon-nucleon centre-of-mass energy of 5.02 TeV with the ATLAS detector at the Large Hadron Collider is presented. The Pb+Pb data were collected in 2015 and 2018 with sampled integrated luminosities of $208~\mathrm{\mu b}^{-1}$ and $38~\mathrm{\mu b^{-1}}$, respectively, and $pp$ data with a sampled integrated luminosity of $1.17~\mathrm{pb}^{-1}$ were collected in 2017. Muons from heavy-flavour semileptonic decays are separated from the light-flavour hadronic background using the momentum imbalance between the inner detector and muon spectrometer measurements, and muons originating from charm and bottom decays are further separated via the muon track's transverse impact parameter. Differential yields in Pb+Pb collisions and differential cross sections in $pp$ collisions for such muons are measured as a function of muon transverse momentum from 4 GeV to 30 GeV in the absolute pseudorapidity interval $|\eta| < 2$. Nuclear modification factors for charm and bottom muons are presented as a function of muon transverse momentum in intervals of Pb+Pb collision centrality. The measured nuclear modification factors quantify a significant suppression of the yields of muons from decays of charm and bottom hadrons, with stronger effects for muons from charm hadron decays.

A search for quantum black holes in electron+jet and muon+jet invariant mass spectra is performed with 140 fb$^{-1}$ of data collected by the ATLAS detector in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider. The observed invariant mass spectrum of lepton+jet pairs is consistent with Standard Model expectations. Upper limits are set at 95% confidence level on the production cross-sections times branching fractions for quantum black holes decaying into a lepton and a quark in a search region with invariant mass above 2.0 TeV. The resulting quantum black hole lower mass threshold limit is 9.2 TeV in the Arkani-Hamed-Dimopoulos-Dvali model, and 6.8 TeV in the Randall-Sundrum model.

A search for Higgs boson pair production via vector-boson fusion (VBF) in the $b\bar{b}b\bar{b}$ final state is carried out with the ATLAS experiment using 126 fb$^{-1}$ of proton-proton collision data delivered at $\sqrt{s} = 13$ TeV by the Large Hadron Collider. This search is sensitive to VBF production of additional heavy bosons that may decay into Higgs boson pairs, and in a non-resonant topology it can constrain the quartic coupling between the Higgs bosons and vector bosons. No significant excess relative to the Standard Model expectation is observed, and limits on the production cross-section are set at the 95 % confidence level for a heavy scalar resonance in the context of an extended Higgs sector, and for non-resonant Higgs boson pair production. Interpretation in terms of the coupling between a Higgs boson pair and two vector bosons is also provided: coupling values normalised to the Standard Model expectation of $\kappa_{2V} < -0.43$ and $\kappa_{2V} > 2.56$ are excluded at the 95 % confidence level in data.

A search for charged Higgs bosons decaying into $W^\pm W^\pm$ or $W^\pm Z$ bosons is performed, involving experimental signatures with two leptons of the same charge, or three or four leptons with a variety of charge combinations, missing transverse momentum and jets. A data sample of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider between 2015 and 2018 is used. The data correspond to a total integrated luminosity of 139 fb$^{-1}$. The search is guided by a type-II seesaw model that extends the scalar sector of the Standard Model with a scalar triplet, leading to a phenomenology that includes doubly and singly charged Higgs bosons. Two scenarios are explored, corresponding to the pair production of doubly charged $H^{\pm\pm}$ bosons, or the associated production of a doubly charged $H^{\pm\pm}$ boson and a singly charged $H^\pm$ boson. No significant deviations from the Standard Model predictions are observed. $H^{\pm\pm}$ bosons are excluded at 95% confidence level up to 350 GeV and 230 GeV for the pair and associated production modes, respectively.

A search for pair production of scalar leptoquarks, each decaying into either an electron or a muon and a top quark, is presented. This is the first leptoquark search using ATLAS data to investigate top-philic cross-generational couplings that could provide explanations for recently observed anomalies in $B$ meson decays. This analysis targets high leptoquark masses which cause the decay products of each resultant top quark to be contained within a single high-$p_{\mathrm{T}}$ large-radius jet. The full Run 2 dataset is exploited, consisting of 139 fb$^{-1}$ of data collected from proton-proton collisions at $\sqrt{s}=13$ TeV from 2015 to 2018 with the ATLAS detector at the CERN Large Hadron Collider. In the absence of any significant deviation from the background expectation, lower limits on the leptoquark masses are set at 1480 GeV and 1470 GeV for the electron and muon channel, respectively.

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.

Studies of the fragmentation of jets into charged particles in heavy-ion collisions can provide information about the mechanism of jet-quenching by the hot and dense QCD matter created in such collisions, the quark-gluon plasma. This paper presents a measurement of the angular distribution of charged particles around the jet axis in $\sqrt{s_{\mathrm{NN}}}=$ 5.02 TeV Pb+Pb and $pp$ collisions, using the ATLAS detector at the LHC. The Pb+Pb and $pp$ data sets have integrated luminosities of 0.49 nb$^{-1}$ and 25 pb$^{-1}$, respectively. The measurement is performed for jets reconstructed with the anti-$k_{t}$ algorithm with radius parameter $R = 0.4$ and is extended to an angular distance of $r= 0.8$ from the jet axis. Results are presented as a function of Pb+Pb collision centrality and distance from the jet axis for charged particles with transverse momenta in the 1$-$63 GeV range, matched to jets with transverse momenta in the 126$-$316 GeV range and an absolute value of jet rapidity of less than 1.7. Modifications to the measured distributions are quantified by taking a ratio to the measurements in $pp$ collisions. Yields of charged particles with transverse momenta below 4 GeV are observed to be increasingly enhanced as a function of angular distance from the jet axis, reaching a maximum at $r=0.6$. Charged particles with transverse momenta above 4 GeV have an enhanced yield in Pb+Pb collisions in the jet core for angular distances up to $r = 0.05$ from the jet axis, with a suppression at larger distances.