Understanding the equilibrium properties and out of equilibrium dynamics of quantum field theories are key aspects of fundamental problems in theoretical particle physics and cosmology. However, their classical simulation is highly challenging. In this work, we introduce a tensor network method for the classical simulation of continuous quantum field theories that is suitable for the study of low-energy eigenstates and out-of-equilibrium time evolution. The method is built on Hamiltonian truncation and tensor network techniques, bridging the gap between two successful approaches. One of the key developments is the exact construction of matrix product state representations of global projectors, crucial for the implementation of interacting theories. Despite featuring a relatively high computational effort, our method dramatically improves predictive precision compared to exact diagonalisation-based Hamiltonian truncation, allowing the study of so far unexplored parameter regimes and dynamical effects. We corroborate trust in the accuracy of the method by comparing it with exact theoretical results for ground state properties of the sine-Gordon model. We then proceed with discussing $(1+1)$-dimensional quantum electrodynamics, the massive Schwinger model, for which we accurately locate its critical point and study the growth and saturation of momentum-space entanglement in sudden quenches.
Understanding the equilibrium properties and out of equilibrium dynamics of quantum field theories are key aspects of fundamental problems in theoretical particle physics and cosmology. However, their classical simulation is highly challenging. In this work, we introduce a tensor network method for the classical simulation of continuous quantum field theories that is suitable for the study of low-energy eigenstates and out-of-equilibrium time evolution. The method is built on Hamiltonian truncation and tensor network techniques, bridging the gap between two successful approaches. One of the key developments is the exact construction of matrix product state representations of global projectors, crucial for the implementation of interacting theories. Despite featuring a relatively high computational effort, our method dramatically improves predictive precision compared to exact diagonalisation-based Hamiltonian truncation, allowing the study of so far unexplored parameter regimes and dynamical effects. We corroborate trust in the accuracy of the method by comparing it with exact theoretical results for ground state properties of the sine-Gordon model. We then proceed with discussing (1+1)-dimensional quantum electrodynamics, the massive Schwinger model, for which we accurately locate its critical point and study the growth and saturation of momentum-space entanglement in sudden quenches.
Entanglement is the fundamental difference between classical and quantum systems and has become one of the guiding principles in the exploration of high-and low-energy physics.The calculation of entanglement entropies in interacting quantum field theories, however, remains challenging.Here, we present the first method for the explicit computation of reduced density matrices of interacting quantum field theories using truncated Hamiltonian methods.The method is based on constructing an isomorphism between the Hilbert space of the full system and the tensor product of Hilbert spaces of subintervals.This naturally enables the computation of the von Neumann and arbitrary Rényi entanglement entropies as well as mutual information, logarithmic negativity, and other measures of entanglement.Our method is applicable to equilibrium states and nonequilibrium evolution in real time.It is model independent and can be applied to any Hamiltonian truncation method that uses a free basis expansion.We benchmark the method on the free Klein-Gordon theory finding excellent agreement with the analytic results.We further demonstrate its potential on the interacting sine-Gordon model, studying the scaling of von Neumann entropy in ground states and real-time dynamics following quenches of the model.
[eng] Quantum many-body dynamics is a vast field of physics at the intersection of statistical mechanics, condensed matter, quantum information, and computational physics. Understanding the time evolution of observable correlations—such as those involved in material responses—is both a milestone for technological advancement and a corner-stone for generating fundamental knowledge. At the intersection of emergent collective behavior in many-body systems, quantum physics, and information theory lies the classical simulation of out-of-equilibrium phenomena. Without numerical tools, progress in disciplines that require heavy computations with high-dimensional vectors and matrices would be limited to purely analytical predictions and constrained by the available expertise and resources needed to build benchmarking experiments. The goal of this thesis is to advance the classical simulation of quantum many-body systems using Tensor Networks. This technique has proven fruitful for predicting equilibrium ground states and thermal properties of one-dimensional systems, yet it still faces challenges when applied to out-of-equilibrium scenarios. In brief, the dynamics of experimentally accessible expectation values can be represented as a 2D Tensor Network, which must be contracted using iterative methods such as Time-Evolving Block Decimation. However, the memory and computational time required for these calculations—which reflect the methods’ complexity—are known to grow exponentially with the simulated physical time in generic systems. This limitation is often referred to in the literature as the entanglement barrier. This manuscript tackles the problem through the three inequivalent ways in which the 2D Tensor Network can be contracted. In Chapter 3, we study how Schrödinger-picture contraction from the state edge of the network can be modified to improve the convergence of iterative contraction and reduce truncation errors in small systems. To this end, we explore physical decoherence in optimized bases and artificial generalized decoherence, finding that both alternatives show promise for new protocols that overcome the entanglement barrier. In Chapter 4, we explore Heisenberg-picture contraction from the observable edge, combining the concept of generalized decoherence with Pauli Weight truncation. By thoroughly characterizing the Pauli weight requirements to accurately describe the time evolution of local expectation values from initial product states, we conclude that truncating off-diagonal Pauli Weights is an effective scheme for simulating long-time dynamics—and we successfully implement it. We conclude with a spatial characterization of Pauli Weight spreading, bridging Tensor Network simulations and the field of operator hydrodynamics. Finally, in Chapter 5, we focus on the transverse contraction of the 2D Tensor Net-work using the light-cone algorithm, analyzing its variants in terms of truncation strategies and convergence across different models. This leads to a systematic characterization of the method based on the integrability of the underlying Hamiltonians and single-spin noise models.
This thesis develops advanced Tensor Network (TN) methods to address Hamiltonian Lattice Gauge Theories (LGTs), overcoming limitations in real-time dynamics and finite-density regimes. A novel dressed-site formalism is introduced, enabling efficient truncation of gauge fields while preserving gauge invariance for both Abelian and non-Abelian theories. This formalism is successfully applied to SU(2) Yang-Mills LGTs in two dimensions, providing the first TN simulations of this system and revealing critical aspects of its phase diagram and non-equilibrium behavior, such as a Quantum Many-Body (QMB) scarring dynamics. A generalization of the dressed-site formalism is proposed through a new fermion-to-qubit mapping for general lattice fermion theories, revealing powerful for classical and quantum simulations. Numerical innovations, including the use of optimal space-filling curves such as the Hilbert curve to preserve locality in high-dimensional simulations, further enhance the efficiency of these methods. Together with high-performance computing techniques, these advances open current and future development pathways toward optimized, efficient, and faster simulations on scales comparable to Monte Carlo state-of-the-art.
Abstract Tensor network methods as presented in our open source matrix product states software have opened up the possibility to study many-body quantum physics in one and quasi-one-dimensional systems in an easily accessible package similar to density functional theory codes but for strongly correlated dynamics. Here, we address methods which allow one to capture the full entanglement without truncation of the Hilbert space. Such methods are suitable for validation of and comparisons to tensor network algorithms, but especially useful in the case of new kinds of quantum states with high entanglement violating the truncation in tensor networks. Quantum cellular automata are one example for such a system, characterized by tunable complexity, entanglement, and a large spread over the Hilbert space. Beyond the evolution of pure states as a closed system, we adapt the techniques for open quantum systems simulated via the Lindblad master equation. We present three algorithms for solving closed-system many-body time evolution without truncation of the Hilbert space. Exact diagonalization methods have the advantage that they not only keep the full entanglement but also require no approximations to the propagator. Seeking the limits of a maximal number of qubits on a single core, we use Trotter decompositions or Krylov approximation to the exponential of the Hamiltonian. All three methods are also implemented for open systems represented via the Lindblad master equation built from local channels. We show their convergence parameters and focus on efficient schemes for their implementations including Abelian symmetries, e.g. symmetry used for number conservation in the Bose–Hubbard model or discrete symmetries in the quantum Ising model. We present the thermalization timescale in the long-range quantum Ising model as a key example of how exact diagonalization contributes to novel physics.
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 devices, which have a limited number of qubits and short coherence times. We validate our approach by running simulations on a noisy-intermediate scale quantum device, showcasing strong agreement with theoretical predictions. We highlight the potential of HT for simulating QFTs on quantum hardware. Published by the American Physical Society 2024
With the aim of studying nonperturbative out-of-equilibrium dynamics of high-energy particle collisions on quantum simulators, we investigate the scattering dynamics of lattice quantum electrodynamics in 1+1 dimensions. Working in the bosonized formulation of the model and in the thermodynamic limit, we use uniform-matrix-product-state tensor networks to construct multiparticle wave-packet states, evolve them in time, and detect outgoing particles post collision. This facilitates the numerical simulation of scattering experiments in both confined and deconfined regimes of the model at different energies, giving rise to rich phenomenology, including inelastic production of quark and meson states, meson disintegration, and dynamical string formation and breaking. We obtain elastic and inelastic scattering cross sections, together with time-resolved momentum and position distributions of the outgoing particles. Furthermore, we propose an analog circuit-QED implementation of the scattering process that is native to the platform, requires minimal ingredients and approximations, and enables practical schemes for particle wave-packet preparation and evolution. This study highlights the role of classical and quantum simulation in enhancing our understanding of scattering processes in quantum field theories in real time.