The second-order ($v_2$) and third-order ($v_3$) Fourier coefficients describing the azimuthal anisotropy of prompt and nonprompt (from b-hadron decays) J/$\psi$, as well as prompt $\psi$(2S) mesons are measured in lead-lead collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The analysis uses a data set corresponding to an integrated luminosity of 1.61 nb$^{-1}$ recorded with the CMS detector. The J/$\psi$ and $\psi$(2S) mesons are reconstructed using their dimuon decay channel. The $v_2$ and $v_3$ coefficients are extracted using the scalar product method and studied as functions of meson transverse momentum and collision centrality. The measured $v_2$ values for prompt J/$\psi$ mesons are found to be larger than those for nonprompt J/$\psi$ mesons. The prompt J/$\psi$$v_2$ values at high $p_\mathrm{T}$ are found to be underpredicted by a model incorporating only parton energy loss effects in a quark-gluon plasma medium. Prompt and nonprompt J/$\psi$ meson $v_3$ and prompt $\psi$(2S) $v_2$ and $v_3$ values are also reported for the first time, providing new information about heavy quark interactions in the hot and dense medium created in heavy ion collisions.

A search for resonances in events with the $\gamma$+jet final state has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC. The total data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$. Models of excited quarks and quantum black holes are considered. Using a wide-jet reconstruction for the candidate jet, the $\gamma$+jet invariant mass spectrum measured in data is examined for the presence of resonances over the standard model continuum background. The background is estimated by fitting the mass distribution with a functional form. The data exhibit no statistically significant deviations from the expected standard model background. Exclusion limits at 95% confidence level on the resonance mass and other parameters are set. Excited light-flavor quarks (excited bottom quarks) are excluded up to a mass of 6.0 (3.8) TeV. Quantum black hole production is excluded for masses up to 7.5 (5.2) TeV in the Arkani-Hamed-Dimopoulos-Dvali (Randall-Sundrum) model. These lower mass bounds are the most stringent to date among those obtained in the $\gamma$+jet final state.

The transverse momentum balance of pairs of back-to-back b quark jets in PbPb and pp collisions recorded with the CMS detector at the LHC is reported. The center-of-mass energy in both collision systems is 5.02 TeV per nucleon pair. Compared to the pp collision baseline, b quark jets have a larger imbalance in the most central PbPb collisions, as expected from the jet quenching effect. The data are also compared to the corresponding measurement with inclusive dijets. In the most central collisions, the imbalance of b quark dijets is comparable to that of inclusive dijets.

A search for the production of pairs of heavy Majorana neutrinos (N$_\ell$) from the decays of Z' bosons is performed using the CMS detector at the LHC. The data were collected in proton-proton collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, with an integrated luminosity of 138 fb$^{-1}$. The signature for the search is an excess in the invariant mass distribution of the final-state objects, two same-flavor leptons (e or $\mu$) and at least two jets. No significant excess of events beyond the expected background is observed. Upper limits at 95% confidence level are set on the product of the Z' production cross section and its branching fraction to a pair of N$_\ell$, as functions of N$_\ell$ and Z' boson masses ($m_{\mathrm{N}_\ell}$ and $m_\mathrm{Z'}$, respectively) for $m_\mathrm{Z'}$ from 0.4 to 4.6 TeV and $m_{\mathrm{N}_\ell}$ from 0.1 TeV to $m_\mathrm{Z'}$/2. In the theoretical framework of a left-right symmetric model, exclusion bounds in the $m_{\mathrm{N}_\ell}-m_\mathrm{Z'}$ plane are presented in both the electron and muon channels. The observed upper limit on $m_\mathrm{Z'}$ reaches up to 4.42 TeV. These are the most restrictive limits to date on the mass of N$_\ell$ as a function of the Z' boson mass.

The production cross sections of B$^0_\mathrm{s}$ mesons and charge conjugates are measured in proton-proton (pp) and PbPb collisions via the exclusive decay channel B$^0_\mathrm{s}$ $\to$ $\mu^+\mu^-$K$^+$K$^-$ at a centre-of-mass energy of 5.02 TeV per nucleon pair and within the rapidity range $|y|$ $\lt$ 2.4 using the CMS detector at the LHC. The pp measurement is performed as a function of transverse momentum (p$_\mathrm{T}$) of the B$^0_\mathrm{s}$ meson in the range of 7 to 50 GeV/$c$ and is compared to the predictions of perturbative QCD calculations. The B$^0_\mathrm{s}$ production yield in PbPb collisions is measured in two p$_\mathrm{T}$ intervals, 7 to 15 and 15 to 50 GeV/$c$, and compared to the yield in pp collisions in the same kinematic region. The nuclear modification factor ($R_\mathrm{AA}$) is found to be 1.5 $\pm$ 0.6 (stat) $\pm$ 0.5 (syst) for 7-15 GeV/$c$, and 0.87 $\pm$ 0.30 (stat) $\pm$ 0.17 (syst) for 15-50 GeV/$c$, respectively. Within current uncertainties, the results are consistent with models of strangeness enhancement and a suppression as observed for the B$^+$ mesons.

The J/$\psi$$\to$$\mu^+\mu^-\mu^+\mu^-$ decay has been observed with a statistical significance in excess of five standard deviations. The analysis is based on an event sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment in 2018 and corresponding to an integrated luminosity of 33.6 fb${-1}$. Normalizing to the J/$\psi$$\to$$\mu^+\mu^-$ decay mode leads to a branching fraction [10.1$^{+3.3}_{-2.7}$ (stat) $\pm$ 0.4 (syst) ]$\times$ 10$^{-7}$, a value that is consistent with the standard model prediction.

A search for the standard model Higgs boson decaying to a charm quark-antiquark pair, H $\to$$\mathrm{c\bar{c}}$, produced in association with a leptonically decaying V (W or Z) boson is presented. The search is performed with proton-proton collisions at $\sqrt{s}$ = 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb$^{-1}$. Novel charm jet identification and analysis methods using machine learning techniques are employed. The analysis is validated by searching for Z $\to$$\mathrm{c\bar{c}}$ in VZ events, leading to its first observation at a hadron collider with a significance of 5.7 standard deviations. The observed (expected) upper limit on $\sigma$(VH)$ \mathcal{B}$(H $\to$$\mathrm{c\bar{c}}$) is 0.94 (0.50 $^{+0.22}_{-0.15}$) pb at 95% confidence level (CL), corresponding to 14 (7.6 $^{+3.4}_{-2.3}$) times the standard model prediction. For the Higgs-charm Yukawa coupling modifier, $\kappa_\mathrm{c}$, the observed (expected) 95% CL interval is 1.1 $\lt$ $\vert\kappa_\mathrm{C}\vert$ $\lt$ 5.5 ($\vert\kappa_\mathrm{c}\vert$ $\lt$ 3.4), the most stringent constraint to date.

A combination of the results of several searches for the electroweak production of the supersymmetric partners of standard model bosons, and of charged leptons, is presented. All searches use proton-proton collision data at $\sqrt{s}$ = 13 TeV recorded with the CMS detector at the LHC in 2016-2018. The analyzed data correspond to an integrated luminosity of up to 137 fb$^{-1}$. The results are interpreted in terms of simplified models of supersymmetry. Two new interpretations are added with this combination: a model spectrum with the bino as the lightest supersymmetric particle together with mass-degenerate higgsinos decaying to the bino and a standard model boson, and the compressed-spectrum region of a previously studied model of slepton pair production. Improved analysis techniques are employed to optimize sensitivity for the compressed spectra in the wino and slepton pair production models. The results are consistent with expectations from the standard model. The combination provides a more comprehensive coverage of the model parameter space than the individual searches, extending the exclusion by up to 125 GeV, and also targets some of the intermediate gaps in the mass coverage.

A search for long-lived particles (LLPs) decaying in the CMS muon detectors is presented. A data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded at the LHC in 2016-2018, is used. The decays of LLPs are reconstructed as high multiplicity clusters of hits in the muon detectors. In the context of twin Higgs models, the search is sensitive to LLP masses from 0.4 to 55 GeV and a broad range of LLP decay modes, including decays to hadrons, $\tau$ leptons, electrons, or photons. No excess of events above the standard model background is observed. The most stringent limits to date from LHC data are set on the branching fraction of the Higgs boson decay to a pair of LLPs with masses below 10 GeV. This search also provides the best limits for various intervals of LLP proper decay length and mass. Finally, this search sets the first limits at the LHC on a dark quantum chromodynamic sector whose particles couple to the Higgs boson through gluon, Higgs boson, photon, vector, and dark-photon portals, and is sensitive to branching fractions of the Higgs boson to dark quarks as low as 2$\times$10$^{-3}$.

The first measurement of the dependence of $\gamma\gamma$$\to$$\mu^{+}\mu^{-}$ production on the multiplicity of neutrons emitted very close to the beam direction in ultraperipheral heavy ion collisions is reported. Data for lead-lead interactions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV, with an integrated luminosity of approximately 1.5 nb$^{-1}$, were collected using the CMS detector at the LHC. The azimuthal correlations between the two muons in the invariant mass region 8 $\lt$$m_{\mu\mu}$$\lt$ 60 GeV are extracted for events including 0, 1, or at least 2 neutrons detected in the forward pseudorapidity range $|\eta|$$\gt$ 8.3. The back-to-back correlation structure from leading-order photon-photon scattering is found to be significantly broader for events with a larger number of emitted neutrons from each nucleus, corresponding to interactions with a smaller impact parameter. This observation provides a data-driven demonstration that the average transverse momentum of photons emitted from relativistic heavy ions has an impact parameter dependence. These results provide new constraints on models of photon-induced interactions in ultraperipheral collisions. They also provide a baseline to search for possible final-state effects on lepton pairs caused by traversing a quark-gluon plasma produced in hadronic heavy ion collisions.