A measurement of jet substructure variables is presented using data collected in 2016 by the ATLAS experiment at the LHC with proton-proton collisions at $\sqrt{s}=13$ TeV. Large-radius jets groomed with the trimming and soft-drop algorithms are studied. Dedicated event selections are used to study jets produced by light quarks or gluons, and hadronically decaying top quarks and $W$ bosons. The variables measured are sensitive to pronged substructure, and therefore are typically used for tagging jets from boosted massive particles. These include the energy correlation functions and the $N$-subjettiness variables. The number of subjets and the Les Houches angularity are also considered. The distributions of the substructure variables, corrected for detector effects, are compared to the predictions of various Monte Carlo event generators. They are also compared between the large-radius jets originating from light quarks or gluons, and hadronically decaying top quarks and $W$ bosons.

A search is presented for four-top-quark production using an integrated luminosity of 139 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the LHC. Events are selected if they contain a same-sign lepton pair or at least three leptons (electrons or muons). Jet multiplicity, jet flavour and event kinematics are used to separate signal from the background through a multivariate discriminant, and dedicated control regions are used to constrain the dominant backgrounds. The four-top-quark production cross section is measured to be 24$^{+7}_{-6}$ fb. This corresponds to an observed (expected) significance with respect to the background-only hypothesis of 4.3 (2.4) standard deviations and provides evidence for this process.

A generic search for resonances is performed with events containing a $Z$ boson with transverse momentum greater than 100 GeV, decaying into $e^+e^-$ or $\mu^+\mu^-$. The analysed data collected with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider correspond to an integrated luminosity of 139 fb$^{-1}$. Two invariant mass distributions are examined for a localised excess relative to the expected Standard Model background in six independent event categories (and their inclusive sum) to increase the sensitivity. No significant excess is observed. Exclusion limits at 95% confidence level are derived for two cases: a model-independent interpretation of Gaussian-shaped resonances with the mass width between 3% and 10% of the resonance mass, and a specific heavy vector triplet model with the decay mode $W'\to ZW \to \ell\ell qq$.

The elliptic ($v_2$) and triangular ($v_3$) azimuthal anisotropy coefficients in central $^{3}$He+Au, $d$+Au, and $p$+Au collisions at $\mbox{$\sqrt{s_{\mathrm{NN}}}$}$ = 200 GeV are measured as a function of transverse momentum ($p_{\mathrm{T}}$) at mid-rapidity ($|\eta|<$0.9), via the azimuthal angular correlation between two particles both at $|\eta|<$0.9. While the $v_2(p_{\mathrm{T}})$ values depend on the colliding systems, the $v_3(p_{\mathrm{T}})$ values are system-independent within the uncertainties, suggesting an influence on eccentricity from sub-nucleonic fluctuations in these small-sized systems. These results also provide stringent constraints for the hydrodynamic modeling of these systems.

A search for chargino$-$neutralino pair production in three-lepton final states with missing transverse momentum is presented. The study is based on a dataset of $\sqrt{s} = 13$ TeV $pp$ collisions recorded with the ATLAS detector at the LHC, corresponding to an integrated luminosity of 139 fb$^{-1}$. No significant excess relative to the Standard Model predictions is found in data. The results are interpreted in simplified models of supersymmetry, and statistically combined with results from a previous ATLAS search for compressed spectra in two-lepton final states. Various scenarios for the production and decay of charginos ($\tilde\chi^\pm_1$) and neutralinos ($\tilde\chi^0_2$) are considered. For pure higgsino $\tilde\chi^\pm_1\tilde\chi^0_2$ pair-production scenarios, exclusion limits at 95% confidence level are set on $\tilde\chi^0_2$ masses up to 210 GeV. Limits are also set for pure wino $\tilde\chi^\pm_1\tilde\chi^0_2$ pair production, on $\tilde\chi^0_2$ masses up to 640 GeV for decays via on-shell $W$ and $Z$ bosons, up to 300 GeV for decays via off-shell $W$ and $Z$ bosons, and up to 190 GeV for decays via $W$ and Standard Model Higgs bosons.

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.

This paper presents a search for pair production of higgsinos, the supersymmetric partners of the Higgs bosons, in scenarios with gauge-mediated supersymmetry breaking. Each higgsino is assumed to decay into a Higgs boson and a nearly massless gravitino. The search targets events where each Higgs boson decays into $b\bar{b}$, leading to a reconstructed final state with at least three energetic $b$-jets and missing transverse momentum. Two complementary analysis channels are used, with each channel specifically targeting either low or high values of the higgsino mass. The low-mass (high-mass) channel exploits 126 (139) fb$^{-1}$ of $\sqrt{s}=13$ TeV data collected by the ATLAS detector during Run 2 of the Large Hadron Collider. No significant excess above the Standard Model prediction is found. At 95% confidence level, masses between 130 GeV and 940 GeV are excluded for higgsinos decaying exclusively into Higgs bosons and gravitinos. Exclusion limits as a function of the higgsino decay branching ratio to a Higgs boson are also reported.

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.

The production of a $W$ boson in association with a single charm quark is studied using 140 $\mathrm{fb}^{-1}$ of $\sqrt{s} = 13\,\mathrm{TeV}$ proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider. The charm quark is tagged by a charmed hadron, reconstructed with a secondary-vertex fit. The $W$ boson is reconstructed from an electron/muon decay and the missing transverse momentum. The mesons reconstructed are $D^{\pm} \to K^\mp \pi^\pm \pi^\pm$ and $D^{*\pm} \to D^{0} \pi^\pm \to (K^\mp \pi^\pm) \pi^\pm$, where $p_{\text{T}}(e, \mu) > 30\,\mathrm{GeV}$, $|\eta(e, \mu)| < 2.5$, $p_{\text{T}}(D) > 8\,\mathrm{GeV}$, and $|\eta(D)| < 2.2$. The integrated and normalized differential cross-sections as a function of the pseudorapidity of the lepton from the $W$ boson decay, and of the transverse momentum of the meson, are extracted from the data using a profile likelihood fit. The measured fiducial cross-sections are $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{-}{+}D^{+}) = 50.2\pm0.2\,\mathrm{(stat.)}\,^{+2.4}_{-2.3}\,\mathrm{(syst.)}\,\mathrm{pb}$, $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{+}{+}D^{-}) = 48.5\pm0.2\,\mathrm{(stat.)}\,^{+2.3}_{-2.2}\,\mathrm{(syst.)}\,\mathrm{pb}$, $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{-}{+}D^{*+}) = 51.1\pm0.4\,\mathrm{(stat.)}\,^{+1.9}_{-1.8}\,\mathrm{(syst.)}\,\mathrm{pb}$, and $\sigma^{\mathrm{OS-SS}}_{\mathrm{fid}}(W^{+}{+}D^{*-}) = 50.0\pm0.4\,\mathrm{(stat.)}\,^{+1.9}_{-1.8}\,\mathrm{(syst.)}\,\mathrm{pb}$. Results are compared with the predictions of next-to-leading-order quantum chromodynamics calculations performed using state-of-the-art parton distribution functions. The ratio of charm to anti-charm production cross-sections is studied to probe the $s$-$\bar{s}$ quark asymmetry and is found to be $R_c^\pm = 0.971\pm0.006\,\mathrm{(stat.)}\pm0.011\,\mathrm{(syst.)}$.

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