Center for Theoretical Physics Publications

The Hamiltonian limit of lattice gauge theories can be found by extrapolating the results of anisotropic lattice computations, i.e., computations using lattice actions with different temporal and spatial lattice spacings (\(a_t\neq a_s\)), to the limit of \(a_t\to 0\). In this work, we present a study of this Hamiltonian limit for a Euclidean \(U(1)\) gauge theory in 2+1 dimensions (QED3), regularized on a toroidal lattice. The limit is found using the renormalized anisotropy \(\xi_R=a_t/a_s\), by sending \(\xi_R \to 0\) while keeping the spatial lattice spacing constant. We compute \(\xi_R\) in \(3\) different ways: using both the ``normal'' and the ``sideways'' static quark potential, as well as the gradient flow evolution of gauge fields. The latter approach will be particularly relevant for future investigations of combining quantum computations with classical Monte Carlo computations, which requires the matching of lattice results obtained in the Hamiltonian and Lagrangian formalisms.

The Heavy Quark Expansion (HQE) gives an expansion of inclusive decay rates of \(b\)-hadrons as a simultaneous series in \(\alpha_s\) and \(1/m_b\), in terms of perturbatively defined coefficients and non-perturbative matrix elements. Spectator effects arise from the dimension-\(6\) operators in the HQE in which along with the heavy quark, a light spectator quark from the hadron participates in the weak decay. In this work, we provide a lattice determination of the bare matrix elements that contribute to the spectator effects for the \(B^+\)-meson, \(B_d\)-meson and the \(\Lambda_b\)-baryon. The computations were performed on two separate lattice spacings of the RBC-UKQCD (2+1)-flavor Domain Wall Fermion ensembles. Renormalization and continuum extrapolation of the matrix elements have not yet been computed, and are left to future work.

Dark matter (DM) in the Milky Way halo may annihilate or decay to photons, producing monochromatic gamma rays. We search for DM-induced spectral lines using 14 years of data from the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope (\(\textit{Fermi}\)-LAT) between \(10\,\mathrm{GeV}\) and \(2\,\mathrm{TeV}\) in the inner Milky Way leveraging both the spatial and spectral morphology of an expected signal. We present new constraints as strong as \(\langle \sigma v \rangle \lesssim 6\times 10^{-30}\, \mathrm{cm}^3/\mathrm{s}\) for the two-to-two annihilations and \(\tau \gtrsim 10^{30}\,\mathrm{s}\) for one-to-two decays, representing leading sensitivity between \(10\,\mathrm{GeV}\) and \(\sim\)\(500\,\mathrm{GeV}\). We consider the implications of our line-constraints on the Galactic Center Excess (GCE), which is a previously-observed excess of continuum \(\sim\)GeV gamma-rays that may be explained by DM annihilation. The Higgs portal and neutralino-like DM scenarios, which have been extensively discussed as possible origins of the GCE, are constrained by our work because of the lack of observed one-loop decays to two photons. More generally, we interpret our null results in a variety of annihilating and decaying DM models, such as neutralinos, gravitinos, and glueballs, showing that in many cases the line search is more powerful than the continuum, despite the continuum annihilation being at tree level.

We present a neural flow wavefunction, Gauge-Fermion FlowNet, and use it to simulate 2+1D lattice compact quantum electrodynamics with finite density dynamical fermions. The gauge field is represented by a neural network which parameterizes a discretized flow-based transformation of the amplitude while the fermionic sign structure is represented by a neural net backflow. This approach directly represents the \(U(1)\) degree of freedom without any truncation, obeys Guass's law by construction, samples autoregressively avoiding any equilibration time, and variationally simulates Gauge-Fermion systems with sign problems accurately. In this model, we investigate confinement and string breaking phenomena in different fermion density and hopping regimes. We study the phase transition from the charge crystal phase to the vacuum phase at zero density, and observe the phase seperation and the net charge penetration blocking effect under magnetic interaction at finite density. In addition, we investigate a magnetic phase transition due to the competition effect between the kinetic energy of fermions and the magnetic energy of the gauge field. With our method, we further note potential differences on the order of the phase transitions between a continuous \(U(1)\) system and one with finite truncation. Our state-of-the-art neural network approach opens up new possibilities to study different gauge theories coupled to dynamical matter in higher dimensions.

This paper is essentially a verbatim reconstruction of lectures that I gave at the Los Alamos School on Relativistic Dynamics and Quark Nuclear Physics in 1985. They were published in the school proceedings, but the book is not widely available. The Los Alamos School took place at the height of the first wave of interest in the quark substructure of nuclei, stimulated by the 1983 discovery of the EMC Effect. Interest in this subject has been increasing for years and the prospect of a dedicated Electron Ion Collider within the decade guarantees even greater attention to quarks and gluons in nuclei among both theorists and experimentalists. Recently, to my surprise, I learned that copies of my old lectures have been circulating and been found useful by the relatively few people who know about them. The are, of course, dated: experiments have far outstripped what was available 37 years ago and theory has progressed too. However, the rest frame derivation of the parton model, the derivation and discussion of the convolution formalism for nucleons, nucleon correlations, and other, virtual, constituents of nuclei, and sections on scaling violation and the operator product expansion have aged pretty well and seem to still be useful. With the help and encouragement of Richard Milner, I have recreated the LaTeX files necessary to post the 1985 Lectures on the arXiv, making them available to the nuclear and particle physics community. Apart from correcting some typographical errors, I have made no attempt to edit, improve, or update these lectures. I hope readers will nevertheless find them useful.

In axion-Maxwell theory at the minimal axion-photon coupling, we find non-invertible 0- and 1-form global symmetries arising from the naive shift and center symmetries. Since the Gauss law is anomalous, there is no conserved, gauge-invariant, and quantized electric charge. Rather, using half higher gauging, we find a non-invertible Gauss law associated with a non-invertible 1-form global symmetry, which is related to the Page charge. These symmetries act invertibly on the axion field and Wilson line, but non-invertibly on the monopoles and axion strings, leading to selection rules related to the Witten effect. We also derive various crossing relations between the defects. Finally, we show that the non-invertible 0- and 1-form global symmetries mix with other invertible symmetries in a way reminiscent of a higher-group symmetry. Using this non-invertible higher symmetry structure, we derive universal inequalities on the energy scales where different infrared symmetries emerge in any renormalization group flow to the axion-Maxwell theory.

A central challenge in analog quantum simulation is to characterize desirable physical properties of quantum states produced in experiments. However, in conventional approaches, the extraction of arbitrary information requires performing measurements in many different bases, which necessitates a high level of control that present-day quantum devices may not have. Here, we propose and analyze a scalable protocol that leverages the ergodic nature of generic quantum dynamics, enabling the efficient extraction of many physical properties. The protocol does not require sophisticated controls and can be generically implemented in analog quantum simulation platforms today. Our protocol involves introducing ancillary degrees of freedom in a predetermined state to a system of interest, quenching the joint system under Hamiltonian dynamics native to the particular experimental platform, and then measuring globally in a single, fixed basis. We show that arbitrary information of the original quantum state is contained within such measurement data, and can be extracted using a classical data-processing procedure. We numerically demonstrate our approach with a number of examples, including the measurements of entanglement entropy, many-body Chern number, and various superconducting orders in systems of neutral atom arrays, bosonic and fermionic particles on optical lattices, respectively, only assuming existing technological capabilities. Our protocol excitingly promises to overcome limited controllability and, thus, enhance the versatility and utility of near-term quantum technologies.

Physicists dating back to Feynman have lamented the difficulties of applying the variational principle to quantum field theories. In non-relativistic quantum field theories, the challenge is to parameterize and optimize over the infinitely many \(n\)-particle wave functions comprising the state's Fock space representation. Here we approach this problem by introducing neural-network quantum field states, a deep learning ansatz that enables application of the variational principle to non-relativistic quantum field theories in the continuum. Our ansatz uses the Deep Sets neural network architecture to simultaneously parameterize all of the \(n\)-particle wave functions comprising a quantum field state. We employ our ansatz to approximate ground states of various field theories, including an inhomogeneous system and a system with long-range interactions, thus demonstrating a powerful new tool for probing quantum field theories.

The presence of symmetries can convert an otherwise computationally hard task into an efficient one. In light of recently proposed quantum algorithms that incorporate symmetries in the hope of quantum advantage, we show that with symmetries that are restrictive enough, classical algorithms can efficiently emulate their quantum counterparts. As an example, we show that simulating and finding the ground state of permutation invariant Hamiltonians can be efficiently performed in polynomial classical time, further restricting the possible set of symmetries that provide quantum advantage.

The pion light-cone distribution amplitude (LCDA) is a central non-perturbative object of interest for the calculation of high-energy exclusive processes in quantum chromodynamics. This article describes the progress in the lattice QCD calculation of the fourth Mellin moment of the pion LCDA using a heavy-quark operator product expansion (HOPE).

We propose a variational quantum eigensolver suitable for exploring the phase structure of the multi-flavor Schwinger model in the presence of a chemical potential. The parametric ansatz circuit we design is capable of incorporating the symmetries of the model, present in certain parameter regimes, which allows for reducing the number of variational parameters substantially. Moreover, the ansatz circuit can be implementated on both measurement-based and circuit-based quantum hardware. We numerically demonstrate that our ansatz circuit is able to capture the phase structure of the model and allows for faithfully approximating the ground state. Our results show that our approach is suitable for current intermediate-scale quantum hardware and can be readily implemented on existing quantum devices.

Lattice computations in the Hamiltonian formulation have so far mainly focused on staggered fermions. In these proceedings, we study Wilson fermions in the Hamiltonian formulation and propose a new method to determine the resulting mass shift. As a benchmark study, we examine the one-flavour Schwinger model with Wilson fermions and a topological \(\theta\)-term using matrix product states. Wilson fermions explicitly break chiral symmetry; thus, the bare mass of the lattice model receives an additive renormalization. In order to measure this mass shift directly, we develop a method that is suitable for the Hamiltonian formulation, which relies on the fact that the vacuum expectation value of the electric field density vanishes when the renormalized mass is zero. We examine the dependence of the mass shift on the lattice spacing, the lattice volume, the \(\theta\)-parameter, and the Wilson parameter. Using the mass shift, we then perform the continuum extrapolation of the electric field density and compare the resulting mass dependence to the analytical predictions of mass perturbation theory. We demonstrate that incorporating the mass shift significantly improves the continuum extrapolation. Finally, we apply our method to the same model using staggered fermions instead of Wilson fermions and compare the resulting mass shift to recent theoretical predictions.

We study the massless charged spectrum of U(1) gauge fields in F-theory that arise from flux breaking of a nonabelian group. The U(1) charges that arise in this way can be very large. In particular, using vertical flux breaking, we construct an explicit 4D F-theory model with a U(1) decoupled from other gauge sectors, in which the massless/light fields have charges as large as 657. This result greatly exceeds prior results in the literature. We argue heuristically that this result may provide an upper bound on charges for light fields under decoupled U(1) factors in the F-theory landscape. We also show that the charges can be even larger when the U(1) is coupled to other gauge groups.

The infinite-component Chern-Simons-Maxwell (iCSM) theory is a 3+1D generalization of the 2+1D Chern-Simons-Maxwell theory by including an infinite number of coupled gauge fields. It can be used to describe interesting 3+1D systems. In Phys. Rev. B 105, 195124 (2022), it was used to construct gapped fracton models both within and beyond the foliation framework. In this paper, we study the nontrivial features of gapless iCSM theories. In particular, we find that while gapless 2+1D Maxwell theories are confined and not robust due to monopole effect, gapless iCSM theories are deconfined and robust against all local perturbation and hence represent a robust 3+1D deconfined gapless order. The gaplessness of the gapless iCSM theory can be understood as a consequence of the spontaneous breaking of an exotic one-form symmetry. Moreover, for a subclass of the gapless iCSM theories, we find interesting topological features in the correlation and response of the system. Finally, for this subclass of theories, we propose a fully continuous field theory description of the model that captures all these features.

We explore the relationship between renormalization group (RG) flow and error correction by constructing quantum algorithms that exactly recognize 1D symmetry-protected topological (SPT) phases protected by finite internal Abelian symmetries. For each SPT phase, our algorithm runs a quantum circuit which emulates RG flow: an arbitrary input ground state wavefunction in the phase is mapped to a unique minimally-entangled reference state, thereby allowing for efficient phase identification. This construction is enabled by viewing a generic input state in the phase as a collection of coherent `errors' applied to the reference state, and engineering a quantum circuit to efficiently detect and correct such errors. Importantly, the error correction threshold is proven to coincide exactly with the phase boundary. We discuss the implications of our results in the context of condensed matter physics, machine learning, and near-term quantum algorithms.

The geometry of 4-string contact interaction of closed string field theory is characterized using machine learning. We obtain Strebel quadratic differentials on 4-punctured spheres as a neural network by performing unsupervised learning with a custom-built loss function. This allows us to solve for local coordinates and compute their associated mapping radii numerically. We also train a neural network distinguishing vertex from Feynman region. As a check, 4-tachyon contact term in the tachyon potential is computed and a good agreement with the results in the literature is observed. We argue that our algorithm is manifestly independent of number of punctures and scaling it to characterize the geometry of \(n\)-string contact interaction is feasible.

We prove novel speed limits on the growth of entanglement, equal time correlators, and spacelike Wilson loops in spatially uniform time-evolving states in strongly coupled CFTs with holographic duals. These bounds can also be viewed as quantum weak energy conditions. Several of the speed limits are valid for regions of arbitrary size and with multiple connected components, and our findings imply new bounds on the effective entanglement velocity of small subregions. In 2d CFT, our results prove a conjecture by Liu and Suh for a large class of states. Key to our findings is a momentum-entanglement correspondence, showing that entanglement growth is computed by the momentum crossing the HRT surface. In our setup, we prove a number of general features of boundary-anchored extremal surfaces, such as a sharp bound on the smallest radius that a surface can probe, and that the tips of extremal surfaces cannot lie in trapped regions. Our methods rely on novel global GR techniques, including a delicate interplay between Lorentzian and Riemannian Hawking masses. While our proofs assume the dominant energy condition in the bulk, we provide numerical evidence that our bounds are true under less restrictive assumptions.

Recent applications of machine-learned normalizing flows to sampling in lattice field theory suggest that such methods may be able to mitigate critical slowing down and topological freezing. However, these demonstrations have been at the scale of toy models, and it remains to be determined whether they can be applied to state-of-the-art lattice quantum chromodynamics calculations. Assessing the viability of sampling algorithms for lattice field theory at scale has traditionally been accomplished using simple cost scaling laws, but as we discuss in this work, their utility is limited for flow-based approaches. We conclude that flow-based approaches to sampling are better thought of as a broad family of algorithms with different scaling properties, and that scalability must be assessed experimentally.

We study the properties of three-body resonances using a lattice complex scalar \(\varphi^4\) theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy-plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory.

Gauge Theory plays a crucial role in many areas in science, including high energy physics, condensed matter physics and quantum information science. In quantum simulations of lattice gauge theory, an important step is to construct a wave function that obeys gauge symmetry. In this paper, we have developed gauge equivariant neural network wave function techniques for simulating continuous-variable quantum lattice gauge theories in the Hamiltonian formulation. We have applied the gauge equivariant neural network approach to find the ground state of 2+1-dimensional lattice gauge theory with U(1) gauge group using variational Monte Carlo. We have benchmarked our approach against the state-of-the-art complex Gaussian wave functions, demonstrating improved performance in the strong coupling regime and comparable results in the weak coupling regime.