MIT Center for Theoretical Physics — a Leinweber Institute Publications

Gapped two-dimensional gauge theories with massless fermions generically have rich vacuum structures consisting of many degenerate vacua related by the action of topological line operators. The algebra of such operators has been used to calculate ratios of vacuum expectation values of local operators and to predict nontrivial particle-soliton degeneracies. In this paper, we use recently-developed tensor network methods to study several examples of such theories via their Hamiltonian lattice descriptions. Our lattice results agree with all previously-made predictions. Furthermore, we identify the lattice strong-coupling states that can be adiabatically continued to the degenerate vacua in the continuum limit. We conjecture a procedure, referred to as a lattice decay rule, for how this identification works in general. This rule allows us to compute the continuum vacuum degeneracy by studying the lattice Hamiltonian in the strong-coupling limit.

We present cosmological parameters measurements from the full combination of DESI DR1 galaxy clustering data described with large-scale structure effective field theory. By incorporating additional datasets (photometric galaxies and CMB lensing cross-correlations) and extending the bispectrum likelihood to smaller scales using a consistent one-loop theory computation, we achieve substantial gains in constraining power relative to previous analyses. Combining with the latest DESI baryon acoustic oscillation data and using cosmic microwave background (CMB) priors on the power spectrum tilt and baryon density, we obtain tight constraints on the $Λ$CDM model, finding the Hubble constant $H_0=69.08\pm 0.37~\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, the matter density fraction $Ω_m=0.2973\pm 0.0050$, and the mass fluctuation amplitude $σ_8 = 0.815\pm 0.016$ (or the lensing parameter $S_8\equivσ_8\sqrt{Ω_m/0.3}=0.811\pm 0.016$), corresponding to $0.6\%$, $1.7\%$, and $2\%$ precision respectively. Adding the Pantheon+ supernova sample (SNe), we find a preference of $2.6σ$ for the $w_0w_a$ dynamical dark energy model from low-redshift data alone, which increases to $2.8σ$ when exchanging the SNe with Planck CMB data. Combining full-shape data with BAO, CMB, and SNe likelihoods, we improve the dark energy figure-of-merit by $18\%$ and bound the sum of the neutrino masses to $M_ν<0.057$ eV in $Λ$CDM and $M_ν<0.095$ eV in the $w_0w_a$ dynamical dark energy model (both at 95% CL). This represents an improvement of $25\%$ over the background expansion constraints and the strongest bound on neutrino masses in $w_0w_a$CDM to date. Our results suggest that the preference for the normal ordering of neutrino mass states holds regardless of the cosmological background model, and is robust in light of tensions between cosmological datasets.

Primordial black holes (PBHs) are theorized to form from the collapse of overdensities in the very early Universe. PBHs in the asteroid-mass range $10^{17} \, {\rm g}\lesssim M \lesssim 10^{23} \, {\rm g}$ could serve as all or most of the dark matter today, but are particularly difficult to detect due to their modest rates of Hawking emission and sub-micron Schwarzschild radii. We consider whether the steep gradients of a PBH’s gravitational field could generate tidal forces strong enough to disrupt atoms and nuclei. Such phenomena may yield new observables that could uniquely distinguish a PBH from a macroscopic object of the same mass. We first consider the gravitational ionization of ambient neutral hydrogen and evaluate prospects for detecting photon radiation from the recombination of ionized atoms. During the present epoch, this effect would be swamped by Hawking radiation — which would itself be difficult to detect for PBHs at the upper end of the asteroid-mass window. We then consider the gravitational ionization and heating of neutral hydrogen immediately following recombination at $z\simeq1090$, and identify a broad class of PBH distributions with typical mass $5\times10^{21}\,{\rm g}\lesssim M \lesssim 10^{23}\, {\rm g}$ within which gravitational interactions would have been the dominant form of energy deposition to the medium. We also identify conditions under which tidal forces from a transiting PBH could overcome the strong nuclear force, either by dissociating deuterons, which would be relevant during big bang nucleosynthesis (BBN), or by inducing fission of heavy nuclei. We find that gravitational dissociation of deuterons dominates photodissociation rates due to Hawking radiation for PBHs with masses $10^{14}\,{\rm g}\lesssim M \lesssim 10^{16}\,{\rm g}$. We additionally identify the phenomenon of gravitationally induced fission of heavy nuclei via tidal deformation.

We investigate the sensitivity of the Electron-Ion Collider (EIC) to invisible final states in coherent exclusive electroproduction. The characteristic signal is a forward proton with reduced energy and little additional detector activity. Using the excellent particle detection capabilities and kinematics reconstruction at the EIC, we argue that backgrounds can be strongly suppressed. While our analysis applies to various states, we specifically focus on pseudoscalar particles: (i) neutral mesons ($π^0,η^{(\prime)}$), whose invisible Standard Model decays are extremely suppressed, and (ii) gluon-coupled axion-like particles (ALPs) decaying invisibly to a dark sector. Depending on the meson species and the achievable background rejection, the EIC could strengthen existing bounds on invisible decays of pseudoscalar mesons by up to four orders of magnitude, probing branching ratios as small as ${\rm BR}(η^{(\prime)}\to{\rm inv})\sim 10^{-8}$. In addition, the EIC would directly probe invisibly decaying ALPs with the couplings up to $f_a\sim 10^5\,\text{GeV}$ and masses in the range $0.1$-$2\,\text{GeV}$.

We develop a quadratic-in-Riemann worldline action for a Kerr black hole at infinite spin orders by matching to a proposed tree-level Kerr Compton amplitude, originally obtained from higher-spin QFT considerations. A worldline action is an effective theory, and as such the tree-level matching needs to be corrected by loop effects, including UV counter terms, renormalization, and higher-order matching to general relativity. However, we anticipate that many features of the Wilson coefficients of the proposed tree-level action will remain unchanged even after a loop-level matching. While the worldline action is given in closed form, it contains an infinite number of quadratic-in-Riemann operators $R^2$, even for the same-helicity sector. We argue that in the same-helicity sector the $R^2$ operators have no intrinsic meaning, as they merely remove unwanted terms produced by the linear-in-Riemann operators, which are well-established in the literature. The opposite-helicity sector is somewhat more complicated, it contains both $R^2$ operators that removes unwanted terms, and $R^2$ operators that add new needed terms to the Compton amplitude. We discuss and classify all independent $R^2$ operators that can feature in the worldline action.

Gapless quantum phases can become distinct when internal symmetries are enforced, in analogy with gapped symmetry-protected topological (SPT) phases. However, this distinction does not always lead to protected edge modes, raising the question of how the bulk-boundary correspondence is generalized to gapless cases. We propose that the spatial interface between gapless phases — rather than their boundaries — provides a more robust fingerprint. We show that whenever two 1+1d conformal field theories (CFTs) differ in symmetry charge assignments of local operators or twisted sectors, any symmetry-preserving spatial interface between the theories must flow to a non-invertible defect. We illustrate this general result for different versions of the Ising CFT with $\mathbb{Z}_2 \times \mathbb{Z}_2^T$ symmetry, obtaining a complete classification of allowed conformal interfaces. When the Ising CFTs differ by nonlocal operator charges, the interface hosts 0+1d symmetry-breaking phases with finite-size splittings scaling as $1/L^3$, as well as continuous phase transitions between them. For general gapless phases differing by an SPT entangler, the interfaces between them can be mapped to conformal defects with a certain defect ‘t Hooft anomaly. This classification also gives implications for higher-dimensional examples, including symmetry-enriched variants of the 2+1d Ising CFT. Our results establish a physical indicator for symmetry-enriched criticality through symmetry-protected interfaces, giving a new handle on the interplay between topology and gapless phases.

It was recently shown that black holes decohere any quantum superpositions in their vicinity. This decoherence is mediated by soft radiation through the horizon, and can be understood as the result of the fact that quantum states in the exterior source distinguishable states of long-range fields in the interior. To study this phenomenon and others, we extend Tomita-Takesaki theory to accommodate states of soft radiation such as arise in the electromagnetic and gravitational memory effects, and provide a general framework for computing the distinguishability of general coherent states. Applying these tools, we use the methods of unambiguous state discrimination and approximate quantum error correction to prove some new relations regarding the distinguishability of quantum states, and the quantum information content of soft radiation, and thereby show that a black hole (or any horizon) decoheres its environment as though its interior were full of optimal observers.

In this paper, we study pole-skipping in non-maximally quantum chaotic systems. Using Rindler conformal field theories and the large-$q$ SYK chain as illustrative examples, we argue that the pole-skipping points of few-body operators organize into trajectories in the complex frequency-momentum plane, with the leading trajectory encoding the quantum Lyapunov exponent. We further propose that these trajectories admit a natural interpretation as Regge trajectories of stringy excitations in a dual stringy black hole geometry. From this perspective, pole-skipping for an individual operator can be viewed as tracking the stringy horizon through the response of a single excitation. Our results suggest that pole-skipping reflects intrinsic properties of quantum chaotic systems and may be deeply connected to the structure of horizons in the stringy regime.

Despite its successes, the large-$N$ holographic dictionary remains incomplete. Key features of gravitational path integrals–most notably Euclidean wormholes and the associated failure of factorization–lack a clear interpretation in the standard large-$N$ framework. A related challenge is the possibility of erratic $N$-dependence in CFT observables, behavior with no evident semiclassical gravitational counterpart. We argue that these puzzles point to a missing ingredient in the dictionary: a large-$N$ filter. This filter projects out the erratic $N$-dependence of CFT quantities when mapping them to semiclassical bulk physics, providing an intrinsic boundary definition of gravitational “averages.” It also offers a boundary explanation of wormhole contributions and a boundary prediction of their amplitudes, thereby giving a natural resolution of the factorization puzzle. In addition, we derive an infinite tower of inequalities constraining wormhole amplitudes and argue that internal wormholes do not induce random couplings in the low-energy effective theory. Beyond resolving factorization, the large-$N$ filter supplies a generalized framework from which richer Lorentzian spacetime structures can emerge, including closed universes and black hole interiors. We argue that, as a consequence of erratic large-$N$ behavior, both closed universes and black hole interiors are quantum volatile, and that an AdS spacetime entangled with a baby universe is likewise quantum volatile. This volatility may allow an observer in AdS to infer the existence of the baby universe, whereas for an infalling observer, the ability to make measurements near a black hole horizon may become fundamentally limited–even if they may not live long enough to notice.

The Lyman-alpha (Lya) forest is a unique probe of cosmology and the intergalactic medium at high redshift and small scales. The statistical power of the ongoing Dark Energy Spectroscopic Instrument (DESI) demands precise theoretical tools to model the Lya forest. We present a hybrid effective field theory (HEFT) forward model in redshift space that leverages the accuracy of non-linear particle displacements computed using the N-body simulation suite AbacusSummit with the predictive power of an analytical, perturbative bias forward model in the framework of the effective field theory (EFT). The residual noise between the model and the simulated Lya field has a nearly white (scale-and orientation-independent) power spectrum on quasi-linear scales, substantially simplifying its modeling compared to a purely perturbative description. As a consequence of the improved control over the 3D Lya forest stochasticity, we find agreement between the modeled and the true power spectra at the 5 per cent level down to scales of k <= 1 h/Mpc. This procedure offers a promising path toward constructing efficient and accurate emulators to predict large-scale clustering summary statistics for full-shape cosmological analyses of Lya forest data from both DESI and its successor, DESI-II.

The optical dump at the LUXE experiment has the potential to create a large flux of $\mathcal{O}({\rm GeV})$ photons that can be used to look for new physics when directed at a solid material dump. The LUXE-NPOD extension of LUXE, which focuses the hard photons onto a slab of tungsten, offers two interaction points (laser-electron and photon-tungsten), making it well-suited to test theories containing two or more new particles. We examine the dark axion portal, a scenario involving both dark photons (DPs) and axion-like particles (ALPs) and their mutual interactions, and its implications on the phenomenology at LUXE-NPOD. To simulate the spectra of particle populations generated at the electron-laser interaction point, we solve a set of extended 1D cascade equations. We recover a photon spectrum consistent with previous analyses and present previously unconsidered DP and ALP spectra. We derive the overall sensitivity of LUXE-NPOD to various parameters of the new particles, and show that it is capable of probing previously uncharted regions in the dark axion parameter space. For ALPs in the $10$-$500$ MeV mass range and DPs either heavier ($\sim{\rm GeV}$) or significantly lighter $\mathcal{O}(10^{-18}$-$10^{-13})$ eV, we obtain novel constraints on DP kinetic mixing parameters smaller than $10^{-2}$ and on $\mathcal{O}(10^{-4}$-$10^{-3})$ GeV$^{-1}$ ALP-photon couplings. We find that restrictions on $\mathcal{O}\left(10^{-3}\right)$ kinetic mixing can be extracted for arbitrarily small DP masses. Our discussion aims to be systematic and demonstrates a practical method of analyzing constraints on multi-dimensional parameter spaces.

We consider the propagation of heavy quarks through a dense, hydrodynamically flowing QCD medium, representative of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Working in the high-energy limit, we identify two novel mass-dependent effects arising from the heavy quark coupling to the local medium flow. The first is the emergence of a tensorial jet transport coefficient, $\hat{q}_{ij}$, which encodes the directional structure of transverse-momentum broadening. The second, named the directional dead-cone effect, corresponds to an anisotropic suppression of medium-induced radiation aligned with the hydrodynamic flow. We discuss how these effects manifest in jet observables and identify distinctive signature of heavy quark dynamics in an evolving medium.

Motivated by the requirement of the LSZ reduction formula to remove self-energy contributions on external legs, we examine quark self-energy contributions in soft-collinear effective (SCET) theory. We examine an operator basis that follows directly from full quantum chromodynamics (QCD) (upon application of the SCET equations of motion to express small Dirac components in terms of large Dirac components). We find that, for this basis, the self-energy contributions can be identified from their diagrammatic topologies, as in full QCD. However, for an alternative operator basis that is obtained from the direct-QCD basis by an application of Wilson-line identities, interactions are shifted from a covariant derivative to a Wilson line. Consequently, some self-energy contributions are hidden in diagrams involving Wilson lines, making their identification subtle.

We use a combination of perturbation theory, holography, supersymmetric localization, integrability, and numerical conformal bootstrap methods to constrain the energy-energy correlator in $\text{SU}(N_c)$${\mathcal N}=4$ SYM at finite coupling. For finite $N_c$, we derive lower bounds on the second and fourth multipoles of the energy-energy correlator at different couplings, along with a smeared energy-energy correlator as a function of the angle between the two detectors. We present evidence that our lower bounds on the multipoles are nearly saturated by the ${\cal N} = 4$ SYM theory. In the planar limit, we further use dispersive functionals to obtain tight two-sided bounds on both the first three non-trivial multipoles and on the angular dependence of the energy-energy correlator. As the coupling is varied from weak to strong, the energy-energy correlator exhibits a transition from single-trace to double-trace operator dominance in the collinear limit, which we characterize quantitatively. A similar phenomenon occurs in QCD, where a parton-hadron transition is observed as detectors are brought closer together.

Jet substructure observables serve as essential tools for probing the quark-gluon plasma produced in relativistic heavy-ion collisions. Their interpretation, however, is often complicated by edge effects, which arise when correlated particles fall outside the reconstructed jet radius, introducing distortions that obscure the underlying QCD dynamics. In this work, we present a comprehensive phenomenological study of edge effects in soft-insensitive angular observables, taking the two-point energy correlator (EEC) as a representative example. We argue that these distortions scale linearly with the average angular separation between the winner-take-all and $E$-scheme axes $\langle ϕ\rangle$, and validate this behavior across proton-proton (p-p) simulations with Pythia8 and Herwig7, as well as lead-lead (Pb-Pb) simulations using JEWEL and CoLBT. In p-p collisions, edge effects are strongly suppressed, scaling as $(R_L/R)^4$, whereas medium-modified jets can exhibit larger distortions, with contributions scaling as $(R_L/R)^2$ and $(R_L/R)^4$. Taking Pb-Pb/p-p ratios of the EEC substantially reduces, but does not completely eliminate, these distortions, highlighting the need of accounting for edge effects in the interpretation of heavy-ion jet substructure measurements. Since edge effects are largely governed by the $\langle ϕ\rangle$ distribution, studying this distribution provides a new handle for benchmarking and constraining the modeling of edge effects in heavy-ion event generators.

We consider planar $\mathcal{N}=4$ super Yang-Mills at finite temperature with chemical potentials that couple either to the $R$-charges or the spins of the operators. We find expressions for the Hagedorn temperatures at both zero coupling by explicitly counting states, and at strong coupling using the string theory dual. We then apply the quantum spectral curve (QSC) to this problem, which adds additional twists to the $Q$-functions. For a single chemical potential $μ$ coupled to one of the $R$-charges, we find the analytic weak-coupling Hagedorn temperature to one-loop order for any value of $μ$, and to two-loop order for $μ=1/2$. We then solve the QSC numerically, showing that at strong coupling there is good agreement with the string theory prediction to order $1/λ^{1/4}$. This provides further evidence for a recent conjecture of Harmark for the form of the world-sheet zero-point shift. We also use the QSC to find the analytic one-loop correction to the Hagedorn temperature with non-zero chemical potentials coupled to the spins.

Large superconducting quantum circuits have a number of important applications in quantum computing. Accurately predicting the performance of these devices from first principles is challenging, as it requires solving the many-body Schrödinger equation. This work introduces a new, general ab-initio method for analyzing large quantum circuits based on lattice field theory, a tool commonly applied in nuclear and particle physics. This method is competitive with state-of-the-art techniques such as tensor networks, but avoids introducing systematic errors due to truncation of the infinite-dimensional Hilbert space associated with superconducting phases. The approach is applied to fluxonium, a specific many-component superconducting qubit with favorable qualities for quantum computation. A systematic study of the influence of impedance on fluxonium is conducted that parallels previous experimental studies, and ground capacitance effects are explored. The qubit frequency and charge noise dephasing rate are extracted from statistical analyses of charge noise, where thousands of instantiations of charge disorder in the Josephson junction array of a fixed fluxonium qubit are explicitly averaged over at the microscopic level. This is difficult to achieve with any other existing method.

Supersymmetric localization and Ward identities have been used in the past several years to derive two integral constraints on the four-point function of the stress-tensor multiplet in $\mathcal{N} = 4$ super-Yang-Mills theory. These constraints are powerful tools for studying the theory, especially when used in tandem with analytic and/or numerical bootstrap techniques. In this paper, we consider three additional integral constraints that can be derived starting from the $\mathcal{N} = 4$ super-Yang-Mills theory placed on a squashed four-sphere. These constraints are technically much more challenging to derive than the ones in the literature, and much of the paper is devoted to developing techniques to make this computation tractable. Our end result is that these three constraints are implied by the two constraints already appearing in the literature.

Models of cosmic inflation generically predict a weak but potentially detectable amount of primordial non-Gaussianity (PNG), which can be used to obtain insights into the degrees of freedom during inflation and their interactions. The simplest types of PNG are the local and non-local (equilateral and orthogonal) shapes of the primordial three-point correlators, which are predicted by models with multiple light fields and derivative interactions in single-field inflation, respectively. In this paper we place constraints on local, equilateral, and orthogonal non-Gaussianities using the power spectrum and bispectrum extracted from first public release of the Dark Energy Spectroscopic Instrument (DESI). Our analysis makes use of higher-order clustering information through a consistent effective field theory (EFT) model for both the power spectrum and bispectrum at one-loop order. Using robust scale cuts where the EFT description is valid, we find the following constraints on PNG amplitudes: $f^{\rm loc}_{\rm NL}=-0.1\pm 7.4$, $f^{\rm equil}_{\rm NL}=719\pm 390$, $f^{\rm orth}_{\rm NL}=-200\pm 100$ (at $68\%$ CL). Non-local PNG constraints can be further improved by combining high-redshift DESI with legacy BOSS data and using simulation-based priors on bias parameters, yielding the strongest large-scale structure constraints to date $f^{\rm equil}_{\rm NL}=200\pm 230$, $f^{\rm orth}_{\rm NL}=-24\pm 86$. Our constraint on $f^{\rm loc}_{\rm NL}$ is competitive with the cosmic microwave background (CMB) limit; the combination gives $f^{\rm loc}_{\rm NL}=-0.0\pm 4.1$, $18\%$ stronger than the CMB only result, which represents the strongest bound on multi-field inflation yet obtained.

How is the production of primordial black holes (PBHs) in single-field models of inflation impacted by the presence of additional scalar fields? We consider the effect of a spectator field – a free scalar field with sub-Hubble mass, no direct coupling to the inflaton, and which makes a subdominant contribution to the total energy density – in the context of single-field models of inflation featuring a transient phase of ultra-slow roll (USR) evolution. Despite the modest title, a spectator field can have a dramatic impact: the slow-roll evolution of the spectator prevents the combined inflaton-and-spectator system from entering into USR, which naively might be expected to preclude the production of PBHs. However, we demonstrate that the growth of perturbations is maintained or enhanced by the spectator, through the rich interplay of curvature and isocurvature perturbations. We show in a model-independent way that the single-field phase of ultra-slow-roll is replaced by two turns in field space encompassing a phase of tachyonic instability for the isocurvature perturbations and a transfer of power from isocurvature to curvature modes. Furthermore, we highlight a degeneracy between the fine-tuning of the feature in the inflaton potential and the parameters of the spectator, leading to an overall resilience of model predictions to parameter variations. This makes it easier for the underlying PBH model to accommodate both high-precision CMB constraints and production of PBHs in the asteroid-mass range.