The ATLAS Collaboration has measured the inclusive production of $Z$ bosons via their decays into electron and muon pairs in $p+$Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV at the Large Hadron Collider. The measurements are made using data corresponding to integrated luminosities of 29.4 nb$^{-1}$ and 28.1 nb$^{-1}$ for $Z \rightarrow ee$ and $Z \rightarrow \mu\mu$, respectively. The results from the two channels are consistent and combined to obtain a cross section times the $Z \rightarrow \ell\ell$ branching ratio, integrated over the rapidity region $|y^{*}_{Z}|<3.5$, of 139.8 $\pm$ 4.8 (stat.) $\pm$ 6.2 (syst.) $\pm$ 3.8 (lumi.) nb. Differential cross sections are presented as functions of the $Z$ boson rapidity and transverse momentum, and compared with models based on parton distributions both with and without nuclear corrections. The centrality dependence of $Z$ boson production in $p+$Pb collisions is measured and analyzed within the framework of a standard Glauber model and the model's extension for fluctuations of the underlying nucleon-nucleon scattering cross section.

Two-particle pseudorapidity correlations are measured in $\sqrt{s_{\rm{NN}}}$ = 2.76 TeV Pb+Pb, $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV $p$+Pb, and $\sqrt{s}$ = 13 TeV $pp$ collisions at the LHC, with total integrated luminosities of approximately 7 $\mu\mathrm{b}^{-1}$, 28 $\mathrm{nb}^{-1}$, and 65 $\mathrm{nb}^{-1}$, respectively. The correlation function $C_{\rm N}(\eta_1,\eta_2)$ is measured as a function of event multiplicity using charged particles in the pseudorapidity range $|\eta|<2.4$. The correlation function contains a significant short-range component, which is estimated and subtracted. After removal of the short-range component, the shape of the correlation function is described approximately by $1+\langle{a_1^2}\rangle \eta_1\eta_2$ in all collision systems over the full multiplicity range. The values of $\sqrt{\langle{a_1^2}\rangle}$ are consistent between the opposite-charge pairs and same-charge pairs, and for the three collision systems at similar multiplicity. The values of $\sqrt{\langle{a_1^2}\rangle}$ and the magnitude of the short-range component both follow a power-law dependence on the event multiplicity. The $\eta$ distribution of the short-range component, after symmetrizing the proton and lead directions in $p$+Pb collisions, is found to be smaller than that in $pp$ collisions with comparable multiplicity.

Measurements of charged-particle fragmentation functions of jets produced in ultra-relativistic nuclear collisions can provide insight into the modification of parton showers in the hot, dense medium created in the collisions. ATLAS has measured jets in $\sqrt{s_{NN}} = 2.76$ TeV Pb+Pb collisions at the LHC using a data set recorded in 2011 with an integrated luminosity of 0.14 nb$^{-1}$. Jets were reconstructed using the anti-$k_{t}$ algorithm with distance parameter values $R$ = 0.2, 0.3, and 0.4. Distributions of charged-particle transverse momentum and longitudinal momentum fraction are reported for seven bins in collision centrality for $R=0.4$ jets with $p_{{T}}^{\mathrm{jet}}> 100$ GeV. Commensurate minimum $p_{\mathrm{T}}$ values are used for the other radii. Ratios of fragment distributions in each centrality bin to those measured in the most peripheral bin are presented. These ratios show a reduction of fragment yield in central collisions relative to peripheral collisions at intermediate $z$ values, $0.04 \lesssim z \lesssim 0.2$ and an enhancement in fragment yield for $z \lesssim 0.04$. A smaller, less significant enhancement is observed at large $z$ and large $p_{\mathrm{T}}$ in central collisions.

Measurements of the centrality and rapidity dependence of inclusive jet production in $\sqrt{s_\mathrm{NN}} = 5.02$ TeV proton--lead ($p$+Pb) collisions and the jet cross-section in $\sqrt{s} = 2.76$ TeV proton--proton collisions are presented. These quantities are measured in datasets corresponding to an integrated luminosity of 27.8 nb$^{-1}$ and 4.0 pb$^{-1}$, respectively, recorded with the ATLAS detector at the Large Hadron Collider in 2013. The $p$+Pb collision centrality was characterised using the total transverse energy measured in the pseudorapidity interval $-4.9 < \eta < -3.2$ in the direction of the lead beam. Results are presented for the double-differential per-collision yields as a function of jet rapidity and transverse momentum ($p_\mathrm{T}$) for minimum-bias and centrality-selected $p$+Pb collisions, and are compared to the jet rate from the geometric expectation. The total jet yield in minimum-bias events is slightly enhanced above the expectation in a $p_\mathrm{T}$-dependent manner but is consistent with the expectation within uncertainties. The ratios of jet spectra from different centrality selections show a strong modification of jet production at all $p_\mathrm{T}$ at forward rapidities and for large $p_\mathrm{T}$ at mid-rapidity, which manifests as a suppression of the jet yield in central events and an enhancement in peripheral events. These effects imply that the factorisation between hard and soft processes is violated at an unexpected level in proton-nucleus collisions. Furthermore, the modifications at forward rapidities are found to be a function of the total jet energy only, implying that the violations may have a simple dependence on the hard parton-parton kinematics.

Measurements of two-particle correlation functions and the first five azimuthal harmonics, $v_1$ to $v_5$, are presented, using 28 $\mathrm{nb}^{-1}$ of $p$+Pb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV measured with the ATLAS detector at the LHC. Significant long-range "ridge-like" correlations are observed for pairs with small relative azimuthal angle ($|\Delta\phi|<\pi/3$) and back-to-back pairs ($|\Delta\phi| > 2\pi/3$) over the transverse momentum range $0.4 < p_{\rm T} < 12$ GeV and in different intervals of event activity. The event activity is defined by either the number of reconstructed tracks or the total transverse energy on the Pb-fragmentation side. The azimuthal structure of such long-range correlations is Fourier decomposed to obtain the harmonics $v_n$ as a function of $p_{\rm T}$ and event activity. The extracted $v_n$ values for $n=2$ to 5 decrease with $n$. The $v_2$ and $v_3$ values are found to be positive in the measured $p_{\rm T}$ range. The $v_1$ is also measured as a function of $p_{\rm T}$ and is observed to change sign around $p_{\rm T}\approx 1.5$-2.0 GeV and then increase to about 0.1 for $p_{\rm T}>4$ GeV. The $v_2(p_{\rm T})$, $v_3(p_{\rm T})$ and $v_4(p_{\rm T})$ are compared to the $v_n$ coefficients in Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}} =2.76$ TeV with similar event multiplicities. Reasonable agreement is observed after accounting for the difference in the average $p_{\rm T}$ of particles produced in the two collision systems.

Correlations between the elliptic or triangular flow coefficients $v_m$ ($m$=2 or 3) and other flow harmonics $v_n$ ($n$=2 to 5) are measured using $\sqrt{s_{NN}}=2.76$ TeV Pb+Pb collision data collected in 2010 by the ATLAS experiment at the LHC, corresponding to an integrated lumonisity of 7 $\mu$b$^{-1}$. The $v_m$-$v_n$ correlations are measured in midrapidity as a function of centrality, and, for events within the same centrality interval, as a function of event ellipticity or triangularity defined in a forward rapidity region. For events within the same centrality interval, $v_3$ is found to be anticorrelated with $v_2$ and this anticorrelation is consistent with similar anticorrelations between the corresponding eccentricities $\epsilon_2$ and $\epsilon_3$. On the other hand, it is observed that $v_4$ increases strongly with $v_2$, and $v_5$ increases strongly with both $v_2$ and $v_3$. The trend and strength of the $v_m$-$v_n$ correlations for $n$=4 and 5 are found to disagree with $\epsilon_m$-$\epsilon_n$ correlations predicted by initial-geometry models. Instead, these correlations are found to be consistent with the combined effects of a linear contribution to $v_n$ and a nonlinear term that is a function of $v_2^2$ or of $v_2v_3$, as predicted by hydrodynamic models. A simple two-component fit is used to separate these two contributions. The extracted linear and nonlinear contributions to $v_4$ and $v_5$ are found to be consistent with previously measured event-plane correlations.

The centrality dependence of the mean charged-particle multiplicity as a function of pseudorapidity is measured in approximately 1 $\mu$b$^{-1}$ of proton--lead collisions at a nucleon--nucleon centre-of-mass energy of $\sqrt{s_{_{\rm{NN}}}} = 5.02$ TeV using the ATLAS detector at the Large Hadron Collider. Charged particles with absolute pseudorapidity less than 2.7 are reconstructed using the ATLAS pixel detector. The $p$+Pb collision centrality is characterised by the total transverse energy measured in the Pb-going direction of the forward calorimeter. The charged-particle pseudorapidity distributions are found to vary strongly with centrality, with an increasing asymmetry between the proton-going and Pb-going directions as the collisions become more central. Three different estimations of the number of nucleons participating in the $p$+Pb collision have been carried out using the Glauber model as well as two Glauber--Gribov inspired extensions to the Glauber model. Charged-particle multiplicities per participant pair are found to vary differently for these three models, highlighting the importance of including colour fluctuations in nucleon--nucleon collisions in the modelling of the initial state of $p$+Pb collisions.

The relationship between jet production in the central region and the underlying-event activity in a pseudorapidity-separated region is studied in 4.0 pb$^{-1}$ of $\sqrt{s} = 2.76$ TeV $pp$ collision data recorded with the ATLAS detector at the LHC. The underlying event is characterised through measurements of the average value of the sum of the transverse energy at large pseudorapidity downstream of one of the protons, which are reported here as a function of hard-scattering kinematic variables. The hard scattering is characterised by the average transverse momentum and pseudorapidity of the two highest transverse momentum jets in the event. The dijet kinematics are used to estimate, on an event-by-event basis, the scaled longitudinal momenta of the hard-scattered partons in the target and projectile beam-protons moving toward and away from the region measuring transverse energy, respectively. Transverse energy production at large pseudorapidity is observed to decrease with a linear dependence on the longitudinal momentum fraction in the target proton and to depend only weakly on that in the projectile proton. The results are compared to the predictions of various Monte Carlo event generators, which qualitatively reproduce the trends observed in data but generally underpredict the overall level of transverse energy at forward pseudorapidity.

A measurement of $\textit{W}$ boson production in lead-lead collisions at $\sqrt{s_{\mathrm{NN}}}=$2.76 TeV is presented. It is based on the analysis of data collected with the ATLAS detector at the LHC in 2011 corresponding to an integrated luminosity of 0.14 $\mathrm{nb}^{-1}$ and 0.15 $\mathrm{nb}^{-1}$ in the muon and electron decay channels, respectively. The differential production yields and lepton charge asymmetry are each measured as a function of the average number of participating nucleons $< N_{\mathrm{part}} >$ and absolute pseudorapidity of the charged lepton. The results are compared to predictions based on next-to-leading-order QCD calculations. These measurements are, in principle, sensitive to possible nuclear modifications to the parton distribution functions and also provide information on scaling of $\textit{W}$ boson production in multi-nucleon systems.

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.