High-energy nuclear collisions create a quark-gluon plasma, whose initial condition and subsequent expansion vary from event to event, impacting the distribution of the event-wise average transverse momentum ($P([p_{\mathrm{T}}])$). Distinguishing between contributions from fluctuations in the size of the nuclear overlap area (geometrical component) and other sources at fixed size (intrinsic component) presents a challenge. Here, these two components are distinguished by measuring the mean, variance, and skewness of $P([p_{\mathrm{T}}])$ in $^{208}$Pb+$^{208}$Pb and $^{129}$Xe+$^{129}$Xe collisions at $\sqrt{s_{{\mathrm{NN}}}} = 5.02$ and 5.44 TeV, respectively, using the ATLAS detector at the LHC. All observables show distinct changes in behavior in ultra-central collisions, where the geometrical variations are suppressed as the overlap area reaches its maximum. These results demonstrate a new technique to disentangle geometrical and intrinsic fluctuations, enabling constraints on initial condition and properties of the quark-gluon plasma, such as the speed of sound.

A search for decays of the Higgs boson into a $Z$ boson and a light resonance, with a mass of 0.5-3.5 GeV, is performed using the full 140 fb$^{-1}$ dataset of 13 TeV proton-proton collisions recorded by the ATLAS detector during Run~2 of the LHC. Leptonic decays of the $Z$ boson and hadronic decays of the light resonance are considered. The resonance can be interpreted as a $J/\psi$ or $\eta_c$ meson, an axion-like particle, or a light pseudoscalar in two-Higgs-doublet models. Due to its low mass, it would be produced with high boost and reconstructed as a single small-radius jet of hadrons. A neural network is used to correct the Monte Carlo simulation of the background in a data-driven way. Two additional neural networks are used to distinguish signal from background. A binned profile-likelihood fit is performed on the final-state invariant mass distribution. No significant excess of events relative to the expected background is observed, and upper limits at 95% confidence level are set on the Higgs boson's branching fraction to a $Z$ boson and a light resonance. The exclusion limit is 10% for the lower masses, and increases for higher masses. Upper limits on the effective coupling $C^\text{eff}_{ZH}/\Lambda$ of an axion-like particle to a Higgs boson and $Z$ boson are also set at 95% confidence level, and range from 0.9 to 2 TeV$^{-1}$.

A search for Higgs boson pair production via vector-boson fusion is performed in the Lorentz-boosted regime, where a Higgs boson candidate is reconstructed as a single large-radius jet, using 140 fb$^{-1}$ of proton-proton collision data at $\sqrt{s} = 13$ TeV recorded by the ATLAS detector at the Large Hadron Collider. Only Higgs boson decays into bottom quark pairs are considered. The search is particularly sensitive to the quartic coupling between two vector bosons and two Higgs bosons relative to its Standard Model prediction, $\kappa_{2V}$. This study constrains $\kappa_{2V}$ to $0.55 < \kappa_{2V} < 1.49$ at 95% confidence level. The value $\kappa_{2V} = 0$ is excluded with a significance of 3.8 standard deviations with other Higgs boson couplings fixed to their Standard Model values. A search for new heavy spin-0 resonances that would mediate Higgs boson pair production via vector-boson fusion is carried out in the mass range of 1-5 TeV for the first time under several model and decay-width assumptions. No significant deviation from the Standard Model hypothesis is observed and exclusion limits at 95% confidence level are derived.

This paper presents a first measurement of the cross-section for the charged-current Drell-Yan process $pp\rightarrow W^{\pm} \rightarrow \ell^{\pm} \nu$ above the resonance region, where $\ell$ is an electron or muon. The measurement is performed for transverse masses, $m_{\text{T}}^{\text{W}}$, between 200 GeV and 5000 GeV, using a sample of 140~fb$^{-1}$ of $pp$ collision data at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV collected by the ATLAS detector at the LHC during 2015-2018. The data are presented single differentially in transverse mass and double differentially in transverse mass and absolute lepton pseudorapidity. A test of lepton flavour universality shows no significant deviations from the Standard Model. The electron and muon channel measurements are combined to achieve a total experimental precision of 3% at low $m_{\text{T}}^{\text{W}}$. The single- and double differential $W$-boson charge asymmetries are evaluated from the measurements. A comparison to next-to-next-to-leading-order perturbative QCD predictions using several recent parton distribution functions and including next-to-leading-order electroweak effects indicates the potential of the data to constrain parton distribution functions. The data are also used to constrain four fermion operators in the Standard Model Effective Field Theory formalism, in particular the lepton-quark operator Wilson coefficient $c_{\ell q}^{(3)}.$

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-sections for the production of a $Z$ boson in association with two photons are measured in proton$-$proton collisions at a centre-of-mass energy of 13 TeV. The data used correspond to an integrated luminosity of 139 fb$^{-1}$ recorded by the ATLAS experiment during Run 2 of the LHC. The measurements use the electron and muon decay channels of the $Z$ boson, and a fiducial phase-space region where the photons are not radiated from the leptons. The integrated $Z(\rightarrow\ell\ell)\gamma\gamma$ cross-section is measured with a precision of 12% and differential cross-sections are measured as a function of six kinematic variables of the $Z\gamma\gamma$ system. The data are compared with predictions from MC event generators which are accurate to up to next-to-leading order in QCD. The cross-section measurements are used to set limits on the coupling strengths of dimension-8 operators in the framework of an effective field theory.

Differential cross-sections are measured for the production of four charged leptons in association with two jets. These measurements are sensitive to final states in which the jets are produced via the strong interaction as well as to the purely-electroweak vector boson scattering process. The analysis is performed using proton-proton collision data collected by ATLAS at $\sqrt{s}=13$ TeV and with an integrated luminosity of 140 fb$^{-1}$. The data are corrected for the effects of detector inefficiency and resolution and are compared to state-of-the-art Monte Carlo event generator predictions. The differential cross-sections are used to search for anomalous weak-boson self-interactions that are induced by dimension-six and dimension-eight operators in Standard Model effective field theory.

Studies of the correlations of the two highest transverse momentum (leading) jets in individual Pb+Pb collision events can provide information about the mechanism of jet quenching by the hot and dense matter created in such collisions. In Pb+Pb and pp collisions at $\sqrt{s_{_\text{NN}}}$ = 5.02 TeV, measurements of the leading dijet transverse momentum ($p_{\mathrm{T}}$) correlations are presented. Additionally, measurements in Pb+Pb collisions of the dijet pair nuclear modification factors projected along leading and subleading jet $p_{\mathrm{T}}$ are made. The measurements are performed using the ATLAS detector at the LHC with 260 pb$^{-1}$ of pp data collected in 2017 and 2.2 nb$^{-1}$ of Pb+Pb data collected in 2015 and 2018. An unfolding procedure is applied to the two-dimensional leading and subleading jet $p_{\mathrm{T}}$ distributions to account for experimental effects in the measurement of both jets. Results are provided for dijets with leading jet $p_{\mathrm{T}}$ greater than 100 GeV. Measurements of the dijet-yield-normalized $x_{\mathrm{J}}$ distributions in Pb+Pb collisions show an increased fraction of imbalanced jets compared to pp collisions; these measurements are in agreement with previous measurements of the same quantity at 2.76 TeV in the overlapping kinematic range. Measurements of the absolutely-normalized dijet rate in Pb+Pb and pp collisions are also presented, and show that balanced dijets are significantly more suppressed than imbalanced dijets in Pb+Pb collisions. It is observed in the measurements of the pair nuclear modification factors that the subleading jets are significantly suppressed relative to leading jets with $p_{\mathrm{T}}$ between 100 and 316 GeV for all centralities in Pb+Pb collisions.

Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

The structure of nucleons is multidimensional and depends on the transverse momenta, spatial geometry, and polarization of the constituent partons. Such a structure can be studied using high-energy photons produced in ultraperipheral heavy-ion collisions. The first measurement of the azimuthal angular correlations of exclusively produced events with two jets in photon-lead interactions at large momentum transfer is presented, a process that is considered to be sensitive to the underlying nuclear gluon polarization. This study uses a data sample of ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, corresponding to an integrated luminosity of 0.38 nb$^{-1}$, collected with the CMS experiment at the LHC. The measured second harmonic of the correlation between the sum and difference of the two jet momenta is found to be positive, and rising, as the dijet momentum increases. A well-tuned model that has been successful at describing a wide range of proton scattering data from the HERA experiments fails to describe the observed correlations, suggesting the presence of gluon polarization effects.