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.

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 measurement of the inclusive cross section of top quark pair production in association with a Z boson using proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC is performed. The data sample corresponds to an integrated luminosity of 77.5 fb$^{-1}$, collected by the CMS experiment during 2016 and 2017. The measurement is performed using final states containing three or four charged leptons (electrons or muons), and the Z boson is detected through its decay to an oppositely charged lepton pair. The production cross section is measured to be $\sigma(\mathrm{t\bar{t}Z})$ $=$ 0.95 $\pm$ 0.05 (stat) $\pm$ 0.06 (syst) pb. For the first time, differential cross sections are measured as functions of the transverse momentum of the Z boson and the angular distribution of the negatively charged lepton from the Z boson decay. The most stringent direct limits to date on the anomalous couplings of the top quark to the Z boson are presented, including constraints on the Wilson coefficients in the framework of the standard model effective field theory.

This paper presents a search for a new Z' vector gauge boson with the ATLAS experiment at the Large Hadron Collider using pp collision data collected at $\sqrt{s} = 13$ TeV, corresponding to an integrated luminosity of 139 fb$^{-1}$. The new gauge boson Z' is predicted by $L_{\mu}-L_{\tau}$ models to address observed phenomena that can not be explained by the Standard Model. The search examines the four-muon (4$\mu$) final state, using a deep learning neural network classifier to separate the Z' signal from the Standard Model background events. The di-muon invariant masses in the $4\mu$ events are used to extract the Z' resonance signature. No significant excess of events is observed over the predicted background. Upper limits at a 95% confidence level on the Z' production cross-section times the decay branching fraction of $pp \rightarrow Z'\mu\mu \rightarrow 4\mu$ are set from 0.31 to 4.3 fb for the Z' mass ranging from 5 to 81 GeV. The corresponding common coupling strengths, $g_{Z'}$, of the Z' boson to the second and third generation leptons above 0.003 - 0.2 have been excluded.

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

The measurement of hard scatterings in proton-nucleus collisions has resulted in a greater understanding of both the proton and nuclear structure. ATLAS measured the centrality dependence of the dijet yield using 165 nb$^{-1}$ of $p$+Pb data collected at $\sqrt{s_{_\text{NN}}}$ = 8.16 TeV in 2016. The event centrality, which reflects the $p$+Pb impact parameter, is characterized by the total transverse energy registered in the Pb-going side of the forward calorimeter. The central-to-peripheral ratio of the scaled dijet yields, $R_\mathrm{CP}$, is evaluated, and the results are presented as a function of variables that reflect the kinematics of the initial hard parton scattering process. The $R_\mathrm{CP}$ shows a scaling with the Bjorken-$x$ of the parton originating from the proton, $x_p$, while no such trend is observed as a function of $x_\mathrm{Pb}$. This analysis provides unique input to understanding the role of small proton spatial configurations in $p$+Pb collisions by covering parton momentum fractions from the valence region down to $x_p \sim 10^{-3}$ and $x_\mathrm{Pb}\sim 4\cdot10^{-4}$.

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.

A search for nonresonant Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is presented. The analysis uses 126 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}={13}$ TeV collected with the ATLAS detector at the Large Hadron Collider, and targets both the gluon-gluon fusion and vector-boson fusion production modes. No evidence of the signal is found and the observed (expected) upper limit on the cross-section for nonresonant Higgs boson pair production is determined to be 5.4 (8.1) times the Standard Model predicted cross-section at 95% confidence level. Constraints are placed on modifiers to the $HHH$ and $HHVV$ couplings. The observed (expected) $2\sigma$ constraints on the $HHH$ coupling modifier, $\kappa_\lambda$, are determined to be $[-3.5, 11.3]$ ($[-5.4, 11.4]$), while the corresponding constraints for the $HHVV$ coupling modifier, $\kappa_{2V}$, are $[-0.0, 2.1]$ ($[-0.1, 2.1]$). In addition, constraints on relevant coefficients are derived in the context of the Standard Model effective field theory and Higgs effective field theory, and upper limits on the $HH$ production cross-section are placed in seven Higgs effective field theory benchmark scenarios.

A search is conducted for a new scalar boson $S$, with a mass distinct from that of the Higgs boson, decaying into four leptons ($\ell =$$e$, $\mu$) via an intermediate state containing two on-shell, promptly decaying new spin-1 bosons $Z_\text{d}$: $S \rightarrow Z_\text{d}Z_\text{d} \rightarrow 4\ell$, where the $Z_\text{d}$ boson has a mass between 15 and 300 GeV, and the $S$ boson has a mass between either 30 and 115 GeV or 130 and 800 GeV. The search uses proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider with an integrated luminosity of 139 fb$^{-1}$ at a centre-of-mass energy of $\sqrt{s}=13$ TeV. No significant excess above the Standard Model background expectation is observed. Upper limits at 95% confidence level are set on the production cross-section times branching ratio, $\sigma(gg \to S) \times \mathcal{B}(S\rightarrow Z_\text{d}Z_\text{d} \rightarrow 4\ell)$, as a function of the mass of both particles, $m_S$ and $m_{Z\text{d}}$.