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.

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

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 letter presents a combination of searches for Higgs boson pair production using up to 36.1 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy $\sqrt{s} = 13$ TeV recorded with the ATLAS detector at the LHC. The combination is performed using six analyses searching for Higgs boson pairs decaying into the bbbb, bbWW, bb$\tau\tau$, WWWW, bb$\gamma \gamma$ and WW$\gamma\gamma$ final states. Results are presented for non-resonant and resonant Higgs boson pair production modes. No statistically significant excess in data above the Standard Model predictions is found. The combined observed (expected) limit at 95% confidence level on the non-resonant Higgs boson pair production cross-section is 6.9 (10) times the predicted Standard Model cross-section. Limits are also set on the ratio ($ \kappa_{\lambda} $) of the Higgs boson self-coupling to its Standard Model value. This ratio is constrained at 95% confidence level in observation (expectation) to $ -5.0 < \kappa_{\lambda} <12.0 $ ($ -5.8 < \kappa_{\lambda} <12.0 $). In addition, limits are set on the production of narrow scalar resonances and spin-2 Kaluza-Klein Randall-Sundrum gravitons. Exclusion regions are also provided in the parameter space of the habemus Minimal Supersymmetric Standard Model and the Electroweak Singlet Model.

Correlations of two flow harmonics $v_n$ and $v_m$ via three- and four-particle cumulants are measured in 13 TeV $pp$, 5.02 TeV $p$+Pb, and 2.76 TeV peripheral Pb+Pb collisions with the ATLAS detector at the LHC. The goal is to understand the multi-particle nature of the long-range collective phenomenon in these collision systems. The large non-flow background from dijet production present in the standard cumulant method is suppressed using a method of subevent cumulants involving two, three and four subevents separated in pseudorapidity. The results show a negative correlation between $v_2$ and $v_3$ and a positive correlation between $v_2$ and $v_4$ for all collision systems and over the full multiplicity range. However, the magnitudes of the correlations are found to depend strongly on the event multiplicity, the choice of transverse momentum range and collision system. The relative correlation strength, obtained by normalisation of the cumulants with the $\langle v_n^2\rangle$ from a two-particle correlation analysis, is similar in the three collision systems and depends weakly on the event multiplicity and transverse momentum. These results based on the subevent methods provide strong evidence of a similar long-range multi-particle collectivity in $pp$, $p$+Pb and peripheral Pb+Pb collisions.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

Cross-section measurements for a $Z$ boson produced in association with high-transverse-momentum jets ($p_{\mathrm{T}} \geq 100$ GeV) and decaying into a charged-lepton pair ($e^+e^-,\mu^+\mu^-$) are presented. The measurements are performed using proton-proton collisions at $\sqrt{s}=13$ TeV corresponding to an integrated luminosity of $139$ fb$^{-1}$ collected by the ATLAS experiment at the LHC. Measurements of angular correlations between the $Z$ boson and the closest jet are performed in events with at least one jet with $p_{\mathrm{T}} \geq 500$ GeV. Event topologies of particular interest are the collinear emission of a $Z$ boson in dijet events and a boosted $Z$ boson recoiling against a jet. Fiducial cross sections are compared with state-of-the-art theoretical predictions. The data are found to agree with next-to-next-to-leading-order predictions by NNLOjet and with the next-to-leading-order multi-leg generators MadGraph5_aMC@NLO and Sherpa.