MIT Center for Theoretical Physics — a Leinweber Institute Publications

In the $1+1\mathrm{D}$ ultralocal lattice Hamiltonian for staggered fermions with a finite-dimensional Hilbert space, there are two conserved, integer-valued charges that flow in the continuum limit to the vector and axial charges of a massless Dirac fermion with a perturbative anomaly. Each of the two lattice charges generates an ordinary U(1) global symmetry that acts locally on operators and can be gauged individually. Interestingly, they do not commute on a finite lattice and generate the Onsager algebra, but their commutator goes to zero in the continuum limit. The chiral anomaly is matched by this non-Abelian algebra, which is consistent with the Nielsen-Ninomiya theorem. We further prove that the presence of these two conserved lattice charges forces the low-energy phase to be gapless, reminiscent of the consequence from perturbative anomalies of continuous global symmetries in continuum field theory. Upon bosonization, these two charges lead to two exact U(1) symmetries in the XX model that flow to the momentum and winding symmetries in the free boson conformal field theory.

We explore full-shape analysis with simulation-based priors, which is the simplest approach to galaxy clustering data analysis that combines effective field theory (EFT) on large scales and numerical simulations on small scales. The core ingredient of our approach is the prior density of EFT parameters which we extract from a suite of 10,500 galaxy simulations based on the halo occupation distribution (HOD) model. We measure the EFT parameters with the field-level forward model, which enables us to cancel cosmic variance. On the theory side, we develop a new efficient approach to calculate field-level transfer functions using time-sliced perturbation theory and the logarithmic fast Fourier transform. We find that the cosmology dependence of EFT parameters of galaxies is approximately degenerate with the HOD parameters, and hence it can be ignored for the purpose of prior generation. We use neural density estimation to model the measured distribution of EFT parameters. Our distribution model is then used as a prior in a reanalysis of the Baryon Oscillation Spectroscopic Survey (BOSS) full-shape galaxy power spectrum data. Assuming the $\Lambda \mathrm{CDM}$ model, we find significant ($\approx 30\%$ and $\approx 60\%$) improvements for the matter density fraction and the mass fluctuation amplitude, which are constrained to ${\Omega }_m=0.315\pm 0.010$ and ${\sigma }_8=0.671\pm 0.027$. The value of the Hubble constant does not change, $H_0=68.7\pm 1.1\mathrm{km}/s/\mathrm{Mpc}$. This reaffirms earlier reports of the structure growth tension from the BOSS data. Finally, we use the measured EFT parameters to constrain the galaxy-dark matter connection.

Cosmological models are often motivated and formulated in the language of particle physics, using quantities such as the axion decay constant, but tested against data using ostensibly physical quantities, such as energy density ratios, assuming uniform priors on the latter. This approach neglects priors on the model from fundamental theory, including from particle physics and string theory, such as the preference for sub-Planckian axion decay constants. We introduce a novel approach to learning theory-informed priors for Bayesian inference using normalizing flows (NF), a flexible generative machine learning technique that generates priors on model parameters when analytic expressions are unavailable or difficult to compute. As a test case, we focus on early dark energy (EDE), a model designed to address the Hubble tension. Rather than using uniform priors on the phenomenological EDE parameters $f_{\mathrm{E}\mathrm{D}\mathrm{E}}$ and $z_c$, we train a NF on EDE cosmologies informed by theory expectations for axion masses and decay constants. Our method recovers known constraints in this representation while being $\sim 300,000$ times more efficient in terms of total CPU compute time. Applying our NF to $\text{Planck}$ and BOSS data, we obtain the first theory-informed constraints on EDE, finding $f_{\mathrm{E}\mathrm{D}\mathrm{E}}\lesssim 0.02$ at $95\mathrm{\%}$ confidence with an $H_0$ consistent with Planck, but in $\sim 6\sigma$tension with SH0ES. This yields the strongest constraints on EDE to date, additionally challenging its role in resolving the Hubble tension.

This note serves as a companion to a Letter, where we introduce a new energy correlator-based observable designed to minimize the impact of selection bias due to energy loss in inclusive jets in heavy-ion collisions. Here, we apply the method outlined in the Letter to the first-ever measurement of energy correlators in heavy-ion collisions, recently released by the CMS Collaboration.

The first-ever measurement of energy correlators within inclusive jets produced in heavy-ion collisions, revealed by the CMS Collaboration, shows a clear enhancement at large angles relative to the proton-proton (p-p) baseline. However, interpreting this enhancement is complicated due to selection bias from energy loss, which also distorts the energy correlator heavy-ion to p-p ratio in the hadronization region, hindering our understanding of parton/hadron dynamics in a colored medium. In this Letter, we introduce a new ratio of energy correlator observables that removes the leading effects of selection bias from the two-point energy correlator spectrum (E2C). pythia and herwig simulations show that the impact of selection bias in the E2C is reduced by an order of magnitude, while sensitivity to any other medium modifications is retained. This quantity can be obtained directly from the experimental measurements presented by CMS, as illustrated in an accompanying note [C. Andres and J. Holguin, ].

Quantum thermalization describes how closed quantum systems can effectively reach thermal equilibrium, resolving the apparent incongruity between the reversibility of Schrödinger’s equation and the irreversible entropy growth dictated by the second law of thermodynamics. Despite its ubiquity and conceptual significance, a complete proof of quantum thermalization has remained elusive for several decades. Here, we prove that quantum thermalization must occur in any qubit system with local interactions satisfying three conditions: (i) high effective temperature, (ii) translation invariance, and (iii) no perfect resonances in the energy spectrum. Specifically, we show that a typical, initially unentangled pure state drawn from any ensemble with maximum entropy becomes locally indistinguishable from a Gibbs state upon unitary evolution. Our proof relies on a recent breakthrough in quantum information theory proving the separability of high-temperature thermal states, as well as on new technical results identifying sufficient conditions for quantum thermalization, whose applicability extends beyond our main result. Our work illustrates that statistical physics can be understood as an emergent phenomenon, explicitly derived from the first principles of quantum mechanics.

Quantum signal processing (QSP) provides a systematic framework for implementing a polynomial transformation of a linear operator, and unifies nearly all known quantum algorithms. In parallel, recent works have developed randomized compiling, a technique that promotes a unitary gate to a quantum channel and enables a quadratic suppression of error (i.e., ϵ → O(ϵ${}^{}$)) at little to no overhead. Here we integrate randomized compiling into QSP through Stochastic Quantum Signal Processing. Our algorithm implements a probabilistic mixture of polynomials, strategically chosen so that the average evolution converges to that of a target function, with an error quadratically smaller than that of an equivalent individual polynomial. Because nearly all QSP-based algorithms exhibit query complexities scaling as $O(\log (1/ϵ))$—stemming from a result in functional analysis—this error suppression reduces their query complexity by a factor that asymptotically approaches 1/2. By the unifying capabilities of QSP, this reduction extends broadly to quantum algorithms, which we demonstrate on algorithms for real and imaginary time evolution, phase estimation, ground state preparation, and matrix inversion.

We derive an analog of Mirzakhani’s recursion relation for hyperbolic string vertices and investigate its implications for closed string field theory. Central to our construction are systolic volumes: the Weil-Petersson volumes of regions in moduli spaces of Riemann surfaces whose elements have systoles L ≥ 0. These volumes can be shown to satisfy a recursion relation through a modification of Mirzakhani’s recursion as long as L ≤ 2 sinh${}^{}$ 1. Applying the pants decomposition of Riemann surfaces to off-shell string amplitudes, we promote this recursion to hyperbolic string field theory and demonstrate the higher order vertices are determined by the cubic vertex iteratively for any background. Such structure implies the solutions of closed string field theory obey a quadratic integral equation. We illustrate the utility of our approach in an example of a stubbed scalar theory.

We point out that dark matter can illuminate planetary skies via ultraviolet airglow. Dark matter annihilation products can excite molecular hydrogen, which then deexcites to produce ultraviolet emission in the Lyman and Werner bands. We search for this new effect by analyzing nightside ultraviolet radiation data from Voyager 1, Voyager 2, and New Horizons flybys of Neptune, Uranus, Saturn, and Jupiter. Our findings set new constraints on the dark matter-nucleon scattering cross section down to about ${10}^{-40}$cm${}^2$. We highlight that future ultraviolet airglow measurements of Solar System planets or other worlds provide a new dark matter discovery avenue.

The temperature and polarization anisotropies of the cosmic microwave background have been used to set constraints on decaying dark matter models down to keV masses. In this work, we extend these limits down to the sub-keV mass range. Using principal component analysis, we estimate the lower bound on the decay lifetime for a basis of different dark matter masses and Standard Model final states, from which the bound on an arbitrary model can be calculated. We validate our principal component analysis using Markov chain Monte Carlo methods and Planck 2018 data. We perform a separate analysis for models decaying into photons below the hydrogen ionization threshold. We demonstrate that, for these models, the effect of energy deposition can be captured approximately by a single parameter, but the redshift dependence of the effect is very different from higher-energy injections; in particular, the perturbations to cosmic microwave background anisotropies are more sensitive to energy deposited around the time of recombination.

A system with a one-form global symmetry at finite temperature can be viewed as a dissipative fluid of string-like objects. In this work, we classify and construct the most general effective field theories for hydrodynamics of such string fluids, in a probe limit where the one-form charge density is decoupled from the energy-momentum tensor. We show that at leading order in the derivative expansion, there are two distinct types of diffusive transport in a string fluid depending on the discrete spacetime symmetries present in it. One particular application of interest is magnetohydrodynamics (MHD), where the effective field theories describe the diffusion of magnetic field lines. Due to the distinction between the effective field theories for different discrete symmetries, we show that the MHD of a fluid with charge conjugation symmetry is qualitatively different from that of a neutron star, which we previously discussed in arXiv:2207.01636. The explicit effective actions that we write down can be used to obtain the dispersion relations ω(κ) up to cubic order in momenta for each of the different discrete symmetry choices. As another application of this formalism, we show that when the one-form symmetry is spontaneously broken, the effective action reduces to the Maxwell theory. This confirms the interpretation of the photon as a Goldstone boson arising from the spontaneous breaking of a one-form symmetry.

We take a first step towards developing a new language to describe causal structure, event horizons, and quantum extremal surfaces (QES) for the bulk description of holographic systems beyond the standard Einstein gravity regime. By considering the structure of boundary operator algebras, we introduce a stringy
causal depth parameter”, which quantifies the depth of the emergent radial direction in the bulk, and a certain notion of ergodicity on the boundary. We define stringy event horizons in terms of the half-sided inclusion property, which is related to a stronger notion of boundary ergodic or quantum chaotic behavior. Using our definition, we argue that above the Hawking–Page temperature, there is an emergent sharp horizon structure in the large N limit of N=4 Super-Yang–Mills at finite nonzero ‘t Hooft coupling. In contrast, some previously considered toy models of black hole information loss do not have a stringy horizon. Our methods can also be used to probe violations of the equivalence principle for the bulk gravitational system, and to explore aspects of stringy nonlocality.

The so-called metastability bound asserts that an unnaturally small Higgs mass is a necessary condition for electroweak vacuum metastability, offering a new approach toward solving the hierarchy problem. So far, this result relies on the assumption of a negative Higgs mass parameter, or equivalently, on electroweak spontaneous symmetry breaking. We derive a new, corresponding bound for the case of a positive mass parameter. When the new bound is significantly more restrictive than or comparable to its established counterpart, it may offer an explanation for the sign of the Higgs mass parameter, and thus, spontaneous symmetry breaking itself. New physics at scales δ(1–10) TeV can lower these bounds as far as the TeV scale. As an illustration, we consider vacuum stability in the presence of additional TeV-scale fermions with Yukawa couplings to the Higgs, as well as a dimension-six term parametrizing new physics in the UV. This scenario requires new physics that couples strongly to the Higgs, and can potentially be probed at future colliders. Finally, to allow for comparison with concrete mechanisms predicting metastability, we provide the mass-dependent lifetime of the electroweak vacuum for this model.

The abundance of massive primordial black holes has historically been constrained by dynamical probes. Since these objects can participate in hard few-body scattering processes, they can readily transfer energy to stellar systems and, in particular, disrupt wide binaries. However, disruption is not the only possible outcome of such few-body processes. Primordial black holes could also participate in exchange processes, in which one component of a binary system is ejected and replaced by the black hole itself. In this case, the remaining object in the binary would dynamically appear to have an invisible companion. We study the rate of exchange processes for primordial black holes as a component of dark matter and evaluate possible mechanisms for detecting such binaries. We find that many such binaries plausibly exist in the Solar neighborhood and show that this process can account for observed binary systems whose properties run counter to the predictions of isolated binary evolution.

The abundance of massive primordial black holes has historically been constrained by dynamical probes. Since these objects can participate in hard few-body scattering processes, they can readily transfer energy to stellar systems and, in particular, disrupt wide binaries. However, disruption is not the only possible outcome of such few-body processes. Primordial black holes could also participate in exchange processes, in which one component of a binary system is ejected and replaced by the black hole itself. In this case, the remaining object in the binary would dynamically appear to have an invisible companion. We study the rate of exchange processes for primordial black holes as a component of dark matter and evaluate possible mechanisms for detecting such binaries. We find that many such binaries plausibly exist in the Solar neighborhood and show that this process can account for observed binary systems whose properties run counter to the predictions of isolated binary evolution.

The black hole entropy has been observed to generically turn negative at exponentially low temperatures \(T\sim e^{-S_0}\) in the extremal Bekenstein-Hawking entropy \(S_0\), a seeming pathology often attributed to missing non-perturbative effects. In fact, we show that this negativity must happen for any effective theory of quantum gravity with an ensemble description. To do so, we identify the usual gravitational entropy as an annealed entropy \(S_a\), and prove that this quantity gives \(S_0\) at extremality if and only if the ground-state energy is protected by supersymmetry, and diverges negatively otherwise. The actual thermodynamically-behaved quantity is the average or quenched entropy \(S_q\), whose calculation is poorly understood in gravity: it involves replica wormholes in a regime where the topological expansion breaks down. Using matrix integrals we find new instanton saddles that dominate gravitational correlators at \(T\sim e^{-S_0}\) and are dual to semiclassical wormholes involving dynamical branes. These brane solutions give the leading contribution to any black hole very near extremality, and a duality with matrix ensembles would not make sense without them. In the non-BPS case, they are required to make \(S_q\) non-negative and also enhance the negativity of \(S_a\), both effects consistent with matrix integrals evaluated exactly. Our instanton results are tested against the on-shell action of D3-branes dual to multiply wrapped Wilson loops in \(\mathcal{N}=4\) super-YM, and a precise match is found. Our analysis of low-energy random matrix spectra also explains the origin of spectral gaps in supersymmetric theories, not only when there are BPS states at zero energy, but also for purely non-BPS supermultiplets. In the former, our prediction for the gap in terms of the degeneracy of BPS states agrees with the R-charge scaling in gapped multiplets of \(\mathcal{N}=2\) super-JT gravity.

As the partons in a jet propagate through the quark-gluon plasma (QGP) produced in a heavy-ion collision, they lose energy to, kick, and are kicked by the medium. The resulting modifications to the parton shower encode information about the microscopic nature of QGP. The momentum and energy lost by the parton shower are gained by the medium and, since QGP is a strongly coupled liquid, this means that the jet excites a wake in the droplet of QGP. After freezeout, this wake becomes soft hadrons with net momentum in the jet direction meaning that reconstructed jets include hadrons originating from both the modified parton shower and its wake. This makes it challenging to find an unambiguous experimental view of the response of a droplet of QGP to a jet. Recent years have seen significant advances in the understanding of the substructure of jets using correlation functions of the energy flux operator. So far, such studies have focused primarily on the two-point correlator, which serves to identify the angular scale of the underlying dynamics. Higher-point correlators hold the promise of mapping out the dynamics themselves. We perform the first study of the shape-dependent three-point energy-energy-energy correlator in heavy-ion collisions. Using the Hybrid Model to simulate the interactions of high energy jets with QGP, we show that hadrons originating from wakes are the dominant contribution to the three-point correlator in the regime where the three points are well-separated in angle, forming a roughly equilateral triangle. This equilateral region of the correlator is far from the region populated by collinear vacuum emissions, making it a canvas on which jet wakes can be imaged. Our work is a key step towards the systematic use of energy correlators to image and unravel the dynamical response of a droplet of QGP to a passing jet, and motivates many experimental and theoretical studies.

Deconvolving (“unfolding”) detector distortions is a critical step in the comparison of cross section measurements with theoretical predictions in particle and nuclear physics. However, most existing approaches require histogram binning while many theoretical predictions are at the level of statistical moments. We develop a new approach to directly unfold distribution moments as a function of another observable without having to first discretize the data. Our Moment Unfolding technique uses machine learning and is inspired by Generative Adversarial Networks (GANs). We demonstrate the performance of this approach using jet substructure measurements in collider physics. With this illustrative example, we find that our Moment Unfolding protocol is more precise than bin-based approaches and is as or more precise than completely unbinned methods.