MIT Center for Theoretical Physics — a Leinweber Institute Publications

We study the stochastic dynamics of heavy quarks propagating through the strongly coupled plasma of $\mathcal{N}=4$ supersymmetric Yang-Mills (SYM) theory at nonzero temperature in terms of the corresponding Kolmogorov equation, which correctly describes their kinetic equilibration and the non-Gaussian fluctuations in their momenta without having to restrict their velocity to the non-relativistic regime. Leveraging the heavy quark limit, we show that the evolution of the momentum space distribution function can be reformulated as a Hamilton-Jacobi problem, and therefore can be solved in terms of first-order ordinary differential equations. We solve these evolution equations in an infinite thermal plasma with a constant temperature for initial conditions specified by spherically symmetric heavy quark momentum distributions with a phenomenologically motivated shape that is steeply falling at large momentum. To highlight the distinctive features of the kinetic equilibration process, we compare their solutions with Fokker-Planck dynamics with the same drag coefficient and fluctuations chosen by hand to guarantee equilibration. We find qualitatively similar dynamics at small momentum, and very different dynamics at large momentum, where, much like in jet quenching phenomena, the steepness of the momentum distribution gives a larger relevance to unlikely events in which a heavy quark loses little momentum in a given time step. Such events are much less unlikely in the Kolmogorov evolution than in Fokker-Planck evolution with the same mean energy loss, meaning that equilibration at large momentum is significantly delayed. Our results provide a systematic description of heavy quarks propagating through strongly coupled plasma from the ultra-relativistic to the non-relativistic regime and point the way towards implementation in phenomenological studies.

We present the next-to-next-to-leading logarithmic (NNLL) collinear resummation of projected $N$-point energy correlators (ENCs) up to $N=6$, matched to fixed-order predictions at NLO, in both electron-positron annihilation and Higgs decay to gluons. The key new ingredient is the two-loop jet function for $N=4,5,6$, which we compute semi-analytically using Integration-by-Parts and differential equations. We further include the leading non-perturbative corrections for ENCs, described by two universal soft matrix elements $\overlineΩ_{1q},\overlineΩ_{1g}$ of order $Λ_{\rm QCD}$, whose evolution is governed by anomalous dimensions for $(N-1)$-point correlators. The matched distributions are compared with parton-shower simulations from Pythia8 and Herwig7, and we study the sensitivity of both the absolute spectra and their ratios to the two-point energy correlator under variations of $α_s$ and $\overlineΩ_{1q,1g}$. Our results show that higher-point projected energy correlators are now under quantitative control at NNLL accuracy, opening the door to future $α_s$ extractions with complementary systematics.

The reliable treatment of highly boosted hadrons is crucial for many lattice QCD applications. For all these cases the kinematically enhanced interpolators promise very significant improvements and are, therefore, ever more often used in recent calculations, especially for highly boosted mesons like the pion. Motivated by, e.g., the physics program of the future Electron-Ion Collider (EIC) in the US and Electron-Ion collider in China (EIcC), we systematically benchmark our code for the unpolarized isovector nucleon quark matrix elements extracted at large source-sink separations, where excited state artifacts are significantly suppressed. We find that the precision of the renormalized nucleon matrix elements is typically improved by an order of magnitude at momentum $P_z\sim2.5$ GeV. By comparing the results from three CLS ensembles with different lattice spacings $a$ but the same pion mass, we observe no statistically significant dependence on $a$ in the renormalized matrix elements at nearly identical values of $P_z$. These encouraging results suggest that the use of kinematically improved operators is highly advantageous for parton physics calculations and can be extended to a broader class of baryon observables, making them a promising candidate for a standard component of modern lattice QCD.

Primordial black holes (PBHs) in the asteroid-mass range, with typical masses $10^{17}\,{\rm g}\lesssim M \lesssim 10^{23}\,{\rm g}$, have drawn significant recent attention as a viable dark matter candidate. The peak frequencies of photons emitted via Hawking radiation from asteroid-mass PBHs range from infrared to $γ$-ray bands. We calculate expected local transit rates for extended PBH mass distributions which could comprise all the dark matter. We evaluate prospects for detecting Hawking-radiated photons from local PBH transits through the inner Solar System and from PBH explosions in the far outer edges of the Solar System. We consider several existing and proposed ground-based and space-based instruments sensitive to photons from the radio band to ultrahigh energy $γ$-rays. We find that proposed instruments, such as the AMEGO-X satellite, could reliably detect PBH transits within ${\it O} (0.1 \, {\rm AU})$ of the Earth, while the HAWC and LHAASO observatories are both sensitive to PBH explosions out to ${\it O}(0.1 \, {\rm pc})$ and ${\it O}(0.5 \, {\rm pc})$ respectively. We conclude by specifically considering potential companion electromagnetic signatures in the case of a PBH explosion about $10^3\,{\rm AU}$ from Earth, which has been suggested as a potential source for the $\sim 220 \, {\rm PeV}$ ultrahigh-energy KM3-230213A neutrino event observed by the KM3NeT collaboration in 2023. Whereas we find that the recent KM3NeT event would not have yielded detectable electromagnetic signals — due to its location on the sky, proposed distance from Earth, and the offline status of the HAWC observatory at that time — we demonstrate that future PBH explosions at comparable distances could yield measurable electromagnetic signals at Earth, depending on alignment of the PBH burst with detector fields of view.

We investigate the interplay between the thermodynamic (Gibbs) entropy of quark-gluon plasma (QGP) and the quantum entanglement entropy characteristic of confined hadronic states across the quark-hadron phase transition. In the deconfined regime, entropy is well described by the statistical mechanics of colored quarks and gluons. Upon hadronization, however, the macroscopic Gibbs entropy of the plasma cannot simply vanish; instead, it is reorganized into the configurational entropy of a gas of colorless hadrons together with quantum correlations among the confined partons within each hadron. We show that the entanglement entropy of the internal partonic wave functions inside hadrons provides a natural repository for this “converted” thermodynamic entropy, reconciling the apparent reduction of macroscopic entropy with the second law of thermodynamics. Either by extrapolating from known facts about deep inelastic scattering, or starting from model descriptions of the proton wave function, or starting from the Hagedorn spectrum of its resonances, we provide three estimates of the magnitude of the entanglement entropy carried by a proton, with very different uncertainties. All three estimates indicate that the internal entanglement entropy of the proton is similar in magnitude to the Gibbs entropy of the QGP droplet from which the proton formed as QGP cools through the quark-hadron transition, as for example throughout the universe microseconds after the Big Bang. These results suggest that entanglement entropy offers a bridge between the quantum information content of hadronic states and the thermodynamic entropy of the quark-gluon plasma, providing a new lens on the microscopic mechanism of confinement and the nature of the QCD phase transition.

We revisit the Hartle-Hawking wave function in AdS spacetime, where natural spatial slices are open and require an additional spacetime boundary. This leads to two constructions: a fully gravitational wave function, in which the boundary configuration is integrated over, and a partially frozen one, in which it is fixed, as in AdS/CFT. To illustrate the fully gravitational construction, we explicitly analyze it in AdS$_3$ Einstein gravity and AdS$_2$ Jackiw-Teitelboim gravity. We then evaluate the one-loop correction to the hyperbolic-ball partition function in $D$-dimensional AdS Einstein gravity, expected to give the leading contribution to the wave-function norm. We demonstrate that the fully gravitational hyperbolic ball partition function, where the boundary fluctuates, develops a nontrivial one-loop phase of $(\mp i)^{D+1}$, analogous to that of the sphere partition function in dS gravity. By contrast, the partially frozen partition function, where the boundary is fixed, remains real and positive. Motivated by this AdS comparison, we conversely investigate a partially frozen dS sphere partition function where the metric on an equator is fixed, finding that its one-loop phase cancels nontrivially. Our results suggest that the phase problem is controlled by whether the gravitational path integral is fully dynamical or partially frozen.

Tidal stripping is a key feature of the evolution of dark matter (DM) halos, and has major implications for the population of low-mass galaxies. In the case of ultralight DM, tidal stripping proceeds not only classically, at the tidal radius, but also via a process analogous to quantum tunneling by long-wavelength particles out of the potential of a subhalo. This modified tidal stripping behavior leads to tight constraints on the particle mass as a function of subhalo and host properties. As many models of ultralight DM predict several independent species, it is crucial to understand how these constraints can be generalized to multifield halos with different particle masses. However, numerical challenges make it difficult to directly study the tunneling process in all but the simplest multifield scenarios. We introduce a simplified approach based on semiclassical methods that entirely sidesteps the most difficult aspects of the numerical problem, and we apply this to the study of tunneling in multifield halos. Our results significantly clarify the physics of tidal stripping for ultralight DM halos even in the single-field case: we provide first-principles derivations of features of the tunneling rate previously suggested by empirical fits. We then evaluate stability bounds on two-field halos for the first time, for a wide range of density and particle mass ratios. We show that for particular parameter combinations, the stability bounds in the two-field case can be somewhat relaxed relative to the single-field case, but for much of the parameter space, the constraints become more stringent. We discuss the path towards probing realistic multifield ultralight DM halos.

Today, the observable cosmos exhibits a remarkable degree of isotropy and plausibly began in a nearly isotropic initial state. The properties of the Lorentzian Chern-Simons-Kodama (CSK) functional can provide an understanding of this initial state. In gravity with a positive cosmological constant, the Chern-Simons-Kodama (CSK) wavefunctional is an exact, chiral solution of the quantum gravitational constraints. We suggest that the normalizability and other issues with this functional, if interpreted as a proper state of quantum gravity, instead suggest an embedding into a larger quantum gravitational completion, and recast the CSK functional as a gravitational sphaleron with observationally desirable properties. By perturbing around the dominant de Sitter saddle of the wavefunctional with appropriate quantum gravitational boundary conditions, we find that for a closed universe the system is dynamically driven to spatial isotropy, while all anisotropic modes acquire positive quadratic curvature and are Gaussian-suppressed. The decay of this sphaleron therefore proceeds along an isotropic channel, providing an intrinsic quantum-gravitational mechanism for dynamical isotropization. This isotropization effect is robust under the inclusion of a slow-roll inflaton, and no analogous isotropic sphaleron exists for spatially flat or hyperbolic geometries. Taken together, these results recast the Lorentzian CSK functional as a chiral sphaleron that naturally prepares an approximately isotropic de Sitter background for inflation. Beyond this phenomenological study, we further suggest that the CSK functional can be understood as a boundary functional for a class of anomaly-free objects, including a complexified generalization of the Hartle-Hawking state.

In his 1976 paper “More about the massive Schwinger model”, Coleman introduced $1+1$-dimensional Quantum Electrodynamics coupled to two charged massive fermions. By applying Abelian bosonization, he elucidated much of the physics of this two-flavor Schwinger model, but he listed three puzzles at the end of his paper. We present new analytical and numerical calculations to solve Coleman’s three puzzles and thereby deepen our understanding of this model. These puzzles pertain to the theory with equal fermion masses at $θ= 0$ and at $θ= π$, as well as the size of isospin-breaking effects when the fermion masses are unequal. For the puzzle at $θ= π$, the solution is related to the structure of the zero-temperature phase diagram arXiv:2305.04437: for equal fermion masses $m$, the model exhibits spontaneous breaking of charge conjugation symmetry and absence of confinement for any value of the gauge coupling $g$, so that there is a smooth interpolation from weak to strong coupling. Using two-loop Renormalization Group and integrability methods, we show that the mass gap behaves as $\sim m e^{-0.111 g^2/m^2}$ in the strong coupling regime $m\ll g$. Our numerical results using the lattice Hamiltonian are in good agreement with this behavior. For the puzzle at $θ= 0$, the solution is related to a level crossing between two isosinglet particles with different discrete quantum numbers; we demonstrate the necessity of such a crossing by comparing integrability and weak coupling calculations, and we also exhibit the crossing numerically. Finally, we provide a new estimate for the size of isospin-breaking effects caused by different fermion masses at strong coupling.

Recent developments have shown that some semiclassical spacetimes cannot emerge from a traditional application of the rules of holography, prompting proposals for restoring their emergence with “observer rules”. In this paper, we propose a general semiclassical diagnostic of such failures of emergence, and of the extent to which observer rules can fix them. Our diagnostic is the presence of certain “evanescent” quantum extremal surfaces, which are distinguished by an upper bound on their area rather than their generalized entropy. In particular, the generalized entropy of an evanescent QES may be large: even though its area term must be small, its bulk entanglement term is unconstrained. This feature is explained by an operational distinction between classical and quantum connectivity in semiclassical gravity, or equivalently between the two summands of the generalized entropy.

We prove a singularity theorem in semiclassical gravity without assuming global hyperbolicity or the null energy/curvature condition; the former is replaced by the weaker causality conditions of stable causality and past reflectivity, and the latter is replaced as is standard by the Generalized Second Law. This establishes in particular that the standard models of evaporating black holes are singular – i.e. they are null geodesically incomplete.

The annihilation of self-interacting dark matter with long-range interactions can be significantly enhanced at low velocities through the Sommerfeld effect. At special points in parameter space, where near-zero-energy resonances exist in the spectrum of the theory, the standard calculation of this enhancement appears to violate unitarity. Recently, several approaches have been proposed to regulate this behavior, some introducing explicit ultraviolet (UV) scales and others not, raising the question of whether these prescriptions are consistent. In this paper, we compare these approaches and show that even in nominally cutoff-dependent methods, the regulated Sommerfeld-enhanced cross sections are independent of the UV regulator to a good approximation, and that when the unitarity-preserving corrections are large, the regulation schemes coincide to leading order. We use these insights to write down a regulator-independent prescription for unitarization applicable to multi-state systems, where the modified enhancement can be written solely in terms of the standard enhancement factor, the hard annihilation amplitude, and the $S$-matrix for scattering in the long-range potential.

The universality of the large momentum expansion allows computing parton distribution functions (PDFs) starting from any Euclidean correlator with appropriate large momentum Fourier Components. Here we consider current-current correlators which have been used in short-distance expansion to obtain moments of PDFs. The advantage of such correlators is that they have simple renormalization properties and do not have linear power divergences as in quasi-PDF. However, in lattice calculations, four-point functions are needed. Here we present an expansion formula with current-current correlators up to the next-to-leading order, and preliminary numerical calculations with four-point functions.

Jet production at the LHC is a powerful probe of QCD, making it ideal for precision tests and determinations of QCD parameters such as parton distribution functions and the strong coupling constant. To make the most of the abundant jet production data collected at the LHC, precise calculations are required. While state-of-the-art calculations reach next-to-next-to-leading order (NNLO) QCD accuracy, a critical assessment of the remaining uncertainties arising from non-perturbative effects and missing higher orders remains crucial for correctly interpreting comparisons between theory and data. Scale variation is nearly always used to determine effects from missing higher orders. In this article, we reassess this method in the context of inclusive jet production by performing NNLO QCD calculations supplemented by small-jet-radius resummation through next-to-next-to-leading-logarithmic accuracy (NNLL). We find that NNLL resummation can have an appreciable impact on the scale uncertainty for inclusive jet cross sections, and, for some scale choices, can lead to sizeable shifts of the central cross section. We conclude that scale variations in fixed-order and resummed calculations can drastically underestimate the impact of higher orders for commonly used jet radius parameters, and that missing higher-order estimates obtained via scale variations should be considered unreliable. Our findings add further evidence to the importance of going beyond scale variations in jet and jet substructure calculations.

We present an amplitude analysis of high-energy polarized photoproduction of $π^-Δ^{++}$ within a Regge exchange framework. A Regge amplitude model incorporating $π$, $ρ$, $b_1$, and $a_2$ trajectory exchanges is fit simultaneously to spin density matrix elements measured by the GlueX experiment at photon energies of $E_γ= 8.2$–$8.8$ GeV and differential cross section data from SLAC. By including SDME data, the fit constrains not only the magnitudes but also the relative phases of the helicity amplitudes. The results confirm the dominance of pion exchange at small momentum transfer, while natural parity exchanges become significant at larger $t$. We analytically continue the $s$-channel amplitude to the $t$-channel, taking care of the kinematical singularities, and isolate the dynamical residues at the meson poles. The extracted $πNΔ$ coupling constant is found to be consistent with the value obtained from the decay width of the $Δ(1232)$. For the $ρNΔ$, $b_1 NΔ$, and $a_2 NΔ$ vertices, first extractions of the relevant coupling constants are provided.

We find that heavy quark transport beyond leading logarithm at weak coupling is intrinsically non-Gaussian: the longitudinal momentum transfer distribution has asymmetric exponential tails that are crucial for equilibration dynamics. We show this by computing the leading-order momentum transfer kernel for relativistic heavy quarks in weakly coupled non-Abelian plasmas, matching perturbative momentum transfer on the thermal scale to hard-thermal-loop-resummed soft physics. This is the same structure previously found in strongly coupled holographic plasmas, showing that it is not peculiar to weak or strong coupling, conformality, or supersymmetry. We therefore expect that this is a robust feature that physical quark-gluon plasma should also exhibit.

Quantum error correction is widely believed to be essential for large-scale quantum computation, but the required qubit overhead remains a central challenge. Quantum low-density parity-check codes can substantially reduce this overhead through high-rate encodings, yet finite-size instances with practical logical error rates often achieve encoding rates only around or below $1/10$. Here, building on a recent ultra-high-rate construction by Kasai, we identify new structural conditions on the underlying affine permutation matrices that make encoding rates exceeding $1/2$ compatible with efficient implementation on reconfigurable neutral atom arrays. These conditions define a co-designed family of ultra-high-rate quantum codes that supports efficient syndrome extraction and atom rearrangement under realistic parallel control constraints. Using a hierarchical decoder with high accuracy and good throughput, we study the performance under a circuit-level noise model with $p=0.1\%$, achieving per-logical-per-round error rates of $1.3_{-0.9}^{+3.0} \times 10^{-13}$ with a $[[2304,1156,\leq 14]]$ code and $2.9_{-1.5}^{+3.1} \times 10^{-11}$ with a $[[1152,580,\leq 12]]$ code. These results approach the teraquop regime, highlighting the promise of this code family for practical ultra-high-rate quantum error correction.

Energy correlators offer a clean probe of quantum chromodynamics, serving as an ideal laboratory to rigorously investigate non-perturbative power corrections. The recent discovery that linear corrections exhibit a universal anomalous scaling points to a deep, underlying theoretical structure. We uncover the quantum field-theoretic origin of this phenomenon in the energy-energy correlator using light-ray operators. Through an explicit loop calculation, we derive the one-loop anomalous dimension, revealing that the dijet operator must be combined with a specific triple-jet component. This provides a first-principles framework that connects operator theory with high-precision collider phenomenology.

The 1D flux power spectrum ($P_{\mathrm{1D}}$) of the Ly$α$ forest provides an exceptionally high-resolution probe of structure formation down to small scales ($k\approx1-10~\text{$h~$Mpc$^{-1}$}$). These scales carry the imprints of massive neutrinos, warm dark matter, and the running of the primordial power spectrum spectral index. The effective field theory (EFT) is a promising perturbative approach to systematically and efficiently describe the Ly$α$ forest, but it faces challenges in its application to $P_{\mathrm{1D}}$, as many EFT parameters become degenerate when projected along the line of sight. In addition, this projection generates new stochastic terms from the integration over small-scale modes. In this work, we address these issues by compressing the EFT model space using the Fisher matrix formalism and linearizing the resulting compression directions, enabling analytic template marginalization and significantly reducing the computational cost of likelihood evaluation. We use hydrodynamical simulations to obtain a baseline estimate of EFT parameters, and use the DESI DR1 $P_{\mathrm{1D}}$ measurements to derive compression directions. We then marginalize over deviations from the baseline using these compression directions and forecast the constraining power of our formalism. We find that even in conservative scenarios where each data redshift bin requires its own set of EFT parameters, the cosmological constraints saturate with the linear bias, two leading-order 1D stochastic terms, and three principal combinations of the remaining EFT templates. In this case, our forecasted precision of the amplitude ($Δ^2_p$) and the logarithmic slope ($n_p$) of the linear matter power spectrum at the pivot scale ($k_p=0.7~\text{Mpc}^{-1}$) is $10\%$ and $2.0\%$, respectively, which is similar to emulator-based analyses that include observational data systematics.

We present a solvable Hamiltonian that realizes an exact lattice chiral $U(1)_V \times U(1)_A$ symmetry. Nielsen-Ninomiya-type no-go theorems are evaded by using lattice bosons rather than fermions. The continuum limit is a compact boson field theory with an axion-like coupling. The $U(1)_V$ symmetry shifts the scalar, while $U(1)_A$ acts on local operators associated with short axion strings and is transmuted into a higher-form symmetry in the continuum limit. We demonstrate the chiral anomaly by showing that the lattice theta angle is shifted by an axial rotation when $U(1)_V$ is gauged. Gauging either $U(1)_V$ or $U(1)_A$ leads to lattice non-invertible and 2-group symmetries, respectively, matching the continuum picture.