MIT Center for Theoretical Physics — a Leinweber Institute Publications

We study the large \(N_\text{c}\) scaling of meson-meson scattering amplitudes in a theory with \(N_\text{f}=4\) degenerate quark flavors. We focus on two different scattering channels, one having the same quantum numbers as some recently found tetraquark states at LHCb. Using Lüscher’s formalism, we study the \(N_\text{c}\) dependence of the scattering phase shift and investigate the presence of exotic resonances in the scattering amplitude. We analyze the impact of including two-vector-meson and tetraquark-like operators to extract the finite-volume energies.

We review recent progress in understanding quarkonium dynamics inside the quark-gluon plasma as an open quantum system with a focus on the definition and nonperturbative calculations of relevant transport coefficients and generalized gluon distributions.

The high energy (Regge) limit provides a playground for understanding all loop structures of scattering amplitudes, and plays an important role in the description of many phenomenologically relevant cross-sections. While well understood in the planar limit, the structure of non-planar corrections introduces many fascinating complexities, for which a general organizing principle is still lacking. We study the structure of multi-reggeon exchanges in the context of the effective field theory for forward scattering, and derive their factorization into collinear operators (impact factors) and soft operators. We derive the structure of the renormalization group consistency equations in the effective theory, showing how the anomalous dimensions of the soft operators are related to those of the collinear operators, allowing us to derive renormalization group equations in the Regge limit purely from a collinear perspective. The rigidity of the consistency equations provides considerable insight into the all orders organization of Regge amplitudes in the effective theory, as well as its relation to other approaches. Along the way we derive a number of technical results that improve the understanding of the effective theory. We illustrate this collinear perspective by re-deriving all the standard BFKL equations for two-Glauber exchange from purely collinear calculations, and we show that this perspective provides a number of conceptual and computational advantages as compared to the standard view from soft or Glauber physics. We anticipate that this formulation in terms of collinear operators will enable a better understanding of the relation between BFKL and DGLAP in gauge theories, and facilitate the analysis of renormalization group evolution equations describing Reggeization beyond next-to-leading order.

We consider phenomenological aspects of a natural class of Standard Model-like supersymmetric F-theory vacua realized through flux breaking of rigid \(E_7\) gauge factors. Three generations of Standard Model matter are realized in many of these vacua. We further find that many other Standard Model-like features are naturally compatible with these constructions. For example, dimension-4 and 5 terms associated with proton decay are ubiquitously suppressed. Many of these features are due to the group theoretical structure of \(E_7\) and associated F-theory geometry. In particular, a set of approximate global symmetries descends from the \(E_7\) group, leading to exponential suppression of undesired couplings.

Primordial black holes (PBHs) remain a viable dark matter candidate in the asteroid-mass range. We point out that in this scenario, the PBH abundance would be large enough for at least one object to cross through the inner Solar System per decade. Since Solar System ephemerides are modeled and measured to extremely high precision, such close encounters could produce detectable perturbations to orbital trajectories with characteristic features. We evaluate this possibility with a suite of simple Solar System simulations, and we argue that the abundance of asteroid-mass PBHs can plausibly be probed by existing and near-future data.

We present an operator based factorization formula for the transverse energy-energy correlator in the back-to-back (dijet) region, and uncover its remarkable perturbative simplicity and relation to transverse momentum dynamics. This simplicity enables us to achieve next-to-next-to-next-to leading logarithmic (N\(^3\)LL) accuracy for a hadron collider dijet event shape for the first time. Our factorization formula applies to color singlet, \(W/Z/\gamma\) + jet, and dijet production, providing a natural generalization of transverse momentum observables to one- and two-jet final states. This provides a laboratory for precision studies of QCD and transverse momentum dynamics at hadron colliders, as well as an opportunity for understanding factorization and its violation in a perturbatively well controlled setting.

Early matter-dominated eras (EMDEs) are a natural feature arising in many models of the early universe and can generate a stochastic gravitational wave background (SGWB) during the transition from an EMDE to the radiation-dominated universe required by the time of Big Bang Nucleosynthesis. While there are calculations of the SGWB generated in the linear regime, no detailed study has been made of the nonlinear regime. We perform the first comprehensive calculation of GW production in the nonlinear regime, using a hybrid \(N\)-body and lattice simulation to study GW production from both a metastable matter species and the radiation produced in its decay. We find that nonlinearities significantly enhance GW production up to frequencies at least as large as the inverse light-crossing time of the largest halos that form prior to reheating. The resulting SGWB is within future observational reach for curvature perturbations as small as those probed in the cosmic microwave background, depending on the reheating temperature. Out-of-equilibrium dynamics could further boost the induced SGWB, while a fully relativistic gravitational treatment is required to resolve the spectrum at even higher frequencies.

We report on progress in the nonperturbative understanding of quarkonium dynamics inside a thermal plasma. The time evolution of small-size quarkonium is governed by two-point correlation functions of chromoelectric fields dressed with an adjoint Wilson line, known in this context as generalized gluon distributions (GGDs). The GGDs have been calculated in both weakly and strongly coupled plasmas by using perturbative and holographic methods. Strikingly, the results of our calculations for a strongly coupled plasma indicate that the quarkonium dissociation and recombination rates vanish in the transport descriptions that assume quarkonium undergoes Markovian dynamics. However, this does not imply that the dynamics is trivial. As a starting point to explore the phenomenological consequences of the result at strong coupling, we show a calculation of the \(\Upsilon(1S)\) formation probability in time-dependent perturbation theory. This is a first step towards the development of a transport formalism that includes non-Markovian effects, which, depending on how close the as of yet undetermined nonperturbative QCD result of the GGDs is to the strongly coupled \(\mathcal{N}=4\) SYM result, could very well dominate over the Markovian ones in quark-gluon plasma produced at RHIC and the LHC.

The 21-cm signal provides a novel avenue to measure the thermal state of the universe during cosmic dawn and reionization (redshifts \(z\sim 5-30\)), and thus to probe energy injection from decaying or annihilating dark matter (DM). These DM processes are inherently inhomogeneous: both decay and annihilation are density dependent, and furthermore the fraction of injected energy that is deposited at each point depends on the gas ionization and density, leading to further anisotropies in absorption and propagation. In this work, we develop a new framework for modeling the impact of spatially inhomogeneous energy injection and deposition during cosmic dawn, accounting for ionization and baryon density dependence, as well as the attenuation of propagating photons. We showcase how this first completely inhomogeneous treatment affects the predicted 21-cm power spectrum in the presence of exotic sources of energy injection, and forecast the constraints that upcoming HERA measurements of the 21-cm power spectrum will set on DM decays to photons and to electron/positron pairs. These projected constraints considerably surpass those derived from CMB and Lyman-\(\alpha\) measurements, and for decays to electron/positron pairs they exceed all existing constraints in the sub-GeV mass range, reaching lifetimes of \(\sim 10^{28}\,\mathrm{s}\). Our analysis demonstrates the unprecedented sensitivity of 21-cm cosmology to exotic sources of energy injection during the cosmic dark ages. Our code, \(\mathtt{DM21cm}\), includes all these effects and is publicly available in an accompanying release.

Motivated by bit threads, we introduce a new prescription for computing entropy vectors outside the holographic entropy cone. By utilizing cycle flows on directed graphs, we show that the maximum cycle flow associated to any subset of vertices, which corresponds to a subsystem, manifestly obeys purification symmetry. Furthermore, we prove that the maximum cycle flow obeys both subadditivity and strong subadditivity, thereby establishing it as a viable candidate for the entropy associated to the subsystem. Finally, we demonstrate how our model generalizes the entropy vectors obtainable via conventional flows in undirected graphs and hypergraphs.

The adiabatic hydrodynamization framework is a promising framework within which to describe and characterize pre-hydrodynamic attractors in a model-independent fashion. Using this framework, we define a procedure to identify a time-dependent change in coordinates which reveals a dynamical reduction in the number of active degrees of freedom. Applying this procedure to the kinetic theory of a Bjorken-expanding gas of gluons in the small angle elastic scattering limit, we are able to intuitively explain the self-similar evolution of the gluon distribution function long before the applicability of hydrodynamics, as well as the loss of memory of its initial condition.

The rapid growth in size of quantum devices demands efficient ways to control them, which is challenging for systems with thousands of qubits or more. Here, we present a simple yet powerful solution: robust, site-dependent control of an arbitrary number of qubits in parallel with only minimal local tunability of the driving field. Inspired by recent experimental advances, we consider access to only one of three constrained local control capabilities: local control of either the phase or amplitude of the beam at each qubit, or individual Z rotations. In each case, we devise parallelizable composite pulse sequences to realize arbitrary single-qubit unitaries robust against quasistatic amplitude and frequency fluctuations. Numerical demonstration shows that our approach outperforms existing sequences such as BB1 and CORPSE in almost all regimes considered, achieving average fidelity \(>0.999\) under a decoherence rate \(\sim 10^{-5}\), even with a few percent amplitude and frequency error. Our results indicate that even for very large qubit ensembles, accurate, individual manipulation can be achieved despite substantial control inhomogeneity.

It is difficult to construct a post-inflation QCD axion model that solves the axion quality problem (and hence the Strong CP problem) without introducing a cosmological disaster. In a post-inflation axion model, the axion field value is randomized during the Peccei-Quinn phase transition, and axion domain walls form at the QCD phase transition. We emphasize that the gauge equivalence of all minima of the axion potential (i.e., domain wall number one) is insufficient to solve the cosmological domain wall problem. The axion string on which a domain wall ends must exist as an individual object (as opposed to a multi-string state), and it must be produced in the early universe. These conditions are often not satisfied in concrete models. Post-inflation axion models also face a potential problem from fractionally charged relics; solving this problem often leads to low-energy Landau poles for Standard Model gauge couplings, reintroducing the quality problem. We study several examples, finding that models that solve the quality problem face cosmological problems, and vice versa. This is not a no-go theorem; nonetheless, we argue that it is much more difficult than generally appreciated to find a viable post-inflation QCD axion model. Successful examples may have a nonstandard cosmological history (e.g., multiple types of cosmic axion strings of different tensions), undermining the widespread expectation that the post-inflation QCD axion scenario predicts a unique mass for axion dark matter.

To maximize the discovery potential of high-energy colliders, experimental searches should be sensitive to unforeseen new physics scenarios. This goal has motivated the use of machine learning for unsupervised anomaly detection. In this paper, we introduce a new anomaly detection strategy called FORCE: factorized observables for regressing conditional expectations. Our approach is based on the inductive bias of factorization, which is the idea that the physics governing different energy scales can be treated as approximately independent. Assuming factorization holds separately for signal and background processes, the appearance of non-trivial correlations between low- and high-energy observables is a robust indicator of new physics. Under the most restrictive form of factorization, a machine-learned model trained to identify such correlations will in fact converge to the optimal new physics classifier. We test FORCE on a benchmark anomaly detection task for the Large Hadron Collider involving collimated sprays of particles called jets. By teasing out correlations between the kinematics and substructure of jets, our method can reliably extract percent-level signal fractions. This strategy for uncovering new physics adds to the growing toolbox of anomaly detection methods for collider physics with a complementary set of assumptions.

We introduce a diffusion-based generative model to describe the distribution of galaxies in our Universe directly as a collection of points in 3-D space (coordinates) optionally with associated attributes (e.g., velocities and masses), without resorting to binning or voxelization. The custom diffusion model can be used both for emulation, reproducing essential summary statistics of the galaxy distribution, as well as inference, by computing the conditional likelihood of a galaxy field. We demonstrate a first application to massive dark matter haloes in the Quijote simulation suite. This approach can be extended to enable a comprehensive analysis of cosmological data, circumventing limitations inherent to summary statistic — as well as neural simulation-based inference methods.

We investigate Wilsonian effective field theory as a model for the construction of the tachyon potential and nonperturbative vacua in closed string field theory. In a number of cases we are able to find the effective potential exactly, and observe what appear to be universal features. We find that the effective field theory contains the same nonperturbative vacuum structure as the bare Lagrangian, though this information is encoded less efficiently as the distance scale of the effective field theory is increased. The implication is that closed string field theory plausibly contains information about the nonperturbative vacuum structure of string theory, in spite of its similarities to effective field theory. We also truncate the effective potential at a fixed power of the field and investigate how the global structure of the effective potential may be approximated via Padé resummation. Qualitative comparisons suggest that computation of the eighth to sixteenth order closed string vertex should be enough to obtain reliable results for the closed string field theory action evaluated on the tachyon field.

We investigate the transverse energy-energy correlators (TEEC) in the small-\(x\) regime at the upcoming Electron-Ion Collider (EIC). Focusing on the back-to-back production of electron-hadron pairs in both \(ep\) and \(eA\) collisions, we establish a factorization theorem given in terms of the hard function, quark distributions, soft functions, and TEEC jet functions, where the gluon saturation effect is incorporated. Numerical results for TEEC in both \(ep\) and \(eA\) collisions are presented, together with the nuclear modification factor \(R_A\). Our analysis reveals that TEEC observables in deep inelastic scattering provide a valuable approach for probing gluon saturation phenomena. Our findings underscore the significance of measuring TEEC at the EIC, emphasizing its efficacy in advancing our understanding of gluon saturation and nuclear modifications in high-energy collisions.

Accelerator-based experiments reliant on charged pion and kaon decays to produce muon-neutrino beams also deliver an associated powerful flux of muons. Therefore, these experiments can additionally be sensitive to light new particles that preferentially couple to muons and decay to visible final states on macroscopic length scales. Such particles are produced through rare 3-body meson decays in the decay pipe or via muon scattering in the beam dump, and decay in a downstream detector. To demonstrate the potential of this search strategy, we recast existing MiniBooNE and MicroBooNE studies of neutral pion production in neutrino-induced neutral-current scattering (\(\nu_\mu N \to \nu_\mu N \pi^0,~\pi^0\rightarrow \gamma\gamma\)) to place new leading limits on light (\(< 2m_\mu\)) muon-philic scalar particles that decay to diphotons through loops of virtual muons. Our results exclude scalars of mass between 10 and 60 MeV in which this scenario resolves the muon \(g-2\) anomaly. We also make projections for the sensitivity of SBND to these models and provide a road map for future neutrino experiments to perform dedicated searches for muon-philic forces.

Spacetime inversion symmetries such as parity and time reversal play a central role in physics, but they are usually treated as global symmetries. In quantum gravity there are no global symmetries, so any spacetime inversion symmetries must be gauge symmetries. In particular this includes \(\mathcal{CRT}\) symmetry (in even dimensions usually combined with a rotation to become \(\mathcal{CPT}\)), which in quantum field theory is always a symmetry and seems likely to be a symmetry of quantum gravity as well. In this article we discuss what it means to gauge a spacetime inversion symmetry, and we explain some of the more unusual consequences of doing this. In particular we argue that the gauging of \(\mathcal{CRT}\) is automatically implemented by the sum over topologies in the Euclidean gravity path integral, that in a closed universe the Hilbert space of quantum gravity must be a real vector space, and that in Lorentzian signature manifolds which are not time-orientable must be included as valid configurations of the theory. In particular we give an example of an asymptotically-AdS time-unorientable geometry which must be included to reproduce computable results in the dual CFT.

In gravity, spacelike separated regions can be dependent on each other due to the constraint equations. In this paper, we give a natural definition of subsystem independence and gravitational dressing of perturbations in classical gravity. We find that extremal surfaces, non-perturbative lumps of matter, and generic trapped surfaces are structures that enable dressing and subregion independence. This leads to a simple intuitive picture for why extremal surfaces tend to separate independent subsystems. The underlying reason is that localized perturbations on one side of an extremal surface contribute negatively to the mass on the other side, making the gravitational constraints behave as if there exist both negative and positive charges. Our results support the consistency of islands in massless gravity, shed light on the Python’s lunch, and provide hints on the nature of the split property in perturbatively quantized general relativity. We also prove a theorem bounding the area of certain surfaces in spherically symmetric asymptotically de Sitter spacetimes from above and below in terms of the horizon areas of de Sitter and Nariai. This theorem implies that it is impossible to deform a single static patch without also deforming the opposite patch, provided we assume spherical symmetry and an energy condition.