MIT Center for Theoretical Physics — a Leinweber Institute Publications

We study the scaling of meson-meson scattering amplitudes with the number of colors, $N_{\text{c}}$. We use lattice calculations in a theory with $N_{\text{f}}=4$ degenerate flavors, with $N_{\text{c}}=3-6$ and pion mass $M_{\pi }\approx 560$ MeV. We focus on three different scattering channels, two of which have the same quantum numbers as some tetraquark candidates recently found at LHCb: the $T_{cs0}^0\left(2900\right)$, $T_{c\overline{s}0}^{++}\left(2900\right)$, $T_{c\overline{s}0}^0\left(2900\right)$ and $T_{c\overline{s}1}^0\left(2900\right)$ states. Finite-volume energies are extracted using a large set of operators, containing two-particle operators with the form of two pions or two vector mesons, and local tetraquark operators. The resulting energy spectra is used to constrain the infinite-volume scattering amplitude by means of Lüscher's quantization condition. We consider polynomial parametrizations of the phase shift, as well as one-loop chiral perturbation theory (ChPT) predictions. We find that our lattice results follow the expected $N_{\text{c}}$ scaling and are sensitive to subleading $N_{\text{c}}$ corrections. In addition, we constrain the scaling of different combinations of low-energy constants from matching to large $N_{\text{c}}$ ChPT. The results for the channel corresponding to a $({\pi }^{+}{D}_s^{+}-K^{+}{D}^{+})$ state show evidence of a virtual bound state with energy $E_{\text{virtual}}=1.63\left(10\right)M_{\pi }$ for $N_{\text{c}}=3$, while this pole disappears at $N_{\text{c}}>3$. This may be connected to the exotic states found in experiment.

The currently observing Dark Energy Spectroscopic Instrument (DESI) places sub-percent constraints on the Baryon Acoustic Oscillations (BAO) scaling parameters from the Lyman-$\alpha$ (Ly-$\alpha$) forest. However, no systematic error budget stemming from non-linearities in the 3D clustering of the Ly-$\alpha$ forest is included in the DESI-Ly-$\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly-$\alpha$ forest mocks produced on the $N$-body simulation suite \textsc{AbacusSummit}. Specifically, we measure the Ly-$\alpha$ auto-correlation and the Ly-$\alpha$-quasar cross-correlation functions. We use the DESI Ly-$\alpha$ forest fitting pipeline, \textsc{Vega}, with the publicly available covariance matrix from eBOSS DR16. To mitigate the noise, we adopt a linear control variates (LCV) technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the auto-correlation, we detect a small positive shift in radial direction of $\Delta {\alpha }_{\parallel }=0.35\%$ at the 3$\sigma$ level and virtually no shift in the transverse direction, ${\alpha }_{\perp }$. From the cross-correlation, we see a similar shift to $\Delta {\alpha }_{\parallel }$, albeit with larger error bars, and a small negative shift, $\Delta {\alpha }_{\perp }=\sim$0.25%, at the 2$\sigma$ level. We also make a connection with the Ly-$\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly-$\alpha$ auto and the Ly-$\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly-$\alpha$ forest analyses based on EFT. This work paves the way for the full-shape analysis of DESI and future surveys.

The pseudoscalar meson light-cone distribution amplitudes (LCDAs) are essential non-perturbative inputs for a range of high-energy exclusive processes in quantum chromodynamics. In this proceedings, progress towards a determination of the low Mellin moments of the pion and kaon LCDAs by the HOPE Collaboration is reported.

Effective field theory (EFT)-based full-shape analysis with simulation-based priors (SBPs) is a novel approach to galaxy clustering data analysis, which significantly boosts the constraining power by efficiently incorporating field-level simulation information from small scales. So far, SBPs have been mostly extracted from a large set of mock catalogs generated with the Halo Occupation Distribution (HOD) approach. We show that given a halo mass function model and assuming that the EFT parameters of halos depend only on the peak height, HOD-based priors can be computed analytically for the standard 7-parameter HOD model. We derive the relevant expressions for deterministic EFT parameters from the halo model. The halo model, however, fails to accurately describe stochastic EFT parameters, for which we use a physically motivated phenomenological prescription. We compare our analytic priors with simulations and find an excellent agreement. Our approach provides analytic insights into the physics behind the EFT parameters of galaxies, and allows one to reduce the computation time of SBPs to virtually nothing. As an application, we produce a set of 100,000 EFT parameters by sampling both HOD and cosmological models, and explicitly demonstrate that any cosmology dependence can be completely absorbed by relatively minor shifts of the HOD parameters. Our approach can be generalized to other variants of the HOD and cosmological models beyond $\Lambda$CDM.

The so-called metastability bound on the Higgs mass suggests that the smallness of the Higgs mass may be a byproduct of the metastability of the electroweak vacuum. A significantly strong bound requires new physics capable of lowering the scale where the Higgs quartic coupling turns negative through renormalization group effects, without destabilizing the electroweak vacuum entirely. We analyze in this context the low-scale Majoron model of neutrino masses, which automatically contains two key elements for a viable scenario: heavy fermions to lower the instability scale and a extended scalar sector to stabilize the potential and achieve realistic lifetimes for the electroweak vacuum. We show how the metastability bound can be generalized to theories with multiple scalars and present an efficient way of calculating the tunneling rate in such potentials. We also demonstrate that FCC will probe regions of the parameter space relevant for metastability: large regions of the fermionic sector at FCC-ee and some reach to the scalar sector at FCC-hh.

We provide a general prescription for gauging finite non-invertible symmetries in 1+1d lattice Hamiltonian systems. Our primary example is the Rep(D${}_8$) fusion category generated by the Kennedy-Tasaki transformation, which is the simplest anomaly-free non-invertible symmetry on a spin chain of qubits. We explicitly compute its lattice F-symbols and illustrate our prescription for a particular (non-maximal) gauging of this symmetry. In our gauging procedure, we introduce two qubits around each link, playing the role of "gauge fields" for the non-invertible symmetry, and impose novel Gauss's laws. Similar to the Kramers-Wannier transformation for gauging an ordinary $Z_2$, our gauging can be summarized by a gauging map, which is part of a larger, continuous non-invertible cosine symmetry.

Axion-like particles can couple to Standard Model gluons, electroweak gauge bosons, and massive fermions. A future multi-TeV muon collider provides a favorable environment to probe axion-like particles through multiple production channels, including vector boson fusion via electroweak gauge boson couplings and the top-associated production mediated by direct fermionic couplings. Motivated by the quality issue of the QCD axion, we focus on axion-like particles with masses and decay constants around the TeV scale. We explore how different axion-like particle couplings shape its production and decay modes, revealing a rich and intricate phenomenological landscape.

We study systems of two and three mesons composed of pions and kaons at maximal isospin using four CLS ensembles with $a\approx 0.063$fm, including one with approximately physical quark masses. Using the stochastic Laplacian-Heaviside method, we determine the energy spectrum of these systems including many levels in different momentum frames and irreducible representations. Using the relativistic two- and three-body finite-volume formalism, we constrain the two and three-meson K matrices, including not only the leading $s$ wave, but also $p$ and $d$ waves. By solving the three-body integral equations, we determine, for the first time, the physical-point scattering amplitudes for $3{\pi }^{+}$, $3K^{+}$, ${\pi }^{+}{\pi }^{+}{K}^{+}$ and $K^{+}{K}^{+}{\pi }^{+}$ systems. These are determined for total angular momentum $J^P=0^{-}$, $1^{+}$, and $2^{-}$. We also obtain accurate results for $2{\pi }^{+}$, ${\pi }^{+}{K}^{+}$, and $2K^{+}$ phase shifts. We compare our results to Chiral Perturbation Theory, and to phenomenological fits.

We use lattice QCD calculations of the finite-volume spectra of systems of two and three mesons to determine, for the first time, three-particle scattering amplitudes with physical quark masses. Our results are for combinations of ${\pi }^{+}$ and $K^{+}$, at a lattice spacing $a=0.063$fm, and in the isospin-symmetric limit. We also obtain accurate results for maximal-isospin two-meson amplitudes, with those for ${\pi }^{+}{K}^{+}$ and $2K^{+}$ being the first determinations at the physical point. Dense lattice spectra are obtained using the stochastic Laplacian-Heaviside method, and the analysis leading to scattering amplitudes is done using the relativistic finite-volume formalism. Results are compared to chiral perturbation theory and to phenomenological fits to experimental data, finding good agreement.

We present the first study of anti-isolated Upsilon decays to two muons (${\rm Y}\to {\mu }^{+}{\mu }^{-}$) in proton-proton collisions at the Large Hadron Collider. Using a machine learning (ML)-based anomaly detection strategy, we "rediscover" the ${\rm Y}$ in 13 TeV CMS Open Data from 2016, despite overwhelming anti-isolated backgrounds. We elevate the signal significance to $6.4\sigma$ using these methods, starting from $1.6\sigma$ using the dimuon mass spectrum alone. Moreover, we demonstrate improved sensitivity from using an ML-based estimate of the multi-feature likelihood compared to traditional "cut-and-count" methods. Our work demonstrates that it is possible and practical to find real signals in experimental collider data using ML-based anomaly detection, and we distill a readily-accessible benchmark dataset from the CMS Open Data to facilitate future anomaly detection developments.

The energy created in deuterium-tritium fusion reactors originates from a high-intensity neutron flux interacting with the reactor's inner walls. The neutron flux can also be used to produce a self-sustaining reaction by lining the walls with lithium-rich `breeding blankets', in which a fraction of neutrons interacts with lithium, creating the tritium fuel. The high-intensity neutron flux can also result in the production of dark sector particles, feebly interacting light scalars or pseudoscalars, via nuclear transitions within the breeding blanket. We estimate the potential size of such dark sector flux outside the reactor, taking into account all current constraints, and consider possible detection methods at current and future thermonuclear fusion reactors. As a by-product, we also recast the SNO axion bound for a CP even scalar. We find that year-long searches at current and future reactors can set leading constraints on dark scalar -- and dark pseudoscalar -- nucleon couplings.

A global fit for ${\alpha }_s\left(m_Z\right)$ is performed on available $e^{+}{e}^{-}$ data for the heavy jet mass distribution. The state-of-the-art theory prediction includes $O\left({\alpha }_s^3\right)$ fixed-order results, N${}^{}$LL${}^{}$ dijet resummation, N${}^{}$LL Sudakov shoulder resummation, and a first-principles treatment of power corrections in the dijet region. Theoretical correlations are incorporated through a flat random-scan covariance matrix. The global fit results in ${0.1145}_{-0.0019}^{+0.0021}$, compatible with similar determinations from thrust and $C$-parameter. Dijet resummation is essential for a robust fit, as it engenders insensitivity to the fit-range lower cutoff; without resummation the fit-range sensitivity is overwhelming. In addition, we find evidence for a negative power correction in the trijet region if and only if Sudakov shoulder resummation is included.

We present a new solution to the Higgs hierarchy problem based on a dynamical vacuum selection mechanism in a landscape which scans the Higgs mass. In patches where the Higgs mass parameter takes a natural value, the Higgs potential only admits a minimum with a large and negative energy density. This causes a cosmological crunch, removing such patches from the landscape. Conversely, in patches where the Higgs mass parameter is smaller than a critical value, the Higgs potential admits a metastable minimum with the standard cosmological history. This critical value is determined by the instability scale, where the quartic coupling turns negative due to its running. The ability of this mechanism to explain the observed Higgs mass hinges on new physics at the TeV scale, such as vector-like fermions. We study two simple realizations of this scenario in a heavy neutral lepton model and in the singlet-doublet model, the latter mimicking a Higgsino-bino system. We show that the relevant parts of their parameter spaces can be probed by proposed future colliders, such as the FCC-ee or a muon collider.

Composite Higgs models are a class of Beyond the Standard Model (BSM) models proposed to address the hierarchy and naturalness problems associated with the Standard Model (SM) Higgs. A new QCD-like strongly interacting sector based on $SU\left(2\right)$ with two fundamental flavours can be used to build a composite Higgs model which is not yet ruled out by experiment. The role of the singlet scalar resonance will affect Higgs phenomenology at the LHC. In this project our goal is to understand the properties of the singlet scalar state in the new strongly interacting sector in isolation as a first step to understanding the role of this state in composite Higgs models. We present here the first lattice results for the mass of the $\sigma$ in $SU\left(2\right)$ with two fundamental flavours using exponential clover Wilson fermions.

The maximum entropy principle is foundational for statistical analyses of complex dynamics. This principle has been challenged by the findings of a previous work [Phys. Rev. X 7, 031034 (2017)], where it was argued that a quantum system driven in time by a certain aperiodic sequence without any explicit symmetries, dubbed the Thue-Morse drive, gives rise to emergent nonergodic steady states which are underpinned by effective conserved quantities. Here, we resolve this apparent tension. We rigorously prove that the Thue-Morse drive achieves a very strong notion of quantum ergodicity in the long-time limit: the time evolution of any initial state uniformly visits every corner of its Hilbert space. On the other hand, we find the dynamics also approximates a Floquet drive for arbitrarily long albeit finite periods of time with no characteristic timescale, resulting in a scale-free ergodic dynamics we call critically slow complete Hilbert-space ergodicity (CS-CHSE). Furthermore, numerical studies reveal that CS-CHSE is not specific to the Thue-Morse drive and is in fact exhibited by many other aperiodic drives derived from morphic sequences, i.e., words derived from repeatedly applying substitution rules on basic characters. Our work presents a new class of dynamics in time-dependent quantum systems where full ergodicity is eventually attained, but only after astronomically long times.

We explore interacting dark matter (DM) models that allow DM and baryons to scatter off of each other with a cross section that scales with relative particle velocity. Using the effective field theory of large-scale structure, we perform the first analysis of BOSS full-shape galaxy clustering data for velocity-dependent DM-baryon interactions. We determine that while the addition of BOSS full-shape data visibly modifies the shape of the posterior distribution, it does not significantly alter the 95% confidence level intervals for the interaction cross section obtained from an analysis of the cosmic microwave (CMB) anisotropy from Planck measurements alone. Moreover, in agreement with previous findings, we note that the DM-baryon interacting model presents a good fit to both large-scale structure (LSS) data and CMB data and alleviates the $S_8$ tension between the two data sets. After combining LSS and CMB data with weak lensing data from the Dark Energy Survey, we find a $\gtrsim 2\sigma$ preference for non-zero interactions between DM and baryons in a velocity-independent model. We also explore a scenario where only a fraction of DM undergoes scattering with baryons; we find a similar $\gtrsim 2\sigma$ preference for the presence of interactions. Our results suggest that a suppression of the linear matter power spectrum at small scales may be needed to resolve certain discrepancies between LSS and CMB data that are found in the cold DM (CDM) scenario.

We compute scalar static response coefficients (Love numbers) of non-dilatonic black $p$-brane solutions in higher dimensional supergravity. This calculation revels a fine-tuning behavior similar to that of higher dimensional black holes, which we explain by hidden'' near-zone Love symmetries. In general, these symmetries act on equations for perturbations but they are not background isometries. The Love symmetry of charged $p=0$ branes is described by the usual $SL(2,R)$ algebra. For $p=1$ the Love symmetry has an algebraic structure $SL(2,R)\times SL(2,R)$. The $p=0,1$ Love symmetries reduce to isometries of the near-horizon Schwarzschild-AdS${}_{p+2}$ metric in the near-extremal finite temperature limit. They further reduce to the AdS${}_{p+2}$isometries in the extremal zero-temperature limit. We call this process geometrization. In contrast, for the $p>1$ cases, the Love symmetry is always an $SL(2,R)$, and there is no limit in which it becomes geometric. We interpret geometrization and its absence as a consequence of the local equivalence between the Schwarzschild-AdS${}_{p+2}$ and pure AdS${}_{p+2}$ spaces for $p=0,1$, which does not hold for $p>1$. We also show that the static Love numbers of extremal $p$-branes are always zero regardless of spacetime dimensionality, which contrasts starkly with the non-extremal case. Overall, our results suggest that the Love symmetry is hidden by nature, and it can acquire a geometric meaning only if the background has an AdS${}_2$ or AdS${}_3$ limit.

In this paper, we explore a taxonomy of connectivity for space-like structures. It is inspired by isolating posets of connected pieces of a space and examining its embedding in the ambient space. The taxonomy includes in its scope all standard notions of connectivity in point-set and point-free contexts, such as connectivity in graphs and hypergraphs (as well as k-connectivity in graphs), connectivity and path-connectivity in topology, and connectivity of elements in a frame.

We present results on the large $N_{\text{c}}$ scaling of meson-meson scattering amplitudes. We work in a theory with $N_{\text{f}}=4$ degenerate quark flavors and run lattice simulations with $N_{\text{c}}=3-6$ and pion mass $M_{\pi }\approx 590$ MeV. We focus on three different scattering channels, two of which have the same quantum numbers as some tetraquark candidates recently found at LHCb. Finite-volume energies are extracted using a large set of operators, containing two-particle operators corresponding to two pions or two vector mesons, and local tetraquark operators. Using Lüscher's quantization condition, we constrain the infinite-volume scattering amplitudes and investigate subleading $N_{\text{c}}$ corrections to the large $N_{\text{c}}$ limit. For one of the channels, we find indications of a virtual bound state at $N_{\text{c}}=3$, which may be related to one of the aforementioned exotic states.

Normalizing flows have arisen as a tool to accelerate Monte Carlo sampling for lattice field theories. This work reviews recent progress in applying normalizing flows to 4-dimensional nonabelian gauge theories, focusing on two advancements: an architectural improvement referred to as learned active loops, and the application of correlated ensemble methods to QCD with $N_f=2$ dynamical fermions.