The first measurement of the CP structure of the Yukawa coupling between the Higgs boson and $\tau$ leptons is presented. The measurement is based on data collected in proton-proton collisions at $\sqrt{s} =$ 13 TeV by the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. The analysis uses the angular correlation between the decay planes of $\tau$ leptons produced in Higgs boson decays. The effective mixing angle between CP-even and CP-odd $\tau$ Yukawa couplings is found to be $-$1 $\pm$ 19$^\circ$, compared to an expected value of 0 $\pm$ 21$^\circ$ at the 68.3% confidence level. The data disfavour the pure CP-odd scenario at 3.0 standard deviations. The results are compatible with predictions for the standard model Higgs boson.

Studies of $CP$ violation and anomalous couplings of the Higgs boson to vector bosons and fermions are presented. The data were acquired by the CMS experiment at the LHC and correspond to an integrated luminosity of 137 fb$^{-1}$ at a proton-proton collision energy of 13 TeV. The kinematic effects in the Higgs boson's four-lepton decay H $\to$ 4$\ell$ and its production in association with two jets, a vector boson, or top quarks are analyzed, using a full detector simulation and matrix element techniques to identify the production mechanisms and to increase sensitivity to the Higgs boson tensor structure of the Higgs boson interactions. A simultaneous measurement is performed of up to five Higgs boson couplings to electroweak vector bosons (HVV), two couplings to gluons (Hgg), and two couplings to top quarks (Htt). The $CP$ measurement in the Htt interaction is combined with the recent measurement in the H $\to$$\gamma\gamma$ channel. The results are presented in the framework of anomalous couplings and are also interpreted in the framework of effective field theory, including the first study of $CP$ properties of the Htt and effective Hgg couplings from a simultaneous analysis of the gluon fusion and top-associated processes. The results are consistent with the standard model of particle physics.

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.

A search is reported for heavy resonances and quantum black holes decaying into e$\mu$, e$\tau$, and $\mu\tau$ final states in proton-proton collision data recorded by the CMS experiment at the CERN LHC during 2016-2018 at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The e$\mu$, e$\tau$, and $\mu\tau$ invariant mass spectra are reconstructed, and no evidence is found for physics beyond the standard model. Upper limits are set at 95% confidence level on the product of the cross section and branching fraction for lepton flavor violating signals. Three benchmark signals are studied: resonant $\tau$ sneutrino production in $R$ parity violating supersymmetric models, heavy Z' gauge bosons with lepton flavor violating decays, and nonresonant quantum black hole production in models with extra spatial dimensions. Resonant $\tau$ sneutrinos are excluded for masses up to 4.2 TeV in the e$\mu$ channel, 3.7 TeV in the e$\tau$ channel, and 3.6 TeV in the $\mu\tau$ channel. A Z' boson with lepton flavor violating couplings is excluded up to a mass of 5.0 TeV in the e$\mu$ channel, up to 4.3 TeV in the e$\tau$ channel, and up to 4.1 TeV in the $\mu\tau$ channel. Quantum black holes in the benchmark model are excluded up to the threshold mass of 5.6 TeV in the e$\mu$ channel, 5.2 TeV in the e$\tau$ channel, and 5.0 TeV in the $\mu\tau$ channel. In addition, model-independent limits are extracted to allow comparisons with other models for the same final states and similar event selection requirements. The results of these searches provide the most stringent limits available from collider experiments for heavy particles that undergo lepton flavor violating decays.

A search is reported for high-mass hadronic resonances that decay to a parton and a Lorentz-boosted resonance, which in turn decays into a pair of partons. The search is based on data collected with the CMS detector at the LHC in proton-proton collisions at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The boosted resonance is reconstructed as a single wide jet with substructure consistent with a two-body decay. The high-mass resonance is thus considered as a dijet system. The jet substructure information and the kinematic properties of cascade resonance decays are exploited to disentangle the signal from the large quantum chromodynamics multijet background. The dijet mass spectrum is analyzed for the presence of new high-mass resonances, and is found to be consistent with the standard model background predictions. Results are interpreted in a warped extra dimension model where the high-mass resonance is a Kaluza-Klein gluon, the boosted resonance is a radion, and the final state partons are all gluons. Limits on the production cross section are set as a function of the Kaluza-Klein gluon and radion masses. These limits exclude at 95% confidence level models with Kaluza-Klein gluon masses in the range from 2.0 to 4.3 TeV and radion masses in the range from 0.20 to 0.74 TeV. By exploring a novel experimental signature, the observed limits on the Kaluza-Klein gluon mass are extended by up to about 1 TeV compared to previous searches.

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 measurement of event-shape variables is presented, using a data sample produced in a special run with approximately one inelastic proton-proton collision per bunch crossing. The data were collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 64 $\mu$b$^{-1}$. A number of observables related to the overall distribution of charged particles in the collisions are corrected for detector effects and compared with simulations. Inclusive event-shape distributions, as well as differential distributions of event shapes as functions of charged-particle multiplicity, are studied. None of the models investigated is able to satisfactorily describe the data. Moreover, there are significant features common amongst all generator setups studied, particularly showing data being more isotropic than any of the simulations. Multidimensional unfolded distributions are provided, along with their correlations.

The results of a search for Higgs boson pair (HH) production in the WW*WW*, WW*$\tau\tau$, and $\tau\tau\tau\tau$ decay modes are presented. The search uses 138 fb$^{-1}$ of proton-proton collision data recorded by the CMS experiment at the LHC at a center-of-mass energy of 13 TeV from 2016 to 2018. Analyzed events contain two, three, or four reconstructed leptons, including electrons, muons, and hadronically decaying tau leptons. No evidence for a signal is found in the data. Upper limits are set on the cross section for nonresonant HH production, as well as resonant production in which a new heavy particle decays to a pair of Higgs bosons. For nonresonant production, the observed (expected) upper limit on the cross section at 95% confidence level (CL) is 21.3 (19.4) times the standard model (SM) prediction. The observed (expected) ratio of the trilinear Higgs boson self-coupling to its value in the SM is constrained to be within the interval $-$6.9 to 11.1 ($-$6.9 to 11.7) at 95% CL, and limits are set on a variety of new-physics models using an effective field theory approach. The observed (expected) limits on the cross section for resonant HH production range from 0.18 to 0.90 (0.08 to 1.06) pb at 95% CL for new heavy-particle masses in the range 250-1000 GeV.

An inclusive search for long-lived exotic particles decaying to a pair of muons is presented. The search uses data collected by the CMS experiment at the CERN LHC in proton-proton collisions at $\sqrt{s}$ = 13 TeV in 2016 and 2018 and corresponding to an integrated luminosity of 97.6 fb$^{-1}$. The experimental signature is a pair of oppositely charged muons originating from a common secondary vertex spatially separated from the pp interaction point by distances ranging from several hundred $\mu$m to several meters. The results are interpreted in the frameworks of the hidden Abelian Higgs model, in which the Higgs boson decays to a pair of long-lived dark photons Z$_\mathrm{D}$, and of a simplified model, in which long-lived particles are produced in decays of an exotic heavy neutral scalar boson. For the hidden Abelian Higgs model with $m_\mathrm{Z_D}$ greater than 20 GeV and less than half the mass of the Higgs boson, they provide the best limits to date on the branching fraction of the Higgs boson to dark photons for $c\tau$(Z$_\mathrm{D}$) (varying with $m_\mathrm{Z_D}$) between 0.03 and ${\approx}$ 0.5 mm, and above ${\approx}$ 0.5 m. Our results also yield the best constraints on long-lived particles with masses larger than 10 GeV produced in decays of an exotic scalar boson heavier than the Higgs boson and decaying to a pair of muons.

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}$.