Lattice Field Theory Research Articles

Tensor renormalization group and the volume independence in 2D U(N) and SU(N) gauge theories

The tensor renormalization group method is a promising approach to lattice field theories, which is free from the sign problem unlike standard Monte Carlo methods. One of the remaining issues is the application to gauge theories, which is so far limited to U(1) and SU(2) gauge groups. In the case of higher rank, it becomes highly nontrivial to restrict the number of representations in the character expansion to be used in constructing the fundamental tensor. We propose a practical strategy to accomplish this and demonstrate it in 2D U(N) and SU(N) gauge theories, which are exactly solvable. Using this strategy, we obtain the singular-value spectrum of the fundamental tensor, which turns out to have a definite profile in the large-N limit. For the U(N) case, in particular, we show that the large-N behavior of the singular-value spectrum changes qualitatively at the critical coupling of the Gross-Witten-Wadia phase transition. As an interesting consequence, we find a new type of volume independence in the large-N limit of the 2D U(N) gauge theory with the θ term in the strong coupling phase, which goes beyond the Eguchi-Kawai reduction.

Lattice renormalization of quantum simulations

With advances in quantum computing, new opportunities arise to tackle challenging calculations in quantum field theory. We show that trotterized time-evolution operators can be related by analytic continuation to the Euclidean transfer matrix on an anisotropic lattice. In turn, trotterization entails renormalization of the temporal and spatial lattice spacings. Based on the tools of Euclidean lattice field theory, we propose two schemes to determine Minkowski lattice spacings, using Euclidean data and thereby overcoming the demands on quantum resources for scale setting. In addition, we advocate using a fixed-anisotropy approach to the continuum to reduce both circuit depth and number of independent simulations. We demonstrate these methods with Qiskit noiseless simulators for a $2+1$D discrete non-Abelian $D_4$ gauge theory with two spatial plaquettes.

Quantum algorithms for open lattice field theory

Certain aspects of some unitary quantum systems are well-described by evolution via a non-Hermitian effective Hamiltonian, as in the Wigner-Weisskopf theory for spontaneous decay. Conversely, any non-Hermitian Hamiltonian evolution can be accommodated in a corresponding unitary system + environment model via a generalization of Wigner-Weisskopf theory. This demonstrates the physical relevance of novel features such as exceptional points in quantum dynamics, and opens up avenues for studying many body systems in the complex plane of coupling constants. In the case of lattice field theory, sparsity lends these channels the promise of efficient simulation on standardized quantum hardware. We thus consider quantum operations that correspond to Suzuki-Lee-Trotter approximation of lattice field theories undergoing non-Hermitian time evolution, with potential applicability to studies of spin or gauge models at finite chemical potential, with topological terms, to quantum phase transitions - a range of models with sign problems. We develop non-Hermitian quantum circuits and explore their promise on a benchmark, the quantum one-dimensional Ising model with complex longitudinal magnetic field, showing that observables can probe the Lee-Yang edge singularity. The development of attractors past critical points in the space of complex couplings indicates a potential for study on near-term noisy hardware.

Rapid organic dye degradation and wavelength dependent sensing study in Cu1-xFexO

Cu1-xFexO (x = 0, 0.0156, 0.0234 and 0.0312) nanoparticles have been synthesized by using a simple sol-gel route. A pure monoclinic phase (space group c2/c) was observed from XRD studies. A Rietveld analysis of the XRD spectra using GSAS software helped evaluate the variations of lattice constant, lattice strain and bond length/angles. With Fe doping, while the indirect bandgap could be nominally tuned from ∼1.40 eV to ∼1.46 eV, the Fe-doping seems to improve carrier concentration, mobility, and lattice regularity than pure CuO. The photocatalytic degradation of toxic methylene blue (MB) dye could be enhanced from 0.91%/min for the pure CuO to 1.48%/min for 3.125%-Fe-doped CuO. This property has been correlated to the carrier properties of the materials. The UV and visible photoresponse of Cu1-xFexO nanoparticles was investigated. A drastic reduction in both UV and visible light-sensitivity with Fe-doping was observed and correlated with carrier properties.

Generalization capabilities of translationally equivariant neural networks

The rising adoption of machine learning in high energy physics and lattice field theory necessitates the re-evaluation of common methods that are widely used in computer vision, which, when applied to problems in physics, can lead to significant drawbacks in terms of performance and generalizability. One particular example for this is the use of neural network architectures that do not reflect the underlying symmetries of the given physical problem. In this work, we focus on complex scalar field theory on a two-dimensional lattice and investigate the benefits of using group equivariant convolutional neural network architectures based on the translation group. For a meaningful comparison, we conduct a systematic search for equivariant and non-equivariant neural network architectures and apply them to various regression and classification tasks. We demonstrate that in most of these tasks our best equivariant architectures can perform and generalize significantly better than their non-equivariant counterparts, which applies not only to physical parameters beyond those represented in the training set, but also to different lattice sizes.

Quantum simulation of cosmic inflation

In this paper, we generalize Jordan-Lee-Preskill, an algorithm for simulating flat-space quantum field theories, to 3+1 dimensional inflationary spacetime. The generalized algorithm contains the encoding treatment, the initial state preparation, the inflation process, and the quantum measurement of cosmological observables at late time. The algorithm is helpful for obtaining predictions of cosmic non-Gaussianities, serving as useful benchmark problems for quantum devices, and checking assumptions made about interacting vacuum in the inflationary perturbation theory. Components of our work also include a detailed discussion about the lattice regularization of the cosmic perturbation theory, a detailed discussion about the in-in formalism, a discussion about encoding using the HKLL-type formula that might apply for both dS and AdS spacetimes, a discussion about bounding curvature perturbations, a description of the three-party Trotter simulation algorithm for time-dependent Hamiltonians, a ground state projection algorithm for simulating gapless theories, a discussion about the quantum-extended Church-Turing Thesis, and a discussion about simulating cosmic reheating in quantum devices.

RI/MOM renormalization of the parton quasidistribution functions in lattice regularization

We analyze the lattice spacing dependence for the pion unpolarized matrix element of a quark bilinear operator with a Wilson link (parton quasidistribution functions operator, quasi-PDFs) in the rest frame, using 13 lattice spacings ranging from 0.032 to 0.121 fm. We compare results for three different fermion actions with or without good chiral symmetry on dynamical gauge ensembles from three collaborations. This investigation is motivated by the fact that the gauge link generates a $1/a$ divergence, the cancellation of which in many ratios can be numerically tricky. Indeed, our results show that this cancellation deteriorates with decreasing lattice spacing, and that the RI/MOM method leaves a linearly divergent residue for quasi-PDFs. We also show that in the Landau gauge the interaction between the Wilson link and the external state results in a linear divergence which depends on the discretized fermion action.

Scale symmetry breaking, quantum anomalous energy and proton mass decomposition

We study the anomalous scale symmetry breaking effects on the proton mass in QCD due to quantum fluctuations at ultraviolet scales. We confirm that a novel contribution naturally arises as a part of the proton mass, which we call the quantum anomalous energy (QAE). We discuss the QAE origins in both lattice and dimensional regularizations and demonstrate its role as a scheme-and-scale independent component in the mass decomposition. We further argue that QAE role in the proton mass resembles a dynamical Higgs mechanism, in which the anomalous scale symmetry breaking field generates mass scales through its vacuum condensate, as well as its static and dynamical responses to the valence quarks. We demonstrate some of our points in two simpler but closely related quantum field theories, namely the 1+1 dimensional non-linear sigma model in which QAE is non-perturbative and scheme-independent, and QED where the anomalous energy effect is perturbative calculable.

Rotated twisted-mass: a convenient regularization scheme for isospin breaking QCD and QED lattice calculations

We propose a scheme of lattice twisted-mass fermion regularization which is particularly convenient for application to isospin breaking (IB) QCD and QED calculations, based in particular on the so called RM123 approach, in which the IB terms of the action are treated as a perturbation. The main, practical advantage of this scheme is that it allows the calculation of IB effects on some mesonic observables, like e.g. the pi ^+ - pi ^0 mass splitting, using lattice correlation functions in which the quark and antiquark fields in the meson are regularized with opposite values of the Wilson parameter r. These correlation functions are found to be affected by much smaller statistical fluctuations, with respect to the analogous functions in which quark and antiquark fields are regularized with the same value of r. Two numerical application of this scheme, that we call rotated twisted-mass, within pure QCD and QCD + QED respectively, are also provided for illustration.

Generalized Gelfand–Yaglom formula for a discretized quantum mechanic system

The Gelfand–Yaglom formula relates the regularized determinant of a differential operator to the solution of an initial value problem. Here, we develop a generalized Gelfand–Yaglom formula for a Hamiltonian system with Lagrangian boundary conditions in the discrete and continuous settings. Later, we analyze the convergence of the discretized Hamilton–Jacobi operator and propose a lattice regularization for the determinant.

Investigating a (3+1)D topological θ -term in the Hamiltonian formulation of lattice gauge theories for quantum and classical simulations

Quantum technologies offer the prospect to efficiently simulate sign-problem afflicted regimes in lattice field theory, such as the presence of topological terms, chemical potentials, and out-of-equilibrium dynamics. In this work, we derive the (3+1)D topological $\theta$-term for Abelian and non-Abelian lattice gauge theories in the Hamiltonian formulation, paving the way towards Hamiltonian-based simulations of such terms on quantum and classical computers. We further study numerically the zero-temperature phase structure of a (3+1)D U(1) lattice gauge theory with the $\theta$-term via exact diagonalization for a single periodic cube. In the strong coupling regime, our results suggest the occurrence of a phase transition at constant values of $\theta$, as indicated by an avoided level-crossing and abrupt changes in the plaquette expectation value, the electric energy density, and the topological charge density. These results could in principle be cross-checked by the recently developed (3+1)D tensor network methods and quantum simulations, once sufficient resources become available.

Application of the dynamic mean field theory to fluid transport in slit pores.

We explore the applicability of the lattice model and dynamic mean field theory as a computationally efficient tool to study transport across heterogeneous porous media, such as mixed matrix membranes. As a starting point and to establish some basic definitions of properties analogous to those in the off-lattice systems, we consider transport across simple models of porous materials represented by a slit pore in a chemical potential gradient. Using this simple model, we investigate the distribution of density and flux under steady state conditions, define the permeability across the system, and explore how this property depends on the length of the pore and the solid-fluid interactions. Among other effects, we observe that the flux in the system goes through a maximum as the solid-fluid interaction is varied from weak to strong. This effect is dominated by the behavior of the fluid near the walls and is also confirmed by off-lattice molecular dynamics simulations. We further extend this study to explore transport across heterogeneous slit pore channels composed of two solids with different values of solid-fluid interaction strengths. We demonstrate that the lattice models and dynamic mean field theory provide a useful framework to pose questions on the accuracy and applicability of the classical theories of transport across heterogeneous porous systems.

Scaling Limits of Lattice Quantum Fields by Wavelets

We present a rigorous renormalization group scheme for lattice quantum field theories in terms of operator algebras. The renormalization group is considered as an inductive system of scaling maps between lattice field algebras. We construct scaling maps for scalar lattice fields using Daubechies’ wavelets, and show that the inductive limit of free lattice ground states exists and the limit state extends to the familiar massive continuum free field, with the continuum action of spacetime translations. In particular, lattice fields are identified with the continuum field smeared with Daubechies’ scaling functions. We compare our scaling maps with other renormalization schemes and their features, such as the momentum shell method or block-spin transformations.

Testing the dark SU(N) Yang-Mills theory confined landscape: From the lattice to gravitational waves

We pave the way for future gravitational-wave detection experiments, such as the big bang observer and DECIGO, to constraint dark sectors made of $SU(N)$ Yang-Mills confined theories. We go beyond the state-of-the-art by combining first principle lattice results and effective field theory approaches to infer essential information about the nonperturbative dark deconfinement phase transition driving the generation of gravitational-waves in the early Universe, such as the order, duration and energy budget of the phase transition which are essential in establishing the strength of the resulting gravitational-wave signal.

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

Lattice regularisation and entanglement structure of the Gross-Neveu model

We construct a Hamiltonian lattice regularisation of the N-flavour Gross-Neveu model that manifestly respects the full O(2N) symmetry, preventing the appearance of any unwanted marginal perturbations to the quantum field theory. In the context of this lattice model, the dynamical mass generation is intimately related to the Coleman-Mermin-Wagner and Lieb-Schultz-Mattis theorems. In particular, the model can be interpreted as lying at the first order phase transition line between a trivial and symmetry-protected topological (SPT) phase, which explains the degeneracy of the elementary kink excitations. We show that our Hamiltonian model can be solved analytically in the large N limit, producing the correct expression for the mass gap. Furthermore, we perform extensive numerical matrix product state simulations for N = 2, thereby recovering the emergent Lorentz symmetry and the proper non-perturbative mass gap scaling in the continuum limit. Finally, our simulations also reveal how the continuum limit manifests itself in the entanglement spectrum. As expected from conformal field theory we find two conformal towers, one tower spanned by the linear representations of O(4), corresponding to the trivial phase, and the other by the projective (i.e. spinor) representations, corresponding to the SPT phase.

Instanton solution for Schwinger production of ’t Hooft-Polyakov monopoles

We present the results of an explicit numerical computation of a novel instanton in Georgi-Glashow SU(2) theory. The instanton is physically relevant as a mediator of Schwinger production of 't Hooft-Polyakov magnetic monopoles from strong magnetic fields. In weak fields, the pair production rate has previously been computed using the worldline approximation, which breaks down in strong fields due to the effects of finite monopole size. Using lattice field theory we have overcome this limit, including finite monopole size effects to all orders. We demonstrate that a full consideration of the internal monopole structure results in an enhancement to the pair production rate, and confirm earlier results that monopole production becomes classical at the Ambjorn-Olesen critical field strength.

Equilibrium chiral magnetic effect: Spatial inhomogeneity, finite temperature, interactions

We discuss equilibrium relativistic fermionic systems in lattice regularization, and extend the consideration of chiral magnetic effect to systems with spatial inhomogeneity and finite temperature. Besides, we take into account interactions due to exchange by gauge bosons. We find that the equilibrium chiral magnetic conductivity remains equal to zero.

Normalizing flows and the real-time sign problem

Normalizing flows have recently been applied to the problem of accelerating Markov chains in lattice field theory. We propose a generalization of normalizing flows that allows them to applied to theories with a sign problem. These complex normalizing flows are closely related to contour deformations (i.e. the generalized Lefschetz thimble method), which been applied to sign problems in the past. We discuss the question of the existence of normalizing flows: they do not exist in the most general case, but we argue that exact normalizing flows are likely to exist for many physically interesting problems, including cases where the Lefschetz thimble decomposition has an intractable sign problem. Finally, normalizing flows can be constructed in perturbation theory. We give numerical results on their effectiveness across a range of couplings for the Schwinger-Keldysh sign problem associated to a real scalar field in $0+1$ dimensions.

Time-dependent Entanglement Entropy in Dissipative Conformal Theories: TFD Approach

In this work the TFD formalism is explored in order to study a dissipative time-dependent thermal vacuum. This state is a consequence of a particular interaction between two theories, which can be interpreted as two conformal theories defined at the two asymptotic boundaries of an AdS black hole. The initial state is prepared to be the equilibrium TFD thermal vacuum. The interaction causes dissipation from the point of view of observers who measure observables in one of the boundaries. We show that the vacuum evolves as an entangled state at finite temperature and the dissipative dynamics is controlled by the time-dependent entropy operator, defined in the non-equilibrium TFD framework. We use lattice field theory techniques to calculate the non-equilibrium thermodynamic entropy and the finite temperature entanglement entropy. We show that both grow linearly with time.