This paper presents studies of Bose-Einstein correlations (BEC) in proton-proton collisions at a centre-of-mass energy of 13 TeV, using data from the ATLAS detector at the CERN Large Hadron Collider. Data were collected in a special low-luminosity configuration with a minimum-bias trigger and a high-multiplicity track trigger, accumulating integrated luminosities of 151 $\mu$b$^{-1}$ and 8.4 nb$^{-1}$ respectively. The BEC are measured for pairs of like-sign charged particles, each with $|\eta|$ < 2.5, for two kinematic ranges: the first with particle $p_T$ > 100 MeV and the second with particle $p_T$ > 500 MeV. The BEC parameters, characterizing the source radius and particle correlation strength, are investigated as functions of charged-particle multiplicity (up to 300) and average transverse momentum of the pair (up to 1.5 GeV). The double-differential dependence on charged-particle multiplicity and average transverse momentum of the pair is also studied. The BEC radius is found to be independent of the charged-particle multiplicity for high charged-particle multiplicity (above 100), confirming a previous observation at lower energy. This saturation occurs independent of the transverse momentum of the pair.

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.

The exclusive production of pion pairs in the process $pp\to pp\pi^+\pi^-$ has been measured at $\sqrt{s}$ = 7 TeV with the ATLAS detector at the LHC, using 80 $\mu$b$^{-1}$ of low-luminosity data. The pion pairs were detected in the ATLAS central detector while outgoing protons were measured in the forward ATLAS ALFA detector system. This represents the first use of proton tagging to measure an exclusive hadronic final state at the LHC. A cross-section measurement is performed in two kinematic regions defined by the proton momenta, the pion rapidities and transverse momenta, and the pion-pion invariant mass. Cross section values of $4.8 \pm 1.0 \text{(stat.)} + {}^{+0.3}_{-0.2} \text{(syst.)}\mu$b and $9 \pm 6 \text{(stat.)} + {}^{+2}_{-2}\text{(syst.)}\mu$b are obtained in the two regions; they are compared with theoretical models and provide a demonstration of the feasibility of measurements of this type.

In a special run of the LHC with $\beta^\star = 2.5~$km, proton-proton elastic-scattering events were recorded at $\sqrt{s} = 13~$TeV with an integrated luminosity of $340~\mu \textrm{b}^{-1}$ using the ALFA subdetector of ATLAS in 2016. The elastic cross section was measured differentially in the Mandelstam $t$ variable in the range from $-t = 2.5 \cdot 10^{-4}~$GeV$^{2}$ to $-t = 0.46~$GeV$^{2}$ using 6.9 million elastic-scattering candidates. This paper presents measurements of the total cross section $\sigma_{\textrm{tot}}$, parameters of the nuclear slope, and the $\rho$-parameter defined as the ratio of the real part to the imaginary part of the elastic-scattering amplitude in the limit $t \rightarrow 0$. These parameters are determined from a fit to the differential elastic cross section using the optical theorem and different parameterizations of the $t$-dependence. The results for $\sigma_{\textrm{tot}}$ and $\rho$ are \begin{equation*} \sigma_{\textrm{tot}}(pp\rightarrow X) = \mbox{104.7} \pm 1.1 \; \mbox{mb} , \; \; \; \rho = \mbox{0.098} \pm 0.011 . \end{equation*} The uncertainty in $\sigma_{\textrm{tot}}$ is dominated by the luminosity measurement, and in $\rho$ by imperfect knowledge of the detector alignment and by modelling of the nuclear amplitude.

Measurements of the differential production cross-sections of prompt and non-prompt $J/\psi$ and $\psi(2$S$)$ mesons with transverse momenta between 8 and 360 GeV and rapidity in the range $|y|<2$ are reported. Furthermore, measurements of the non-prompt fractions of $J/\psi$ and $\psi(2$S$)$, and the prompt and non-prompt $\psi(2$S$)$-to-$J/\psi$ production ratios, are presented. The analysis is performed using 140 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data recorded by the ATLAS detector at the LHC during the years 2015-2018.

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.

The associated production of a Higgs boson and a top-quark pair is measured in events characterised by the presence of one or two electrons or muons. The Higgs boson decay into a $b$-quark pair is used. The analysed data, corresponding to an integrated luminosity of 139 fb$^{-1}$, were collected in proton-proton collisions at the Large Hadron Collider between 2015 and 2018 at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The measured signal strength, defined as the ratio of the measured signal yield to that predicted by the Standard Model, is $0.35^{+0.36}_{-0.34}$. This result is compatible with the Standard Model prediction and corresponds to an observed (expected) significance of 1.0 (2.7) standard deviations. The signal strength is also measured differentially in bins of the Higgs boson transverse momentum in the simplified template cross-section framework, including a bin for specially selected boosted Higgs bosons with transverse momentum above 300 GeV.

A search for heavy neutral Higgs bosons is performed using the LHC Run 2 data, corresponding to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}=13$ TeV recorded with the ATLAS detector. The search for heavy resonances is performed over the mass range 0.2-2.5 TeV for the $\tau^+\tau^-$ decay with at least one $\tau$-lepton decaying into final states with hadrons. The data are in good agreement with the background prediction of the Standard Model. In the $M_{h}^{125}$ scenario of the Minimal Supersymmetric Standard Model, values of $\tan\beta>8$ and $\tan\beta>21$ are excluded at the 95% confidence level for neutral Higgs boson masses of 1.0 TeV and 1.5 TeV, respectively, where $\tan\beta$ is the ratio of the vacuum expectation values of the two Higgs doublets.

The production cross-sections for $W^{\pm}$ and $Z$ bosons are measured using ATLAS data corresponding to an integrated luminosity of 4.0 pb$^{-1}$ collected at a centre-of-mass energy $\sqrt{s}=2.76$ TeV. The decay channels $W \rightarrow \ell \nu$ and $Z \rightarrow \ell \ell $ are used, where $\ell$ can be an electron or a muon. The cross-sections are presented for a fiducial region defined by the detector acceptance and are also extrapolated to the full phase space for the total inclusive production cross-section. The combined (average) total inclusive cross-sections for the electron and muon channels are: \begin{eqnarray} \sigma^{\text{tot}}_{W^{+}\rightarrow \ell \nu}& = & 2312 \pm 26\ (\text{stat.})\ \pm 27\ (\text{syst.}) \pm 72\ (\text{lumi.}) \pm 30\ (\text{extr.})\text{pb} \nonumber, \\ \sigma^{\text{tot}}_{W^{-}\rightarrow \ell \nu}& = & 1399 \pm 21\ (\text{stat.})\ \pm 17\ (\text{syst.}) \pm 43\ (\text{lumi.}) \pm 21\ (\text{extr.})\text{pb} \nonumber, \\ \sigma^{\text{tot}}_{Z \rightarrow \ell \ell}& = & 323.4 \pm 9.8\ (\text{stat.}) \pm 5.0\ (\text{syst.}) \pm 10.0\ (\text{lumi.}) \pm 5.5 (\text{extr.}) \text{pb} \nonumber. \end{eqnarray} Measured ratios and asymmetries constructed using these cross-sections are also presented. These observables benefit from full or partial cancellation of many systematic uncertainties that are correlated between the different measurements.

The production of a prompt photon in association with a $Z$ boson is studied in proton-proton collisions at a centre-of-mass energy $\sqrt{s} =$ 13 TeV. The analysis uses a data sample with an integrated luminosity of 139 fb$^{-1}$ collected by the ATLAS detector at the LHC from 2015 to 2018. The production cross-section for the process $pp \rightarrow \ell^+\ell^-\gamma+X$ ($\ell = e, \mu$) is measured within a fiducial phase-space region defined by kinematic requirements on the photon and the leptons, and by isolation requirements on the photon. An experimental precision of 2.9% is achieved for the fiducial cross-section. Differential cross-sections are measured as a function of each of six kinematic variables characterising the $\ell^+\ell^-\gamma$ system. The data are compared with theoretical predictions based on next-to-leading-order and next-to-next-to-leading-order perturbative QCD calculations. The impact of next-to-leading-order electroweak corrections is also considered.