A search for physics beyond the standard model (SM) in the final state with a hadronically decaying tau lepton and a neutrino is presented. This analysis is based on data recorded by the CMS experiment from proton-proton collisions at a center-of-mass energy of 13 TeV at the LHC, corresponding to a total integrated luminosity of 138 fb$^{=1}$. The transverse mass spectrum is analyzed for the presence of new physics. No significant deviation from the SM prediction is observed. Limits are set on the production cross section of a W' boson decaying into a tau lepton and a neutrino. Lower limits are set on the mass of the sequential SM-like heavy charged vector boson and the mass of a quantum black hole. Upper limits are placed on the couplings of a new boson to the SM fermions. Constraints are put on a nonuniversal gauge interaction model and an effective field theory model. For the first time, upper limits on the cross section of $t$-channel leptoquark (LQ) exchange are presented. These limits are translated into exclusion limits on the LQ mass and on its coupling in the $t$-channel. The sensitivity of this analysis extends into the parameter space of LQ models that attempt to explain the anomalies observed in B meson decays. The limits presented for the various interpretations are the most stringent to date. Additionally, a model-independent limit is provided.

A search for low-mass dilepton resonances in Higgs boson decays is conducted in the four-lepton final state. The decay is assumed to proceed via a pair of beyond the standard model particles, or one such particle and a Z boson. The search uses proton-proton collision data collected with the CMS detector at the CERN LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$, at a center-of-mass energy $\sqrt{s} =$ 13 TeV. No significant deviation from the standard model expectation is observed. Upper limits at 95% confidence level are set on model-independent Higgs boson decay branching fractions. Additionally, limits on dark photon and axion-like particle production, based on two specific models, are reported.

Proton-proton interactions resulting in final states with two photons are studied in a search for the signature of flavor-changing neutral current interactions of top quarks (t) and Higgs bosons (H). The analysis is based on data collected at a center-of-mass energy of 13 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. No significant excess above the background prediction is observed. Upper limits on the branching fractions ($\mathcal{B}$) of the top quark decaying to a Higgs boson and an up (u) or charm quark (c) are derived through a binned fit to the diphoton invariant mass spectrum. The observed (expected) 95% confidence level upper limits are found to be 0.019 (0.031)% for $\mathcal B$(t $\to$ Hu) and 0.073 (0.051)% for $\mathcal{B}$(t $\to$ Hc). These are the strictest upper limits yet determined.

A search is presented for long-lived particles produced in pairs in proton-proton collisions at the LHC operating at a center-of-mass energy of 13 TeV. The data were collected with the CMS detector during the period from 2015 through 2018, and correspond to a total integrated luminosity of 140 fb$^{-1}$. This search targets pairs of long-lived particles with mean proper decay lengths between 0.1 and 100 mm, each of which decays into at least two quarks that hadronize to jets, resulting in a final state with two displaced vertices. No significant excess of events with two displaced vertices is observed. In the context of $R$-parity violating supersymmetry models, the pair production of long-lived neutralinos, gluinos, and top squarks is excluded at 95% confidence level for cross sections larger than 0.08 fb, masses between 800 and 3000 GeV, and mean proper decay lengths between 1 and 25 mm.

A search for a heavy resonance decaying into a top quark and a W boson in proton-proton collisions at $\sqrt{s} =$ 13 TeV is presented. The data analyzed were recorded with the CMS detector at the LHC and correspond to an integrated luminosity of 138 fb$^{-1}$. The top quark is reconstructed as a single jet and the W boson, from its decay into an electron or muon and the corresponding neutrino. A top quark tagging technique based on jet clustering with a variable distance parameter and simultaneous jet grooming is used to identify jets from the collimated top quark decay. The results are interpreted in the context of two benchmark models, where the heavy resonance is either an excited bottom quark b$^*$ or a vector-like quark B. A statistical combination with an earlier search by the CMS Collaboration in the all-hadronic final state is performed to place upper cross section limits on these two models. The new analysis extends the lower range of resonance mass probed from 1.4 down to 0.7 TeV. For left-handed, right-handed, and vector-like couplings, b$^*$ masses up to 3.0, 3.0, and 3.2 TeV are excluded at 95% confidence level, respectively. The observed upper limits represent the most stringent constraints on the b$^*$ model to date.

The collective behavior of K$^0_\mathrm{S}$ and $\Lambda/\bar{\Lambda}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($v_2$) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20\GeV is present. The strange hadron $v_2$ values extracted in \pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

A combination of measurements of the inclusive top-quark pair production cross-section performed by ATLAS and CMS in proton-proton collisions at centre-of-mass energies of 7 and 8 TeV at the LHC is presented. The cross-sections are obtained using top-quark pair decays with an opposite-charge electron-muon pair in the final state and with data corresponding to an integrated luminosity of about 5 fb$^{-1}$ at $\sqrt{s}=7$ TeV and about 20 fb$^{-1}$ at $\sqrt{s}=8$ TeV for each experiment. The combined cross-sections are determined to be $178.5 \pm 4.7$ pb at $\sqrt{s}=7$ TeV and $243.3^{+6.0}_{-5.9}$ pb at $\sqrt{s}=8$ TeV with a correlation of 0.41, using a reference top-quark mass value of 172.5 GeV. The ratio of the combined cross-sections is determined to be $R_{8/7}= 1.363\pm 0.032$. The combined measured cross-sections and their ratio agree well with theory calculations using several parton distribution function (PDF) sets. The values of the top-quark pole mass (with the strong coupling fixed at 0.118) and the strong coupling (with the top-quark pole mass fixed at 172.5 GeV) are extracted from the combined results by fitting a next-to-next-to-leading-order plus next-to-next-to-leading-log QCD prediction to the measurements. Using a version of the NNPDF3.1 PDF set containing no top-quark measurements, the results obtained are $m_t^\text{pole} = 173.4^{+1.8}_{-2.0}$ GeV and $\alpha_\text{s}(m_Z)= 0.1170^{+ 0.0021}_{-0.0018}$.

Pseudorapidity (eta) distributions of charged particles produced in proton-proton collisions at a centre-of-mass energy of 8 TeV are measured in the ranges abs(eta) < 2.2 and 5.3 < abs(eta) < 6.4 covered by the CMS and TOTEM detectors, respectively. The data correspond to an integrated luminosity of 45 inverse microbarns. Measurements are presented for three event categories. The most inclusive category is sensitive to 91-96% of the total inelastic proton-proton cross section. The other two categories are disjoint subsets of the inclusive sample that are either enhanced or depleted in single diffractive dissociation events. The data are compared to models used to describe high-energy hadronic interactions. None of the models considered provide a consistent description of the measured distributions.

A search for physics beyond the standard model in events with at least three charged leptons (electrons or muons) is presented. The data sample corresponds to an integrated luminosity of 137 fb$^{-1}$ of proton-proton collisions at $\sqrt{s} =$ 13 TeV, collected with the CMS detector at the LHC in 2016-2018. The two targeted signal processes are pair production of type-III seesaw heavy fermions and production of a light scalar or pseudoscalar boson in association with a pair of top quarks. The heavy fermions may be manifested as an excess of events with large values of leptonic transverse momenta or missing transverse momentum. The light scalars or pseudoscalars may create a localized excess in the dilepton mass spectra. The results exclude heavy fermions of the type-III seesaw model for masses below 880 GeV at 95% confidence level in the scenario of equal branching fractions to each lepton flavor. This is the most restrictive limit on the flavor-democratic scenario of the type-III seesaw model to date. Assuming a Yukawa coupling of unit strength to top quarks, branching fractions of new scalar (pseudoscalar) bosons to dielectrons or dimuons above 0.004 (0.03) and 0.04 (0.03) are excluded at 95% confidence level for masses in the range 15-75 and 108-340 GeV, respectively. These are the first limits in these channels on an extension of the standard model with scalar or pseudoscalar particles.

A search for narrow and broad resonances with masses greater than 1.8 TeV decaying to a pair of jets is presented. The search uses proton-proton collision data at $\sqrt{s} =$ 13 TeV collected at the LHC, corresponding to an integrated luminosity of 137 fb$^{-1}$. The background arising from standard model processes is predicted with the fit method used in previous publications and with a new method. The dijet invariant mass spectrum is well described by both data-driven methods, and no significant evidence for the production of new particles is observed. Model independent upper limits are reported on the production cross sections of narrow resonances, and broad resonances with widths up to 55% of the resonance mass. Limits are presented on the masses of narrow resonances from various models: string resonances, scalar diquarks, axigluons, colorons, excited quarks, color-octet scalars, W' and Z' bosons, Randall-Sundrum gravitons, and dark matter mediators. The limits on narrow resonances are improved by 200 to 800 GeV relative to those reported in previous CMS dijet resonance searches. The limits on dark matter mediators are presented as a function of the resonance mass and width, and on the associated coupling strength as a function of the mediator mass. These limits exclude at 95% confidence level a dark matter mediator with a mass of 1.8 TeV and width 1% of its mass or higher, up to one with a mass of 4.8 TeV and a width 45% of its mass or higher.

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 a narrow resonance with a mass between 350 and 700 GeV, and decaying into a pair of jets, is performed using proton-proton collision events containing at least three jets. The data sample corresponds to an integrated luminosity of 18.3 fb$^{-1}$ recorded at $\sqrt{s} =$ 13 TeV with the CMS detector. Data are collected with a technique known as "data scouting", in which the events are reconstructed, selected, and recorded at a high rate in a compact form by the high-level trigger. The three-jet final state provides sensitivity to lower resonance masses than in previous searches using the data scouting technique. The spectrum of the dijet invariant mass, calculated from the two jets with the largest transverse momenta in the event, is used to search for a resonance. No significant excess over a smoothly falling background is found. Limits at 95% confidence level are set on the production cross section of a narrow dijet resonance and compared with the cross section of a vector dark matter mediator coupling to dark matter particles and quarks. Translating to a model where the narrow vector resonance interacts only with quarks, upper limits on this coupling range between 0.10 and 0.15, depending on the resonance mass. These results represent the most stringent upper limits in the mass range between 350 and 450 GeV obtained with a flavor-inclusive dijet resonance search.

A search is presented for a narrow resonance decaying to a pair of oppositely charged muons using $\sqrt{s} =$ 13 TeV proton-proton collision data recorded at the LHC. In the 45-75 and 110-200 GeV resonance mass ranges, the search is based on conventional triggering and event reconstruction techniques. In the 11.5-45 GeV mass range, the search uses data collected with dimuon triggers with low transverse momentum thresholds, recorded at high rate by storing a reduced amount of trigger-level information. The data correspond to integrated luminosities of 137 fb$^{-1}$ and 96.6 fb$^{-1}$ for conventional and high-rate triggering, respectively. No significant resonant peaks are observed in the probed mass ranges. The search sets the most stringent constraints to date on a dark photon in the $\sim$30-75 and 110-200 GeV mass ranges.

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 measurement of the mass of the Higgs boson in the diphoton decay channel is presented. This analysis is based on 35.9 fb$^{-1}$ of proton-proton collision data collected during the 2016 LHC running period, with the CMS detector at a center-of-mass energy of 13 TeV. A refined detector calibration and new analysis techniques have been used to improve the precision of this measurement. The Higgs boson mass is measured to be $m_\mathrm{H} =$ 125.78 $\pm$ 0.26 GeV. This is combined with a measurement of $m_\mathrm{H}$ already performed in the H $\to$ ZZ $\to$ 4$\ell$ decay channel using the same data set, giving $m_\mathrm{H} =$ 125.46 $\pm$ 0.16 GeV. This result, when further combined with an earlier measurement of $m_\mathrm{H}$ using data collected in 2011 and 2012 with the CMS detector, gives a value for the Higgs boson mass of $m_\mathrm{H} =$ 125.38 $\pm$ 0.14 GeV. This is currently the most precise measurement of the mass of the Higgs boson.

The parton-level top quark (t) forward-backward asymmetry and the anomalous chromoelectric ($\hat{d}_\mathrm{t}$) and chromomagnetic ($\hat{\mu}_\mathrm{t}$) moments have been measured using LHC pp collisions at a center-of-mass energy of 13 TeV, collected in the CMS detector in a data sample corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The linearized variable $A_\mathrm{FB}^{(1)}$ is used to approximate the asymmetry. Candidate $\mathrm{t\bar{t}}$ events decaying to a muon or electron and jets in final states with low and high Lorentz boosts are selected and reconstructed using a fit of the kinematic distributions of the decay products to those expected for $\mathrm{t\bar{t}}$ final states. The values found for the parameters are $A_\mathrm{FB}^{(1)} =$ 0.048 $^{+0.095}_{-0.087}$ (stat) $^{+0.020}_{-0.029}$ (syst), $\hat{\mu}_\mathrm{t} =-$ 0.024 $^{+0.013}_{-0.009}$ (stat) $^{+0.016}_{-0.011}$ (syst), and a limit is placed on the magnitude of $|\hat{d}_\mathrm{t}|$ $<$ 0.03 at 95% confidence level.

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 search for beyond the standard model spin-0 bosons, $\phi$, that decay into pairs of electrons, muons, or tau leptons is presented. The search targets the associated production of such bosons with a W or Z gauge boson, or a top quark-antiquark pair, and uses events with three or four charged leptons, including hadronically decaying tau leptons. The proton-proton collision data set used in the analysis was collected at the LHC from 2016 to 2018 at a center-of-mass energy of 13 TeV, and corresponds to an integrated luminosity of 138 fb$^{-1}$. The observations are consistent with the predictions from standard model processes. Upper limits are placed on the product of cross sections and branching fractions of such new particles over the mass range of 15 to 350 GeV with scalar, pseudoscalar, or Higgs-boson-like couplings, as well as on the product of coupling parameters and branching fractions. Several model-dependent exclusion limits are also presented. For a Higgs-boson-like $\phi$ model, limits are set on the mixing angle of the Higgs boson with the $\phi$ boson. For the associated production of a $\phi$ boson with a top quark-antiquark pair, limits are set on the coupling to top quarks. Finally, limits are set for the first time on a fermiophilic dilaton-like model with scalar couplings and a fermiophilic axion-like model with pseudoscalar couplings.

Measurements of the charge-dependent two-particle angular correlation function in proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at$\sqrt{s_\mathrm{NN}}$ = 5.02 TeV are reported. The pPb and PbPb datasets correspond to integrated luminosities of 186\nbinv and 0.607 nb$^{-1}$, respectively, and were collected using the CMS detector at the CERN LHC. The charge-dependent correlations are characterized by balance functions of same- and opposite-sign particle pairs. The balance functions, which contain information about the creation time of charged particle pairs and the development of collectivity, are studied as functions of relative pseudorapidity ($\Delta \eta$) and relative azimuthal angle ($\Delta \phi$), for various multiplicity and transverse momentum ($p_\mathrm{T}$) intervals. A multiplicity dependence of the balance function is observed in $\Delta \eta$ and $\Delta \phi$ for both systems. The width of the balance functions decreases towards high-multiplicity collisions in the momentum region $\lt$2 GeV, for pPb and PbPb results. No multiplicity dependence is observed at higher transverse momentum. The data are compared with HYDJET, HIJING and AMPT generator predictions, none of which capture completely the multiplicity dependence seen in the data.

Quasireal photons exchanged in relativistic heavy ion interactions are powerful probes of the gluonic structure of nuclei. The coherent J/$\psi$ photoproduction cross section in ultraperipheral lead-lead collisions is measured as a function of photon-nucleus center-of-mass energies per nucleon (W$^\text{Pb}_{\gamma\text{N}}$), over a wide range of 40 $\lt$ W$^\text{Pb}_{\gamma\text{N}}$$\lt$ 400 GeV. Results are obtained using data at the nucleon-nucleon center-of-mass energy of 5.02 TeV collected by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of 1.52 nb$^{-1}$. The cross section is observed to rise rapidly at low W$^\text{Pb}_{\gamma\text{N}}$, and plateau above W$^\text{Pb}_{\gamma\text{N}}$$\approx$ 40 GeV, up to 400 GeV, a new regime of small Bjorken-$x$ ($\approx$ 6 $\times$ 10$^{-5}$) gluons being probed in a heavy nucleus. The observed energy dependence is not predicted by current quantum chromodynamic models.

The first measurement of the top quark pair ($\mathrm{t\bar{t}}$) production cross section in proton-proton collisions at $\sqrt{s}$ = 13.6 TeV is presented. Data recorded with the CMS detector at the CERN LHC in Summer 2022, corresponding to an integrated luminosity of 1.21 fb$^{-1}$, are analyzed. Events are selected with one or two charged leptons (electrons or muons) and additional jets. A maximum likelihood fit is performed in event categories defined by the number and flavors of the leptons, the number of jets, and the number of jets identified as originating from b quarks. An inclusive $\mathrm{t\bar{t}}$ production cross section of 881 $\pm$ 23 (stat+syst) $\pm$ 20 (lumi) pb is measured, in agreement with the standard model prediction of 924 $^{+32}_{-40}$ pb.

The production of prompt $\Lambda^+_\mathrm{c}$ baryons is measured via the exclusive decay channel $\Lambda^+_\mathrm{c}$$\to$ pK$^-\pi^+$ at a center-of-mass energy per nucleon pair of 5.02 TeV, using proton-proton (pp) and lead-lead (PbPb) collision data collected by the CMS experiment at the CERN LHC. The pp and PbPb data were obtained in 2017 and 2018 with integrated luminosities of 252 and 0.607 nb$^{-1}$, respectively. The measurements are performed within the $\Lambda^+_\mathrm{c}$ rapidity interval $\vert y \vert$$\lt$ 1 with transverse momentum ($p_\mathrm{T}$) ranges of 3-30 and 6-40 GeV/$c$ for pp and PbPb collisions, respectively. Compared to the yields in pp collisions scaled by the expected number of nucleon-nucleon interactions, the observed yields of $\Lambda^+_\mathrm{c}$ with $p_\mathrm{T}$$\gt$ 10 GeV/$c$ are strongly suppressed in PbPb collisions. The level of suppression depends significantly on the collision centrality. The $\Lambda^+_\mathrm{c}$ / D$^0$ production ratio is similar in PbPb and pp collisions at $p_\mathrm{T}$$\gt$ 10 GeV/$c$, suggesting that the coalescence process does not play a dominant role in prompt $\Lambda^+_\mathrm{c}$ baryon production at higher $p_\mathrm{T}$.

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.

A search for resonances in events with the $\gamma$+jet final state has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC. The total data analyzed correspond to an integrated luminosity of 138 fb$^{-1}$. Models of excited quarks and quantum black holes are considered. Using a wide-jet reconstruction for the candidate jet, the $\gamma$+jet invariant mass spectrum measured in data is examined for the presence of resonances over the standard model continuum background. The background is estimated by fitting the mass distribution with a functional form. The data exhibit no statistically significant deviations from the expected standard model background. Exclusion limits at 95% confidence level on the resonance mass and other parameters are set. Excited light-flavor quarks (excited bottom quarks) are excluded up to a mass of 6.0 (3.8) TeV. Quantum black hole production is excluded for masses up to 7.5 (5.2) TeV in the Arkani-Hamed-Dimopoulos-Dvali (Randall-Sundrum) model. These lower mass bounds are the most stringent to date among those obtained in the $\gamma$+jet final state.

Ratios of cross sections, $\sigma$(Z+c jets)/$\sigma$(Z+jets), $\sigma$(Z+b jets)/$\sigma$(Z+jets), and $\sigma$(Z+c jets)/$\sigma$(Z+b jets) in the associated production of a Z boson with at least one charm or bottom quark jet are measured in proton-proton collisions at $\sqrt{s}=$ 13 TeV. The data sample, collected by the CMS experiment at the CERN LHC, corresponds to an integrated luminosity of 35.9 fb$^{-1}$, with a fiducial volume of $p_\mathrm{T}>$ 30 GeV and $|\eta|<$ 2.4 for the jets, where $p_\mathrm{T}$ and $\eta$ represent transverse momentum and pseudorapidity, respectively. The Z boson candidates come from leptonic decays into electrons or muons with $p_\mathrm{T}>$ 25 GeV and $|\eta|<$ 2.4, and the dilepton mass satisfies 71 $\lt m_\mathrm{Z}\lt$ 111 GeV. The measured values are $\sigma$(Z+c jets)/$\sigma$(Z+jets) = 0.102 $\pm$ 0.002 $\pm$ 0.009, $\sigma$(Z+b jets)/$\sigma$(Z+jets) = 0.0633 $\pm$ 0.0004 $\pm$ 0.0015, and $\sigma$(Z+c jets)/$\sigma$(Z+b jets) = 1.62 $\pm$ 0.03 $\pm$ 0.15. Results on the inclusive and differential cross section ratios as functions of jet and Z boson transverse momentum are compared with predictions from leading and next-to-leading order perturbative quantum chromodynamics calculations. These are the first measurements of the cross section ratios at 13 TeV, and the measurement precision is better than that of the current theoretical predictions.