MIT Center for Theoretical Physics — a Leinweber Institute Publications

Quantum algorithms enable Heisenberg-limited scaling across many tasks in sensing and learning, achieving precision $varepsilon$ using $O(1/varepsilon)$ total evolution time. However, in Hamiltonian property testing, the best known algorithms have $O(1/varepsilon^2)$ total evolution time, unless it is possible to $textit{reverse}$ the flow of time under the Hamiltonian. In this work, we prove that time reversal is necessary to close this gap. We establish evolution time lower bounds for Hamiltonian certification and locality testing, for a continuous family of error metrics ranging from average- to worst-case. We also prove matching upper bounds for Hamiltonian certification in all cases, both with and without time reversal, showing that the advantage from time reversal persists. Our results follow from a new Fourier-analytic approach to quantum lower bounds in continuous time, which is sensitive to time reversal and applicable beyond Hamiltonian property testing. This gives an alternative proof that inverses are required for amplitude amplification and estimation, improving the lower bounds of Tang and Wright and extending to fractional queries. We also show that, for unstructured search, uncertainty in the phase of the oracle destroys the quantum speedup without inverse queries, even when the oracle is fixed and coherent. Our work rigorously establishes time reversal as a resource governing both the precision of quantum measurement and the robustness of quantum algorithms.

Relativistic heavy quarks equilibrate even when their microscopic drag and momentum diffusion violate the Einstein relation. We investigate how equilibration proceeds as a function of the coupling and compared to simplified Fokker-Planck models. To do this, we present a detailed calculation of the momentum-transfer kernel in weakly coupled non-Abelian plasmas through strict $mathcal{O}(g^4)$, and study the resulting dynamics at leading order in the inverse heavy-quark mass. Beyond leading logarithm, its non-Gaussian structure plays a crucial role in equilibration. We examine closely the analytic structure of the evolution kernel, determine the large-order asymptotics of the momentum-transfer cumulants, and show how the kernel's bounded analytic domain fixes asymmetric exponential tails at the level of the momentum transfer probability. Strongly coupled $mathcal{N}=4$ SYM also features an approximately Gaussian core and exponential tails, but differs in the nature of the kernel's singularities and the velocity dependence of the transport coefficients and tail exponents. These differences have distinct dynamical consequences: for the parameters studied, the persistence of the relativistic population in steeply falling spectra is largely explained by drag alone at very weak coupling, whereas fluctuation-induced survivor bias produces a large relative enhancement at strong coupling. Nevertheless, the shared qualitative structure between weak and strong coupling points the way towards a practical phenomenological description of heavy quarks anchored in first-principles field theory.

We present a technique for computing Drell-Yan lepton pair production at particle-level that achieves N$^3$LO and approximate N$^4$LL$'$ accuracy at small to moderate transverse momentum, tree-level precision for $Z+4$, $5$, and $6$ jets, and approximate NLO EW precision at high transverse momentum. A crucial non-perturbative input, the Collins-Soper kernel, is obtained from continuum-limit lattice QCD results. High-precision theory results are recast into strictly positive event weights using an information-theoretic approach based on the maximum entropy principle.

Discretizing electronic degrees of freedom into atomic orbitals is the default choice to accurately model electrons in molecules in a compact and flexible way. We consider the task of simulating such electronic structure models on a quantum computer for $N$ orbitals, keeping the number of orbitals per atom constant. Atomic orbitals give rise to Hamiltonians with potentially up to $mathcal{O}(N^4)$ terms and little structure, which made it challenging to match the complexity of methods based on plane waves and real-space grids. In this work, we show that a Trotter step for atomic orbitals can be implemented in depth $mathcal{O}(mathrm{polylog}(N))$ without ancillas by combining a rigorous treatment of orbital localization, a hierarchical Hamiltonian decomposition based on the Fast Multipole Method, and shallow quantum Fourier arithmetic circuits. By bounding the Trotter error, we obtain sublinear simulation depth, and our total gate count $N^{5/3 + o(1)}$ matches the best scaling known for any basis.

We prove that the Gibbs state of any spin chain can be exactly decomposed into a mixture of product states over contiguous blocks of qubits with bounded length, at any constant temperature. As a consequence, any such Gibbs state can be prepared by local unitary circuits of constant depth, with efficient classical preprocessing. We place strong limits on the spatial structure of thermal entanglement: both the entanglement depth and entanglement width (measures of multiparticle entanglement) are bounded by constants, while the localizable entanglement (the bipartite entanglement that can be generated with the aid of local measurements and classical communication) vanishes exactly beyond a constant threshold distance. Our result also has implications for quantum thermalization: for any translation-invariant spin-chain Hamiltonian with nondegenerate spectral gaps, typical low-complexity pure states drawn from a maximally entropic ensemble at any constant effective temperature become locally indistinguishable from the Gibbs state upon unitary evolution. Additionally, we obtain an analogous decomposition for fermionic systems: the Gibbs state of any 1D local fermionic Hamiltonian is exactly a mixture of constant-depth local fermionic circuits acting on Gaussian states.

Learning and sensing are among the most promising applications of quantum technology, often with provable quantum advantages. In any realistic experiment, however, noise threatens to erase these advantages as the problem size grows. We introduce oracle distillation, a procedure that distills many queries to a noisy oracle into a single high-fidelity oracle. Here the oracle, an unknown unitary from a known family, models the interface through which a quantum device learns about an unknown system. We construct an explicit protocol that distills any Boolean oracle, both under adversarial noise that corrupts a constant fraction of the input qubits and under i.i.d. depolarizing noise at a constant rate on every qubit. At the heart of the protocol is the weak query, in which the device deliberately queries the oracle only weakly, trading response strength for the ability to correct errors. Our protocol is near-optimal among all protocols that distill the oracle while keeping its errors correctable, revealing a fundamental tension between these two requirements. An important implication of the protocol is a threshold theorem: every Boolean oracle problem whose quantum advantage is polynomial in the domain size, including Grover search, $k$-forrelation, and Simon's problem, retains a quantum advantage whenever the per-qubit depolarizing rate is below a constant threshold. In particular, the quadratic quantum advantage of Grover search, long believed fragile under noise, survives nearly intact at low error rates even when almost every query involves errors on one or more qubits. We further extend the protocol to fractional and continuous-time Boolean oracles. Quantum error correction can thus protect not only prescribed operations but also unknown dynamics. Our results make a large family of quantum advantages in learning robust against noise and open a path toward robust computational sensing.

The quantum circuit complexity of an evolving many-body quantum system is believed to exhibit a sustained growth, maintained for timescales much longer than the onset of thermalization. Most previous works have focused on models which violate energy conservation, such as random unitary circuits. Here we study generic, local time-independent Hamiltonian dynamics. We unconditionally prove that for generic local Hamiltonians and typical initial product states at high effective temperature, the robust quantum circuit complexity must grow over a very long period of time, attaining an exponentially large value at late times. Our approach relies on two structural properties that we prove rigorously: (i) generic local Hamiltonians satisfy generalized spectral no-resonance conditions of arbitrary order, and (ii) typical high-temperature product states are effectively supported in exponentially many energy eigenstates. These properties have been widely assumed without proof in prior works. As corollaries of this result, we show the late-time state displays robust volume-law entanglement that is irremovable by polynomial-size circuits, and we establish a no fast-forwarding result for generic local Hamiltonians. We also lower bound the complexity of preparing sufficiently large subsystems with local quantum channels, in contrast with sufficiently small thermalizing regions which we show retain low complexity at late times.

Decoherence in quantum systems is conventionally modeled as the effect of interactions with an external environment. However, such a prescription excludes isolated many-body systems, which are also expected to display classical behavior at macroscopic scales. In isolated systems, decoherence must emerge internally from microscopic degrees of freedom that are invisible to the macroscopic description. Here we explicitly show that classicality can emerge in such a fashion. We consider a weakly disordered, $3$-local chaotic kicked top of $N$ qubits, where the collective spin sector serves as the macroscopic description, while the microscopic permutation sector acts as an internal bath, decohering the collective spin sector. Starting from closed unitary dynamics, we derive and numerically confirm an effective Lindblad equation for the collective spin variables. In the thermodynamic limit these reduced dynamics converge to a classical chaotic Fokker--Planck equation with vanishingly small diffusion on the spherical phase space, producing a quantum-classical correspondence beyond the Ehrenfest time. The chaotic dynamics evolve the pure many-body wavefunction into continuously branching components associated with distinct classical trajectories. These branches acquire nearly orthogonal microscopic records in the permutation sector, preventing quantum interferences and ensuring the corresponding histories remain consistent.

Charge-dependent directed flow provides a sensitive probe of early electromagnetic fields and baryon stopping in relativistic heavy-ion collisions. Recent STAR measurements show a centrality-dependent sign change in the directed flow splitting (the difference between the directed flow of protons and antiprotons), indicating that electromagnetic effects alone are not sufficient to describe this observable and that the baryon stopping $-$ in particular the component of the stopped proton distribution that is odd in rapidity and odd under reflection in the impact parameter direction $-$ must also be included. We develop a semi-analytic hydrodynamic framework that combines spectator-induced electromagnetic fields with a Glauber-based description of baryon stopping, built upon an analytic solution for the background hydrodynamic flow due to Gubser together with the simplifying assumption of a constant electrical conductivity. For Au+Au collisions at $sqrt{s_{NN}}=200$ GeV, we find that baryon stopping gives a positive contribution to the directed flow splitting that decreases for more peripheral collisions, while electromagnetic fields give a negative contribution that is larger for more peripheral collisions. The competition between these two effects naturally reproduces the observed sign change in Au+Au collisions as a function of centrality, describes the observed rapidity dependence in the $50$--$80%$ centrality interval, and reproduces trends seen in U+U collisions. Although the simplifying assumptions that we have made regarding the analytic background and constant conductivity limit our ability to make quantitative comparisons, our model provides a transparent explanation of how transported baryon number and spectator-induced electromagnetic fields jointly shape charge-dependent directed flow.

The frequency spectrum of the cosmic microwave background is the most precise blackbody ever measured in nature, with deviations constrained at the level of almost one part per million from the COBE satellite. Nevertheless, departures away from a perfect blackbody are present in standard $Λ$CDM cosmology, lurking just beneath the surface of our current observational bounds. These spectral distortions provide invaluable information on our thermal history in both the post- and pre-recombination epochs, allowing us to peer beyond last scattering and into the primordial Universe. Here, we present an overview of the underlying physics responsible for generating CMB spectral distortions at all epochs. As an illustration of this rich physics, we review a comprehensive set of mechanisms capable of generating distortions both within and beyond the standard $Λ$CDM paradigm. We also discuss the information that can be gleaned by going beyond the monopole (sky-averaged) spectrum and exploiting the spatial information present in anisotropic spectral distortions. To supplement our discussion of the diverse science of spectral distortions, we provide an overview of the upcoming and proposed experimental landscape. We highlight that the combination of the TMS, COSMO, and BISOU experiments will provide the first discovery of a monopole $y$-type distortion within the coming decade. From space, the proposed FOSSIL experiment is forecasted to improve upon the original COBE/FIRAS measurement by roughly three orders of magnitude in sensitivity, bringing with it the detection of the $Λ$CDM $μ$-type distortion sourced by the dissipation of small scale acoustic modes in the pre-recombination plasma. With transformational measurements on the horizon, CMB spectral distortions offer a uniquely sensitive probe of the thermal history of the Universe at redshifts $z lesssim 2 times 10^6$.

Recent progress in collider physics has reformulated phenomenological questions in terms of detector correlation functions and advanced their theoretical understanding. These advances have primarily focused on correlators of the average null energy operator, $mathcal{E}(vec{n})$, known as energy correlators. However, colliders have access to a much broader class of detector operators, $mathcal{E}^n_R(vec{n})$, which measure powers of the energy on a subset of hadrons $R$, such as charged hadrons. These generalized detectors are generically not infrared and collinear safe, and their description requires nonperturbative matching between the hadronic detectors measured in the infrared and the partonic detectors used in ultraviolet calculations. We develop a framework for computing their multi-point correlation functions. We introduce universal nonperturbative matching coefficients, termed "detector functions", that implement this infrared-ultraviolet matching. For the $mathcal{E}_R^n$ operators studied here, these coefficients are represented by energy-weighted moments of single- and multi-hadron fragmentation functions. We present their renormalization group structure, derive QCD factorization theorems for the projected correlators, and compute jet functions through next-to-leading order. In the fixed-coupling pure Yang-Mills limit, we derive the light-ray OPE of hadronic detectors and connect it to the QCD factorization framework. We also identify universal nonperturbative power corrections generated by soft radiation, which are enhanced in the collinear limit and can substantially modify the perturbative angular scaling. A parton shower study finds qualitative agreement with the predicted perturbative and nonperturbative scaling behaviors. Our work significantly broadens the space of detector operators under theoretical control, with potential phenomenological applications.

Motivated by searches for weak coherent drives, we formulate repeated quantum sensing with a fixed shot budget as an asymmetric composite hypothesis test. Taking Rabi sensing as a concrete example, we benchmark detection power, sensitivity, and Type-II error exponents in the resonant case with unknown signal amplitude and phase. We compare non-adaptive population and transverse readouts with a myopic Bayesian policy that selects each projective axis by maximizing the expected information gain in one step. The common decision statistic is a log Bayes factor, with a policy specific threshold calibrated under the null to enforce a common Type-I error. A weak signal expansion shows that population readout is phase independent but quadratic in amplitude, giving $n^{-1/4}$ sensitivity, whereas transverse readout is linear in amplitude and permits $n^{-1/2}$ sensitivity without adaptation, but is phase-dependent. In Monte Carlo pseudoexperiments, the adaptive policy exploits posterior information about the unknown direction to guide subsequent readouts; its sensitivity is consistent with $n^{-1/2}$ over the simulated large $n$ range and, at the largest simulated shot counts, outperforms the fixed transverse schedules. Its phase-averaged effective Type-II exponent also exceeds the non-adaptive references over the simulated range. These results demonstrate the finite budget value of exploiting nuisance parameter information under calibrated false positive control.

The kinetic Sunyaev-Zeldovich (kSZ) effect is a sensitive probe of the free-electron momentum field, tapping into spatial correlations relevant for both large-scale cosmology and smaller-scale astrophysics. Precision kSZ measurements will provide a robust channel for constraining inflationary physics and equivalence principle violation, while simultaneously characterizing the small-scale gas distribution and associated baryonic feedback. However, robust kSZ measurement for cosmological applications is challenging, as the primary source of large-scale information is convolved with unknown small-scale dynamics in a degenerate combination. In this paper we address this challenge using the Effective Field Theory (EFT) framework that provides a principled way of treating the effects of small-scale dynamics on large scale-structure. Extracting the kSZ signal can be reformulated as a measurement of the one-loop galaxy-galaxy-electron momentum bispectrum. We compute here a pure dark matter version of this bispectrum in EFT and compare it to simulations. We derive the EFT counterterms that capture small-scale backreaction and find that two of them acquire a transverse direction. One transverse counterterm gets generated entirely by the vorticity of the velocity field. Our computation of the bispectrum dipole agrees with simulation data to 5% up to k_max = 0.23 h/Mpc, which is significantly larger than the reach of tree level theory, which breaks down beyond k_max = 0.07 h/Mpc (at redshift z = 0.5). For the bispectrum computed with the projected momentum field, relevant to measuring the kSZ field, we find a similar k_max, and make the first detection of the transverse counterterms from simulations. Our results suggest that one-loop bispectrum computations can play an important role in advancing kSZ cosmology.

We carry out a search for the relative baryon-cold dark matter (CDM) density and velocity modes in the DESI full-shape (FS) data, and investigate whether DESI's preference for evolving dark energy can be affected by these relative perturbations, which are a known contaminant of the baryon acoustic oscillations. For the first time, we use a holistic effective field theory (EFT) treatment of the power spectrum and bispectrum both for the standard terms (without relative perturbations), modeled in EFT to one-loop order, and for the relative perturbations. For the latter we implement an efficient computation of the one-loop power spectrum that allows us to scan over different cosmologies. We obtain robust constraints on the relative density and velocity bias parameters $b_{δ_{bc}}$, and $b_{v^2_{bc}}$ for all the DESI DR1 tracers, finding e.g. $b_{δ_{bc}}=5 pm 6$, $b_{v^2_{bc}}=-0.031pm 0.023$ (at 68% CL) for the LRG2 sample in a DESI FS-only analysis within $Λ$CDM, supplemented by CMB priors on the baryon density and spectral tilt. Assuming instead a $w_0w_a$CDM cosmological model, and combining DESI FS with the CMB primary anisotropies, CMB lensing, and Pantheon+ supernovae, our LRG2 constraints read $b_{δ_{bc}}=2 pm 5$, $b_{v^2_{bc}}=-0.032pm 0.022$. Adding the relative velocity divergence biases we find $b_{δ_{bc}}=-17_{-17}^{+19}$, $b_{θ_{bc}}=-760_{-680}^{+700}$, $b_{v^2_{bc}}=-0.032pm 0.024$ for the same sample. The constraints on the dark energy equation of state parameters are virtually unchanged when marginalizing over the baryon-CDM modes. We conclude that DESI's preference for evolving dark energy is robust with respect to the baryon-CDM relative perturbations.

The Sommerfeld enhancement (SE) from long-range interactions, and the related bound state formation rate, can be important non-perturbative inputs to dark matter (DM) annihilation signals. In the presence of singular interaction potentials, the conventional boundary conditions of the Schrödinger equation fail, obfuscating the computation of SE, its physical origin, and its relation to bound states in such potentials. In this work, we clarify the origin of SE in singular potentials in a two-fold manner: using the framework of velocity power counting in non-relativistic effective field theory (NREFT), and via position-space regularization of singular potentials at short distances to compute SE and bound states. We illustrate our findings through the case of pseudoscalar mediators interacting with massive Dirac DM. We find that when such a system arises from a UV-complete theory at weak coupling, no SE is generated. However, in the case of a derivatively coupled pseudoscalar, where the interaction is described by a higher-dimension operator in an effective theory, SE (and bound states) can occur for weak couplings if there is a large hierarchy between the cutoff scale of the effective theory and the dark matter mass. This SE is sensitive to the choice of the UV completion of the potential at short distances, but a non-negligible SE can persist even when the UV physics alone would not generate any SE; we elucidate the interplay of UV and IR physics in this case.

In a broad class of inflationary models, the leading non-Gaussian bispectrum arises from one-loop rather than tree-level processes. However, realistic one-loop contributions remain largely unexplored for phenomenologically relevant masses beyond the simplest bubble diagrams. This is insufficient for observational purposes because, for generic masses and couplings, triangle contributions need not be suppressed relative to their bubble counterparts. We present a numerical method for evaluating scalar one-loop in-in diagrams contributing to the inflationary bispectrum at fully general external momenta and masses. First, the Witten-Feynman parameterization reduces de Sitter loop integrals to generalized Euler-Mellin integrals governed by Symanzik graph polynomials. Unfortunately, the exponents of the polynomials in the integrand become complex for sufficiently heavy masses, leading to sign problems when evaluating the integral using standard numerical techniques. To remedy this problem, we introduce the reduced Schwinger method, which evaluates a highly oscillatory subintegral analytically to obtain a Gauss hypergeometric kernel, leaving the remaining integrals to standard numerical quadrature. We validate the algorithm by reproducing known analytic results for the tree-level bispectrum in terms of ${}_3F_2$ functions, and then apply it to the one-loop bubble and triangle contributions at general kinematics. This yields the first direct numerical evaluation of the complete one-loop bispectrum, valid across the full kinematic range. Using these numerical results, we construct bispectrum templates and compare them with CMB data.

We solve in closed form a one-dimensional relativistic system: two masses interacting only through elastic collisions with a massless mediator bouncing between them. Momenta, times, and positions are hyperbolic functions of the collision index. The mediator energy, interpreted as the pair's effective potential, obeys an exact discrete Coulomb law, $Vpropto 1/r$, with a Lorentz-invariant action as coupling. A massive Newtonian mediator instead transmits a $1/r^{3}$ force; one adiabatic invariant traces both laws to the mediator's dispersion relation. Continued to negative mediator energy, the closed forms turn trigonometric, binding a one-dimensional mechanical analog of the Coulomb atom.

We study the growth of massive $SO(9)$ super-representations in type II and type I string theories. We do this first directly by finding an empirical formula that specifies which representations appear at any level, and then compute the multiplicities from a refined partition function evaluated over finite fields up to level 501 for individual representations, and level 226 for all representations. We then derive asymptotic formulae for the growth of any representation, which are constructed by integrating about the peaks of the refined partition function. Using a Rademacher sum we can find very accurate approximations for the multiplicities of the representations. Unlike the superstring partition function, which only receives contributions from the odd Rademacher terms, the multiplicities for any representation have contributions from both even and odd terms. Finally, we apply the same techniques to the ordinary representations, where a simplified refined partition function allows for better computational speed and simpler expressions for the asymptotic approximations.

Quantum thermalization explains how an isolated many-body system naturally evolves towards a thermal state, rendering information about the initial conditions inaccessible to local measurements. This is precisely the mechanism utilized in quantum error correction, where information is protected by design through a nonlocal encoding. In this work, we leverage this connection to port the rigorous framework of (approximate) quantum error correction to the study of quantum thermalization. Treating typical late-time states as codewords, we characterize the error-correcting properties of generic thermalizing dynamics. We numerically uncover a universal relationship between the encoding rate, distance, and thermal entropy density of the emergent code. At infinite temperature, this universal curve saturates the quantum Singleton bound, achieving the same optimal limit as Haar-random codes. At finite temperature, we introduce a code family based on the Scrooge ensemble, the natural thermal analogue of the Haar ensemble, and prove it saturates the entropic quantum Singleton bound, establishing this family as optimal within entropic constraints. Our extracted universal curve independently saturates this same bound, revealing that finite-temperature thermalization is itself optimal. Finally, we show how conserved quantities limit the error-correcting behavior of thermalization: codewords with differing energies, or other conserved charges, leak only classical information, and correctability persists until the difference reaches the scale of thermal fluctuations. Our results reveal a universal optimal coding structure in thermalizing dynamics, while introducing new optimal codes that achieve fundamental limits of approximate quantum error correction.

We present predictions for the collinear energy-energy correlator (EEC) measured inside $Z$-tagged jets in proton-proton collisions across the confinement transition. Our framework combines transverse-momentum-dependent factorization for back-to-back $Z+$jet production with an exclusive differential EEC jet function incorporating perturbative evolution and flavor-dependent non-perturbative fragmentation. The EEC jet function is evaluated at $mathrm{NLL}^{prime}$ accuracy, with separate quark and gluon non-perturbative scales constrained by existing $e^+e^-$ and inclusive-jet data. We provide predictions for forward $Z+$jet production at $sqrt{s}=8~mathrm{TeV}$ and $13~mathrm{TeV}$ in the jet transverse-momentum intervals $(20, 30)$, $(30, 50)$ and $(50, 100) , mathrm{GeV}$. The EEC distributions as functions of $z_χ$ exhibit finite plateaus at small $z_χ$, while the corresponding angular distributions display the characteristic turnover associated with the confinement transition. A flavor decomposition shows that this transition becomes increasingly dominated by quark-initiated jets as the jet transverse momentum increases. Forward $Z$-tagged jets therefore provide a clean probe of non-perturbative quark fragmentation and offer a direct test of whether the characteristic confinement scale extracted from $e^+e^-$ collisions remains universal in a hadronic environment.