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 (below 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 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.

A search for long-lived heavy neutrinos (N) in the decays of B mesons produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data sample corresponds to an integrated luminosity of 41.6 fb$^{-1}$ collected in 2018 by the CMS experiment at the CERN LHC, using a dedicated data stream that enhances the number of recorded events containing B mesons. The search probes heavy neutrinos with masses in the range 1 $\lt$$m_\mathrm{N}$$\lt$ 3 GeV and decay lengths in the range 10$^{-2}$$\lt$$c\tau_\mathrm{N}$$\lt$ 10$^{4}$ mm, where $\tau_\mathrm{N}$ is the N proper mean lifetime. Signal events are defined by the signature B $\to$$\ell_\mathrm{B}$NX; N $\to$$\ell^{\pm} \pi^{\mp}$, where the leptons $\ell_\mathrm{B}$ and $\ell$ can be either a muon or an electron, provided that at least one of them is a muon. The hadronic recoil system, X, is treated inclusively and is not reconstructed. No significant excess of events over the standard model background is observed in any of the $\ell^{\pm} \pi^{\mp}$ invariant mass distributions. Limits at 95% confidence level on the sum of the squares of the mixing amplitudes between heavy and light neutrinos, $\vert V_\mathrm{N}\vert^2$, and on $c\tau$ are obtained in different mixing scenarios for both Majorana and Dirac-like N particles. The most stringent upper limit $\vert V_\mathrm{N}\vert^2$ $\lt$ 2.0 $\times$ 10$^{-5}$ is obtained at $m_\mathrm{N}$ = 1.95 GeV for the Majorana case where N mixes exclusively with muon neutrinos. The limits on $\vert V_\mathrm{N}\vert^2$ for masses 1 $\lt$ $m_\mathrm{N}$ $\lt$ 1.7 GeV are the most stringent from a collider experiment to date.

A search for the pair production of the lightest supersymmetric partner of the top quark, the top squark ($\tilde{\mathrm{t}}_1$), is presented. The search targets the four-body decay of the $\tilde{\mathrm{t}}_1$, which is preferred when the mass difference between the top squark and the lightest supersymmetric particle is smaller than the mass of the W boson. This decay mode consists of a bottom quark, two other fermions, and the lightest neutralino ($\tilde{\chi}^0_1$), which is assumed to be the lightest supersymmetric particle. The data correspond to an integrated luminosity of 138 fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13 TeV collected by the CMS experiment at the CERN LHC. Events are selected using the presence of a high-momentum jet, an electron or muon with low transverse momentum, and a significant missing transverse momentum. The signal is selected based on a multivariate approach that is optimized for the difference between $m(\tilde{\mathrm{t}}_1)$ and $m(\tilde{\chi}^0_1)$. The contribution from leading background processes is estimated from data. No significant excess is observed above the expectation from standard model processes. The results of this search exclude top squarks at 95% confidence level for masses up to 480 and 700 GeV for $m(\tilde{\mathrm{t}}_1) - m(\tilde{\chi}^0_1$) = 10 and 80 GeV, respectively.

The structure of nucleons is multidimensional and depends on the transverse momenta, spatial geometry, and polarization of the constituent partons. Such a structure can be studied using high-energy photons produced in ultraperipheral heavy-ion collisions. The first measurement of the azimuthal angular correlations of exclusively produced events with two jets in photon-lead interactions at large momentum transfer is presented, a process that is considered to be sensitive to the underlying nuclear gluon polarization. This study uses a data sample of ultraperipheral lead-lead collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, corresponding to an integrated luminosity of 0.38 nb$^{-1}$, collected with the CMS experiment at the LHC. The measured second harmonic of the correlation between the sum and difference of the two jet momenta is found to be positive, and rising, as the dijet momentum increases. A well-tuned model that has been successful at describing a wide range of proton scattering data from the HERA experiments fails to describe the observed correlations, suggesting the presence of gluon polarization effects.

A search for a heavy pseudoscalar Higgs boson, A, decaying to a 125 GeV Higgs boson h and a Z boson is presented. The h boson is identified via its decay to a pair of tau leptons, while the Z boson is identified via its decay to a pair of electrons or muons. The search targets the production of the A boson via the gluon-gluon fusion process, gg $\to$ A, and in association with bottom quarks, $\mathrm{b\bar{b}}$A. The analysis uses a data sample corresponding to an integrated luminosity of 138 fb$^{-1}$ collected with the CMS detector at the CERN LHC in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV. Constraints are set on the product of the cross sections of the A production mechanisms and the A $\to$ Zh decay branching fraction. The observed (expected) upper limit at 95% confidence level ranges from 0.049 (0.060) pb to 1.02 (0.79) pb for the gg $\to$ A process and from 0.053 (0.059) pb to 0.79 (0.61) pb for the $\text{b}\bar{\text{b}}$A process in the probed range of the A boson mass, $m_\text{A}$, from 225 GeV to 1 TeV. The results of the search are used to constrain parameters within the ${\text{M}_{\text{h,EFT}}^{\text{125}}}$ benchmark scenario of the minimal supersymmetric extension of the standard model. Values of $\tanβ$ below 2.2 are excluded in this scenario at 95% confidence level for all $m_\text{A}$ values in the range from 225 to 350 GeV.

Measurements of light-by-light scattering (LbL, $γγ$$\to$$γγ$) and the Breit-Wheeler process (BW, $γγ$$\to$$\mathrm{e^+e^-}$) are reported in ultraperipheral PbPb collisions at a centre-of-mass energy per nucleon pair of 5.02 TeV. The data sample, corresponding to an integrated luminosity of 1.7 nb$^{-1}$, was collected by the CMS experiment at the CERN LHC in 2018. Events with an exclusively produced $γγ$ or $\mathrm{e^+e^-}$ pair with invariant masses $m^{γγ,\mathrm{ee}}$$>$ 5 GeV, along with other fiducial criteria, are selected. The measured BW fiducial production cross section, $σ_\text{fid}$($γγ$$\to$$\mathrm{e^+e^-}$) = 263.5 $\pm$ 1.8 (stat) $\pm$ 17.8 (syst) $μ$b, as well as the differential distributions for various kinematic observables, are in agreement with leading-order quantum electrodynamics predictions complemented with final-state photon radiation. The measured differential BW cross sections allow discriminating between different theoretical descriptions of the photon flux of the lead ion. In the LbL final state, 26 exclusive diphoton candidate events are observed compared with 12.0 $\pm$ 2.9 expected for the background. Combined with previous results, the observed significance of the LbL signal with respect to the background-only hypothesis is above five standard deviations. The measured fiducial LbL scattering cross section, $σ_\text{fid} (γγ$$\to$$γγ)$ = 107 $\pm$ 24 (stat) $\pm$ 13 (syst) nb, is in agreement with next-to-leading-order predictions. Limits on the production of axion-like particles coupled to photons are set over the mass range 5-100 GeV, including the most stringent limits to date in the range of 5-10 GeV.

The production cross sections of B$^0_\mathrm{s}$ and B$^+$ mesons are reported in proton-proton (pp) collisions recorded by the CMS experiment at the CERN LHC with a center-of-mass energy of 5.02 TeV. The data sample corresponds to an integrated luminosity of 302 pb$^{-1}$. The cross sections are based on measurements of the B$^0_\mathrm{s}$$\to$ J/$ψ(μ^+μ^-)ϕ$(1020) (K$^+$K$^-$) and B$^+$$\to$ J/$ψ(μ^+μ^-)$K$^+$ decay channels. Results are presented in the transverse momentum ($p_\mathrm{T}$) range 7-50 GeV/$c$ and the rapidity interval $\lvert y \rvert$$\lt$ 2.4 for the B mesons. The measured $p_\mathrm{T}$-differential cross sections of B$^+$ and B$^0_\mathrm{s}$ in pp collisions are well described by fixed-order plus next-to-leading logarithm perturbative quantum chromodynamics calculations. Using previous PbPb collision measurements at the same nucleon-nucleon center-of-mass energy, the nuclear modification factors, $R_\mathrm{AA}$, of the B mesons are determined. For $p_\mathrm{T}$$\lt$ 10 GeV/$c$, both mesons are found to be suppressed in PbPb collisions (with $R_\mathrm{AA}$ values significantly below unity), with less suppression observed for the B$^0_\mathrm{s}$ mesons. In this $p_\mathrm{T}$ range, the $R_\mathrm{AA}$ values for the B$^+$ mesons are consistent with those for inclusive charged hadrons and D$^0$ mesons. Below 10 GeV/$c$, both B$^+$ and B$^0_\mathrm{s}$ are found to be less suppressed than either inclusive charged hadrons or D$^0$ mesons, with the B$^0_\mathrm{s}$$R_\mathrm{AA}$ value consistent with unity. The $R_\mathrm{AA}$ values found for the B$^+$ and B$^0_\mathrm{s}$ are compared to theoretical calculations, providing constraints on the mechanism of bottom quark energy loss and hadronization in the quark-gluon plasma, the hot and dense matter created in ultrarelativistic heavy ion collisions.

Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

A new algorithm is presented to discriminate reconstructed hadronic decays of tau leptons ($\tau_\mathrm{h}$) that originate from genuine tau leptons in the CMS detector against $\tau_\mathrm{h}$ candidates that originate from quark or gluon jets, electrons, or muons. The algorithm inputs information from all reconstructed particles in the vicinity of a $\tau_\mathrm{h}$ candidate and employs a deep neural network with convolutional layers to efficiently process the inputs. This algorithm leads to a significantly improved performance compared with the previously used one. For example, the efficiency for a genuine $\tau_\mathrm{h}$ to pass the discriminator against jets increases by 10-30% for a given efficiency for quark and gluon jets. Furthermore, a more efficient $\tau_\mathrm{h}$ reconstruction is introduced that incorporates additional hadronic decay modes. The superior performance of the new algorithm to discriminate against jets, electrons, and muons and the improved $\tau_\mathrm{h}$ reconstruction method are validated with LHC proton-proton collision data at $\sqrt{s} =$ 13 TeV.