The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

This article reports measurements characterizing the Underlying Event (UE) associated with hard scatterings at midrapidity in pp collisions at $\sqrt{s}=13$ TeV. The hard scatterings are identified by the leading particle, the charged particle with the highest transverse momentum ($p_{\rm T}^{\rm leading}$) in the event. Charged-particle number and summed transverse-momentum densities are measured in different azimuthal regions defined with respect to the leading particle direction: Toward, Transverse, and Away. The Toward and Away regions contain the fragmentation products of the hard scatterings in addition to the UE contribution, whereas particles in the Transverse region are expected to originate predominantly from the UE. The study is performed as a function of $p_{\rm T}^{\rm leading}$ with three different $p_{\rm T}$ thresholds for the associated particles, $p_{\rm T}^{\rm min} >$ 0.15, 0.5, and 1.0 GeV/$c$. The charged-particle density in the Transverse region rises steeply for low values of $p_{\rm T}^{\rm leading}$ and reaches a plateau. The results confirm the trend that the charged-particle density in the Transverse region shows a stronger increase with $\sqrt{s}$ than the inclusive charged-particle density at midrapidity. The UE activity is increased by approximately 20% when going from 7 to 13 TeV. The plateau in the Transverse region ($5 < p_{\rm T}^{\rm leading} < ~ 40$ GeV/$c$ ) is further characterized by the probability distribution of its charged-particle multiplicity normalized to its average value (relative transverse activity, $R_{T}$) and the mean transverse momentum as a function of $R_{T}$. Experimental results are compared to model calculations using PYTHIA 8 and EPOS LHC. The overall agreement between models and data is within 30%. These measurements provide new insights on the interplay between hard scatterings and the associated UE in pp collisions.

This Letter presents the first direct investigation of the p-$\Sigma^{0}$ interaction, using the femtoscopy technique in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV measured by the ALICE detector. The $\Sigma^{0}$ is reconstructed via the decay channel to $\Lambda \gamma$, and the subsequent decay of $\Lambda$ to p$\pi^-$. The photon is detected via the conversion in material to e$^{+}$e$^{-}$ pairs exploiting the unique capability of the ALICE detector to measure electrons at low transverse momenta. The measured p-$\Sigma^{0}$ correlation indicates a shallow strong interaction. The comparison of the data to several theoretical predictions obtained employing the $Correlation~Analysis~Tool~using~the~Schr\"odinger~Equation$ (CATS) and the Lednick\'y-Lyuboshits approach shows that the current experimental precision does not yet allow to discriminate between different models, as it is the case for the available scattering and hypernuclei data. Nevertheless, the p-$\Sigma^{0}$ correlation function is found to be sensitive to the strong interaction, and driven by the interplay of the different spin and isospin channels. This pioneering study demonstrates the feasibility of a femtoscopic measurement in the p-$\Sigma^{0}$ channel and with the expected larger data samples in LHC Run 3 and Run 4, the p-$\Sigma^{0}$ interaction will be constrained with high precision.

Deuteron production in high-energy collisions is sensitive to the space-time evolution of the collision system, and is typically described by a coalescence mechanism. For the first time, we present results on jet-associated deuteron production in \pp collisions at $\sqrt{s}=13$ TeV, providing an opportunity to test the established picture for deuteron production in events with a hard scattering. Using a trigger particle with high transverse-momentum ($p_{\rm T}>5$ GeV/$c$) as a proxy for the presence of a jet at midrapidity, we observe a measurable population of deuterons being produced around the jet proxy. The associated deuteron yield measured in a narrow angular range around the trigger particle differs by 2.4-4.8 standard deviations from the uncorrelated background. The data are described by PYTHIA model calculations featuring baryon coalescence.

The cross section of jets reconstructed from charged particles is measured in the transverse momentum range of $5<p_\mathrm{T}<100\ \mathrm{GeV}/c$ in pp collisions at the center-of-mass energy of $\sqrt{s} = 5.02\ \mathrm{TeV}$ with the ALICE detector. The jets are reconstructed using the anti-$k_\mathrm{T}$ algorithm with resolution parameters $R=0.2$, $0.3$, $0.4$, and $0.6$ in the pseudorapidity range $|\eta|< 0.9-R$. The charged jet cross sections are compared with the leading order (LO) and to next-to-leading order (NLO) perturbative Quantum ChromoDynamics (pQCD) calculations. It was found that the NLO calculations agree better with the measurements. The cross section ratios for different resolution parameters were also measured. These ratios increase from low $p_\mathrm{T}$ to high $p_\mathrm{T}$ and saturate at high $p_\mathrm{T}$, indicating that jet collimation is larger at high $p_\mathrm{T}$ than at low $p_\mathrm{T}$. These results provide a precision test of pQCD predictions and serve as a baseline for the measurement in Pb$-$Pb collisions at the same energy to quantify the effects of the hot and dense medium created in heavy-ion collisions at the LHC.

A measurement of the W$^+$W$^-$ boson pair production cross section in proton-proton collisions at $\sqrt{s} =$ 13 TeV is presented. The data used in this study are collected with the CMS detector at the CERN LHC and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The W$^+$W$^-$ candidate events are selected by requiring two oppositely charged leptons (electrons or muons). Two methods for reducing background contributions are employed. In the first one, a sequence of requirements on kinematic quantities is applied allowing a measurement of the total production cross section: 117.6 $\pm$ 6.8 pb, which agrees well with the theoretical prediction. Fiducial cross sections are also reported for events with zero or one jet, and the change in the zero-jet fiducial cross section with the jet transverse momentum threshold is measured. Normalized differential cross sections are reported within the fiducial region. A second method for suppressing background contributions employs two random forest classifiers. The analysis based on this method includes a measurement of the total production cross section and also a measurement of the normalized jet multiplicity distribution in W$^+$W$^-$ events. Finally, a dilepton invariant mass distribution is used to probe for physics beyond the standard model in the context of an effective field theory, and constraints on the presence of dimension-6 operators are derived.

The first experimental information on the strong interaction between $\Lambda$ and $\Xi^-$ strange baryons is presented in this Letter. The correlation function of $\Lambda-\Xi^-$ and $\overline{\Lambda}-\overline{\Xi}^{+}$ pairs produced in high-multiplicity proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the LHC is measured as a function of the relative momentum of the pair. The femtoscopy method is used to calculate the correlation function, which is then compared with theoretical expectations obtained using a meson exchange model, chiral effective field theory, and Lattice QCD calculations close to the physical point. Data support predictions of small scattering parameters while discarding versions with large ones, thus suggesting a weak $\Lambda-\Xi^{-}$ interaction. The limited statistical significance of the data does not yet allow one to constrain the effects of coupled channels like $\Sigma-\Xi$ and N$-\Omega$.

The study of the strength and behaviour of the antikaon-nucleon ($\mathrm{\overline{K}N}$) interaction constitutes one of the key focuses of the strangeness sector in low-energy Quantum Chromodynamics (QCD). In this letter a unique high-precision measurement of the strong interaction between kaons and protons, close and above the kinematic threshold, is presented. The femtoscopic measurements of the correlation function at low pair-frame relative momentum of (K$^+$ p $\oplus$ K$^-$ $\overline{\mathrm{p}}$) and (K$^-$ p $\oplus$ K$^+$ $\overline{\mathrm{p}}$) pairs measured in pp collisions at $\sqrt{s}$ = 5, 7 and 13 TeV are reported. A structure observed around a relative momentum of 58 MeV/$c$ in the measured correlation function of (K$^-$ p $\oplus$ K$^+$ $\overline{\mathrm{p}}$) with a significance of 4.4. $\sigma$ constitutes the first experimental evidence for the opening of the $(\mathrm{\overline{K}^0 n} \oplus \mathrm{K^0 \overline{n}})$ isospin breaking channel due to the mass difference between charged and neutral kaons. The measured correlation functions have been compared to J\"{u}lich and Kyoto models in addition to the Coulomb potential. The high-precision data at low relative momenta presented in this work prove femtoscopy to be a powerful complementary tool to scattering experiments and provide new constraints above the $\mathrm{\overline{K}N}$ threshold for low-energy QCD chiral models.

A combination of four searches for new physics involving signatures with at least one photon and large missing transverse momentum, motivated by generalized models of gauge-mediated supersymmetry (SUSY) breaking, is presented. All searches make use of proton-proton collision data at $\sqrt{s}=$ 13 TeV, which were recorded with the CMS detector at the LHC in 2016, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. Signatures with at least one photon and large missing transverse momentum are categorized into events with two isolated photons, events with a lepton and a photon, events with additional jets, and events with at least one high-energy photon. No excess of events is observed beyond expectations from standard model processes, and limits are set in the context of gauge-mediated SUSY. Compared to the individual searches, the combination extends the sensitivity to gauge-mediated SUSY in both electroweak and strong production scenarios by up to 100 GeV in neutralino and chargino masses, and yields the first CMS result combining various SUSY searches in events with photons at $\sqrt{s}=$ 13 TeV.

The inclusive production of the charm-strange baryon $\rm \Omega_c^0$ is measured for the first time via its hadronic decay into $\rm \Omega^-\pi^+$ at midrapidity ($|y|<0.5$) in proton-proton (pp) collisions at the centre-of-mass energy $\sqrt{s}=13$ TeV with the ALICE detector at the LHC. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval $2<p_{\rm T}<12~{\rm GeV}/c$. The $p_{\rm T}$ dependence of the $\rm \Omega_c^0$-baryon production relative to the prompt $\rm D^0$-meson and to the prompt $\rm \Xi_c^0$-baryon production is compared to various models that take different hadronisation mechanisms into consideration. In the measured $p_{\rm T}$ interval, the ratio of the $p_{\rm T}$-integrated cross sections of $\rm \Omega_c^0$ and prompt $\Lambda_{\rm c}^{+}$ baryons multiplied by the $\rm \Omega^-\pi^+$ branching ratio is found to be larger by a factor of about 20 with a significance of about 4$\sigma$ when compared to $\rm e^+e^-$ collisions.

Measurements of the inclusive J/$\psi$ yield as a function of charged-particle pseudorapidity density ${\rm d}N_{\rm ch}/\rm{d}\eta$ in pp collisions at $\sqrt{s}$ = 13 TeV with ALICE at the LHC are reported. The J/$\psi$ meson yield is measured at midrapidity ($|y|<0.9$) in the dielectron channel, for events selected based on the charged-particle multiplicity at midrapidity ($|\eta|<1$) and at forward rapidity ($-3.7<\eta<-1.7$ and $2.8<\eta<5.1$); both observables are normalized to their corresponding averages in minimum bias events. The increase of the normalized J/$\psi$ yield with normalized ${\rm d}N_{\rm ch}/\rm{d}\eta$ is significantly stronger than linear and dependent on the transverse momentum. The data are compared to theoretical predictions, which describe the observed trends well, albeit not always quantitatively.

The production rates and the transverse momentum distribution of strange hadrons at mid-rapidity ($\ |y\ | < 0.5$) are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV as a function of the charged particle multiplicity, using the ALICE detector at the LHC. The production rates of $\rm{K}^{0}_{S}$, $\Lambda$, $\Xi$, and $\Omega$ increase with the multiplicity faster than what is reported for inclusive charged particles. The increase is found to be more pronounced for hadrons with a larger strangeness content. Possible auto-correlations between the charged particles and the strange hadrons are evaluated by measuring the event-activity with charged particle multiplicity estimators covering different pseudorapidity regions. When comparing to lower energy results, the yields of strange hadrons are found to depend only on the mid-rapidity charged particle multiplicity. Several features of the data are reproduced qualitatively by general purpose QCD Monte Carlo models that take into account the effect of densely-packed QCD strings in high multiplicity collisions. However, none of the tested models reproduce the data quantitatively. This work corroborates and extends the ALICE findings on strangeness production in proton-proton collisions at 7 TeV.

This paper presents the measurements of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{p}$ and $\bar{\rm{p}}$ transverse momentum ($p_{\rm{T}}$) spectra as a function of charged-particle multiplicity density in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector at the LHC. Such study allows us to isolate the center-of-mass energy dependence of light-flavour particle production. The measurements reported here cover a $p_{\rm{T}}$ range from 0.1 GeV/$c$ to 20 GeV/$c$ and are done in the rapidity interval $|y|<0.5$. The $p_{\rm{T}}$-differential particle ratios exhibit an evolution with multiplicity, similar to that observed in pp collisions at $\sqrt{s}$ = 7 TeV, which is qualitatively described by some of the hydrodynamical and pQCD-inspired models discussed in this paper. Furthermore, the $p_{\rm{T}}$-integrated hadron-to-pion yield ratios measured in pp collisions at two different center-of-mass energies are consistent when compared at similar multiplicities. This also extends to strange and multistrange hadrons, suggesting that, at LHC energies, particle hadrochemistry scales with particle multiplicity the same way under different collision energies and colliding systems.

We present results on transverse momentum ($p_{\rm T}$) and rapidity ($y$) differential production cross sections, mean transverse momentum and mean transverse momentum square of inclusive J/$\psi$ and $\psi(2S)$ at forward rapidity ($2.5<y<4$) as well as $\psi(2S)$-to-J/$\psi$ cross section ratios. These quantities are measured in pp collisions at center of mass energies $\sqrt{s}=5.02$ and 13 TeV with the ALICE detector. Both charmonium states are reconstructed in the dimuon decay channel, using the muon spectrometer. A comprehensive comparison to inclusive charmonium cross sections measured at $\sqrt{s}=2.76$, 7 and 8 TeV is performed. A comparison to non-relativistic quantum chromodynamics and fixed-order next-to-leading logarithm calculations, which describe prompt and non-prompt charmonium production respectively, is also presented. A good description of the data is obtained over the full $p_{\rm T}$ range, provided that both contributions are summed. In particular, it is found that for $p_{\rm T}>15$ GeV/$c$ the non-prompt contribution reaches up to 50% of the total charmonium yield.

Measurements are presented of the differential cross sections for Z bosons produced in proton-proton collisions at $\sqrt{s} =$ 13 TeV and decaying to muons and electrons. The data analyzed were collected in 2016 with the CMS detector at the LHC and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The measured fiducial inclusive product of cross section and branching fraction agrees with next-to-next-to-leading order quantum chromodynamics calculations. Differential cross sections of the transverse momentum $p_\mathrm{T}$, the optimized angular variable $\phi^*_\eta$, and the rapidity of lepton pairs are measured. The data are corrected for detector effects and compared to theoretical predictions using fixed order, resummed, and parton shower calculations. The uncertainties of the measured normalized cross sections are smaller than 0.5% for $\phi^*_\eta <$ 0.5 and for $p_\mathrm{T}^\mathrm{Z} <$ 50 GeV.

The production of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{K}^{0}_{S}$, $\rm{K}^{*}(892)^{0}$, $\rm{p}$, $\phi(1020)$, $\Lambda$, $\Xi^{-}$, $\Omega^{-}$, and their antiparticles was measured in inelastic proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV at midrapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm{T}}$) using the ALICE detector at the CERN LHC. Furthermore, the single-particle $p_{\rm{T}}$ distributions of $\rm{K}^{0}_{S}$, $\Lambda$, and $\overline{\Lambda}$ in inelastic pp collisions at $\sqrt{s}$ = 7 TeV are reported here for the first time. The $p_{\rm{T}}$ distributions are studied at midrapidity within the transverse momentum range $0\leq p_{\rm{T}}\leq20$ GeV/$c$, depending on the particle species. The $p_{\rm{T}}$ spectra, integrated yields, and particle yield ratios are discussed as a function of collision energy and compared with measurements at lower $\sqrt{s}$ and with results from various general-purpose QCD-inspired Monte Carlo models. A hardening of the spectra at high $p_{\rm{T}}$ with increasing collision energy is observed, which is similar for all particle species under study. The transverse mass and $x_{\rm{T}}\equiv2p_{\rm{T}}/\sqrt{s}$ scaling properties of hadron production are also studied. As the collision energy increases from $\sqrt{s}$ = 7 to 13 TeV, the yields of non- and single-strange hadrons normalized to the pion yields remain approximately constant as a function of $\sqrt{s}$, while ratios for multi-strange hadrons indicate enhancements. The $p_{\rm{T}}$-differential cross sections of $\pi^{\pm}$, $\rm{K}^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) are compared with next-to-leading order perturbative QCD calculations, which are found to overestimate the cross sections for $\pi^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) at high $p_{\rm{T}}$.

A search is performed for neutral non-standard-model Higgs bosons decaying to two muons in the context of the minimal supersymmetric standard model (MSSM). Proton-proton collision data recorded by the CMS experiment at the CERN Large Hadron Collider at a center-of-mass energy of 13 TeV were used, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search is sensitive to neutral Higgs bosons produced via the gluon fusion process or in association with a $\mathrm{b\overline{b}}$ quark pair. No significant deviations from the standard model expectation are observed. Upper limits at 95% confidence level are set in the context of the $m_\mathrm{h}^{\text{mod+}}$ and phenomenological MSSM scenarios on the parameter $\tan\beta$ as a function of the mass of the pseudoscalar A boson, in the range from 130 to 600 GeV. The results are also used to set a model-independent limit on the product of the branching fraction for the decay into a muon pair and the cross section for the production of a scalar neutral boson, either via gluon fusion, or in association with b quarks, in the mass range from 130 to 1000 GeV.

This work presents new constraints on the existence and the binding energy of a possible $\Lambda$-$\Lambda$ bound state, the H-dibaryon, derived from $\Lambda$-$\Lambda$ femtoscopic measurements by the ALICE collaboration. The results are obtained from a new measurement using the femtoscopy technique in pp collisions at $\sqrt{s}=13$ TeV and p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV, combined with previously published results from p-Pb collisions at $\sqrt{s}=7$ TeV. The $\Lambda$-$\Lambda$ scattering parameter space, spanned by the inverse scattering length $f_0^{-1}$ and the effective range $d_0$, is constrained by comparing the measured $\Lambda$-$\Lambda$ correlation function with calculations obtained within the Lednicky model. The data are compatible with hypernuclei results and lattice computations, both predicting a shallow attractive interaction, and permit to test different theoretical approaches describing the $\Lambda$-$\Lambda$ interaction. The region in the $(f_0^{-1},d_0)$ plane which would accommodate a $\Lambda$-$\Lambda$ bound state is substantially restricted compared to previous studies. The binding energy of the possible $\Lambda$-$\Lambda$ bound state is estimated within an effective-range expansion approach and is found to be $B_{\Lambda\Lambda}=3.2^{+1.6}_{-2.4}\mathrm{(stat)}^{+1.8}_{-1.0}\mathrm{(syst)}$ MeV.

A search for long-lived particles decaying into jets is presented. Data were collected with the CMS detector at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV in 2016, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The search examines the distinctive topology of displaced tracks and secondary vertices. The selected events are found to be consistent with standard model predictions. For a simplified model in which long-lived neutral particles are pair produced and decay to two jets, pair production cross sections larger than 0.2 fb are excluded at 95% confidence level for a long-lived particle mass larger than 1000 GeV and proper decay lengths between 3 and 130 mm. Several supersymmetry models with gauge-mediated supersymmetry breaking or $R$-parity violation, where pair-produced long-lived gluinos or top squarks decay to several final-state topologies containing displaced jets, are also tested. For these models, in the mass ranges above 200 GeV, gluino masses up to 2300-2400 GeV and top squark masses up to 1350-1600 GeV are excluded for proper decay lengths approximately between 10 and 100 mm. These are the most restrictive limits to date on these models.

Central exclusive and semiexclusive production of $\pi^+\pi^-$ pairs is measured with the CMS detector in proton-proton collisions at the LHC at center-of-mass energies of 5.02 and 13 TeV. The theoretical description of these nonperturbative processes, which have not yet been measured in detail at the LHC, poses a significant challenge to models. The two pions are measured and identified in the CMS silicon tracker based on specific energy loss, whereas the absence of other particles is ensured by calorimeter information. The total and differential cross sections of exclusive and semiexclusive central $\pi^+\pi^-$ production are measured as functions of invariant mass, transverse momentum, and rapidity of the $\pi^+\pi^-$ system in the fiducial region defined as transverse momentum $p_\mathrm{T}(\pi)$ $>$ 0.2 GeV and pseudorapidity $|\eta(\pi)|$ $<$ 2.4. The production cross sections for the four resonant channels f$_0(500)$, $\rho^0(770)$, f$_0(980)$, and f$_2(1270)$ are extracted using a simple model. These results represent the first measurement of this process at the LHC collision energies of 5.02 and 13 TeV.

Measurements of identified hadrons as a function of the charged-particle multiplicity in pp collisions enable a search for the onset of collective effects in small collision systems. With such measurements, it is possible to study the mechanisms that determine the shapes of hadron transverse momentum ($p_{\rm{T}}$) spectra, to search for possible modifications of the yields of short-lived hadronic resonances due to scattering effects in the hadron-gas phase, and to investigate different explanations for the multiplicity evolution of strangeness production provided by phenomenological models. In this paper, these topics are addressed through measurements of the $\rm{K}^{*}(892)^{0}$ and $\phi(1020)$ mesons at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV as a function of the charged-particle multiplicity. The results include the $p_{\rm{T}}$ spectra, $p_{\rm{T}}$-integrated yields, mean transverse momenta, and the ratios of the yields of these resonances to those of longer-lived hadrons. Comparisons with results from other collision systems and energies, as well as predictions from phenomenological models, are also discussed.

Results are reported from a search for supersymmetric particles in the final state with multiple jets and large missing transverse momentum. The search uses a sample of proton-proton collisions at $\sqrt{s} =$ 13 TeV collected with the CMS detector in 2016-2018, corresponding to an integrated luminosity of 137 fb$^{-1}$, representing essentially the full LHC Run 2 data sample. The analysis is performed in a four-dimensional search region defined in terms of the number of jets, the number of tagged bottom quark jets, the scalar sum of jet transverse momenta, and the magnitude of the vector sum of jet transverse momenta. No significant excess in the event yield is observed relative to the expected background contributions from standard model processes. Limits on the pair production of gluinos and squarks are obtained in the framework of simplified models for supersymmetric particle production and decay processes. Assuming the lightest supersymmetric particle to be a neutralino, lower limits on the gluino mass as large as 2000 to 2310 GeV are obtained at 95% confidence level, while lower limits on the squark mass as large as 1190 to 1630 GeV are obtained, depending on the production scenario.

A first observation is presented for the electroweak production of a W boson, a photon, and two jets in proton-proton collisions. The W boson decays are selected by requiring one identified electron or muon and an imbalance in transverse momentum. The two jets are required to have a high dijet mass and a large separation in pseudorapidity. The measurement is based on data collected with the CMS detector at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The observed (expected) significance for this process is 4.9 (4.6) standard deviations. After combining with previously reported CMS results at 8 TeV, the observed (expected) significance is 5.3 (4.8) standard deviations. The cross section for the electroweak W$\gamma$jj production in a restricted fiducial region is measured as 20.4 $\pm$ 4.5 fb and the total cross section for W$\gamma$ production in association with 2 jets in the same fiducial region is 108 $\pm$ 16 fb. All results are in good agreement with recent theoretical predictions. Constraints are placed on anomalous quartic gauge couplings in terms of dimension-8 effective field theory operators.

Multiplicity ($N_{\rm ch}$) distributions and transverse momentum ($p_{\rm T}$) spectra of inclusive primary charged particles in the kinematic range of $|\eta| < 0.8$ and 0.15 GeV/$c$$< p_{T} <$ 10 GeV/$c$ are reported for pp, p-Pb, Xe-Xe and Pb-Pb collisions at centre-of-mass energies per nucleon pair ranging from $\sqrt{s_{\rm NN}} = 2.76$ TeV up to $13$ TeV. A sequential two-dimensional unfolding procedure is used to extract the correlation between the transverse momentum of primary charged particles and the charged-particle multiplicity of the corresponding collision. This correlation sharply characterises important features of the final state of a collision and, therefore, can be used as a stringent test of theoretical models. The multiplicity distributions as well as the mean and standard deviation derived from the $p_{\rm T}$ spectra are compared to state-of-the-art model predictions. Providing these fundamental observables of bulk particle production consistently across a wide range of collision energies and system sizes can serve as an important input for tuning Monte Carlo event generators.

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.