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:quant-ph/0310089 · Physical Review Letters

Efficient Simulation of One-Dimensional Quantum Many-Body Systems

Guifré Vidal

Published Jul 19, 2004• 1,771 citations
arXivHTMLPDF

Abstract

We present a numerical method to simulate the time evolution, according to a generic Hamiltonian made of local interactions, of quantum spin chains and systems alike. The efficiency of the scheme depends on the amount of entanglement involved in the simulated evolution. Numerical analysis indicates that this method can be used, for instance, to efficiently compute time-dependent properties of low-energy dynamics in sufficiently regular but otherwise arbitrary one-dimensional quantum many-body systems. As by-products, we describe two alternatives to the density matrix renormalization group method.

DOI
10.1103/physrevlett.93.040502

Related papers

Effcient simulation of the adaptive time-dependent density-matrix renormalization-group with periodic boundary conditions

Jie Ren

We introduce a numerical method of the adaptive time-dependent density-matrix renormalization-group to compute one-dimensional quantum spin systems with periodic boundary condition. We check our algorithm to study the dynamic correlation in spin-1/2 Heisenberg XX chain at zero temperature, and the numerical analysis of errors indicates that this method could be used to efficiently simulate the time-dependent properties of low-energy dynamics in an arbitrary one-dimensional quantum many-body systems with the nearest-neighbor interaction.

arXiv:1105.6023•May 30, 2011•0 citations

Classical simulation of quantum many-body systems

Yichen Huang

Classical simulation of quantum many-body systems is in general a challenging problem for the simple reason that the dimension of the Hilbert space grows exponentially with the system size. In particular, merely encoding a generic quantum many-body state requires an exponential number of bits. However, condensed matter physicists are mostly interested in local Hamiltonians and especially their ground states, which are highly non-generic. Thus, we might hope that at least some physical systems allow efficient classical simulation.Starting with one-dimensional (1D) quantum systems (i.e., the simplest nontrivial case), the first basic question is: Which classes of states have efficient classical representations? It turns out that this question is quantitatively related to the amount of entanglement in the state, for states with ``little entanglement'' are well approximated by matrix product states (a data structure that can be manipulated efficiently on a classical computer). At a technical level, the mathematical notion for ``little entanglement'' is area law, which has been proved for unique ground states in 1D gapped systems. We establish an area law for constant-fold degenerate ground states in 1D gapped systems and thus explain the effectiveness of matrix-product-state methods in (e.g.) symmetry breaking phases. This result might not be intuitively trivial as degenerate ground states in gapped systems can be long-range correlated.Suppose an efficient classical representation exists. How can one find it efficiently? The density matrix renormalization group is the leading numerical method for computing ground states in 1D quantum systems. However, it is a heuristic algorithm and the possibility that it may fail in some cases cannot be completely ruled out. Recently, a provably efficient variant of the density matrix renormalization group has been developed for frustration-free 1D gapped systems. We generalize this algorithm to all (i.e., possibly frustrated) 1D gapped systems. Note that the ground-state energy of 1D gapless Hamiltonians is computationally intractable even in the presence of translational invariance.It is tempting to extend methods and tools in 1D to two and higher dimensions (2+D), e.g., matrix product states are generalized to tensor network states. Since an area law for entanglement (if formulated properly) implies efficient matrix product state representations in 1D, an interesting question is whether a similar implication holds in 2+D. Roughly speaking, we show that an area law for entanglement (in any reasonable formulation) does not always imply efficient tensor network representations of the ground states of 2+D local Hamiltonians even in the presence of translational invariance. It should be emphasized that this result does not contradict with the common sense that in practice quantum states with more entanglement usually require more space to be stored classically; rather, it demonstrates that the relationship between entanglement and efficient classical representations is still far from being well understood.Excited eigenstates participate in the dynamics of quantum systems and are particularly relevant to the phenomenon of many-body localization (absence of transport at finite temperature in strongly correlated systems). We study the entanglement of excited eigenstates in random spin chains and expect that its singularities coincide with dynamical quantum phase transitions. This expectation is confirmed in the disordered quantum Ising chain using both analytical and numerical methods.Finally, we study the problem of generating ground states (possibly with topological order) in 1D gapped systems using quantum circuits. This is an interesting problem both in theory and in practice. It not only characterizes the essential difference between the entanglement patterns that give rise to trivial and nontrivial topological order, but also quantifies the difficulty of preparing quantum states with a quantum computer (in experiments).

•Jan 1, 2015•19 citations

Efficient time-evolution of matrix product states using average Hamiltonians

Belal Abouraya, Jirawat Saiphet, Fedor Jelezko +1 more

Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.

quant-pharXiv:2602.04955•Feb 6, 2026•0 citations

Time evolution of one-dimensional Quantum Many Body Systems

Salvatore R. Manmana

The level of current understanding of the physics of time‐dependent strongly correlated quantum systems is far from complete, principally due to the lack of effective controlled approaches. Recently, there has been progress in the development of approaches for one‐dimensional systems. We describe recent developments in the construction of numerical schemes for general (one‐dimensional) Hamiltonians: in particular, schemes based on exact diagonalization techniques and on the density matrix renormalization group method (DMRG). We present preliminary results for spinless fermions with nearest‐neighbor‐interaction and investigate their accuracy by comparing with exact results.

arXiv:cond-mat/0502396•AIP conference proceedings•Jan 1, 2005•84 citations

Is the efficiency of classical simulations of quantum dynamics related to integrability?

Tomaž Prosen, Marko Žnidarič

Efficiency of time evolution of quantum observables, and thermal states of quenched Hamiltonians, is studied using time-dependent density-matrix renormalization group method in a family of generic quantum spin chains which undergo a transition from integrable to nonintegrable-quantum chaotic case as control parameters are varied. Quantum states (observables) are represented in terms of matrix-product operators with rank D(epsilon)(t), such that evolution of a long chain is accurate within fidelity error epsilon up to time t. We found that the rank generally increases exponentially D(epsilon)(t) proportional to exp(const t), unless the model system was integrable in which case we found polynomial increase.

arXiv:quant-ph/0608057•Physical Review E•Jan 11, 2007•138 citations

Mixed-State Dynamics in One-Dimensional Quantum Lattice Systems: A Time-Dependent Superoperator Renormalization Algorithm

Michael P. Zwolak, Guifré Vidal

We present an algorithm to study mixed-state dynamics in one-dimensional quantum lattice systems. The algorithm can be used, e.g., to construct thermal states or to simulate real time evolution given by a generic master equation. Its two main ingredients are (i) a superoperator renormalization scheme to efficiently describe the state of the system and (ii) the time evolving block decimation technique to efficiently update the state during a time evolution. The computational cost of a simulation increases significantly with the amount of correlations between subsystems, but it otherwise depends only linearly on the system size. We present simulations involving quantum spins and fermions in one spatial dimension.

arXiv:cond-mat/0406440•Physical Review Letters•Nov 12, 2004•655 citations
We present a numerical method to simulate the time evolution, according to a generic Hamiltonian made of local interactions, of quantum spin chains and systems alike. The efficiency of the scheme depends on the amount of entanglement involved in the simulated evolution. Numerical analysis indicates that this method can be used, for instance, to efficiently compute time-dependent properties of low-energy dynamics in sufficiently regular but otherwise arbitrary one-dimensional quantum many-body systems. As by-products, we describe two alternatives to the density matrix renormalization group method.
We introduce a numerical method of the adaptive time-dependent density-matrix renormalization-group to compute one-dimensional quantum spin systems with periodic boundary condition. We check our algorithm to study the dynamic correlation in spin-1/2 Heisenberg XX chain at zero temperature, and the numerical analysis of errors indicates that this method could be used to efficiently simulate the time-dependent properties of low-energy dynamics in an arbitrary one-dimensional quantum many-body systems with the nearest-neighbor interaction.
Classical simulation of quantum many-body systems is in general a challenging problem for the simple reason that the dimension of the Hilbert space grows exponentially with the system size. In particular, merely encoding a generic quantum many-body state requires an exponential number of bits. However, condensed matter physicists are mostly interested in local Hamiltonians and especially their ground states, which are highly non-generic. Thus, we might hope that at least some physical systems allow efficient classical simulation.Starting with one-dimensional (1D) quantum systems (i.e., the simplest nontrivial case), the first basic question is: Which classes of states have efficient classical representations? It turns out that this question is quantitatively related to the amount of entanglement in the state, for states with ``little entanglement'' are well approximated by matrix product states (a data structure that can be manipulated efficiently on a classical computer). At a technical level, the mathematical notion for ``little entanglement'' is area law, which has been proved for unique ground states in 1D gapped systems. We establish an area law for constant-fold degenerate ground states in 1D gapped systems and thus explain the effectiveness of matrix-product-state methods in (e.g.) symmetry breaking phases. This result might not be intuitively trivial as degenerate ground states in gapped systems can be long-range correlated.Suppose an efficient classical representation exists. How can one find it efficiently? The density matrix renormalization group is the leading numerical method for computing ground states in 1D quantum systems. However, it is a heuristic algorithm and the possibility that it may fail in some cases cannot be completely ruled out. Recently, a provably efficient variant of the density matrix renormalization group has been developed for frustration-free 1D gapped systems. We generalize this algorithm to all (i.e., possibly frustrated) 1D gapped systems. Note that the ground-state energy of 1D gapless Hamiltonians is computationally intractable even in the presence of translational invariance.It is tempting to extend methods and tools in 1D to two and higher dimensions (2+D), e.g., matrix product states are generalized to tensor network states. Since an area law for entanglement (if formulated properly) implies efficient matrix product state representations in 1D, an interesting question is whether a similar implication holds in 2+D. Roughly speaking, we show that an area law for entanglement (in any reasonable formulation) does not always imply efficient tensor network representations of the ground states of 2+D local Hamiltonians even in the presence of translational invariance. It should be emphasized that this result does not contradict with the common sense that in practice quantum states with more entanglement usually require more space to be stored classically; rather, it demonstrates that the relationship between entanglement and efficient classical representations is still far from being well understood.Excited eigenstates participate in the dynamics of quantum systems and are particularly relevant to the phenomenon of many-body localization (absence of transport at finite temperature in strongly correlated systems). We study the entanglement of excited eigenstates in random spin chains and expect that its singularities coincide with dynamical quantum phase transitions. This expectation is confirmed in the disordered quantum Ising chain using both analytical and numerical methods.Finally, we study the problem of generating ground states (possibly with topological order) in 1D gapped systems using quantum circuits. This is an interesting problem both in theory and in practice. It not only characterizes the essential difference between the entanglement patterns that give rise to trivial and nontrivial topological order, but also quantifies the difficulty of preparing quantum states with a quantum computer (in experiments).
Simulating quantum many-body systems (QMBS) is one of the long-standing, highly non-trivial challenges in condensed matter physics and quantum information due to the exponentially growing size of the system's Hilbert space. To date, tensor networks have been an essential tool for studying such quantum systems, owing to their ability to efficiently capture the entanglement properties of the systems they represent. One of the well-known tensor network architectures, namely matrix product states (MPS), is the standard method for simulating one-dimensional QMBS. Here, we propose a simple, yet efficient, method to augment the already available MPS algorithms to simulate the dynamics of time-dependent Hamiltonians with better accuracy and a faster convergence rate, giving a second-order convergence compared to the first-order convergence of the standard method. We apply our proposed method to simulate the dynamics of a chain of single spins associated with nitrogen-vacancy color centers in diamonds, which has potential applications for practical and scalable quantum technologies, and find that our method improves the average error for a system of few NV centers by a factor of about 1000 for moderate step sizes. Our work paves the way for efficient simulation of QMBS under the influence of time-dependent Hamiltonians.
The level of current understanding of the physics of time‐dependent strongly correlated quantum systems is far from complete, principally due to the lack of effective controlled approaches. Recently, there has been progress in the development of approaches for one‐dimensional systems. We describe recent developments in the construction of numerical schemes for general (one‐dimensional) Hamiltonians: in particular, schemes based on exact diagonalization techniques and on the density matrix renormalization group method (DMRG). We present preliminary results for spinless fermions with nearest‐neighbor‐interaction and investigate their accuracy by comparing with exact results.
Efficiency of time evolution of quantum observables, and thermal states of quenched Hamiltonians, is studied using time-dependent density-matrix renormalization group method in a family of generic quantum spin chains which undergo a transition from integrable to nonintegrable-quantum chaotic case as control parameters are varied. Quantum states (observables) are represented in terms of matrix-product operators with rank D(epsilon)(t), such that evolution of a long chain is accurate within fidelity error epsilon up to time t. We found that the rank generally increases exponentially D(epsilon)(t) proportional to exp(const t), unless the model system was integrable in which case we found polynomial increase.
We present an algorithm to study mixed-state dynamics in one-dimensional quantum lattice systems. The algorithm can be used, e.g., to construct thermal states or to simulate real time evolution given by a generic master equation. Its two main ingredients are (i) a superoperator renormalization scheme to efficiently describe the state of the system and (ii) the time evolving block decimation technique to efficiently update the state during a time evolution. The computational cost of a simulation increases significantly with the amount of correlations between subsystems, but it otherwise depends only linearly on the system size. We present simulations involving quantum spins and fermions in one spatial dimension.