The ATLAS Collaboration has measured the inclusive production of $Z$ bosons via their decays into electron and muon pairs in $p+$Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV at the Large Hadron Collider. The measurements are made using data corresponding to integrated luminosities of 29.4 nb$^{-1}$ and 28.1 nb$^{-1}$ for $Z \rightarrow ee$ and $Z \rightarrow \mu\mu$, respectively. The results from the two channels are consistent and combined to obtain a cross section times the $Z \rightarrow \ell\ell$ branching ratio, integrated over the rapidity region $|y^{*}_{Z}|<3.5$, of 139.8 $\pm$ 4.8 (stat.) $\pm$ 6.2 (syst.) $\pm$ 3.8 (lumi.) nb. Differential cross sections are presented as functions of the $Z$ boson rapidity and transverse momentum, and compared with models based on parton distributions both with and without nuclear corrections. The centrality dependence of $Z$ boson production in $p+$Pb collisions is measured and analyzed within the framework of a standard Glauber model and the model's extension for fluctuations of the underlying nucleon-nucleon scattering cross section.

Measurements of $ZZ$ production in the $\ell^{+}\ell^{-}\ell^{\prime +}\ell^{\prime -}$ channel in proton-proton collisions at 13 TeV center-of-mass energy at the Large Hadron Collider are presented. The data correspond to 36.1 $\mathrm{fb}^{-1}$ of collisions collected by the ATLAS experiment in 2015 and 2016. Here $\ell$ and $\ell'$ stand for electrons or muons. Integrated and differential $ZZ \to \ell^{+}\ell^{-}\ell^{\prime +}\ell^{\prime -}$ cross sections with $Z \to \ell^+\ell^-$ candidate masses in the range of 66 GeV to 116 GeV are measured in a fiducial phase space corresponding to the detector acceptance and corrected for detector effects. The differential cross sections are presented in bins of twenty observables, including several that describe the jet activity. The integrated cross section is also extrapolated to a total phase space and to all Standard-Model decays of $Z$ bosons with mass between 66 GeV and 116 GeV, resulting in a value of $17.3 \pm 0.9$ [$\pm 0.6$ (stat.) $\pm 0.5$ (syst.) $\pm 0.6$ (lumi.)] pb. The measurements are found to be in good agreement with the Standard-Model predictions. A search for neutral triple gauge couplings is performed using the transverse momentum distribution of the leading $Z$-boson candidate. No evidence for such couplings is found and exclusion limits are set on their parameters.

The Standard Model of particle physics describes the known fundamental particles and forces that make up our universe, with the exception of gravity. One of the central features of the Standard Model is a field that permeates all of space and interacts with fundamental particles. The quantum excitation of this field, known as Higgs field, manifests itself as the Higgs boson, the only fundamental particle with no spin. In 2012, a particle with properties consistent with the Higgs boson of the Standard Model was observed by the ATLAS and CMS experiments at the Large Hadron Collider at CERN. Since then, more than 30 times as many Higgs bosons have been recorded by the ATLAS experiment, allowing much more precise measurements and new tests of the theory. Here, on the basis of this larger dataset, we combine an unprecedented number of production and decay processes of the Higgs boson to scrutinize its interactions with elementary particles. Interactions with gluons, photons, and $W$ and $Z$ bosons -- the carriers of the strong, electromagnetic, and weak forces -- are studied in detail. Interactions with three third-generation matter particles (bottom ($b$) and top ($t$) quarks, and tau leptons ($\tau$)) are well measured and indications of interactions with a second-generation particle (muons, $\mu$) are emerging. These tests reveal that the Higgs boson discovered ten years ago is remarkably consistent with the predictions of the theory and provide stringent constraints on many models of new phenomena beyond the Standard Model.

Jet substructure quantities are measured using jets groomed with the soft-drop grooming procedure in dijet events from 32.9 fb$^{-1}$ of $pp$ collisions collected with the ATLAS detector at $\sqrt{s} = 13$ TeV. These observables are sensitive to a wide range of QCD phenomena. Some observables, such as the jet mass and opening angle between the two subjets which pass the soft-drop condition, can be described by a high-order (resummed) series in the strong coupling constant $\alpha_S$. Other observables, such as the momentum sharing between the two subjets, are nearly independent of $\alpha_S$. These observables can be constructed using all interacting particles or using only charged particles reconstructed in the inner tracking detectors. Track-based versions of these observables are not collinear safe, but are measured more precisely, and universal non-perturbative functions can absorb the collinear singularities. The unfolded data are directly compared with QCD calculations and hadron-level Monte Carlo simulations. The measurements are performed in different pseudorapidity regions, which are then used to extract quark and gluon jet shapes using the predicted quark and gluon fractions in each region. All of the parton shower and analytical calculations provide an excellent description of the data in most regions of phase space.

The ratio of production cross sections of the W and Z bosons with exactly one associated jet is presented as a function of jet transverse momentum threshold. The measurement has been designed to maximise cancellation of experimental and theoretical uncertainties, and is reported both within a particle-level kinematic range corresponding to the detector acceptance and as a total cross-section ratio. Results are obtained with the ATLAS detector at the LHC in pp collisions at a centre-of-mass energy of 7 TeV using an integrated luminosity of 33 pb^-1. The results are compared with perturbative leading-order, leading-log, and next-to-leading-order QCD predictions, and are found to agree within experimental and theoretical uncertainties. The ratio is measured for events with a single jet with p_T > 30 GeV to be 8.73 +/- 0.30 (stat) +/- 0.40 (syst) in the electron channel, and $ 8.49 +/- 0.23 (stat) +/- 0.33 (syst) in the muon channel.

The ratio of the production cross sections for W and Z bosons in association with jets has been measured in proton-proton collisions at $\sqrt{s}$ = 7 TeV with the ATLAS experiment at the Large Hadron Collider. The measurement is based on the entire 2011 dataset, corresponding to an integrated luminosity of 4.6 fb$^{-1}$. Inclusive and differential cross-section ratios for massive vector bosons decaying to electrons and muons are measured in association with jets with transverse momentum $p_T$ > 30 GeV and jet rapidity $|y|$ < 4.4. The measurements are compared to next-to-leading-order perturbative QCD calculations and to predictions from different Monte Carlo generators implementing leading-order matrix elements supplemented by parton showers.

Jet substructure observables have significantly extended the search program for physics beyond the Standard Model at the Large Hadron Collider. The state-of-the-art tools have been motivated by theoretical calculations, but there has never been a direct comparison between data and calculations of jet substructure observables that are accurate beyond leading-logarithm approximation. Such observables are significant not only for probing the collinear regime of QCD that is largely unexplored at a hadron collider, but also for improving the understanding of jet substructure properties that are used in many studies at the Large Hadron Collider. This Letter documents a measurement of the first jet substructure quantity at a hadron collider to be calculated at next-to-next-to-leading-logarithm accuracy. The normalized, differential cross-section is measured as a function of log$_{10}\rho^2$, where $\rho$ is the ratio of the soft-drop mass to the ungroomed jet transverse momentum. This quantity is measured in dijet events from 32.9 fb$^{-1}$ of $\sqrt{s} = 13$ TeV proton-proton collisions recorded by the ATLAS detector. The data are unfolded to correct for detector effects and compared to precise QCD calculations and leading-logarithm particle-level Monte Carlo simulations.

A search is presented for the direct pair production of the stop, the supersymmetric partner of the top quark, that decays through an $R$-parity-violating coupling to a final state with two leptons and two jets, at least one of which is identified as a $b$-jet. The dataset corresponds to an integrated luminosity of 36.1 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of $\sqrt{s} = 13$ TeV, collected in 2015 and 2016 by the ATLAS detector at the LHC. No significant excess is observed over the Standard Model background, and exclusion limits are set on stop pair production at a 95% confidence level. Lower limits on the stop mass are set between 600 GeV and 1.5 TeV for branching ratios above 10% for decays to an electron or muon and a $b$-quark.

A search for new particles that decay into top quark pairs is reported. The search is performed with the ATLAS experiment at the LHC using an integrated luminosity of 20.3 fb$^{-1}$ of proton-proton collision data collected at a centre-of-mass energy of $\sqrt{s}=8$ TeV. The lepton-plus-jets final state is used, where the top pair decays to $W^+bW^-\bar{b}$, with one $W$ boson decaying leptonically and the other hadronically. The invariant mass spectrum of top quark pairs is examined for local excesses or deficits that are inconsistent with the Standard Model predictions. No evidence for a top quark pair resonance is found, and 95% confidence-level limits on the production rate are determined for massive states in benchmark models. The upper limits on the cross-section times branching ratio of a narrow $Z'$ boson decaying to top pairs range from 4.2 pb to 0.03 pb for resonance masses from 0.4 TeV to 3.0 TeV. A narrow leptophobic topcolour $Z'$ boson with mass below 1.8 TeV is excluded. Upper limits are set on the cross-section times branching ratio for a broad colour-octet resonance with $\Gamma/m =$ 15% decaying to $t\bar{t}$. These range from 2.5 pb to 0.03 pb for masses from 0.4 TeV to 3.0 TeV. A Kaluza-Klein excitation of the gluon in a Randall-Sundrum model is excluded for masses below 2.2 TeV.

This search, a type not previously performed at ATLAS, uses a comparison of the production cross sections for $e^+ \mu^-$ and $e^- \mu^+$ pairs to constrain physics processes beyond the Standard Model. It uses $139 \text{fb}^{-1}$ of proton$-$proton collision data recorded at $\sqrt{s} = 13$ TeV at the LHC. Targeting sources of new physics which prefer final states containing $e^{+}\mu^{-}$ to $e^{-}\mu^{+}$, the search contains two broad signal regions which are used to provide model-independent constraints on the ratio of cross sections at the 2% level. The search also has two special selections targeting supersymmetric models and leptoquark signatures. Observations using one of these selections are able to exclude, at 95% confidence level, singly produced smuons with masses up to 640 GeV in a model in which the only other light sparticle is a neutralino when the $R$-parity-violating coupling $\lambda'_{231}$ is close to unity. Observations using the other selection exclude scalar leptoquarks with masses below 1880 GeV when $g_{\text{1R}}^{eu}=g_{\text{1R}}^{\mu c}=1$, at 95% confidence level. The limit on the coupling reduces to $g_{\text{1R}}^{eu}=g_{\text{1R}}^{\mu c}=0.46$ for a mass of 1420 GeV.

A search for heavy Higgs bosons produced in association with a vector boson and decaying into a pair of vector bosons is performed in final states with two leptons (electrons or muons) of the same electric charge, 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 observed data are in agreement with Standard Model background expectations. The results are interpreted using higher-dimensional operators in an effective field theory. Upper limits on the production cross-section are calculated at 95% confidence level as a function of the heavy Higgs boson's mass and coupling strengths to vector bosons. Limits are set in the Higgs boson mass range from 300 to 1500 GeV, and depend on the assumed couplings. The highest excluded mass for a heavy Higgs boson with the coupling combinations explored is 900 GeV. Limits on coupling strengths are also provided.

A search for high-mass resonances decaying to $\tau\nu$ using proton-proton collisions at $\sqrt{s}$ = 13 TeV produced by the Large Hadron Collider is presented. Only $\tau$-lepton decays with hadrons in the final state are considered. The data were recorded with the ATLAS detector and correspond to an integrated luminosity of 36.1 fb$^{-1}$. No statistically significant excess above the Standard Model expectation is observed; model-independent upper limits are set on the visible $\tau\nu$ production cross section. Heavy $W^{\prime}$ bosons with masses less than 3.7 TeV in the Sequential Standard Model and masses less than 2.2-3.8 TeV depending on the coupling in the non-universal G(221) model are excluded at the 95% credibility level.

A search for high-mass resonances decaying into $\tau^{+}\tau^{-}$ final states using proton-proton collisions at $\sqrt{s}= 8$ TeV produced by the Large Hadron Collider is presented. The data were recorded with the ATLAS detector and correspond to an integrated luminosity of 19.5-20.3 fb$^{-1}$. No statistically significant excess above the Standard Model expectation is observed; 95% credibility upper limits are set on the cross section times branching fraction of $Z^{\prime}$ resonances decaying into $\tau^+\tau^-$ pairs as a function of the resonance mass. As a result, $Z^{\prime}$ bosons of the Sequential Standard Model with masses less than 2.02 TeV are excluded at 95% credibility. The impact of the fermionic couplings on the $Z^{\prime}$ acceptance is investigated and limits are also placed on a $Z^{\prime}$ model that exhibits enhanced couplings to third-generation fermions.

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

A search for resonant and non-resonant pair production of Higgs bosons in the $b\bar{b}\tau^+\tau^-$ final state is presented. The search uses 36.1 fb$^{-1}$ of $pp$ collision data with $\sqrt{s}= 13$ TeV recorded by the ATLAS experiment at the LHC in 2015 and 2016. The semileptonic and fully hadronic decays of the $\tau$-lepton pair are considered. No significant excess above the expected background is observed in the data. The cross-section times branching ratio for non-resonant Higgs boson pair production is constrained to be less than 30.9 fb, 12.7 times the Standard Model expectation, at 95% confidence level. The data are also analyzed to probe resonant Higgs boson pair production, constraining a model with an extended Higgs sector based on two doublets and a Randall-Sundrum bulk graviton model. Upper limits are placed on the resonant Higgs boson pair production cross-section times branching ratio, excluding resonances $X$ in the mass range $305~{\rm GeV} < m_X < 402~{\rm GeV}$ in the simplified hMSSM minimal supersymmetric model for $\tan\beta=2$ and excluding bulk Randall-Sundrum gravitons $G_{\mathrm{KK}}$ in the mass range $325~{\rm GeV} < m_{G_{\mathrm{KK}}} < 885~{\rm GeV}$ for $k/\overline{M}_{\mathrm{Pl}} = 1$.

A search for the $Z\gamma$ decay of the Higgs boson, with $Z$ boson decays into pairs of electrons or muons is presented. The analysis uses proton$-$proton collision data at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS detector at the Large Hadron Collider. The observed data are consistent with the expected background with a $p$-value of 1.3%. An upper limit at 95% confidence level on the production cross-section times the branching ratio for $pp\to H\to Z\gamma$ is set at 3.6 times the Standard Model prediction while 2.6 times is expected in the presence of the Standard Model Higgs boson. The best-fit value for the signal yield normalised to the Standard Model prediction is $2.0^{+1.0}_{-0.9}$ where the statistical component of the uncertainty is dominant.

Statistical combinations of searches for charginos and neutralinos using various decay channels are performed using $139\,$fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13\,$TeV with the ATLAS detector at the Large Hadron Collider. Searches targeting pure-wino chargino pair production, pure-wino chargino-neutralino production, or higgsino production decaying via Standard Model $W$, $Z$, or $h$ bosons are combined to extend the mass reach to the produced SUSY particles by 30-100 GeV. The depth of the sensitivity of the original searches is also improved by the combinations, lowering the 95% CL cross-section upper limits by 15%-40%.

This paper reviews and extends searches for the direct pair production of the scalar supersymmetric partners of the top and bottom quarks in proton-proton collisions collected by the ATLAS collaboration during the LHC Run 1. Most of the analyses use 20 fb$^{-1}$ of collisions at a centre-of-mass energy of $\sqrt{s}$ = 8 TeV, although in some case an additional 4.7 fb$^{-1}$ of collision data at $\sqrt{s}$ = 7 TeV are used. New analyses are introduced to improve the sensitivity to specific regions of the model parameter space. Since no evidence of third-generation squarks is found, exclusion limits are derived by combining several analyses and are presented in both a simplified model framework, assuming simple decay chains, as well as within the context of more elaborate phenomenological supersymmetric models.

Measurements are presented of the properties of high transverse momentum jets, produced in proton-proton collisions at a center-of-mass energy of sqrt(s) = 7 TeV. The data correspond to an integrated luminosity of 35 pb^-1 and were collected with the ATLAS detector in 2010. Jet mass, width, eccentricity, planar flow and angularity are measured for jets reconstructed using the anti-kt algorithm with distance parameters R = 0.6 and 1.0, with transverse momentum pT > 300 GeV and pseudorapidity |eta| < 2. The measurements are compared to the expectations of Monte Carlo generators that match leading-logarithmic parton showers to leading-order, or next-to-leading-order, matrix elements. The generators describe the general features of the jets, although discrepancies are observed in some distributions.

An analysis is presented of events containing jets including at least one $b$-tagged jet, sizeable missing transverse momentum, and at least two leptons including a pair of the same electric charge, with the scalar sum of the jet and lepton transverse momenta being large. A data sample with an integrated luminosity of 20.3 fb$^{-1}$ of $pp$ collisions at $\sqrt{s} = 8$ TeV recorded by the ATLAS detector at the Large Hadron Collider is used. Standard Model processes rarely produce these final states, but there are several models of physics beyond the Standard Model that predict an enhanced rate of production of such events; the ones considered here are production of vector-like quarks, enhanced four-top-quark production, pair production of chiral $b^\prime$-quarks, and production of two positively charged top quarks. Eleven signal regions are defined; subsets of these regions are combined when searching for each class of models. In the three signal regions primarily sensitive to positively charged top quark pair production, the data yield is consistent with the background expectation. There are more data events than expected from background in the set of eight signal regions defined for searching for vector-like quarks and chiral $b^\prime$-quarks, but the significance of the discrepancy is less than two standard deviations. The discrepancy reaches 2.5 standard deviations in the set of five signal regions defined for searching for four-top-quark production. The results are used to set 95% CL limits on various models.

A search is presented for a heavy resonance $Y$ decaying into a Standard Model Higgs boson $H$ and a new particle $X$ in a fully hadronic final state. The full Large Hadron Collider Run 2 dataset of proton-proton collisions at $\sqrt{s}= 13$ TeV collected by the ATLAS detector from 2015 to 2018 is used, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets the high $Y$-mass region, where the $H$ and $X$ have a significant Lorentz boost in the laboratory frame. A novel signal region is implemented using anomaly detection, where events are selected solely because of their incompatibility with a learned background-only model. It is defined using a jet-level tagger for signal-model-independent selection of the boosted $X$ particle, representing the first application of fully unsupervised machine learning to an ATLAS analysis. Two additional signal regions are implemented to target a benchmark $X$ decay into two quarks, covering topologies where the $X$ is reconstructed as either a single large-radius jet or two small-radius jets. The analysis selects Higgs boson decays into $b\bar{b}$, and a dedicated neural-network-based tagger provides sensitivity to the boosted heavy-flavor topology. No significant excess of data over the expected background is observed, and the results are presented as upper limits on the production cross section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) for signals with $m_Y$ between 1.5 and 6 TeV and $m_X$ between 65 and 3000 GeV.

Measurements of the centrality and rapidity dependence of inclusive jet production in $\sqrt{s_\mathrm{NN}} = 5.02$ TeV proton--lead ($p$+Pb) collisions and the jet cross-section in $\sqrt{s} = 2.76$ TeV proton--proton collisions are presented. These quantities are measured in datasets corresponding to an integrated luminosity of 27.8 nb$^{-1}$ and 4.0 pb$^{-1}$, respectively, recorded with the ATLAS detector at the Large Hadron Collider in 2013. The $p$+Pb collision centrality was characterised using the total transverse energy measured in the pseudorapidity interval $-4.9 < \eta < -3.2$ in the direction of the lead beam. Results are presented for the double-differential per-collision yields as a function of jet rapidity and transverse momentum ($p_\mathrm{T}$) for minimum-bias and centrality-selected $p$+Pb collisions, and are compared to the jet rate from the geometric expectation. The total jet yield in minimum-bias events is slightly enhanced above the expectation in a $p_\mathrm{T}$-dependent manner but is consistent with the expectation within uncertainties. The ratios of jet spectra from different centrality selections show a strong modification of jet production at all $p_\mathrm{T}$ at forward rapidities and for large $p_\mathrm{T}$ at mid-rapidity, which manifests as a suppression of the jet yield in central events and an enhancement in peripheral events. These effects imply that the factorisation between hard and soft processes is violated at an unexpected level in proton-nucleus collisions. Furthermore, the modifications at forward rapidities are found to be a function of the total jet energy only, implying that the violations may have a simple dependence on the hard parton-parton kinematics.

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.

Measurements of distributions of charged particles produced in proton-proton collisions with a centre-of-mass energy of 13 TeV are presented. The data were recorded by the ATLAS detector at the LHC and correspond to an integrated luminosity of 151 $\mu$b$^{-1}$. The particles are required to have a transverse momentum greater than 100 MeV and an absolute pseudorapidity less than 2.5. The charged-particle multiplicity, its dependence on transverse momentum and pseudorapidity and the dependence of the mean transverse momentum on multiplicity are measured in events containing at least two charged particles satisfying the above kinematic criteria. The results are corrected for detector effects and compared to the predictions from several Monte Carlo event generators.

Charged-particle distributions are measured in proton-proton collisions at a centre-of-mass energy of 13 TeV, using a data sample of nearly 9 million events, corresponding to an integrated luminosity of 170 $\mu$b$^{-1}$, recorded by the ATLAS detector during a special Large Hadron Collider fill. The charged-particle multiplicity, its dependence on transverse momentum and pseudorapidity and the dependence of the mean transverse momentum on the charged-particle multiplicity are presented. The measurements are performed with charged particles with transverse momentum greater than 500 MeV and absolute pseudorapidity less than 2.5, in events with at least one charged particle satisfying these kinematic requirements. Additional measurements in a reduced phase space with absolute pseudorapidity less than 0.8 are also presented, in order to compare with other experiments. The results are corrected for detector effects, presented as particle-level distributions and are compared to the predictions of various Monte Carlo event generators.