MIT Center for Theoretical Physics — a Leinweber Institute Publications

We propose a symplectic structure for the phase space of a generic Lagrangian field theory expressed in the framework of $L_\infty$ algebras. The symplectic structure does not require explicit knowledge of the derivative content of the Lagrangian, and therefore is applicable to nonlocal models, such as string field theory, where traditional constructions are difficult to apply. We test our proposal in a number of examples ranging from general relativity to $p$-adic string theory.

We perform the first high-throughput search and evaluation of materials that can serve as excellent low-mass dark matter detectors. Using properties of close to one thousand materials from the Materials Project database, we project the sensitivity in dark matter parameter space for experiments constructed from each material, including both absorption and scattering processes between dark matter and electrons. Using the anisotropic materials in the dataset, we further compute the level of daily modulation in interaction rates and the resulting directional sensitivities, highlighting materials with prospects to detect the dark matter wind. Our methods provide the basic tools for the data-driven design of dark matter detectors, and our findings lay the groundwork for the next generation of highly optimized direct searches for dark matter as light as the keV scale.

We propose the use of transition-edge sensor (TES) single-photon detectors as a simultaneous target and sensor for direct dark matter searches, and report results from the first search of this kind. We perform a 489 h science run with a TES device optimized for the detection of 1064 nm photons, with a mass of ~0.2 ng and an energy threshold of ~0.3 eV, and set new limits on dark matter interactions with both electrons and nucleons for dark matter with mass below the MeV scale. With their excellent energy resolution, TESs enable search strategies that are complementary to recent results from superconducting nanowire single-photon detectors and kinetic inductance detectors. We show that next-generation TES arrays hold promise to probe new regions of light dark matter parameter space.

Non-relativistic QCD (NRQCD) is used to calculate the relativistic correction to the amplitude for exclusive photoproduction of vector Quarkonia in the near-threshold region within the generalized parton distribution (GPD) framework. The relativistic corrections are found to be large for $J/ψ$. Cross-sections for both $J/ψ$ and $Υ$ are calculated with the former being compared to the data. We also demonstrate the presence of endpoint divergences for the relativistic correction away from the near-threshold regime.

We consider the gauge algebra of closed string field theory with a focus on diffeomorphisms. This algebra contains off-shell information in two ways. The first way is geometric, through the choice of three-punctured sphere defining the three-string vertex. We establish that to leading order in derivatives the superstring algebra is universal: identical for any choice of vertex. For bosonic strings, however, some off-shell dependence remains for vertices that require symmetrization. Off-shell information also appears because field-dependent redefinition of the gauge parameters can alter the algebra. We analyze this dependence in the language of $L_\infty$ algebras, looking at the role of trivial gauge transformations.

We present a construction of the double copy for gauge theories that exhibit off-shell color-kinematics duality within the Batalin-Vilkovitsky (BV) formalism. As illustrative examples, we consider the double copies of Chern-Simons theory, four-dimensional BF theory, and two-dimensional Yang-Mills theory, and we discuss possible gravity interpretations for these cases. We formalize the construction and demonstrate that Kodaira-Spencer gravity, Kähler gravity, and their generalizations, fit naturally within this framework. In particular, Kodaira-Spencer gravity emerges as a gauge theory describing deformations of generalized complex structures, while the double copy of Chern-Simons theory captures the deformation of the Courant bracket.

We derive high-precision results for the $e^+e^-$ heavy jet mass (HJM) $d \sigma/d \rho$ and dihemisphere mass (DHM) $d^2\sigma/(d s_1 d s_2)$ distributions, for $s_1\sim s_2$, in the dijet region. New results include: i) the N$^3$LL resummation for HJM of large logarithms $\ln^n(\rho)$ at small $\rho$ including the exact two-loop non-global hemisphere soft function, the 4-loop cusp anomalous dimension and the 3-loop hard and jet functions, ii) N$^3$LL results for DHM with resummation of logarithms $\ln(s_{1,2}/Q^2)$ when there is no large separation between $s_1$ and $s_2$, iii) profile functions for HJM to give results simultaneously valid in the peak and tail regions, iv) a complete two-dimensional basis of non-perturbative functions which can be used for double differential observables, that are needed for both HJM and DHM in the peak region, and v) an implementation of renormalon subtractions for large-angle soft radiation to ${\cal O}(\alpha_s^3)$ together with a resummation of the additional large $\ln(Q\rho/\Lambda_{QCD})$ logarithms. Here $Q$ is the $e^+e^-$ center-of-mass energy. Our resummation results are combined with known fixed-order ${\cal O}(\alpha_s^3)$ results and we discuss the convergence and remaining perturbative uncertainty in the cross section. We also prove that, at order $1/Q$, the first moment of the HJM distribution involves an additional non-perturbative parameter compared to the power correction that shifts the tail of the spectrum (where $1\gg \rho\gg \Lambda_{QCD}/Q$). This differs from thrust where a single non-perturbative parameter at order $1/Q$ describes both the first moment and the tail, and it disfavors models of power corrections employing a single non-perturbative parameter, such as the low-scale effective coupling model. In this paper we focus only on the dijet region, not the far-tail distribution for $\rho \gtrsim 0.2$.

Dark matter consisting of ultralight bosons can form a macroscopic Bose-Einstein condensate with distinctive observational signatures. While this possibility has been extensively studied for axions and axion-like particles $-$ pseudoscalars with masses protected by shift symmetry $-$ realistic models from string theory and other higher-dimensional theories predict more complex structures. Here we investigate a two-field generalization where an axion couples to a moduli field through its kinetic term, representing the phase and radial modes of a complex scalar field. We demonstrate that when this system forms a gravitationally bound Bose-Einstein condensate, the kinetic coupling produces dramatic modifications to cosmological evolution compared to the canonical single-field case. Most notably, the axion Jeans scale becomes dynamically dependent on the moduli field's evolution, fundamentally altering structure formation. By mapping existing observational constraints from canonical axion models to our two-field scenario, we identify regions of parameter space that are already excluded by current observations. In particular, consistency with observations requires that the moduli field must take on small field values, $\chi/M_{\rm pl} \ll 1$, throughout most of cosmic history for this class of axions to remain a viable description of all dark matter.

We argue that collider observables such as hadron number flux can be matched onto a linear combination of detectors/light-ray operators in perturbative QCD. The spectrum of detectors in QCD is subtle, due to recombination between the DGLAP and BFKL trajectories. We explain how to define and renormalize these trajectories at one-loop, systematically incorporating their recombination. The leading and subleading soft gluon theorems play an important role, and our analysis suggests the presence of an infinite series of further subleading soft theorems for squared-amplitudes/form factors. Combined with our light-ray matching hypothesis, the anomalous dimensions of recombined DGLAP/BFKL detectors yield a prediction for the energy dependence of the number of particles in a jet, as well as other predictions for more general energy-weighted hadron measurements. We compare these predictions to Monte-Carlo simulations, finding good agreement.

We study the theory of massless fields of type II strings arising from the string field theory that uses two string fields, a physical one and an extra one that allows the writing of an action, but whose degrees of freedom ultimately decouple. The mechanism allowing the description of the self-dual five-form of type IIB, anticipated by Sen, is used by the SFT to describe all Ramond-Ramond forms in type IIB and IIA in a manifestly duality-invariant way. We find explicit expressions for the leading terms in the gauge transformation of the RR fields and focus on diffeomorphisms, which are exotic for both the physical and the extra fields, perhaps as needed to describe propagating degrees of freedom that do not gravitate. The algebra of diffeomorphisms includes field-dependent structure constants and only closes on-shell, as predicted by the type II SFT gauge algebra.

We introduce wifi ensembles as a novel framework to obtain asymptotic frequentist uncertainties on density ratios, with a particular focus on neural ratio estimation in the context of high-energy physics. When the density ratio of interest is a likelihood ratio conditioned on parameters, wifi ensembles can be used to perform simulation-based inference on those parameters. After training the basis functions f_i(x), uncertainties on the weights w_i can be straightforwardly propagated to the estimated parameters without requiring extraneous bootstraps. To demonstrate this approach, we present an application in quantum chromodynamics at the Large Hadron Collider, using wifi ensembles to estimate the likelihood ratio between generated quark and gluon jets. We use this learned likelihood ratio to estimate the quark fraction in a synthetic mixed quark/gluon sample, showing that the resultant uncertainties empirically satisfy the desired coverage properties.

We introduce a new mechanism for the simultaneous generation of baryon and dark matter asymmetries through ultraviolet-dominated freeze-in scatterings. The mechanism relies on heavy Majorana neutrinos that connect the visible Standard Model sector to a dark sector through the neutrino portal. Following reheating of the visible sector to a temperature well below the heavy neutrino masses, we show that 2-to-2 scattering processes can populate the dark sector and generate both baryon and dark matter asymmetries. In some parameter regions, the dominant source of baryon asymmetry can be charge transfer from the dark sector, a process we call dark wash-in. We also demonstrate that annihilation of the dark matter to massless states within the dark sector can deplete the symmetric population without destroying the net baryon charge to leave only an asymmetric dark matter abundance today. Depending on the specific model parameters, the observed baryon and dark matter abundances can be attained with heavy neutrino masses $M_N \gtrsim 10^{10}$ GeV, and dark matter masses in the range 0.1 GeV $\lesssim m_\chi \lesssim 10^3$ GeV if the dark matter relic abundance is mainly asymmetric and even lower masses if it is symmetric.

We revisit perturbative unitarity in scalar field inflation with a nonminimal coupling, with Higgs inflation serving as the most prominent example. Although such models are phenomenologically successful, it is critical to examine whether or not unitarity violations spoil their theoretical self-consistency. The analysis of these issues has so far typically relied on order-of-magnitude estimates of scattering amplitudes, which are appropriate for generic parameters. It is not evident that these methods apply to scenarios relying on a near-critical inflationary potential, for which an interplay of both small scalar self-couplings and nonminimal couplings could partially alleviate the unitarity issues. To allow for an exploration of this possibility, we consider the full $S$-matrix for the relevant scattering processes, taking into account important phase space volume factors, leading to a precise evaluation of the cut-off scale. In the single-field case, we demonstrate that near-criticality raises the cut-off scale considerably, compared to previous estimates. In the multifield case, momentum-dependent self-interactions in the kinetic sector lower the cut-off compared to the single-field case to a value comparable to but slightly larger than previous estimates. We carefully study both the single-field and multifield cases in metric and metric-affine (Palatini) formulations of gravity, as well as introduce a new phenomenologically viable model with a canonical kinetic term and a significantly raised cut-off, and discuss the importance of background field effects.

Simulating generic quantum states and dynamics is practically intractable using classical computers. However, certain special classes -- namely Clifford and matchgate circuits -- permit efficient computation. They provide invaluable tools for studying many-body physics, quantum chemistry, and quantum computation. While both play foundational roles across multiple disciplines, the origins of their tractability seem disparate, and their relationship remain unclear. A deeper understanding of such tractable classes could expand their scope and enable a wide range of new applications. In this work, we make progress toward the unified understanding of the Clifford and matchgate -- these two classes are, in fact, distinct special cases of a single underlying structure. Specifically, we introduce the 2D Quon language, which combines Majorana worldlines with their underlying spacetime topology to diagrammatically represent quantum processes and tensor networks. In full generality, the 2D Quon language is universal -- capable of representing arbitrary quantum states, dynamics, or tensor networks -- yet they become especially powerful in describing Clifford and matchgate classes. Each class can be efficiently characterized in a visually recognizable manner using the Quon framework. This capability naturally gives rise to several families of efficiently computable tensor networks introduced in this work: punctured matchgates, hybrid Clifford-matchgate-MPS, and ansatze generated from factories of tractable networks. All of these exhibit high non-Cliffordness, high non-matchgateness, and large bipartite entanglement entropy. We discuss a range of applications of our approach, from recovering well-known results such as the Kramers-Wannier duality and the star-triangle relation of the Ising model, to enabling variational optimization with novel ansatz states.

We show that the honeycomb tight-binding model hosts an exact microscopic avatar of its low-energy SU(2) valley symmetry and parity anomaly. Specifically, the SU(2) valley symmetry arises from a collection of conserved, integer-quantized charge operators that obey the Onsager algebra. Along with lattice reflection and time-reversal symmetries, this Onsager symmetry has a Lieb-Schultz-Mattis (LSM) anomaly that matches the parity anomaly in the IR. Indeed, we show that any local Hamiltonian commuting with these symmetries cannot have a trivial unique gapped ground state. We study the phase diagram of the simplest symmetric model and survey various deformations, including Haldane's mass term, which preserves only the Onsager symmetry. Our results place the parity anomaly in 2+1D alongside Schwinger's anomaly in 1+1D and Witten's SU(2) anomaly in 3+1D as 't Hooft anomalies that can arise from the Onsager symmetry on the lattice.

Effective Field Theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. In this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and Astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at $z=3$. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find non-perturbative velocity contributions (Fingers of God) of a similar size for LBGs to the familiar case of Luminous Red Galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to $\sim 30\%$. Our results indicate that the effective momentum reach $k_{\rm{Max}}$ at $z=3$ for LAEs (LBGs) will be in the range $0.3-0.6 ~h\rm{Mpc}^{-1}$ ($0.2-0.8~h\rm{Mpc}^{-1}$), suggesting that EFT will perform well for high redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

To date, holographic calculations in strongly coupled plasma have provided separate descriptions for the rates of energy loss either for ultrarelativistic massless quarks and gluons or for infinitely massive quarks, with the latter calculation valid for $\sqrt{\gamma} < M/(\sqrt{\lambda}T)$, where $\gamma$ is the Lorentz boost factor for a heavy quark with velocity $v$ and mass $M$ moving through plasma with 't Hooft coupling $\lambda$ and temperature $T$. These two calculations should apply sequentially in the description of the energy loss of a heavy quark that starts out ultrarelativistic, loses energy, slows down, becomes non-relativistic at later times, and ultimately comes to rest and diffuses in the strongly coupled plasma. We provide an ansatz for uniquely incorporating both regimes to give an approximate but unified description of how a heavy quark that is initially ultrarelativistic loses energy all the way until it comes to rest. We implement this ansatz in the Hybrid Strong/Weak Coupling Model. With this new, consistent, treatment of heavy quark energy loss at strong coupling, we confront our predictions for the suppression and azimuthal anisotropies of D- and B-mesons, as well as B-tagged jets, with available experimental data.

Scalar fields are ubiquitous in theories of high-energy physics. In the context of cosmic inflation, this suggests the existence of spectator fields, which provide a subdominant source of energy density. We show that spectator fields boost the inflationary production of primordial black holes, with single-field ultra-slow roll evolution supplanted by a phase of evolution along the spectator direction, and primordial perturbations amplified by the resulting multifield dynamics. This generic mechanism is largely free from the severe fine-tuning that afflicts single-field inflationary PBH models.

The far-from equilibrium dynamics of the pre-hydrodynamic quark-gluon plasma (QGP) formed in heavy ion collisions can be characterized by distinct stages, during each of which the system loses some memory of its initial condition, until only the hydrodynamic modes remain. However, even though it has been repeatedly observed, finding intuitive physical explanations of how and why attractor behavior occurs has remained a challenge. The Adiabatic Hydrodynamization (AH) framework provides exactly such an explanation, showing that the attractor solution can be thought of as the ground state of an analog to quantum mechanical adiabatic evolution, provided we identify appropriate coordinate rescalings. Using the example of a simplified QCD kinetic theory in the small-angle scattering limit, we show how AH can explain both the early pre-hydrodynamic attractor and the later hydrodynamizing attractor in a longitudinally expanding gluon gas in a unified framework. By doing this, we provide a unified description of, and intuition for, all the stages of what in QCD would be bottom-up thermalization, starting from a pre-hydrodynamic attractor and ending with hydrodynamization.