The differential and total cross sections for kaon pair production in the pp->ppK+K- reaction have been measured at three beam energies of 2.65, 2.70, and 2.83 GeV using the ANKE magnetic spectrometer at the COSY-Juelich accelerator. These near-threshold data are separated into pairs arising from the decay of the phi-meson and the remainder. For the non-phi selection, the ratio of the differential cross sections in terms of the K-p and K+p invariant masses is strongly peaked towards low masses. This effect can be described quantitatively by using a simple ansatz for the K-p final state interaction, where it is seen that the data are sensitive to the magnitude of an effective K-p scattering length. When allowance is made for a small number of phi events where the K- rescatters from the proton, the phi region is equally well described at all three energies. A very similar phenomenon is discovered in the ratio of the cross sections as functions of the K-pp and K+pp invariant masses and the identical final state interaction model is also very successful here. The world data on the energy dependence of the non-phi total cross section is also reproduced, except possibly for the results closest to threshold.

This letter describes the observation of the light-by-light scattering process, $\gamma\gamma\rightarrow\gamma\gamma$, in Pb+Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The analysis is conducted using a data sample corresponding to an integrated luminosity of 1.73 nb$^{-1}$, collected in November 2018 by the ATLAS experiment at the LHC. Light-by-light scattering candidates are selected in events with two photons produced exclusively, each with transverse energy $E_{\textrm{T}}^{\gamma} > 3$ GeV and pseudorapidity $|\eta_{\gamma}| < 2.37$, diphoton invariant mass above 6 GeV, and small diphoton transverse momentum and acoplanarity. After applying all selection criteria, 59 candidate events are observed for a background expectation of 12 $\pm$ 3 events. The observed excess of events over the expected background has a significance of 8.2 standard deviations. The measured fiducial cross section is 78 $\pm$ 13 (stat.) $\pm$ 7 (syst.) $\pm$ 3 (lumi.) nb.

A search for the production of three massive vector bosons in proton--proton collisions is performed using data at $\sqrt{s}=13\,TeV$ recorded with the ATLAS detector at the Large Hadron Collider in the years 2015--2017, corresponding to an integrated luminosity of $79.8\,\text{fb}^{-1}$. Events with two same-sign leptons $\ell$ (electrons or muons) and at least two reconstructed jets are selected to search for $WWW\to\ell\nu\ell\nu qq$. Events with three leptons without any same-flavour opposite-sign lepton pairs are used to search for $WWW\to\ell\nu\ell\nu\ell\nu$, while events with three leptons and at least one same-flavour opposite-sign lepton pair and one or more reconstructed jets are used to search for $WWZ\to\ell\nu qq \ell\ell$. Finally, events with four leptons are analysed to search for $WWZ\to\ell\nu\ell\nu\ell\ell$ and $WZZ\to qq \ell\ell\ell\ell$. Evidence for the joint production of three massive vector bosons is observed with a significance of 4.0 standard deviations, where the expectation is 3.1 standard deviations.

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 prevalence of hadronic jets at the LHC requires that a deep understanding of jet formation and structure is achieved in order to reach the highest levels of experimental and theoretical precision. There have been many measurements of jet substructure at the LHC and previous colliders, but the targeted observables mix physical effects from various origins. Based on a recent proposal to factorize physical effects, this Letter presents a double-differential cross-section measurement of the Lund jet plane using 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collision data collected with the ATLAS detector using jets with transverse momentum above 675 GeV. The measurement uses charged particles to achieve a fine angular resolution and is corrected for acceptance and detector effects. Several parton shower Monte Carlo models are compared with the data. No single model is found to be in agreement with the measured data across the entire plane.

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.

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.

A search for magnetic monopoles and high-electric-charge objects is presented using 34.4 fb$^{-1}$ of 13 TeV $pp$ collision data collected by the ATLAS detector at the LHC during 2015 and 2016. The considered signature is based upon high ionization in the transition radiation tracker of the inner detector associated with a pencil-shape energy deposit in the electromagnetic calorimeter. The data were collected by a dedicated trigger based on the tracker high-threshold hit capability. The results are interpreted in models of Drell-Yan pair production of stable particles with two spin hypotheses (0 and 1/2) and masses ranging from 200 GeV to 4000 GeV. The search improves by approximately a factor of five the constraints on the direct production of magnetic monopoles carrying one or two Dirac magnetic charges and stable objects with electric charge in the range $20\le|z|\le60$ and extends the charge range to $60<|z|\le100$.

The dynamics of isolated-photon plus two-jet production in $pp$ collisions at a centre-of-mass energy of 13 TeV are studied with the ATLAS detector at the LHC using a dataset corresponding to an integrated luminosity of 36.1 fb$^{-1}$. Cross sections are measured as functions of a variety of observables, including angular correlations and invariant masses of the objects in the final state, $\gamma+jet+jet$. Measurements are also performed in phase-space regions enriched in each of the two underlying physical mechanisms, namely direct and fragmentation processes. The measurements cover the range of photon (jet) transverse momenta from 150 GeV (100 GeV) to 2 TeV. The tree-level plus parton-shower predictions from SHERPA and PYTHIA as well as the next-to-leading-order QCD predictions from SHERPA are compared with the measurements. The next-to-leading-order QCD predictions describe the data adequately in shape and normalisation except for regions of phase space such as those with high values of the invariant mass or rapidity separation of the two jets, where the predictions overestimate the data.

This paper presents a measurement of quantities related to the formation of jets from high-energy quarks and gluons (fragmentation). Jets with transverse momentum 100 GeV $<p_T<$ 2.5 TeV and pseudorapidity $|\eta| < 2.1$ from an integrated luminosity of 33 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collisions are reconstructed with the ATLAS detector at the Large Hadron Collider. Charged-particle tracks with $p_T > 500$ MeV and $|\eta| < 2.5$ are used to probe the detailed structure of the jet. The fragmentation properties of the more forward and the more central of the two leading jets from each event are studied. The data are unfolded to correct for detector resolution and acceptance effects. Comparisons with parton shower Monte Carlo generators indicate that existing models provide a reasonable description of the data across a wide range of phase space, but there are also significant differences. Furthermore, the data are interpreted in the context of quark- and gluon-initiated jets by exploiting the rapidity dependence of the jet flavor fraction. A first measurement of the charged-particle multiplicity using model-independent jet labels (topic modeling) provides a promising alternative to traditional quark and gluon extractions using input from simulation. The simulations provide a reasonable description of the quark-like data across the jet $p_T$ range presented in this measurement, but the gluon-like data have systematically fewer charged particles than the simulations.

This paper reports on a search for heavy resonances decaying into $WW$, $ZZ$ or $WZ$ using proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The data, corresponding to an integrated luminosity of 139 $\mathrm{fb^{-1}}$, were recorded with the ATLAS detector from 2015 to 2018 at the Large Hadron Collider. The search is performed for final states in which one $W$ or $Z$ boson decays leptonically, and the other $W$ boson or $Z$ boson decays hadronically. The data are found to be described well by expected backgrounds. Upper bounds on the production cross sections of heavy scalar, vector or tensor resonances are derived in the mass range 300-5000 GeV within the context of Standard Model extensions with warped extra dimensions or including a heavy vector triplet. Production through gluon-gluon fusion, Drell-Yan or vector-boson fusion are considered, depending on the assumed model.

A dedicated sample of Large Hadron Collider proton-proton collision data at centre-of-mass energy $\sqrt{s}=8$ TeV is used to study inclusive single diffractive dissociation, $pp \rightarrow Xp$. The intact final-state proton is reconstructed in the ATLAS ALFA forward spectrometer, while charged particles from the dissociated system $X$ are measured in the central detector components. The fiducial range of the measurement is $-4.0 < \log_{10} \xi < -1.6$ and $0.016 < |t| < 0.43 \ {\rm GeV^2}$, where $\xi$ is the proton fractional energy loss and $t$ is the squared four-momentum transfer. The total cross section integrated across the fiducial range is $1.59 \pm 0.13 \ {\rm mb}$. Cross sections are also measured differentially as functions of $\xi$, $t$, and $\Delta \eta$, a variable that characterises the rapidity gap separating the proton and the system $X$. The data are consistent with an exponential $t$ dependence, ${\rm d} \sigma / {\rm d} t \propto \text{e}^{Bt}$ with slope parameter $B = 7.65 \pm 0.34 \ {\rm GeV^{-2}}$. Interpreted in the framework of triple Regge phenomenology, the $\xi$ dependence leads to a pomeron intercept of $\alpha(0) = 1.07 \pm 0.09$.

A measurement of $W^\pm$ boson production in lead-lead collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV is reported using data recorded by the ATLAS experiment at the LHC in 2015, corresponding to a total integrated luminosity of $0.49\;\mathrm{nb^{-1}}$. The $W^\pm$ bosons are reconstructed in the electron or muon leptonic decay channels. Production yields of leptonically decaying $W^\pm$ bosons, normalised by the total number of minimum-bias events and the nuclear thickness function, are measured within a fiducial region defined by the detector acceptance and the main kinematic requirements. These normalised yields are measured separately for $W^+$ and $W^-$ bosons, and are presented as a function of the absolute value of pseudorapidity of the charged lepton and of the collision centrality. The lepton charge asymmetry is also measured as a function of the absolute value of lepton pseudorapidity. In addition, nuclear modification factors are calculated using the $W^\pm$ boson production cross-sections measured in $pp$ collisions. The results are compared with predictions based on next-to-leading-order calculations with CT14 parton distribution functions as well as with predictions obtained with the EPPS16 and nCTEQ15 nuclear parton distribution functions. No dependence of normalised production yields on centrality and a good agreement with predictions are observed for mid-central and central collisions. For peripheral collisions, the data agree with predictions within 1.7 (0.9) standard deviations for $W^-$ ($W^+$) bosons.

This Letter presents direct searches for lepton flavour violation in Higgs boson decays, $H\rightarrow e\tau$ and $H\rightarrow\mu\tau$, performed with the ATLAS detector at the LHC. The searches are based on a data sample of proton-proton collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV, corresponding to an integrated luminosity of $36.1\,\mathrm{fb}^{-1}$. No significant excess is observed above the expected background from Standard Model processes. The observed (median expected) 95 % confidence-level upper limits on the lepton-flavour-violating branching ratios are $0.47\%$ ($0.34^{+0.13}_{-0.10}\,\%$) and $0.28\%$ ($0.37^{+0.14}_{-0.10}\,\%$) for $H\to e\tau$ and $H\to\mu\tau$, respectively.

The azimuthal anisotropy of charged particles produced in $\sqrt{s_{\mathrm{NN}}}=8.16$ TeV $p$+Pb collisions is measured with the ATLAS detector at the LHC. The data correspond to an integrated luminosity of $165$ $\mathrm{nb}^{-1}$ that was collected in 2016. Azimuthal anisotropy coefficients, elliptic $v_2$ and triangular $v_3$, extracted using two-particle correlations with a non-flow template fit procedure, are presented as a function of particle transverse momentum ($p_\mathrm{T}$) between $0.5$ and $50$ GeV. The $v_2$ results are also reported as a function of centrality in three different particle $p_\mathrm{T}$ intervals. The results are reported from minimum-bias events and jet-triggered events, where two jet $p_\mathrm{T}$ thresholds are used. The anisotropies for particles with $p_\mathrm{T}$ less than about $2$ GeV are consistent with hydrodynamic flow expectations, while the significant non-zero anisotropies for $p_\mathrm{T}$ in the range $9$-$50$ GeV are not explained within current theoretical frameworks. In the $p_\mathrm{T}$ range $2$-$9$ GeV, the anisotropies are larger in minimum-bias than in jet-triggered events. Possible origins of these effects, such as the changing admixture of particles from hard scattering and the underlying event, are discussed.

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.

Differential cross-section measurements are presented for the electroweak production of two jets in association with a $Z$ boson. These measurements are sensitive to the vector-boson fusion production mechanism and provide a fundamental test of the gauge structure of the Standard Model. The analysis is performed using proton-proton collision data collected by ATLAS at $\sqrt{s}$=13 TeV and with an integrated luminosity of 139 fb$^{-1}$. The differential cross-sections are measured in the $Z\rightarrow \ell^+\ell^-$ decay channel ($\ell=e,\mu$) as a function of four observables: the dijet invariant mass, the rapidity interval spanned by the two jets, the signed azimuthal angle between the two jets, and the transverse momentum of the dilepton pair. The data are corrected for the effects of detector inefficiency and resolution and are sufficiently precise to distinguish between different state-of-the-art theoretical predictions calculated using Powheg+Pythia8, Herwig7+Vbfnlo and Sherpa 2.2. The differential cross-sections are used to search for anomalous weak-boson self-interactions using a dimension-six effective field theory. The differential cross-section as a function of the signed azimuthal angle between the two jets is found to be particularly sensitive to the interference between the Standard Model and dimension-six scattering amplitudes and provides a direct test of charge-conjugation and parity invariance in the weak-boson self-interactions.

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.

This paper presents results of searches for electroweak production of supersymmetric particles in models with compressed mass spectra. The searches use 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collision data collected by the ATLAS experiment at the Large Hadron Collider. Events with missing transverse momentum and two same-flavor, oppositely charged, low transverse momentum leptons are selected, and are further categorized by the presence of hadronic activity from initial-state radiation or a topology compatible with vector-boson fusion processes. The data are found to be consistent with predictions from the Standard Model. The results are interpreted using simplified models of $R$-parity-conserving supersymmetry in which the lightest supersymmetric partner is a neutralino with a mass similar to the lightest chargino, the second-to-lightest neutralino or the slepton. Lower limits on the masses of charginos in different simplified models range from 193 GeV to 240 GeV for moderate mass splittings, and extend down to mass splittings of 1.5 GeV to 2.4 GeV at the LEP chargino bounds (92.4 GeV). Similar lower limits on degenerate light-flavor sleptons extend up to masses of 251 GeV and down to mass splittings of 550 MeV. Constraints on vector-boson fusion production of electroweak SUSY states are also presented.

A search for dark matter is conducted in final states containing a photon and missing transverse momentum in proton$-$proton collisions at $\sqrt{s}$ = 13 TeV. The data, collected during 2015$-$2018 by the ATLAS experiment at the CERN LHC, correspond to an integrated luminosity of 139 fb$^{-1}$. No deviations from the predictions of the Standard Model are observed and 95% confidence-level upper limits between 2.45 fb and 0.5 fb are set on the visible cross section for contributions from physics beyond the Standard Model, in different ranges of the missing transverse momentum. The results are interpreted as 95% confidence-level limits in models where weakly interacting dark-matter candidates are pair-produced via an s-channel axial-vector or vector mediator. Dark-matter candidates with masses up to 415 (580) GeV are excluded for axial-vector (vector) mediators, while the maximum excluded mass of the mediator is 1460 (1470) GeV. In addition, the results are expressed in terms of 95% confidence-level limits on the parameters of a model with an axion-like particle produced in association with a photon, and are used to constrain the coupling $g_{aZ\gamma}$ of an axion-like particle to the electroweak gauge bosons.

Fiducial and differential measurements of $W^+W^-$ production in events with at least one hadronic jet are presented. These cross-section measurements are sensitive to the properties of electroweak-boson self-interactions and provide a test of perturbative quantum chromodynamics and the electroweak theory. The analysis is performed using proton$-$proton collision data collected at $\sqrt{s}=13~$TeV with the ATLAS experiment, corresponding to an integrated luminosity of 139$~$fb$^{-1}$. Events are selected with exactly one oppositely charged electron$-$muon pair and at least one hadronic jet with a transverse momentum of $p_{\mathrm{T}}>30~$GeV and a pseudorapidity of $|\eta|<4.5$. After subtracting the background contributions and correcting for detector effects, the jet-inclusive $W^+W^-+\ge 1~$jet fiducial cross-section and $W^+W^-+$ jets differential cross-sections with respect to several kinematic variables are measured, thus probing a previously unexplored event topology at the LHC. These measurements include leptonic quantities, such as the lepton transverse momenta and the transverse mass of the $W^+W^-$ system, as well as jet-related observables such as the leading jet transverse momentum and the jet multiplicity. Limits on anomalous triple-gauge-boson couplings are obtained in a phase space where interference between the Standard Model amplitude and the anomalous amplitude is enhanced.

The differential cross section for isolated-photon production in $pp$ collisions is measured at a centre-of-mass energy of 13 TeV with the ATLAS detector at the LHC using an integrated luminosity of 36.1 fb$^{-1}$. The differential cross section is presented as a function of the photon transverse energy in different regions of photon pseudorapidity. The differential cross section as a function of the absolute value of the photon pseudorapidity is also presented in different regions of photon transverse energy. Next-to-leading-order QCD calculations from JETPHOX and SHERPA as well as next-to-next-to-leading-order QCD calculations from NNLOJET are compared with the measurement, using several parameterisations of the proton parton distribution functions. The predictions provide a good description of the data within the experimental and theoretical uncertainties.

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.

A search for the direct production of the supersymmetric partners of $\tau$-leptons (staus) in final states with two hadronically decaying $\tau$-leptons is presented. The analysis uses a dataset of $pp$ collisions corresponding to an integrated luminosity of $139$ fb$^{-1}$, recorded with the ATLAS detector at the Large Hadron Collider at a center-of-mass energy of 13 TeV. No significant deviation from the expected Standard Model background is observed. Limits are derived in scenarios of direct production of stau pairs with each stau decaying into the stable lightest neutralino and one $\tau$-lepton in simplified models where the two stau mass eigenstates are degenerate. Stau masses from 120 GeV to 390 GeV are excluded at 95% confidence level for a massless lightest neutralino.

A search for the pair production of heavy leptons as predicted by the type-III seesaw mechanism is presented. The search uses proton-proton collision data at a centre-of-mass energy of 13 TeV, corresponding to 139 fb$^{-1}$ of integrated luminosity recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. The analysis focuses on the final state with two light leptons (electrons or muons) of different flavour and charge combinations, with at least two jets and large missing transverse momentum. No significant excess over the Standard Model expectation is observed. The results are translated into exclusion limits on heavy-lepton masses, and the observed lower limit on the mass of the type-III seesaw heavy leptons is 790 GeV at 95% confidence level.

This paper presents a search for new heavy particles decaying into a pair of top quarks using 139 fb$^{-1}$ of proton--proton collision data recorded at a centre-of-mass energy of $\sqrt{s}=13$ TeV with the ATLAS detector at the Large Hadron Collider. The search is performed using events consistent with pair production of high-transverse-momentum top quarks and their subsequent decays into the fully hadronic final states. The analysis is optimized for resonances decaying into a $t\bar{t}$ pair with mass above 1.4 TeV, exploiting a dedicated multivariate technique with jet substructure to identify hadronically decaying top quarks using large-radius jets and evaluating the background expectation from data. No significant deviation from the background prediction is observed. Limits are set on the production cross-section times branching fraction for the new $Z'$ boson in a topcolor-assisted-technicolor model. The $Z'$ boson masses below 3.9 and 4.7 TeV are excluded at 95% confidence level for the decay widths of 1% and 3%, respectively.

The elliptic flow of muons from the decay of charm and bottom hadrons is measured in $pp$ collisions at $\sqrt{s}=13$ TeV using a data sample with an integrated luminosity of 150 pb$^{-1}$ recorded by the ATLAS detector at the LHC. The muons from heavy-flavor decay are separated from light-hadron decay muons using momentum imbalance between the tracking and muon spectrometers. The heavy-flavor decay muons are further separated into those from charm decay and those from bottom decay using the distance-of-closest-approach to the collision vertex. The measurement is performed for muons in the transverse momentum range 4-7 GeV and pseudorapidity range $|\eta|<2.4$. A significant non-zero elliptic anisotropy coefficient $v_{2}$ is observed for muons from charm decays, while the $v_{2}$ value for muons from bottom decays is consistent with zero within uncertainties.

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.

Narrow resonances decaying into $WW$, $WZ$ or $ZZ$ boson pairs are searched for in 139 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV recorded with the ATLAS detector at the Large Hadron Collider from 2015 to 2018. The diboson system is reconstructed using pairs of high transverse momentum, large-radius jets. These jets are built from a combination of calorimeter- and tracker-inputs compatible with the hadronic decay of a boosted $W$ or $Z$ boson, using jet mass and substructure properties. The search is performed for diboson resonances with masses greater than 1.3 TeV. No significant deviations from the background expectations are observed. Exclusion limits at the 95% confidence level are set on the production cross-section times branching ratio into dibosons for resonances in a range of theories beyond the Standard Model, with the highest excluded mass of a new gauge boson at 3.8 TeV in the context of mass-degenerate resonances that couple predominantly to gauge bosons.

This paper describes precision measurements of the transverse momentum $p_\mathrm{T}^{\ell\ell}$ ($\ell=e,\mu$) and of the angular variable $\phi^{*}_{\eta}$ distributions of Drell-Yan lepton pairs in a mass range of 66-116 GeV. The analysis uses data from 36.1 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV collected by the ATLAS experiment at the LHC in 2015 and 2016. Measurements in electron-pair and muon-pair final states are performed in the same fiducial volumes, corrected for detector effects, and combined. Compared to previous measurements in proton-proton collisions at $\sqrt{s}=$7 and 8 TeV, these new measurements probe perturbative QCD at a higher centre-of-mass energy with a different composition of initial states. They reach a precision of 0.2% for the normalized spectra at low values of $p_\mathrm{T}^{\ell\ell}$. The data are compared with different QCD predictions, where it is found that predictions based on resummation approaches can describe the full spectrum within uncertainties.

A search for new-physics resonances decaying into a lepton and a jet performed by the ATLAS experiment is presented. Scalar leptoquarks pair-produced in $pp$ collisions at $\sqrt{s}=13$ TeV at the Large Hadron Collider are considered using an integrated luminosity of 139 fb$^{-1}$, corresponding to the full Run 2 dataset. They are searched for in events with two electrons or two muons and two or more jets, including jets identified as arising from the fragmentation of $c$- or $b$-quarks. The observed yield in each channel is consistent with the Standard Model background expectation. Leptoquarks with masses below 1.8 TeV and 1.7 TeV are excluded in the electron and muon channels, respectively, assuming a branching ratio into a charged lepton and a quark of 100%, with minimal dependence on the quark flavour. Upper limits on the aforementioned branching ratio are also given as a function of the leptoquark mass.

Measurements of joint-polarisation states of $W$ and $Z$ gauge bosons in $W^{\pm}Z$ production are presented. The data set used corresponds to an integrated luminosity of $139$ fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of $13$ TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. The $W^{\pm}Z$ candidate events are reconstructed using leptonic decay modes of the gauge bosons into electrons and muons. The simultaneous pair-production of longitudinally polarised vector bosons is measured for the first time with a significance of $7.1$ standard deviations. The measured joint helicity fractions integrated over the fiducial region are $f_{\mathrm{00}} = 0.067 \pm 0.010$, $f_{\mathrm{0T}} = 0.110 \pm 0.029$, $f_{\mathrm{T0}} = 0.179 \pm 0.023$ and $f_{\mathrm{TT}} = 0.644 \pm 0.032$, in agreement with the next-to-leading-order Standard Model predictions. Individual helicity fractions of the $W$ and $Z$ bosons are also measured and found to be consistent with joint helicity fractions within the expected amount of correlations. Both the joint and individual helicity fractions are also measured separately in $W^+Z$ and $W^-Z$ events. Inclusive and differential cross sections for several kinematic observables sensitive to polarisation are presented.

A search for quantum black holes in electron+jet and muon+jet invariant mass spectra is performed with 140 fb$^{-1}$ of data collected by the ATLAS detector in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider. The observed invariant mass spectrum of lepton+jet pairs is consistent with Standard Model expectations. Upper limits are set at 95% confidence level on the production cross-sections times branching fractions for quantum black holes decaying into a lepton and a quark in a search region with invariant mass above 2.0 TeV. The resulting quantum black hole lower mass threshold limit is 9.2 TeV in the Arkani-Hamed-Dimopoulos-Dvali model, and 6.8 TeV in the Randall-Sundrum model.

This paper reports cross-section measurements of $ZZ$ production in $pp$ collisions at $\sqrt{s}=13.6$ TeV at the Large Hadron Collider. The data were collected by the ATLAS detector in 2022, and correspond to an integrated luminosity of 29 fb$^-1$. Events in the $ZZ\rightarrow4\ell$ ($\ell = e$, $\mu$) final states are selected and used to measure the inclusive and differential cross-sections in a fiducial region defined close to the analysis selections. The inclusive cross-section is further extrapolated to the total phase space with a requirement of 66 $< m_Z <$ 116 GeV for both $Z$ bosons, yielding $16.8 \pm 1.1$ pb. The results are well described by the Standard Model predictions.

Higgsinos with masses near the electroweak scale can solve the hierarchy problem and provide a dark matter candidate, while detecting them at the LHC remains challenging if their mass-splitting is $\mathcal{O}$(1 GeV). This Letter presents a novel search for nearly mass-degenerate higgsinos in events with an energetic jet, missing transverse momentum, and a low-momentum track with a significant transverse impact parameter using 140 fb$^{-1}$ of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment. For the first time since LEP, a range of mass-splittings between the lightest charged and neutral higgsinos from 0.3 GeV to 0.9 GeV is excluded at 95% confidence level, with a maximum reach of approximately 170 GeV in the higgsino mass.

This letter presents a search for narrow, high-mass resonances in the $Z\gamma$ final state with the $Z$ boson decaying into a pair of electrons or muons. The $\sqrt{s}=13$ TeV $pp$ collision data were recorded by the ATLAS detector at the CERN Large Hadron Collider and have an integrated luminosity of 140 fb$^{-1}$. The data are found to be in agreement with the Standard Model background expectation. Upper limits are set on the resonance production cross section times the decay branching ratio into $Z\gamma$. For spin-0 resonances produced via gluon-gluon fusion, the observed limits at 95% confidence level vary between 65.5 fb and 0.6 fb, while for spin-2 resonances produced via gluon-gluon fusion (or quark-antiquark initial states) limits vary between 77.4 (76.1) fb and 0.6 (0.5) fb, for the mass range from 220 GeV to 3400 GeV.

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 paper presents measurements of charged-particle distributions sensitive to the properties of the underlying event in events containing a $Z$ boson decaying into a muon pair. The data were obtained using the ATLAS detector at the LHC in proton-proton collisions at a centre-of-mass energy of 13 TeV with an integrated luminosity of 3.2 fb$^{-1}$. Distributions of the charged-particle multiplicity and of the charged-particle transverse momentum are measured in regions of the azimuth defined relative to the $Z$ boson direction. The measured distributions are compared with the predictions of various Monte Carlo generators which implement different underlying-event models. The Monte Carlo model predictions qualitatively describe the data well, but with some significant discrepancies.

Measurements of the Standard Model Higgs boson decaying into a $b\bar{b}$ pair and produced in association with a $W$ or $Z$ boson decaying into leptons, using proton-proton collision data collected between 2015 and 2018 by the ATLAS detector, are presented. The measurements use collisions produced by the Large Hadron Collider at a centre-of-mass energy of $\sqrt{s} = $13 TeV, corresponding to an integrated luminosity of 139 fb$^{-1}$. The production of a Higgs boson in association with a $W$ or $Z$ boson is established with observed (expected) significances of 4.0 (4.1) and 5.3 (5.1) standard deviations, respectively. Cross-sections of associated production of a Higgs boson decaying into bottom quark pairs with an electroweak gauge boson, $W$ or $Z$, decaying into leptons are measured as a function of the gauge boson transverse momentum in kinematic fiducial volumes. The cross-section measurements are all consistent with the Standard Model expectations, and the total uncertainties vary from 30% in the high gauge boson transverse momentum regions to 85% in the low regions. Limits are subsequently set on the parameters of an effective Lagrangian sensitive to modifications of the $WH$ and $ZH$ processes as well as the Higgs boson decay into $b\bar{b}$.

The problems of neutrino masses, matter-antimatter asymmetry, and dark matter could be successfully addressed by postulating right-handed neutrinos with Majorana masses below the electroweak scale. In this work, leptonic decays of $W$ bosons extracted from 32.9 fb$^{-1}$ to 36.1 fb$^{-1}$ of 13 TeV proton-proton collisions at the LHC are used to search for heavy neutral leptons (HNLs) that are produced through mixing with muon or electron neutrinos. The search is conducted using the ATLAS detector in both prompt and displaced leptonic decay signatures. The prompt signature requires three leptons produced at the interaction point (either $\mu\mu e$ or $e e\mu$) with a veto on same-flavour opposite-charge topologies. The displaced signature comprises a prompt muon from the $W$ boson decay and the requirement of a dilepton vertex (either $\mu\mu$ or $\mu e$) displaced in the transverse plane by 4-300 mm from the interaction point. The search sets constraints on the HNL mixing to muon and electron neutrinos for HNL masses in the range 4.5-50 GeV.

The production of single top quarks and top antiquarks via the $t$-channel exchange of a virtual $W$ boson is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV at the LHC using $140\,\mathrm{fb^{-1}}$ of ATLAS data. The total cross-sections are determined to be $\sigma(tq)=137^{+8}_{-8}\,\mathrm{pb}$ and $\sigma(\bar{t}q)=84^{+6}_{-5}\,\mathrm{pb}$ for top-quark and top-antiquark production, respectively. The combined cross-section is found to be $\sigma(tq+\bar{t}q)=221^{+13}_{-13}\,\mathrm{pb}$ and the cross-section ratio is $R_{t}=\sigma(tq)/\sigma(\bar{t}q)=1.636^{+0.036}_{-0.034}$. The predictions at next-to-next-to-leading-order in quantum chromodynamics are in good agreement with these measurements. The predicted value of $R_{t}$ using different sets of parton distribution functions is compared with the measured value, demonstrating the potential to further constrain the functions when using this result in global fits. The measured cross-sections are interpreted in an effective field theory approach, setting limits at the 95% confidence level on the strength of a four-quark operator and an operator coupling the third quark generation to the Higgs boson doublet: $-0.37 < C_{Qq}^{3,1}/\Lambda^2 < 0.06$ and $-0.87 < C_{\phi Q}^{3}/\Lambda^2 < 1.42$. The constraint $|V_{tb}|>0.95$ at the 95% confidence level is derived from the measured value of $\sigma(tq+\bar{t}q)$. In a more general approach, pairs of CKM matrix elements involving top quarks are simultaneously constrained, leading to confidence contours in the corresponding two-dimensional parameter spaces.

A search for new resonances decaying into a pair of jets is reported using the dataset of proton-proton collisions 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}$. The distribution of the invariant mass of the two leading jets is examined for local excesses above a data-derived estimate of the Standard Model background. In addition to an inclusive dijet search, events with jets identified as containing $b$-hadrons are examined specifically. No significant excess of events above the smoothly falling background spectra is observed. The results are used to set cross-section upper limits at 95% confidence level on a range of new physics scenarios. Model-independent limits on Gaussian-shaped signals are also reported. The analysis looking at jets containing $b$-hadrons benefits from improvements in the jet flavour identification at high transverse momentum, which increases its sensitivity relative to the previous analysis beyond that expected from the higher integrated luminosity.

A search for long-lived particles decaying into an oppositely charged lepton pair, $\mu\mu$, $ee$, or $e\mu$, is presented using 32.8 fb$^{-1}$ of $pp$ collision data collected at $\sqrt{s}=13$ TeV by the ATLAS detector at the LHC. Candidate leptons are required to form a vertex, within the inner tracking volume of ATLAS, displaced from the primary $pp$ interaction region. No lepton pairs with an invariant mass greater than 12 GeV are observed, consistent with the background expectations derived from data. The detection efficiencies for generic resonances with lifetimes ($c\tau$) of 100-1000 mm decaying into a dilepton pair with masses between 0.1-1.0 TeV are presented as a function of $p_T$ and decay radius of the resonances to allow the extraction of upper limits on the cross sections for theoretical models. The result is also interpreted in a supersymmetric model in which the lightest neutralino, produced via squark-antisquark production, decays into $\ell^{+}\ell^{'-}\nu$ ($\ell, \ell^{'} = e$, $\mu$) with a finite lifetime due to the presence of R-parity violating couplings. Cross-section limits are presented for specific squark and neutralino masses. For a 700 GeV squark, neutralinos with masses of 50-500 GeV and mean proper lifetimes corresponding to $c\tau$ values between 1 mm to 6 m are excluded. For a 1.6 TeV squark, $c\tau$ values between 3 mm to 1 m are excluded for 1.3 TeV neutralinos.

Results of a search for new physics in final states with an energetic jet and large missing transverse momentum are reported. The search uses proton-proton collision data corresponding to an integrated luminosity of 139 fb$^{-1}$ at a center-of-mass energy of 13 TeV collected in the period 2015-2018 with the ATLAS detector at the Large Hadron Collider. Compared to previous publications, in addition to an increase of almost a factor of four in the data size, the analysis implements a number of improvements in the signal selection and the background determination leading to enhanced sensitivity. Events are required to have at least one jet with transverse momentum above 150 GeV and no reconstructed leptons ($e$, $\mu$ or $\tau$) or photons. Several signal regions are considered with increasing requirements on the missing transverse momentum starting at 200 GeV. Overall agreement is observed between the number of events in data and the Standard Model predictions. Model-independent 95 % confidence-level limits on visible cross sections for new processes are obtained in the range between 736 fb and 0.3 fb with increasing missing transverse momentum. Results are also translated into improved exclusion limits in models with pair-produced weakly interacting dark-matter candidates, large extra spatial dimensions, supersymmetric particles in several compressed scenarios, axion-like particles, and new scalar particles in dark-energy-inspired models. In addition, the data are translated into bounds on the invisible branching ratio of the Higgs boson.

A search for the Higgs boson decaying into a photon and a pair of electrons or muons with an invariant mass $m_{\ell\ell} < 30$ GeV is presented. The analysis is performed using 139 fb$^{-1}$ of proton-proton collision data, produced by the LHC at a centre-of-mass energy of 13 TeV and collected by the ATLAS experiment. Evidence for the $H \rightarrow \ell \ell \gamma$ process is found with a significance of 3.2$\sigma$ over the background-only hypothesis, compared to an expected significance of 2.1$\sigma$. The best-fit value of the signal strength parameter, defined as the ratio of the observed signal yield to the one expected in the Standard Model, is $\mu = 1.5 \pm 0.5$. The Higgs boson production cross-section times the $H \rightarrow\ell\ell\gamma$ branching ratio for $m_{\ell\ell} <$ 30 GeV is determined to be 8.7 $^{+2.8}_{-2.7}$ fb.

Measurements of $K_S^0$ and $\Lambda^0$ production in $t\bar{t}$ final states have been performed. They are based on a data sample with integrated luminosity of 4.6 $\mathrm{fb}^{-1}$ from proton-proton collisions at a centre-of-mass energy of 7 TeV, collected in 2011 with the ATLAS detector at the Large Hadron Collider. Neutral strange particles are separated into three classes, depending on whether they are contained in a jet, with or without a $b$-tag, or not associated with a selected jet. The aim is to look for differences in their main kinematic distributions. A comparison of data with several Monte Carlo simulations using different hadronisation and fragmentation schemes, colour reconnection models and different tunes for the underlying event has been made. The production of neutral strange particles in $t\bar{t}$ dileptonic events is found to be well described by current Monte Carlo models for $K_S^0$ and $\Lambda^0$ production within jets, but not for those produced outside jets.

To assess the properties of the quark-gluon plasma formed in heavy-ion collisions, the ATLAS experiment at the LHC measures a correlation between the mean transverse momentum and the magnitudes of the flow harmonics. The analysis uses data samples of lead-lead and proton-lead collisions obtained at the centre-of-mass energy per nucleon pair of 5.02 TeV, corresponding to total integrated luminosities of $22 ~\mu b^{-1}$ and $28~nb^{-1}$, respectively. The measurement is performed using a modified Pearson correlation coefficient with the charged-particle tracks on an event-by-event basis. The modified Pearson correlation coefficients for the $2^{nd}$-, 3$^{rd}$-, and 4$^{th}$-order harmonics are measured as a function of event centrality quantified as the number of charged particles or the number of nucleons participating in the collision. The measurements are performed for several intervals of the charged-particle transverse momentum. The correlation coefficients for all studied harmonics exhibit a strong centrality evolution in the lead-lead collisions, which only weakly depends on the charged-particle momentum range. In the proton-lead collisions, the modified Pearson correlation coefficient measured for the second harmonics shows only weak centrality dependence. The data is qualitatively described by the predictions based on the hydrodynamical model.

This Letter presents a search for the production of new heavy resonances decaying into a Higgs boson and a photon using proton-proton collision data at $\sqrt{s}=13$ TeV collected by the ATLAS detector at the LHC. The data correspond to an integrated luminosity of 139 fb$^{-1}$. The analysis is performed by reconstructing hadronically decaying Higgs boson $(H\to b\bar{b})$ candidates as single large-radius jets. A novel algorithm using information about the jet constituents in the center-of-mass frame of the jet is implemented to identify the two $b$-quarks in the single jet. No significant excess of events is observed above the expected background. Upper limits are set on the production cross-section times branching fraction for narrow spin-1 resonances decaying into a Higgs boson and a photon in the resonance mass range from 0.7 to 4 TeV, cross-sections times branching fraction are excluded between 11.6 fb and 0.11 fb at a 95% confidence level.

Azimuthal anisotropies of muons from charm and bottom hadron decays are measured in Pb+Pb collisions at $\sqrt{s_\mathrm{NN}}= 5.02$ TeV. The data were collected with the ATLAS detector at the Large Hadron Collider in 2015 and 2018 with integrated luminosities of $0.5~\mathrm{nb}^{-1}$ and $1.4~\mathrm{nb^{-1}}$, respectively. The kinematic selection for heavy-flavor muons requires transverse momentum $4 < p_\mathrm{T} < 30$ GeV and pseudorapidity $|\eta|<2.0$. The dominant sources of muons in this $p_\mathrm{T}$ range are semi-leptonic decays of charm and bottom hadrons. These heavy-flavor muons are separated from light-hadron decay muons and punch-through hadrons using the momentum imbalance between the measurements in the tracking detector and in the muon spectrometers. Azimuthal anisotropies, quantified by flow coefficients, are measured via the event-plane method for inclusive heavy-flavor muons as a function of the muon $p_\mathrm{T}$ and in intervals of Pb+Pb collision centrality. Heavy-flavor muons are separated into contributions from charm and bottom hadron decays using the muon transverse impact parameter with respect to the event primary vertex. Non-zero elliptic ($v_{2}$) and triangular ($v_{3}$) flow coefficients are extracted for charm and bottom muons, with the charm muon coefficients larger than those for bottom muons for all Pb+Pb collision centralities. The results indicate substantial modification to the charm and bottom quark angular distributions through interactions in the quark-gluon plasma produced in these Pb+Pb collisions, with smaller modifications for the bottom quarks as expected theoretically due to their larger mass.