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.

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

Infrared and collinear safe event shape distributions and their mean values are determined in e+e- collisions at centre-of-mass energies between 45 and 202 GeV. A phenomenological analysis based on power correction models including hadron mass effects for both differential distributions and mean values is presented. Using power corrections, alpha_s is extracted from the mean values and shapes. In an alternative approach, renormalisation group invariance (RGI) is used as an explicit constraint, leading to a consistent description of mean values without the need for sizeable power corrections. The QCD beta-function is precisely measured using this approach. From the DELPHI data on Thrust, including data from low energy experiments, one finds beta_0 = 7.86 +/- 0.32 for the one loop coefficient of the beta-function or, assuming QCD, n_f = 4.75 +/- 0.44 for the number of active flavours. These values agree well with the QCD expectation of beta_0=7.67 and n_f=5. A direct measurement of the full logarithmic energy slope excludes light gluinos with a mass below 5 GeV.

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.

We present an angular analysis of the $B^{+}\rightarrow K^{\ast+}(\rightarrow K_{S}^{0}\pi^{+})\mu^{+}\mu^{-}$ decay using 9$\,\mbox{fb}^{-1}$ of $pp$ collision data collected with the LHCb experiment. For the first time, the full set of CP-averaged angular observables is measured in intervals of the dimuon invariant mass squared. Local deviations from Standard Model predictions are observed, similar to those in previous LHCb analyses of the isospin-partner $B^{0}\rightarrow K^{\ast0}\mu^{+}\mu^{-}$ decay. The global tension is dependent on which effective couplings are considered and on the choice of theory nuisance parameters.

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

Inclusive distributions of charged particles in hadronic W decays are experimentally investigated using the statistics collected by the DELPHI experiment at LEP during 1997, 1998 and 1999, at centre-of-mass energies from 183 to around 200 GeV. The possible effects of interconnection between the hadronic decays of two Ws are not observed. Measurements of the average multiplicity for charged and identified particles in q qbar and WW events at centre-of-mass energies from 130 to 200 GeV and in W decays are presented. The results on the average multiplicity of identified particles and on the position xi^* of the maximum of the xi_p = -log(2p/sqrt(s)) distribution are compared with predictions of JETSET and MLLA calculations.

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.

This paper presents measurements of distributions of charged particles which are produced in proton--proton collisions at a centre-of-mass energy of $\sqrt{s} = 8$ TeV and recorded by the ATLAS detector at the LHC. A special dataset recorded in 2012 with a small number of interactions per beam crossing (below 0.004) and corresponding to an integrated luminosity of $160 \mathrm{\mu b^{-1}}$ was used. A minimum-bias trigger was utilised to select a data sample of more than 9 million collision events. The multiplicity, pseudorapidity, and transverse momentum distributions of charged particles are shown in different regions of kinematics and charged-particle multiplicity, including measurements of final states at high multiplicity. The results are corrected for detector effects and are compared to the predictions of various Monte Carlo event generator models which simulate the full hadronic final state.

Measurements are presented from proton-proton collisions at centre-of-mass energies of sqrt(s) = 0.9, 2.36 and 7 TeV recorded with the ATLAS detector at the LHC. Events were collected using a single-arm minimum-bias trigger. The charged-particle multiplicity, its dependence on transverse momentum and pseudorapidity and the relationship between the mean transverse momentum and charged-particle multiplicity are measured. Measurements in different regions of phase-space are shown, providing diffraction-reduced measurements as well as more inclusive ones. The observed distributions are corrected to well-defined phase-space regions, using model-independent corrections. The results are compared to each other and to various Monte Carlo models, including a new AMBT1 PYTHIA 6 tune. In all the kinematic regions considered, the particle multiplicities are higher than predicted by the Monte Carlo models. The central charged-particle multiplicity per event and unit of pseudorapidity, for tracks with pT >100 MeV, is measured to be 3.483 +- 0.009 (stat) +- 0.106 (syst) at sqrt(s) = 0.9 TeV and 5.630 +- 0.003 (stat) +- 0.169 (syst) at sqrt(s) = 7 TeV.

A combination of the searches for pair-produced vector-like partners of the top and bottom quarks in various decay channels ($T$$\rightarrow$$Zt/Wb/Ht$, $B$$\rightarrow$$Zb/Wt/Hb$) is performed using 36.1 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. The observed data are found to be in good agreement with the Standard Model background prediction in all individual searches. Therefore, combined 95% confidence-level upper limits are set on the production cross-section for a range of vector-like quark scenarios, significantly improving upon the reach of the individual searches. Model-independent limits are set assuming the vector-like quarks decay to Standard Model particles. A singlet $T$ is excluded for masses below 1.31 TeV and a singlet $B$ is excluded for masses below 1.22 TeV. Assuming a weak isospin $(T,B)$ doublet and $|V_{Tb}| \ll |V_{tB}|$, $T$ and $B$ masses below 1.37 TeV are excluded.

This article presents measurements of the $t$-channel single top-quark ($t$) and top-antiquark ($\bar{t}$) total production cross sections $\sigma(tq)$ and $\sigma(\bar{t}q)$, their ratio $R_{t}=\sigma(tq)/\sigma(\bar{t}q)$, and a measurement of the inclusive production cross section $\sigma(tq + \bar{t}q)$ in proton--proton collisions at $\sqrt{s} = 7$ TeV at the LHC. Differential cross sections for the $tq$ and $\bar{t}q$ processes are measured as a function of the transverse momentum and the absolute value of the rapidity of $t$ and $\bar{t}$, respectively. The analyzed data set was recorded with the ATLAS detector and corresponds to an integrated luminosity of 4.59 fb$^{-1}$. Selected events contain one charged lepton, large missing transverse momentum, and two or three jets. The cross sections are measured by performing a binned maximum-likelihood fit to the output distributions of neural networks. The resulting measurements are $\sigma(tq)= 46\pm 6\; \mathrm{pb}$, $\sigma(\bar{t}q)= 23 \pm 4\; \mathrm{pb}$, $R_{t}=2.04\pm 0.18$, and $\sigma(tq + \bar{t}q)= 68 \pm 8\; \mathrm{pb}$, consistent with the Standard Model expectation. The uncertainty on the measured cross sections is dominated by systematic uncertainties, while the uncertainty on $R_{t}$ is mainly statistical. Using the ratio of $\sigma(tq + \bar{t}q)$ to its theoretical prediction, and assuming that the top-quark-related CKM matrix elements obey the relation $|V_{tb}|\gg |V_{ts}|, |V_{td}|$, we determine $|V_{tb}|=1.02 \pm 0.07$.

An updated analysis using about 1.5 million events recorded at $\sqrt{s} = M_Z$ with the DELPHI detector in 1994 is presented. Eighteen infrared and collinear safe event shape observables are measured as a function of the polar angle of the thrust axis. The data are compared to theoretical calculations in ${\cal O} (\alpha_s^2)$ including the event orientation. A combined fit of $\alpha_s$ and of the renormalization scale $x_{\mu}$ in $\cal O(\alpha_s^2$) yields an excellent description of the high statistics data. The weighted average from 18 observables including quark mass effects and correlations is $\alpha_s(M_Z^2) = 0.1174 \pm 0.0026$. The final result, derived from the jet cone energy fraction, the observable with the smallest theoretical and experimental uncertainty, is $\alpha_s(M_Z^2) = 0.1180 \pm 0.0006 (exp.) \pm 0.0013 (hadr.) \pm 0.0008 (scale) \pm 0.0007 (mass)$. Further studies include an $\alpha_s$ determination using theoretical predictions in the next-to-leading log approximation (NLLA), matched NLLA and $\cal O(\alpha_s^2$) predictions as well as theoretically motivated optimized scale setting methods. The influence of higher order contributions was also investigated by using the method of Pad\'{e} approximants. Average $\alpha_s$ values derived from the different approaches are in good agreement.

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.

Charge-dependent azimuthal correlations of same- and opposite-sign pairs with respect to the second- and third-order event planes have been measured in pPb collisions at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV and PbPb collisions at 5.02 TeV with the CMS experiment at the LHC. The measurement is motivated by the search for the charge separation phenomenon predicted by the chiral magnetic effect (CME) in heavy ion collisions. Three- and two-particle azimuthal correlators are extracted as functions of the pseudorapidity difference, the transverse momentum ($p_\mathrm{t}$) difference, and the $p_\mathrm{t}$ average of same- and opposite-charge pairs in various event multiplicity ranges. The data suggest that the charge-dependent three-particle correlators with respect to the second- and third-order event planes share a common origin, predominantly arising from charge-dependent two-particle azimuthal correlations coupled with an anisotropic flow. The CME is expected to lead to a $v_2$-independent three-particle correlation when the magnetic field is fixed. Using an event shape engineering technique, upper limits on the $v_2$-independent fraction of the three-particle correlator are estimated to be 13% for pPb and 7% for PbPb collisions at 95% confidence level. The results of this analysis, both the dominance of two-particle correlations as a source of the three-particle results and the similarities seen between PbPb and pPb, provide stringent constraints on the origin of charge-dependent three-particle azimuthal correlations and challenge their interpretation as arising from a chiral magnetic effect in heavy ion collisions.

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.