Searches for the exclusive decays of Higgs and $Z$ bosons into a vector quarkonium state and a photon are performed in the $\mu^+\mu^- \gamma$ final state with a proton$-$proton collision data sample corresponding to an integrated luminosity of $139$ fb$^{-1}$ collected at $\sqrt{s}=13$ TeV with the ATLAS detector at the CERN Large Hadron Collider. The observed data are compatible with the expected backgrounds. The 95% confidence-level upper limits on the branching fractions of the Higgs boson decays into $J/\psi \gamma$, $\psi(2S) \gamma$, and $\Upsilon(1S,2S,3S) \gamma$ are found to be $2.0\times10^{-4}$, $10.5\times10^{-4}$, and $(2.5,4.2,3.4)\times10^{-4}$, respectively, assuming Standard Model production of the Higgs boson. The corresponding 95% CL upper limits on the branching fractions of the $Z$ boson decays are $1.2\times10^{-6}$, $2.4\times10^{-6}$, and $(1.1,1.3,2.4)\times10^{-6}$. An observed 95% CL interval of $(-133,175)$ is obtained for the $\kappa_c/\kappa_\gamma$ ratio of Higgs boson coupling modifiers, and a 95% CL interval of $(-37,40)$ is obtained for $\kappa_b/\kappa_\gamma$.

This paper presents measurements of charged-hadron spectra obtained in $pp$, $p$+Pb, and Pb+Pb collisions at $\sqrt{s}$ or $\sqrt{s_{_\text{NN}}}=5.02$ TeV, and in Xe+Xe collisions at $\sqrt{s_{_\text{NN}}}=5.44$ TeV. The data recorded by the ATLAS detector at the LHC have total integrated luminosities of 25 pb${}^{-1}$, 28 nb${}^{-1}$, 0.50 nb${}^{-1}$, and 3 $\mu$b${}^{-1}$, respectively. The nuclear modification factors $R_{p\text{Pb}}$ and $R_\text{AA}$ are obtained by comparing the spectra in heavy-ion and $pp$ collisions in a wide range of charged-particle transverse momenta and pseudorapidity. The nuclear modification factor $R_{p\text{Pb}}$ shows a moderate enhancement above unity with a maximum at $p_{\mathrm{T}} \approx 3$ GeV; the enhancement is stronger in the Pb-going direction. The nuclear modification factors in both Pb+Pb and Xe+Xe collisions feature a significant, centrality-dependent suppression. They show a similar distinct $p_{\mathrm{T}}$-dependence with a local maximum at $p_{\mathrm{T}} \approx 2$ GeV and a local minimum at $p_{\mathrm{T}} \approx 7$ GeV. This dependence is more distinguishable in more central collisions. No significant $|\eta|$-dependence is found. A comprehensive comparison with several theoretical predictions is also provided. They typically describe $R_\text{AA}$ better in central collisions and in the $p_{\mathrm{T}}$ range from about 10 to 100 GeV.

The azimuthal variation of jet yields in heavy-ion collisions provides information about the path-length dependence of the energy loss experienced by partons passing through the hot, dense nuclear matter known as the quark-gluon plasma. This paper presents the azimuthal anisotropy coefficients $v_2$, $v_3$, and $v_4$ measured for jets in Pb+Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV using the ATLAS detector at the LHC. The measurement uses data collected in 2015 and 2018, corresponding to an integrated luminosity of 2.2 nb$^{-1}$. The $v_n$ values are measured as a function of the transverse momentum of the jets between 71 GeV and 398 GeV and the event centrality. A nonzero value of $v_2$ is observed in all but the most central collisions. The value of $v_2$ is largest for jets with lower transverse momentum, with values up to 0.05 in mid-central collisions. A smaller, nonzero value of $v_3$ of approximately 0.01 is measured with no significant dependence on jet $p_T$ or centrality, suggesting that fluctuations in the initial state play a small but distinct role in jet energy loss. No significant deviation of $v_4$ from zero is observed in the measured kinematic region.

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.

The ATLAS detector at the Large Hadron Collider has been used to measure jet substructure modification and suppression in Pb+Pb collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_\mathrm{NN}}}=5.02~\mathrm{TeV}$ in comparison with $pp$ collisions at $\sqrt{s}=5.02~\mathrm{TeV}$. The Pb+Pb data, collected in 2018, have an integrated luminosity of $1.72~\mathrm{nb^{-1}}$, while the $pp$ data, collected in 2017, have an integrated luminosity of $260~\mathrm{pb}^{-1}$. Jets used in this analysis are clustered using the anti-$k_{t}$ algorithm with a radius parameter $R=0.4$. The jet constituents, defined by both tracking and calorimeter information, are used to determine the angular scale $r_\mathrm{g}$ of the first hard splitting inside the jet by reclustering them using the Cambridge-Aachen algorithm and employing the soft-drop grooming technique. The nuclear modification factor, $R_\mathrm{AA}$, used to characterize jet suppression in Pb+Pb collisions, is presented differentially in $r_\mathrm{g}$, jet transverse momentum, and in intervals of collision centrality. The $R_\mathrm{AA}$ value is observed to depend significantly on jet $r_\mathrm{g}$. Jets produced with the largest measured $r_\mathrm{g}$ are found to be twice as suppressed as those with the smallest $r_\mathrm{g}$ in central Pb+Pb collisions. The $R_\mathrm{AA}$ values do not exhibit a strong variation with jet $p_\mathrm{T}$ in any of the $r_\mathrm{g}$ intervals. The $r_\mathrm{g}$ and $p_\mathrm{T}$ dependence of jet $R_\mathrm{AA}$ is qualitatively consistent with a picture of jet quenching arising from coherence and provides the most direct evidence in support of this approach.

A measurement of four-top-quark production using proton-proton collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the Large Hadron Collider corresponding to an integrated luminosity of 139 fb$^{-1}$ is presented. Events are selected if they contain a single lepton (electron or muon) or an opposite-sign lepton pair, in association with multiple jets. The events are categorised according to the number of jets and how likely these are to contain $b$-hadrons. A multivariate technique is then used to discriminate between signal and background events. The measured four-top-quark production cross section is found to be 26$^{+17}_{-15}$ fb, with a corresponding observed (expected) significance of 1.9 (1.0) standard deviations over the background-only hypothesis. The result is combined with the previous measurement performed by the ATLAS Collaboration in the multilepton final state. The combined four-top-quark production cross section is measured to be 24$^{+7}_{-6}$ fb, with a corresponding observed (expected) signal significance of 4.7 (2.6) standard deviations over the background-only predictions. It is consistent within 2.0 standard deviations with the Standard Model expectation of 12.0$\pm$2.4 fb.

A search for charged-lepton-flavour violating $\mu\tau qt$ ($q=u,c$) interactions is presented, considering both top-quark production and decay. The data analysed correspond to 140 $\textrm{fb}^{-1}$ of proton-proton collisions at a centre-of-mass energy of $\sqrt{s}= $13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The analysis targets events containing two muons with the same electric charge, a hadronically decaying $\tau$-lepton and at least one jet, with exactly one $b$-tagged jet, produced by a $\mu\tau qt$ interaction. Agreement with the Standard Model expectation within $1.6\sigma$ is observed, and limits are set at the 95% CL on the charged-lepton-flavour violation branching ratio of $\mathcal{B}(t \to \mu\tau q) < 8.7 \times 10^{-7}$. An Effective Field Theory interpretation is performed yielding 95% CL limits on Wilson coefficients, dependent on the flavour of the associated light quark and the Lorentz structure of the coupling. These range from $|c_{\mathsf{lequ}}^{3(2313)}| / \Lambda^{2} < 0.10\textrm{ TeV}^{-2}$ for $\mu\tau ut$ to $|c_{\mathsf{ lequ}}^{1(2323)}| / \Lambda^{2} < 1.8\textrm{ TeV}^{-2}$ for $\mu\tau ct$. An additional interpretation is performed for scalar leptoquark production inducing charged lepton flavour violation, with fixed inter-generational couplings. Upper limits on leptoquark coupling strengths are set at the 95% CL, ranging from $\lambda^{\textrm{LQ}} = $1.3 to $\lambda^{\textrm{LQ}} = $3.7 for leptoquark masses between 0.5 and 2.0 TeV.

A search is presented for flavour-changing neutral-current interactions involving the top quark, the Higgs boson and an up-type quark ($q=u,c$) with the ATLAS detector at the Large Hadron Collider. The analysis considers leptonic decays of the top quark along with Higgs boson decays into two $W$ bosons, two $Z$ bosons or a $\tau^{+}\tau^{-}$ pair. It focuses on final states containing either two leptons (electrons or muons) of the same charge or three leptons. The considered processes are $t\bar{t}$ and $Ht$ production. For the $t\bar{t}$ production, one top quark decays via $t\to Hq$. The proton-proton collision data set analysed amounts to 140 fb$^{-1}$ at $\sqrt{s}=13$ TeV. No significant excess beyond Standard Model expectations is observed and upper limits are set on the $t\to Hq$ branching ratios at 95% confidence level, amounting to observed (expected) limits of $\mathcal{B}(t\to Hu)<2.8\,(3.0) \times 10^{-4}$ and $\mathcal{B}(t\to Hc)<3.3\,(3.8) \times 10^{-4}$. Combining this search with other searches for $tHq$ flavour-changing neutral-current interactions previously conducted by ATLAS, considering $H\to b\bar{b}$ and $H\to\gamma\gamma$ decays, as well as $H\to\tau^{+}\tau^{-}$ decays with one or two hadronically decaying $\tau$-leptons, yields observed (expected) upper limits on the branching ratios of $\mathcal{B}(t\to Hu)<2.6\,(1.8) \times 10^{-4}$ and $\mathcal{B}(t\to Hc)<3.4\,(2.3) \times 10^{-4}$.

The electroweak production of $Z(\nu\bar{\nu})\gamma$ in association with two jets is studied in a regime with a photon of high transverse momentum above 150 GeV using proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider. The analysis uses a data sample with an integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector during the 2015-2018 LHC data-taking period. This process is an important probe of the electroweak symmetry breaking mechanism in the Standard Model and is sensitive to quartic gauge boson couplings via vector-boson scattering. The fiducial $Z(\nu\bar{\nu})\gamma jj$ cross section for electroweak production is measured to be 0.77$^{+0.34}_{-0.30}$ fb and is consistent with the Standard Model prediction. Evidence of electroweak $Z(\nu\bar{\nu})\gamma jj$ production is found with an observed significance of 3.2$\sigma$ for the background-only hypothesis, compared with an expected significance of 3.7$\sigma$. The combination of this result with the previously published ATLAS observation of electroweak $Z(\nu\bar{\nu})\gamma jj$ production yields an observed (expected) signal significance of 6.3$\sigma$ (6.6$\sigma$). Limits on anomalous quartic gauge boson couplings are obtained in the framework of effective field theory with dimension-8 operators.

A search for pair-produced scalar or vector leptoquarks decaying into a $b$-quark and a $\tau$-lepton is presented using the full LHC Run 2 (2015-2018) data sample of 139 fb$^{-1}$ collected with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} =13$ TeV. Events in which at least one $\tau$-lepton decays hadronically are considered, and multivariate discriminants are used to extract the signals. No significant deviations from the Standard Model expectation are observed and 95% confidence-level upper limits on the production cross-section are derived as a function of leptoquark mass and branching ratio $B$ into a $\tau$-lepton and $b$-quark. For scalar leptoquarks, masses below 1460 GeV are excluded assuming $B=100$%, while for vector leptoquarks the corresponding limit is 1650 GeV (1910 GeV) in the minimal-coupling (Yang-Mills) scenario.