A search is presented for a heavy resonance $Y$ decaying into a Standard Model Higgs boson $H$ and a new particle $X$ in a fully hadronic final state. The full Large Hadron Collider Run 2 dataset of proton-proton collisions at $\sqrt{s}= 13$ TeV collected by the ATLAS detector from 2015 to 2018 is used, and corresponds to an integrated luminosity of 139 fb$^{-1}$. The search targets the high $Y$-mass region, where the $H$ and $X$ have a significant Lorentz boost in the laboratory frame. A novel signal region is implemented using anomaly detection, where events are selected solely because of their incompatibility with a learned background-only model. It is defined using a jet-level tagger for signal-model-independent selection of the boosted $X$ particle, representing the first application of fully unsupervised machine learning to an ATLAS analysis. Two additional signal regions are implemented to target a benchmark $X$ decay into two quarks, covering topologies where the $X$ is reconstructed as either a single large-radius jet or two small-radius jets. The analysis selects Higgs boson decays into $b\bar{b}$, and a dedicated neural-network-based tagger provides sensitivity to the boosted heavy-flavor topology. No significant excess of data over the expected background is observed, and the results are presented as upper limits on the production cross section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) for signals with $m_Y$ between 1.5 and 6 TeV and $m_X$ between 65 and 3000 GeV.

A search for leptoquark pair production decaying into $te^- \bar{t}e^+$ or $t\mu^- \bar{t}\mu^+$ in final states with multiple leptons is presented. The search is based on a dataset of $pp$ collisions at $\sqrt{s}$=13 TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 139~fb$^{-1}$. Events are selected with two or more light leptons (electron or muon) and at least two jets out of which at least one jet is identified as coming from a $b$-hadron. Four signal regions, with the requirement of at least three light leptons, are considered based on the number of leptons of a given flavour. The main background processes are estimated using dedicated control regions in a simultaneous fit with the signal regions to data. No excess above the Standard Model background prediction is observed and 95% confidence level limits on the production cross section times branching ratio are derived as a function of the leptoquark mass. Under the assumption of exclusive decays into $te^{-}$ ($t\mu^{-}$), the corresponding lower limit on the scalar mixed-generation leptoquark mass $m_{\mathrm{LQ}_{\mathrm{mix}}^{\mathrm{d}}}$ is at 1.58 (1.59) TeV and on the vector leptoquark mass $m_{\tilde{U}_1}$ at 1.67 (1.67) TeV in the minimal coupling scenario and at 1.95 (1.95) TeV in the Yang-Mills scenario.

A search is presented for the pair production of heavy vector-like $T$ and $B$ quarks, primarily targeting the $T$ quark decays to a $W$ boson and a $b$-quark. The search is based on $36.1$ fb$^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 13 TeV recorded in 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider. Data are analysed in the lepton-plus-jets final state, including at least one $b$-tagged jet and a large-radius jet identified as originating from the hadronic decay of a high-momentum $W$ boson. No significant deviation from the Standard Model expectation is observed in the reconstructed $T$ mass distribution. The observed 95% confidence level lower limit on the $T$ mass are 1350 GeV assuming 100% branching ratio to $Wb$. In the SU(2) singlet scenario, the lower mass limit is 1170 GeV. This search is also sensitive to a heavy vector-like $B$ quark decaying to $Wt$ and other final states. The results are thus reinterpreted to provide a 95% confidence level lower limit on the $B$ quark mass at 1250 GeV assuming 100% branching ratio to $Wt$; in the SU(2) singlet scenario, the limit is 1080 GeV. Mass limits on both $T$ and $B$ production are also set as a function of the decay branching ratios. The 100% branching ratio limits are found to be applicable to heavy vector-like $Y$ and $X$ production that decay to $Wb$ and $Wt$, respectively.

A measurement of the $ZZ$ production in the $\ell^{-}\ell^{+}\ell^{\prime -}\ell^{\prime +}$ and $\ell^{-}\ell^{+}\nu\bar{\nu}$ channels $(\ell = e, \mu)$ in proton--proton collisions at $\sqrt{s} = 8$ TeV at the Large Hadron Collider at CERN, using data corresponding to an integrated luminosity of 20.3 fb$^{-1}$ collected by the ATLAS experiment in 2012 is presented. The fiducial cross sections for $ZZ\to\ell^{-}\ell^{+}\ell^{\prime -}\ell^{\prime +}$ and $ZZ\to \ell^{-}\ell^{+}\nu\bar{\nu}$ are measured in selected phase-space regions. The total cross section for $ZZ$ events produced with both $Z$ bosons in the mass range 66 to 116 GeV is measured from the combination of the two channels to be $7.3\pm0.4\textrm{(stat)}\pm0.3\textrm{(syst)}\pm0.2\textrm{(lumi)}$ pb, which is consistent with the Standard Model prediction of $6.6^{+0.7}_{-0.6}$ pb. The differential cross sections in bins of various kinematic variables are presented. The differential event yield as a function of the transverse momentum of the leading $Z$ boson is used to set limits on anomalous neutral triple gauge boson couplings in $ZZ$ production.

A search for heavy resonances decaying to a $W$ or $Z$ boson and a Higgs boson in the $q\bar{q}^{(\prime)}b\bar{b}$ final state is described. The search uses 36.1 fb$^{-1}$ of proton-proton collision data at $\sqrt{s} =$ 13 TeV collected by the ATLAS detector at the CERN Large Hadron Collider in 2015 and 2016. The data are in agreement with the Standard Model expectations, with the largest excess found at a resonance mass of 3.0 TeV with a local (global) significance of 3.3 (2.1) $\sigma$. The results are presented in terms of constraints on a simplified model with a heavy vector triplet. Upper limits are set on the production cross-section times branching ratio for resonances decaying to a $W$ ($Z$) boson and a Higgs boson, itself decaying to $b\bar{b}$, in the mass range between 1.1 and 3.8 TeV; the limits range between 83 and 1.6 fb (77 and 1.1 fb) at 95% confidence level.

This paper presents a study of the production of $WW$ or $WZ$ boson pairs, with one $W$ boson decaying to $e\nu$ or $\mu\nu$ and one $W$ or $Z$ boson decaying hadronically. The analysis uses 20.2 fb$^{-1}$ of $\sqrt{s}=8$ TeV $pp$ collision data, collected by the ATLAS detector at the Large Hadron Collider. Cross-sections for $WW/WZ$ production are measured in high-$p_{T}$ fiducial regions defined close to the experimental event selection. The cross-section is measured for the case where the hadronically decaying boson is reconstructed as two resolved jets, and the case where it is reconstructed as a single jet. The transverse momentum distribution of the hadronically decaying boson is used to search for new physics. Observations are consistent with the Standard Model predictions, and $95\%$ confidence intervals are calculated for parameters describing anomalous triple gauge-boson couplings.

A search is conducted for new resonances decaying into a $W$ or $Z$ boson and a 125 GeV Higgs boson in the $\nu\bar{\nu}b\bar{b}$, $\ell^{\pm}{\nu}b\bar{b}$, and $\ell^+\ell^-b\bar{b}$ final states, where $\ell ^{\pm}= e^{\pm}$ or $\mu^{\pm}$, in $pp$ collisions at $\sqrt s = 13$ TeV. The data used correspond to a total integrated luminosity of 36.1 fb$^{-1}$ collected with the ATLAS detector at the Large Hadron Collider during the 2015 and 2016 data-taking periods. The search is conducted by examining the reconstructed invariant or transverse mass distributions of $Wh$ and $Zh$ candidates for evidence of a localised excess in the mass range of 220 GeV up to 5 TeV. No significant excess is observed and the results are interpreted in terms of constraints on the production cross-section times branching fraction of heavy $W^\prime$ and $Z^\prime$ resonances in heavy-vector-triplet models and the CP-odd scalar boson $A$ in two-Higgs-doublet models. Upper limits are placed at the 95 % confidence level and range between $9.0\times 10^{-4}$ pb and $8.1\times 10^{-1}$ pb depending on the model and mass of the resonance.

A search for events with two displaced vertices from long-lived particles (LLP) pairs using data collected by the ATLAS detector at the LHC is presented. This analysis uses 139~fb$^{-1}$ of proton-proton collision data at $\sqrt{s}=13$ TeV recorded in 2015-2018. The search employs techniques for reconstructing vertices of LLPs decaying to jets in the muon spectrometer displaced between 3 m and 14 m with respect to the primary interaction vertex. The observed numbers of events are consistent with the expected background and limits for several benchmark signals are determined. For the Higgs boson with a mass of 125 GeV, the paper reports the first exclusion limits for branching fractions into neutral long-lived particles below 0.1%, while branching fractions above 10% are excluded at 95% confidence level for LLP proper lifetimes ranging from 4 cm to 72.4 m. In addition, the paper present the first results for the decay of LLPs into into $t\bar{t}$ in the ATLAS muon spectrometer.

A search is presented for the pair production of heavy vector-like $B$ quarks, primarily targeting $B$ quark decays into a $W$ boson and a top quark. The search is based on $36.1$ $fb^{-1}$ of $pp$ collisions at $\sqrt{s}$ = 13 TeV recorded in 2015 and 2016 with the ATLAS detector at the CERN Large Hadron Collider. Data are analysed in the lepton-plus-jets final state, characterised by a high-transverse-momentum isolated electron or muon, large missing transverse momentum, and multiple jets, of which at least one is $b$-tagged. No significant deviation from the Standard Model expectation is observed. The 95% confidence level lower limit on the $B$ mass is 1350 GeV assuming a 100% branching ratio to $Wt$. In the SU(2) singlet scenario, the lower mass limit is 1170 GeV. This search is also sensitive to a heavy vector-like $B$ quark decaying into other final states ($Zb$ and $Hb$) and thus mass limits on $B$ production are set as a function of the decay branching ratios. The 100% branching ratio limits are found to be also applicable to heavy vector-like $X$ production, with charge $+$5/3, that decay into $Wt$.

This paper presents a search for decays of the Higgs boson with a mass of 125 GeV into a pair of new pseudoscalar particles, $H\rightarrow aa$, where one $a$-boson decays into a $b$-quark pair and the other into a muon pair. The search uses 139 fb$^{-1}$ of proton-proton collision data at a center-of-mass energy of $\sqrt{s}=13$ TeV recorded between 2015 and 2018 by the ATLAS experiment at the LHC. A narrow dimuon resonance is searched for in the invariant mass spectrum between 16 GeV and 62 GeV. The largest excess of events above the Standard Model backgrounds is observed at a dimuon invariant mass of 52 GeV and corresponds to a local (global) significance of $3.3 \sigma$ ($1.7 \sigma$). Upper limits at 95% confidence level are placed on the branching ratio of the Higgs boson to the $bb\mu\mu$ final state, $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$, and are in the range $\text{(0.2-4.0)} \times 10^{-4}$, depending on the signal mass hypothesis.

This search, a type not previously performed at ATLAS, uses a comparison of the production cross sections for $e^+ \mu^-$ and $e^- \mu^+$ pairs to constrain physics processes beyond the Standard Model. It uses $139 \text{fb}^{-1}$ of proton$-$proton collision data recorded at $\sqrt{s} = 13$ TeV at the LHC. Targeting sources of new physics which prefer final states containing $e^{+}\mu^{-}$ to $e^{-}\mu^{+}$, the search contains two broad signal regions which are used to provide model-independent constraints on the ratio of cross sections at the 2% level. The search also has two special selections targeting supersymmetric models and leptoquark signatures. Observations using one of these selections are able to exclude, at 95% confidence level, singly produced smuons with masses up to 640 GeV in a model in which the only other light sparticle is a neutralino when the $R$-parity-violating coupling $\lambda'_{231}$ is close to unity. Observations using the other selection exclude scalar leptoquarks with masses below 1880 GeV when $g_{\text{1R}}^{eu}=g_{\text{1R}}^{\mu c}=1$, at 95% confidence level. The limit on the coupling reduces to $g_{\text{1R}}^{eu}=g_{\text{1R}}^{\mu c}=0.46$ for a mass of 1420 GeV.

A search for resonant Higgs boson pair production in the $b\bar{b}b\bar{b}$ final state is presented. The analysis uses 126-139 fb$^{-1}$ of $pp$ collision data at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector at the Large Hadron Collider. The analysis is divided into two channels, targeting Higgs boson decays which are reconstructed as pairs of small-radius jets or as individual large-radius jets. Spin-0 and spin-2 benchmark signal models are considered, both of which correspond to resonant $HH$ production via gluon$-$gluon fusion. The data are consistent with Standard Model predictions. Upper limits are set on the production cross-section times branching ratio to Higgs boson pairs of a new resonance in the mass range from 251 GeV to 5 TeV.

A search for charginos and neutralinos at the Large Hadron Collider is reported using fully hadronic final states and missing transverse momentum. Pair-produced charginos or neutralinos are explored, each decaying into a high-$p_{\text{T}}$ Standard Model weak boson. Fully-hadronic final states are studied to exploit the advantage of the large branching ratio, and the efficient background rejection by identifying the high-$p_{\text{T}}$ bosons using large-radius jets and jet substructure information. An integrated luminosity of 139 fb$^{-1}$ of proton-proton collision data collected by the ATLAS detector at a center-of-mass energy of 13 TeV is used. No significant excess is found beyond the Standard Model expectation. The 95% confidence level exclusion limits are set on wino or higgsino production with varying assumptions in the decay branching ratios and the type of the lightest supersymmetric particle. A wino (higgsino) mass up to 1060 (900) GeV is excluded when the lightest SUSY particle mass is below 400 (240) GeV and the mass splitting is larger than 400 (450) GeV. The sensitivity to high-mass wino and higgsino is significantly extended compared with the previous LHC searches using the other final states.

Results of a search for the electroweak associated production of charginos and next-to-lightest neutralinos, pairs of charginos or pairs of tau sleptons are presented. These processes are characterised by final states with at least two hadronically decaying tau leptons, missing transverse momentum and low jet activity. The analysis is based on an integrated luminosity of 20.3 fb$^{-1}$ of proton--proton collisions at $\sqrt{s}=8$ TeV recorded with the ATLAS experiment at the Large Hadron Collider. No significant excess is observed with respect to the predictions from Standard Model processes. Limits are set at 95% confidence level on the masses of the lighter chargino and next-to-lightest neutralino for various hypotheses for the lightest neutralino mass in simplified models. In the scenario of direct production of chargino pairs, with each chargino decaying into the lightest neutralino via an intermediate tau slepton, chargino masses up to 345 GeV are excluded for a massless lightest neutralino. For associated production of mass-degenerate charginos and next-to-lightest neutralinos, both decaying into the lightest neutralino via an intermediate tau slepton, masses up to 410 GeV are excluded for a massless lightest neutralino.

Measurements of jet activity in top-quark pair events produced in proton--proton collisions are presented, using 3.2 fb$^{-1}$ of $pp$ collision data at a centre-of-mass energy of 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Events are chosen by requiring an opposite-charge $e\mu$ pair and two $b$-tagged jets in the final state. The normalised differential cross-sections of top-quark pair production are presented as functions of additional-jet multiplicity and transverse momentum, $p_{\mathrm T}$. The fraction of signal events that do not contain additional jet activity in a given rapidity region, the gap fraction, is measured as a function of the $p_{\mathrm T}$ threshold for additional jets, and is also presented for different invariant mass regions of the $e\mu b\bar{b}$ system. All measurements are corrected for detector effects and presented as particle-level distributions compared to predictions with different theoretical approaches for QCD radiation. While the kinematics of the jets from top-quark decays are described well, the generators show differing levels of agreement with the measurements of observables that depend on the production of additional jets.

A search for a heavy neutral Higgs boson, $A$, decaying into a $Z$ boson and another heavy Higgs boson, $H$, is performed using a data sample corresponding to an integrated luminosity of 36.1 fb$^{-1}$ from proton-proton collisions at $\sqrt{s} = 13$ TeV recorded in 2015 and 2016 by the ATLAS detector at the Large Hadron Collider. The search considers the $Z$ boson decaying to electrons or muons and the $H$ boson into a pair of $b$-quarks. No evidence for the production of an $A$ boson is found. Considering each production process separately, the 95% confidence-level upper limits on the $pp\rightarrow A\rightarrow ZH$ production cross-section times the branching ratio $H\rightarrow bb$ are in the range of 14-830 fb for the gluon-gluon fusion process and 26-570 fb for the $b$-associated process for the mass ranges 130-700 GeV of the $H$ boson and process for the mass ranges 130-700 GeV of the $H$ boson and 230-800 GeV of the $A$ boson. The results are interpreted in the context of the two-Higgs-doublet model.

A search optimized for new heavy particles decaying to two $b$-quarks and produced in association with additional $b$-quarks is reported. The sensitivity is improved by $b$-tagging at least one lower-$p_{\rm{T}}$ jet in addition to the two highest-$p_{\rm{T}}$ jets. The data used in this search correspond to an integrated luminosity of 103 $\text{fb}^{-1}$ collected with a dedicated trijet trigger during the 2017 and 2018 $\sqrt{s} = 13$ TeV proton-proton collision runs with the ATLAS detector at the LHC. The search looks for resonant peaks in the $b$-tagged dijet invariant mass spectrum over a smoothly falling background. The background is estimated with an innovative data-driven method based on orthonormal functions. The observed $b$-tagged dijet invariant mass spectrum is compatible with the background-only hypothesis. Upper limits at 95% confidence level on a heavy vector-boson production cross section times branching ratio to a pair of $b$-quarks are derived.

The production rates of prompt and non-prompt $J/\psi$ and $\psi(2\mathrm{S})$ mesons are measured using 2.1 $fb^{-1}$ and 11.4 $fb^{-1}$ of data collected with the ATLAS experiment at the LHC, in proton-proton collisions at $\sqrt{s}=7$ and 8 TeV respectively. Production cross-sections for both prompt and non-prompt production sources, ratios of $\psi(2\mathrm{S})$ to $J/\psi$ production, and fractions of non-prompt to inclusive production for $J/\psi$ and $\psi(2\mathrm{S})$ are measured double-differentially as a function of meson $p_{T}$ and rapidity. These measurements are made in a restricted fiducial volume and also corrected for geometrical acceptance after which they are compared to a variety of theoretical predictions.

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

Searches are performed for nonresonant and resonant di-Higgs boson production in the $b\bar{b}\gamma\gamma$ final state. 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. No excess above the expected background is found and upper limits on the di-Higgs boson production cross sections are set. A 95% confidence-level upper limit of 4.2 times the cross section predicted by the Standard Model is set on $pp \rightarrow HH$ nonresonant production, where the expected limit is 5.7 times the Standard Model predicted value. The expected constraints are obtained for a background hypothesis excluding $pp \rightarrow HH$ production. The observed (expected) constraints on the Higgs boson trilinear coupling modifier $\kappa_{\lambda}$ are determined to be $[-1.5, 6.7]$ $([-2.4, 7.7])$ at 95% confidence level, where the expected constraints on $\kappa_{\lambda}$ are obtained excluding $pp \rightarrow HH$ production from the background hypothesis. For resonant production of a new hypothetical scalar particle $X$ ($X \rightarrow HH \rightarrow b\bar{b}\gamma\gamma$), limits on the cross section for $pp \to X \to HH$ are presented in the narrow-width approximation as a function of $m_{X}$ in the range $251 \leq m_{X} \leq 1000$ GeV. The observed (expected) limits on the cross section for $pp \to X \to HH$ range from 640 fb to 44 fb (391 fb to 46 fb) over the considered mass range.

A search for Supersymmetry involving the pair production of gluinos decaying via third-generation squarks to the lightest neutralino is reported. It uses an LHC proton--proton dataset at a center-of-mass energy $\sqrt{s} = 13$ TeV with an integrated luminosity of 3.2 fb$^{-1}$ collected with the ATLAS detector in 2015. The signal is searched for in events containing several energetic jets, of which at least three must be identified as $b$-jets, large missing transverse momentum and, potentially, isolated electrons or muons. Large-radius jets with a high mass are also used to identify highly boosted top quarks. No excess is found above the predicted background. For neutralino masses below approximately 700 GeV, gluino masses of less than 1.78 TeV and 1.76 TeV are excluded at the 95% CL in simplified models of the pair production of gluinos decaying via sbottom and stop, respectively. These results significantly extend the exclusion limits obtained with the $\sqrt{s} = 8$ TeV dataset.

High transverse momentum jets produced in pp collisions at a centre of mass energy of 7 TeV are used to measure the transverse energy-energy correlation function and its associated azimuthal asymmetry. The data were recorded with the ATLAS detector at the LHC in the year 2011 and correspond to an integrated luminosity of 158 $\mathrm{pb}^{-1}$. The selection criteria demand the average transverse momentum of the two leading jets in an event to be larger than 250 GeV. The data at detector level are well described by Monte Carlo event generators. They are unfolded to the particle level and compared with theoretical calculations at next-to-leading-order accuracy. The agreement between data and theory is good and provides a precision test of perturbative Quantum Chromodynamics at large momentum transfers. From this comparison, the strong coupling constant given at the $Z$ boson mass is determined to be $\alpha_{\mathrm{s}}(m_Z) = 0.1173 \pm 0.0010 \mbox{ (exp.) }^{+0.0065}_{-0.0026} \mbox{ (theo.)}$.

A search for new phenomena in final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair, jets, and large missing transverse momentum is presented. This analysis makes use of proton--proton collision data with an integrated luminosity of $36.1 \; \mathrm{fb}^{-1}$, collected during 2015 and 2016 at a centre-of-mass energy $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. The search targets the pair production of supersymmetric coloured particles (squarks or gluinos) and their decays into final states containing an $e^+e^-$ or $\mu^+\mu^-$ pair and the lightest neutralino ($\tilde{\chi}_1^0$) via one of two next-to-lightest neutralino ($\tilde{\chi}_2^0$) decay mechanisms: $\tilde{\chi}_2^0 \rightarrow Z \tilde{\chi}_1^0$, where the $Z$ boson decays leptonically leading to a peak in the dilepton invariant mass distribution around the $Z$ boson mass; and $\tilde{\chi}_2^0 \rightarrow \ell^+\ell^- \tilde{\chi}_1^0$ with no intermediate $\ell^+\ell^-$ resonance, yielding a kinematic endpoint in the dilepton invariant mass spectrum. The data are found to be consistent with the Standard Model expectation. Results are interpreted using simplified models, and exclude gluinos and squarks with masses as large as 1.85 TeV and 1.3 TeV at 95% confidence level, respectively.

This paper reports the results of a search for strong production of supersymmetric particles in 20.1 fb$^{-1}$ of proton-proton collisions at a centre-of-mass energy of 8 TeV using the ATLAS detector at the LHC. The search is performed separately in events with either zero or at least one high $p_\mathrm{T}$ lepton (electron or muon), large missing transverse momentum, high jet multiplicity and at least three jets identified as originated from the fragmentation of a b-quark. No excess is observed with respect to the Standard Model predictions. The results are interpreted in the context of several supersymmetric models involving gluinos and scalar top and bottom quarks, as well as a mSUGRA/CMSSM model. Gluino masses up to 1340 GeV are excluded, depending on the model, significantly extending the previous ATLAS limits.

The production of a $Z$ boson in association with a $J/\psi$ meson in proton-proton collisions probes the production mechanisms of quarkonium and heavy flavour in association with vector bosons, and allows studies of multiple parton scattering. Using $20.3\, \mathrm{fb^{-1}}$ of data collected with the ATLAS experiment at the LHC in $pp$ collisions at $\sqrt{s}=8\, \mathrm{TeV}$, the first measurement of associated $Z + J/\psi$ production is presented for both prompt and non-prompt $J/\psi$ production, with both signatures having a significance in excess of $5\sigma$. The inclusive production cross-sections for $Z$ boson production (analysed in $\mu^+\mu^-$ or $e^+e^-$ decay modes) in association with prompt and non-prompt $J/\psi(\to\mu^+\mu^-)$ are measured relative to the inclusive production rate of $Z$ bosons in the same fiducial volume to be $(36.8\pm 6.7\pm 2.5) \times 10^{-7}$ and $(65.8\pm 9.2\pm 4.2) \times 10^{-7}$ respectively. Normalised differential production cross-section ratios are also determined as a function of the $J/\psi$ transverse momentum. The fraction of signal events arising from single and double parton scattering is estimated, and a lower limit of $5.3\ (3.7)\, \mathrm{mb}$ at $68\ (95)$% confidence level is placed on the effective cross-section regulating double parton interactions.