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.

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.

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.

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