Experimental particle physics seeks to understand the universe by probing its fundamental particles and forces and exploring how they govern the large-scale processes that shape cosmic evolution. This whitepaper presents a vision for how Artificial Intelligence (AI) can accelerate discovery in this field. We outline grand challenges that must be addressed to enable transformative breakthroughs and describe how current and planned experimental facilities can implement this vision to advance our understanding of the vast and complex physical world from the smallest to the largest scales. We show how facilities currently under construction, such as the HL-LHC, DUNE and soon EIC, can both benefit from and serve as proving grounds for this vision, while also enabling a longer-term goal for how future experiments – like FCC-ee at CERN, IceCube-Gen2, a Muon Collider in the U.S., and smaller to mid-scale projects – can be fully AI-native. We describe how a truly national-scale collaboration, jointly managed across large funding partners, and involving both DOE laboratories and universities, can make this happen.

The limits on lepton radii and contact interaction scales are used to derive limits on the deviations of lepton anomalous magnetic dipole moments (g-2)/2 from the Standard Model predictions. In the case of quadratic dependence of the deviations on the lepton mass the limits for electrons are better compared to low energy measurements, for muons somewhat weaker than the latest precision experiments, and for taus by far the best.

Modern distributed computing systems produce large amounts of monitoring data. For these systems to operate smoothly, underperforming or failing components must be identified quickly, and preferably automatically, enabling the system managers to react accordingly. In this contribution, we analyze jobs and transfer data collected in the running of the LHC computing infrastructure. The monitoring data is harvested from the Elasticsearch database and converted to formats suitable for further processing. Based on various machine and deep learning techniques, we develop automatic tools for continuous monitoring of the health of the underlying systems. Our initial implementation is based on publicly available deep learning tools, PyTorch or TensorFlow packages, running on state-of-the-art GPU systems.

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

The B$^0_s$ and B$^+$ production yields are measured in PbPb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data sample, collected with the CMS detector at the LHC, corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The mesons are reconstructed in the exclusive decay channels B$^0_s$ $\to$ J/$\psi(\mu^+\mu^-)\phi($K$^+$K$^-)$ and B$^+$ $\to$ J/$\psi(\mu^+\mu^-)$K$^+$, in the transverse momentum range 7-50 GeV/c and absolute rapidity 0-2.4. The B$^0_s$ meson is observed with a statistical significance in excess of five standard deviations for the first time in nucleus-nucleus collisions. The measurements are performed as functions of the transverse momentum of the B mesons and of the PbPb collision centrality. The ratio of production yields of B$^0_s$ and B$^+$ is measured and compared to theoretical models that include quark recombination effects.

A combination of searches for top squark pair production using proton-proton collision data at a center-of-mass energy of 13 TeV at the CERN LHC, corresponding to an integrated luminosity of 137 fb(-1) collected by the CMS experiment, is presented. Signatures with at least 2 jets and large missing transverse momentum are categorized into events with 0, 1, or 2 leptons. New results for regions of parameter space where the kinematical properties of top squark pair production and top quark pair production are very similar are presented. Depending on themodel, the combined result excludes a top squarkmass up to 1325 GeV for amassless neutralino, and a neutralinomass up to 700 GeV for a top squarkmass of 1150 GeV. Top squarks with masses from 145 to 295 GeV, for neutralino masses from 0 to 100 GeV, with a mass difference between the top squark and the neutralino in a window of 30 GeV around the mass of the top quark, are excluded for the first time with CMS data. The results of theses searches are also interpreted in an alternative signal model of dark matter production via a spin-0 mediator in association with a top quark pair. Upper limits are set on the cross section for mediator particle masses of up to 420 GeV.

The second-order Fourier coefficients (v(2)) characterizing the azimuthal distributions of Y(1S) and Y(2S) mesons produced in PbPb collisions at root s(NN) = 5.02 TeV are studied. The Y mesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The collected data set corresponds to an integrated luminosity of 1.7 nb(-1). The scalar product method is used to extract the v2 coefficients of the azimuthal distributions. Results are reported for the rapidity range vertical bar y vertical bar < 2.4, in the transverse momentum interval 0 < pT < 50 GeV/c, and in three centrality ranges of 10-30%, 30-50% and 50-90%. In contrast to the J/psi mesons, the measured v(2) values for the Y mesons are found to be consistent with zero. (C) 2021 The Author(s). Published by Elsevier B.V.

The first measurement of the dependence of $\\gamma\\gamma \\to \\mu^{+}\\mu^{-}$ production on the multiplicity of neutrons emitted very close to the beam direction in ultraperipheral heavy ion collisions is reported. Data for lead-lead interactions at $\\sqrt{s_\\mathrm{NN}} =$ 5.02 TeV, with an integrated luminosity of approximately 1.5 nb$^{-1}$, were collected using the CMS detector at the LHC. The azimuthal correlations between the two muons in the invariant mass region 8 $\\lt m_{\\mu\\mu} \\lt$ 60 GeV are extracted for events including zero, one, or at least two neutrons detected in the forward pseudorapidity range $|\\eta| \\lt$ 8.3. The back-to-back correlation structure from leading-order photon-photon scattering is found to be significantly broader for events with a larger number of emitted neutrons from each nucleus, corresponding to interactions with a smaller impact parameter. This observation provides a data-driven demonstration that the average transverse momentum of photons emitted from relativistic heavy ions has an impact parameter dependence. These results provide new constraints on models of photon-induced interactions in ultraperipheral collisions. They also provide a baseline to search for possible final-state effects on lepton pairs caused by traversing a quark-gluon plasma produced in hadronic heavy ion collisions.

A measurement is presented of electroweak (EW) production of a W boson in association with two jets in proton-proton collisions at s = 13 Te . The data sample was recorded by the CMS Collaboration at the LHC and corresponds to an integrated luminosity of 35.9 fb - 1 . The measurement is performed for the l ν jj final state (with l ν indicating a lepton-neutrino pair, and j representing the quarks produced in the hard interaction) in a kinematic region defined by invariant mass m jj \u003e 120 Ge and transverse momenta p T j \u003e 25 Ge . The cross section of the process is measured in the electron and muon channels yielding σ EW ( W jj ) = 6.23 ± 0.12 (stat) ± 0.61 (syst) pb per channel, in agreement with leading-order standard model predictions. The additional hadronic activity of events in a signal-enriched region is studied, and the measurements are compared with predictions. The final state is also used to perform a search for anomalous trilinear gauge couplings. Limits on anomalous trilinear gauge couplings associated with dimension-six operators are given in the framework of an effective field theory. The corresponding 95% confidence level intervals are - 2.3 c W W W / Λ 2 2.5 Te - 2 , - 8.8 c W / Λ 2 16 Te - 2 , and - 45 c B / Λ 2 46 Te - 2 . These results are combined with the CMS EW Zjj analysis, yielding the constraint on the c W W W coupling: - 1.8 c W W W / Λ 2 2.0 Te - 2 .

A search is conducted for a low-mass charged Higgs boson produced in a top quark decay and subsequently decaying into a charm and a strange quark. The data sample was recorded in proton-proton collisions at $\\sqrt{s}=$ 13 TeV by the CMS experiment at the LHC and corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The search is performed in the process of top quark pair production, where one top quark decays to a bottom quark and a charged Higgs boson, and the other to a bottom quark and a W boson. With the W boson decaying to a charged lepton (electron or muon) and a neutrino, the final state comprises an isolated lepton, missing transverse momentum, and at least four jets, of which two are tagged as b jets. To enhance the search sensitivity, one of the jets originating from the charged Higgs boson is required to satisfy a charm tagging selection. No significant excess beyond standard model predictions is found in the dijet invariant mass distribution. An upper limit in the range 1.68-0.25% is set on the branching fraction of the top quark decay to the charged Higgs boson and bottom quark for a charged Higgs boson mass between 80 and 160 GeV.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in collisions at nucleon-nucleon center-of-mass energies and 5.02 with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

Events with no charged particles produced between the two leading jets are studied in proton-proton collisions at

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\\sqrt{S}$$_{NN}$=2.76 and 5.02TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

Abstract Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in $$\mathrm {PbPb}$$ PbPb collisions at nucleon-nucleon center-of-mass energies $$\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 2.76$$ sNN=2.76 and 5.02$$\,\text {TeV}$$ TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

Central exclusive and semiexclusive production of [see pdf] pairs is measured with the CMS detector in proton-proton collisions at the LHC at center-of-mass energies of 5.02 and 13TeV. 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 [see pdf] production are measured as functions of invariant mass, transverse momentum, and rapidity of the [see pdf] system in the fiducial region defined as transverse momentum [see pdf] and pseudorapidity [see pdf] . The production cross sections for the four resonant channels , [see pdf] , , and are extracted using a simple model. These results represent the first measurement of this process at the LHC collision energies of 5.02 and 13TeV.

Two related searches for phenomena beyond the standard model (BSM) are performed using events with hadronic jets and significant transverse momentum imbalance. The results are based on a sample of proton-proton collisions at a center-of-mass energy of 13 Te , collected by the CMS experiment at the LHC in 2016-2018 and corresponding to an integrated luminosity of 137 fb - 1 . The first search is inclusive, based on signal regions defined by the hadronic energy in the event, the jet multiplicity, the number of jets identified as originating from bottom quarks, and the value of the kinematic variable M T 2 for events with at least two jets. For events with exactly one jet, the transverse momentum of the jet is used instead. The second search looks in addition for disappearing tracks produced by BSM long-lived charged particles that decay within the volume of the tracking detector. No excess event yield is observed above the predicted standard model background. This is used to constrain a range of BSM models that predict the following: the pair production of gluinos and squarks in the context of supersymmetry models conserving R-parity, with or without intermediate long-lived charginos produced in the decay chain; the resonant production of a colored scalar state decaying to a massive Dirac fermion and a quark; or the pair production of scalar and vector leptoquarks each decaying to a neutrino and a top, bottom, or light-flavor quark. In most of the cases, the results obtained are the most stringent constraints to date.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies s NN = 2.76 and 5.02 TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.