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.

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.

Many measurements at the LHC require efficient identification of heavy-flavour jets, i.e. jets originating from bottom (b) or charm (c) quarks. An overview of the algorithms used to identify c jets is described and a novel method to calibrate them is presented. This new method adjusts the entire distributions of the outputs obtained when the algorithms are applied to jets of different flavours. It is based on an iterative approach exploiting three distinct control regions that are enriched with either b jets, c jets, or light-flavour and gluon jets. Results are presented in the form of correction factors evaluated using proton-proton collision data with an integrated luminosity of 41.5 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected by the CMS experiment in 2017. The closure of the method is tested by applying the measured correction factors on simulated data sets and checking the agreement between the adjusted simulation and collision data. Furthermore, a validation is performed by testing the method on pseudodata, which emulate different miscalibration conditions. The calibrated results enable the use of the full distributions of heavy-flavour identification algorithm outputs, e.g. as inputs to machine-learning models. Thus, they are expected to increase the sensitivity of future physics analyses.

A search for decays to invisible particles of Higgs bosons produced in association with a top-antitop quark pair or a vector boson, which both decay to a fully hadronic final state, has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. The 95% confidence level upper limit set on the branching fraction of the 125 GeV Higgs boson to invisible particles, $\mathcal{B}$(H $\to$ inv), is 0.54 (0.39 expected), assuming standard model production cross sections. The results of this analysis are combined with previous $\mathcal{B}$(H $\to$ inv) searches carried out at $\sqrt{s}$ = 7, 8, and 13 TeV in complementary production modes. The combined upper limit at 95% confidence level on $\mathcal{B}$(H $\to$ inv) is 0.15 (0.08 expected).

The observation of the production of four top quarks in proton-proton collisions is reported, based on a data sample collected by the CMS experiment at a center-of-mass energy of 13 TeV in 2016-2018 at the CERN LHC and corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with two same-sign, three, or four charged leptons (electrons and muons) and additional jets are analyzed. Compared to previous results in these channels, updated identification techniques for charged leptons and jets originating from the hadronization of b quarks, as well as a revised multivariate analysis strategy to distinguish the signal process from the main backgrounds, lead to an improved expected signal significance of 4.9 standard deviations above the background-only hypothesis. Four top quark production is observed with a significance of 5.6 standard deviations, and its cross section is measured to be 17.7 $^{+3.7}_{-3.5}$ (stat) $^{+2.3}_{-1.9}$ (syst) fb, in agreement with the available standard model predictions.

A search is performed for charged-lepton flavor violating processes in top quark (t) production and decay. The data were collected by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV and correspond to an integrated luminosity of 138 fb$^{-1}$. The selected events are required to contain one opposite-sign electron-muon pair, a third charged lepton (electron or muon), and at least one jet of which no more than one is associated with a bottom quark. Boosted decision trees are used to distinguish signal from background, exploiting differences in the kinematics of the final states particles. The data are consistent with the standard model expectation. Upper limits at 95% confidence level are placed in the context of effective field theory on the Wilson coefficients, which range between 0.024-0.424 TeV$^{-2}$ depending on the flavor of the associated light quark and the Lorentz structure of the interaction. These limits are converted to upper limits on branching fractions involving up (charm) quarks, t$\to$e$\mu$u (t$\to$e$\mu$c), of 0.032 (0.498)$\times$10$^{-6}$, 0.022 (0.369)$\times$10$^{-6}$, and 0.012 (0.216)$\times$10$^{-6}$ for tensor-like, vector-like, and scalar-like interactions, respectively.

A search is presented for fractionally charged particles with charge below 1$e$, using their small energy loss in the tracking detector as a key variable to observe a signal. The analyzed data set corresponds to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions collected at $\sqrt{s}$ = 13 TeV in 2016-2018 at the CERN LHC. This is the first search at the LHC for new particles with charges between $e/$3 and $e$. Masses up to 640 GeV and charges as low as $e/$3 are excluded at 95% confidence level. These are the most stringent limits to date for the considered Drell-Yan-like production mode.

The inclusive jet cross section is measured as a function of jet transverse momentum $p_\mathrm{T}$ and rapidity $y$. The measurement is performed using proton-proton collision data at $\sqrt{s}$ = 5.02 TeV, recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 27.4 pb$^{-1}$. The jets are reconstructed with the anti-$k_\mathrm{T}$ algorithm using a distance parameter of $R$ = 0.4, within the rapidity interval $\lvert y\rvert$$\lt$ 2, and across the kinematic range 0.06 $\lt$$p_\mathrm{T}$$\lt$ 1 TeV. The jet cross section is unfolded from detector to particle level using the determined jet response and resolution. The results are compared to predictions of perturbative quantum chromodynamics, calculated at both next-to-leading order and next-to-next-to-leading order. The predictions are corrected for nonperturbative effects, and presented for a variety of parton distribution functions and choices of the renormalization/factorization scales and the strong coupling $\alpha_\mathrm{S}$.

A search for an exotic decay of the Higgs boson to a Z boson and a light pseudoscalar particle (a), decaying to a pair of leptons and a pair of photons, respectively, is presented. The search is based on proton-proton collision data at a center-of-mass energy of $\sqrt{s}$ = 13 TeV, collected with the CMS detector and corresponding to an integrated luminosity of 138 fb$^{-1}$. The analysis probes pseudoscalar masses $m_\mathrm{a}$ between 1 and 30 GeV, leading to two pairs of well-isolated leptons and photons. Upper limits at 95% confidence level are set on the Higgs boson production cross section times its branching fraction to two leptons and two photons. The observed (expected) limits are in the range of 1.1-17.8 (1.7-17.9) fb within the probed $m_\mathrm{a}$ interval. An excess of data above the expected standard model background with a local (global) significance of 2.6 (1.3) standard deviations is observed for a mass hypothesis of $m_\mathrm{a}$ = 3 GeV. Limits on models involving axion-like particles, formulated as an effective field theory, are also reported.

A search for W' bosons decaying to a top and a bottom quark in final states including an electron or a muon is performed with the CMS detector at the LHC. The analyzed data correspond to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 Tev. Good agreement with the standard model expectation is observed and no evidence for the existence of the W' boson is found over the mass range examined. The largest observed deviation from the standard model expectation is found for a W' boson mass ($m_\mathrm{W'}$) hypothesis of 3.8 TeV with a relative decay width of 1%, with a local (global) significance of 2.6 (2.0) standard deviations. Upper limits on the production cross sections of W' bosons decaying to a top and a bottom quark are set. Left- and right-handed W' bosons with $m_\mathrm{W'}$ below 3.9 and 4.3 TeV, respectively, are excluded at the 95% confidence level, under the assumption that the new particle has a narrow decay width. Limits are also set for relative decay widths up to 30%. These are the most stringent limits to date on this W' boson decay channel.

Measurements of inclusive and normalized differential cross sections of the associated production of top quark-antiquark and bottom quark-antiquark pairs, ttbb, are presented. The results are based on data from proton-proton collisions collected by the CMS detector at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. The cross sections are measured in the lepton+jets decay channel of the top quark pair, using events containing exactly one isolated electron or muon and at least five jets. Measurements are made in four fiducial phase space regions, targeting different aspects of the ttbb process. Distributions are unfolded to the particle level through maximum likelihood fits, and compared with predictions from several event generators. The inclusive cross section measurements of this process in the fiducial phase space regions are the most precise to date. In most cases, the measured inclusive cross sections exceed the predictions with the chosen generator settings. The only exception is when using a particular choice of dynamic renormalization scale, $\mu_\mathrm{R}=\frac{1}{2} \prod_{i=\mathrm{t, \bar{t}, b, \bar{b}}} m_{\mathrm{T},i}^{1/4}$, where $m_{\mathrm{T},i}^2=m_i^2+p^2_{\mathrm{T},i}$ are the transverse masses of top and bottom quarks. The differential cross sections show varying degrees of compatibility with the theoretical predictions, and none of the tested generators with the chosen settings simultaneously describe all the measured distributions.

The first evidence for the standard model production of a top quark in association with a W boson and a Z boson is reported. The measurement is performed in multilepton final states, where the Z boson is reconstructed via its decays to electron or muon pairs and the W boson decays either to leptons or hadrons. The analysed data were recorded by the CMS experiment at the CERN LHC in 2016-2018 in proton-proton collisions at $\sqrt{s}$ = 13 TeV, and correspond to an integrated luminosity of 138 fb$^{-1}$. The measured cross section is 354 $\pm$ 54 (stat) $\pm$ 95 (syst) fb, and corresponds to a statistical significance of 3.4 standard deviations.

A search is presented for a third-generation leptoquark (LQ) coupled exclusively to a $\tau$ lepton and a b quark. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV recorded with the CMS detector, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events with $\tau$ leptons and a varying number of jets originating from b quarks are considered, targeting the single and pair production of LQs, as well as nonresonant $t$-channel LQ exchange. An excess is observed in the data with respect to the background expectation in the combined analysis of all search regions. For a benchmark LQ mass of 2 TeV and an LQ-b-$\tau$ coupling strength of 2.5, the excess reaches a local significance of up to 2.8 standard deviations. Upper limits at the 95% confidence level are placed on the LQ production cross section in the LQ mass range 0.5-2.3 TeV, and up to 3 TeV for $t$-channel LQ exchange. Leptoquarks are excluded below masses of 1.22-1.88 TeV for different LQ models and varying coupling strengths up to 2.5. The study of nonresonant $\tau\tau$ production through $t$-channel LQ exchange allows lower limits on the LQ mass of up to 2.3 TeV to be obtained.

The first search for narrow resonances decaying to three well-separated hadronic jets is presented. The search uses proton-proton collision data corresponding to an integrated luminosity of 138 fb$^{-1}$ at $\sqrt{s}$ = 13 TeV, collected at the CERN LHC. No significant deviations from the background predictions are observed between 1.75-9.00 TeV. The results provide the first mass limits on a right-handed boson Z$_{\mathrm{R}}$ decaying to three gluons, an excited quark decaying via a vector boson to three quarks, as well as updated limits on a Kaluza-Klein gluon decaying via a radion to three gluons.

A measurement is presented of the primary Lund jet plane (LJP) density in inclusive jet production in proton-proton collisions. The analysis uses 138 fb$^{-1}$ of data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV. The LJP, a representation of the phase space of emissions inside jets, is constructed using iterative jet declustering. The transverse momentum $k_\mathrm{T}$ and the splitting angle $\Delta R$ of an emission relative to its emitter are measured at each step of the jet declustering process. The average density of emissions as function of $\ln(k_\mathrm{T}$/GeV) and $\ln(R/\Delta R)$ is measured for jets with distance parameters $R$ = 0.4 or 0.8, transverse momentum $p_\mathrm{T} \gt$ 700 GeV, and rapidity $\vert y\vert \lt $ 1.7. The jet substructure is measured using the charged-particle tracks of the jet. The measured distributions, unfolded to the level of stable particles, are compared with theoretical predictions from simulations and with perturbative quantum chromodynamics calculations. Due to the ability of the LJP to factorize physical effects, these measurements can be used to improve different aspects of the physics modeling in event generators.

Diboson production in association with jets is studied in the fully leptonic final states, pp $\to$ (Z$\gamma^*$)(Z/$\gamma^*$)+jets $\to$ 2$\ell$2$\ell'$+jets, ($\ell,\ell'$ = e or $\mu$) in proton-proton collisions at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the LHC. Differential distributions and normalized differential cross sections are measured as a function of jet multiplicity, transverse momentum $p_\mathrm{T}$, pseudorapidity $\eta$, invariant mass and $\Delta\eta$ of the highest-$p_\mathrm{T}$ and second-highest-$p_\mathrm{T}$ jets, and as a function of invariant mass of the four-lepton system for events with various jet multiplicities. These differential cross sections are compared with theoretical predictions that mostly agree with the experimental data. However, in a few regions we observe discrepancies between the predicted and measured values. Further improvement of the predictions is required to describe the ZZ+jets production in the whole phase space.

Using proton-proton collision data corresponding to an integrated luminosity of 140 fb$^{-1}$ collected by the CMS experiment at $\sqrt{s}$ = 13 TeV, the $\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$ decay is observed for the first time, with a statistical significance exceeding 5 standard deviations. The relative branching fraction, with respect to the $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$ decay, is measured to be $\mathcal{B}$($\Lambda_\text{b}^0$$\to$ J/$\psi\Xi^-$K$^+$)/$\mathcal{B}$( $\Lambda_\text{b}^0$$\to$$\psi$(2S)$\Lambda$) = [3.38 $\pm$ 1.02 $\pm$ 0.61 $\pm$ 0.03]%, where the first uncertainty is statistical, the second is systematic, and the third is related to the uncertainties in $\mathcal{B}$($\psi$(2S) $\to$ J/$\psi\pi^+\pi^-$) and $\mathcal{B}$($\Xi^-$ $\to$ $\Lambda\pi^-$).

A search for Higgs boson pair (HH) production with one Higgs boson decaying to two bottom quarks and the other to two W bosons are presented. The search is done using proton-proton collisions data at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded by the CMS detector at the LHC from 2016 to 2018. The final states considered include at least one leptonically decaying W boson. No evidence for the presence of a signal is observed and corresponding upper limits on the HH production cross section are derived. The limit on the inclusive cross section of the nonresonant HH production, assuming that the distributions of kinematic observables are as expected in the standard model (SM), is observed (expected) to be 14 (18) times the value predicted by the SM, at 95% confidence level. The limits on the cross section are also presented as functions of various Higgs boson coupling modifiers, and anomalous Higgs boson coupling scenarios. In addition, limits are set on the resonant HH production via spin-0 and spin-2 resonances within the mass range 250-900 GeV.

A search for Higgs boson pair (HH) production in association with a vector boson V (W or Z boson) is presented. The search is based on proton-proton collision data at a center-of-mass energy of 13 TeV, collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. All hadronic and leptonic decays of V bosons are used. The leptons considered are electrons, muons, and neutrinos. The HH production is searched for in the $\mathrm{b\bar{b}b\bar{b}}$ decay channel. An observed (expected) upper limit at 95% confidence level of VHH production cross section is set at 294 (124) times the standard model prediction. Constraints are also set on the modifiers of the Higgs boson trilinear self-coupling, $\kappa_{\lambda}$, assuming $\kappa_{2\mathrm{V}}$ = 1 and vice versa on the coupling of two Higgs bosons with two vector bosons, $\kappa_{2\mathrm{V}}$. The observed (expected) 95% confidence intervals of these coupling modifiers are -37.7 $\lt$ $\kappa_{\lambda}$ $\lt$ 37.2 (-30.1 $\lt$ $\kappa_{\lambda}$ $\lt$ 8.9) and -12.2 $\lt$ $\kappa_{2\mathrm{V}}$ $\lt$ 13.5 (-7.2 $\lt$ $\kappa_{2\mathrm{V}}$ $\lt$ 8.9), respectively.

A search for exotic decays of the Higgs boson (H) with a mass of 125 GeV to a pair of light pseudoscalars $\mathrm{a}_1$ is performed in final states where one pseudoscalar decays to two b quarks and the other to a pair of muons or $\tau$ leptons. A data sample of proton-proton collisions at $\sqrt{s}$ = 13 TeV corresponding to an integrated luminosity of 138 fb$^{-1}$ recorded with the CMS detector is analyzed. No statistically significant excess is observed over the standard model backgrounds. Upper limits are set at 95% confidence level (CL) on the Higgs boson branching fraction to $\mu\mu$bb and to $\tau\tau$bb, via a pair of $\mathrm{a}_1$s. The limits depend on the pseudoscalar mass $m_{\mathrm{a}_1}$ and are observed to be in the range (0.17-3.3) $\times$ 10$^{-4}$ and (1.7-7.7) $\times$ 10$^{2}$ in the $\mu\mu$bb and $\tau\tau$bb final states, respectively. In the framework of models with two Higgs doublets and a complex scalar singlet (2HDM+S), the results of the two final states are combined to determine model-independent upper limits on the branching fraction $\mathcal{B}$(H $\to$ $\mathrm{a}_1\mathrm{a}_1$ $\to$ $\ell\ell$bb) at 95% CL, with $\ell$ being a muon or a $\tau$ lepton. For different types of 2HDM+S, upper bounds on the branching fraction $\mathcal{B}$(H $\to$ $\mathrm{a}_1\mathrm{a}_1$) are extracted from the combination of the two channels. In most of the Type II 2HDM+S parameter space, $\mathcal{B}($H $\to$ $\mathrm{a}_1\mathrm{a}_1$) values above 0.23 are excluded at 95% CL for $m_{\mathrm{a}_1}$ values between 15 and 60 GeV.

This Letter presents a search for direct production of charginos and neutralinos via electroweak interactions. The results are based on data from proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. The search considers final states with large missing transverse momentum and pairs of hadronically decaying bosons WW, WZ, and WH, where H is the Higgs boson. These bosons are identified using novel algorithms. No significant excess of events is observed relative to the expectations from the standard model. Limits at the 95% confidence level are placed on the cross section for production of mass-degenerate wino-like supersymmetric particles $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$, and mass-degenerate higgsino-like supersymmetric particles $\tilde{\chi}_1^\pm$, $\tilde{\chi}_2^0$, and $\tilde{\chi}_3^0$. In the limit of a nearly-massless lightest supersymmetric particle $\tilde{\chi}_1^0$, wino-like particles with masses up to 870 and 960 GeV are excluded in the cases of $\tilde{\chi}_2^0$ $\to$ Z$\tilde{\chi}_1^0$ and $\tilde{\chi}_2^0$ $\to$ H$\tilde{\chi}_1^0$, respectively, and higgsino-like particles are excluded between 300 and 650 GeV.