Measurements of prompt $\psi$(2S) meson production cross sections in proton-lead (pPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_{_\mathrm{NN}}}=$ 5.02 TeV are reported. The results are based on pPb and pp data collected by the CMS experiment at the LHC, corresponding to integrated luminosities of 34.6 nb$^{-1}$ and 28.0 pb$^{-1}$, respectively. The nuclear modification factor $R_\mathrm{pPb}$ is measured for prompt $\psi$(2S) in the transverse momentum range 4 $<$ p$_\mathrm{T}$ $<$ 30 GeV$/c$ and the center-of-mass rapidity range $-$2.4 $< y_\mathrm{cm} <$ 1.93. The results on $\psi$(2S) $R_\mathrm{pPb}$ are compared to the corresponding modification factor for prompt J$/\psi$ mesons and are found to be more suppressed than the J$/\psi$ states over the entire kinematic range studied.

A measurement is presented of the charged hadron multiplicity in hadronic PbPb collisions, as a function of pseudorapidity and centrality, at a collision energy of 2.76 TeV per nucleon pair. The data sample is collected using the CMS detector and a minimum-bias trigger, with the CMS solenoid off. The number of charged hadrons is measured both by counting the number of reconstructed particle hits and by forming hit doublets of pairs of layers in the pixel detector. The two methods give consistent results. The charged hadron multiplicity density dN(ch)/d eta, evaluated at eta=0 for head-on collisions, is found to be 1612 +/- 55, where the uncertainty is dominated by systematic effects. Comparisons of these results to previous measurements and to various models are also presented.

Evidence for the light-by-light scattering process, $\gamma\gamma$ $\to$ $\gamma\gamma$, in ultraperipheral PbPb collisions at a centre-of-mass energy per nucleon pair of 5.02 TeV is reported. The analysis is conducted using a data sample corresponding to an integrated luminosity of 390 $\mu$b$^{-1}$ recorded by the CMS experiment at the LHC. Light-by-light scattering processes are selected in events with two photons exclusively produced, each with transverse energy E$_\mathrm{T}^{\gamma}$ $>$ 2 GeV, pseudorapidity $|\eta^{\gamma}|$ $\lt$ 2.4, diphoton invariant mass $m^{\gamma\gamma}$ $\gt$ 5 GeV, diphoton transverse momentum $p_\mathrm{T}^{\gamma\gamma}$ $\lt$ 1 GeV, and diphoton acoplanarity below 0.01. After all selection criteria are applied, 14 events are observed, compared to expectations of 9.0 $\pm$ 0.9 (theo) events for the signal and 4.0 $\pm$ 1.2 (stat) for the background processes. The excess observed in data relative to the background-only expectation corresponds to a significance of 3.7 standard deviations, and has properties consistent with those expected for the light-by-light scattering signal. The measured fiducial light-by-light scattering cross section, $\sigma_\mathrm{fid} (\gamma\gamma$ $\to$ $\gamma\gamma) =$ 120 $\pm$ 46 (stat) $\pm$ 28 (syst) $\pm$ 12 (theo) nb, is consistent with the standard model prediction. The $m^{\gamma\gamma}$ distribution is used to set new exclusion limits on the production of pseudoscalar axion-like particles, via the $\gamma\gamma$ $\to$ a $\to$ $\gamma\gamma$ process, in the mass range $m_{\mathrm{a}} =$ 5-90 GeV.

Measurements of charged-particle production in pp, p$-$Pb, and Pb$-$Pb collisions in the toward, away, and transverse regions with the ALICE detector are discussed. These regions are defined event-by-event relative to the azimuthal direction of the charged trigger particle, which is the reconstructed particle with the largest transverse momentum ($p_{\mathrm{T}}^{\rm trig}$) in the range $8<p_{\mathrm{T}}^{\rm trig}<15$ GeV$/c$. The toward and away regions contain the primary and recoil jets, respectively; both regions are accompanied by the underlying event (UE). In contrast, the transverse region perpendicular to the direction of the trigger particle is dominated by the so-called UE dynamics, and includes also contributions from initial- and final-state radiation. The relative transverse activity classifier, $R_{\mathrm{T}}=N_{\mathrm{ch}}^{\mathrm{T}}/\langle N_{\mathrm{ch}}^{\mathrm{T}}\rangle$, is used to group events according to their UE activity, where $N_{\mathrm{ch}}^{\mathrm{T}}$ is the charged-particle multiplicity per event in the transverse region and $\langle N_{\mathrm{ch}}^{\mathrm{T}}\rangle$ is the mean value over the whole analysed sample. The energy dependence of the $R_{\mathrm{T}}$ distributions in pp collisions at $\sqrt{s}=2.76$, 5.02, 7, and 13 TeV is reported, exploring the Koba-Nielsen-Olesen (KNO) scaling properties of the multiplicity distributions. The first measurements of charged-particle $p_{\rm T}$ spectra as a function of $R_{\mathrm{T}}$ in the three azimuthal regions in pp, p$-$Pb, and Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV are also reported. Data are compared with predictions obtained from the event generators PYTHIA 8 and EPOS LHC. This set of measurements is expected to contribute to the understanding of the origin of collective-like effects in small collision systems (pp and p$-$Pb).

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\sqrt{s_\mathrm{NN}} =$ 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

The cross sections for $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) production in lead-lead (PbPb) and proton-proton (pp) collisions at $\sqrt{s_{_\mathrm{NN}}} =$ 5.02 TeV have been measured using the CMS detector at the LHC. The nuclear modification factors, R$_\mathrm{AA}$, derived from the PbPb-to-pp ratio of yields for each state, are studied as functions of meson rapidity and transverse momentum, as well as PbPb collision centrality. The yields of all three states are found to be significantly suppressed, and compatible with a sequential ordering of the suppression, R$_\mathrm{AA}$($\Upsilon$(1S)) $>$ R$_\mathrm{AA}$($\Upsilon$(2S)) $>$ R$_\mathrm{AA}$($\Upsilon$(3S)) . The suppression of $\Upsilon$(1S) is larger than that seen at $\sqrt{s_{_\mathrm{NN}}} =$ 2.76 TeV, although the two are compatible within uncertainties. The upper limit on the R$_\mathrm{AA}$ of $\Upsilon$(3S) integrated over $p_\mathrm{T}$ and rapidity is 0.094 at 95% confidence level, which is the strongest suppression observed for any hadron species in heavy ion collisions to date.

The B$^0_s$ and B$^+$ production yields are measured in PbPb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data sample, collected with the CMS detector at the LHC, corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The mesons are reconstructed in the exclusive decay channels B$^0_s$$\to$ J/$\psi(\mu^+\mu^-)\phi($K$^+$K$^-)$ and B$^+$$\to$ J/$\psi(\mu^+\mu^-)$K$^+$, in the transverse momentum range 7-50 GeV/c and absolute rapidity 0-2.4. The B$^0_s$ meson is observed with a statistical significance in excess of five standard deviations for the first time in nucleus-nucleus collisions. The measurements are performed as functions of the transverse momentum of the B mesons and of the PbPb collision centrality. The ratio of production yields of B$^0_s$ and B$^+$ is measured and compared to theoretical models that include quark recombination effects.

The transverse-momentum ($p_{\rm T}$) spectra of ${\rm K}^*(892)^0$ and $\phi(1020)$ measured with the ALICE detector up to $p_{\rm T}$ = 16 GeV/$c$ in the rapidity range $-1.2 < y < 0.3$, in p-Pb collisions at the center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured $p_{\rm T}$ distributions show a dependence on both multiplicity and rapidity at low $p_{\rm T}$ whereas no significant dependence is observed at high $p_{\rm T}$. A rapidity dependence is observed in the $p_{\rm T}$-integrated yield (d$N$/d$y$), whereas the mean transverse momentum ($\left< p_{\rm T} \right>$) shows a flat behavior as a function of rapidity. The rapidity asymmetry ($Y_{\rm asym}$) at low $p_{\rm T}$ ( < 5 GeV/$c$) is more significant for higher multiplicity classes. At high $p_{\rm T}$, no significant rapidity asymmetry is observed in any of the multiplicity classes. Both ${\rm K}^*(892)^0$ and $\phi(1020)$ show similar $Y_{\rm asym}$. The nuclear modification factor ($Q_{\rm CP}$) as a function of $p_{\rm T}$ shows a Cronin-like enhancement at intermediate $p_{\rm T}$, which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high $p_{\rm T}$ (> 5 GeV/$c$), the $Q_{\rm CP}$ values are greater than unity and no significant rapidity dependence is observed.

The azimuthal ($\Delta\varphi$) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p$-$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV. Results are reported for electrons with transverse momentum $4<p_{\rm T}<16$ GeV/$c$ and pseudorapidity $|\eta|<0.6$. The associated charged particles are selected with transverse momentum $1<p_{\rm T}<7$ GeV/$c$, and relative pseudorapidity separation with the leading electron $|\Delta\eta| < 1$. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p$-$Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The $\Delta\varphi$ distribution and the peak observables in pp and p$-$Pb collisions are compared with calculations from various Monte Carlo event generators.

The production of prompt charged particles in proton-lead collisions and in proton-proton collisions at the nucleon-nucleon centre-of-mass energy ${\sqrt{s_{\scriptscriptstyle\mathrm{NN}}}=5\,\mathrm{TeV}}$ is studied at LHCb as a function of pseudorapidity ($\eta$) and transverse momentum ($p_{\mathrm{T}}$) with respect to the proton beam direction. The nuclear modification factor for charged particles is determined as a function of $\eta$ between ${-4.8<\eta<-2.5}$ (backward region) and ${2.0<\eta<4.8}$ (forward region), and $p_{\mathrm{T}}$ between ${0.2<p_{\mathrm{T}}<8.0\,\mathrm{GeV}/c}$. The results show a suppression of charged particle production in proton-lead collisions relative to proton-proton collisions in the forward region and an enhancement in the backward region for $p_{\mathrm{T}}$ larger than $1.5\,\mathrm{GeV}/c$. This measurement constrains nuclear PDFs and saturation models at previously unexplored values of the parton momentum fraction down to $10^{-6}$.