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.

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.

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.

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.

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.

A search for supersymmetric partners of top quarks decaying as $\tilde{t}_1\to c\tilde\chi^0_1$ and supersymmetric partners of charm quarks decaying as $\tilde{c}_1\to c\tilde\chi^0_1$, where $\tilde\chi^0_1$ is the lightest neutralino, is presented. The search uses 36.1 ${\rm fb}^{-1}$ $pp$ collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS experiment at the Large Hadron Collider and is performed in final states with jets identified as containing charm hadrons. Assuming a 100% branching ratio to $c\tilde\chi^0_1$, top and charm squarks with masses up to 850 GeV are excluded at 95% confidence level for a massless lightest neutralino. For $m_{\tilde{t}_1,\tilde{c}_1}-m_{\tilde\chi^0_1}

Multi-particle cumulants and corresponding Fourier harmonics are measured for azimuthal angle distributions of charged particles in $pp$ collisions at $\sqrt{s}$ = 5.02 and 13 TeV and in $p$+Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, and compared to the results obtained for low-multiplicity Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. These measurements aim to assess the collective nature of particle production. The measurements of multi-particle cumulants confirm the evidence for collective phenomena in $p$+Pb and low-multiplicity Pb+Pb collisions. On the other hand, the $pp$ results for four-particle cumulants do not demonstrate collective behaviour, indicating that they may be biased by contributions from non-flow correlations. A comparison of multi-particle cumulants and derived Fourier harmonics across different collision systems is presented as a function of the charged-particle multiplicity. For a given multiplicity, the measured Fourier harmonics are largest in Pb+Pb, smaller in $p$+Pb and smallest in $pp$ collisions. The $pp$ results show no dependence on the collision energy, nor on the multiplicity.

Inclusive jet and dijet cross-sections are measured in proton-proton collisions at a centre-of-mass energy of 13 TeV. The measurement uses a dataset with an integrated luminosity of 3.2 fb$^{-1}$ recorded in 2015 with the ATLAS detector at the Large Hadron Collider. Jets are identified using the anti-${k_t}$ algorithm with a radius parameter value of $R=0.4$. The inclusive jet cross-sections are measured double-differentially as a function of the jet transverse momentum, covering the range from 100 GeV to 3.5 TeV, and the absolute jet rapidity up to $|y|=3$. The double-differential dijet production cross-sections are presented as a function of the dijet mass, covering the range from 300 GeV to 9 TeV, and the half absolute rapidity separation between the two leading jets within $|y|<3$, $y*$, up to $y*=3$. Next-to-leading-order, and next-to-next-to-leading-order for the inclusive jet measurement, perturbative QCD calculations corrected for non-perturbative and electroweak effects are compared to the measured cross-sections.

A search for the dimuon decay of the Higgs boson was performed using data corresponding to an integrated luminosity of 36.1 fb$^{-1}$ collected with the ATLAS detector in $pp$ collisions at $\sqrt{s}=13$ TeV at the Large Hadron Collider. No significant excess is observed above the expected background. The observed (expected) upper limit on the cross section times branching ratio is 3.0 (3.1) times the Standard Model prediction at the 95% confidence level for a Higgs boson mass of 125 GeV. When combined with the $pp$ collision data at $\sqrt{s}=7$ TeV and $\sqrt{s}=8$ TeV, the observed (expected) upper limit is 2.8 (2.9) times the Standard Model prediction.

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.