MIT Center for Theoretical Physics — a Leinweber Institute Publications

As of today there exist consistent, gauge-invariant string field theories describing all string theories: bosonic open and closed strings, open superstrings, heterotic strings and type II strings. The construction of these theories require algebraic ingredients, such as \(A_\infty\) and \(L_\infty\) homotopy algebras, geometric ingredients, relevant to the building of moduli spaces of Riemann surfaces and the distribution of picture changing operators, and field-theoretic ingredients, involving two-dimensional CFT's and BCFT's and Batalin-Vilkovisky quantization. Applications of string field theory include the description of non-perturbative phenomena such as tachyon condensation and classical solutions, and the resolution of a number of ambiguities that bedevil the world-sheet formulation of perturbative string theory. It also allows, given a proper definition of contours of integration for momenta, for a proof of unitarity and a clear understanding of the ultraviolet finiteness of the theory. In this article we review these developments. This is an expanded version of the review written for the ``Handbook of Quantum Gravity'', eds. C. Bambi, L. Modesto and I.Shapiro.

Energy correlators provide a powerful observable to study fragmentation dynamics in QCD. We demonstrate that the leading nonperturbative corrections for projected \(N\)-point energy correlators are described by the same universal parameter for any \(N\), which has already been determined from other event shape fits. Including renormalon-free nonperturbative corrections substantially improves theoretical predictions of energy correlators, notably the transition into the confining region at small angles. Nonperturbative corrections are shown to have a significant impact on \(\alpha_s\) extractions.

"Attractor" solutions for the pre-hydrodynamic, far-from-equilibrium, evolution of the matter produced in relativistic heavy ion collisions have emerged as crucial descriptors of the rapid hydrodynamization of quark-gluon plasma (QGP). Adiabatic Hydrodynamization (AH) has been proposed as a framework with which to describe, explain, and predict attractor behavior that draws upon an analogy to the adiabatic approximation in quantum mechanics. In this work, we systematize the description of pre-hydrodynamic attractors in kinetic theory by showing how to use the AH framework to identify these long-lived solutions to which varied initial conditions rapidly evolve, demonstrating the robustness of this framework. In a simplified QCD kinetic theory in the small-angle scattering limit, we use AH to explain both the early- and late-time scaling behavior of a longitudinally expanding gluon gas in a unified framework. In this context, we show that AH provides a unified description of, and intuition for, all the stages of what in QCD would be bottom-up thermalization, starting from a pre-hydrodynamic attractor and ending with hydrodynamization. We additionally discuss the connection between the notions of scaling behavior and adiabaticity and the crucial role of time-dependent coordinate redefinitions in identifying the degrees of freedom of kinetic theories that give rise to attractor solutions. The tools we present open a path to the intuitive explanation of how attractor behavior arises and how the attractor evolves in all stages of the hydrodynamization of QGP in heavy ion collisions.

Extracting scientific understanding from particle-physics experiments requires solving diverse learning problems with high precision and good data efficiency. We propose the Lorentz Geometric Algebra Transformer (L-GATr), a new multi-purpose architecture for high-energy physics. L-GATr represents high-energy data in a geometric algebra over four-dimensional space-time and is equivariant under Lorentz transformations, the symmetry group of relativistic kinematics. At the same time, the architecture is a Transformer, which makes it versatile and scalable to large systems. L-GATr is first demonstrated on regression and classification tasks from particle physics. We then construct the first Lorentz-equivariant generative model: a continuous normalizing flow based on an L-GATr network, trained with Riemannian flow matching. Across our experiments, L-GATr is on par with or outperforms strong domain-specific baselines.

We update and extend a previous N\(^3\)LL\(^\prime\)+\({\cal O}(\alpha_s^3)\) strong coupling determination from thrust data. In particular, we carry out a fit with data fully restricted to the dijet region seeking to minimize the potential impact of power corrections that go beyond dijet configurations. In addition, we parametrize deviations from the dijet power correction in order to add an additional source of uncertainty in the result for \(\alpha_s(m_Z)\). We also show that the inclusion of resummation is important to achieve stability with respect to varying the fit region.

We demonstrate that "natural inflation", also known as "axion inflation", can be compatible with Planck 2018 measurements of the cosmic microwave background, while predicting an exponentially small tensor-to-scalar ratio, e.g., \(r\sim 10^{-15}\). The strong suppression of \(r\) arises from dynamics of the radial component of the complex scalar field, whose phase is the axion. Such tiny values of \(r\) remain well below the threshold for detection by CMB-S4 or Simons Observatory B-mode searches. The model is testable with the running \(\alpha_s\) of the spectral index, which is within reach of next-generation CMB and large-scale structure experiments, motivating the running as a primary science goal for future experiments.

We present an effective field theory (EFT) approach to extract fundamental cosmological parameters from the Lyman-alpha forest flux fluctuations as an alternative to the standard simulation-based techniques. As a first application, we re-analyze the publicly available one-dimensional Lyman-alpha flux power spectrum (FPS) data from the Sloan Digital Sky Survey (SDSS). Our analysis relies on informative priors on EFT parameters which we extract from a combination of public hydrodynamic simulation and emulator data. Assuming the concordance cosmological model, our one-parameter analysis yields a \(2\%\) measurement of the late time mass fluctuation amplitude \(\sigma_8 = 0.841\pm 0.017\), or equivalently,the structure growth parameter \(S_8 = 0.852\pm 0.017\), consistent with the standard cosmology. Combining our EFT likelihood with Planck + baryon acoustic oscillation data, we find a new constraint on the total neutrino mass, \(\sum m_\nu<\) 0.08 eV (at 95% CL). Our study defines priorities for the development of EFT methods and sets the benchmark for cosmological analyses of the Lyman-alpha forest data from the Dark Energy Spectroscopic Instrument.

Sen's formalism is a mechanism for eliminating constraints on the dynamical fields that are imposed independently from equations of motion by employing spurious free fields. In this note a cyclic homotopy associative algebra underlying Sen's formalism is developed. The novelty lies in the construction of a symplectic form and cyclic \(A_\infty\) maps on an extended algebra that combines the dynamical and spurious fields. This algebraic presentation makes the gauge invariance of theories using Sen's formulation manifest.

We use the factorization constraints of open-closed string field theory to determine the signs and normalizations of general string amplitudes with both open and closed string external states. The normalization of all amplitudes is controlled by the genus, the number of boundaries, the number of open and closed string insertions, the string coupling and the D-brane tension. The challenge with signs arises because the relevant moduli spaces are not complex manifolds and have no obvious orientation. We deal with this by fixing a specific convention for the sign of the integration measure over the moduli space and adopting a consistent prescription for the ordering of operators and ghost insertions inside correlators.

The properties of the Higgs potential are determined by three parameters: the mass parameter, the quartic self-coupling, and a constant term. Remarkably, all three of these parameters seem subject to a significant amount of fine-tuning. All these tunings can be seen as their corresponding parameters being close to critical values marking quantum phase transitions. While such behavior is surprising from a conventional particle physics perspective, it is a common feature of dynamical systems. This has motivated the conjecture that the values of the Higgs' parameters are the result of some dynamical mechanism. This possibility suggests the construction of mechanisms dynamically choosing sets of Higgs parameters. In these notes, I discuss a complementary approach. Taking seriously that such a mechanism could exist in nature, it is plausible to assume that it also influences Beyond-Standard-Model physics. This suggests considering near-criticality in any model of interest and investigating its consequences more generally, in particular independent of a concrete mechanism responsible for it. I first explain what it means for the parameters of the Higgs potential to be near-critical. This includes a discussion of the recently discovered "metastability bound" on the Higgs mass, which can be understood through a critical point. I then review two concrete examples of mechanisms in which the parameters of the Higgs potential are dynamically driven towards critical values. These mechanisms also serve as an important proof on concept for the feasibility of the assumptions at the foundation of these notes. Using a simple example for concreteness, the final part of these notes then explicitly demonstrates how to approach a given model in the light of the near-criticality conjecture.

It has been shown that entropy differences between certain states of perturbative quantum gravity can be computed without specifying an ultraviolet completion. This is analogous to the situation in classical statistical mechanics, where entropy differences are defined but absolute entropy is not. Unlike in classical statistical mechanics, however, the entropy differences computed in perturbative quantum gravity do not have a clear physical interpretation. Here we construct a family of perturbative black hole states for which the entropy difference can be interpreted as a relative counting of states. Conceptually, this paper begins with the algebra of mass fluctuations around a fixed black hole background, and points out that while this is a type I algebra, it is not a factor and therefore has no canonical definition of entropy. As in previous work, coupling the mass fluctuations to quantum matter embeds the mass algebra within a type II factor, in which entropy differences (but not absolute entropies) are well defined. It is then shown that for microcanonical wavefunctions of mass fluctuation, the type II entropy difference equals the logarithm of the dimension of the extra Hilbert space that is needed to map one microcanonical window to another using gauge-invariant unitaries. The paper closes with comments on type II entropy difference in a more general class of states, where the von Neumann entropy difference does not have a physical interpretation, but "one-shot" entropy differences do.

X-space schemes are gauge-invariant, regulator-independent renormalization schemes that are defined by requiring position-space correlation functions of gauge invariant operators to be equal to their noninteracting values at particular kinematic points. These schemes can be used to nonperturbatively renormalize composite operators in Lattice Quantum Chromodynamics (LQCD), and by computing matching coefficients between the X-space scheme and MSbar in the dimensionally-regulated continuum, matrix elements calculated with LQCD can be converted to MSbar-renormalized matrix elements. Using X-space schemes for Heavy Quark Effective Theory (HQET) operators has the additional benefit that appropriate ratios of position-space correlation functions cancel the power divergent static-quark self-energy of Lattice HQET nonperturbatively. This work presents the O(\(\alpha_S\)) matching coefficients between X-space renormalized four-quark flavor-nonsinglet HQET operators relevant for the lifetimes of charm- and bottom-hadrons, and four-quark HQET operators relevant for mixing between neutral mesons containing a heavy quark, such as B-Bbar mixing.

We present the first constraints on tidal heating for the binary systems detected in the LIGO-Virgo-KAGRA (LVK) gravitational wave data. Tidal heating, also known as tidal dissipation, characterizes the viscous nature of an astrophysical body and provides a channel for exchanging energy and angular momentum with the tidal environment. Using the worldline effective field theory formalism, we introduce a physically motivated and easily interpretable parametrization of tidal heating valid for an arbitrary compact astrophysical object. We then derive the imprints of the spin-independent and linear-in-spin tidal heating effects of generic binary components on the waveform phases and amplitudes of quasi-circular orbits. Notably, the mass-weighted spin-independent tidal heating coefficient derived in this work, \(\mathcal{H}_0\), is the dissipative analog of the tidal Love number. We constrain the tidal heating coefficients using the public LVK O1-O3 data. Our parameter estimation study includes two separate analyses: the first treats the catalog of binary events as binary black holes (BBH), while the second makes no assumption about the nature of the binary constituents and can therefore be interpreted as constraints for exotic compact objects. In the former case, we combine the posterior distributions of the individual BBH events and obtain a joint constraint of \(-13 < \mathcal{H}_0 < 20\) at the \(90\%\) credible interval for the BBH population. This translates into a bound on the fraction of the emitted gravitational wave energy lost due to tidal heating (or gained due to radiation enhancement effects) at \(|\Delta E_H/\Delta E_{\infty}|\lesssim 3\cdot 10^{-3}\). Our work provides the first robust framework for deriving and measuring tidal heating effects in merging binary systems, demonstrating its potential as a powerful probe of the nature of binary constituents and tests of new physics.

The low-energy finite-volume spectrum of the two-nucleon system at a pion mass of \(m_\pi \approx 806\) MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Phys.Rev.D 107 (2023) 9, 094508] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, these results provide compelling evidence for the presence of an additional state in the low-energy spectrum not present in a non-interacting two-nucleon system. Within the space of operators that are considered in this work, the strongest overlap of this state is onto hexaquark operators, and variational analyses based on operator sets containing only dibaryon operators are insensitive to this state over the temporal extent studied here. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated.

Normalizing flows are machine-learned maps between different lattice theories which can be used as components in exact sampling and inference schemes. Ongoing work yields increasingly expressive flows on gauge fields, but it remains an open question how flows can improve lattice QCD at state-of-the-art scales. We discuss and demonstrate two applications of flows in replica exchange (parallel tempering) sampling, aimed at improving topological mixing, which are viable with iterative improvements upon presently available flows.

Scale separation is an important physical principle that has previously enabled algorithmic advances such as multigrid solvers. Previous work on normalizing flows has been able to utilize scale separation in the context of scalar field theories, but the principle has been largely unexploited in the context of gauge theories. This work gives an overview of a new method for generating gauge fields using hierarchical normalizing flow models. This method builds gauge fields from the outside in, allowing different parts of the model to focus on different scales of the problem. Numerical results are presented for \(U(1)\) and \(SU(3)\) gauge theories in 2, 3, and 4 spacetime dimensions.

We study the general framework of effective theories obtained via exact path integration of a complete theory over some sector of its degrees of freedom. Theories constructed this way contain multi-integrals which couple fields arbitrarily far apart, and in certain settings even on path-disconnected components of the space. These are not just entanglement, but genuine non-local interactions that we dub quantum wormholes. Any state the path integral of such an effective theory prepares is shown to be a partial trace of a state of the complete theory over the integrated-out sector. The resulting reduced density operator is generally mixed due to bra-ket wormholes. An infinite family of ensembles of pure states of the complete theory giving the same effective state is identified. These allow one to equivalently interpret any effective state as being prepared by an ensemble of theories. When computing entropic quantities, bra-ket wormholes give rise to replica wormholes. This causes replica path integrals for the effective theory to not factorize even when the underlying manifold does, as expected from mixing. In contrast, effective theories obtained by derivative expansions have no quantum wormholes and prepare pure states. There exist operators in the algebra of effective theories which can distinguish mixed from pure states, implying a breakdown of non-exact effective theories for sufficiently complex observables. This framework unifies and provides new insights into much of the phenomena observed in quantum gravity, including the interplay between wormholes and unitarity, the breakdown of bulk effective theory, the factorization puzzle, state ensembles, theory ensembles, quantum error correction, and baby universes. Some interesting lessons are drawn accounting also for characteristic aspects of gravity concerning IR/UV mixing and Kaluza-Klein reductions.

The transverse momentum-dependent fragmentation functions (TMD FFs) of heavy (bottom and charm) quarks, which we recently introduced, are universal building blocks that enter predictions for a large number of observables involving final-state heavy quarks or hadrons. They enable the extension of fixed-order subtraction schemes to quasi-collinear limits, and are of particular interest in their own right as probes of the nonperturbative dynamics of hadronization. In this paper we calculate all TMD FFs involving heavy quarks and the associated TMD matrix element in heavy-quark effective theory (HQET) to next-to-leading order in the strong interaction. Our results confirm the renormalization properties, large-mass, and small-mass consistency relations predicted in our earlier work. We also derive and confirm a prediction for the large-\(z\) behavior of the heavy-quark TMD FF by extending, for the first time, the formalism of joint resummation to capture quark mass effects in heavy-quark fragmentation. Our final results in position space agree with those of a recent calculation by another group that used a highly orthogonal organization of singularities in the intermediate momentum-space steps, providing a strong independent cross check. As an immediate application, we present the complete quark mass dependence of the energy-energy correlator (EEC) in the back-to-back limit at \(\mathcal{O}(\alpha_s)\).

We study future lepton collider prospects for testing predictive models of leptophilic dark matter candidates with a thermal origin. We calculate experimental milestones for testing the parameter space compatible with freeze-out and the associated collider signals at past, present, and future facilities. This analysis places new limits on such models by leveraging the utility of lepton colliders. At \(e^+e^-\) machines, we make projections using precision \(Z\)-pole observables from \(e^+e^-\to l^+l^- + \) missing energy signatures at LEP and future projections for FCC-ee in these channels. Additionally, a muon collider could also probe new thermal relic parameter space in this scenario via \(\mu^+\mu^- \to X + \) missing energy where \(X\) is any easy identifiable SM object. Collectively, these processes can probe much all of the parameter space for which DM direct annihilation to \(l^+l^-\) yields the observed relic density in Higgs-like models with mass-proportional couplings to charged leptons.

Massive particles leave imprints on primordial non-Gaussianity via couplings to the inflaton, even despite their exponential dilution during inflation: practically, the Universe acts as a Cosmological Collider. We present the first dedicated search for spin-zero particles using BOSS redshift-space galaxy power spectrum and bispectrum multipoles, as well as Planck CMB non-Gaussianity data. We demonstrate that some Cosmological Collider models are well approximated by the standard equilateral and orthogonal parametrization; assuming negligible inflaton self-interactions, this facilitates us translating Planck non-Gaussianity constraints into bounds on Collider models. Many models have signatures that are not degenerate with equilateral and orthogonal non-Gaussianity and thus require dedicated searches. Here, we constrain such models using BOSS three-dimensional redshift-space galaxy clustering data, focusing on spin-zero particles in the principal series and constraining their couplings to the inflaton at varying speed and mass, marginalizing over the unknown inflaton self-interactions. This is made possible through an improvement in Cosmological Bootstrap techniques and the combination of perturbation theory and halo occupation distribution models for galaxy clustering. Our work sets the standard for inflationary spectroscopy with cosmological observations, providing the ultimate link between physics on the largest and smallest scales.