We obtain a well-behaved continuum limit of projected entangled pair states (PEPS) that provides an abstract class of quantum field states with natural symmetries. Making use of the recently introduced path integral representation of one-dimensional con- tinuous matrix product states (cMPS) for quantum fields, we demonstrate how symmetries of the physical field state are encoded within the dynamics of an auxiliary field system of one dimension less. In particular, the imposition of euclidean symmetries on the physical system requires that the auxiliary system involved in the class' definition must be Lorentz invariant. The physical field states automatically inherit entropy area laws from the PEPS class, and are fully described by the dissipative dynamics of the lower dimensional auxiliary field system.
We argue that the natural way to generalise a tensor network variational class to a continuous quantum system is to use the Feynman path integral to implement a continuous tensor contraction. This approach is illustrated for the case of a recently introduced class of quantum field states known as continuous matrix-product states (cMPS). As an example of the utility of the path-integral representation we argue that the state of a dynamically evolving quantum field admits a natural representation as a cMPS. An argument that all states in Fock space admit a cMPS representation when the number of variational parameters tends to infinity is also provided.
We introduce a new class of states for bosonic quantum fields which extend tensor network states to the continuum and generalize continuous matrix product states (cMPS) to spatial dimensions
d≥2. By construction, they are Euclidean invariant, and are genuine continuum limits of discrete tensor network states. Admitting both a functional integral and an operator representation, they share the important properties of their discrete counterparts: expressiveness, invariance under gauge transformations, simple rescaling flow, and compact expressions for the
N-point functions of local observables. While we discuss mostly the continuous tensor network states extending Projected Entangled Pair States (PEPS), we propose a generalization bearing similarities with the continuum Multi-scale Entanglement Renormalization Ansatz (cMERA).
These notes provide an introduction to tensor-network methods in quantum many-body physics, with an emphasis on matrix-product states (MPS). They develop the basic tensor-network language, including graphical notation, virtual indices, bond dimensions, gauge freedom, canonical forms, QR and singular-value decompositions, and the role of entanglement in controlling the efficiency of the representation. The main MPS algorithms are then introduced, including contractions, correlation functions, matrix-product operators, DMRG, and time-evolution methods. The notes also briefly discuss projected entangled-pair states (PEPS) as a higher-dimensional generalization of MPS, together with the basic ideas behind approximate PEPS contraction. Finally, tensor-network representations of mixed states, quantum channels, and Lindblad dynamics are presented, with applications to thermal states and open quantum systems. The presentation is accompanied by short Julia code examples based on ITensor, ITensorMPS, and TensorMixedStates. These notes were written for the 9th Les Houches Summer School on Computational Physics: Open Quantum Systems, held in June 2026.
We introduce a new class of states for bosonic quantum fields which extend tensor network states to the continuum and generalize continuous matrix product states to spatial dimensions d ≥ 2. By construction, they are Euclidean invariant and are genuine continuum limits of discrete tensor network states.Admitting both a functional integral and an operator representation, they share the important properties of their discrete counterparts: expressiveness, invariance under gauge transformations, simple rescaling flow, and compact expressions for the N-point functions of local observables.While we discuss mostly the continuous tensor network states extending projected entangled-pair states, we propose a generalization bearing similarities with the continuum multiscale entanglement renormalization ansatz.
Tensor network methods are playing a central role in multiple disciplines of modern quantum physics, such as condensed matter physics, quantum information and quantum gravity. Based on renormalization group (RG) ideas, they provide a large class of variational ansatz for quantum many-body systems. There are various tensor network structures which are suited for different systems based on their entanglement structures, such as matrix product states (MPS), projected entangled pair states (PEPS), multi-scale entanglement renormalization ansatz (MERA). Numerous techniques involving coarse-graining tensor networks have been developed, such as tensor renormalization group (TRG), tensor network renormalization (TNR) and its variants. \n \nFirst, we study the generalization of MERA to continuous systems, or cMERA, which is expected to become a powerful variational ansatz for the ground state of strongly interacting quantum field theories. We investigate, in the simpler context of Gaussian cMERA for free theories, the extent to which the cMERA state
ketPsicM with finite UV cut-off
Lambda can capture the spacetime symmetries of the ground state
ketPsi. For a free boson conformal field theory (CFT) in 1+1 dimensions as a concrete example, we build a quasi-local unitary transformation
V that maps
ketPsi into
ketPsicM and show two main results. (i) Any spacetime symmetry of the ground state
ketPsi is also mapped by
V into a spacetime symmetry of the cMERA
ketPsicM. However, while in the CFT the stress-energy tensor
Tmunu(x) (in terms of which all the spacetime symmetry generators are expressed) is local, the corresponding cMERA stress-energy tensor
TmunucM(x)=VTmunu(x)Vdagger is quasi-local. (ii) From the cMERA, we can extract quasi-local scaling operators
OalphacM(x) characterized by the exact same scaling dimensions
Deltaalpha, conformal spins
salpha, operator product expansion coefficients
Calphabetagamma, and central charge
c as the original CFT. We argue that these results should also apply to interacting theories. \n \nSecond we extend TNR to field theories in the continuum. A short-distance length scale
1/Lambda is introduced in the continuum partition function by smearing the fields. The resulting object is still defined in the continuum but has no fluctuations at distances shorter than
1/Lambda. An infinitesimal coarse-graining step is then generated by the combined action of a \\textit{rescaling} operator
mathfrakL and a \\textit{disentangling} operator
mathfrakK that implements a quasilocal field redefinition. As demonstrated for a free boson in two dimensions, continuous TNR exactly preserves translation and rotation symmetries and can generate a proper RG flow. Moreover, from a critical fixed point of this RG flow one can then extract the conformal data of the underlying conformal field theory. \n \nLeaving behind tensor network methods, we explore the universal structure that emerges from the eigenvectors of the reduced density matrix (RDM) for 1+1 dimensional critical lattice systems. Using the Ising model as an example, we demonstrate that the large-weight eigenvectors are related to the low-lying states of the underlying CFT with suitable boundary conditions by a conformal map, and they transform accordingly under the Virasoro generators. \n \nFinally, using the quantum renormalization group (QRG), we derive the bulk geometry that emerges in the holographic dual of the fermionic U(
N) vector model at a nonzero charge density. The obstruction that prohibits the metallic state from being smoothly deformable to the direct product state under the renormalization group flow gives rise to a horizon at a finite radial coordinate in the bulk. The region outside the horizon is described by the Lifshitz geometry with a higher-spin hair determined by microscopic details of the boundary theory. On the other hand, the interior of the horizon is not described by any Riemannian manifold, as it exhibits an algebraic non-locality. The non-local structure inside the horizon carries the information on the shape of the filled Fermi sea.