MIT Center for Theoretical Physics — a Leinweber Institute Publications

We review briefly the characteristic topological data of Calabi--Yau threefolds and focus on the question of when two threefolds are equivalent through related topological data. This provides an interesting test case for machine learning methodology in discrete mathematics problems motivated by physics.

The LUXE experiment (LASER Und XFEL Experiment) is a new experiment in planning at DESY Hamburg, which will study Quantum Electrodynamics (QED) at the strong-field frontier. In this regime, QED is non-perturbative. This manifests itself in the creation of physical electron-positron pairs from the QED vacuum. LUXE intends to measure the positron production rate in this unprecedented regime by using, among others, a silicon tracking detector. The large number of expected positrons traversing the sensitive detector layers results in an extremely challenging combinatorial problem, which can become computationally very hard for classical computers. This paper presents a preliminary study to explore the potential of quantum computers to solve this problem and to reconstruct the positron trajectories from the detector energy deposits. The reconstruction problem is formulated in terms of a quadratic unconstrained binary optimisation. Finally, the results from the quantum simulations are discussed and compared with traditional classical track reconstruction algorithms.

There is great potential to apply machine learning in the area of numerical lattice quantum field theory, but full exploitation of that potential will require new strategies. In this white paper for the Snowmass community planning process, we discuss the unique requirements of machine learning for lattice quantum field theory research and outline what is needed to enable exploration and deployment of this approach in the future.

Finite-volume pionless effective field theory provides an efficient framework for the extrapolation of nuclear spectra and matrix elements calculated at finite volume in lattice QCD to infinite volume, and to nuclei with larger atomic number. In this work, it is demonstrated how this framework may be implemented via a set of correlated Gaussian wavefunctions optimised using differentiable programming and via solution of a generalised eigenvalue problem. This approach is shown to be significantly more efficient than a stochastic implementation of the variational method based on the same form of correlated Gaussian wavefunctions, yielding comparably accurate representations of the ground-state wavefunctions with an order of magnitude fewer terms. The efficiency of representation allows such calculations to be extended to larger systems than in previous work. The method is demonstrated through calculations of the binding energies of nuclei with atomic number \(A\in\{2,3,4\}\) in finite volume, matched to lattice QCD calculations at quark masses corresponding to \(m_\pi=806\) MeV, and infinite-volume effective field theory calculations of \(A\in\{2,3,4,5,6\}\) systems based on this matching.

We present the first lattice study of pion-pion scattering with varying number of colors, \(N_\text{c}\). We use lattice simulations with four degenerate quark flavors, \(N_\text{f}=4\), and \(N_\text{c}=3-6\). We focus on two scattering channels that do not involve vacuum diagrams. These correspond to two irreducible representations of the SU(4) flavor group: the fully symmetric one, \(SS\), and the fully antisymmetric one, \(AA\). The former is a repulsive channel equivalent to the isospin-2 channel of SU(2). By contrast, the latter is attractive and only exists for \(N_\text{f} \geq 4\). A representative state is \(\left( |D_s^+ \pi^+\rangle - |D^+ K^+\rangle\right)/\sqrt{2}\). Using Lüscher's formalism, we extract the near-threshold scattering amplitude and we match our results to Chiral Perturbation Theory (ChPT) at large \(N_\text{c}\). For this, we compute the analytical U\((N_\text{f})\) ChPT prediction for two-pion scattering, and use the lattice results to constrain the \(N_\text{c}\) scaling of the relevant low-energy couplings.

In a chain of \(N\) spins governed by a local Hamiltonian, we consider eigenstates in a microcanonical ensemble in the middle of the energy spectrum. We prove that if the bandwidth of the ensemble is \(C\ln^2N\) for a large constant \(C\), then the average entanglement entropy of these eigenstates is bounded away from the maximum entropy. This result highlights the difference between the entanglement of mid-spectrum eigenstates of (chaotic) local Hamiltonians and that of random states, for the latter is nearly maximal. We also prove that the former is different from the thermodynamic entropy at the same energy.

Transverse-momentum-dependent parton distribution functions (TMDs) can be studied from first principles by a perturbative matching onto lattice-calculable quantities: so-called lattice TMDs, which are a class of equal-time correlators that includes quasi-TMDs and TMDs in the Lorentz-invariant approach. We introduce a general correlator that includes as special cases these two Lattice TMDs and continuum TMDs, like the Collins scheme. Then, to facilitate the derivation of a factorization relation between lattice and continuum TMDs, we construct a new scheme, the Large Rapidity (LR) scheme, intermediate between the Collins and quasi-TMDs. The LR and Collins schemes differ only by an order of limits, and can be matched onto one another by a multiplicative kernel. We show that this same matching also holds between quasi and Collins TMDs, which enables us to prove a factorization relation between these quantities to all orders in \(\alpha_s\). Our results imply that there is no mixing between various quark flavors or gluons when matching Collins and quasi TMDs, making the lattice calculation of individual flavors and gluon TMDs easier than anticipated. We cross-check these results explicitly at one loop and discuss implications for other physical-to-lattice scheme factorizations.

We show that the canonical purification of an evaporating black hole after the Page time consists of a short, connected, Lorentzian wormhole between two asymptotic boundaries, one of which is unitarily related to the radiation. This provides a quantitative and general realization of the predictions of ER=EPR in an evaporating black hole after the Page time; this further gives a standard AdS/CFT calculation of the entropy of the radiation (without modifications of the homology constraint). Before the Page time, the canonical purification consists of two disconnected, semiclassical black holes. From this, we construct two bipartite entangled holographic CFT states, with equal (and large) amount of entanglement, where the semiclassical dual of one has a connected ERB and the other does not. From this example, we speculate that that measures of multipartite entanglement may offer a more complete picture into the emergence of spacetime.

Jets of hadrons produced at high-energy colliders provide experimental access to the dynamics of asymptotically free quarks and gluons and their confinement into hadrons. In this paper, we show that the high energies of the Large Hadron Collider (LHC), together with the exceptional resolution of its detectors, allow multipoint correlation functions of energy flow operators to be directly measured within jets for the first time. Using Open Data from the CMS experiment, we show that reformulating jet substructure in terms of these correlators provides new ways of probing the dynamics of QCD jets, which enables direct imaging of the confining transition to free hadrons as well as precision measurements of the scaling properties and interactions of quarks and gluons. This opens a new era in our understanding of jet substructure and illustrates the immense unexploited potential of high-quality LHC data sets for elucidating the dynamics of QCD.

Quarkonium suppression in relativistic heavy ion collisions has been studied experimentally for decades to probe the properties of the quark-gluon plasma. For this purpose, complete theoretical understanding of the time evolution of quarkonium inside the quark-gluon plasma is needed but challenging. Here I review recent progress in applying the open quantum system framework to describe the real-time dynamics of quarkonium, with a focus on the gauge-invariant chromoelectric field correlators of the plasma that control the dynamics.

The extraction of nonperturbative TMD physics is made challenging by prescriptions that shield the Landau pole, which entangle long- and short-distance contributions in momentum space. The use of different prescriptions then makes the comparison of fit results for underlying nonperturbative contributions not meaningful on their own. We propose a model-independent method to restrict momentum-space observables to the perturbative domain. This method is based on a set of integral functionals that act linearly on terms in the conventional position-space operator product expansion (OPE). Artifacts from the truncation of the integral can be systematically pushed to higher powers in \(\Lambda_{\rm QCD}/k_T\). We demonstrate that this method can be used to compute the cumulative integral of TMD PDFs over \(k_T \le k_T^\mathrm{cut}\) in terms of collinear PDFs, accounting for both radiative corrections and evolution effects. This yields a systematic way of correcting the naive picture where the TMD PDF integrates to a collinear PDF, and for unpolarized quark distributions we find that when renormalization scales are chosen near \(k_T^\mathrm{cut}\), such corrections are a percent-level effect. We also show that, when supplemented with experimental data and improved perturbative inputs, our integral functionals will enable model-independent limits to be put on the nonperturbative OPE contributions to the Collins-Soper kernel and intrinsic TMD distributions.

We study non-perturbative superpotential generated by D(-1)-branes in type IIB compactifications on orientifolds of Calabi-Yau threefold hypersurfaces. To compute D-instanton superpotential, we study F-theory compactification on toric complete intersection elliptic Calabi-Yau fourfolds. We take the Sen-limit, but with finite \(g_s,\) in F-theory compactification with a restriction that all D7-branes are carrying SO(8) gauge groups, which we call the global Sen-limit. In the global Sen-limit, the axio-dilaton is not varying in the compactification manifold. We compute the Picard-Fuchs equations of elliptic Calabi-Yau fourfolds in the global Sen-limit, and show that the Picard-Fuchs equations of the elliptic fourfolds split into that of the underlying Calabi-Yau threefolds and of the elliptic fiber. We then demonstrate that this splitting property of the Picard-Fuchs equation implies that the fourform period of the elliptic Calabi-Yau fourfolds in the global Sen-limit does not contain exponentially suppressed terms \(\mathcal{O}(e^{-\pi/g_s})\). With this result, we finally show that in the global Sen-limit, superpotential of the underlying type IIB compactification does not receive D-instanton contributions.

The computational cost required to calculate nuclear correlation functions grows factorially in the number of quarks, making the study of large nuclei inaccessible to ab initio study using lattice QCD at the present time. However, the tensor expressions corresponding to many of these correlation functions exhibit a high degree of permutation symmetry that can be exploited to reduce computational work. We present promising speed-ups for certain choices of interpolating operators using two new algorithms for computing nuclear correlation functions.

Tracking the dynamics of physical systems in real time is a prime application of digital quantum computers. Using a trapped-ion system with up to six qubits, we simulate the real-time dynamics of a lattice gauge theory in 1+1 dimensions, i.e., the lattice Schwinger model, and demonstrate non-perturbative effects such as pair creation for times much longer than previously accessible. We study the gate requirement of two formulations of the model using the Suzuki-Trotter product formula, as well as the trade-off between errors from the ordering of the Hamiltonian terms, the Trotter step size, and experimental imperfections. To mitigate experimental errors, a recent symmetry-protection protocol for suppressing coherent errors and a symmetry-inspired post-selection scheme are applied. This work demonstrates the integrated theoretical, algorithmic, and experimental approach that is essential for efficient simulation of lattice gauge theories and other complex physical systems.

Nuclei with equal number of baryons but varying proton number (isobars) have many commonalities, but differ in both electric charge and nuclear structure. Relativistic collisions of such isobars provide unique opportunities to study the variation of the magnetic field, provided the nuclear structure is well understood. In this Letter we simulate collisions using several state-of-the-art parametrizations of the \(^{96}_{40}\)Zr and \(^{96}_{44}\)Ru isobars and show that a comparison with the exciting STAR measurement arXiv:2109.00131 of ultrarelativistic collisions can uniquely identify the structure of both isobars. This not only provides an urgently needed understanding of the structure of the Zirconium and Ruthenium isobars, but also paves the way for more detailed studies of nuclear structure using relativistic heavy ion collisions.

In holographic duality an eternal AdS black hole is described by two copies of the boundary CFT in the thermal field double state. In this paper we provide explicit constructions in the boundary theory of infalling time evolutions which can take bulk observers behind the horizon. The constructions also help to illuminate the boundary emergence of the black hole horizons, the interiors, and the associated causal structure. A key element is the emergence, in the large \(N\) limit of the boundary theory, of a type III\(_1\) von Neumann algebraic structure from the type I boundary operator algebra and the half-sided modular translation structure associated with it.

Motivated by recent experimental realizations of exotic phases of matter on programmable quantum simulators, we carry out a comprehensive theoretical study of quantum phase transitions in a Rydberg atom array on a square lattice, with both open and periodic boundary conditions. In the bulk, we identify several first-order and continuous phase transitions by performing large-scale quantum Monte Carlo simulations and develop an analytical understanding of the nature of these transitions using the framework of Landau-Ginzburg-Wilson theory. Remarkably, we find that under open boundary conditions, the boundary itself undergoes a second-order quantum phase transition, independent of the bulk. These results explain recent experimental observations and provide important new insights into both the adiabatic state preparation of novel quantum phases and quantum optimization using Rydberg atom array platforms.

A new population of millisecond pulsars is a long-standing proposed explanation for the excess of GeV-scale gamma rays emanating from the region surrounding the center of the Milky Way (the "Galactic Center excess"). We examine several simple parameterizations of possible luminosity functions for this population, as well as several benchmark luminosity functions proposed in the literature, and compare the predicted populations of resolved point sources to the Fermi 4FGL-DR2 point source catalog and a sub-population recently identified using wavelet-based methods. We provide general results that can be used to translate upper limits on the number of resolved point sources associated with the excess, and the fraction of the flux in the excess that can be attributed to resolved sources, into limits on the luminosity function parameter space. We discuss a number of important systematic uncertainties, including in the detection threshold model and the total flux attributed to the excess. We delineate regions of parameter space (containing existing benchmark models) where there is no apparent tension with current data, and the number of total pulsars needed to explain the excess is in the range of O(10,000-100,000). Forecasting the effects of lowered point source detection thresholds, we show that novel analysis methods that probe sub-threshold point source populations can hope to resolve more than 30% of the flux of the excess.

Differential measurements of the semi-inclusive deep inelastic scattering (SIDIS) process with polarized beams provide important information on the three-dimensional structure of hadrons. Among the various observables are azimuthal asymmetries that start at subleading power, and which give access to novel transverse momentum dependent distributions (TMDs). Theoretical predictions for these distributions are currently based on the parton model rather than a rigorous factorization based analysis. Working under the assumption that leading power Glauber interactions do not spoil factorization at this order, we use the Soft Collinear Effective Theory to derive a complete factorization formula for power suppressed hard scattering effects in SIDIS. This yields generalized definitions of the TMDs that depend on two longitudinal momentum fractions (one of them only relevant beyond tree level), and a complete proof that only the same leading power soft function appears and can be absorbed into the TMD distributions at this order. We also show that perturbative corrections can be accounted for with only one new hard coefficient. Factorization formulae are given for all spin dependent structure functions which start at next-to-leading power. Prospects for improved subleading power predictions that include resummation are discussed.

We study the mass-deformed sphere free energy of three-dimensional \(\mathcal{N} = 2\) superconformal field theories with holographic duals. Building on previous observations, we conjecture a proportionality relation between the sphere free energy on the boundary and the prepotential of the four-dimensional \(\mathcal{N} = 2\) supergravity theory in the bulk. We verify this formula by explicit computation in several examples of supergravity theories with vector multiplets and hypermultiplets.