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

A search for the rare $\eta$$\to$$\mu^+\mu^-\mu^+\mu^-$ double-Dalitz decay is performed using a sample of proton-proton collisions, collected by the CMS experiment at the CERN LHC with high-rate muon triggers in 2017-2018 and corresponding to an integrated luminosity of 101 fb$^{-1}$. A signal having a statistical significance well in excess of 5 standard deviations is observed. Using the $\eta$$\to$$\mu^+ \mu^-$ decay as normalization, the branching fraction $\mathcal{B}(\eta$$\to$$\mu^+\mu^-\mu^+\mu^-)$ = [5.0 $\pm$ 0.8 (stat) $\pm$ 0.7 (syst) $\pm$ 0.7 ($\mathcal{B}_{2\mu}$)] $\times$ 10$^{-9}$ is measured, where the last term is the uncertainty in the normalization channel branching fraction. This work achieves an improved precision of over five orders of magnitude compared to previous results, leading to the first measurement of this branching fraction, which is found to agree with theoretical predictions.

A search is conducted for a new scalar boson $S$, with a mass distinct from that of the Higgs boson, decaying into four leptons ($\ell =$$e$, $\mu$) via an intermediate state containing two on-shell, promptly decaying new spin-1 bosons $Z_\text{d}$: $S \rightarrow Z_\text{d}Z_\text{d} \rightarrow 4\ell$, where the $Z_\text{d}$ boson has a mass between 15 and 300 GeV, and the $S$ boson has a mass between either 30 and 115 GeV or 130 and 800 GeV. The search uses proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider with an integrated luminosity of 139 fb$^{-1}$ at a centre-of-mass energy of $\sqrt{s}=13$ TeV. No significant excess above the Standard Model background expectation is observed. Upper limits at 95% confidence level are set on the production cross-section times branching ratio, $\sigma(gg \to S) \times \mathcal{B}(S\rightarrow Z_\text{d}Z_\text{d} \rightarrow 4\ell)$, as a function of the mass of both particles, $m_S$ and $m_{Z\text{d}}$.

A search for the resonant production of a heavy scalar $X$ decaying into a Higgs boson and a new lighter scalar $S$, through the process $X \to S(\to bb) H(\to \gamma\gamma)$, where the two photons are consistent with the Higgs boson decay, is performed. The search is conducted using an integrated luminosity of 140 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider. The search is performed over the mass range 170 $\leq$$m_{X}$$\leq$ 1000 GeV and 15 $\leq$$m_{S}$$\leq$ 500 GeV. Parameterised neural networks are used to enhance the signal purity and to achieve continuous sensitivity in a domain of the ($m_{X}$, $m_{S}$) plane. No significant excess above the expected background is found and 95% CL upper limits are set on the cross section times branching ratio, ranging from 39 fb to 0.09 fb. The largest deviation from the background-only expectation occurs for ($m_{X}$, $m_{S}$) = (575, 200) GeV with a local (global) significance of 3.5 (2.0) standard deviations.

A search for Higgs boson pair production via vector-boson fusion is performed in the Lorentz-boosted regime, where a Higgs boson candidate is reconstructed as a single large-radius jet, using 140 fb$^{-1}$ of proton-proton collision data at $\sqrt{s} = 13$ TeV recorded by the ATLAS detector at the Large Hadron Collider. Only Higgs boson decays into bottom quark pairs are considered. The search is particularly sensitive to the quartic coupling between two vector bosons and two Higgs bosons relative to its Standard Model prediction, $\kappa_{2V}$. This study constrains $\kappa_{2V}$ to $0.55 < \kappa_{2V} < 1.49$ at 95% confidence level. The value $\kappa_{2V} = 0$ is excluded with a significance of 3.8 standard deviations with other Higgs boson couplings fixed to their Standard Model values. A search for new heavy spin-0 resonances that would mediate Higgs boson pair production via vector-boson fusion is carried out in the mass range of 1-5 TeV for the first time under several model and decay-width assumptions. No significant deviation from the Standard Model hypothesis is observed and exclusion limits at 95% confidence level are derived.

This paper presents a search for a new $Z^\prime$ resonance decaying into a pair of dark quarks which hadronise into dark hadrons before promptly decaying back as Standard Model particles. This analysis is based on proton-proton collision data recorded at $\sqrt{s}=13$ TeV with the ATLAS detector at the Large Hadron Collider between 2015 and 2018, corresponding to an integrated luminosity of 139 fb$^{-1}$. After selecting events containing large-radius jets with high track multiplicity, the invariant mass distribution of the two highest-transverse-momentum jets is scanned to look for an excess above a data-driven estimate of the Standard Model multijet background. No significant excess of events is observed and the results are thus used to set 95 % confidence-level upper limits on the production cross-section times branching ratio of the $Z^\prime$ to dark quarks as a function of the $Z^\prime$ mass for various dark-quark scenarios.

A search is reported for near-threshold structures in the J/$\psi$J/$\psi$ invariant mass spectrum produced in proton-proton collisions at $\sqrt{s}$ = 13 TeV from data collected by the CMS experiment, corresponding to an integrated luminosity of 135 fb$^{-1}$. Three structures are found, and a model with quantum interference among these structures provides a good description of the data. A new structure is observed with a significance above 5 standard deviations at a mass of 6638 $^{+43}_{-38}$ (stat) $^{+16}_{-31}$ (syst) MeV. Another structure with even higher significance is found at a mass of 6847 $^{+44}_{-28}$ (stat) $^{+48}_{-20}$ (syst) MeV, which is consistent with the X(6900) resonance reported by the LHCb experiment and confirmed by the ATLAS experiment. Evidence for another new structure, with a local significance of 4.7 standard deviations, is found at a mass of 7134 $^{+48}_{-25}$ (stat) $^{+41}_{-15}$ (syst) MeV. Results are also reported for a model without interference, which does not fit the data as well and shows mass shifts up to 150 MeV relative to the model with interference.