A search for supersymmetry in events with four or more charged leptons (electrons, muons and $\tau$-leptons) is presented. The analysis uses a data sample corresponding to $139\,\mbox{fb\(^{-1}\)}$ of proton-proton collisions delivered by the Large Hadron Collider at $\sqrt{s}=13$ TeV and recorded by the ATLAS detector. Four-lepton signal regions with up to two hadronically decaying $\tau$-leptons are designed to target several supersymmetric models, while a general five-lepton signal region targets any new physics phenomena leading to a final state with five charged leptons. Data yields are consistent with Standard Model expectations and results are used to set upper limits on contributions from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of general gauge-mediated supersymmetry, excluding higgsino masses up to $540$ GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of $1.6$ TeV, $1.2$ TeV, and $2.5$ TeV are placed on wino, slepton and gluino masses, respectively.

A search for dark-matter particles in events with large missing transverse momentum and a Higgs boson candidate decaying into two photons is reported. The search uses $139$ fb$^{-1}$ of proton-proton collision data collected at $\sqrt{s}=13$ TeV with the ATLAS detector at the CERN LHC between 2015 and 2018. No significant excess of events over the Standard Model predictions is observed. The results are interpreted by extracting limits on three simplified models that include either vector or pseudoscalar mediators and predict a final state with a pair of dark-matter candidates and a Higgs boson decaying into two photons.

Using 1.8 fb-1 of pp collisions at a center-of-mass energy of 7 TeV recorded by the ATLAS detector at the Large Hadron Collider, we present measurements of the production cross sections of Upsilon(1S,2S,3S) mesons. Upsilon mesons are reconstructed using the di-muon decay mode. Total production cross sections for p_T<70 GeV and in the rapidity interval |Upsilon|<2.25 are measured to be 8.01+-0.02+-0.36+-0.31 nb, 2.05+-0.01+-0.12+-0.08 nb, 0.92+-0.01+-0.07+-0.04 nb respectively, with uncertainties separated into statistical, systematic, and luminosity measurement effects. In addition, differential cross section times di-muon branching fractions for Upsilon(1S), Upsilon(2S), and Upsilon(3S) as a function of Upsilon transverse momentum p_T and rapidity are presented. These cross sections are obtained assuming unpolarized production. If the production polarization is fully transverse or longitudinal with no azimuthal dependence in the helicity frame the cross section may vary by approximately +-20%. If a non-trivial azimuthal dependence is considered, integrated cross sections may be significantly enhanced by a factor of two or more. We compare our results to several theoretical models of Upsilon meson production, finding that none provide an accurate description of our data over the full range of Upsilon transverse momenta accessible with this dataset.

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.

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.

Exclusive dimuon production in ultraperipheral collisions (UPC), resulting from photon-photon interactions in the strong electromagnetic fields of colliding high-energy lead nuclei, $\mathrm{PbPb}(\gamma\gamma) \rightarrow \mu^+\mu^- (\mathrm{Pb}^{(\star)}\mathrm{Pb}^{(\star)} )$, is studied using $\mathcal{L}_{\mathrm{int}} = 0.48$ nb$^{-1}$ of $\sqrt{s_\mathrm{NN}}=5.02$ TeV lead-lead collision data at the LHC with the ATLAS detector. Dimuon pairs are measured in the fiducial region $p_{\mathrm{T}\mu} > 4$ GeV, $|\eta_{\mu}| < 2.4$, invariant mass $m_{\mu\mu} > 10$ GeV, and $p_{\mathrm{T,\mu\mu}} < 2$ GeV. The primary background from single-dissociative processes is extracted from the data using a template fitting technique. Differential cross sections are presented as a function of $m_{\mu\mu}$, absolute pair rapidity ($|y_{\mu\mu}|$), scattering angle in the dimuon rest frame ($|\cos \vartheta^{\star}_{\mu\mu}|$) and the colliding photon energies. The total cross section of the UPC $\gamma \gamma \rightarrow \mu^{+}\mu^{-}$ process in the fiducial volume is measured to be $\sigma_{\mathrm{fid}}^{\mu\mu} = 34.1 \! \pm \! 0.3 \mathrm{(stat.)} \! \pm \! 0.7 \mathrm{(syst.)}$ $\mu\mathrm{b}$. Generally good agreement is found with calculations from STARlight, which incorporate the leading-order Breit-Wheeler process with no final-state effects, albeit differences between the measurements and theoretical expectations are observed. In particular, the measured cross sections at larger $|y_{\mu\mu}|$ are found to be about 10-20% larger in data than in the calculations, suggesting the presence of larger fluxes of photons in the initial state. Modification of the dimuon cross sections in the presence of forward and/or backward neutron production is also studied and is found to be associated with a harder incoming photon spectrum, consistent with expectations.

An inclusive search for a new-physics signature of lepton-jet resonances has been performed by the ATLAS experiment. Scalar leptoquarks, pair-produced in $pp$ collisions at $\sqrt{s}$ = 13 TeV at the Large Hadron Collider, have been considered. An integrated luminosity of 3.2 fb$^{-1}$, corresponding to the full 2015 dataset was used. First (second) generation leptoquarks were sought in events with two electrons (muons) and two or more jets. The observed event yield in each channel is consistent with Standard Model background expectations. The observed (expected) lower limits on the leptoquark mass at 95% confidence level are 1100 GeV and 1050 GeV (1160 GeV and 1040 GeV) for first and second generation leptoquarks, respectively, assuming a branching ratio into a charged lepton and a quark of 100%. Upper limits on the aforementioned branching ratio are also given as a function of leptoquark mass. Compared with the results of earlier ATLAS searches, the sensitivity is increased for leptoquark masses above 860 GeV, and the observed exclusion limits confirm and extend the published results.

The prompt and non-prompt production cross-sections for the chi_c1 and chi_c2 charmonium states are measured in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector at the LHC using 4.5 fb^-1 of integrated luminosity. The chi_c states are reconstructed through the radiative decay chi_c -> J/psi gamma (with J/psi -> mu+mu-) where photons are reconstructed from gamma -> e+e- conversions. The production rate of the chi_c2 state relative to the chi_c1 state is measured for prompt and non-prompt chi_c as a function of J/psi transverse momentum. The prompt chi_c cross-sections are combined with existing measurements of prompt J/psi production to derive the fraction of prompt J/psi produced in feed-down from chi_c decays. The fractions of chi_c1 and chi_c2 produced in b-hadron decays are also measured.

This Letter presents a search for magnetic monopoles with the ATLAS detector at the CERN Large Hadron Collider using an integrated luminosity of 2.0 fb^-1 of pp collisions recorded at a center-of-mass energy of sqrt(s)=7 TeV. No event is found in the signal region, leading to an upper limit on the production cross section at 95% confidence level of 1.6/epsilon fb for Dirac magnetic monopoles with the minimum unit magnetic charge and with mass between 200 GeV and 1500 GeV, where epsilon is the monopole reconstruction efficiency. The efficiency epsilon is high and uniform in the fiducial region given by pseudorapidity |eta|<1.37 and transverse kinetic energy 600-700<E^kin sin(theta)<1400 GeV. The minimum value of 700 GeV is for monopoles of mass 200 GeV, whereas the minimum value of 600 GeV is applicable for higher mass monopoles. Therefore, the upper limit on the production cross section at 95% confidence level is 2 fb in this fiducial region. Assuming the kinematic distributions from Drell-Yan pair production of spin-1/2 Dirac magnetic monopoles, the efficiency is in the range 1%-10%, leading to an upper limit on the cross section at 95% confidence level that varies from 145 fb to 16 fb for monopoles with mass between 200 GeV and 1200 GeV. This limit is weaker than the fiducial limit because most of these monopoles lie outside the fiducial region.

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.

Measurements of four-lepton differential and integrated fiducial cross-sections in events with two same-flavour, opposite-charge electron or muon pairs are presented. The data correspond to 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collisions, collected by the ATLAS detector during Run 2 of the Large Hadron Collider (2015-2018). The final state has contributions from a number of interesting Standard Model processes that dominate in different four-lepton invariant mass regions, including single $Z$ boson production, Higgs boson production and on-shell $ZZ$ production, with a complex mix of interference terms, and possible contributions from physics beyond the Standard Model. The differential cross-sections include the four-lepton invariant mass inclusively, in slices of other kinematic variables, and in different lepton flavour categories. Also measured are dilepton invariant masses, transverse momenta, and angular correlation variables, in four regions of four-lepton invariant mass, each dominated by different processes. The measurements are corrected for detector effects and are compared with state-of-the-art Standard Model calculations, which are found to be consistent with the data. The $Z\rightarrow 4\ell$ branching fraction is extracted, giving a value of $\left(4.41 \pm 0.30\right) \times 10^{-6}$. Constraints on effective field theory parameters and a model based on a spontaneously broken $B-L$ gauge symmetry are also evaluated. Further reinterpretations can be performed with the provided information.

Light-by-light scattering ($\gamma\gamma\rightarrow\gamma\gamma$) is a quantum-mechanical process that is forbidden in the classical theory of electrodynamics. This reaction is accessible at the Large Hadron Collider thanks to the large electromagnetic field strengths generated by ultra-relativistic colliding lead (Pb) ions. Using 480 $\mu$b$^{-1}$ of Pb+Pb collision data recorded at a centre-of-mass energy per nucleon pair of 5.02 TeV by the ATLAS detector, the ATLAS Collaboration reports evidence for the $\gamma\gamma\rightarrow\gamma\gamma$ reaction. A total of 13 candidate events are observed with an expected background of 2.6$\pm$0.7 events. After background subtraction and analysis corrections, the fiducial cross section of the process $\textrm{Pb+Pb}\,(\gamma\gamma)\rightarrow \textrm{Pb}^{(\ast)}\textrm{+}\textrm{Pb}^{(\ast)}\,\gamma\gamma$, for photon transverse energy $E_{\mathrm{T}}>$3 GeV, photon absolute pseudorapidity $|\eta|<$2.4, diphoton invariant mass greater than 6 GeV, diphoton transverse momentum lower than 2 GeV and diphoton acoplanarity below 0.01, is measured to be 70 $\pm$ 24 (stat.) $\pm$ 17 (syst.) nb, which is in agreement with Standard Model predictions.

A search is presented for photonic signatures, motivated by generalized models of gauge-mediated supersymmetry breaking. This search makes use of proton-proton collision data at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ recorded by the ATLAS detector at the LHC, and it explores models dominated by both strong and electroweak production of supersymmetric partner states. Experimental signatures incorporating an isolated photon and significant missing transverse momentum are explored. These signatures include events with an additional photon or additional jet activity not associated with any specific underlying quark flavor. No significant excess of events is observed above the Standard Model prediction, and 95% confidence-level upper limits of between 0.083 fb and 0.32 fb are set on the visible cross section of contributions from physics beyond the Standard Model. These results are interpreted in terms of lower limits on the masses of gluinos, squarks, and gauginos in the context of generalized models of gauge-mediated supersymmetry, which reach as high as 2.3 TeV for strongly produced and 1.3 TeV for weakly produced supersymmetric partner pairs.

This paper presents a measurement of the triple-differential cross section for the Drell--Yan process $Z/\gamma^*\rightarrow \ell^+\ell^-$ where $\ell$ is an electron or a muon. The measurement is performed for invariant masses of the lepton pairs, $m_{\ell\ell}$, between $46$ and $200$ GeV using a sample of $20.2$ fb$^{-1}$ of $pp$ collisions data at a centre-of-mass energy of $\sqrt{s}=8$ TeV collected by the ATLAS detector at the LHC in 2012. The data are presented in bins of invariant mass, absolute dilepton rapidity, $|y_{\ell\ell}|$, and the angular variable $\cos\theta^{*}$ between the outgoing lepton and the incoming quark in the Collins--Soper frame. The measurements are performed in the range $|y_{\ell\ell}|<2.4$ in the muon channel, and extended to $|y_{\ell\ell}|<3.6$ in the electron channel. The cross sections are used to determine the $Z$ boson forward-backward asymmetry as a function of $|y_{\ell\ell}|$ and $m_{\ell\ell}$. The measurements achieve high-precision, below the percent level in the pole region, excluding the uncertainty in the integrated luminosity, and are in agreement with predictions. These precision data are sensitive to the parton distribution functions and the effective weak mixing angle.

A search is performed for the electroweak pair production of charginos and associated production of a chargino and neutralino, each of which decays through an $R$-parity-violating coupling into a lepton and a $W$, $Z$, or Higgs boson. The trilepton invariant-mass spectrum is constructed from events with three or more leptons, targeting chargino decays that include an electron or muon and a leptonically decaying $Z$ boson. The analyzed dataset corresponds to an integrated luminosity of 139 fb$^{-1}$ of proton-proton collision data produced by the Large Hadron Collider at a center-of-mass energy of $\sqrt{s}$ = 13 TeV and collected by the ATLAS experiment between 2015 and 2018. The data are found to be consistent with predictions from the Standard Model. The results are interpreted as limits at 95% confidence level on model-independent cross sections for processes beyond the Standard Model. Limits are also set on the production of charginos and neutralinos for a Minimal Supersymmetric Standard Model with an approximate $B$-$L$ symmetry. Charginos and neutralinos with masses between 100 GeV and 1100 GeV are excluded depending on the assumed decay branching fractions into a lepton (electron, muon, or $\tau$-lepton) plus a boson ($W$, $Z$, or Higgs).

A study of WZ production in proton-proton collisions at sqrt(s) = 7 TeV is presented using data corresponding to an integrated luminosity of 4.6 fb^-1 collected with the ATLAS detector at the Large Hadron Collider in 2011. In total, 317 candidates, with a background expectation of 68+/-10 events, are observed in double-leptonic decay final states with electrons, muons and missing transverse momentum. The total cross-section is determined to be sigma_WZ(tot) = 19.0+1.4/-1.3(stat.)+/-0.9(syst.)+/-0.4(lumi.) pb, consistent with the Standard Model expectation of 17.6+1.1/-1.0 pb. Limits on anomalous triple gauge boson couplings are derived using the transverse momentum spectrum of Z bosons in the selected events. The cross section is also presented as a function of Z boson transverse momentum and diboson invariant mass.

A measurement of event-plane correlations involving two or three event planes of different order is presented as a function of centrality for 7 ub-1 Pb+Pb collision data at sqrt(s_NN)=2.76 TeV, recorded by the ATLAS experiment at the LHC. Fourteen correlators are measured using a standard event-plane method and a scalar-product method, and the latter method is found to give a systematically larger correlation signal. Several different trends in the centrality dependence of these correlators are observed. These trends are not reproduced by predictions based on the Glauber model, which includes only the correlations from the collision geometry in the initial state. Calculations that include the final-state collective dynamics are able to describe qualitatively, and in some cases also quantitatively, the centrality dependence of the measured correlators. These observations suggest that both the fluctuations in the initial geometry and non-linear mixing between different harmonics in the final state are important for creating these correlations in momentum space.

We present a measurement of two-particle angular correlations in proton-proton collisions at sqrt(s) = 900 GeV and 7 TeV. The collision events were collected during 2009 and 2010 with the ATLAS detector at the Large Hadron Collider using a single-arm minimum bias trigger. Correlations are measured for charged particles produced in the kinematic range of transverse momentum pT > 100 MeV and pseudorapidity |eta| < 2.5. A complex structure in pseudorapidity and azimuth is observed at both collision energies. Results are compared to Pythia 8 and Herwig++ as well as to the AMBT2B, DW and Perugia 2011 tunes of Pythia 6. The data are not satisfactorily described by any of these models.

A $6.8 \ {\rm nb^{-1}}$ sample of $pp$ collision data collected under low-luminosity conditions at $\sqrt{s} = 7$ TeV by the ATLAS detector at the Large Hadron Collider is used to study diffractive dijet production. Events containing at least two jets with $p_\mathrm{T} > 20$ GeV are selected and analysed in terms of variables which discriminate between diffractive and non-diffractive processes. Cross sections are measured differentially in $\Delta\eta^F$, the size of the observable forward region of pseudorapidity which is devoid of hadronic activity, and in an estimator, $\tilde{\xi}$, of the fractional momentum loss of the proton assuming single diffractive dissociation ($pp \rightarrow pX$). Model comparisons indicate a dominant non-diffractive contribution up to moderately large $\Delta\eta^F$ and small $\tilde{\xi}$, with a diffractive contribution which is significant at the highest $\Delta\eta^F$ and the lowest $\tilde{\xi}$. The rapidity-gap survival probability is estimated from comparisons of the data in this latter region with predictions based on diffractive parton distribution functions.

A search is presented for photonic signatures motivated by generalised models of gauge-mediated supersymmetry breaking. This search makes use of $20.3{\rm fb}^{-1}$ of proton-proton collision data at $\sqrt{s}=8$ TeV recorded by the ATLAS detector at the LHC, and explores models dominated by both strong and electroweak production of supersymmetric partner states. Four experimental signatures incorporating an isolated photon and significant missing transverse momentum are explored. These signatures include events with an additional photon, lepton, $b$-quark jet, or jet activity not associated with any specific underlying quark flavor. No significant excess of events is observed above the Standard Model prediction and model-dependent 95% confidence-level exclusion limits are set.

The four-lepton ($4\ell$, $\ell = e, \mu$) production cross section is measured in the mass range from 80 to 1000 GeV using 20.3 fb$^{-1}$ of data in $pp$ collisions at $\sqrt{s}=8$ TeV collected with the ATLAS detector at the LHC. The $4\ell$ events are produced in the decays of resonant $Z$ and Higgs bosons and the non-resonant $ZZ$ continuum originating from $q\bar q$, $gg$, and $qg$ initial states. A total of 476 signal candidate events are observed with a background expectation of $26.2 \pm 3.6$ events, enabling the measurement of the integrated cross section and the differential cross section as a function of the invariant mass and transverse momentum of the four-lepton system. In the mass range above $180$ GeV, assuming the theoretical constraint on the $q\bar q$ production cross section calculated with perturbative NNLO QCD and NLO electroweak corrections, the signal strength of the gluon-fusion component relative to its leading-order prediction is determined to be $\mu_{gg}=2.4 \pm 1.0 (stat.) \pm 0.5 (syst.)\pm 0.8 (theory)$.

Pseudorapidity gap distributions in proton-proton collisions at sqrt(s) = 7 TeV are studied using a minimum bias data sample with an integrated luminosity of 7.1 inverse microbarns. Cross sections are measured differentially in terms of Delta eta F, the larger of the pseudorapidity regions extending to the limits of the ATLAS sensitivity, at eta = +/- 4.9, in which no final state particles are produced above a transverse momentum threshold p_T Cut. The measurements span the region 0 < Delta eta F < 8 for 200 < p_T Cut < 800 MeV. At small Delta eta F, the data test the reliability of hadronisation models in describing rapidity and transverse momentum fluctuations in final state particle production. The measurements at larger gap sizes are dominated by contributions from the single diffractive dissociation process (pp -> Xp), enhanced by double dissociation (pp -> XY) where the invariant mass of the lighter of the two dissociation systems satisfies M_Y <~ 7 GeV. The resulting cross section is d sigma / d Delta eta F ~ 1 mb for Delta eta F >~ 3. The large rapidity gap data are used to constrain the value of the pomeron intercept appropriate to triple Regge models of soft diffraction. The cross section integrated over all gap sizes is compared with other LHC inelastic cross section measurements.

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.

The result of a search for pair production of the supersymmetric partner of the Standard Model bottom quark ($\tilde{b}_1$) is reported. The search uses 3.2 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}=$13 TeV collected by the ATLAS experiment at the Large Hadron Collider in 2015. Bottom squarks are searched for in events containing large missing transverse momentum and exactly two jets identified as originating from $b$-quarks. No excess above the expected Standard Model background yield is observed. Exclusion limits at 95% confidence level on the mass of the bottom squark are derived in phenomenological supersymmetric $R$-parity-conserving models in which the $\tilde{b}_1$ is the lightest squark and is assumed to decay exclusively via $\tilde{b}_1 \rightarrow b \tilde{\chi}_1^0$, where $\tilde{\chi}_1^0$ is the lightest neutralino. The limits significantly extend previous results; bottom squark masses up to 800 (840) GeV are excluded for the $\tilde{\chi}_1^0$ mass below 360 (100) GeV whilst differences in mass above 100 GeV between the $\tilde{b}_1$ and the $\tilde{\chi}_1^0$ are excluded up to a $\tilde{b}_1$ mass of 500 GeV.

The modification of the production of $J/\psi$, $\psi(\mathrm{2S})$, and $\mit{\Upsilon}(n\mathrm{S})$ ($n = 1, 2, 3$) in $p$+Pb collisions with respect to their production in $pp$ collisions has been studied. The $p$+Pb and $pp$ datasets used in this paper correspond to integrated luminosities of $28$ $\mathrm{nb}^{-1}$ and $25$ $\mathrm{pb}^{-1}$ respectively, collected in 2013 and 2015 by the ATLAS detector at the LHC, both at a centre-of-mass energy per nucleon pair of 5.02 TeV. The quarkonium states are reconstructed in the dimuon decay channel. The yields of $J/\psi$ and $\psi(\mathrm{2S})$ are separated into prompt and non-prompt sources. The measured quarkonium differential cross sections are presented as a function of rapidity and transverse momentum, as is the nuclear modification factor, $R_{p\mathrm{Pb}}$ for $J/\psi$ and $\mit{\Upsilon}(\mathrm{1S})$. No significant modification of the $J/\psi$ production is observed while $\mit{\Upsilon}(\mathrm{1S})$ production is found to be suppressed at low transverse momentum in $p$+Pb collisions relative to $pp$ collisions. The production of excited charmonium and bottomonium states is found to be suppressed relative to that of the ground states in central $p$+Pb collisions.