Multiparticle azimuthal correlations of prompt D$^0$ mesons are measured in PbPb collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. For the first time, a four-particle cumulant method is used to extract the second Fourier coefficient of the azimuthal distribution ($v_2$) of D$^0$ mesons as a function of event centrality and the D$^0$ transverse momentum. The ratios of the four-particle $v_2$ values to previously measured two-particle cumulant results provide direct experimental access to event-by-event fluctuations of charm quark azimuthal anisotropies. These ratios are also found to be comparable to those of inclusive charged particles in the event. However, hints of deviations are seen in the most central and peripheral collisions. To investigate the origin of flow fluctuations in the charm sector, these measurements are compared with models implementing fluctuations of charm quark energy loss via collisional or radiative processes in the quark-gluon plasma. These models cannot quantitatively describe the data over the full transverse momentum and centrality ranges, although the calculations with collisional energy loss provide a better description of the data.

The transverse momentum spectra of charged particles have been measured in pp and PbPb collisions at sqrt(sNN) = 2.76 TeV by the CMS experiment at the LHC. In the transverse momentum range pt = 5-10 GeV/c, the charged particle yield in the most central PbPb collisions is suppressed by up to a factor of 5 compared to the pp yield scaled by the number of incoherent nucleon-nucleon collisions. At higher pt, this suppression is significantly reduced, approaching roughly a factor of 2 for particles with pt in the range pt=40-100 GeV/c.

The $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

A search for pairs of Higgs bosons produced via gluon and vector boson fusion is presented, focusing on the four b quark final state. The data sample consists of proton-proton collisions at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC, and corresponds to an integrated luminosity of 138 fb$^{-1}$. No deviation from the background-only hypothesis is observed. A 95% confidence level upper limit on the Higgs boson pair production cross section is observed at 3.9 times the standard model prediction for an expected value of 7.8. Constraints are also set on the modifiers of the Higgs field self-coupling, $\kappa_\lambda$, and of the coupling of two Higgs bosons to two vector bosons, $\kappa_\mathrm{2V}$. The observed (expected) allowed intervals at the 95% confidence level are $-$2.3 $\lt$ $\kappa_\lambda$ $\lt$ 9.4 ($-$5.0 $\lt$ $\kappa_\lambda$ $\lt$ 12.0) and $-$0.1 $\lt$ $\kappa_\mathrm{2V}$ $\lt$ 2.2 ($-$0.4 $\lt$ $\kappa_\mathrm{2V}$ $\lt$ 2.5). These are the most stringent observed constraints to date on the HH production cross section and on the $\kappa_\mathrm{2V}$ coupling.

A search for flavor-changing neutral current interactions of the top quark (t) and the Higgs boson (H) is presented. The search is based on a data sample corresponding to an integrated luminosity of 137 fb$^{-1}$ recorded by the CMS experiment at the LHC in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Events containing exactly one lepton (muon or electron) and at least three jets, among which at least two are identified as originating from the hadronization of a bottom quark, are analyzed. A set of deep neural networks is used for kinematic event reconstruction, while boosted decision trees distinguish the signal from the background events. No significant excess over the background predictions is observed, and upper limits on the signal production cross sections are extracted. These limits are interpreted in terms of top quark decay branching fractions ($\mathcal{B}$) to the Higgs boson and an up (u) or a charm quark (c). Assuming one nonvanishing extra coupling at a time, the observed (expected) upper limits at 95% confidence level are $\mathcal{B}$(t $\to$ Hu) $\lt$ 0.079 (0.11)% and $\mathcal{B}$(t $\to$ Hc) $\lt$ 0.094 (0.086)%.

A search has been performed for heavy resonances decaying to ZZ or ZW and for axion-like particles (ALPs) mediating nonresonant ZZ or ZH production, in final states with two charged leptons ($\ell$ = e, $\mu$) produced by the decay of a Z boson, and two quarks produced by the decay of a Z, W, or Higgs boson H. The analysis is sensitive to resonances with masses in the range 450 to 2000 GeV. Two categories are defined corresponding to the merged or resolved reconstruction of the hadronically decaying boson. The search is based on data collected during 2016-2018 by the CMS experiment at the LHC in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. No significant excess is observed in the data above the standard model background expectation. Upper limits on the production cross section of heavy, narrow spin-2 and spin-1 resonances are derived as functions of the resonance mass, and exclusion limits on the production of bulk graviton particles and W$'$ bosons are calculated in the framework of the warped extra dimensions and heavy vector triplet models, respectively. In addition, upper limits on the ALP-mediated diboson production cross section and ALP couplings to standard model particles are obtained in the framework of linear and chiral effective field theories. These are the first limits on nonresonant ALP-mediated ZZ and ZH production obtained by the LHC experiments.

A search is presented for a right-handed W boson (W$_\mathrm{R}$) and a heavy neutrino (N), in a final state consisting of two same-flavor leptons (ee or $\mu\mu$) and two quarks. The search is performed with the CMS experiment at the CERN LHC using a data sample of proton-proton collisions at a center-of-mass energy of 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$. The search covers two regions of phase space, one where the decay products of the heavy neutrino are merged into a single large-area jet, and one where the decay products are well separated. The expected signal is characterized by an excess in the invariant mass distribution of the final-state objects. No significant excess over the standard model background expectations is observed. The observations are interpreted as upper limits on the product of W$_\mathrm{R}$ production cross sections and branching fractions assuming that couplings are identical to those of the standard model W boson. For N masses $m_\mathrm{N}$ equal to half the W$_\mathrm{R}$ mass $m_\mathrm{W_R}$ ($m_\mathrm{N}$ = 0.2 TeV), $m_\mathrm{W_R}$ is excluded at 95% confidence level up to 4.7 (4.8) and 5.0 (5.4) TeV for the electron and muon channels, respectively. This analysis provides the most stringent limits on the W$_\mathrm{R}$ mass to date.

A search for heavy neutral leptons (HNLs), the right-handed Dirac or Majorana neutrinos, is performed in final states with three charged leptons (electrons or muons) using proton-proton collision data collected by the CMS experiment at $\sqrt{s} =$ 13 TeV at the CERN LHC. The data correspond to an integrated luminosity of 138 fb$^{-1}$. The HNLs could be produced through mixing with standard model neutrinos $\nu$. For small values of the HNL mass ($\lt$ 20 GeV) and the square of the HNL-$\nu$ mixing parameter (10$^{-7}$-10$^{-2}$), the decay length of these particles can be large enough so that the secondary vertex of the HNL decay can be resolved with the CMS silicon tracker. The selected final state consists of one lepton emerging from the primary proton-proton collision vertex, and two leptons forming a displaced, secondary vertex. No significant deviations from the standard model expectations are observed, and constraints are obtained on the HNL mass and coupling strength parameters, excluding previously unexplored regions of parameter space in the mass range 1-20 GeV and squared mixing parameter values as low as 10$^{-7}$.

The measurement of the cross section for the production of a Z boson, decaying to dielectrons or dimuons, in association with at least one bottom quark jet are performed with proton-proton collision data at $\sqrt{s} =$ 13 TeV. The data sample corresponds to an integrated luminosity of 137 fb$^{-1}$, collected by the CMS experiment at the LHC during 2016-2018. The integrated cross sections for Z + $\ge$ 1 b jet and Z + $\ge$ 2 b jets are reported for the electron, muon, and combined channels. The fiducial cross sections in the combined channel are 6.52 $\pm$ 0.04 (stat) $\pm$ 0.40 (syst) $\pm$ 0.14 (theo) pb for Z + $\ge$ 1 b jet and 0.65 $\pm$ 0.03 (stat) $\pm$ 0.07 (syst) $\pm$ 0.02 (theo) pb for Z + $\ge$ 2 b jets. The differential cross section distributions are measured as functions of various kinematic observables that are useful for precision tests of perturbative quantum chromodynamics predictions. The ratios of integrated and differential cross sections for Z + $\ge$ 2 b jets and Z + $\ge$ 1 b jet processes are also determined. The value of the integrated cross section ratio measured in the combined channel is 0.100 $\pm$ 0.005 (stat) $\pm$ 0.007 (syst) $\pm$ 0.003 (theo). All measurements are compared with predictions from various event generators.

Results are presented from a search for physics beyond the standard model in proton-proton collisions at $\sqrt{s} =$ 13 TeV in channels with two Higgs bosons, each decaying via the process H $\to$$\mathrm{b\bar{b}}$, and large missing transverse momentum. The search uses a data sample corresponding to an integrated luminosity of 137 fb$^{-1}$ collected by the CMS experiment at the CERN LHC. The search is motivated by models of supersymmetry that predict the production of neutralinos, the neutral partners of the electroweak gauge and Higgs bosons. The observed event yields in the signal regions are found to be consistent with the standard model background expectations. The results are interpreted using simplified models of supersymmetry. For the electroweak production of nearly mass-degenerate higgsinos, each of whose decay chains yields a neutralino ($\tilde{\chi}^0_1$) that in turn decays to a massless goldstino and a Higgs boson, $\tilde{\chi}^0_1$ masses in the range 175 to 1025 GeV are excluded at 95% confidence level. For the strong production of gluino pairs decaying via a slightly lighter $\tilde{\chi}^0_2$ to H and a light $\tilde{\chi}^0_1$, gluino masses below 2330 GeV are excluded.