Signals consistent with the B$^+_\mathrm{c}$(2S) and B$^{*+}_\mathrm{c}$(2S) states are observed in proton-proton collisions at $\sqrt{s} =$ 13 TeV, in an event sample corresponding to an integrated luminosity of 140 fb$^{-1}$, collected by the CMS experiment during the 2016, 2017, and 2018 LHC running periods. These excited $\bar{\mathrm{b}}$c states are observed in the B$^+_\mathrm{c}\pi^+\pi^-$ invariant mass spectrum, with the ground state B$^+_\mathrm{c}$ reconstructed through its decay to J/$\psi\,\pi^+$. The two states are well resolved from each other and are observed with a significance exceeding five standard deviations. The mass of the B$^+_\mathrm{c}$(2S) meson is measured to be 6871.0 $\pm$ 1.2 (stat) $\pm$ 0.8 (syst) $\pm$ 0.8 (B$^+_\mathrm{c}$) MeV, where the last term corresponds to the uncertainty in the world-average B$^+_\mathrm{c}$ mass.

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

The J/$\psi$$\to$$\mu^+\mu^-\mu^+\mu^-$ decay has been observed with a statistical significance in excess of five standard deviations. The analysis is based on an event 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 33.6 fb${-1}$. Normalizing to the J/$\psi$$\to$$\mu^+\mu^-$ decay mode leads to a branching fraction [10.1$^{+3.3}_{-2.7}$ (stat) $\pm$ 0.4 (syst) ]$\times$ 10$^{-7}$, a value that is consistent with the standard model prediction.

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.

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 the lepton flavor violating $\tau$$\to$ 3$\mu$ decay is performed using proton-proton collision events at a center-of-mass energy of 13 TeV collected by the CMS experiment at the LHC in 2017-2018, corresponding to an integrated luminosity of 97.7 fb$^{-1}$. Tau leptons produced in both heavy-flavor hadron and W boson decays are exploited in the analysis. No evidence for the decay is observed. The results of this search are combined with an earlier null result based on data collected in 2016 to obtain a total integrated luminosity of 131 fb$^{-1}$. The observed (expected) upper limits on the branching fraction $\mathcal{B}$($\tau$$\to$ 3$\mu$) at confidence levels of 90 and 95% are 2.9 $\times$ 10$^{-8}$ (2.4 $\times$ 10$^{-8}$) and 3.6 $\times$ 10$^{-8}$ (3.0 $\times$ 10$^{-8}$), respectively.

The observation of the $\Lambda_\mathrm{b}^0 \to$J/$\psi \Lambda \phi$ decay is reported using proton-proton collision data collected at $\sqrt{s} =$ 13 TeV by the CMS experiment at the LHC in 2018, corresponding to an integrated luminosity of 60 fb$^{-1}$. The ratio of the branching fractions $\mathcal{B}(\Lambda_\mathrm{b}^0 \to$J/$\psi \Lambda \phi)/\mathcal{B}(\Lambda_\mathrm{b}^0\to\psi \Lambda)$ is measured to be (8.26$\pm$0.90 (stat) $\pm$ 0.68 (syst) $\pm$ 0.11 $(\mathcal{B}))\times $10$^{-2}$, where the first uncertainty is statistical, the second is systematic, and the last uncertainty reflects the uncertainties in the world-average branching fractions of $\phi$ and $\psi$(2S) decays to the reconstructed final states.

The $\chi_\mathrm{b1}$(3P) and $\chi_\mathrm{b2}$(3P) states are observed through their $\Upsilon$(3S) $\gamma$ decays, using an event sample of proton-proton collisions collected by the CMS experiment at the CERN LHC. The data were collected at a center-of-mass energy of 13 TeV and correspond to an integrated luminosity of 80.0 fb$^{-1}$. The $\Upsilon$(3S) mesons are identified through their dimuon decay channel, while the low-energy photons are detected after converting to e$^+$e$^-$ pairs in the silicon tracker, leading to a $\chi_\mathrm{b}$(3P) mass resolution of 2.2 MeV. This is the first time that the $J =$ 1 and 2 states are well resolved and their masses individually measured: 10$\,$513.42 $\pm$ 0.41 (stat) $\pm$ 0.18 (syst) MeV and 10$\,$524.02 $\pm$ 0.57 (stat) $\pm$ 0.18 (syst) MeV; they are determined with respect to the world-average value of the $\Upsilon$(3S) mass, which has an uncertainty of 0.5 MeV. The mass splitting is measured to be 10.60 $\pm$ 0.64 (stat) $\pm$ 0.17 (syst) MeV.

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.

Using a data sample of $\sqrt{s} =$ 13 TeV proton-proton collisions collected by the CMS experiment at the LHC in 2017 and 2018 with an integrated luminosity of 103 fb$^{-1}$, the B$^0$$\to$$\psi$(2S)K$^0_\mathrm{S}$ and B$^0_\mathrm{S}$$\to$$\psi$(2S)K$^0_\mathrm{S}\pi^+\pi^-$ decays are observed with significances exceeding 5 standard deviations. The resulting branching fraction ratios, measured for the first time, correspond to $\mathcal{B}$(B$^0_\mathrm{S}$$\to$$\psi$(2S)K$^0_\mathrm{S}$) / $\mathcal{B}$(B$^0$$\to$$\psi$(2S)K$^0_\mathrm{S}$) = (3.33 $\pm$ 0.69 (stat) $\pm$ 0.11 (syst) $\pm$ 0.34 ($f_\mathrm{s} / f_\mathrm{d}$)) $\times$ 10$^{-2}$ and $\mathcal{B}$(B$^0$$\to$$\psi$(2S)K$^0_\mathrm{S}\pi^+\pi^-$) / $\mathcal{B}$(B$^0$$\to$$\psi$(2S)K$^0_\mathrm{S}$) = 0.480 $\pm$ 0.013 (stat) $\pm$ 0.032 (syst), where the last uncertainty in the first ratio is related to the uncertainty in the ratio of production cross sections of B$^0_\mathrm{s}$ and B$^0$ mesons, $f_\mathrm{s} / f_\mathrm{d}$.