MIT Center for Theoretical Physics — a Leinweber Institute Publications

The pion light-cone distribution amplitude (LCDA) is an essential non-perturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of $m_π= 550$ MeV. This calculation finds $\langleξ^2\rangle = 0.202(8)(9)$ and $\langle ξ^4 \rangle = 0.039(28)(11)$ where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the $\overline{\text{MS}}$ scheme at a renormalization scale of $μ= 2$ GeV.

Measuring the energy loss of energetic jet partons may provide experimental opportunities to constrain the initial nonequilibrium stages in heavy-ion collisions, which requires theoretical predictions for the jet-medium potential. In this letter, we go beyond the usual harmonic approximation, extract this potential for the first time using QCD kinetic theory simulations, and provide a simple cutoff-independent small-distance form. We find significant differences to thermal media at early times in momentum and position space, as well as substantial angular dependence. Applying our results to the gluon splitting rates using a novel computational method reveals tensions in the current kinetic theory implementations for early times.

Quantum systems of physical interest are often local, but there are at least three competing perspectives on how “locality” should be formalized: an algebraic framework, a path-integral framework, and a lattice framework. One puzzle in this competition is that systems with higher-form symmetries, which are perfectly local from the path-integral and lattice perspectives, can violate the algebraic principle of “additivity”. In this paper, we propose a resolution to this puzzle by introducing a weaker locality principle, “disjoint additivity”, which together with Haag duality should always be satisfied in local quantum systems. As evidence, we give examples in which disjoint additivity is preserved when ordinary additivity is violated; we show that Haag duality and disjoint additivity are satisfied in rather general lattice systems with local symmetry constraints; we give examples of nonlocal theories in which either disjoint additivity or Haag duality is violated; and finally we give examples of systems with nonlocal symmetry constraints in which disjoint additivity is violated, but can be restored by passing to a local “SymTFT” system in one higher dimension.

We present simulations of the real time dynamics of a fluid in the vicinity of a critical endpoint in the phase diagram. The relevant hydrodynamic theory is known as model H, and it is expected to describe the long-distance dynamics of QCD matter in the vicinity of a possible critical endpoint in the QCD phase diagram.

The spectrum of coherent gluon radiation from a quark-anti-quark pair experiencing multiple scatterings within a coloured medium is central for understanding in-medium parton cascades. Despite its foundational importance, current results are limited by reliance on simplified scattering rates, such as the harmonic oscillator approximation, valid only in restricted phase-space regions. Using the formalism introduced in a previous article, we express the gluon emission spectrum as a set of differential equations that can be solved numerically, circumventing conventional approximations. We present the transverse momentum and energy distributions of emitted gluons for realistic interaction models, illustrating the breakdown of colour coherence across the entire accessible phase-space, and consequently enabling a higher-precision description of jet observables.

This work reports the first calculation of the gravitational form factors (GFFs) of the scalar glueball, performed via lattice field theory in Yang-Mills theory at a single lattice spacing. The glueball GFFs are compared with those of other hadrons as determined in previous lattice calculations, providing strong indications that glueballs have a different gluonic structure than typical hadronic states. A mass radius of 0.263(31) fm is predicted, supporting previous suggestions that the scalar glueball is significantly smaller than other hadrons. These results point towards a potential smoking-gun characteristic to target by experimental glueball searches.

We summarize and update using new NICER measurements the results of arXiv:2401.16253, in which we used various astrophysical neutron-star observations to set an upper bound on the CFL color-superconducting gap in a range of baryon chemical potentials $μ_B \in [2.1,3.2]$, above those reached within neutron stars. We also corroborate the “reasonable” constraint from arXiv:2401.16253 on the maximum value of the color-superconducting gap by performing a new Bayesian analysis using a prior that extends a two-segment Gaussian process connecting the whole density range between CEFT and pQCD.

Two-point correlation functions of systems with baryon number $B \in \{1,2,3,4\}$ are investigated using lattice Quantum Chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the $O(N)$ model. The theoretical O(N) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behaviour of the baryon number $B$ QCD correlation function at large temporal separation is well-reproduced by the $O(N \sim 2/B)$ model distribution.

The detection of compact binary mergers with sub-solar masses at gravitational-wave observatories could mark the groundbreaking discovery of primordial black holes (PBHs). Concurrently, evidence for a nHz stochastic gravitational wave background observed by pulsar timing arrays (PTAs) could suggest a non-astrophysical origin, potentially arising from scalar-induced gravitational waves (SIGW). In this work, we analyze the connection between the two phenomena in the case where they share a common origin: the collapse of large primordial curvature perturbations in the early universe. We focus on sub-solar PBH populations within reach of upcoming experiments, including the current and future runs of LIGO-Virgo-KAGRA as well as the third generation observatories such as the Einstein Telescope and Cosmic Explorer. Using a Bayesian framework with physically motivated priors, we perform a consistent model comparison that incorporates existing astrophysical bounds together with the discovery potential of future detectors. Our analysis lends stronger support for the SIGW interpretation over the astrophysical one, as the narrowed priors place greater weight on the region of highest likelihood. Ultimately, we illustrate that combining PTA data with interferometer searches can deliver correlated evidence for new physics across multiple gravitational-wave bands.

We introduce heavy-flavor energy correlators, $\langle Ψ|\mathcal{E}_{\mathcal{H}}(\vec{n}_1) \mathcal{E}_{\mathcal{H}}(\vec{n}_2)|Ψ\rangle$, as a powerful observable for profiling partonic splittings and characterizing heavy-flavor fragmentation. We present its collinear factorization, perform resummation, and demonstrate the angular distribution’s sensitivity to both the heavy-quark pair splitting and their subsequent fragmentation. We then apply the heavy-flavor EEC to probe medium-induced effects, revealing its sensitivity to medium modifications to heavy-quark pair splitting functions and to the medium’s spatial structure.

A key question about the QCD phase diagram is whether there is a critical point somewhere on the boundary between the hadronic and quark-gluon plasma phases, and if so where. Heavy-ion collisions offer a unique opportunity to search for signatures of such a critical point by analyzing event-by-event fluctuations in particle multiplicities. To draw meaningful conclusions from experimental data, a theoretical framework is needed to link QCD thermodynamics with the particle spectra and correlations observed in detectors. The Equation of State (EoS) of QCD near a critical point can be related to the universal Gibbs free energy of the 3D Ising model using four currently unknown non-universal mapping parameters whose values are determined by the microscopic details of QCD. We utilize the maximum entropy approach to freeze-out the fluctuations in order to make estimates for factorial cumulants of proton multiplicities, assuming thermal equilibrium, for a family of EoS with a 3D Ising-like critical point, varying the microscopic inputs that determine the strength and structure of the critical features. We quantify the effect of the non-universal mapping parameters, and the distance between the critical point and the freeze-out curve, on the factorial cumulants of proton multiplicities.

We study a number of 1+1d lattice models with anti-unitary symmetries that simultaneously reflect space and reverse time. Some of these symmetries are anomalous, leading to Lieb-Schultz-Mattis-type constraints, thus excluding a trivially gapped phase. Examples include a mod 8 anomaly in the Majorana chain and various mod 2 anomalies in the spin chain. In some cases, there is an exact, non-anomalous lattice symmetry that flows in the continuum to CPT. In some other cases, the CPT symmetry of the continuum theory is emergent or absent. Depending on the model, the anomaly of the lattice model is matched in the continuum in different ways. In particular, it can be mapped to an emergent anomaly of an emanant symmetry.

We derive a factorization formula for coherent and incoherent $ep$ diffraction using the soft collinear effective theory, utilizing multiple power expansion parameters to handle different kinematic regions. This goes beyond the known hard-collinear diffractive factorization to address the small-$x$ Regge dynamics and Pomeron exchange from first principles. The effective field theory analysis also uncovers and factorizes an important irreducible incoherent background generated by color-nonsinglet exchange, dubbed “quasi-diffraction”, for which we calculate the associated Sudakov suppression. For unpolarized scattering we show that there are four diffractive structure functions at leading power, and point out the importance of studying $F_{3,4}^D$ through asymmetries, in addition to $F_{2,L}^D$. For the quasi-diffractive background, we make model independent predictions for ratios of the corresponding structure functions in a perturbative kinematic region. Our analysis also makes predictions for six leading-power spin-dependent structure functions. Finally, we provide connections to diffractive parton distributions, and assess the Ingelman-Schlein model. Our work lays a path for further QCD-based studies of diffraction.

We consider the quantum creation of a universe with flat spatial sections and the topology of a 3-torus, taking into account the effect of Casimir energy. We show that the corresponding instantons are singular. Since these instantons describe universes originating in a state of infinite energy, we argue that they cannot be interpreted as quantum creation from `nothing’. If quantum corrections to the energy-momentum tensor are neglected, the spacetime of the toroidal universe reduces to de Sitter space with appropriate periodic identifications. Contrary to previous claims in the literature, this spacetime is geodesically incomplete. We argue that this spacetime describes a classical universe originating at a singularity, and not a quantum origin. We conclude that the quantum creation of a toroidal universe from nothing cannot be described in the context of semiclassical quantum gravity — it is either impossible, or it depends essentially on Planck-scale physics. We therefore see no reasonable way to compare the probability of creation of a toroidal universe, if it is possible at all, with that of a spherical universe.

We present updated constraints on the abundance of primordial black holes (PBHs) dark matter from the high-redshift Lyman-$α$ forest data from MIKE/HIRES experiments. Our analysis leverages an effective field theory (EFT) description of the 1D flux power spectrum, allowing us to analytically predict the Lyman-$α$ fluctuations on quasi-linear scales from first principles. Our EFT-based likelihood enables robust inference across redshifts $z = 4.2-5.4$ and down to scales of 100 kpc, within previously unexplored regions of parameter space for this dataset. We derive new bounds on the PBH fraction with respect to the total dark matter $f_{\text{PBH}}$, excluding populations with $f_{\text{PBH}} \gtrsim 10^{-3}$ for masses $M_{\text{PBH}} \sim 10^{4}-10^{16} M_{\odot}$. This offers the leading constraint for PBHs heavier than $10^{9} M_{\odot}$ and highlights the Lyman-$α$ forest as a uniquely sensitive probe of new physics models that modify the structure formation history of our universe.

Numerical analytic continuation arises frequently in lattice field theory, particularly in spectroscopy problems. This work shows the equivalence of common spectroscopic problems to certain classes of moment problems that have been studied thoroughly in the mathematical literature. Mathematical results due to Kovalishina enable rigorous bounds on smeared matrix-valued spectral functions, which are implemented numerically for the first time. The required input is a positive-definite matrix of Euclidean-time correlation functions; such matrices are routinely computed in variational spectrum studies using lattice quantum chromodynamics. This work connects the moment-problem perspective to recent developments using the Rayleigh–Ritz method and Lanczos algorithm. Possible limitations due to finite numerical precision are discussed.

We investigate dark matter (DM) interactions via spectroscopic signatures of energy injection in planetary environments. We develop a general framework to account for how DM energy injection signals depend on the DM spatial distribution, planetary structure, and DM energy deposition profile. We combine UV airglow data on the Solar System’s gas giants from the Voyager and New Horizons flybys, and ionospheric measurements from AMS-02 and ELFIN CubeSat on Earth, with internal heat flow data from Cassini, Voyager, and terrestrial boreholes, to constrain DM-nucleon scattering across both heavy and light mediator scenarios. We show that Earth, gas giants, and ice giants probe complementary DM masses and mediator properties, and forecast the reach of a free-floating Super-Jupiter. These results establish planetary spectroscopy as a powerful and versatile probe of the dark sector, complementary to direct detection, cosmology, and collider searches.

We present an efficient and accurate pipeline for the analysis of the redshift-space galaxy bispectrum multipoles at one-loop order in effective field theory (EFT). We provide a systematic theory derivation based on power counting, which features the first comprehensive treatment of stochastic EFT contributions — these are found to significantly improve the match to data. Our computational pipeline utilizes the COBRA technique that expands the linear matter power spectrum over a basis of principal components based on a singular value decomposition, allowing the cosmology dependence to be captured to sub-permille accuracy with just eight templates. This transforms the problem of computing the one-loop EFT bispectrum to a simple tensor multiplication, reducing the computation time to around a second per cosmology with negligible loss of accuracy. Using these tools, we study the cosmological information in the bispectrum by analyzing PTChallenge simulations, whose gigantic volume provides the most powerful test of the one-loop EFT bispectrum so far. We find that the one-loop prediction provides an excellent match to the bispectrum data up to $k_{\rm max}=0.15~h$ Mpc$^{-1}$, as evidenced by the precise recovery of the dark matter density $ω_\text{cdm}$, Hubble constant $H_0$, and mass fluctuation amplitude $σ_8$ parameters, and the amplitude of equilateral primordial non-Gaussianity (PNG) $f_{\rm NL}^{\rm equil}$. Combined with the power spectrum, the COBRA-based one-loop bispectrum multipoles yield tighter constraints than the tree-level bispectrum monopole, with the posteriors on $ω_{\text{cdm}}$, $H_0$, and $σ_8$ shrinking by 43%, 31%, and 4%, respectively. This suggests that the COBRA-based bispectrum analysis will be an important tool in the interpretation of data from ongoing redshift surveys such as DESI and Euclid.

We study the process of hydrodynamization in kinetic theories of gluons undergoing boost-invariant expansion using the Adiabatic Hydrodynamization (AH) framework. We study both number-conserving and non-conserving theories, and find that including number non-conserving inelastic scattering processes restores many qualitative features of hydrodynamization in QCD EKT despite the simplicity of our model. In particular, introducing inelastic scattering results in a more realistic hydrodynamization time. With or without number non-conservation, we find that first-order hydrodynamics becomes applicable at the same time that a unique ground state emerges in the dynamical evolution of the one-particle distribution function. Furthermore, we find a set of low-effective-energy attractor modes which evolve adiabatically long before hydrodynamization, and find that the emergence of a gap between these ground state modes and the excited modes coincides with the time at which the system falls onto an attractor surface. Strikingly, this is the case even in the absence of pre-thermal scaling of the gluon distribution function, which has previously been strongly associated with pre-thermal attractor behavior. Finally, motivated by a generic feature we observe in the spectrum, we show that as the system hydrodynamizes, the rapid decoupling of non-hydrodynamic modes in boost-invariant kinetic theory can be understood with the AH framework in a model-independent fashion.

Massive compact objects soften binaries. This process has been used for decades to constrain the population of such objects, particularly as a component of dark matter (DM). The effects of light compact objects, such as those in the unconstrained asteroid-mass range, have generally been neglected. In principle, low-energy perturbers can harden binaries instead of softening them, but the standard lore is that this effect vanishes when the perturber velocities are large compared to the binary’s orbital velocity, as is typical for DM constituents. Here, we revisit the computation of the hardening rate induced by light perturbers. We show that although the perturbations average to zero over many encounters, many scenarios of interest for DM constraints are in the regime where the variance cannot be neglected. We show that a few fast-moving perturbers can leave stochastic perturbations in systems that are measured with great precision, and we use this framework to revisit the constraint potential of systems such as binary pulsars and the Solar System. This opens a new class of dynamical probes with potential applications to asteroid-mass DM candidates.