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.

The cross section for W or Z boson production in association with two photons is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV. The data set corresponds to an integrated luminosity of 137 fb$^{-1}$ collected by the CMS experiment at the LHC. The W $\to$$\ell\nu$ and Z $\to$$\ell\ell$ decay modes (where $\ell =$ e, $\mu$) are used to extract the W$\gamma\gamma$ and Z$\gamma\gamma$ cross sections in a phase space defined by electron (muon) with transverse momentum larger than 30 GeV and photon transverse momentum larger than 20 GeV. All leptons and photons are required to have absolute pseudorapidity smaller than 2.5. The measured cross sections in this phase space are $\sigma$(W$\gamma\gamma$) = 13.6 $^{+1.9}_{-1.9}$ (stat) ${}^{+4.0}_{-4.0}$ (syst) $\pm$ 0.08 (PDF + scale) fb and $\sigma$(Z$\gamma\gamma$) = 5.41 $^{+0.58}_{-0.55}$ (stat) ${}^{+0.64}_{-0.70}$ (syst) $\pm$ 0.06 (PDF + scale) fb. Limits on anomalous quartic gauge couplings are set in the framework of an effective field theory with dimension-8 operators.

Double-parton scattering is investigated using events with a Z boson and jets. The Z boson is reconstructed using only the dimuon channel. The measurements are performed with proton-proton collision data recorded by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$ collected in the year 2016. Differential cross sections of Z + $\geq$ 1 jet and Z + $\geq$ 2 jets are measured with transverse momentum of the jets above 20 GeV and pseudorapidity $|\eta|$$\lt$ 2.4. Several distributions with sensitivity to double-parton scattering effects are measured as functions of the angle and the transverse momentum imbalance between the Z boson and the jets. The measured distributions are compared with predictions from several event generators with different hadronization models and different parameter settings for multiparton interactions. The measured distributions show a dependence on the hadronization and multiparton interaction simulation parameters, and are important input for future improvements of the simulations.

The production of a top quark-antiquark pair in association with a W boson ($\mathrm{t\bar{t}}$W) is measured in proton-proton collisions at a center-of-mass energy of 13 TeV. The analyzed data was recorded by the CMS experiment at the CERN LHC and corresponds to an integrated luminosity of 138 fb$^{-1}$. Events with two or three leptons (electrons and muons) and additional jets are selected. In events with two leptons, a multiclass neural network is used to distinguish between the signal and background processes. Events with three leptons are categorized based on the number of jets and of jets originating from b quark hadronization, and the lepton charges. The inclusive $\mathrm{t\bar{t}}$W production cross section in the full phase space is measured to be 868 $\pm$ 40 (stat) $\pm$ 51 (syst) fb. The $\mathrm{t\bar{t}}$W$^+$ and $\mathrm{t\bar{t}}$W$^-$ cross sections are also measured as 553 $\pm$ 30 (stat) $\pm$ 30 (syst) and 343 $\pm$ 26 (stat) $\pm$ 25 (syst) fb, respectively, and the corresponding ratio of the two cross sections is found to be 1.61 $\pm$ 0.15 (stat) $^{+0.07}_{-0.05}$ (syst). The measured cross sections are larger than but consistent with the standard model predictions within two standard deviations, and represent the most precise measurement of these cross sections to date.

Measurements of Higgs boson production cross sections and couplings in events where the Higgs boson decays into a pair of photons are reported. Events are selected from a sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV collected by the CMS detector at the LHC from 2016 to 2018, corresponding to an integrated luminosity of 137 fb$^{-1}$. Analysis categories enriched in Higgs boson events produced via gluon fusion, vector boson fusion, vector boson associated production, and production associated with top quarks are constructed. The total Higgs boson signal strength, relative to the standard model (SM) prediction, is measured to be 1.12 $\pm$ 0.09. Other properties of the Higgs boson are measured, including SM signal strength modifiers, production cross sections, and its couplings to other particles. These include the most precise measurements of gluon fusion and vector boson fusion Higgs boson production in several different kinematic regions, the first measurement of Higgs boson production in association with a top quark pair in five regions of the Higgs boson transverse momentum, and an upper limit on the rate of Higgs boson production in association with a single top quark. All results are found to be in agreement with the SM expectations.

The hydrodynamic flow-like behavior of charged hadrons in high-energy lead-lead collisions is studied through multiparticle correlations. The elliptic anisotropy values based on different orders of multiparticle cumulants, $v_{2}\{2k\}$, are measured up to the tenth order ($k$ = 5) as functions of the collision centrality at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The data were recorded by the CMS experiment at the LHC and correspond to an integrated luminosity of 0.607 nb$^{-1}$. A hierarchy is observed between the coefficients, with $v_{2}\{2\} > v_{2}\{4\} \gtrsim v_{2}\{6\} \gtrsim v_{2}\{8\} \gtrsim v_{2}\{10\}$. Based on these results, centrality-dependent moments for the fluctuation-driven event-by-event $v_{2}$ distribution are determined, including the skewness, kurtosis and, for the first time, superskewness. Assuming a hydrodynamic expansion of the produced medium, these moments directly probe the initial-state geometry in high-energy nucleus-nucleus collisions.

A search for supersymmetry is presented in events with a single charged lepton, electron or muon, and multiple hadronic jets. The data correspond to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 TeV, recorded by the CMS experiment at the CERN LHC. The search targets gluino pair production, where the gluinos decay into final states with the lightest supersymmetric particle (LSP) and either a top quark-antiquark ($\mathrm{t\bar{t}}$) pair, or a light-flavor quark-antiquark ($\mathrm{q\bar{q}}$) pair and a virtual or on-shell W boson. The main backgrounds, $\mathrm{t\bar{t}}$ pair and W+jets production, are suppressed by requirements on the azimuthal angle between the momenta of the lepton and of its reconstructed parent W boson candidate, and by top quark and W boson identification based on a machine-learning technique. The number of observed events is consistent with the expectations from standard model processes. Limits are evaluated on supersymmetric particle masses in the context of two simplified models of gluino pair production. Exclusions for gluino masses reach up to 2120 (2050) GeV at 95% confidence level for a model with gluino decay to a $\mathrm{t\bar{t}}$ pair (a $\mathrm{q\bar{q}}$ pair and a W boson) and the LSP. For the same models, limits on the mass of the LSP reach up to 1250 (1070) GeV.

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 strange quark content of the proton is probed through the measurement of the production cross section for a W boson and a charm (c) quark in proton-proton collisions at a center-of-mass energy of 13 TeV. The analysis uses a data sample corresponding to a total integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. The W bosons are identified through their leptonic decays to an electron or a muon, and a neutrino. Charm jets are tagged using the presence of a muon or a secondary vertex inside the jet. The W+c production cross section and the cross section ratio $R^\pm_\text{c}$ = $\sigma$(W$^+$+$\bar{\text{c}}$) / $\sigma$(W$^-$+$\text{c}$) are measured inclusively and differentially as functions of the transverse momentum and the pseudorapidity of the lepton originating from the W boson decay. The precision of the measurements is improved with respect to previous studies, reaching 1% in $R^\pm_\text{c}$. The precision of the measurements is improved with respect to previous studies, reaching 1% in $R^\pm_\text{c}$ = 0.950 $\pm$ 0.005 (stat) $\pm$ 0.010 (syst). The measurements are compared with theoretical predictions up to next-to-next-to-leading order in perturbative quantum chromodynamics.

Production cross sections of the standard model Higgs boson decaying to a pair of W bosons are measured in proton-proton collisions at a center-of-mass energy of 13 TeV. The analysis targets Higgs bosons produced via gluon fusion, vector boson fusion, and in association with a W or Z boson. Candidate events are required to have at least two charged leptons and moderate missing transverse momentum, targeting events with at least one leptonically decaying W boson originating from the Higgs boson. Results are presented in the form of inclusive and differential cross sections in the simplified template cross section framework, as well as couplings of the Higgs boson to vector bosons and fermions. The data set collected by the CMS detector during 2016-2018 is used, corresponding to an integrated luminosity of 138 fb$^{-1}$. The signal strength modifier $\mu$, defined as the ratio of the observed production rate in a given decay channel to the standard model expectation, is measured to be $\mu$ = 0.95 $^{+0.10}_{-0.09}$. All results are found to be compatible with the standard model within the uncertainties.