Jet quenching is the process of color-charged partons losing energy via interactions with quark-gluon plasma droplets created in heavy-ion collisions. The collective expansion of such droplets is well described by viscous hydrodynamics. Similar evidence of collectivity is consistently observed in smaller collision systems, including $pp$ and $p$+Pb collisions. In contrast, while jet quenching is observed in Pb+Pb collisions, no evidence has been found in these small systems to date, raising fundamental questions about the nature of the system created in these collisions. The ATLAS experiment at the Large Hadron Collider has measured the yield of charged hadrons correlated with reconstructed jets in 0.36 nb$^{-1}$ of $p$+Pb and 3.6 pb$^{-1}$ of $pp$ collisions at 5.02 TeV. The yields of charged hadrons with $p_\mathrm{T}^\mathrm{ch} >0.5$ GeV near and opposite in azimuth to jets with $p_\mathrm{T}^\mathrm{jet} > 30$ or $60$ GeV, and the ratios of these yields between $p$+Pb and $pp$ collisions, $I_{p\mathrm{Pb}}$, are reported. The collision centrality of $p$+Pb events is categorized by the energy deposited by forward neutrons from the struck nucleus. The $I_{p\mathrm{Pb}}$ values are consistent with unity within a few percent for hadrons with $p_\mathrm{T}^\mathrm{ch} >4$ GeV at all centralities. These data provide new, strong constraints which preclude almost any parton energy loss in central $p$+Pb collisions.

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 suppression and correlations of dijets is performed using 3 $\mu$b$^{-1}$ of Xe+Xe data at $\sqrt{s_{\mathrm{NN}}} = 5.44$ TeV collected with the ATLAS detector at the LHC. Dijets with jets reconstructed using the $R=0.4$ anti-$k_t$ algorithm are measured differentially in jet $p_{\text{T}}$ over the range of 32 GeV to 398 GeV and the centrality of the collisions. Significant dijet momentum imbalance is found in the most central Xe+Xe collisions, which decreases in more peripheral collisions. Results from the measurement of per-pair normalized and absolutely normalized dijet $p_{\text{T}}$ balance are compared with previous Pb+Pb measurements at $\sqrt{s_{\mathrm{NN}}} =5.02$ TeV. The differences between the dijet suppression in Xe+Xe and Pb+Pb are further quantified by the ratio of pair nuclear-modification factors. The results are found to be consistent with those measured in Pb+Pb data when compared in classes of the same event activity and when taking into account the difference between the center-of-mass energies of the initial parton scattering process in Xe+Xe and Pb+Pb collisions. These results should provide input for a better understanding of the role of energy density, system size, path length, and fluctuations in the parton energy loss.

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.

Studies of the correlations of the two highest transverse momentum (leading) jets in individual Pb+Pb collision events can provide information about the mechanism of jet quenching by the hot and dense matter created in such collisions. In Pb+Pb and pp collisions at $\sqrt{s_{_\text{NN}}}$ = 5.02 TeV, measurements of the leading dijet transverse momentum ($p_{\mathrm{T}}$) correlations are presented. Additionally, measurements in Pb+Pb collisions of the dijet pair nuclear modification factors projected along leading and subleading jet $p_{\mathrm{T}}$ are made. The measurements are performed using the ATLAS detector at the LHC with 260 pb$^{-1}$ of pp data collected in 2017 and 2.2 nb$^{-1}$ of Pb+Pb data collected in 2015 and 2018. An unfolding procedure is applied to the two-dimensional leading and subleading jet $p_{\mathrm{T}}$ distributions to account for experimental effects in the measurement of both jets. Results are provided for dijets with leading jet $p_{\mathrm{T}}$ greater than 100 GeV. Measurements of the dijet-yield-normalized $x_{\mathrm{J}}$ distributions in Pb+Pb collisions show an increased fraction of imbalanced jets compared to pp collisions; these measurements are in agreement with previous measurements of the same quantity at 2.76 TeV in the overlapping kinematic range. Measurements of the absolutely-normalized dijet rate in Pb+Pb and pp collisions are also presented, and show that balanced dijets are significantly more suppressed than imbalanced dijets in Pb+Pb collisions. It is observed in the measurements of the pair nuclear modification factors that the subleading jets are significantly suppressed relative to leading jets with $p_{\mathrm{T}}$ between 100 and 316 GeV for all centralities in Pb+Pb collisions.

Differential cross-section measurements of $Z\gamma$ production in association with hadronic jets are presented, using the full 139 fb$^{-1}$ dataset of $\sqrt{s}=13$ TeV proton-proton collisions collected by the ATLAS detector during Run 2 of the LHC. Distributions are measured using events in which the $Z$ boson decays leptonically and the photon is usually radiated from an initial-state quark. Measurements are made in both one and two observables, including those sensitive to the hard scattering in the event and others which probe additional soft and collinear radiation. Different Standard Model predictions, from both parton-shower Monte Carlo simulation and fixed-order QCD calculations, are compared with the measurements. In general, good agreement is observed between data and predictions from MATRIX and MiNNLO$_\text{PS}$, as well as next-to-leading-order predictions from MadGraph5_aMC@NLO and Sherpa.

Jet shapes have been measured in inclusive jet production in proton-proton collisions at sqrt(s) = 7 TeV using 3 pb^{-1} of data recorded by the ATLAS experiment at the LHC. Jets are reconstructed using the anti-kt algorithm with transverse momentum 30 GeV < pT < 600 GeV and rapidity in the region |y| < 2.8. The data are corrected for detector effects and compared to several leading-order QCD matrix elements plus parton shower Monte Carlo predictions, including different sets of parameters tuned to model fragmentation processes and underlying event contributions in the final state. The measured jets become narrower with increasing jet transverse momentum and the jet shapes present a moderate jet rapidity dependence. Within QCD, the data test a variety of perturbative and non-perturbative effects. In particular, the data show sensitivity to the details of the parton shower, fragmentation, and underlying event models in the Monte Carlo generators. For an appropriate choice of the parameters used in these models, the data are well described.

A measurement of the cross section for the inclusive production of isolated prompt photons in pp collisions at a centre-of-mass energy sqrt(s) = 7TeV is presented. The measurement covers the pseudorapidity ranges |eta|<1.37 and 1.52<|eta|<1.81 in the transverse energy range 15 < E_T <100 GeV. The results are based on an integrated luminosity of 880 nb-1, collected with the ATLAS detector at the Large Hadron Collider. Photon candidates are identified by combining information from the calorimeters and from the inner tracker. Residual background in the selected sample is estimated from data based on the observed distribution of the transverse isolation energy in a narrow cone around the photon candidate. The results are compared to predictions from next-to-leading order perturbative QCD calculations.

Dijet angular distributions from the first LHC pp collisions at center-of-mass energy sqrt(s) = 7 TeV have been measured with the ATLAS detector. The dataset used for this analysis represents an integrated luminosity of 3.1 pb-1. Dijet $\chi$ distributions and centrality ratios have been measured up to dijet masses of 2.8 TeV, and found to be in good agreement with Standard Model predictions. Analysis of the $\chi$ distributions excludes quark contact interactions with a compositeness scale $\Lambda$ below 3.4 TeV, at 95% confidence level, significantly exceeding previous limits.

A search for new heavy particles manifested as resonances in two-jet final states is presented. The data were produced in 7 TeV proton-proton collisions by the Large Hadron Collider (LHC) and correspond to an integrated luminosity of 315 nb^-1 collected by the ATLAS detector. No resonances were observed. Upper limits were set on the product of cross section and signal acceptance for excited-quark (q*) production as a function of q* mass. These exclude at the 95% CL the q* mass interval 0.30 < mq* < 1.26 TeV, extending the reach of previous experiments.

A measurement of single top-quark production in the s-channel is performed in proton$-$proton collisions at a centre-of-mass energy of 13 TeV with the ATLAS detector at the CERN Large Hadron Collider. The dataset corresponds to an integrated luminosity of 139 fb$^{-1}$. The analysis is performed on events with an electron or muon, missing transverse momentum and exactly two $b$-tagged jets in the final state. A discriminant based on matrix element calculations is used to separate single-top-quark s-channel events from the main background contributions, which are top-quark pair production and $W$-boson production in association with jets. The observed (expected) signal significance over the background-only hypothesis is 3.3 (3.9) standard deviations, and the measured cross-section is $\sigma=8.2^{+3.5}_{-2.9}$ pb, consistent with the Standard Model prediction of $\sigma^{\mathrm{SM}}=10.32^{+0.40}_{-0.36}$ pb.

Measurements of differential cross sections are presented for inclusive isolated-photon production in $pp$ collisions at a centre-of-mass energy of 13 TeV provided by the LHC and using 139 fb$^{-1}$ of data recorded by the ATLAS experiment. The cross sections are measured as functions of the photon transverse energy in different regions of photon pseudorapidity. The photons are required to be isolated by means of a fixed-cone method with two different cone radii. The dependence of the inclusive-photon production on the photon isolation is investigated by measuring the fiducial cross sections as functions of the isolation-cone radius and the ratios of the differential cross sections with different radii in different regions of photon pseudorapidity. The results presented in this paper constitute an improvement with respect to those published by ATLAS earlier: the measurements are provided for different isolation radii and with a more granular segmentation in photon pseudorapidity that can be exploited in improving the determination of the proton parton distribution functions. These improvements provide a more in-depth test of the theoretical predictions. Next-to-leading-order QCD predictions from JETPHOX and SHERPA and next-to-next-to-leading-order QCD predictions from NNLOJET are compared to the measurements, using several parameterisations of the proton parton distribution functions. The measured cross sections are well described by the fixed-order QCD predictions within the experimental and theoretical uncertainties in most of the investigated phase-space region.

A search for pair-produced vector-like quarks using events with exactly one lepton ($e$ or $\mu$), at least four jets including at least one $b$-tagged jet, and large missing transverse momentum is presented. Data from proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV, recorded by the ATLAS detector at the LHC from 2015 to 2018 and corresponding to an integrated luminosity of 139 fb$^{-1}$, are analysed. Vector-like partners $T$ and $B$ of the top and bottom quarks are considered, as is a vector-like $X$ with charge +5/3, assuming their decay into a $W$, $Z$, or Higgs boson and a third-generation quark. No significant deviations from the Standard Model expectation are observed. Upper limits on the production cross-section of $T$ and $B$ quark pairs as a function of their mass are derived for various decay branching ratio scenarios. The strongest lower limits on the masses are 1.59 TeV assuming mass-degenerate VLQs and branching ratios corresponding to the weak-isospin doublet model, and 1.47 TeV (1.46 TeV) for exclusive $T \rightarrow Zt$ ($B/X \rightarrow Wt$) decays. In addition, lower limits on the $T$ and $B$ quark masses are derived for all possible branching ratios.

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.

Jet cross sections have been measured for the first time in proton-proton collisions at a centre-of-mass energy of 7 TeV using the ATLAS detector. The measurement uses an integrated luminosity of 17 nb-1 recorded at the Large Hadron Collider. The anti-kt algorithm is used to identify jets, with two jet resolution parameters, R = 0.4 and 0.6. The dominant uncertainty comes from the jet energy scale, which is determined to within 7% for central jets above 60 GeV transverse momentum. Inclusive single-jet differential cross sections are presented as functions of jet transverse momentum and rapidity. Dijet cross sections are presented as functions of dijet mass and the angular variable $\chi$. The results are compared to expectations based on next-to-leading-order QCD, which agree with the data, providing a validation of the theory in a new kinematic regime.

A search for a $WZ$ resonance, in the fully leptonic final state (electrons and muons), is performed using 139 fb$^{-1}$ of data collected at a centre-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The results are interpreted in terms of a singly charged Higgs boson of the Georgi$-$Machacek model, produced by $WZ$ fusion, and of a Heavy Vector Triplet, with the resonance produced by $WZ$ fusion or the Drell$-$Yan process. No significant excess over the Standard Model predictions is observed and limits are set on the production cross-section times branching ratio as a function of the resonance mass for these processes.

A study of $B_c^+\to J/\psi D_s^+$ and $B_c^+\to J/\psi D_s^{*+}$ decays using 139 fb$^{-1}$ of integrated luminosity collected with the ATLAS detector from $\sqrt{s} = 13$ TeV $pp$ collisions at the LHC is presented. The ratios of the branching fractions of the two decays to the branching fraction of the $B_c^+\to J/\psi \pi^+$ decay are measured: $\mathcal B(B_c^+\to J/\psi D_s^+)/\mathcal B(B_c^+\to J/\psi \pi^+) = 2.76\pm 0.47$ and $\mathcal B(B_c^+\to J/\psi D_s^{*+})/\mathcal B(B_c^+\to J/\psi \pi^+) = 5.33\pm 0.96$. The ratio of the branching fractions of the two decays is found to be $\mathcal B(B_c^+\to J/\psi D_s^{*+})/\mathcal B(B_c^+\to J/\psi D_s^+) = 1.93\pm0.26$. For the $B_c^+\to J/\psi D_s^{*+}$ decay, the transverse polarization fraction, $\Gamma_{\pm\pm}/\Gamma$, is measured to be $0.70\pm0.11$. The reported uncertainties include both the statistical and systematic components added in quadrature. The precision of the measurements exceeds that in all previous studies of these decays. These results supersede those obtained in the earlier ATLAS study of the same decays with $\sqrt{s} = 7$ and 8 TeV $pp$ collision data. A comparison with available theoretical predictions for the measured quantities is presented.

A search for pair-produced scalar and vector leptoquarks decaying into quarks and leptons of different generations is presented. It uses the full LHC Run 2 (2015-2018) data set of 139 fb$^{-1}$ collected with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV. Scalar leptoquarks with charge -(1/3)e as well as scalar and vector leptoquarks with charge +(2/3)e are considered. All possible decays of the pair-produced leptoquarks into quarks of the third generation (t, b) and charged or neutral leptons of the first or second generation ($e, \mu, \nu$) with exactly one electron or muon in the final state are investigated. No significant deviations from the Standard Model expectation are observed. Upper limits on the production cross-section are provided for eight models as a function of the leptoquark mass and the branching ratio of the leptoquark into the charged or neutral lepton. In addition, lower limits on the leptoquark masses are derived for all models across a range of branching ratios. Two of these models have the goal of providing an explanation for the recent B-anomalies. In both models, a vector leptoquark decays into charged and neutral leptons of the second generation with a similar branching fraction. Lower limits of 1980 GeV and 1710 GeV are set on the leptoquark mass for these two models.

A search for diphoton resonances in the mass range between 10 and 70 GeV with the ATLAS experiment at the Large Hadron Collider (LHC) is presented. The analysis is based on $pp$ collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at a centre-of-mass energy of 13 TeV recorded from 2015 to 2018. Previous searches for diphoton resonances at the LHC have explored masses down to 65 GeV, finding no evidence of new particles. This search exploits the particular kinematics of events with pairs of closely spaced photons reconstructed in the detector, allowing examination of invariant masses down to 10 GeV. The presented strategy covers a region previously unexplored at hadron colliders because of the experimental challenges of recording low-energy photons and estimating the backgrounds. No significant excess is observed and the reported limits provide the strongest bound on promptly decaying axion-like particles coupling to gluons and photons for masses between 10 and 70 GeV.

First measurements of the W -> lnu and Z/gamma* -> ll (l = e, mu) production cross sections in proton-proton collisions at sqrt(s) = 7 TeV are presented using data recorded by the ATLAS experiment at the LHC. The results are based on 2250 W -> lnu and 179 Z/gamma* -> ll candidate events selected from a data set corresponding to an integrated luminosity of approximately 320 nb-1. The measured total W and Z/gamma*-boson production cross sections times the respective leptonic branching ratios for the combined electron and muon channels are $\stotW$ * BR(W -> lnu) = 9.96 +- 0.23(stat) +- 0.50(syst) +- 1.10(lumi) nb and $\stotZg$ * BR(Z/gamma* -> ll) = 0.82 +- 0.06(stat) +- 0.05(syst) +- 0.09(lumi) nb (within the invariant mass window 66 < m_ll < 116 GeV). The W/Z cross-section ratio is measured to be 11.7 +- 0.9(stat) +- 0.4(syst). In addition, measurements of the W+ and W- production cross sections and of the lepton charge asymmetry are reported. Theoretical predictions based on NNLO QCD calculations are found to agree with the measurements.

This paper presents direct searches for lepton flavour violation in Higgs boson decays, $H\rightarrow e\tau$ and $H\rightarrow\mu\tau$, performed using data collected 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 138 fb$^{-1}$. Leptonic ($\tau \rightarrow \ell \nu_\ell \nu_\tau$) and hadronic ($\tau \rightarrow $ hadrons $ \nu_\tau$) decays of the $\tau$-lepton are considered. Two background estimation techniques are employed: the MC-template method, based on data-corrected simulation samples, and the Symmetry method, based on exploiting the symmetry between electrons and muons in the Standard Model backgrounds. No significant excess of events is observed and the results are interpreted as upper limits on lepton-flavour-violating branching ratios of the Higgs boson. The observed (expected) upper limits set on the branching ratios at 95% confidence level, $\mathcal{B}(H\rightarrow e\tau)<0.20\%$ (0.12%) and $\mathcal{B}(H\rightarrow \mu\tau)<0.18\%$ (0.09%), are obtained with the MC-template method from a simultaneous measurement of potential $H \rightarrow e\tau$ and $H \rightarrow\mu\tau$ signals. The best-fit branching ratio difference, $\mathcal{B}(H\rightarrow \mu\tau)- \mathcal{B}(H\rightarrow e\tau)$, measured with the Symmetry method in the channel where the $\tau$-lepton decays to leptons, is (0.25 $\pm$ 0.10)%, compatible with a value of zero within 2.5$\sigma$.

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.

A generic search for resonances is performed with events containing a $Z$ boson with transverse momentum greater than 100 GeV, decaying into $e^+e^-$ or $\mu^+\mu^-$. The analysed data collected with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider correspond to an integrated luminosity of 139 fb$^{-1}$. Two invariant mass distributions are examined for a localised excess relative to the expected Standard Model background in six independent event categories (and their inclusive sum) to increase the sensitivity. No significant excess is observed. Exclusion limits at 95% confidence level are derived for two cases: a model-independent interpretation of Gaussian-shaped resonances with the mass width between 3% and 10% of the resonance mass, and a specific heavy vector triplet model with the decay mode $W'\to ZW \to \ell\ell qq$.

A search for charged leptons with large impact parameters using 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV $pp$ collision data from the ATLAS detector at the LHC is presented, addressing a long-standing gap in coverage of possible new physics signatures. Results are consistent with the background prediction. This search provides unique sensitivity to long-lived scalar supersymmetric lepton-partners (sleptons). For lifetimes of 0.1 ns, selectron, smuon and stau masses up to 720 GeV, 680 GeV, and 340 GeV are respectively excluded at 95% confidence level, drastically improving on the previous best limits from LEP.

Two-particle long-range azimuthal correlations are measured in photonuclear collisions using 1.7 nb$^{-1}$ of 5.02 TeV Pb+Pb collision data collected by the ATLAS experiment at the LHC. Candidate events are selected using a dedicated high-multiplicity photonuclear event trigger, a combination of information from the zero-degree calorimeters and forward calorimeters, and from pseudorapidity gaps constructed using calorimeter energy clusters and charged-particle tracks. Distributions of event properties are compared between data and Monte Carlo simulations of photonuclear processes. Two-particle correlation functions are formed using charged-particle tracks in the selected events, and a template-fitting method is employed to subtract the non-flow contribution to the correlation. Significant nonzero values of the second- and third-order flow coefficients are observed and presented as a function of charged-particle multiplicity and transverse momentum. The results are compared with flow coefficients obtained in proton-proton and proton-lead collisions in similar multiplicity ranges, and with theoretical expectations. The unique initial conditions present in this measurement provide a new way to probe the origin of the collective signatures previously observed only in hadronic collisions.