This paper presents studies of Bose-Einstein correlations (BEC) in proton-proton collisions at a centre-of-mass energy of 13 TeV, using data from the ATLAS detector at the CERN Large Hadron Collider. Data were collected in a special low-luminosity configuration with a minimum-bias trigger and a high-multiplicity track trigger, accumulating integrated luminosities of 151 $\mu$b$^{-1}$ and 8.4 nb$^{-1}$ respectively. The BEC are measured for pairs of like-sign charged particles, each with $|\eta|$ < 2.5, for two kinematic ranges: the first with particle $p_T$ > 100 MeV and the second with particle $p_T$ > 500 MeV. The BEC parameters, characterizing the source radius and particle correlation strength, are investigated as functions of charged-particle multiplicity (up to 300) and average transverse momentum of the pair (up to 1.5 GeV). The double-differential dependence on charged-particle multiplicity and average transverse momentum of the pair is also studied. The BEC radius is found to be independent of the charged-particle multiplicity for high charged-particle multiplicity (above 100), confirming a previous observation at lower energy. This saturation occurs independent of the transverse momentum of the pair.

Jet shapes have been measured in inclusive jet production in proton-proton collisions at sqrt(s) = 7 TeV using 3 pb^{-1} of data recorded by the ATLAS experiment at the LHC. Jets are reconstructed using the anti-kt algorithm with transverse momentum 30 GeV < pT < 600 GeV and rapidity in the region |y| < 2.8. The data are corrected for detector effects and compared to several leading-order QCD matrix elements plus parton shower Monte Carlo predictions, including different sets of parameters tuned to model fragmentation processes and underlying event contributions in the final state. The measured jets become narrower with increasing jet transverse momentum and the jet shapes present a moderate jet rapidity dependence. Within QCD, the data test a variety of perturbative and non-perturbative effects. In particular, the data show sensitivity to the details of the parton shower, fragmentation, and underlying event models in the Monte Carlo generators. For an appropriate choice of the parameters used in these models, the data are well described.

A study of $B_c^+\to J/\psi D_s^+$ and $B_c^+\to J/\psi D_s^{*+}$ decays using 139 fb$^{-1}$ of integrated luminosity collected with the ATLAS detector from $\sqrt{s} = 13$ TeV $pp$ collisions at the LHC is presented. The ratios of the branching fractions of the two decays to the branching fraction of the $B_c^+\to J/\psi \pi^+$ decay are measured: $\mathcal B(B_c^+\to J/\psi D_s^+)/\mathcal B(B_c^+\to J/\psi \pi^+) = 2.76\pm 0.47$ and $\mathcal B(B_c^+\to J/\psi D_s^{*+})/\mathcal B(B_c^+\to J/\psi \pi^+) = 5.33\pm 0.96$. The ratio of the branching fractions of the two decays is found to be $\mathcal B(B_c^+\to J/\psi D_s^{*+})/\mathcal B(B_c^+\to J/\psi D_s^+) = 1.93\pm0.26$. For the $B_c^+\to J/\psi D_s^{*+}$ decay, the transverse polarization fraction, $\Gamma_{\pm\pm}/\Gamma$, is measured to be $0.70\pm0.11$. The reported uncertainties include both the statistical and systematic components added in quadrature. The precision of the measurements exceeds that in all previous studies of these decays. These results supersede those obtained in the earlier ATLAS study of the same decays with $\sqrt{s} = 7$ and 8 TeV $pp$ collision data. A comparison with available theoretical predictions for the measured quantities is presented.

A search for nonresonant Higgs boson pair production in the $b\bar{b}\gamma\gamma$ final state is performed using 140 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. This analysis supersedes and expands upon the previous nonresonant ATLAS results in this final state based on the same data sample. The analysis strategy is optimised to probe anomalous values not only of the Higgs ($H$) boson self-coupling modifier $\kappa_\lambda$ but also of the quartic $HHVV$ ($V=W,Z$) coupling modifier $\kappa_{2V}$. No significant excess above the expected background from Standard Model processes is observed. An observed upper limit $\mu_{HH}<4.0$ is set at 95% confidence level on the Higgs boson pair production cross-section normalised to its Standard Model prediction. The 95% confidence intervals for the coupling modifiers are $-1.4<\kappa_\lambda<6.9$ and $-0.5<\kappa_{2V}<2.7$, assuming all other Higgs boson couplings except the one under study are fixed to the Standard Model predictions. The results are interpreted in the Standard Model effective field theory and Higgs effective field theory frameworks in terms of constraints on the couplings of anomalous Higgs boson (self-)interactions.

Searches for new phenomena inspired by supersymmetry in final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair, jets, and missing transverse momentum are presented. These searches make use of proton-proton collision data with an integrated luminosity of 139 $\text{fb}^{-1}$, collected during 2015-2018 at a centre-of-mass energy $\sqrt{s}=13 $TeV by the ATLAS detector at the Large Hadron Collider. Two searches target the pair production of charginos and neutralinos. One uses the recursive-jigsaw reconstruction technique to follow up on excesses observed in 36.1 $\text{fb}^{-1}$ of data, and the other uses conventional event variables. The third search targets pair production of coloured supersymmetric particles (squarks or gluinos) decaying through the next-to-lightest neutralino $(\tilde\chi_2^0)$ via a slepton $(\tilde\ell)$ or $Z$ boson into $\ell^+\ell^-\tilde\chi_1^0$, resulting in a kinematic endpoint or peak in the dilepton invariant mass spectrum. The data are found to be consistent with the Standard Model expectations. Results are interpreted using simplified models and exclude masses up to 900 GeV for electroweakinos, 1550 GeV for squarks, and 2250 GeV for gluinos.

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

A search for a new massive charged gauge boson, $W'$, is performed with the ATLAS detector at the LHC. The dataset used in this analysis was collected from proton-proton collisions at a centre-of-mass energy of $\sqrt{s} =13$ TeV, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The reconstructed $tb$ invariant mass is used to search for a $W'$ boson decaying into a top quark and a bottom quark. The result is interpreted in terms of a $W'$ boson with purely right-handed or left-handed chirality in a mass range of 0.5-6 TeV. Different values for the coupling of the $W'$ boson to the top and bottom quarks are considered, taking into account interference with single-top-quark production in the $s$-channel. No significant deviation from the background prediction is observed. The results are expressed as upper limits on the $W' \rightarrow tb$ production cross-section times branching ratio as a function of the $W'$-boson mass and in the plane of the coupling vs the $W'$-boson mass.

A search for the pair production of heavy leptons as predicted by the type-III seesaw mechanism is presented. The search uses proton-proton collision data at a centre-of-mass energy of 13 TeV, corresponding to 139 fb$^{-1}$ of integrated luminosity recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. The analysis focuses on final states with three or four electrons or muons from the possible decays of new heavy leptons via intermediate electroweak bosons. No significant deviations above the Standard Model expectation are observed; upper and lower limits on the heavy lepton production cross-section and masses are derived respectively. These results are then combined for the first time with the ones already published by ATLAS using the channel with two leptons in the final state. The observed lower limit on the mass of the type-III seesaw heavy leptons combining two, three and four lepton channels together is 910 GeV at the 95% confidence level.

A search for supersymmetry involving the pair production of gluinos decaying via off-shell third-generation squarks into the lightest neutralino ($\tilde\chi^0_1$) is reported. It exploits LHC proton$-$proton collision data at a centre-of-mass energy $\sqrt{s} = 13$ TeV with an integrated luminosity of 139 fb$^{-1}$ collected with the ATLAS detector from 2015 to 2018. The search uses events containing large missing transverse momentum, up to one electron or muon, and several energetic jets, at least three of which must be identified as containing $b$-hadrons. Both a simple kinematic event selection and an event selection based upon a deep neural-network are used. No significant excess above the predicted background is found. In simplified models involving the pair production of gluinos that decay via off-shell top (bottom) squarks, gluino masses less than 2.44 TeV (2.35 TeV) are excluded at 95% CL for a massless $\tilde\chi^0_1$. Limits are also set on the gluino mass in models with variable branching ratios for gluino decays to $b\bar{b}\tilde\chi^0_1$, $t\bar{t}\tilde\chi^0_1$ and $t\bar{b}\tilde\chi^-_1$ / $\bar{t}b\tilde\chi^+_1$.

A search for supersymmetry in events with four or more charged leptons (electrons, muons and $\tau$-leptons) is presented. The analysis uses a data sample corresponding to $139\,\mbox{fb\(^{-1}\)}$ of proton-proton collisions delivered by the Large Hadron Collider at $\sqrt{s}=13$ TeV and recorded by the ATLAS detector. Four-lepton signal regions with up to two hadronically decaying $\tau$-leptons are designed to target several supersymmetric models, while a general five-lepton signal region targets any new physics phenomena leading to a final state with five charged leptons. Data yields are consistent with Standard Model expectations and results are used to set upper limits on contributions from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of general gauge-mediated supersymmetry, excluding higgsino masses up to $540$ GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of $1.6$ TeV, $1.2$ TeV, and $2.5$ TeV are placed on wino, slepton and gluino masses, respectively.