Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

We have studied the dependence of azimuthal anisotropy $v_2$ for inclusive and identified charged hadrons in Au$+$Au and Cu$+$Cu collisions on collision energy, species, and centrality. The values of $v_2$ as a function of transverse momentum $p_T$ and centrality in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV and 62.4 GeV are the same within uncertainties. However, in Cu$+$Cu collisions we observe a decrease in $v_2$ values as the collision energy is reduced from 200 to 62.4 GeV. The decrease is larger in the more peripheral collisions. By examining both Au$+$Au and Cu$+$Cu collisions we find that $v_2$ depends both on eccentricity and the number of participants, $N_{\rm part}$. We observe that $v_2$ divided by eccentricity ($\varepsilon$) monotonically increases with $N_{\rm part}$ and scales as ${N_{\rm part}^{1/3}}$. The Cu$+$Cu data at 62.4 GeV falls below the other scaled $v_{2}$ data. For identified hadrons, $v_2$ divided by the number of constituent quarks $n_q$ is independent of hadron species as a function of transverse kinetic energy $KE_T=m_T-m$ between $0.1<KE_T/n_q<1$ GeV. Combining all of the above scaling and normalizations, we observe a near-universal scaling, with the exception of the Cu$+$Cu data at 62.4 GeV, of $v_2/(n_q\cdot\varepsilon\cdot N^{1/3}_{\rm part})$ vs $KE_T/n_q$ for all measured particles.

Constraints on the Higgs boson self-coupling are set by combining double-Higgs boson analyses in the $b\bar{b}b\bar{b}$, $b\bar{b}\tau^+\tau^-$ and $b\bar{b} \gamma \gamma$ decay channels with single-Higgs boson analyses targeting the $\gamma \gamma$, $ZZ^*$, $WW^*$, $\tau^+ \tau^-$ and $b\bar{b}$ decay channels. The data used in these analyses were recorded by the ATLAS detector at the LHC in proton$-$proton collisions at $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 126$-$139 fb$^{-1}$. The combination of the double-Higgs analyses sets an upper limit of $\mu_{HH} < 2.4$ at 95% confidence level on the double-Higgs production cross-section normalised to its Standard Model prediction. Combining the single-Higgs and double-Higgs analyses, with the assumption that new physics affects only the Higgs boson self-coupling ($\lambda_{HHH}$), values outside the interval $-0.4< \kappa_{\lambda}=(\lambda_{HHH}/\lambda_{HHH}^{\textrm{SM}})< 6.3$ are excluded at 95% confidence level. The combined single-Higgs and double-Higgs analyses provide results with fewer assumptions, by adding in the fit more coupling modifiers introduced to account for the Higgs boson interactions with the other Standard Model particles. In this relaxed scenario, the constraint becomes $-1.4 < \kappa_{\lambda} < 6.1$ at 95% CL.

The nuclear modification factors of J/$\psi$ and $\psi$(2S) mesons are measured in PbPb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC in 2015, corresponding to integrated luminosities of 464 $\mu$b$^{-1}$ and 28 pb$^{-1}$, respectively. The measurements are performed in the dimuon rapidity range of $|y| <$ 2.4 as a function of centrality, rapidity, and transverse momentum (p$_\mathrm{T}$) from p$_\mathrm{T}=$ 3 GeV/$c$ in the most forward region and up to 50 GeV/$c$. Both prompt and nonprompt (coming from b hadron decays) mesons are observed to be increasingly suppressed with centrality, with a magnitude similar to the one observed at $\sqrt{s_{\mathrm{NN}}}=$ 2.76 TeV for the two J/$\psi$ meson components. No dependence on rapidity is observed for either prompt or nonprompt J/$\psi$ mesons. An indication of a lower prompt J/$\psi$ meson suppression at p$_\mathrm{T} >$ 25 GeV/$c$ is seen with respect to that observed at intermediate p$_\mathrm{T}$. The prompt $\psi$(2S) meson yield is found to be more suppressed than that of the prompt J/$\psi$ mesons in the entire p$_\mathrm{T}$ range.

Measurements of the $\pi^{\pm}$, $K^{\pm}$, and proton double differential yields emitted from the surface of the 90-cm-long carbon target (T2K replica) were performed for the incoming 31 GeV/c protons with the NA61/SHINE spectrometer at the CERN SPS using data collected during 2010 run. The double differential $\pi^{\pm}$ yields were measured with increased precision compared to the previously published NA61/SHINE results, while the $K^{\pm}$ and proton yields were obtained for the first time. A strategy for dealing with the dependence of the results on the incoming proton beam profile is proposed. The purpose of these measurements is to reduce significantly the (anti)neutrino flux uncertainty in the T2K long-baseline neutrino experiment by constraining the production of (anti)neutrino ancestors coming from the T2K target.

Measurements of particle emission from a replica of the T2K 90 cm-long carbon target were performed in the NA61/SHINE experiment at CERN SPS, using data collected during a high-statistics run in 2009. An efficient use of the long-target measurements for neutrino flux predictions in T2K requires dedicated reconstruction and analysis techniques. Fully-corrected differential yields of $\pi^\pm$-mesons from the surface of the T2K replica target for incoming 31 GeV/c protons are presented. A possible strategy to implement these results into the T2K neutrino beam predictions is discussed and the propagation of the uncertainties of these results to the final neutrino flux is performed.

Studies of the fragmentation of jets into charged particles in heavy-ion collisions can provide information about the mechanism of jet-quenching by the hot and dense QCD matter created in such collisions, the quark-gluon plasma. This paper presents a measurement of the angular distribution of charged particles around the jet axis in $\sqrt{s_{\mathrm{NN}}}=$ 5.02 TeV Pb+Pb and $pp$ collisions, using the ATLAS detector at the LHC. The Pb+Pb and $pp$ data sets have integrated luminosities of 0.49 nb$^{-1}$ and 25 pb$^{-1}$, respectively. The measurement is performed for jets reconstructed with the anti-$k_{t}$ algorithm with radius parameter $R = 0.4$ and is extended to an angular distance of $r= 0.8$ from the jet axis. Results are presented as a function of Pb+Pb collision centrality and distance from the jet axis for charged particles with transverse momenta in the 1$-$63 GeV range, matched to jets with transverse momenta in the 126$-$316 GeV range and an absolute value of jet rapidity of less than 1.7. Modifications to the measured distributions are quantified by taking a ratio to the measurements in $pp$ collisions. Yields of charged particles with transverse momenta below 4 GeV are observed to be increasingly enhanced as a function of angular distance from the jet axis, reaching a maximum at $r=0.6$. Charged particles with transverse momenta above 4 GeV have an enhanced yield in Pb+Pb collisions in the jet core for angular distances up to $r = 0.05$ from the jet axis, with a suppression at larger distances.

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 azimuthal anisotropy of charged particles produced in $\sqrt{s_{\mathrm{NN}}}=8.16$ TeV $p$+Pb collisions is measured with the ATLAS detector at the LHC. The data correspond to an integrated luminosity of $165$ $\mathrm{nb}^{-1}$ that was collected in 2016. Azimuthal anisotropy coefficients, elliptic $v_2$ and triangular $v_3$, extracted using two-particle correlations with a non-flow template fit procedure, are presented as a function of particle transverse momentum ($p_\mathrm{T}$) between $0.5$ and $50$ GeV. The $v_2$ results are also reported as a function of centrality in three different particle $p_\mathrm{T}$ intervals. The results are reported from minimum-bias events and jet-triggered events, where two jet $p_\mathrm{T}$ thresholds are used. The anisotropies for particles with $p_\mathrm{T}$ less than about $2$ GeV are consistent with hydrodynamic flow expectations, while the significant non-zero anisotropies for $p_\mathrm{T}$ in the range $9$-$50$ GeV are not explained within current theoretical frameworks. In the $p_\mathrm{T}$ range $2$-$9$ GeV, the anisotropies are larger in minimum-bias than in jet-triggered events. Possible origins of these effects, such as the changing admixture of particles from hard scattering and the underlying event, are discussed.

Correlations of two flow harmonics $v_n$ and $v_m$ via three- and four-particle cumulants are measured in 13 TeV $pp$, 5.02 TeV $p$+Pb, and 2.76 TeV peripheral Pb+Pb collisions with the ATLAS detector at the LHC. The goal is to understand the multi-particle nature of the long-range collective phenomenon in these collision systems. The large non-flow background from dijet production present in the standard cumulant method is suppressed using a method of subevent cumulants involving two, three and four subevents separated in pseudorapidity. The results show a negative correlation between $v_2$ and $v_3$ and a positive correlation between $v_2$ and $v_4$ for all collision systems and over the full multiplicity range. However, the magnitudes of the correlations are found to depend strongly on the event multiplicity, the choice of transverse momentum range and collision system. The relative correlation strength, obtained by normalisation of the cumulants with the $\langle v_n^2\rangle$ from a two-particle correlation analysis, is similar in the three collision systems and depends weakly on the event multiplicity and transverse momentum. These results based on the subevent methods provide strong evidence of a similar long-range multi-particle collectivity in $pp$, $p$+Pb and peripheral Pb+Pb collisions.