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.

The ratio of production cross sections of the W and Z bosons with exactly one associated jet is presented as a function of jet transverse momentum threshold. The measurement has been designed to maximise cancellation of experimental and theoretical uncertainties, and is reported both within a particle-level kinematic range corresponding to the detector acceptance and as a total cross-section ratio. Results are obtained with the ATLAS detector at the LHC in pp collisions at a centre-of-mass energy of 7 TeV using an integrated luminosity of 33 pb^-1. The results are compared with perturbative leading-order, leading-log, and next-to-leading-order QCD predictions, and are found to agree within experimental and theoretical uncertainties. The ratio is measured for events with a single jet with p_T > 30 GeV to be 8.73 +/- 0.30 (stat) +/- 0.40 (syst) in the electron channel, and $ 8.49 +/- 0.23 (stat) +/- 0.33 (syst) in the muon channel.

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.

A search for new particles that decay into top quark pairs is reported. The search is performed with the ATLAS experiment at the LHC using an integrated luminosity of 20.3 fb$^{-1}$ of proton-proton collision data collected at a centre-of-mass energy of $\sqrt{s}=8$ TeV. The lepton-plus-jets final state is used, where the top pair decays to $W^+bW^-\bar{b}$, with one $W$ boson decaying leptonically and the other hadronically. The invariant mass spectrum of top quark pairs is examined for local excesses or deficits that are inconsistent with the Standard Model predictions. No evidence for a top quark pair resonance is found, and 95% confidence-level limits on the production rate are determined for massive states in benchmark models. The upper limits on the cross-section times branching ratio of a narrow $Z'$ boson decaying to top pairs range from 4.2 pb to 0.03 pb for resonance masses from 0.4 TeV to 3.0 TeV. A narrow leptophobic topcolour $Z'$ boson with mass below 1.8 TeV is excluded. Upper limits are set on the cross-section times branching ratio for a broad colour-octet resonance with $\Gamma/m =$ 15% decaying to $t\bar{t}$. These range from 2.5 pb to 0.03 pb for masses from 0.4 TeV to 3.0 TeV. A Kaluza-Klein excitation of the gluon in a Randall-Sundrum model is excluded for masses below 2.2 TeV.

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 search for high-mass resonances decaying into $\tau^{+}\tau^{-}$ final states using proton-proton collisions at $\sqrt{s}= 8$ TeV produced by the Large Hadron Collider is presented. The data were recorded with the ATLAS detector and correspond to an integrated luminosity of 19.5-20.3 fb$^{-1}$. No statistically significant excess above the Standard Model expectation is observed; 95% credibility upper limits are set on the cross section times branching fraction of $Z^{\prime}$ resonances decaying into $\tau^+\tau^-$ pairs as a function of the resonance mass. As a result, $Z^{\prime}$ bosons of the Sequential Standard Model with masses less than 2.02 TeV are excluded at 95% credibility. The impact of the fermionic couplings on the $Z^{\prime}$ acceptance is investigated and limits are also placed on a $Z^{\prime}$ model that exhibits enhanced couplings to third-generation fermions.

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.

An analysis is presented of events containing jets including at least one $b$-tagged jet, sizeable missing transverse momentum, and at least two leptons including a pair of the same electric charge, with the scalar sum of the jet and lepton transverse momenta being large. A data sample with an integrated luminosity of 20.3 fb$^{-1}$ of $pp$ collisions at $\sqrt{s} = 8$ TeV recorded by the ATLAS detector at the Large Hadron Collider is used. Standard Model processes rarely produce these final states, but there are several models of physics beyond the Standard Model that predict an enhanced rate of production of such events; the ones considered here are production of vector-like quarks, enhanced four-top-quark production, pair production of chiral $b^\prime$-quarks, and production of two positively charged top quarks. Eleven signal regions are defined; subsets of these regions are combined when searching for each class of models. In the three signal regions primarily sensitive to positively charged top quark pair production, the data yield is consistent with the background expectation. There are more data events than expected from background in the set of eight signal regions defined for searching for vector-like quarks and chiral $b^\prime$-quarks, but the significance of the discrepancy is less than two standard deviations. The discrepancy reaches 2.5 standard deviations in the set of five signal regions defined for searching for four-top-quark production. The results are used to set 95% CL limits on various models.

A search is presented for a heavy resonance $Y$ decaying into a Standard Model Higgs boson $H$ and a new particle $X$ in a fully hadronic final state. The full Large Hadron Collider Run 2 dataset of proton-proton collisions at $\sqrt{s}= 13$ TeV collected by the ATLAS detector from 2015 to 2018 is used, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets the high $Y$-mass region, where the $H$ and $X$ have a significant Lorentz boost in the laboratory frame. A novel signal region is implemented using anomaly detection, where events are selected solely because of their incompatibility with a learned background-only model. It is defined using a jet-level tagger for signal-model-independent selection of the boosted $X$ particle, representing the first application of fully unsupervised machine learning to an ATLAS analysis. Two additional signal regions are implemented to target a benchmark $X$ decay into two quarks, covering topologies where the $X$ is reconstructed as either a single large-radius jet or two small-radius jets. The analysis selects Higgs boson decays into $b\bar{b}$, and a dedicated neural-network-based tagger provides sensitivity to the boosted heavy-flavor topology. No significant excess of data over the expected background is observed, and the results are presented as upper limits on the production cross section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) for signals with $m_Y$ between 1.5 and 6 TeV and $m_X$ between 65 and 3000 GeV.

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.

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.

A combination of searches for resonant Higgs boson pair production is presented, using up to 139 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the LHC. The combination includes searches performed in three decay channels: $b\bar{b}b\bar{b}$, $bb\tau^+\tau^-$ and $bb\gamma\gamma$. No excess above the expected background is observed and upper limits are set at the 95% confidence level on the production cross section of Higgs boson pairs originating from the decay of a narrow scalar resonance with mass in the range 251 GeV-5 TeV. The observed (expected) limits are in the range 0.96-600 fb (1.2-390 fb). The limits are interpreted in the Type-I Two-Higgs-Doublet Model and the Minimimal Supersymmetric Standard Model, and constrain parameter space not previously excluded by other searches.

Parton energy loss in the quark-gluon plasma (QGP) is studied with a measurement of photon-tagged jet production in 1.7 nb$^{-1}$ of Pb+Pb data and 260 pb$^{-1}$ of $pp$ data, both at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV, with the ATLAS detector. The process $pp \to \gamma$+jet+$X$ and its analogue in Pb+Pb collisions is measured in events containing an isolated photon with transverse momentum ($p_\mathrm{T}$) above $50$ GeV and reported as a function of jet $p_\mathrm{T}$. This selection results in a sample of jets with a steeply falling $p_\mathrm{T}$ distribution that are mostly initiated by the showering of quarks. The $pp$ and Pb+Pb measurements are used to report the nuclear modification factor, $R_\mathrm{AA}$, and the fractional energy loss, $S_\mathrm{loss}$, for photon-tagged jets. In addition, the results are compared with the analogous ones for inclusive jets, which have a significantly smaller quark-initiated fraction. The $R_\mathrm{AA}$ and $S_\mathrm{loss}$ values are found to be significantly different between those for photon-tagged jets and inclusive jets, demonstrating that energy loss in the QGP is sensitive to the colour-charge of the initiating parton. The results are also compared with a variety of theoretical models of colour-charge-dependent energy loss.

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 article presents the results of two studies of Higgs boson properties using the $WW^*(\rightarrow e\nu\mu\nu)jj$ final state, based on a dataset corresponding to 36.1/fb of $\sqrt{s}$=13 TeV proton$-$proton collisions recorded by the ATLAS experiment at the Large Hadron Collider. The first study targets Higgs boson production via gluon$-$gluon fusion and constrains the CP properties of the effective Higgs$-$gluon interaction. Using angular distributions and the overall rate, a value of $\tan(\alpha) = 0.0 \pm 0.4$ stat. $ \pm 0.3$ syst is obtained for the tangent of the mixing angle for CP-even and CP-odd contributions. The second study exploits the vector-boson fusion production mechanism to probe the Higgs boson couplings to longitudinally and transversely polarised $W$ and $Z$ bosons in both the production and the decay of the Higgs boson; these couplings have not been directly constrained previously. The polarisation-dependent coupling-strength scale factors are defined as the ratios of the measured polarisation-dependent coupling strengths to those predicted by the Standard Model, and are determined using rate and kinematic information to be $a_L=0.91^{+0.10}_{-0.18}$(stat.)$^{+0.09}_{-0.17}$(syst.) and $a_{T}=1.2 \pm 0.4 $(stat.)$ ^{+0.2}_{-0.3} $(syst.). These coupling strengths are translated into pseudo-observables, resulting in $\kappa_{VV}= 0.91^{+0.10}_{-0.18}$(stat.)$^{+0.09}_{-0.17}$(syst.) and $\epsilon_{VV} =0.13^{+0.28}_{-0.20}$ (stat.)$^{+0.08}_{-0.10}$(syst.). All results are consistent with the Standard Model predictions.

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.