The inclusive production rates of isolated, prompt photons in $p$+Pb collisions at $\sqrt{s_\mathrm{NN}} = 8.16$ TeV are studied with the ATLAS detector at the Large Hadron Collider using a dataset with an integrated luminosity of 165 nb$^{-1}$ recorded in 2016. The cross-section and nuclear modification factor $R_{p\mathrm{Pb}}$ are measured as a function of photon transverse energy from 20 GeV to 550 GeV and in three nucleon-nucleon centre-of-mass pseudorapidity regions, (-2.83,-2.02), (-1.84,0.91), and (1.09,1.90). The cross-section and $R_{p\mathrm{Pb}}$ values are compared with the results of a next-to-leading-order perturbative QCD calculation, with and without nuclear parton distribution function modifications, and with expectations based on a model of the energy loss of partons prior to the hard scattering. The data disfavour a large amount of energy loss and provide new constraints on the parton densities in nuclei.

Multi-particle azimuthal cumulants are measured as a function of centrality and transverse momentum using 470 $\mu$b$^{-1}$ of Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ATLAS detector at the LHC. These cumulants provide information on the event-by-event fluctuations of harmonic flow coefficients $v_n$ and correlated fluctuations between two harmonics $v_n$ and $v_m$. For the first time, a non-zero four-particle cumulant is observed for dipolar flow, $v_1$. The four-particle cumulants for elliptic flow, $v_2$, and triangular flow, $v_3$, exhibit a strong centrality dependence and change sign in ultra-central collisions. This sign change is consistent with significant non-Gaussian fluctuations in $v_2$ and $v_3$. The four-particle cumulant for quadrangular flow, $v_4$, is found to change sign in mid-central collisions. Correlations between two harmonics are studied with three- and four-particle mixed-harmonic cumulants, which indicate an anti-correlation between $v_2$ and $v_3$, and a positive correlation between $v_2$ and $v_4$. These correlations decrease in strength towards central collisions and either approach zero or change sign in ultra-central collisions. To investigate the possible flow fluctuations arising from intrinsic centrality or volume fluctuations, the results are compared between two different event classes used for centrality definitions. In peripheral and mid-central collisions where the cumulant signals are large, only small differences are observed. In ultra-central collisions, the differences are much larger and transverse momentum dependent. These results provide new information to disentangle flow fluctuations from the initial and final states, as well as new insights on the influence of centrality fluctuations.

A search for heavy long-lived multi-charged particles is performed using the ATLAS detector at the LHC. Data with an integrated luminosity of 36.1 fb$^{-1}$ collected in 2015 and 2016 from proton-proton collisions at $\sqrt{s}$ = 13 TeV are examined. Particles producing anomalously high ionization, consistent with long-lived massive particles with electric charges from |q|=2e to |q|=7e, are searched for. No events are observed, and 95% confidence level cross-section upper limits are interpreted as lower mass limits for a Drell-Yan production model. Multi-charged particles with masses between 50 GeV and 980-1220 GeV (depending on their electric charge) are excluded.

A search for excited electrons produced in $pp$ collisions at $\sqrt{s} = 13$ TeV via a contact interaction $q\bar{q} \to ee^*$ is presented. The search uses 36.1 fb$^{-1}$ of data collected in 2015 and 2016 by the ATLAS experiment at the Large Hadron Collider. Decays of the excited electron via a contact interaction into an electron and a pair of quarks ($eq\bar{q}$) are targeted in final states with two electrons and two hadronic jets, and decays via a gauge interaction into a neutrino and a $W$ boson ($\nu W$) are probed in final states with an electron, missing transverse momentum, and a large-radius jet consistent with a hadronically decaying $W$ boson. No significant excess is observed over the expected backgrounds. Upper limits are calculated for the $pp \to ee^* \to eeq\bar{q}$ and $pp \to ee^* \to e\nu W$ production cross sections as a function of the excited electron mass $m_{e^*}$ at 95% confidence level. The limits are translated into lower bounds on the compositeness scale parameter $\Lambda$ of the model as a function of $m_{e^*}$. For $m_{e^*} < 0.5$ TeV, the lower bound for $\Lambda$ is 11 TeV. In the special case of $m_{e^*} = \Lambda$, the values of $m_{e^*} < 4.8$ TeV are excluded. The presented limits on $\Lambda$ are more stringent than those obtained in previous searches.

The fragmentation of high-energy gluons at small opening angles is largely unconstrained by present measurements. Gluon splitting to $b$-quark pairs is a unique probe into the properties of gluon fragmentation because identified $b$-tagged jets provide a proxy for the quark daughters of the initial gluon. In this study, key differential distributions related to the $g\rightarrow b\bar{b}$ process are measured using 33 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data recorded by the ATLAS experiment at the LHC in 2016. Jets constructed from charged-particle tracks, clustered with the anti-$k_t$ jet algorithm with radius parameter $R = 0.2$, are used to probe angular scales below the $R=0.4$ jet radius. The observables are unfolded to particle level in order to facilitate direct comparisons with predictions from present and future simulations. Multiple significant differences are observed between the data and parton shower Monte Carlo predictions, providing input to improve these predictions of the main source of background events in analyses involving boosted Higgs bosons decaying into $b$-quarks.

This paper presents a measurement of forward-forward and forward-central dijet azimuthal angular correlations and conditional yields in proton-proton ($pp$) and proton-lead ($p$+Pb) collisions as a probe of the nuclear gluon density in regions where the fraction of the average momentum per nucleon carried by the parton entering the hard scattering is low. In these regions, gluon saturation can modify the rapidly increasing parton distribution function of the gluon. The analysis utilizes 25 pb$^{-1}$ of $pp$ data and 360 $\mu \mathrm{b}^{-1}$ of $p$+Pb data, both at $\sqrt{s_{\rm NN}}$ = 5.02 TeV, collected in 2015 and 2016, respectively, with the ATLAS detector at the LHC. The measurement is performed in the center-of-mass frame of the nucleon-nucleon system in the rapidity range between $-$4.0 and 4.0 using the two highest transverse momentum jets in each event, with the highest transverse momentum jet restricted to the forward rapidity range. No significant broadening of azimuthal angular correlations is observed for forward-forward or forward-central dijets in $p$+Pb compared to $pp$ collisions. For forward-forward jet pairs in the proton-going direction, the ratio of conditional yields in $p$+Pb collisions to those in $pp$ collisions is suppressed by approximately 20%, with no significant dependence on the transverse momentum of the dijet system. No modification of conditional yields is observed for forward-central dijets.

This paper describes a search for pairs of neutral, long-lived particles decaying in the ATLAS calorimeter. Long-lived particles occur in many extensions to the Standard Model and may elude searches for new promptly decaying particles. The analysis considers neutral, long-lived scalars with masses between 5 GeV and 400 GeV, produced from decays of heavy bosons with masses between 125 GeV and 1000 GeV, where the long-lived scalars decay into Standard Model fermions. The analysis uses either 10.8 fb$^{-1}$ or 33.0 fb$^{-1}$ of data (depending on the trigger) recorded in 2016 at the LHC with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of 13 TeV. No significant excess is observed, and limits are reported on the production cross section times branching ratio as a function of the proper decay length of the long-lived particles.

This paper presents measurements of the $W^+ \rightarrow \mu^+\nu$ and $W^- \rightarrow \mu^-\nu$ cross-sections and the associated charge asymmetry as a function of the absolute pseudorapidity of the decay muon. The data were collected in proton--proton collisions at a centre-of-mass energy of 8 TeV with the ATLAS experiment at the LHC and correspond to a total integrated luminosity of $20.2~\mbox{fb$^{-1}$}$. The precision of the cross-section measurements varies between 0.8% to 1.5% as a function of the pseudorapidity, excluding the 1.9% uncertainty on the integrated luminosity. The charge asymmetry is measured with an uncertainty between 0.002 and 0.003. The results are compared with predictions based on next-to-next-to-leading-order calculations with various parton distribution functions and have the sensitivity to discriminate between them.

The ratio of the cross sections for inclusive isolated-photon production in $pp$ collisions at centre-of-mass energies of 13 and 8 TeV is measured using the ATLAS detector at the LHC. The integrated luminosities of the 13 TeV and 8 TeV datasets are 3.2 fb$^{-1}$ and 20.2 fb$^{-1}$, respectively. The ratio is measured as a function of the photon transverse energy in different regions of the photon pseudorapidity. The predictions from next-to-leading-order perturbative QCD calculations are compared with the measured ratio. The experimental systematic uncertainties as well as the uncertainties affecting the predictions are evaluated taking into account the correlations between the two centre-of-mass energies, resulting in a reduction of up to a factor of $2.5$ ($5$) in the experimental (theoretical) systematic uncertainties. The predictions based on several parameterisations of the proton parton distribution functions agree with the data within the reduced experimental and theoretical uncertainties. In addition, this ratio to that of the fiducial cross sections for $Z$ boson production at 13 and 8 TeV using the decay channels $Z \rightarrow e^+e^-$ and $Z \rightarrow \mu^+\mu^-$ is made and compared with the theoretical predictions. In this double ratio, a further reduction of the experimental uncertainty is obtained because the uncertainties arising from the luminosity measurement cancel out. The predictions describe the measurements of the double ratio within the theoretical and experimental uncertainties.

A measurement of jet substructure variables is presented using data collected in 2016 by the ATLAS experiment at the LHC with proton-proton collisions at $\sqrt{s}=13$ TeV. Large-radius jets groomed with the trimming and soft-drop algorithms are studied. Dedicated event selections are used to study jets produced by light quarks or gluons, and hadronically decaying top quarks and $W$ bosons. The variables measured are sensitive to pronged substructure, and therefore are typically used for tagging jets from boosted massive particles. These include the energy correlation functions and the $N$-subjettiness variables. The number of subjets and the Les Houches angularity are also considered. The distributions of the substructure variables, corrected for detector effects, are compared to the predictions of various Monte Carlo event generators. They are also compared between the large-radius jets originating from light quarks or gluons, and hadronically decaying top quarks and $W$ bosons.

A measurement of the associated production of a top-quark pair ($t\bar{t}$) with a vector boson ($W$, $Z$) in proton-proton collisions at a center-of-mass energy of 13 TeV is presented, using $36.1$ fb$^{-1}$ of integrated luminosity collected by the ATLAS detector at the Large Hadron Collider. Events are selected in channels with two same- or opposite-sign leptons (electrons or muons), three leptons or four leptons, and each channel is further divided into multiple regions to maximize the sensitivity of the measurement. The $t\bar{t}Z$ and $t\bar{t}W$ production cross sections are simultaneously measured using a combined fit to all regions. The best-fit values of the production cross sections are $\sigma_{t\bar{t}Z} = 0.95 \pm 0.08_{\mathrm{stat.}} \pm 0.10_{\mathrm{syst.}}$ pb and $\sigma_{t\bar{t}W} = 0.87 \pm 0.13_{\mathrm{stat.}} \pm 0.14_{\mathrm{syst.}}$ pb in agreement with the Standard Model predictions. The measurement of the $t\bar{t}Z$ cross section is used to set constraints on effective field theory operators which modify the $t\bar{t}Z$ vertex.

A search for new particles decaying into a pair of top quarks is performed using proton-proton collision data recorded with the ATLAS detector at the Large Hadron Collider at a center-of-mass energy of $\sqrt{s} = $13 TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$. Events consistent with top-quark pair production and the fully hadronic decay mode of the top quarks are selected by requiring multiple high transverse momentum jets including those containing $b$-hadrons. Two analysis techniques, exploiting dedicated top-quark pair reconstruction in different kinematic regimes, are used to optimize the search sensitivity to new hypothetical particles over a wide mass range. The invariant mass distribution of the two reconstructed top-quark candidates is examined for resonant production of new particles with various spins and decay widths. No significant deviation from the Standard Model prediction is observed and limits are set on the production cross-section times branching fraction for new hypothetical $Z'$ bosons, dark-matter mediators, Kaluza-Klein gravitons and Kaluza-Klein gluons. By comparing with the predicted production cross-sections, the $Z'$ boson in the topcolor-assisted-technicolor model is excluded for masses up to 3.1$-$3.6 TeV, the dark-matter mediators in a simplified framework are excluded in the mass ranges from 0.8 TeV to 0.9 TeV and from 2.0 TeV to 2.2 TeV, and the Kaluza-Klein gluon is excluded for masses up to 3.4 TeV, depending on the decay widths of the particles.

A search is performed for localised excesses in dijet mass distributions of low-dijet-mass events produced in association with a high transverse energy photon. The search uses up to 79.8 fb$^{-1}$ of LHC proton-proton collisions collected by the ATLAS experiment at a centre-of-mass energy of 13 TeV during 2015-2017. Two variants are presented: one which makes no jet flavour requirements and one which requires both jets to be tagged as $b$-jets. The observed mass distributions are consistent with multi-jet processes in the Standard Model. The data are used to set upper limits on the production cross-section for a benchmark $Z^\prime$ model and, separately, on generic Gaussian-shape contributions to the mass distributions, extending the current ATLAS constraints on dijet resonances to the mass range between 225 and 1100 GeV.

A search for a narrow scalar resonance decaying into an opposite-sign muon pair produced in events with and without $b$-tagged jets is presented in this paper. The search uses 36.1 fb$^{-1}$ of $\sqrt{s}= 13$ TeV proton-proton collision data recorded by the ATLAS experiment at the LHC. No significant excess of events above the expected Standard Model background is observed in the investigated mass range of 0.2 to 1.0 TeV. The observed upper limits at 95$\%$ confidence level on the cross section times branching ratio for $b$-quark associated production and gluon-gluon fusion are between 1.9 and 41 fb and 1.6 and 44 fb respectively, which is consistent with expectations.

A search for a right-handed gauge boson $W_{\mathrm{R}}$, decaying into a boosted right-handed heavy neutrino $N_{\mathrm{R}}$, in the framework of Left-Right Symmetric Models is presented. It is based on data from proton-proton collisions with a centre-of-mass energy of 13 TeV collected by the ATLAS detector at the Large Hadron Collider during the years 2015, 2016 and 2017, corresponding to an integrated luminosity of 80 fb$^{-1}$. The search is performed separately for electrons and muons in the final state. A distinguishing feature of the search is the use of large-radius jets containing electrons. Selections based on the signal topology result in smaller background compared with to expected signal. No significant deviation from the Standard Model prediction is observed and lower limits are set in the $W_{\mathrm{R}}$ and $N_{\mathrm{R}}$ mass plane. Mass values of the $W_{\mathrm{R}}$ smaller than 3.8-5 TeV are excluded for $N_{\mathrm{R}}$ in the mass range 0.1-1.8 TeV.

This paper reports on a search for the electroweak diboson ($WW/WZ/ZZ$) production in association with a high-mass dijet system, using data from proton-proton collisions at a center-of-mass energy of $\sqrt{s}=13$ TeV. The data, corresponding to an integrated luminosity of 35.5 fb$^{-1}$, were recorded with the ATLAS detector in 2015 and 2016 at the Large Hadron Collider. The search is performed in final states in which one boson decays leptonically, and the other boson decays hadronically. The hadronically decaying $W/Z$ boson is reconstructed as either two small-radius jets or one large-radius jet using jet substructure techniques. The electroweak production of $WW/WZ/ZZ$ in association with two jets is measured with an observed (expected) significance of 2.7 (2.5) standard deviations, and the fiducial cross section is measured to be $45.1 \pm 8.6(\mathrm{stat.}) ^{+15.9} _{-14.6} (\mathrm{syst.})$ fb.

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.

This paper presents studies of Bose-Einstein correlations (BEC) in proton-proton collisions at a centre-of-mass energy of 13 TeV, using data from the ATLAS detector at the CERN Large Hadron Collider. Data were collected in a special low-luminosity configuration with a minimum-bias trigger and a high-multiplicity track trigger, accumulating integrated luminosities of 151 $\mu$b$^{-1}$ and 8.4 nb$^{-1}$ respectively. The BEC are measured for pairs of like-sign charged particles, each with $|\eta|$ < 2.5, for two kinematic ranges: the first with particle $p_T$ > 100 MeV and the second with particle $p_T$ > 500 MeV. The BEC parameters, characterizing the source radius and particle correlation strength, are investigated as functions of charged-particle multiplicity (up to 300) and average transverse momentum of the pair (up to 1.5 GeV). The double-differential dependence on charged-particle multiplicity and average transverse momentum of the pair is also studied. The BEC radius is found to be independent of the charged-particle multiplicity for high charged-particle multiplicity (above 100), confirming a previous observation at lower energy. This saturation occurs independent of the transverse momentum of the pair.

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.

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$.

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.

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.