MIT Center for Theoretical Physics — a Leinweber Institute Publications

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked (CFL) quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory (CEFT) to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J0348+0432, PSR J1624-2230, PSR J0740+6620, GW170817). At baryon number chemical potential \(\mu = 2.6~\text{GeV}\) we find a 95% upper limit on the CFL pairing gap \(\Delta\) of \(457~\text{MeV}\) using overly conservative assumptions and \(216~\text{MeV}\) with more reasonable assumptions. This constraint may be strengthened by future astrophysical measurements as well as by future advances in high density QCD calculations.

The three-particle \(K\)-matrix, \(\mathcal{K}_{\mathrm{df},3}\), is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three-particle finite-volume formalism. In this work, we compute its value for systems of three pions in all isospin channels through next-to-leading order in Chiral Perturbation Theory, generalizing previous work done at maximum isospin. We obtain analytic expressions through quadratic order (or cubic order, in the case of zero isospin) in the expansion about the three-pion threshold.

Machine-learned normalizing flows can be used in the context of lattice quantum field theory to generate statistically correlated ensembles of lattice gauge fields at different action parameters. This work demonstrates how these correlations can be exploited for variance reduction in the computation of observables. Three different proof-of-concept applications are demonstrated using a novel residual flow architecture: continuum limits of gauge theories, the mass dependence of QCD observables, and hadronic matrix elements based on the Feynman-Hellmann approach. In all three cases, it is shown that statistical uncertainties are significantly reduced when machine-learned flows are incorporated as compared with the same calculations performed with uncorrelated ensembles or direct reweighting.

We present a framework to compute amplitudes for the gravitational analog of the Raman process, a quasi-elastic scattering of waves off compact objects, in worldline effective field theory (EFT). As an example, we calculate third post-Minkowskian (PM) order (\(\mathcal{O}(G^3)\)), or two-loop, phase shifts for the scattering of a massless scalar field including all tidal effects and dissipation. Our calculation unveils two sources of the classical renormalization-group flow of dynamical Love numbers: a universal running independent of the nature of the compact object, and a running self-induced by tides. Restricting to the black hole case, we find that our EFT phase shifts agree exactly with those from general relativity, provided that the relevant static Love numbers are set to zero. In addition, we carry out a complete matching of the leading scalar dynamical Love number required to renormalize a universal short scale divergence in the S-wave. Our results pave the way for systematic calculations of gravitational Raman scattering at higher PM orders.

We report on the ongoing effort of improving the determination of the gradient flow scale on the (2+1+1)-flavor HISQ ensembles generated by the MILC collaboration. We compute the scales \(\sqrt{t_0}/a\) and \(w_0/a\) with the Wilson and Symanzik flow using three discretizations for the action density: clover, Wilson and tree-level Symanzik improved. For the absolute scale setting, we intend to employ the \(\Omega\)-baryon mass, but are also using the pion decay constant while the \(\Omega\)-mass calculations are in progress.

We generalize the relativistic field-theoretic three-particle finite-volume scattering formalism to describe generic \(DD\pi\) systems in the charm \(C=2\) sector. This includes the isospin-0 channel, in which the recently discovered doubly-charmed tetraquark \(T_{cc}(3875)^+\) is expected to manifest as a pole in the \(DD \pi \to DD \pi\) scattering amplitude. The formalism presented here can also be applied to lattice QCD settings in which the \(D^*\) is bound and, in particular, remains valid below the left-hand cut in \(D D^*\) scattering, thus resolving an issue in previous analyses of lattice-determined finite-volume energies.

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.