High Energy Physics - Lattice (hep-lat)

Topologically-ordered phases are stable to local perturbations, and topological quantum error-correcting codes enjoy thresholds to local errors. We connect the two notions of stability by constructing classical statistical mechanics models for decoding general CSS codes and classical linear codes. Our construction encodes correction success probabilities under uncorrelated bit-flip and phase-flip errors, and simultaneously describes a generalized Z2 lattice gauge theory with quenched disorder. We observe that the clean limit of the latter is precisely the discretized imaginary time path integral of the corresponding quantum code Hamiltonian when the errors are turned into a perturbative X or Z magnetic field. Motivated by error correction considerations, we define general order parameters for all such generalized Z2 lattice gauge theories, and show that they are generally lower bounded by success probabilities of error correction. For CSS codes satisfying the LDPC condition and with a sufficiently large code distance, we prove the existence of a low temperature ordered phase of the corresponding lattice gauge theories, particularly for those lacking Euclidean spatial locality and/or when there is a nonzero code rate. We further argue that these results provide evidence to stable phases in the corresponding perturbed quantum Hamiltonians, obtained in the limit of continuous imaginary time. To do so, we distinguish space- and time-like defects in the lattice gauge theory. A high free-energy cost of space-like defects corresponds to a successful "memory experiment" and suppresses the energy splitting among the ground states, while a high free-energy cost of time-like defects corresponds to a successful "stability experiment" and points to a nonzero gap to local excitations.

We develop a hybrid oscillator-qubit processor framework for quantum simulation of strongly correlated fermions and bosons that avoids the boson-to-qubit mapping overhead encountered in qubit hardware. This framework gives exact decompositions of particle interactions such as density-density terms and gauge-invariant hopping, as well as approximate methods based on the Baker-Campbell Hausdorff formulas including the magnetic field term for the $U(1)$ quantum link model in $(2+1)$D. We use this framework to show how to simulate dynamics using Trotterisation, perform ancilla-free partial error detection using Gauss's law, measure non-local observables, estimate ground state energies using a oscillator-qubit variational quantum eigensolver as well as quantum signal processing, and we numerically study the influence of hardware errors in circuit QED experiments. To show the advantages over all-qubit hardware, we perform an end-to-end comparison of the gate complexity for the gauge-invariant hopping term and find an improvement of the asymptotic scaling with the boson number cutoff $S$ from $\mathcal{O}(\log(S)^2)$ to $\mathcal{O}(1)$ in our framework as well as, for bosonic matter, a constant factor improvement of better than $10^4$. We also find an improvement from $\mathcal{O}(\log(S))$ to $\mathcal{O}(1)$ for the $U(1)$ magnetic field term. While our work focusses on an implementation in superconducting hardware, our framework can also be used in trapped ion, and neutral atom hardware. This work establishes digital quantum simulation with hybrid oscillator-qubit hardware as a viable and advantageous method for the study of qubit-boson models in materials science, chemistry, and high-energy physics.

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a different approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as $\widetilde{O}(\lambda N^7 |\Omega|^3/(M^{5/2} \epsilon^{3/2})$ where $\lambda$ is the coupling strength, $N$ is the occupation cutoff $|\Omega|$ is the volume of the spatial lattice, $M$ is the mass of the particles and $\epsilon$ is the uncertainty in the energy calculation used for the $S$-matrix determination. Qubitization in the field basis scales as $\widetilde{O}(|\Omega|^2 (k^2 \Lambda +kM^2)/\epsilon)$ where $k$ is the cutoff in the field and $\Lambda$ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of $4\times 10^6$ physical qubits and $10^{12}$ $T$-gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns, placing simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Block Encoding (BE) is a crucial subroutine in many modern quantum algorithms, including those with near-optimal scaling for simulating quantum many-body systems, which often rely on Quantum Signal Processing (QSP). Currently, the primary methods for constructing BEs are the Linear Combination of Unitaries (LCU) and the sparse oracle approach. In this work, we demonstrate that QSP-based techniques, such as Quantum Singular Value Transformation (QSVT) and Quantum Eigenvalue Transformation for Unitary Matrices (QETU), can themselves be efficiently utilized for BE implementation. Specifically, we present several examples of using QSVT and QETU algorithms, along with their combinations, to block encode Hamiltonians for lattice bosons, an essential ingredient in simulations of high-energy physics. We also introduce a straightforward approach to BE based on the exact implementation of Linear Operators Via Exponentiation and LCU (LOVE-LCU). We find that, while using QSVT for BE results in the best asymptotic gate count scaling with the number of qubits per site, LOVE-LCU outperforms all other methods for operators acting on up to $\lesssim11$ qubits, highlighting the importance of concrete circuit constructions over mere comparisons of asymptotic scalings. Using LOVE-LCU to implement the BE, we simulate the time evolution of single-site and two-site systems in the lattice $\varphi^4$ theory using the Generalized QSP algorithm and compare the gate counts to those required for Trotter simulation.

Simulation of quantum systems of a large number of strongly interacting particles persists as one of the most challenging, and computationally demanding, tasks in classical simulation, involving both non-relativistic applications like condensed matter physics and quantum chemistry, as well as relativistic applications like lattice gauge theory simulation. One of the major motivations for building a fault-tolerant quantum computer is the efficient simulation of many-body systems on such a device. While significant developments have been made in the quantum simulation of non-relativistic systems, the simulation of lattice gauge theories has lagged behind, with state-of-the-art Trotterized simulations requiring many orders of magnitude more resources than non-relativistic simulation, in stark contrast to the similar difficulty of these tasks in classical simulation. In this work, we conduct an in-depth analysis of the cost of simulating Abelian and non-Abelian lattice gauge theories in the Kogut-Susskind formulation using simulation methods with near-optimal scaling in system size, evolution time, and error. We provide explicit circuit constructions, as well as T-gate counts and qubit counts for the entire simulation algorithm. This investigation, the first of its kind, leads to up to 25 orders of magnitude improvement over Trotterization in spacetime volume for non-Abelian simulations. Such a dramatic improvement results largely from our algorithm having polynomial scaling with the number of colors, as opposed to exponential scaling in existing approaches. Our work demonstrates that the use of advanced algorithmic techniques leads to dramatic reductions in the cost of ab initio simulations of fundamental interactions, bringing it in step with resources required for first principles quantum simulation of chemistry and condensed matter physics.

Non-Abelian anyons, a promising platform for fault-tolerant topological quantum computation, adhere to the charge super-selection rule (cSSR), which imposes restrictions on physically allowed states and operations. However, the ramifications of cSSR and fusion rules in anyonic quantum information theory remain largely unexplored. In this study, we unveil that the information-theoretic characteristics of anyons diverge fundamentally from those of non-anyonic systems such as qudits, bosons, and fermions and display intricate structures. In bipartite anyonic systems, pure states may have different marginal spectra, and mixed states may contain pure marginal states. More striking is that in a pure entangled state, parties may lack equal access to entanglement. This entanglement asymmetry is manifested in quantum teleportation employing an entangled anyonic state shared between Alice and Bob, where Alice can perfectly teleport unknown quantum information to Bob, but Bob lacks this capability. These traits challenge conventional understanding, necessitating new approaches to characterize quantum information and correlations in anyons. We expect that these distinctive features will also be present in non-Abelian lattice gauge field theories. Our findings significantly advance the understanding of the information-theoretic aspects of anyons and may lead to realizations of quantum communication and cryptographic protocols where one party holds sway over the other.

One of the most challenging problems in the computational study of localization in quantum manybody systems is to capture the effects of rare events, which requires sampling over exponentially many disorder realizations. We implement an efficient procedure on a quantum processor, leveraging quantum parallelism, to efficiently sample over all disorder realizations. We observe localization without disorder in quantum many-body dynamics in one and two dimensions: perturbations do not diffuse even though both the generator of evolution and the initial states are fully translationally invariant. The disorder strength as well as its density can be readily tuned using the initial state. Furthermore, we demonstrate the versatility of our platform by measuring Renyi entropies. Our method could also be extended to higher moments of the physical observables and disorder learning.

Quantum Fourier transformations are an essential component of many quantum algorithms, from prime factoring to quantum simulation. While the standard abelian QFT is well-studied, important variants corresponding to \emphnonabelian groups of interest have seen less development. In particular, fast nonabelian Fourier transformations are important components for both quantum simulations of field theories as well as approaches to the nonabelian hidden subgroup problem. In this work, we present fast quantum Fourier transformations for a number of nonabelian groups of interest for high energy physics, $\mathbb{BT}$, $\mathbb{BO}$, $\Delta(27)$, $\Delta(54)$, and $\Sigma(36\times3)$. For each group, we derive explicit quantum circuits and estimate resource scaling for fault-tolerant implementations. Our work shows that the development of a fast Fourier transformation can substantively reduce simulation costs by up to three orders of magnitude for the finite groups that we have investigated.

We construct the primitive gate set for the digital quantum simulation of the 108-element $\Sigma(36\times3)$ group. This is the first time a nonabelian crystal-like subgroup of $SU(3)$ has been constructed for quantum simulation. The gauge link registers and necessary primitives -- the inversion gate, the group multiplication gate, the trace gate, and the $\Sigma(36\times3)$ Fourier transform -- are presented for both an eight-qubit encoding and a heterogeneous three-qutrit plus two-qubit register. For the latter, a specialized compiler was developed for decomposing arbitrary unitaries onto this architecture.

We develop a Hamiltonian formalism for simulating interacting chiral fermions on the lattice while preserving unitarity and locality and without breaking the chiral symmetry. The fermion doubling problem is circumvented by constructing a Fock space endowed with a semi-definite norm. When projecting our theory on the the single-particle sector, we recover the framework of Stacey fermions, and we demonstrate that the scaling limit of the free model recovers the chiral fermion field. Technically, we make use of a matrix product operator norm to mimick the boundary of a higher dimensional topological theory. As a proof of principle, we consider a single Weyl fermion on a periodic ring with Hubbard-type nearest-neighbor interactions and construct a variational generalized DMRG code to demonstrate that the ground state for large system sizes can be determined efficiently. As our tensor network approach does not exhibit any sign problem, we can add a chemical potential and study real-time evolution.

Lattice gauge theories (LGTs) can be employed to understand a wide range of phenomena, from elementary particle scattering in high-energy physics to effective descriptions of many-body interactions in materials. Studying dynamical properties of emergent phases can be challenging as it requires solving many-body problems that are generally beyond perturbative limits. We investigate the dynamics of local excitations in a $\mathbb{Z}_2$ LGT using a two-dimensional lattice of superconducting qubits. We first construct a simple variational circuit which prepares low-energy states that have a large overlap with the ground state; then we create particles with local gates and simulate their quantum dynamics via a discretized time evolution. As the effective magnetic field is increased, our measurements show signatures of transitioning from deconfined to confined dynamics. For confined excitations, the magnetic field induces a tension in the string connecting them. Our method allows us to experimentally image string dynamics in a (2+1)D LGT from which we uncover two distinct regimes inside the confining phase: for weak confinement the string fluctuates strongly in the transverse direction, while for strong confinement transverse fluctuations are effectively frozen. In addition, we demonstrate a resonance condition at which dynamical string breaking is facilitated. Our LGT implementation on a quantum processor presents a novel set of techniques for investigating emergent particle and string dynamics.

Tensor network (TN) states, including entanglement renormalization (ER), can encompass a wider variety of entangled states. When the entanglement structure of the quantum state of interest is non-uniform in real space, accurately representing the state with a limited number of degrees of freedom hinges on appropriately configuring the TN to align with the entanglement pattern. However, a proposal has yet to show a structural search of ER due to its high computational cost and the lack of flexibility in its algorithm. In this study, we conducted an optimal structural search of TN, including ER, based on the reconstruction of their local structures with respect to variational energy. Firstly, we demonstrated that our algorithm for the spin-$1/2$ tetramer singlets model could calculate exact ground energy using the multi-scale entanglement renormalization ansatz (MERA) structure as an initial TN structure. Subsequently, we applied our algorithm to the random XY models with the two initial structures: MERA and the suitable structure underlying the strong disordered renormalization group. We found that, in both cases, our algorithm achieves improvements in variational energy, fidelity, and entanglement entropy. The degree of improvement in these quantities is superior in the latter case compared to the former, suggesting that utilizing an existing TN design method as a preprocessing step is important for maximizing our algorithm's performance.

Lattice gauge theories (LGTs) describe a broad range of phenomena in condensed matter and particle physics. A prominent example is confinement, responsible for bounding quarks inside hadrons such as protons or neutrons. When quark-antiquark pairs are separated, the energy stored in the string of gluon fields connecting them grows linearly with their distance, until there is enough energy to create new pairs from the vacuum and break the string. While such phenomena are ubiquitous in LGTs, simulating the resulting dynamics is a challenging task. Here, we report the observation of string breaking in synthetic quantum matter using a programmable quantum simulator based on neutral atom arrays. We show that a (2+1)D LGT with dynamical matter can be efficiently implemented when the atoms are placed on a Kagome geometry, with a local U(1) symmetry emerging from the Rydberg blockade, while long-range Rydberg interactions naturally give rise to a linear confining potential for a pair of charges, allowing us to tune both their masses as well as the string tension. We experimentally map out the corresponding phase diagram by adiabatically preparing the ground state of the atom array in the presence of defects, and observe substructure of the confined phase, distinguishing regions dominated by fluctuating strings or by broken string configurations. Finally, by harnessing local control over the atomic detuning, we quench string states and observe string breaking dynamics exhibiting a many-body resonance phenomenon. Our work paves a way to explore phenomena in high-energy physics using programmable quantum simulators.

We propose and analyze a deformation of the 3+1d lattice $\mathbb Z_2$ gauge theory that preserves the non-invertible Wegner duality symmetry at the self-dual point. We identify a frustration-free point along this deformation where there are nine exactly degenerate ground states (on a periodic cubic lattice) even at finite volume. One of these ground states is a trivial product state and the rest are the topologically-ordered ground states of the 3+1d toric code. We also prove that the frustration-free point is gapped in the thermodynamic limit. Our model, therefore, realizes a gapped phase with spontaneously broken Wegner duality symmetry. Furthermore, by imposing the Gauss law constraints energetically, all the above features can be realized on a tensor product Hilbert space. Finally, we discuss a generalization of this deformation to the 3+1d lattice $\mathbb Z_N$ gauge theory and conjecture the possible phase diagram.

Simulating fermions coupled to spin degrees of freedom, relevant for a range of quantum field theories, represents a promising application for quantum simulators. Mapping fermions to qubits is challenging in $2+1$ and higher spacetime dimensions, and mapping bosons demands substantial quantum-computational overhead. These features complicate the realization of mixed fermion-boson quantum systems in digital quantum computers. We propose a native fermion-(large-)spin analog quantum simulator by utilizing dopant arrays in silicon. Specifically, we show how to use a dynamical lattice of coupled nuclear spins and conduction-band electrons to realize a quantum field theory: an extended Jackiw-Rebbi model involving coupled fermions and quantum rotors. We demonstrate the feasibility of observing dynamical mass generation and a confinement-deconfinement quantum phase transition in 1+1 dimensions on this platform, even in the presence of strong long-range Coulomb interactions. Furthermore, we employ finite-temperature Hartree-Fock-Bogoliubov simulations to investigate the dynamics of mass generation in two-dimensional square and honeycomb arrays, showing that this phenomenon can be simulated with realistic experimental parameters. Our findings reveal two distinct phases, and demonstrate robustness against the addition of Coulomb interactions. Finally, we discuss experimental signatures of the phases through transport and local charge sensing in dopant arrays. This study lays the foundation for quantum simulations of quantum field theories exhibiting fermions coupled to spin degrees of freedom using donors in silicon.

We present a quantum computational framework for SU(2) lattice gauge theory, leveraging continuous variables instead of discrete qubits to represent the infinite-dimensional Hilbert space of the gauge fields. We consider a ladder as well as a two-dimensional grid of plaquettes, detailing the use of gauge fixing to reduce the degrees of freedom and simplify the Hamiltonian. We demonstrate how the system dynamics, ground states, and energy gaps can be computed using the continuous-variable approach to quantum computing. Our results indicate that it is feasible to study non-Abelian gauge theories with continuous variables, providing new avenues for understanding the real-time dynamics of quantum field theories.

We define and solve the $\text{U(1)}$ Chern-Simons-Maxwell theory on spacetime lattice, with an emphasis on the chirality of the theory. Realizing Chern-Simons theory on lattice has been a problem of interest for decades, and over the years it has gradually become clear that there are two key points: 1) Some non-topological term, such as a Maxwell term, is necessary -- this is true even in the continuum, but more manifestly on the lattice; 2) the $\text{U(1)}$ gauge field should be implemented in the Villainized form to retain its topological properties. Putting the two ideas together seriously, we show all interesting properties of a chiral Chern-Simons theory are reproduced in an explicitly regularized manner on the lattice. These include the bosonic and fermionic level quantization, the bulk and chiral edge spectrum, the Wilson loop flux attachment (with point-split framing or geometric framing depending on the Maxwell coupling), the Wilson loop spin, the ground state degeneracy, and, most non-trivially, the chiral gravitational anomaly.

Tensor network renormalization group maps study critical points of 2d lattice models like the Ising model by finding the fixed point of the RG map. In a prior work arXiv:2408.10312 we showed that by adding a rotation to the RG map, the Newton method could be implemented to find an extremely accurate fixed point. For a particular RG map (Gilt-TNR) we studied the spectrum of the Jacobian of the RG map at the fixed point and found good agreement between the eigenvalues corresponding to relevant and marginal operators and their known exact values. In this companion work we use two further methods to extract many more scaling dimensions from this Newton method fixed point, and compare the numerical results with the predictions of conformal field theory (CFT). The first method is the well-known transfer matrix. We introduce some extensions of this method that provide spins of the CFT operators modulo an integer. We find good agreement for the scaling dimensions and spins up to $\Delta=3\frac1 8$. The second method we refer to as the lattice dilatation operator (LDO). This lesser-known method obtains good agreement with CFT up to $\Delta=6$. Moreover, the inclusion of a rotation in the RG map makes it possible to extract spins of the CFT operators modulo 4 from this method. Some of the eigenvalues of the Jacobian of the RG map can come from perturbations associated with total derivative interactions and so are not universal. In some past studies (arXiv:2102.08136, arXiv:2305.09899) such non-universal eigenvalues did not appear in the Jacobian. We explain this surprising result by showing that their RG map has the unusual property that the Jacobian is equivalent to the LDO operator.

In the tensor network approach to statistical physics, properties of the critical point of a 2D lattice model are encoded by a four-legged tensor which is a fixed point of an RG map. The traditional way to find the fixed point tensor consists in iterating the RG map after having tuned the temperature to criticality. Here we develop a different and more direct technique, which solves the fixed point equation via the Newton method. This is challenging due to the existence of marginal deformations -- linear transformations of the coordinate frame, which parametrize a two-dimensional family of fixed points. We address this challenge by including a 90 degree rotation into the RG map. This flips the sign of the problematic marginal eigenvalues, rendering the fixed point isolated and accessible via the Newton method. We demonstrate the power of this technique via explicit computations for the 2D Ising model. Using the Gilt-TNR algorithm at bond dimension $\chi=30$, we find the fixed point tensor with $10^{-9}$ accuracy, much higher than what was previously achieved.

Quantum many-body scarring (QMBS) is an intriguing mechanism of ergodicity breaking that has recently spurred significant attention. Particularly prominent in Abelian lattice gauge theories (LGTs), an open question is whether QMBS nontrivially arises in non-Abelian LGTs. Here, we present evidence of robust QMBS in a non-Abelian SU(2) LGT with dynamical matter. Starting in product states that require little experimental overhead, we show that prominent QMBS arises for certain quenches, facilitated through meson and baryon-antibaryon excitations, highlighting its non-Abelian nature. The uncovered scarred dynamics manifests as long-lived coherent oscillations in experimentally accessible local observables as well as prominent revivals in the state fidelity. Our findings bring QMBS to the realm of non-Abelian LGTs, highlighting the intimate connection between scarring and gauge symmetry, and are amenable for observation in a recently proposed trapped-ion qudit quantum computer.

Spontaneous particle-pair formation is a fundamental phenomenon in nature. It can, for example, appear when the potential energy between two particles increases with separation, as if they were connected by a tense string. Beyond a critical separation, new particle pairs can form, causing the string to break. String-breaking dynamics in quantum chromodynamics play a vital role in high-energy particle collisions and early universe evolution. Simulating string evolution and hadron formation is, therefore, a grand challenge in modern physics. Quantum simulators, well-suited for studying dynamics, are expected to outperform classical computing methods. However, the required experimental capabilities to simulate string-breaking dynamics have not yet been demonstrated, even for simpler models of the strong force. We experimentally probe, for the first time, the spatiotemporal dynamics of string-breaking in a (1+1)-dimensional $\mathbb{Z}_2$ lattice gauge theory using a fully programmable trapped-ion quantum simulator. We emulate external static charges and strings via site-dependent magnetic-field control enabled by a dual array of tightly focused laser beams targeting individual ions. First, we study how confinement affects isolated charges, finding that they freely spread without string tension but exhibit localized oscillations when tension is increased. Then, we observe and characterize string-breaking dynamics of a string stretched between two static charges after an abrupt increase in string tension. Charge pairs appear near the string edges and spread into the bulk, revealing a route to dynamical string-breaking distinct from the conventional Schwinger mechanism. Our work demonstrates that analog quantum simulators have achieved the necessary control to explore string-breaking dynamics, which may ultimately be relevant to nuclear and high-energy physics.

Spatially modulated symmetries have emerged since the discovery of fractons, which characterize unconventional topological phases with mobility-constrained quasiparticle excitations. On the other hand, non-invertible duality defects have attracted substantial attention in communities of high energy and condensed matter physics due to their deep insight into quantum anomalies and exotic phases of matter. However, the connection between these exotic symmetries and defects has not been fully explored. In this paper, we construct concrete lattice models with non-invertible duality defects via gauging spatially modulated symmetries and investigate their exotic fusion rules. Specifically, we construct spin models with subsystem symmetries or dipole symmetries on one, two, and three-dimensional lattices. Gauging subsystem symmetries leads to non-invertible duality defects whose fusion rules involve $0$-form subsystem charges in two dimensions and higher-form operators that correspond to ``lineon'' excitations (excitations which are mobile along one-dimensional line) in three dimensions. Gauging dipole symmetries leads to non-invertible duality defects with dipole algebras that describe a hierarchical structure between global and dipole charges. Notably, the hierarchical structure of the dual dipole charges is inverted compared with the original ones. Our work provides a unified and systematic analytical framework for constructing exotic duality defects by gauging relevant symmetries.

Real-time dynamics of the Schwinger model provide an effective description of quark confinement out of equilibrium, routinely employed to model hadronization processes in particle-physics event generators. Ab-initio simulations of such non-perturbative processes are far beyond the reach of existing computational tools, and remain an outstanding open quest for quantum simulators to date. In this work we develop a general strategy to run Schwinger model dynamics on synthetic quantum spin lattices, such as neutral-atom or superconducting-qubit arrays. Our construction encodes the constrained fermionic and bosonic degrees of freedom of the model into the geometric shape of a magnetic interface. We show that global magnetic field patterns can drive coherent quantum dynamics of the interface equivalent to the lattice Schwinger Hamiltonian. We rigorously establish that the optimal array required for simulating real-time wave packet collisions and string fragmentation processes with accuracy $\epsilon$ in the continuum field-theory limit, is a quasi-one-dimensional ribbon with polynomial length and polylogarithmic width in $\epsilon^{-1}$. We finally discuss a concrete advantageous implementation using a state-of-the-art dual-species Rydberg atom array. This work opens up a path for near-term quantum simulators to address questions of immediate relevance to particle physics.

We demonstrate how to construct a fully gauge-fixed lattice Hamiltonian for a pure SU(2) gauge theory. Our work extends upon previous work, where a formulation of an SU(2) lattice gauge theory was developed that is efficient to simulate at all values of the gauge coupling. That formulation utilized maximal-tree gauge, where all local gauge symmetries are fixed and a residual global gauge symmetry remains. By using the geometric picture of an SU(2) lattice gauge theory as a system of rotating rods, we demonstrate how to fix the remaining global gauge symmetry. In particular, the quantum numbers associated with total charge can be isolated by rotating between the lab and body frames using the three Euler angles. The Hilbert space in this new `sequestered' basis partitions cleanly into sectors with differing total angular momentum, which makes gauge-fixing to a particular total charge sector trivial, particularly for the charge-zero sector. In addition to this sequestered basis inheriting the property of being efficient at all values of the coupling, we show that, despite the global nature of the final gauge-fixing procedure, this Hamiltonian can be simulated using quantum resources scaling only polynomially with the lattice volume.

The simulation of quantum lattice gauge theories faces the major challenge of maintaining gauge invariance, as various errors in the simulation push the state of the system out of the physical subspace of the system's exponentially larger Hilbert space. This paper outlines a method, based off of previous work, that uses gauge transformations in two ways. Firstly, the method exploits the Zeno effect by conducting frequent projections to suppress gauge drift. These projections utilize local gauge transformations to destructively interfere unphysical amplitudes via coupling to an ancillary qubit while the physical amplitudes are left untouched, up to a less than unity normalization factor. Secondly, gauge transformations are conducted throughout the time evolution of the system to hamper the speed of gauge drift. This paper demonstrates this method on a pure 1D SU$(2)$ toy model.

Simulating non-equilibrium phenomena in strongly-interacting quantum many-body systems, including thermalization, is a promising application of near-term and future quantum computation. By performing experiments on a digital quantum computer consisting of fully-connected optically-controlled trapped ions, we study the role of entanglement in the thermalization dynamics of a $Z_2$ lattice gauge theory in 2+1 spacetime dimensions. Using randomized-measurement protocols, we efficiently learn a classical approximation of non-equilibrium states that yields the gap-ratio distribution and the spectral form factor of the entanglement Hamiltonian. These observables exhibit universal early-time signals for quantum chaos, a prerequisite for thermalization. Our work, therefore, establishes quantum computers as robust tools for studying universal features of thermalization in complex many-body systems, including in gauge theories.

A natural definition for instanton density operator in lattice QCD has been long desired. We show this problem is, and has to be, resolved by higher category theory. The problem is resolved by refining at a conceptual level the Yang-Mills theory on lattice, in order to recover the homotopy information in the continuum, which would have been lost if we put the theory on lattice in the traditional way. The refinement needed is a generalization -- through the lens of higher category theory -- of the familiar process of Villainization that captures winding in lattice XY model and Dirac quantization in lattice Maxwell theory. The apparent difference is that Villainization is in the end described by principal bundles, hence familiar, but more general topological operators can only be captured on the lattice by more flexible structures beyond the usual group theory and fibre bundles, hence the language of categories becomes natural and necessary. The key structure we need for our particular problem is called multiplicative bundle gerbe, based upon which we can construct suitable structures to naturally define the 2d Wess-Zumino-Witten term, 3d skyrmion density operator and 4d hedgehog defect for lattice $S^3$ (pion vacua) non-linear sigma model, and the 3d Chern-Simons term, 4d instanton density operator and 5d Yang monopole defect for lattice $SU(N)$ Yang-Mills theory. In a broader perspective, higher category theory enables us to rethink more systematically the relation between continuum quantum field theory and lattice quantum field theory. We sketch a proposal towards a general machinery that constructs the suitably refined lattice degrees of freedom for a given non-linear sigma model or gauge theory in the continuum, realizing the desired topological operators on the lattice.

In recent years, the quantum computing method has been used to address the sign problem in traditional Monte Carlo lattice gauge theory (LGT) simulations. We propose that the Coulomb gauge (CG) should be used in quantum simulations of LGT. This is because the redundant degrees of freedom can be eliminated in CG. Therefore, the Hamiltonian in CG does not need to be gauge invariance, allowing the gauge field to be discretized naively. We point out that discretized gauge fields and fermion fields should be placed on momentum and position lattices, respectively. Under this scheme, the CG condition and Gauss's law can be conveniently preserved by solving algebraic equations of polarization vectors. We also discuss the procedure for mapping gauge fields to qubits, and then demonstrate the polynomial scaling of qubits and the complexity of time evolution. Finally, we calculate the vacuum expectation value (VEV) of the U(1) plaquette operator and the Wilson loop on a classical device to test the performance of our discretization scheme.

A framework for simulating the real-time dynamics of particles in dense matter using quantum computers is developed. As a demonstration, we perform classical simulations of heavy-hadrons propagating through a dense medium in the Schwinger model. Measurements of the time-dependent energy and charge density are used to identify mechanisms responsible for energy loss and hadron production (hadronization). A study of entanglement dynamics highlights the importance of quantum coherence between the particles that make up the dense medium. Throughout this work, care is taken to isolate, and remove, phenomena that arise solely from a finite lattice spacing. It is found that signatures of entanglement are more sensitive to lattice artifacts than other observables. Toward quantum simulations, we present an efficient method and the corresponding quantum circuits for preparing ground states in the presence of heavy mesons. These circuits are used to estimate the resources required to simulate in-medium energy loss and hadronization in the Schwinger model using quantum computers.

Non-Abelian gauge theories provide an accurate description of fundamental interactions, as both perturbation theory and quantum Monte Carlo computations in lattice gauge theory, it when applicable, show remarkable agreement with experimental data from particle colliders and cosmological observations. Complementing these computations, or combining them with quantum-inspired Hamiltonian lattice computations on quantum machines to improve continuum limit predictions with current quantum resources, is a formidable open challenge. Here, we propose a resource-efficient method to compute the running of the coupling in non-Abelian gauge theories beyond one spatial dimension. We first represent the Hamiltonian on periodic lattices in terms of loop variables and conjugate loop electric fields, exploiting the Gauss law to retain the gauge-independent ones. Then, we identify a local basis for small and large loops variationally to minimize the truncation error while computing the running of the coupling on small tori. Our method enables computations at arbitrary values of the bare coupling and lattice spacing with current quantum computers, simulators and tensor-network calculations, in regimes otherwise inaccessible.

Simulating nonequilibirum dynamics of a quantum many-body system is one of the promising applications of quantum computing. We simulate the time evolution of one-dimensional ${\bf Z}_2$ lattice gauge theory on IBM's superconducting 156-qubit device ibm\_fez. We consider the Floquet circuit made of the Trotter decomposition of Hamiltonian evolution and focus on its dynamics toward thermalization. Quantum simulation with the help of error mitigation is successful in running the Floquet circuit made of $38$ and $116$ qubits up to $10$ Trotter steps in the best case. This is enough to reach the early stage of prethermalization. Our work would be a benchmark for the potential power of quantum computing for high-energy physics problems.

Quantum computers can efficiently simulate highly entangled quantum systems, offering a solution to challenges facing classical simulation of Quantum Field Theories (QFTs). This paper presents an alternative to traditional methods for simulating the real-time evolution in QFTs by leveraging Hamiltonian Truncation (HT). As a use case, we study the Schwinger model, systematically reducing the complexity of the Hamiltonian via HT while preserving essential physical properties. For the observables studied in this paper, the HT approach converges quickly with the number of qubits, allowing for the interesting physics processes to be captured without needing many qubits. Identifying the truncated free Hamiltonian's eigenbasis with the quantum device's computational basis avoids the need for complicated and costly state preparation routines, reducing the algorithm's overall circuit depth and required coherence time. As a result, the HT approach to simulating QFTs on a quantum device is well suited to Noisy-Intermediate Scale Quantum (NISQ) devices, which have a limited number of qubits and short coherence times. We validate our approach by running simulations on a NISQ device, showcasing strong agreement with theoretical predictions. We highlight the potential of HT for simulating QFTs on quantum hardware.

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. We present a protocol that constrains quantum states by experimentally measured correlation functions which only scales polynomially with system size. This method enables measurement of a quantum state's entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol's usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest non-trivial expansion quantitatively predicts early time thermalization dynamics, including signaling the on-set of quantum chaos indicated by the entanglement Hamiltonian.

Tensor network methods are a class of numerical tools and algorithms to study many-body quantum systems in and out of equilibrium, based on tailored variational wave functions. They have found significant applications in simulating lattice gauge theories approaching relevant problems in high-energy physics. Compared to Monte Carlo methods, they do not suffer from the sign problem, allowing them to explore challenging regimes such as finite chemical potentials and real-time dynamics. Further development is required to tackle fundamental challenges, such as accessing continuum limits or computations of large-scale quantum chromodynamics. In this work, we review the state-of-the-art of Tensor Network methods and discuss a possible roadmap for algorithmic development and strategies to enhance their capabilities and extend their applicability to open high-energy problems. We provide tailored estimates of the theoretical and computational resource scaling for attacking large-scale lattice gauge theories.

We perform a tensor renormalization group simulation of non-Abelian gauge theory in three dimensions using a formulation based on the `armillary sphere.' In this formulation, matrix indices are completely traced out, eliminating the degeneracy in the singular value spectrum of the initial tensor. We demonstrate the usefulness of this technique by computing the average plaquette at zero temperature and the Polyakov loop susceptibility at finite temperatures for 2+1D SU(2) and SU(3) gauge theories. The deconfinement transition is identified for both gauge groups, with the SU(2) case being consistent with previous Monte Carlo results.

There has been recently considerable progress in understanding the nature of perturbation theory in UV free and gapped $2d$ integrable field theories with renormalon singularities. Thanks to Bethe ansatz and large $N$ techniques, non-perturbative corrections can also be computed and lead to the reconstruction of the trans-series for the free energy in presence of a chemical potential. This is an ideal arena to test resurgence in QFT and determine if and how the exact result can be reconstructed from the knowledge of the perturbative series only. In these notes we give a pedagogical introduction to this subject starting from the basics. In the first lecture we give an overview of applications in QFT of Borel resummations before the advent of resurgence. The second lecture introduces the key concepts of resurgence and finally in the third lecture we discuss a specific application in the context of the principal chiral field model. Extended version of three lectures given at IHES and review talks given at Les Diablerets and Mainz, in 2023.

We use the higher-order tensor renormalization group method to study the two-dimensional generalized XY model that admits integer and half-integer vortices. This model is the deformation of the classical XY model and has a rich phase structure consisting of nematic, ferromagnetic, and disordered phases and three transition lines belonging to the Berezinskii-Kosterlitz-Thouless (BKT) and Ising class. We explore the model for a wide range of temperatures, $T$, and the deformation parameter, $\Delta$, and compute specific heat along with integer and half-integer magnetic susceptibility, finding both BKT-like and Ising-like transitions and the region where they meet.

The circuit complexity for Hamiltonian simulation of the sparsified SYK model with $N$ Majorana fermions and $q=4$ (quartic interactions) which retains holographic features (referred to as `minimal holographic sparsified SYK') with $k\ll N^{3}/24$ (where $k$ is the total number of interaction terms times 1/$N$) using second-order Trotter method and Jordan-Wigner encoding is found to be $\widetilde{\mathcal{O}}(k^{p}N^{3/2} \log N (\mathcal{J}t)^{3/2}\varepsilon^{-1/2})$ where $t$ is the simulation time, $\varepsilon$ is the desired error in the implementation of the unitary $U = \exp(-iHt)$, $\mathcal{J}$ is the disorder strength, and $p < 1$. This complexity implies that with less than a hundred logical qubits and about $10^{6}$ gates, it will be possible to achieve an advantage in this model and simulate real-time dynamics up to scrambling time.

We construct a Grassmann tensor network representing the partition function of (1+1)-dimensional two-color QCD with staggered fermions. The Grassmann path integral is rewritten as the trace of a Grassmann tensor network by introducing two-component auxiliary Grassmann fields on every edge of the lattice. We introduce an efficient initial tensor compression scheme to reduce the size of initial tensors. The Grassmann bond-weighted tensor renormalization group approach is adopted to evaluate the quark number density, fermion condensate, and diquark condensate at different gauge couplings as a function of the chemical potential. Different transition behavior is observed as the quark mass is varied. We discuss the efficiency of our initial tensor compression scheme and the future application toward the corresponding higher-dimensional models.

We discuss the classical continuum limit of the string field theory dual to the $\mathrm{SU}(N)$ lattice gauge theory and investigate various fundamental phenomena in the continuum theory at the mean-field level. Our construction of the continuum theory is based on the concept of \it area derivative, which can be regarded as a generalization of the ordinary derivative $\partial/\partial x^\mu$ to operators acting on functional fields $\phi[C]$ on the loop space. The resultant continuum theory has a $\mathbb{Z}_N^{}$ $1$-form global symmetry, which originates in the $\mathbb{Z}_N^{}$ center symmetry in the gauge theory. We find that the confined and deconfined phases of the gauge theory are identified by the unbroken and broken phases of the $\mathbb{Z}_N^{}$ symmetry respectively by showing the Area/Perimeter law of the classical solution. In the broken phase, the low-energy effective theory is described by a $\mathrm{BF}$-type topological field theory and has a emergent $\mathbb{Z}_N^{}$ $(D-2)$-form global symmetry. The existence of the emergent symmetry is deeply related to $(D-2)$-dimensional topological configurations (i.e. center vortex for $D=4$), and we explicitly construct such a topological defect in the continuum theory. Finally, we also comment on the upper and lower critical dimensions of the gauge theory/string field theory.

Non-equilibrium Monte Carlo simulations based on Jarzynski's equality are a well-understood method to compute differences in free energy and also to sample from a target probability distribution without the need to thermalize the system under study. In each evolution, the system starts from a given base distribution at equilibrium and it is gradually driven out-of-equilibrium while evolving towards the target parameters. If the target distribution suffers from long autocorrelation times, this approach represents a promising candidate to mitigate critical slowing down. Out-of-equilibrium evolutions are conceptually similar to Normalizing Flows and they can be combined into a recently-developed architecture called Stochastic Normalizing Flows (SNFs). In this contribution we first focus on the promising scaling with the volume guaranteed by the purely stochastic approach in the $\mathrm{SU}(3)$ lattice gauge theory in 4 dimensions; then, we define an SNF by introducing gauge-equivariant layers between the out-of-equilibrium Monte Carlo updates, and we analyse the improvement obtained as well as the inherited scaling with the volume. Finally, we discuss how this approach can be systematically improved and how simulations of lattice gauge theories in four dimensions for large volumes and close to criticality can be realistically achieved.

$\mathbb{Z}_3$ lattice gauge theory is the simplest discrete gauge theory with three-quark bound states, i.e., baryons. Since it has a finite-dimensional Hilbert space, it can be used for testing quantum simulation of lattice gauge theory at nonzero baryon density. We discuss global and local gauge symmetries and their importance in quantum simulation. We perform quantum emulator calculation and demonstrate how to study the ground state property of baryonic matter.

Parton distribution functions (PDFs) describe universal properties of bound states and allow us to calculate scattering amplitudes in processes with large momentum transfer. Calculating PDFs involves the evaluation of matrix elements with a Wilson line in a light-cone direction. In contrast to Monte Carlo methods in Euclidean spacetime, these matrix elements can be directly calculated in Minkowski-space using the Hamiltonian formalism. The necessary spatial- and time-evolution can be efficiently applied using established tensor network methods. We present PDFs in the Schwinger model calculated with matrix product states.

The XY-plaquette model is the most straightforward lattice realization of a broad class of fractonic field theories that host quasiparticles with restricted mobility. The plaquette interaction appears naturally as a ring-exchange term in the low-energy description of exciton Bose liquids, cold atomic gases, and quantum dimer models. Using first-principle Monte Carlo simulations, we study the phase diagram and the vortex dynamics in the XY-plaquette model on a square lattice in two spatial dimensions. In its minimal formulation, the model contains a ring-exchange plaquette term in two spatial dimensions and a standard XY-link term in the (imaginary) time direction. We show that the phase diagram of the minimal XY-plaquette model possesses two phases: (i) a disordered vortex-dominated phase in which a single percolating vortex trajectory occupies the whole 3d spacetime; (ii) a partially disordered phase in which the vortices become partially immobile, with their worldlines strictly confined to one or several infinite two-dimensional planes. The spatial positions and spatial orientations (along $x$ or $y$ axis) of these vortex domain walls appear to be spontaneous. Individual vortices form a disordered system within each vortex domain wall, so the fractal spacetime dimension of vortex trajectories approaches $D_f = 2$. We argue that the appearance of the vortex walls could be interpreted as a consequence of the spontaneous breaking of a global internal symmetry in the compact XY-plaquette model.

Lattice gauge theories, the discretized cousins of continuum gauge theories, have become an important platform for the exploration of non-equilibrium phenomena beyond their original scope in the Standard Model. In particular, recent works have reported the possibility of disorder-free localization in the lattice Schwinger model. Using degenerate perturbation theory and numerical simulations based on exact diagonalization and matrix product states, we perform a detailed characterization of thermalization breakdown in the Schwinger model including its spectral properties, the structure of eigenstates, and out-of-equilibrium quench dynamics. We scrutinize the strong-coupling limit of the model, in which an intriguing, double-logarithmic-in-time, growth of entanglement was previously proposed from the initial vacuum state. We identify the origin of this ultraslow growth of entanglement as due to an approximate Hilbert space fragmentation and the emergence of a dynamical constraint on particle hopping, which gives rise to sharp jumps in the entanglement entropy dynamics within individual background charge sectors. Based on the statistics of jump times, we argue that the entanglement growth, averaged over charge sectors, is more naturally explained as either single-logarithmic or a weak power law in time. Our results thus suggest the existence of a single ergodicity-breaking regime due to Hilbert space fragmentation, whose properties are reminiscent of conventional many-body localization within the numerically accessible system sizes.

We study a single exactly massless staggered fermion in the fundamental representation of an $SU(2)$ gauge group. We utilize an nHYP-smeared fermion action supplemented with additional heavy Pauli-Villars fields which serve to decrease lattice artifacts. The phase diagram exhibits a clear two-phase structure with a conformal phase at weak coupling and a novel new phase, the Symmetric Mass Generation (SMG) phase, appearing at strong coupling. The SMG phase is confining with all states gapped and chiral symmetry unbroken. Our finite size scaling analysis provides strong evidence that the phase transition between these two phases is continuous, which would allow for the existence of a continuum SMG phase. Furthermore, the RG flows are consistent with a $\beta$-function that vanishes quadratically at the new fixed point suggesting that the $N_f=4$ flavor SU(2) gauge theory lies at the opening of the conformal window.

We present a thorough investigation of non-stabilizerness - a fundamental quantum resource that quantifies state complexity within the framework of quantum computing - in a one-dimensional U(1) lattice gauge theory. We show how non-stabilizerness is always extensive with volume, and has no direct relation to the presence of critical points. However, its derivatives typically display discontinuities across the latter: This indicates that non-stabilizerness is strongly sensitive to criticality, but in a manner that is very different from entanglement (that, typically, is maximal at the critical point). Our results indicate that error-corrected simulations of lattice gauge theories close to the continuum limit have similar computational costs to those at finite correlation length and provide rigorous lower bounds for quantum resources of such quantum computations.

Quantum computers can efficiently solve problems which are widely believed to lie beyond the reach of classical computers. In the near-term, hybrid quantum-classical algorithms, which efficiently embed quantum hardware in classical frameworks, are crucial in bridging the vast divide in the performance of the purely-quantum algorithms and their classical counterparts. Here, a hybrid quantum-classical algorithm is presented for the computation of non-perturbative characteristics of quantum field theories. The presented algorithm relies on a universal parametrized quantum circuit ansatz based on Euler and Cartan's decompositions of single and two-qubit operators. It is benchmarked by computing the energy spectra of lattice realizations of quantum field theories with both short and long range interactions. Low depth circuits are provided for false vacua as well as highly excited states corresponding to mesonic and baryonic excitations occurring in the analyzed models. The described algorithm opens a hitherto-unexplored avenue for the investigation of mass-ratios, scattering amplitudes and false-vacuum decays in quantum field theories.

In this paper, a Hamiltonian lattice formulation for 2+1D compact Maxwell-Chern-Simons theory is derived. We analytically solve this theory and demonstrate that the mass gap in the continuum limit matches the well-known continuum formula. Our formulation preserves topological features such as the quantization of the Chern-Simons level, the degeneracy of energy eigenstates, the non-trivial properties of Wilson loops, and the mutual and self statistics of anyons. This work lays the groundwork for future Hamiltonian-based simulations of Maxwell-Chern-Simons theory on classical and quantum computers.

This work explores the application of the concurrent variational quantum eigensolver (cVQE) for computing excited states of the Schwinger model. By designing suitable ansatz circuits utilizing universal SO(4) or SO(8) qubit gates, we demonstrate how to efficiently obtain the lowest two, four, and eight eigenstates with one, two, and three ancillary qubits for both vanishing and non-vanishing background electric field cases. Simulating the resulting quantum circuits classically with tensor network techniques, we demonstrate the capability of our approach to compute the two lowest eigenstates of systems with up to $\mathcal{O}(100)$ qubits. Given that our method allows for measuring the low-lying spectrum precisely, we also present a novel technique for estimating the additive mass renormalization of the lattice based on the energy gap. As a proof-of-principle calculation, we prepare the ground and first-excited states with one ancillary and four physical qubits on quantum hardware, demonstrating the practicality of using the cVQE to simulate excited states.