The coherent J/$\psi$ photoproduction cross section is measured for the first time at midrapidity in peripheral to semicentral Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV. The centrality differential cross section ${\rm d} \sigma/ {\rm d}y$ is reported for the centrality range 40-90%, together with the doubly-differential cross section ${\rm d}^2 \sigma /{\rm d}y {\rm d} p_{\rm T}$, extracted in two peripheral centrality classes. The J/$\psi$ mesons are reconstructed in the dielectron channel, in the rapidity interval $|y| <$ 0.9 using the ALICE central barrel detectors. The J/$\psi$ cross section at midrapidity is statistically compatible to the earlier ALICE measurement at forward rapidity and at the same centre-of-mass energy, and shows only a mild centrality dependence over the covered range. Several sets of theoretical calculations taking into account the hadronic overlap in the collisions but ignoring possible final-state effects from a hot expanding medium are found to give a fairly good description of the current measurements within uncertainties.

ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015--2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons ($N_\mu>4$) and in the zenith angle range $0^{\circ}<\theta<50^{\circ}$. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about $ 4 \times 10^{15}<E_\mathrm{prim}< 6 \times 10^{16}$~eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET-II-04 is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL 2.3d and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than $20\%$ and $30\%$, respectively. The rate of high muon multiplicity events ($N_\mu>100$) obtained with QGSJET-II-04 and SIBYLL 2.3d is compatible with the data, while EPOS-LHC produces a significantly lower rate ($55\%$ of the measured rate). For both QGSJET-II-04 and SIBYLL 2.3d, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy $E_\mathrm{prim} \sim 10^{17}$ eV of these events. This result places significant constraints on more exotic production mechanisms.

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

A search for the resonant production of a heavy scalar $X$ decaying into a Higgs boson and a new lighter scalar $S$, through the process $X \to S(\to bb) H(\to \gamma\gamma)$, where the two photons are consistent with the Higgs boson decay, is performed. The search is conducted using an integrated luminosity of 140 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The search is performed over the mass range 170 $\leq$$m_{X}$$\leq$ 1000 GeV and 15 $\leq$$m_{S}$$\leq$ 500 GeV. Parameterised neural networks are used to enhance the signal purity and to achieve continuous sensitivity in a domain of the ($m_{X}$, $m_{S}$) plane. No significant excess above the expected background is found and 95% CL upper limits are set on the cross section times branching ratio, ranging from 39 fb to 0.09 fb. The largest deviation from the background-only expectation occurs for ($m_{X}$, $m_{S}$) = (575, 200) GeV with a local (global) significance of 3.5 (2.0) standard deviations.

This paper presents a search for a new $Z^\prime$ resonance decaying into a pair of dark quarks which hadronise into dark hadrons before promptly decaying back as Standard Model particles. This analysis is based on proton-proton collision data recorded at $\sqrt{s}=13$ TeV with the ATLAS detector at the Large Hadron Collider between 2015 and 2018, corresponding to an integrated luminosity of 139 fb$^{-1}$. After selecting events containing large-radius jets with high track multiplicity, the invariant mass distribution of the two highest-transverse-momentum jets is scanned to look for an excess above a data-driven estimate of the Standard Model multijet background. No significant excess of events is observed and the results are thus used to set 95 % confidence-level upper limits on the production cross-section times branching ratio of the $Z^\prime$ to dark quarks as a function of the $Z^\prime$ mass for various dark-quark scenarios.

Measurements of $W^+W^-\rightarrow e^\pm νμ^\mp ν$ production cross-sections are presented, providing a test of the predictions of perturbative quantum chromodynamics and the electroweak theory. The measurements are based on data from $pp$ collisions at $\sqrt{s}=13$ TeV recorded by the ATLAS detector at the Large Hadron Collider in 2015-2018, corresponding to an integrated luminosity of 140 fb$^{-1}$. The number of events due to top-quark pair production, the largest background, is reduced by rejecting events containing jets with $b$-hadron decays. An improved methodology for estimating the remaining top-quark background enables a precise measurement of $W^+W^-$ cross-sections with no additional requirements on jets. The fiducial $W^+W^-$ cross-section is determined in a maximum-likelihood fit with an uncertainty of 3.1%. The measurement is extrapolated to the full phase space, resulting in a total $W^+W^-$ cross-section of $127\pm4$ pb. Differential cross-sections are measured as a function of twelve observables that comprehensively describe the kinematics of $W^+W^-$ events. The measurements are compared with state-of-the-art theory calculations and excellent agreement with predictions is observed. A charge asymmetry in the lepton rapidity is observed as a function of the dilepton invariant mass, in agreement with the Standard Model expectation. A CP-odd observable is measured to be consistent with no CP violation. Limits on Standard Model effective field theory Wilson coefficients in the Warsaw basis are obtained from the differential cross-sections.

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.

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

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/ψ$ or $η_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}/Λ$ 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 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.