MIT Center for Theoretical Physics — a Leinweber Institute Publications

We present the first study of jet substructure modifications during the bottom-up evolution that describes the early stages of heavy-ion collisions. To this end, we study the bremsstrahlung radiation rate of soft gluons from a hard parton propagating through out-of-equilibrium QCD matter. The gluon spectrum is computed within the Improved Opacity Expansion, which accounts for both multiple soft and single hard momentum exchanges between the hard probe and the medium. The background evolution is obtained from effective kinetic theory simulations that determine the jet quenching parameter, which in turn controls the radiation rate. We compute the radiation rate for initially under- and over-occupied systems, as well as for an expanding system undergoing hydrodynamization, which typically represents the initial stages of heavy-ion collisions. The results for these dynamical backgrounds are compared to static and thermally matched scenarios, allowing to gauge the importance of bulk expansion in the evolution of the jet cascade. Our findings show that the early stages of the bulk matter evolution in heavy-ion collisions leave a sizable imprint on the radiation pattern inside jets. These results establish a basis for incorporating pre-equilibrium dynamics into realistic descriptions of jet quenching and hard-probe evolution.

Applications of effective field theory (EFT) and scattering amplitudes to gravitational problems have recently produced many unique results that advanced our understanding of the dynamics of compact binaries. Many of these results were made possible by comparing gravitational scattering amplitudes in EFT with exact expressions from general relativity. However, the latter expressions are not easily available as they require intricate solution techniques for the Teukolsky master equation, such as the Mano-Suzuki-Takasugi (MST) method. In this paper, we develop and present the first public package that enables computations of gravitational scattering amplitudes in black hole perturbation theory within the MST framework. Our package supports both analytic computations to a given post-Minkowskian (PM) order in the low-frequency limit and exact numerical computations for an arbitrary frequency of the perturbing field. As an application, we compute scattering phase shifts and inelasticity parameters for massless spin – 0, 1, and 2 fields resulting from scattering off a rotating Kerr black hole through the third PM order and compare these results with the exact numerical solutions.

We present a first joint extraction of the Collins-Soper kernel (CSK) combining experimental and lattice QCD data in the context of an analysis of transverse-momentum-dependent distributions (TMDs). Based on a neural-network parametrization, we perform a Bayesian reweighting of an existing fits of TMDs using lattice data, as well as a joint TMD fit to lattice and experimental data. We consistently find that the inclusion of lattice information shifts the central value of the CSK by approximately 10% and reduces its uncertainty by 40-50%, highlighting the potential of lattice inputs to improve TMD extractions.

The 21-cm signal from the epoch of cosmic dawn ($z \sim 10-30$) offers a powerful probe of new physics. One standard mechanism for constraining decaying dark matter from 21-cm observations relies on heating of the intergalactic medium by the decay products, an effect whose observability is entangled with the uncertain Lyman-$α$ fluxes and X-ray heating from the first stars. In this Letter, we explore a novel mechanism, where the Lyman-$α$ photons produced from dark matter decay initiate early Wouthuysen-Field coupling of the spin temperature to the gas temperature, thereby boosting the 21-cm signal. This mechanism provides constraints on dark matter that are less dependent on uncertainties associated with star formation than constraints on exotic heating. We study this effect for decaying dark matter with masses $m_χ\sim20.4-27.2$ eV, where diphoton decay efficiently produces Lyman-series photons. We present forecasts for the Hydrogen Epoch of Reionization Array and the Square Kilometre Array, showing their potential to probe an unconstrained parameter space for light decaying DM, including axion-like particles.

We apply weakly coupled QCD kinetic theory to investigate the thermalization of high-momentum on-shell partons (minijets) in a Quark-Gluon Plasma (QGP). Our approach incorporates isotropic hard thermal loop screening to model soft quark and gluon exchanges, allowing us to verify consistency with established analytic results of jet transport coefficients. We perform kinetic simulations of minijets propagating through a thermal gluon plasma, incorporating both collinear radiation and elastic scatterings. The resulting evolution is compared to predictions from jet transport coefficients, including the longitudinal and transverse jet-quenching parameters $\hat{q}$, energy loss, and the drag coefficient. We find that standard definitions of jet transport coefficients neglect the contributions from recoiling medium particles. Including these contributions restores consistency with the kinetic evolution. Finally, we show that the minijet thermalization time scales remarkably well with $\hat{q}$ and we produce a phenomenological estimate of the minijet quenching time in heavy-ion collisions.

Heavy quarks offer an invaluable hard probe of the droplets of quark gluon plasma (QGP) formed in heavy ion collisions at the LHC and RHIC. Given their large mass, they are predominantly produced in hard scattering processes at the earliest moment of a collision and given their rarity they almost never annihilate with a heavy antiquark subsequently. This means that they experience, and probe, the entire history of the expanding, cooling, droplet of QGP from hydrodynamization through hadronization. Quantitative measurements of heavy quark final state observables therefore give us access to information about the transport properties of QGP as well as about medium modifications of hadronization. To date, the Hybrid strong/weak coupling Model of jet quenching has not included any implementation of the heavy-quark sector, which has made it impossible to confront its predictions with measurements of heavy quark and jet observables together, in a unified fashion. Here, we extend the Hybrid Model to investigate heavy quark observables for the first time. We introduce a strongly-coupled calculation of heavy-quark energy loss with the correct behavior when the heavy quarks are either ultrarelativistic or non-relativistic, Gaussian momentum broadening, and recombination of heavy quarks with medium partons using a local color neutralization model of hadronization. We compare our results for the suppression $R_{\rm AA}$ and azimuthal anisotropies $v_2$ of B- and D-mesons and $Λ_c$ baryons, the $R_{\rm AA}$ of B-tagged jets, as well as baryon-to-meson ratios, with available experimental data from ALICE, ATLAS and CMS.

We investigate quantum-information-theoretic observables in electron-proton scattering at the Electron-Ion Collider (EIC). Our analysis focuses on entanglement and magic, two complementary indicators of non-classicality in quantum states. We show that while unpolarized and longitudinally polarized beams yield unentangled separable outcomes, transverse beam polarization enables the generation of entangled and non-stabilizer states. This result holds for both elastic and deep inelastic electron-proton scattering in QED. In the deep inelastic regime, the degree of quantum correlation is governed by the transversity parton distribution functions, providing a novel perspective on spin dynamics within QCD. These results establish the EIC as a promising environment for generating entangled and non-stabilizer states in high-energy physics, and they highlight opportunities for future lepton-hadron colliders to extend such studies into new kinematic domains.

Cosmological studies of the Lyman-Alpha (Lya) forest typically constrain parameters using two-point statistics. However, higher-order statistics, such as the three-point function (or its Fourier counterpart, the bispectrum) offer additional information and help break the degeneracy between the mean flux and power spectrum amplitude, albeit at a significant computational cost. To address this, we extend an existing highly informative compression of the bispectrum, the skew spectra, to the Lya forest. We derive the tree-level bispectrum of Lya forest fluctuations in the framework of effective field theory (EFT) directly in redshift space and validate our methodology on synthetic Lya forest data. We measure the anisotropic cross-spectra between the transmitted flux fraction and all quadratic operators arising in the bispectrum, yielding a set of 26 skew spectra. Using idealized 3D Gaussian smoothing (R=10 Mpc/h), we find good agreement (1-2 sigma level based on the statistical errors of the mocks) with the theoretical tree-level bispectrum prediction for monopole and quadrupole up to k <= 0.17 h/Mpc. To enable the cosmological analysis of Lya forest data from the currently observing Dark Energy Spectroscopic Instrument (DESI), where we cannot do 3D smoothing, we use a line-of-sight smoothing and introduce a new statistic, the shifted skew spectra. These probe non-squeezed bispectrum triangles and avoid locally applying quadratic operators to the field by displacing one copy of the field in the radial direction. Using a fixed displacement of 40 Mpc/h (and line-of-sight smoothing of 10 Mpc/h) yields a similar agreement with the theory prediction. For the special case of correlating the squared (and displaced) field with the original one, we analytically forward model the window function making this approach readily applicable to DESI data.

Light scalar particles arise naturally in many extensions of the Standard Model and are well-motivated dark matter candidates. Gravitational interactions near black holes can trigger the growth of dense scalar configurations that, if sustained during inspiral, alter binary dynamics and imprint signatures on gravitational-wave signals. Detecting such effects would provide a novel probe of fundamental physics and dark matter. Here we develop a semi-analytic waveform model for binaries in scalar environments, validated against numerical relativity simulations, and apply it in a Bayesian analysis of the LIGO-Virgo-KAGRA catalog. Our results set physically meaningful upper bounds on scalar environments around compact binaries. When superradiance priors are included, we find tentative evidence for such an environment in GW190728 with $\ln B_{\mathrm{vac}}^{\mathrm{env}} \approx 3.5$, which would correspond to the existence of a light scalar field with mass $\sim 10^{-12}\,\mathrm{eV}$.

The observation of hybrid mesons in photoproduction experiments can provide essential insight into the inner workings of quantum chromodynamics in the strong coupling regime. In particular, the study of final $η^{(\prime)}π$ states is of great interest due to the presence of the lowest lying hybrid candidate with manifestly exotic quantum numbers, the $π_1(1600)$. In this work, a double-vector exchange model with Reggeized $ρ$ and $ω$ trajectories is developed to describe the photoproduction of $η^{(\prime)}π$ in the high-mass region. Results are presented for the differential cross sections and forward-backward asymmetries in the energy region of interest to the GlueX experiment. The model contains no free parameters, and reproduces the magnitude and $t$-dependence of existing CLAS data at $E_γ=5$\gev. The model predicts a stronger asymmetry in the $η^\prime π$ channel than in the $ηπ$ channel, consistent with what has previously been observed in pion beam experiments. This suggests that the sizeable production of exotic odd waves in $η’π$ is not necessarily related to the presence of gluon-rich environments. Confirmation of these predicted asymmetries from forthcoming GlueX data would enable further predictions of the low-energy spectrum.

We study heat conduction and momentum transport in the context of stochastic fluid dynamics. We consider a fluid described by model H in the classification of Hohenberg and Halperin. We study both non-critical and critical fluids, and we investigate transport properties in two as well as three dimensions. Our results are based on numerical simulations of model H using a Metropolis algorithm, and we employ Kubo relations to extract transport coefficients. We observe the expected logarithmic divergence of the shear viscosity in a two-dimensional non-critical fluid. At a critical point, we find that the transport coefficients exhibit power-law scaling with the system size $L$. The strongest divergence is seen for the thermal conductivity $κ$ in two dimensions. We find $κ\sim L^{x_κ}$ with $x_κ=1.6\pm 0.1$. The divergence is weaker in three dimensions, $x_κ=1.25 \pm 0.3$, and the scaling exponent for the shear viscosity, $x_η$, is significantly smaller than $x_κ$ in both two and three dimensions.

Gauge theories describe the fundamental forces of nature. However, high-energy dynamics, such as the formation of quark-gluon plasmas, is notoriously difficult to model with classical methods. Quantum simulation offers a promising alternative in this regime, yet experiments have mainly probed low energies. Here, we observe the formation of a ballistic plasma and long-time memory effects in high-energy gauge theory dynamics on a high-precision quantum simulator. Both observations are unexpected, as the initial state – fully filled with particle-antiparticle pairs – was thought to rapidly thermalize. Instead, we find correlations spreading ballistically to long distances and a memory of charge clusters. Our observations cannot be explained by many-body scars, but are captured by a new theory of plasma oscillations between electric field and current operators, persisting all the way to the continuum limit of the (1+1)D Schwinger model, of which we simulate a lattice version. Adapting techniques from quantum optics, we visualize plasma oscillations as rotations of Wigner distributions, leading to a novel set of predictions which we test in experiment and numerics. The new framework encompasses both our scenario and scars, which show up as coherent states of the plasma. The experimental surprises we observe in the high-energy dynamics of a simple gauge theory point to the potential of high-precision quantum simulations of gauge theories for general scientific discovery.

Benchmarking quantum devices is a foundational task for the sustained development of quantum technologies. However, accurate in situ characterization of large-scale quantum devices remains a formidable challenge: such systems experience many different sources of errors, and cannot be simulated on classical computers. Here, we introduce new benchmarking methods based on random circuit sampling (RCS), that substantially extend the scope of conventional approaches. Unlike existing benchmarks that report only a single quantity–the circuit fidelity–our framework extracts rich diagnostic information, including spatiotemporal error profiles, correlated and contextual errors, and biased readout errors, without requiring any modifications of the experiment. Furthermore, we develop techniques that achieve this task without classically intractable simulations of the quantum circuit, by leveraging side information, in the form of bitstring samples obtained from reference quantum devices. Our approach is based on advanced high-dimensional statistical modeling of RCS data. We sharply characterize the information-theoretic limits of error estimation, deriving matching upper and lower bounds on the sample complexity across all regimes of side information. We identify surprising phase transitions in learnability as the amount of side information varies. We demonstrate our methods using publicly available RCS data from a state-of-the-art superconducting processor, obtaining in situ characterizations that are qualitatively consistent yet quantitatively distinct from component-level calibrations. Our results establish both practical benchmarking protocols for current and future quantum computers and fundamental information-theoretic limits on how much can be learned from RCS data.

A fundamental task in quantum information is to approximate a pure quantum state in terms of sparse states or, for a bipartite system, states of bounded Schmidt rank. The optimal deterministic approximation in each case is straightforward, and maximizes the fidelity: keep the largest entries or singular values. On the other hand, random mixtures of sparse states can achieve quadratically improved trace distances, and yield nontrivial bounds on other distance measures like the robustness. In this work, we give efficient algorithms for finding mixtures of sparse states that optimally approximate a given pure state in either trace distance or robustness. These algorithms also yield descriptions of efficiently samplable ensembles of sparse, or less-entangled, states that correspond to these optimal mixed approximations. This can be used for the truncation step of algorithms for matrix product states, improving their accuracy while using no extra memory, and we demonstrate this improvement numerically. Our proofs use basic facts about convex optimization and zero-sum games, as well as rigorous guarantees for computing maximum-entropy distributions.

In this Letter, we initiate a systematic study of the $n$-point correlation functions (CF) in gauge theories in the sequential light-cone (SLC) limit. Focusing on QCD, we formulate a factorization theorem for the CF of four vector currents in this limit using tools from soft-collinear effective field theory (SCET). This result unveils a duality between CF and Wilson loops in non-conformal field theories, according to which the singular structure of the CF is described by a null polygonal Wilson loop dressed by universal jet functions and current form factors, directly related to well-known ingredients in collider physics. We employ this new factorization theorem to obtain the singular terms of the four-point CF in QCD in the SLC limit up to three loops. This constitutes the first determination of terms beyond one loop and in particular determines, for the first time, conformal-symmetry-breaking terms induced by the non-vanishing $β$-function of QCD. As a stringent check, we verify that our result satisfies the anomalous conformal Ward identities and its leading-transcendental term reproduces the known correlation function in $\mathcal{N}=4$ SYM to three loops. Our results are easily generalized to the multi-point case, offering a powerful tool for future applications of SCET to the study of the fundamental structure of gauge theories.

We review recent developments in the use of von Neumann algebras to analyze the entanglement structure of quantum gravity and the emergence of spacetime in the semi-classical limit. Von Neumann algebras provide a natural framework for describing quantum subsystems when standard tensor factorizations are unavailable, capturing both kinematic and dynamical aspects of entanglement. The first part of the review introduces the fundamentals of von Neumann algebras, including their classification, and explains how they can be applied to characterize entanglement. Topics covered include modular and half-sided modular flows and their role in the emergence of time, as well as the crossed-product construction of von Neumann algebras. The second part turns to applications in quantum gravity, including an algebraic formulation of AdS/CFT in the large-$N$ limit, the emergence of bulk spacetime structure through subregion-subalgebra duality, and an operator-algebraic perspective on gravitational entropy. It also discusses simple operator-algebraic models of quantum gravity, which provide concrete settings in which to explore these ideas. In addition, several original conceptual contributions are presented, including a diagnostic of firewalls and an algebraic formulation of entanglement islands. The review concludes with some speculative remarks on the mathematical structures underlying quantum gravity.

We report on our investigation in arXiv:2509.08881 of how recent jet substructure measurements constrain the resolution length $L_{\rm res}$ of the quark-gluon plasma formed in heavy-ion collisions. $L_{\rm res}$ is defined such that high-energy partons within a jet shower are resolved by the medium if and only if they are separated by a distance greater than $L_{\rm res}$. Using the Hybrid Model, we reproduce ALICE data on the scaled Soft Drop angle $θ_g$ for $R=0.2$ charged-particle jets and ATLAS data on the Hard Group angle $ΔR_{12}$ for $R=1$ jets reclustered from skinny $R = 0.2$ inclusive subjets. We find that the narrowing of the $θ_g$-distribution in PbPb collisions observed by ALICE and the suppression of $R=1$ jets with multiple skinny subjets in PbPb collisions observed by ATLAS rule out $L_{\rm res} = \infty$, where each entire parton shower loses energy to the plasma coherently as if it were a single colored object. We then compare Hybrid Model calculations to ATLAS measurements of $R_{\rm AA}$ of $R=1$ jets reclustered from $R=0.2$ subjets, as a function of the Soft Drop angle $dR_{12}$ obtained by grooming all charged-particle tracks associated with each $R=1$ jet. We demonstrate, for the first time, that the ATLAS data is inconsistent with $L_{\rm res} = 0$, where the plasma resolves every splitting in a parton shower. Our results agree best with the data when QGP possesses a finite, nonzero $L_{\rm res}\sim (1-2)/(πT)$.

The effective kinetic theory (EKT) of QCD provides a possible picture of various non-equilibrium processes in heavy- and light-ion collisions. While there have been substantial advances in simulating the EKT in simple systems with enhanced symmetry, eventually, event-by-event simulations will be required for a comprehensive phenomenological modeling. As of now, these simulations are prohibitively expensive due to the numerical complexity of the Monte Carlo evaluation of the collision kernels. In this talk, we show how the evaluation of the collision kernels can be performed using neural networks paving the way to full event-by-event simulations.

We obtain the collision kernel and related dipole cross section during the initial nonequilibrium stages in heavy-ion collisions. These quantities are a crucial input for jet quenching calculations. We further compute the gluon splitting rates in the AMY formalism resulting from this nonequilibrium kernel. Comparing with thermal and commonly used forms, we find that particularly the gluon splitting rate for parton energies of the order of the hard effective temperature significantly differs from these approximations.

In this paper, we present a framework for the analytic bootstrap of three-point energy correlators, a crucial observable in $\mathcal{N}=4$ super Yang-Mills theory and quantum chromodynamics (QCD). Our approach combines spherical contour techniques, general physical constraints such as pole cancellations, and power correction data in the singular limits to determine its analytic expression. In contrast to previous bootstrap studies restricted to scattering amplitudes for supersymmetric theories, our framework makes use of the properties of Feynman integrals, marking a significant step toward bootstrapping realistic QCD observables. Using this method, we derive analytic expressions for leading-order three-point energy correlators with equal and unequal energy weights, where the latter are crucial ingredients for projected $N$-point energy correlators. We also apply the recently developed technique of analytic regression with lattice reduction as a way to bypass needing explicit expressions for the singular limits. Bridging theoretical advances in scattering amplitudes with the renewed interest in weighted cross-sections, our work opens the door to precision tests of QCD dynamics through analytic event-shape predictions.