MIT Center for Theoretical Physics — a Leinweber Institute Publications

We consider the energy loss of a high-energy parton in the early anisotropic plasma in heavy-ion collisions. Working in the harmonic approximation, we compute the change in the mean energy of an emitted gluon in the presence of an anisotropic background, characterized by anisotropic jet quenching parameters $\hat q_{x}\neq \hat q_{y}$. Evaluating the resulting integrals numerically, we compare with isotropic media, and obtain a simple pocket formula to estimate the impact of anisotropy on the mean emitted gluon energy, which is generally small. We then combine our results with the values of the jet quenching parameter extracted from QCD kinetic theory simulations and show that the medium length dependence of this mean energy loss exhibits the characteristics of limiting attractors, which can be obtained by extrapolating to zero and infinite coupling. Our study thus relates energy loss of jet partons to universal dynamics of anisotropic plasmas.

Quark-gluon plasma (QGP) is a strongly coupled liquid when viewed at length scales of order the inverse of its temperature and longer. However, when it is probed at short enough length scales, asymptotic freedom mandates the presence of quark- and gluon-like quasiparticles. Partons in jets can trigger perturbative, high momentum-exchange $2\rightarrow2$ Molière scatterings off quasiparticles in the medium, making jets useful probes of the microscopic structure of QGP. Prior to this work, soft strongly coupled momentum-exchanges between jet partons and the QGP droplet produced in a heavy-ion collision, as well as the wakes that jets excite in the droplet, had been accounted for in the Hybrid Model of jet quenching. Here, we present a full calculation of Molière scattering off a QGP quasiparticle which results in the deflection of the jet parton and the excitation of a parton from the thermal medium that recoils after being kicked, and describe how it is implemented in the Hybrid Model. The scattered jet and recoil partons continue to propagate through the QGP, lose energy and momentum, excite wakes, and may further re-scatter. Using the Hybrid Model, we study how Molière scatterings impact jet shapes and fragmentation functions, the Soft Drop angle $R_g$, jet girth $g$, and observables that focus on the number and angular distribution of subjets within jets. We demonstrate that photon-tagged jets provide a particularly sensitive probe: selecting events by the photon energy mitigates the selection bias inherent in inclusive jet measurements and enhances sensitivity to rare large-angle scatterings. We find that Molière scatterings broaden both the $R_g$ and $g$ distributions when jets significantly softer than the photon are included. Our results point the way towards distinctive model-independent experimental signatures of hard scattering of jet partons off quasiparticles in QGP.

We present a combined analytical and numerical study of the effective action of field theory cosmic strings in the Abelian-Higgs model in flat space. Starting directly from the underlying solitonic field theory description, we provide a systematic derivation of the low energy effective action and present evidence for the absence of nontrivial curvature correction terms when only the translational Goldstone modes are retained. Using this framework, we extend the effective theory to include higher energy fluctuations of the soliton profile, which map to massive degrees of freedom propagating on the worldsheet. We show that the leading curvature contribution enters only through the coupling between these massive modes and the worldsheet Ricci scalar. We validate the resulting effective theory via lattice simulations of the full field theory equations of motion in flat space, implemented with Adaptive Mesh Refinement to capture the string dynamics across different scales. The numerical simulations confirm the dynamics obtained using the effective action in its validity range. Furthermore, they also demonstrate the existence of the predicted parametric instability of excited strings that drives the transfer of energy from massive excitations to the Goldstone sector.

Developing measures of quantum ergodicity and chaos stands as a foundational task in the study of quantum many-body systems. In this work, we propose metrics for these effects based on Hamiltonian learning that unify multiple advantages of existing metrics. In particular, we show how ergodicity and chaos improve the robustness of Hamiltonian learning to small errors and furthermore demonstrate that this robustness can be used as a metric for such phenomena. We analytically and numerically show that our metrics not only distinguish between integrable and ergodic regimes in various spin chains but also quantify chaos and ergodicity, allowing us to locate regions of parameter space displaying maximal ergodicity and maximal sensitivity to local perturbations. Our approach not only provides conceptual ways to study quantum chaos and ergodicity but also presents viable experimental methods for quantum simulators.

The light-cone distribution amplitude (LCDA) is a fundamental non-perturbative quantity for understanding hadron structure and exclusive scattering processes. We report on our calculation of the pion and kaon LCDAs using the heavy-quark operator product expansion (HOPE) framework. This method employs an OPE analysis of hadronic amplitudes through the inclusion of a fictitious valence heavy quark. In these proceedings, we report progress on the determination of the first three nontrivial Mellin moments of the kaon LCDAs from dynamical lattice QCD calculations, and we summarize the recently published continuum-limit result for the pion fourth Mellin moment obtained in the quenched approximation, thereby demonstrating the feasibility of the HOPE method for accessing higher moments.

Normalizing flows can be used to construct unbiased, reduced-variance estimators for lattice field theory observables that are defined by a derivative with respect to action parameters. This work implements the approach for observables involving gluonic operator insertions in the SU(3) Yang-Mills theory and two-flavor Quantum Chromodynamics (QCD) in four space-time dimensions. Variance reduction by factors of $10$-$60$ is achieved in glueball correlation functions and in gluonic matrix elements related to hadron structure, with demonstrated computational advantages. The observed variance reduction is found to be approximately independent of the lattice volume, so that volume transfer can be utilized to minimize training costs.

We analyze the public DESI full-shape clustering data using simulation-based priors (SBPs). Our priors are obtained by fitting normalizing flows to the distribution of EFT parameters measured from field-level simulations, themselves generated using tailored halo occupation distribution (HOD) models for each tracer. Incorporating SBPs in a power spectrum analysis significantly enhances $Λ$CDM cosmological parameter constraints; in combination with BAO information from DESI DR2 and a BBN prior on the baryon density, we find the matter density parameter $Ω_m=0.2987\pm 0.0066$, the Hubble constant $H_0=68.80\pm 0.35\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, and the mass fluctuation amplitude $σ_8 = 0.766\pm 0.015$ (or the lensing parameter $S_8=0.764\pm 0.018$), which are $1\%$, $40\%$ and $50\%$ stronger than the baseline results, though with a notable downwards shift in $σ_8$. The SBPs also have a significant impact in extended models, with the dark energy figure-of-merit improving by $70\%$ ($20\%$) in a $w_0w_a$CDM analysis when combining with the CMB (and supernovae). In the SBP analysis, we do not find statistically significant evidence for dynamical dark energy: the equation of state parameters are consistent with a cosmological constant within $2.2σ$ ($1.4σ$) in analyses without (with) supernovae. The neutrino mass constraints are also enhanced, with the $95\%$ limits $M_ν<0.073\,\mathrm{eV}$ and $M_ν<0.090\,\mathrm{eV}$ in $Λ$CDM and $w_0w_a$CDM respectively. The latter is the strongest constraint obtained to date and reinforces the preference for the normal neutrino mass hierarchy, regardless of the background dynamics. While our results are sensitive to HOD modeling assumptions, they clearly demonstrate that the inclusion of small-scale information can significantly sharpen cosmological parameter constraints.

Sterile neutrinos in the keV mass range have long been considered a well-motivated dark matter (DM) candidate. In this work, we explore a sterile neutrino production mechanism through active neutrino self-interactions in the early universe, assuming that they constitute the full DM abundance. We implement a self-consistent treatment of the sterile-neutrino free streaming and the active-neutrino self-interactions on structure formation, which yield a unique scale-dependent modification to the linear matter power spectrum. We then set bounds on this scenario using a combination of the cosmic microwave background and Ly$α$ forest constraints. Specifically, we utilize the two recent likelihoods derived from eBOSS data: (i) an effective field theory (EFT) based full-shape likelihood and (ii) a compressed likelihood obtained from the PRIYA-simulation emulator. We produce some of the most stringent observational constraints to date on sterile neutrino DM, comparable to the bounds from the most stringent laboratory constraints.

Modern precision experiments often probe unknown classical fields with bosonic sensors in quantum-noise-limited regimes where vacuum fluctuations limit conventional readout. We introduce Quantum Signal Learning (QSL), a sensing framework that extends metrology to a broader property-learning setting, and propose a quantum-enhanced protocol that simultaneously estimates many properties of a classical signal with shot noise suppressed below the vacuum level. Our scheme requires only two-mode squeezing, passive optics, and static homodyne measurements, and enables post-hoc classical estimation of many properties from the same experimental dataset. We prove that our protocol enables a quantum speedup for common classical sensing tasks, including measuring electromagnetic correlations, real-time feedback control of interferometric cavities, and Fourier-domain matched filtering. To establish these separations, we introduce an optimal-transport conditioning method, and show both worst-case exponential separations from all entanglement-free strategies and practical speedups over homodyne and heterodyne baselines. We further show that when squeezing is treated as a resource, a protocol with squeezed light can sense a structured classical background exponentially faster than any coherent classical probe.

Entanglement is the hallmark of quantum physics, yet its characterization in interacting many-body systems at thermal equilibrium remains one of the most important challenges in quantum statistical physics. We prove that the Gibbs state of any quantum spin chain can be exactly decomposed into a mixture of matrix product states with a bond dimension that is independent of the system size, at any finite temperature. As a consequence, the Schmidt number, arguably the most stringent measure of bipartite entanglement, is strictly finite for thermal states, even in the thermodynamic limit. Our decomposition is explicit and is accompanied by an efficient classical algorithm to sample the resulting matrix product states.

Dark sectors consisting of atomic constituents (electrons, protons, and photons) offer a well-motivated extension to the Standard Model while providing multiple avenues for phenomenological study. In this work, we explore the impact of conversions between the dark and Standard Model photons in the primordial CMB spectral distortion epoch ($10^3 \lesssim z \lesssim 10^6$). These conversions are resonantly enhanced when the induced thermal masses of both photonic species are equal, thus leading to the possibility that sizeable distortions can be produced. To this end, we solve the Boltzmann equation at early times to determine the (irreducible) freeze-in or freeze-out abundance of dark photons. This procedure also allows us to update the limits on generic milli-charged dark sectors using the ACT DR6 bound on the number of effective radiative degrees of freedom ($N_{\rm eff}$). By then modeling the evolution of the thermal masses in both sectors, we compute the primordial CMB distortion using the Landau-Zener formalism. We find that when the dark electron and proton are roughly similar in mass (the positronium limit), current spectral distortion data from the COBE/FIRAS instrument is able to rule out novel regions of parameter space. We also forecast bounds from the proposed FOSSIL satellite, finding that spectral distortions can also be used to probe the ultra-low dark electric charge regions of parameter space, which are difficult to investigate by other means.

Recently, Jones et al. [arXiv:2310.12859] claimed strong evidence for the statistical anisotropy of the universe. The claim is based on a joint analysis of four different anomaly tests of the cosmic microwave background data, each of which is known to be anomalous, with a lower level of significance. They reported a combined $p$-value of about $3\times 10^{-8}$, which is more than a $5σ$ level of significance. We observe that statistical anisotropy is not even relevant for two of the four considered tests, which seems sufficient to invalidate the authors’ claim. Furthermore, even if one reinterprets the claim as evidence against $Λ$CDM rather than statistical anisotropy, we argue that this result significantly suffers from the look-elsewhere effect. Assuming a set of independent (i.e., uncorrelated) tests, we show that if the four tests with the smallest $p$-values are cherry-picked from 10 independent tests, the $p$-value reported by Jones et al. corresponds to only $3σ$ significance. If there are 27 independent tests, the significance falls to $2σ$. These numbers, however, overstate our argument, since the four tests used by Jones et al. are slightly correlated. Determining the correlation of Jones et al.’s tests by comparing their joint $p$-value with the product of the four separate $p$-values, we find that about 16 or 50 tests are sufficient to reduce the significance of Jones et al.’s results to 3$σ$ or 2$σ$ significance, respectively. We also provide a list of anomaly tests discussed in the literature (and propose a few generalizations), suggesting that very plausibly 16 (or even 50) independent tests have been published, and possibly many more have been considered but not published. We conclude that the current data is consistent with the $Λ$CDM model and, in particular, with statistical isotropy.

We present a framework for systematic computations of scattering amplitudes for gravitational Raman scattering, — the inelastic scattering of massless fields off compact relativistic objects. We focus on the small-frequency (post-Minkowskian, PM) regime relevant for the study of tidal effects, which can be mapped onto gravitational wave observables during the inspiraling phase of a merger. We demonstrate that this setup is ideal for systematic studies of tidal effects, in a way that is free from coordinate, gauge, and field redefinition ambiguities. We use a combination of worldline effective field theory, the background field method, and advanced scattering amplitude techniques to derive phase shifts for scattering of spin-$0,1,2$ fields off generic compact objects at third PM order. We demonstrate that the inclusion of the recoil of the object is crucial for consistency of this calculation. Focusing on a particular case of black holes, we extract the leading static and dynamical Love numbers of the spin-0 field and the static Love number of the spin-1 field in four dimensions by matching our EFT amplitudes and calculations in General Relativity. We show, fully on-shell, that the leading static Love numbers vanish identically, while the dynamical Love numbers are not zero and run logarithmically. The latter resolves the ambiguities of previous off-shell matching calculations. We also extend our results to seven dimensions, where spin-2 Love numbers undergo a renormalization group running at 2PM, which we compute explicitly. In addition, we extract the leading static Love numbers of spin-0 and spin-1 fields in five dimensions, which also run.

The one-state statement for closed universes has sparked considerable discussion. In this paper we examine its physical meaning in the context of the Hartle-Hawking state and de Sitter space. We argue that the one-state property of closed universes is fully compatible with the finite-dimensional quantum mechanics experienced by observers inside de Sitter space, and that this compatibility requires neither mixing of alpha sectors nor any modification of the rules of the gravitational path integral. The apparent tension is resolved by sharply distinguishing the baby-universe Hilbert space, namely the space of closed universes viewed from the outside, from the bulk Hilbert space that governs quantum mechanics for an observer inside a single de Sitter universe. The baby-universe Hilbert space, together with its commutative operator algebra, is not a quantum-mechanical Hilbert space: it is merely a mathematical repackaging of classical probability theory and does not carry any quantum-mechanical structure at all, a direct consequence of the one-state property of closed universes. Accordingly, attempting to formulate quantum mechanics directly on the baby-universe Hilbert space conflates classical ensemble data with the quantum mechanics experienced by bulk observers and leads to physically incorrect conclusions. By contrast, the quantum mechanics experienced by an observer inside de Sitter space emerges from the classical statistics encoded in the baby-universe Hilbert space, providing a concrete realization of Wheeler’s idea of “It from Bit”. We demonstrate these features by completely solving a topological toy model of one-dimensional de Sitter spacetime. Along the way we clarify the physical meaning of de Sitter entropy, showing that it corresponds to the coarse-grained entropy of the underlying state.

In this paper we extend recent ideas about observers and closed universes to theories where observers can be fluctuated into existence in the Hartle-Hawking state. This introduces a phenomenon that was not considered in these earlier discussions: the dominant transition from one cosmological state to another can go through a fluctuation that annihilates the universe and creates a new one. We nonetheless argue that the observer decoherence rule allows for the third-quantized description of such a theory to emerge from a factorizing holographic theory with a one-dimensional Hilbert space, without any need for $α$-parameters. We also point out a close analogy between the observer rule in this context and the coarse-graining of the spectral form factor at late times for AdS black holes. Along the way we clarify several aspects of the relationship between holography, the gravitational path integral, and $α$-parameters. We also explain why string theory scattering amplitudes do not lead to a one-dimensional Hilbert space on the worldsheet, despite being computed by a gravitational path integral with a sum over topology. Finally we point out that using the path integral to compute integrated local operators conditioned on an observer in the context of a theory with a landscape can lead to rather surprising conclusions. For example we argue that in a landscape with one AdS minimum and one dS minimum, both of which can support observers, an observer almost surely finds themself in dS and not AdS even if the boundary conditions are dual to a state with an observer in AdS.

[Abridged] Some cosmological theories propose that the observable universe is a small part of a much larger universe in which parameters describing the low-energy laws of physics vary from region to region. How can we reasonably assess a theory that describes such a mostly unobservable universe? We propose a Bayesian method based on theory-generated probability distributions for our observations. We focus on basic principles, leaving aside concerns about practicality. (We also leave aside the measure problem, to discuss other issues.) We argue that cosmological theories can be tested by standard Bayesian updating, but we need to use theoretical predictions for “first-person” probabilities — i.e., probabilities for our observations, accounting for all relevant selection effects. These selection effects can depend on the observer, and on time, so in principle first-person probabilities are defined for each observer-instant — an observer at an instant of time. First-person probabilities should take into account everything the observer believes about herself and her surroundings — i.e., her “subjective state”. We advocate a “Principle of Self-Locating Indifference” (PSLI), asserting that any real observer should make predictions as if she were chosen randomly from the theoretically predicted observer-instants that share her subjective state. We believe the PSLI is intuitively very reasonable, but also argue that it maximizes the expected fraction of observers who will make correct predictions. Cosmological theories will in general predict a set of possible universes, each with a probability. To calculate first-person probabilities, we argue that each possible universe should be weighted by the number of observer-instants in the specified subjective state that it contains. We also discuss Boltzmann brains, the humans/Jovians parable of Hartle and Srednicki, and the use of “old evidence”.

Wilson loops are essential objects in QCD and have been pivotal in scale setting and demonstrating confinement. Various generalizations are crucial for computations needed in effective field theories. In lattice gauge theory, Wilson loop calculations face challenges, including excited-state contamination at short times and the signal-to-noise ratio issue at longer times. To address these problems, we develop a new method by using neural networks to parametrize interpolators for the static quark-antiquark pair. We construct gauge-equivariant layers for the network and train it to find the ground state of the system. The trained network itself is then treated as our new observable for the inference. Our results demonstrate a significant improvement in the signal compared to traditional Wilson loops, performing as well as Coulomb-gauge Wilson-line correlators while maintaining gauge invariance. Additionally, we present an example where the optimized ground state is used to measure the static force directly, as well as another example combining this method with the multilevel algorithm. Finally, we extend the formalism to find excited-state interpolators for static quark-antiquark systems. To our knowledge, this work is the first study of neural networks with a physically motivated loss function for Wilson loops.

We propose and analyze a versatile and efficient multiparameter quantum sensing protocol, which simultaneously estimates many non-commuting and time-dependent signals that are coherently or incoherently coupled to sensing particles. Even in the presence of control imperfections and readout errors, our approach can detect exponentially many parameters in the system size while maintaining the optimal scaling of sensitivity. To accomplish this, scrambling dynamics are leveraged to map distinct signals to unique patterns of bitstring measurements, which distinguishes a large number of signals without significant sensitivity loss. Based on this principle, we develop a computationally efficient protocol utilizing random global Clifford unitaries and evaluate its performance both analytically and numerically. Our protocol naturally extends to scrambling dynamics generated by random local Clifford circuits, local random unitary circuits (RUCs), and ergodic Hamiltonian evolution–commonly realized in near-term quantum hardware–and opens the door to applications ranging from precise noise benchmarking of quantum dynamics to learning time-dependent Hamiltonians.

Excited-state effects lead to hard-to-quantify systematic uncertainties in lattice quantum chromodynamics (LQCD) spectroscopy calculations when computationally accessible imaginary times are smaller than inverse excitation gaps, as often arises for multi-hadron systems with signal-to-noise problems. Lanczos residual bounds address this by providing two-sided constraints on energies that do not require assumptions beyond Hermiticity, but often give very conservative systematic uncertainty estimates. Here, a more-constraining set of gap bounds is introduced for hadron spectroscopy. These bounds provide tighter constraints whose validity requires an explicit assumption about an energy gap. Exactly solvable lattice field theory correlators are used to test the utility of residual and gap bounds at finite and infinite statistics. Two-sided bounds and other analysis methods are then applied to a high-statistics LQCD calculation of nucleon-nucleon scattering at $m_π\sim 800$ MeV. Generalized eigenvalue problem (GEVP) and Lanczos energy estimators are compatible when applied to the same correlator data, but analyses including different interpolating operators show statistically significant inconsistencies. However, two-sided bounds from all operators are consistent. Under the assumption that the number of energy levels below $NΔ$ and $ΔΔ$ thresholds is the same as for non-interacting nucleons, gap bounds are sufficient to constrain nucleon-nucleon scattering amplitudes at phenomenologically relevant precision. Lanczos methods further reveal that energy-eigenstate estimates from previously studied asymmetric correlators have not converged over accessible imaginary times. Nevertheless, data-driven examples demonstrate why assumptions are required to draw conclusions about the natures of two-nucleon ground states at these masses.

We study the lattice version of higher-form symmetries on tensor-product Hilbert spaces. Interestingly, at low energies, these symmetries may not flow to the topological higher-form symmetries familiar from relativistic quantum field theories, but instead to non-topological higher-form symmetries. We present concrete lattice models exhibiting this phenomenon. One particular model is an $\mathbb{R}$ generalization of the Kitaev honeycomb model featuring an $\mathbb{R}$ lattice 1-form symmetry. We show that its low-energy effective field theory is a gapless, non-relativistic theory with a non-topological $\mathbb{R}$ 1-form symmetry. In both the lattice model and the effective field theory, we demonstrate that the non-topological $\mathbb{R}$ 1-form symmetry is not robust against local perturbations. In contrast, we also study various modifications of the toric code and their low-energy effective field theories to demonstrate that the compact $\mathbb{Z}_2$ lattice 1-form symmetry does become topological at low energies unless the Hamiltonian is fine-tuned. Along the way, we clarify the rules for constructing low-energy effective field theories in the presence of multiple superselection sectors. Finally, we argue on general grounds that non-compact higher-form symmetries (such as $\mathbb{R}$ and $\mathbb{Z}$ 1-form symmetries) in lattice systems generically remain non-topological at low energies, whereas compact higher-form symmetries (such as $\mathbb{Z}_{n}$ and $U(1)$ 1-form symmetries) generically become topological.