Experimental results are presented on event-by-event net-proton fluctuation measurements in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV, recorded by the ALICE detector at the CERN LHC. These measurements have as their ultimate goal an experimental test of Lattice QCD (LQCD) predictions on second and higher order cumulants of net-baryon distributions to search for critical behavior near the QCD phase boundary. Before confronting them with LQCD predictions, account has to be taken of correlations stemming from baryon number conservation as well as fluctuations of participating nucleons. Both effects influence the experimental measurements and are usually not considered in theoretical calculations. For the first time, it is shown that event-by-event baryon number conservation leads to subtle long-range correlations arising from very early interactions in the collisions.

A measurement is presented of the ratio of the inclusive 3-jet cross section to the inclusive 2-jet cross section as a function of the average transverse momentum, <pT[1,2]>, of the two leading jets in the event. The data sample was collected during 2011 at a proton-proton centre-of-mass energy of 7 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 5.0 inverse femtobarns. The strong coupling constant at the scale of the Z boson mass is determined to be alphaS[MZ] = 0.1148 +/- 0.0014 (exp.) +/- 0.0018 (PDF) +/- 0.0050 (theory), by comparing the ratio in the range 0.42 < <pT[1,2]> < 1.39 TeV to the predictions of perturbative QCD at next-to-leading order. This is the first determination of alphaS[MZ] from measurements at momentum scales beyond 0.6 TeV. The predicted ratio depends only indirectly on the evolution of the parton distribution functions of the proton such that this measurement also serves as a test of the evolution of the strong coupling constant. No deviation from the expected behaviour is observed.

A search for a light charged Higgs boson (H$^+$) decaying to a W boson and a CP-odd Higgs boson (A) in final states with e$\mu\mu$ or $\mu\mu\mu$ is performed using data from pp collisions at $\sqrt{s}=$ 13 TeV, recorded by the CMS detector at the LHC and corresponding to an integrated luminosity of 35.9 fb$^{-1}$. In this search, it is assumed that the H$^+$ boson is produced in decays of top quarks, and the A boson decays to two oppositely charged muons. The presence of signals for H$^+$ boson masses between 100 and 160 GeV and A boson masses between 15 and 75 GeV is investigated. No evidence for the production of the H$^+$ boson is found. Upper limits at 95% confidence level are obtained on the combined branching fraction for the decay chain t $\to$ bH$^+$ $\to$ bW$^+$A $\to$ bW$^+\mu^+\mu^-$, of 1.9 $\times$ 10$^{-6}$ to 8.6 $\times$ 10$^{-6}$, depending on the masses of the H$^+$ and A bosons. These are the first limits for these decay modes of the H$^+$ and A bosons.

The dependence of the ratio between the B$_\mathrm{s}^0$ and B$^+$ hadron production fractions, $f_\mathrm{s} / f_\mathrm{u}$, on the transverse momentum ($p_\mathrm{T}$) and rapidity of the B mesons is studied using the decay channels B$_\mathrm{s}^0$$\to$ J$/\psi\,\phi$ and B$^+$$\to$ J$/\psi$ K$^+$. The analysis uses a data 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 61.6 fb$^{-1}$. The $f_\mathrm{s} / f_\mathrm{u}$ ratio is observed to depend on the B $p_\mathrm{T}$ and to be consistent with becoming asymptotically constant at large $p_\mathrm{T}$. No rapidity dependence is observed. The ratio of the B$^0$ to B$^+$ hadron production fractions, $f_\mathrm{d} / f_\mathrm{u}$, measured using the B$^0$$\to$ J$/\psi$ K$^{*0}$ decay channel, is found to be consistent with unity and independent of $p_\mathrm{T}$ and rapidity, as expected from isospin invariance.

Differential cross sections are measured for the standard model Higgs boson produced in association with vector bosons (W, Z) and decaying to a pair of b quarks. Measurements are performed within the framework of the simplified template cross sections. The analysis relies on the leptonic decays of the W and Z bosons, resulting in final states with 0, 1, or 2 electrons or muons. The Higgs boson candidates are either reconstructed from pairs of resolved b-tagged jets, or from single large distance parameter jets containing the particles arising from two b quarks. Proton-proton collision data at $\sqrt{s}$ = 13 TeV, collected by the CMS experiment in 2016-2018 and corresponding to a total integrated luminosity of 138 fb$^{-1}$, are analyzed. The inclusive signal strength, defined as the product of the observed production cross section and branching fraction relative to the standard model expectation, combining all analysis categories, is found to be $\mu$ = 1.15$^{+0.22}_{-0.20}$. This corresponds to an observed (expected) significance of 6.3 (5.6) standard deviations.

Two-particle correlation functions were measured for $\rm p\overline{p}$, $\rm p\overline{\Lambda}$, $\rm \overline{p}\Lambda$, and $\Lambda\overline{\Lambda}$ pairs in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV and $\sqrt{s_{\rm NN}}=5.02$ TeV recorded by the ALICE detector. From a simultaneous fit to all obtained correlation functions, real and imaginary components of the scattering lengths, as well as the effective ranges, were extracted for combined $\rm p\overline{\Lambda}$ and $\rm \overline{p}\Lambda$ pairs and, for the first time, for $\Lambda\overline{\Lambda}$ pairs. Effective averaged scattering parameters for heavier baryon-antibaryon pairs, not measured directly, are also provided. The results reveal similarly strong interaction between measured baryon-antibaryon pairs, suggesting that they all annihilate in the same manner at the same pair relative momentum $k^{*}$. Moreover, the reported significant non-zero imaginary part and negative real part of the scattering length provide motivation for future baryon-antibaryon bound state searches.

A search is presented for pairs of light pseudoscalar bosons, in the mass range from 4 to 15 GeV, produced from decays of the 125 GeV Higgs boson. The decay modes considered are final states that arise when one of the pseudoscalars decays to a pair of tau leptons, and the other one either into a pair of tau leptons or muons. The search is based on proton-proton collisions collected by the CMS experiment in 2016 at a center-of-mass energy of 13 TeV that correspond to an integrated luminosity of 35.9 fb${-1}$. The 2$\mu$2$\tau$ and 4$\tau$ channels are used in combination to constrain the product of the Higgs boson production cross section and the branching fraction into 4$\tau$ final state, $\sigma\mathcal{B}$, exploiting the linear dependence of the fermionic coupling strength of pseudoscalar bosons on the fermion mass. No significant excess is observed beyond the expectation from the standard model. The observed and expected upper limits at 95% confidence level on $\sigma\mathcal{B}$, relative to the standard model Higgs boson production cross section, are set respectively between 0.022 and 0.23 and between 0.027 and 0.19 in the mass range probed by the analysis.

The Lambda(b) differential production cross section and the cross section ratio anti-Lambda(b)/Lambda(b) are measured as functions of transverse momentum pt(Lambda(b)) and rapidity abs(y(Lambda(b))) in pp collisions at sqrt(s) = 7 TeV using data collected by the CMS experiment at the LHC. The measurements are based on Lambda(b) decays reconstructed in the exclusive final state J/Psi Lambda, with the subsequent decays J/Psi to an opposite-sign muon pair and Lambda to proton pion, using a data sample corresponding to an integrated luminosity of 1.9 inverse femtobarns. The product of the cross section times the branching ratio for Lambda(b) to J/Psi Lambda versus pt(Lambda(b)) falls faster than that of b mesons. The measured value of the cross section times the branching ratio for pt(Lambda(b)) > 10 GeV and abs(y(Lambda(b))) < 2.0 is 1.06 +/- 0.06 +/- 0.12 nb, and the integrated cross section ratio for anti-Lambda(b)/Lambda(b) is 1.02 +/- 0.07 +/- 0.09, where the uncertainties are statistical and systematic, respectively.

A search for exotic decays of the Higgs boson to a pair of light pseudoscalar particles a$_1$ is performed under the hypothesis that one of the pseudoscalars decays to a pair of opposite sign muons and the other decays to b$\overline{\mathrm{b}}$. Such signatures are predicted in a number of extensions of the standard model (SM), including next-to-minimal supersymmetry and two-Higgs-doublet models with an additional scalar singlet. The results are based on a data set of proton-proton collisions corresponding to an integrated luminosity of 35.9 fb$^{-1}$, accumulated with the CMS experiment at the CERN LHC in 2016 at a centre-of-mass energy of 13 TeV. No statistically significant excess is observed with respect to the SM backgrounds in the search region for pseudoscalar masses from 20 GeV to half of the Higgs boson mass. Upper limits at 95% confidence level are set on the product of the production cross section and branching fraction, $\sigma_{\mathrm{h}}\mathcal{B}$(h $\to$ a$_1$ a$_1$ $\to$ $\mu^+\mu^-\mathrm{b}\bar{\mathrm{b}}$), ranging from 5 to 33 fb, depending on the pseudoscalar mass. Corresponding limits on the branching fraction, assuming the SM prediction for $\sigma_{\mathrm{h}}$, are (1$-$7)$\times$ 10$^{-4}$.

A search for new massive particles decaying into a pair of Higgs bosons in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. Data were collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search is performed for resonances with a mass between 0.8 and 3.5 TeV using events in which one Higgs boson decays into a bottom quark pair and the other decays into two W bosons that subsequently decay into a lepton, a neutrino, and a quark pair. The Higgs boson decays are reconstructed with techniques that identify final state quarks as substructure within boosted jets. The data are consistent with standard model expectations. Exclusion limits are placed on the product of the cross section and branching fraction for generic spin-0 and spin-2 massive resonances. The results are interpreted in the context of radion and bulk graviton production in models with a warped extra spatial dimension. These are the best results to date from searches for an HH resonance decaying to this final state, and they are comparable to the results from searches in other channels for resonances with masses below 1.5 TeV.

A measurement is presented of the cross section for electroweak production of a Z boson and a photon in association with two jets (Z$\gamma$jj) in proton-proton collisions. The Z boson candidates are selected through their decay into a pair of electrons or muons. The process of interest, electroweak Z$\gamma$jj production, is isolated by selecting events with a large dijet mass and a large pseudorapidity gap between the two jets. The measurement is based on data collected at the CMS experiment at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The observed significance of the signal is 3.9 standard deviations, where a significance of 5.2 standard deviations is expected in the standard model. These results are combined with published results by CMS at $\sqrt{s} =$ 8 TeV, which leads to observed and expected respective significances of 4.7 and 5.5 standard deviations. From the 13 TeV data, a value is obtained for the signal strength of electroweak Z$\gamma$jj production and bounds are given on quartic vector boson interactions in the framework of dimension-eight effective field theory operators.

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).

A search is presented for the decay of the 125 GeV Higgs boson (H) to a pair of new light pseudoscalar bosons (a), followed by the prompt decay of each a boson to a bottom quark-antiquark pair, H $\to$ aa $\to$$\mathrm{b\bar{b}b\bar{b}}$. The analysis is performed using a data sample of proton-proton collisions collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. To reduce the background from standard model processes, the search requires the Higgs boson to be produced in association with a leptonically decaying W or Z boson. The analysis probes the production of new light bosons in a 15 $\lt$$m_\mathrm{a}$$\lt$ 60 GeV mass range. Assuming the standard model predictions for the Higgs boson production cross sections for pp $\to$ WH and ZH, model independent upper limits at 95% confidence level are derived for the branching fraction $\mathcal{B}$(H $\to$ aa $\to$ $\mathrm{b\bar{b}b\bar{b}}$). The combined WH and ZH observed upper limit on the branching fraction ranges from 1.10 for $m_\mathrm{a} =$ 20 GeV to 0.36 for $m_\mathrm{a} =$ 60 GeV, complementing other measurements in the $\mu\mu\tau\tau$, $\tau\tau\tau\tau$ and bb$\ell\ell$ ($\ell=$ $\mu$,$\tau$) channels.

A measurement of the production cross section of a Higgs boson with transverse momentum greater than 250 GeV is presented where the Higgs boson decays to a pair of $\tau$ leptons. It is based on proton-proton collision data collected by the CMS experiment at the CERN LHC at a center-of-mass energy of 13 TeV. The data sample corresponds to an integrated luminosity of 138 fb$^{-1}$. Because of the large transverse momentum of the Higgs boson the $\tau$ leptons from its decays are boosted and produced spatially close, with their decay products overlapping. Therefore, a dedicated algorithm was developed to reconstruct and identify them. The observed (expected) significance of the measured signal with respect to the standard model background-only hypothesis is 3.5 (2.2) standard deviations. The product of the production cross section and branching fraction is measured to be 1.64$^{+0.68}_{-0.54}$ times the standard model expectation. The fiducial differential production cross section is also measured as functions of the Higgs boson and leading jet transverse momenta. This measurement extends the probed large-transverse-momentum region beyond 600 GeV.

Transverse-momentum ($p_{\rm T}$) differential yields of electrons from semileptonic heavy-flavour hadron decays have been measured in the most central (0-10%) and in semi-central (20-40%) Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV. The corresponding production cross section in pp collisions has been measured at the same energy with substantially reduced systematic uncertainties with respect to previously published results. The modification of the yield in Pb-Pb collisions with respect to the expectation from an incoherent superposition of nucleon-nucleon collisions is quantified at mid-rapidity ($|y|$ $<$ 0.8) in the $p_{\rm T}$ interval 0.5-3 GeV/$c$ via the nuclear modification factor, $R_{\rm AA}$. This paper extends the $p_{\rm T}$ reach of the $R_{\rm AA}$ measurement towards significantly lower values with respect to a previous publication. In Pb-Pb collisions the $p_{\rm T}$-differential measurements of yields at low $p_{\rm T}$ are essential to investigate the scaling of heavy-flavour production with the number of binary nucleon-nucleon collisions. Heavy-quark hadronization, a collective expansion and even initial-state effects, such as the nuclear modification of the Parton Distribution Functions, are also expected to have a significant effect on the measured distribution.

A search for supersymmetric particles produced in the vector boson fusion topology in proton-proton collisions is presented. The search targets final states with one or zero leptons, large missing transverse momentum, and two jets with a large separation in rapidity. The data sample corresponds to an integrated luminosity of 35.9 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV collected in 2016 with the CMS detector at the LHC. The observed dijet invariant mass and lepton-neutrino transverse mass spectra are found to be consistent with the standard model predictions. Upper limits are set on the cross sections for chargino ($\widetilde\chi_1^\pm$) and neutralino ($\widetilde\chi_2^0$) production with two associated jets. For a compressed mass spectrum scenario in which the $\widetilde\chi_1^\pm$ and $\widetilde\chi_2^0$ decays proceed via a light slepton and the mass difference between the lightest neutralino $\widetilde\chi_1^0$ and the mass-degenerate particles $\widetilde\chi_1^\pm$ and $\widetilde\chi_2^0$ is 1 (30) GeV, the most stringent lower limit to date of 112 (215) GeV is set on the mass of these latter two particles.

Two related searches for phenomena beyond the standard model (BSM) are performed using events with hadronic jets and significant transverse momentum imbalance. The results are based on a sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment at the LHC in 2016-2018 and corresponding to an integrated luminosity of 137 fb$^{-1}$. The first search is inclusive, based on signal regions defined by the hadronic energy in the event, the jet multiplicity, the number of jets identified as originating from bottom quarks, and the value of the kinematic variable $M_\mathrm{T2}$ for events with at least two jets. For events with exactly one jet, the transverse momentum of the jet is used instead. The second search looks in addition for disappearing tracks produced by BSM long-lived charged particles that decay within the volume of the tracking detector. No excess event yield is observed above the predicted standard model background. This is used to constrain a range of BSM models that predict the following: the pair production of gluinos and squarks in the context of supersymmetry models conserving $R$-parity, with or without intermediate long-lived charginos produced in the decay chain; the resonant production of a colored scalar state decaying to a massive Dirac fermion and a quark; or the pair production of scalar and vector leptoquarks each decaying to a neutrino and a top, bottom, or light-flavor quark. In most of the cases, the results obtained are the most stringent constraints to date.

Decays of the Higgs boson into a Z boson and a J/$\psi$ or $\psi$(2S) meson are searched for in four-lepton final states with the CMS detector at the LHC. A data set of proton-proton collisions corresponding to an integrated luminosity of 138 fb$^{-1}$ is used. Using the same data set, decays of the Higgs and Z boson into quarkonium pairs are also searched for. An observation of such decays with this sample would indicate the presence of physics beyond the standard model. No evidence for these decays has been observed and upper limits at the 95% confidence level are placed on the corresponding branching fractions ($\mathcal{B}$). Assuming longitudinal polarization of the Higgs boson decay products, 95% confidence level observed upper limits for $\mathcal{B}$(H $\to$ J/$\psi$) and $\mathcal{B}$(H $\to$ Z$\psi$(2S)) are 1.9$\times$10$^{-3}$ and 6.6$\times$10$^{-3}$, respectively.

The production cross section of a top quark pair in association with a photon is measured in proton-proton collisions in the decay channel with two oppositely charged leptons (e$^\pm\mu^\mp$, e$^+$e$^-$, or $\mu^+\mu^-$). The measurement is performed using 138 fb$^{-1}$ of proton-proton collision data recorded by the CMS experiment at $\sqrt{s} =$ 13 TeV during the 2016-2018 data-taking period of the CERN LHC. A fiducial phase space is defined such that photons radiated by initial-state particles, top quarks, or any of their decay products are included. An inclusive cross section of 175.2 $\pm$ 2.5 (stat) $\pm$ 6.3 (syst) fb is measured in a signal region with at least one jet coming from the hadronization of a bottom quark and exactly one photon with transverse momentum above 20 GeV. Differential cross sections are measured as functions of several kinematic observables of the photon, leptons, and jets, and compared to standard model predictions. The measurements are also interpreted in the standard model effective field theory framework, and limits are found on the relevant Wilson coefficients from these results alone and in combination with a previous CMS measurement of the $\mathrm{t\bar{t}}\gamma$ production process using the lepton+jets final state.

Measurements are presented of the B$^0_\mathrm{S}$$\to$$\mu^+\mu^-$ branching fraction and effective lifetime, as well as results of a search for the B$^0$$\to$$\mu^+\mu^-$ decay in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the LHC. The analysis is based on data collected with the CMS detector in 2016-2018 corresponding to an integrated luminosity of 140 fb$^{-1}$. The branching fraction of the B$^0_\mathrm{S}$$\to$$\mu^+\mu^-$ decay and the effective B$^0_\mathrm{S}$ meson lifetime are the most precise single measurements to date. No evidence for the B$^0$$\to$$\mu^+\mu^-$ decay has been found. All results are found to be consistent with the standard model predictions and previous measurements.

A search for a new boson X is presented using CERN LHC proton-proton collision data collected by the CMS experiment at $\sqrt{s}$ = 13 TeV in 2016-2018, and corresponding to an integrated luminosity of 138 fb$^{-1}$. The resonance X decays into either a pair of Higgs bosons HH of mass 125 GeV or an H and a new spin-0 boson Y. One H subsequently decays to a pair of photons, and the second H or Y, to a pair of bottom quarks. The explored mass ranges of X are 260-1000 GeV and 300-1000 GeV, for decays to HH and to HY, respectively, with the Y mass range being 90-800 GeV. For a spin-0 X hypothesis, the 95% confidence level upper limit on the product of its production cross section and decay branching fraction is observed to be within 0.90-0.04 fb, depending on the masses of X and Y. The largest deviation from the background-only hypothesis with a local (global) significance of 3.8 (2.8) standard deviations is observed for X and Y masses of 650 and 90 GeV, respectively. The limits are interpreted using several models of new physics.

A test of lepton flavor universality in B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ and B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$ decays, as well as a measurement of differential and integrated branching fractions of a nonresonant B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$ decay are presented. The analysis is made possible by a dedicated data set of proton-proton collisions at $\sqrt{s}$ = 13 TeV recorded in 2018, by the CMS experiment at the LHC, using a special high-rate data stream designed for collecting about 10 billion unbiased b hadron decays. The ratio of the branching fractions $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) to $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}$e$^+$e$^-$) is determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the B$^\pm$$\to$ J/$\psi$K$^\pm$ with J/$\psi$$\to$$\mu^+\mu^-$ and e$^+$e$^-$ decays, which allow for significant cancellation of systematic uncertainties. The ratio $R$(K) is measured in the range 1.1 $\lt q^2 \lt$ 6.0 GeV$^2$, where $q$ is the invariant mass of the lepton pair, and is found to be $R$(K) = 0.78$^{+0.47}_{-0.23}$, in agreement with the standard model expectation $R$(K) $\approx$ 1. This measurement is limited by the statistical precision of the electron channel. The integrated branching fraction in the same $q^2$ range, $\mathcal{B}$(B$^{\pm}$$\to$ K$^{\pm}\mu^+\mu^-$) = (12.42 $\pm$ 0.68) $\times$ 10$^{-8}$, is consistent with the present world-average value and has a comparable precision.

A data sample of events from proton-proton collisions with at least two jets, and two isolated same-sign or three or more charged leptons, is studied in a search for signatures of new physics phenomena. The data correspond to an integrated luminosity of 137 fb$^{-1}$ at a center-of-mass energy of 13 TeV, collected in 2016-2018 by the CMS experiment at the LHC. The search is performed using a total of 168 signal regions defined using several kinematic variables. The properties of the events are found to be consistent with the expectations from standard model processes. Exclusion limits at 95% confidence level are set on cross sections for the pair production of gluinos or squarks for various decay scenarios in the context of supersymmetric models conserving or violating R parity. The observed lower mass limits are as large as 2.1 TeV for gluinos and 0.9 TeV for top and bottom squarks. To facilitate reinterpretations, model-independent limits are provided in a set of simplified signal regions.

A search for the lepton-flavor violating decay of the Higgs boson and potential additional Higgs bosons with a mass in the range 110-160 GeV to an e$^{\pm}\mu^{\mp}$ pair is presented. The search is performed with a proton-proton collision dataset at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. No excess is observed for the Higgs boson. The observed (expected) upper limit on the e$^{\pm}\mu^{\mp}$ branching fraction for it is determined to be 4.4 (4.7) $\times$ 10$^{-5}$ at 95% confidence level, the most stringent limit set thus far from direct searches. The largest excess of events over the expected background in the full mass range of the search is observed at an e$^{\pm}\mu^{\mp}$ invariant mass of approximately 146 GeV with a local (global) significance of 3.8 (2.8) standard deviations.

The first search is presented for vector-like leptons (VLLs) in the context of the "4321 model", an ultraviolet-complete model with the potential to explain existing B physics measurements that are in tension with standard model predictions. The analyzed data, corresponding to an integrated luminosity of 96.5 fb$^{-1}$, were recorded in 2017 and 2018 with the CMS detector at the LHC in proton-proton collisions at $\sqrt{s}$ = 13 TeV. Final states with ${\geq}$ 3 b-tagged jets and two third-generation leptons ($\tau\tau$, $\tau\nu_\tau$, or $\nu_\tau\nu_\tau$) are considered. Upper limits are derived on the VLL production cross section in the VLL mass range 500-1050 GeV. The maximum likelihood fit prefers the presence of signal at the level of 2.8 standard deviations, for a representative VLL mass point of 600 GeV. As a consequence, the observed upper limits are approximately double the expected limits.