Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1002.1824 · Physical Review Letters

Continuous Matrix Product States for Quantum Fields

Frank Verstraete, J. I. Cirac

Published May 14, 2010• 270 citations
arXivHTMLPDF

Abstract

We define matrix product states in the continuum limit, without any reference to an underlying lattice parameter. This allows us to extend the density matrix renormalization group and variational matrix product state formalism to quantum field theories and continuum models in 1 spatial dimension. We illustrate our procedure with the Lieb-Liniger model.

DOI
10.1103/physrevlett.104.190405

Related papers

Continuous Matrix Product States for Quantum Fields: An Energy Minimization Algorithm

Martin Ganahl, Julián Rincón, Guifré Vidal

The generalization of matrix product states (MPS) to continuous systems, as proposed in the breakthrough Letter of Verstraete and Cirac [Phys. Rev. Lett. 104, 190405 (2010).PRLTAO0031-900710.1103/PhysRevLett.104.190405], provides a powerful variational ansatz for the ground state of strongly interacting quantum field theories in one spatial dimension. A continuous MPS (cMPS) approximation to the ground state can be obtained by simulating a Euclidean time evolution. In this Letter we propose a cMPS optimization algorithm based instead on energy minimization by gradient methods and demonstrate its performance by applying it to the Lieb-Liniger model (an integrable model of an interacting bosonic field) directly in the thermodynamic limit. We observe a very significant computational speed-up, of more than 2 orders of magnitude, with respect to simulating a Euclidean time evolution. As a result, a much larger cMPS bond dimension D can be reached (e.g., D=256 with moderate computational resources), thus helping unlock the full potential of the cMPS representation for ground state studies.

arXiv:1611.03779•Physical Review Letters•Jun 2, 2017•53 citations

Finite-Element Matrix Product States for Continuum Models in One Dimension

Akshay Shankar, Karel Van Acoleyen, Jutho Haegeman

We present a matrix product state framework for simulating one-dimensional quantum many-body systems in the continuum using non-orthogonal single-particle basis sets. By mapping the physical problem to an auxiliary computational space, we show that the resulting many-body overlap operator can be efficiently encoded as a matrix product operator for sufficiently localized orbitals, thereby generalizing a construction that first appeared in [arXiv:2405.10285]. This construction recasts the variational ground-state search into a generalized eigenvalue problem, which can be solved using a generalized density matrix renormalization group algorithm. As a primary application, we employ a first-order finite-element expansion to study the ground state properties of the Lieb-Liniger gas in the presence of inhomogeneities. This approach also provides a natural setting for exactly refining the lattice, thereby enabling multigrid optimization strategies for matrix product states.

quant-pharXiv:2606.14873•Jun 16, 2026•0 citations

Continuous matrix product operators for quantum fields

Erickson Tjoa, J. Ignacio Cirac

In this work we introduce an ansatz for continuous matrix product operators for quantum field theory. We show that (i) they admit a closed-form expression in terms of finite number of matrix-valued functions without reference to any lattice parameter; (ii) they are obtained as a suitable continuum limit of matrix product operators; (iii) they preserve the entanglement area law directly in the continuum, and in particular they map continuous matrix product states (cMPS) to another cMPS. As an application, we use this ansatz to construct several families of continuous matrix product unitaries beyond quantum cellular automata.

quant-pharXiv:2511.04545•Physical Review Research•Nov 7, 2025•0 citations

Continuum Limits of Quantum Lattice Systems

Tobias J. Osborne

We describe a general procedure to give effective continuous descriptions of quantum lattice systems in terms of quantum fields. There are two key novelties of our method: firstly, it is framed in the hamiltonian setting and applies equally to distinguishable quantum spins, bosons, and fermions and, secondly, it works for arbitrary variational tensor network states and can easily produce computable non-gaussian quantum field states. Our construction extends the mean-field fluctuation formalism of Hepp and Lieb (developed later by Verbeure and coworkers) to identify emergent continuous large-scale degrees of freedom - the continuous degrees of freedom are not identified beforehand. We apply the construction to to tensor network states, including, matrix product states and projected entangled-pair states, where we recover their recently introduced continuous counterparts, and also for tree tensor networks and the multi-scale entanglement renormalisation ansatz. Finally, extending the continuum limit to include dynamics we obtain a strict light cone for the propagation of information.

arXiv:1901.06124•Jan 18, 2019•9 citations

Variational Optimization of Continuous Matrix Product States

Benoît Tuybens, Jacopo De Nardis, Jutho Haegeman +1 more

Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case, however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.

arXiv:2006.01801•Physical Review Letters•Jan 10, 2022•17 citations

Density Matrix Renormalization Group for Continuous Quantum Systems

Shovan Dutta, Anton S. Buyskikh, Andrew J. Daley +1 more

We introduce a versatile and practical framework for applying matrix product state techniques to continuous quantum systems. We divide space into multiple segments and generate continuous basis functions for the many-body state in each segment. By combining this mapping with existing numerical density matrix renormalization group routines, we show how one can accurately obtain the ground-state wave function, spatial correlations, and spatial entanglement entropy directly in the continuum. For a prototypical mesoscopic system of strongly interacting bosons we demonstrate faster convergence than standard grid-based discretization. We illustrate the power of our approach by studying a superfluid-insulator transition in an external potential. We outline how one can directly apply or generalize this technique to a wide variety of experimentally relevant problems across condensed matter physics and quantum field theory.

arXiv:2108.05366•Physical Review Letters•Jun 8, 2022•9 citations
We define matrix product states in the continuum limit, without any reference to an underlying lattice parameter. This allows us to extend the density matrix renormalization group and variational matrix product state formalism to quantum field theories and continuum models in 1 spatial dimension. We illustrate our procedure with the Lieb-Liniger model.
The generalization of matrix product states (MPS) to continuous systems, as proposed in the breakthrough Letter of Verstraete and Cirac [Phys. Rev. Lett. 104, 190405 (2010).PRLTAO0031-900710.1103/PhysRevLett.104.190405], provides a powerful variational ansatz for the ground state of strongly interacting quantum field theories in one spatial dimension. A continuous MPS (cMPS) approximation to the ground state can be obtained by simulating a Euclidean time evolution. In this Letter we propose a cMPS optimization algorithm based instead on energy minimization by gradient methods and demonstrate its performance by applying it to the Lieb-Liniger model (an integrable model of an interacting bosonic field) directly in the thermodynamic limit. We observe a very significant computational speed-up, of more than 2 orders of magnitude, with respect to simulating a Euclidean time evolution. As a result, a much larger cMPS bond dimension D can be reached (e.g., D=256 with moderate computational resources), thus helping unlock the full potential of the cMPS representation for ground state studies.
We present a matrix product state framework for simulating one-dimensional quantum many-body systems in the continuum using non-orthogonal single-particle basis sets. By mapping the physical problem to an auxiliary computational space, we show that the resulting many-body overlap operator can be efficiently encoded as a matrix product operator for sufficiently localized orbitals, thereby generalizing a construction that first appeared in [arXiv:2405.10285]. This construction recasts the variational ground-state search into a generalized eigenvalue problem, which can be solved using a generalized density matrix renormalization group algorithm. As a primary application, we employ a first-order finite-element expansion to study the ground state properties of the Lieb-Liniger gas in the presence of inhomogeneities. This approach also provides a natural setting for exactly refining the lattice, thereby enabling multigrid optimization strategies for matrix product states.
In this work we introduce an ansatz for continuous matrix product operators for quantum field theory. We show that (i) they admit a closed-form expression in terms of finite number of matrix-valued functions without reference to any lattice parameter; (ii) they are obtained as a suitable continuum limit of matrix product operators; (iii) they preserve the entanglement area law directly in the continuum, and in particular they map continuous matrix product states (cMPS) to another cMPS. As an application, we use this ansatz to construct several families of continuous matrix product unitaries beyond quantum cellular automata.
We describe a general procedure to give effective continuous descriptions of quantum lattice systems in terms of quantum fields. There are two key novelties of our method: firstly, it is framed in the hamiltonian setting and applies equally to distinguishable quantum spins, bosons, and fermions and, secondly, it works for arbitrary variational tensor network states and can easily produce computable non-gaussian quantum field states. Our construction extends the mean-field fluctuation formalism of Hepp and Lieb (developed later by Verbeure and coworkers) to identify emergent continuous large-scale degrees of freedom - the continuous degrees of freedom are not identified beforehand. We apply the construction to to tensor network states, including, matrix product states and projected entangled-pair states, where we recover their recently introduced continuous counterparts, and also for tree tensor networks and the multi-scale entanglement renormalisation ansatz. Finally, extending the continuum limit to include dynamics we obtain a strict light cone for the propagation of information.
Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case, however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.
We introduce a versatile and practical framework for applying matrix product state techniques to continuous quantum systems. We divide space into multiple segments and generate continuous basis functions for the many-body state in each segment. By combining this mapping with existing numerical density matrix renormalization group routines, we show how one can accurately obtain the ground-state wave function, spatial correlations, and spatial entanglement entropy directly in the continuum. For a prototypical mesoscopic system of strongly interacting bosons we demonstrate faster convergence than standard grid-based discretization. We illustrate the power of our approach by studying a superfluid-insulator transition in an external potential. We outline how one can directly apply or generalize this technique to a wide variety of experimentally relevant problems across condensed matter physics and quantum field theory.