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.

A search for new phenomena is performed in final states containing one or more jets and an imbalance in transverse momentum in pp collisions at a centre-of-mass energy of 13 TeV. The analysed data sample, recorded with the CMS detector at the CERN LHC, corresponds to an integrated luminosity of 2.3 inverse femtobarns. Several kinematic variables are employed to suppress the dominant background, multijet production, as well as to discriminate between other standard model and new physics processes. The search provides sensitivity to a broad range of new-physics models that yield a stable weakly interacting massive particle. The number of observed candidate events is found to agree with the expected contributions from standard model processes, and the result is interpreted in the mass parameter space of fourteen simplified supersymmetric models that assume the pair production of gluinos or squarks and a range of decay modes. For models that assume gluino pair production, masses up to 1575 and 975 GeV are excluded for gluinos and neutralinos, respectively. For models involving the pair production of top squarks and compressed mass spectra, top squark masses up to 400 GeV are excluded.

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

The Fourier coefficients v[2] and v[3] characterizing the anisotropy of the azimuthal distribution of charged particles produced in PbPb collisions at sqrt(s[NN]) = 5.02 TeV are measured with data collected by the CMS experiment. The measurements cover a broad transverse momentum range, 1 < pT < 100 GeV. The analysis focuses on pT > 10 GeV range, where anisotropic azimuthal distributions should reflect the path-length dependence of parton energy loss in the created medium. Results are presented in several bins of PbPb collision centrality, spanning the 60% most central events. The v[2] coefficient is measured with the scalar product and the multiparticle cumulant methods, which have different sensitivities to the initial-state fluctuations. The values of both methods remain positive up to pT of about 60-80 GeV, in all examined centrality classes. The v[3] coefficient, only measured with the scalar product method, tends to zero for pT greater than or equal to 20 GeV. Comparisons between theoretical calculations and data provide new constraints on the path-length dependence of parton energy loss in heavy ion collisions and highlight the importance of the initial-state fluctuations.

Azimuthal correlations between the two jets with the largest transverse momenta $ {p_{\mathrm{T}}} $ in inclusive 2-, 3-, and 4-jet events are presented for several regions of the leading jet $ {p_{\mathrm{T}}} $ up to 4 TeV. For 3- and 4-jet scenarios, measurements of the minimum azimuthal angles between any two of the three or four leading $ {p_{\mathrm{T}}} $ jets are also presented. The analysis is based on data from proton-proton collisions collected by the CMS Collaboration at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Calculations based on leading-order matrix elements supplemented with parton showering and hadronization do not fully describe the data, so next-to-leading-order calculations matched with parton shower and hadronization models are needed to better describe the measured distributions. Furthermore, we show that azimuthal jet correlations are sensitive to details of the parton showering, hadronization, and multiparton interactions. A next-to-leading-order calculation matched with parton showers in the MC@NLO method, as implemented in HERWIG 7, gives a better overall description of the measurements than the POWHEG method.

Measurements of primary charged hadron multiplicity distributions are presented for non-single-diffractive events in proton-proton collisions at centre-of-mass energies of sqrt(s) = 0.9, 2.36, and 7 TeV, in five pseudorapidity ranges from |eta|<0.5 to |eta|<2.4. The data were collected with the minimum-bias trigger of the CMS experiment during the LHC commissioning runs in 2009 and the 7 TeV run in 2010. The multiplicity distribution at sqrt(s) = 0.9 TeV is in agreement with previous measurements. At higher energies the increase of the mean multiplicity with sqrt(s) is underestimated by most event generators. The average transverse momentum as a function of the multiplicity is also presented. The measurement of higher-order moments of the multiplicity distribution confirms the violation of Koba-Nielsen-Olesen scaling that has been observed at lower energies.

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.

Constraints on the Higgs boson self-coupling are set by combining double-Higgs boson analyses in the $b\bar{b}b\bar{b}$, $b\bar{b}\tau^+\tau^-$ and $b\bar{b} \gamma \gamma$ decay channels with single-Higgs boson analyses targeting the $\gamma \gamma$, $ZZ^*$, $WW^*$, $\tau^+ \tau^-$ and $b\bar{b}$ decay channels. The data used in these analyses were recorded by the ATLAS detector at the LHC in proton$-$proton collisions at $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 126$-$139 fb$^{-1}$. The combination of the double-Higgs analyses sets an upper limit of $\mu_{HH} < 2.4$ at 95% confidence level on the double-Higgs production cross-section normalised to its Standard Model prediction. Combining the single-Higgs and double-Higgs analyses, with the assumption that new physics affects only the Higgs boson self-coupling ($\lambda_{HHH}$), values outside the interval $-0.4< \kappa_{\lambda}=(\lambda_{HHH}/\lambda_{HHH}^{\textrm{SM}})< 6.3$ are excluded at 95% confidence level. The combined single-Higgs and double-Higgs analyses provide results with fewer assumptions, by adding in the fit more coupling modifiers introduced to account for the Higgs boson interactions with the other Standard Model particles. In this relaxed scenario, the constraint becomes $-1.4 < \kappa_{\lambda} < 6.1$ at 95% CL.

This paper presents a statistical combination of searches targeting final states with two top quarks and invisible particles, characterised by the presence of zero, one or two leptons, at least one jet originating from a $b$-quark and missing transverse momentum. The analyses are searches for phenomena beyond the Standard Model consistent with the direct production of dark matter in $pp$ collisions at the LHC, using 139 fb$^{-\text{1}}$ of data collected with the ATLAS detector at a centre-of-mass energy of 13 TeV. The results are interpreted in terms of simplified dark matter models with a spin-0 scalar or pseudoscalar mediator particle. In addition, the results are interpreted in terms of upper limits on the Higgs boson invisible branching ratio, where the Higgs boson is produced according to the Standard Model in association with a pair of top quarks. For scalar (pseudoscalar) dark matter models, with all couplings set to unity, the statistical combination extends the mass range excluded by the best of the individual channels by 50 (25) GeV, excluding mediator masses up to 370 GeV. In addition, the statistical combination improves the expected coupling exclusion reach by 14% (24%), assuming a scalar (pseudoscalar) mediator mass of 10 GeV. An upper limit on the Higgs boson invisible branching ratio of 0.38 (0.30$^{+\text{0.13}}_{-\text{0.09}}$) is observed (expected) at 95% confidence level.