The infinite time-evolving block decimation algorithm (iTEBD) of Vidal (2007) approximates ground states of infinite one-dimensional quantum lattice systems by imaginary-time evolution of a matrix product state followed by truncation to a fixed bond dimension, but its validity was justified only for finite, gapped chains, through an overlap bound that becomes vacuous when the initial ground-state overlap vanishes in the thermodynamic limit or when the spectral gap closes. No rigorous convergence theorem or first-principles error analysis is known for this algorithm class in the gapless and/or thermodynamic regimes. The problem asks for that theory: convergence of the truncated non-unitary-gate iteration where the original bound fails, quantitative bounds for the non-unitary small-time-step errors, and a rigorous derivation of the observed finite-entanglement scaling that relates bond dimension, induced correlation length, and the accuracy of long-distance observables and the entanglement spectrum at criticality. Long-range interactions, time-dependent Hamiltonians, and infinite two-dimensional tensor networks are outside the scope of this problem.

Background

Ground states of one-dimensional (1D) quantum lattice systems are the central object of classical many-body simulation. The time-evolving block decimation algorithm (TEBD) of Vidal (2004) represents the time-evolved wavefunction as a matrix product state (MPS) — an ansatz in which the coefficients of a state on nn sites are encoded in a chain of tensors contracted over auxiliary indices, with the maximum index size χ\chi (the bond dimension) bounding how much entanglement the ansatz can carry. Evolution under a nearest-neighbor Hamiltonian H=∑rh[r,r+1]H=\sum_r h[r,r+1] is approximated by a Suzuki–Trotter splitting — a product expansion of e−Hte^{-Ht} into small two-site gates e−ih δte^{-ih\,\delta t} with controlled error in the time step δt\delta t — and after each two-site gate the Schmidt decomposition (the bipartite spectral decomposition ∣Ψ⟩=∑αλα∣Φα⟩⊗∣Φα′⟩|\Psi\rangle=\sum_\alpha \lambda_\alpha |\Phi_\alpha\rangle\otimes|\Phi'_\alpha\rangle) is truncated to the χ\chi largest weights. The infinite TEBD (iTEBD) algorithm of the source paper (Vidal 2007) exploits translation invariance so that the infinite chain is simulated directly with a fixed pair of tensors, and was benchmarked against the exactly solvable transverse-field Ising chain.

Ground states are targeted by imaginary-time projection ∣Ψτ⟩∝e−Hτ∣Ψ0⟩|\Psi_\tau\rangle \propto e^{-H\tau}|\Psi_0\rangle with τ→∞\tau\to\infty. For a finite chain whose Hamiltonian has a spectral gap Δ>0\Delta>0 between the ground state and the rest of the spectrum, this converges: the source's Eq. (17) bounds 1−∣⟨Ψτ∣Ψg⟩∣≲e−2Δτ/δ21-|\langle\Psi_\tau|\Psi_g\rangle| \lesssim e^{-2\Delta\tau}/\delta^2, where δ≡∣⟨Ψ0∣Ψg⟩∣\delta\equiv|\langle\Psi_0|\Psi_g\rangle| is the initial ground-state overlap. The paper states explicitly that this bound may become meaningless in the limit n→∞n\to\infty (where typically δ→0\delta\to0) or for a gapless Hamiltonian (Δ=0\Delta=0), declines to justify the method's validity in those regimes, and refers instead to benchmark calculations. Two further gaps are acknowledged: imaginary time requires non-unitary gates e−h δτe^{-h\,\delta\tau}, which destroy the Schmidt (canonical) structure on which strictly local, parallelizable tensor updates rest, with the induced errors only heuristically controlled by taking δτ→0\delta\tau\to0; and at the critical point of the Ising benchmark the half-chain entanglement spectrum {pα}\{p_\alpha\} (the eigenvalues of the reduced density matrix of half the chain, pα=λα2p_\alpha=\lambda_\alpha^2) fails to converge under the iteration, while two-point correlators remain accurate over thousands of lattice sites — leaving the quantitative relation between χ\chi, achievable correlation length, and observable errors open.

The physics of the critical regime is understood: at a quantum critical point the half-chain entanglement entropy (the von Neumann entropy of the reduced density matrix) grows as S∼(c/6)log⁡ξS\sim(c/6)\log\xi with the correlation length ξ\xi, where cc is the central charge — the effective number of massless degrees of freedom of the conformal field theory describing the continuum limit (Calabrese and Cardy 2004) — so no finite-χ\chi MPS can represent the exact critical ground state, whose ξ=∞\xi=\infty. Subsequent work built the quantitative χ\chi-dependence into the finite-entanglement scaling framework: ξ(χ)∼χκ\xi(\chi)\sim\chi^\kappa with κ\kappa fixed by the central charge (Tagliacozzo et al. 2008; Pollmann et al. 2009), later refined into effective-field-theory and boundary-conformal descriptions. Rigorous results cover the other side of the problem: for gapped 1D systems, area-law theorems (bounded entanglement of bounded regions) and polynomial-time algorithms establish that ground states are approximable by MPS (Hastings 2007; Landau–Vazirani–Vidick 2015; Arad et al. 2017), and local observables of critical states are approximable by MPS whose bond dimension grows polynomially in the inverse target accuracy (Huang 2019). These are existence statements, not convergence theorems for the truncated imaginary-time iteration, and they do not cover long-distance observables or the entanglement spectrum. Successor algorithms — the time-dependent variational principle (TDVP), which evolves within the tangent space of the MPS manifold without Trotter gates, and variational uniform MPS (VUMPS) — remove some error sources and converge faster in practice but likewise come without convergence proofs at criticality.

Problem Statement

Provide the rigorous foundation that the source paper explicitly declined to supply: an a priori convergence and error theory for the truncated imaginary-time matrix product state iteration of the TEBD/iTEBD class, valid in the regimes n→∞n\to\infty and/or Δ=0\Delta=0 in which the source's only convergence argument, Eq. (17), is vacuous. The problem has three linked components, which together constitute one error analysis of one algorithm: component (2) is the error input to the global convergence analysis (1) — the per-step Trotter and truncation errors together with their accumulation over unboundedly many steps — and component (3) is the equilibrium finite-χ\chi error floor of the critical regime that survives after (1) and (2) are controlled, so that (2) and (3) cover disjoint error content.

(1) Convergence. Formulate precisely and prove (or disprove) convergence of the iTEBD iteration — repeated imaginary-time Suzuki–Trotter layers of non-unitary two-site gates e−h δτe^{-h\,\delta\tau}, each followed by truncation of the two-site Schmidt decomposition to bond dimension χ\chi, applied to a suitable initial MPS (translation-invariant in the infinite-chain regimes), in the limits δτ→0\delta\tau\to0 and τ→∞\tau\to\infty — as a ground-state approximation method, in each regime where Eq. (17) fails: (i) the gapped thermodynamic limit, (ii) gapless finite chains, and (iii) the gapless thermodynamic limit. Convergence should be characterized in a suitable local sense (e.g. expectation values of local observables, or energy density), since the global overlap δ\delta vanishes as n→∞n\to\infty, together with rates in δτ\delta\tau, χ\chi, and the model parameters.

(2) Local-update error analysis. Derive quantitative bounds for the error introduced by the non-unitary gate update, which destroys the Schmidt (canonical) form on which strictly local parallelized tensor updates rest, as a function of δτ\delta\tau, χ\chi, and the model parameters — including how small δτ\delta\tau must be chosen at fixed χ\chi. The analysis must separate and jointly control Trotter error, truncation error, and their accumulation over unboundedly many steps.

(3) Finite-entanglement scaling at criticality. Derive from first principles — from properties of the Hamiltonian and its criticality class, not from numerical fits — the relation between χ\chi, the induced correlation length ξ(χ)\xi(\chi) (empirically ξ∼χκ\xi\sim\chi^\kappa with κ\kappa governed by the central charge, approximately κ=6/(c(1+12/c))\kappa=6/(c(1+\sqrt{12/c}))), the achievable accuracy of long-distance observables at distances r≲ξ(χ)r\lesssim\xi(\chi), and the truncation error of the half-chain entanglement spectrum; or establish rigorous conditions under which such scaling laws fail.

Numerical corroboration on exactly solvable models (such as the transverse-field Ising chain on which the source paper was benchmarked) supports an answer but does not by itself constitute one.

References

  1. G. Vidal, 'Classical simulation of infinite-size quantum lattice systems in one spatial dimension', Phys. Rev. Lett. 98, 070201 (2007). DOI: 10.1103/PhysRevLett.98.070201; arXiv:cond-mat/0605597; https://arxiv.org/abs/cond-mat/0605597 (source paper)
  2. G. Vidal, 'Efficient simulation of one-dimensional quantum many-body systems', Phys. Rev. Lett. 93, 040502 (2004). DOI: 10.1103/PhysRevLett.93.040502; arXiv:quant-ph/0310089; https://arxiv.org/abs/quant-ph/0310089
  3. J. Eisert, 'Computational difficulty of global variations in the density matrix renormalization group', Phys. Rev. Lett. 97, 260501 (2006). DOI: 10.1103/PhysRevLett.97.260501; arXiv:quant-ph/0609051; https://arxiv.org/abs/quant-ph/0609051
  4. F. Verstraete and J. I. Cirac, 'Matrix product states represent ground states faithfully', Phys. Rev. B 73, 094423 (2006). DOI: 10.1103/PhysRevB.73.094423; arXiv:cond-mat/0505140; https://arxiv.org/abs/cond-mat/0505140
  5. M. B. Hastings, 'An area law for one-dimensional quantum systems', J. Stat. Mech. (2007) P08024. DOI: 10.1088/1742-5468/2007/08/P08024; arXiv:0705.2024; https://arxiv.org/abs/0705.2024
  6. Z. Landau, U. Vazirani, and T. Vidick, 'A polynomial time algorithm for the ground state of one-dimensional gapped local Hamiltonians', Nature Physics 11, 566-569 (2015). DOI: 10.1038/nphys3345; https://doi.org/10.1038/nphys3345
  7. I. Arad, Z. Landau, U. Vazirani, and T. Vidick, 'Rigorous RG algorithms and area laws for low energy eigenstates in 1D', Commun. Math. Phys. 356, 65-105 (2017). DOI: 10.1007/s00220-017-2973-z; arXiv:1602.08828; https://arxiv.org/abs/1602.08828
  8. A. M. Dalzell and F. G. S. L. Brandão, 'Locally accurate MPS approximations for ground states of one-dimensional gapped local Hamiltonians', Quantum 3, 187 (2019). DOI: 10.22331/q-2019-09-23-187; https://quantum-journal.org/papers/q-2019-09-23-187/
  9. N. Schuch, M. M. Wolf, F. Verstraete, and J. I. Cirac, 'Entropy scaling and simulability by matrix product states', Phys. Rev. Lett. 100, 030504 (2008). DOI: 10.1103/PhysRevLett.100.030504; arXiv:0705.0292; https://arxiv.org/abs/0705.0292
  10. Y. Huang, 'Approximating local properties by tensor network states with constant bond dimension', IEEE Trans. Inf. Theory 72, 4930-4935 (2026). DOI: 10.1109/TIT.2026.3694133; arXiv:1903.10048; https://arxiv.org/abs/1903.10048
  11. L. Tagliacozzo, T. R. de Oliveira, S. Iblisdir, and J. I. Latorre, 'Scaling of entanglement support for matrix product states', Phys. Rev. B 78, 024410 (2008). DOI: 10.1103/PhysRevB.78.024410; arXiv:0712.1976; https://arxiv.org/abs/0712.1976
  12. F. Pollmann, S. Mukerjee, A. M. Turner, and J. E. Moore, 'Theory of finite-entanglement scaling at one-dimensional quantum critical points', Phys. Rev. Lett. 102, 255701 (2009). DOI: 10.1103/PhysRevLett.102.255701; arXiv:0812.2903; https://arxiv.org/abs/0812.2903
  13. B. Pirvu, G. Vidal, F. Verstraete, and L. Tagliacozzo, 'Matrix product states for critical spin chains: finite-size versus finite-entanglement scaling', Phys. Rev. B 86, 075117 (2012). DOI: 10.1103/PhysRevB.86.075117; arXiv:1204.3934; https://arxiv.org/abs/1204.3934
  14. Vid Stojevic, Jutho Haegeman, I. P. McCulloch, Luca Tagliacozzo, and Frank Verstraete, 'Conformal data from finite entanglement scaling', Phys. Rev. B 91, 035120 (2015). DOI: 10.1103/PhysRevB.91.035120; https://doi.org/10.1103/PhysRevB.91.035120
  15. Ho N. Phien, Ian P. McCulloch, and Guifré Vidal, 'Fast convergence of imaginary time evolution tensor network algorithms by recycling the environment', Phys. Rev. B 91, 115137 (2015). DOI: 10.1103/PhysRevB.91.115137; arXiv:1411.0391; https://arxiv.org/abs/1411.0391
  16. J. Haegeman, J. I. Cirac, T. J. Osborne, I. Pižorn, H. Verschelde, and F. Verstraete, 'Time-dependent variational principle for quantum lattices', Phys. Rev. Lett. 107, 070601 (2011). DOI: 10.1103/PhysRevLett.107.070601; arXiv:1103.0936; https://arxiv.org/abs/1103.0936
  17. V. Zauner-Stauber, L. Vanderstraeten, M. T. Fishman, F. Verstraete, and J. Haegeman, 'Variational optimization algorithms for uniform matrix product states', Phys. Rev. B 97, 045145 (2018). DOI: 10.1103/PhysRevB.97.045145; arXiv:1701.07035; https://arxiv.org/abs/1701.07035
  18. Markus Bachmayr, Henrik Eisenmann, Emil Kieri, and André Uschmajew, 'Existence of dynamical low-rank approximations to parabolic problems', Mathematics of Computation 90, 1799-1830 (2021). DOI: 10.1090/mcom/3626; https://doi.org/10.1090/mcom/3626
  19. Markus Bachmayr, Henrik Eisenmann, and André Uschmajew, 'Dynamical low-rank tensor approximations to high-dimensional parabolic problems: existence and convergence of spatial discretizations', Numerische Mathematik 157, 781-822 (2025). DOI: 10.1007/s00211-025-01465-8; https://doi.org/10.1007/s00211-025-01465-8
  20. Y. Zhang and E. Solomonik, 'On stability of tensor networks and canonical forms', arXiv:2001.01191 (2020). https://arxiv.org/abs/2001.01191
  21. A. Ueda and M. Oshikawa, 'Finite-size and finite bond dimension effects of tensor network renormalization', Phys. Rev. B 108, 024413 (2023). DOI: 10.1103/PhysRevB.108.024413; arXiv:2302.06632; https://arxiv.org/abs/2302.06632
  22. N. E. Sherman, A. Avdoshkin, and J. E. Moore, 'Universality of critical dynamics with finite entanglement', Phys. Rev. Lett. 131, 106501 (2023). DOI: 10.1103/PhysRevLett.131.106501; arXiv:2301.09681; https://arxiv.org/abs/2301.09681
  23. Rui-Zhen Huang, Long Zhang, Andreas M. Läuchli, Jutho Haegeman, Frank Verstraete, and Laurens Vanderstraeten, 'Emergent conformal boundaries from finite-entanglement scaling in matrix product states', Phys. Rev. Lett. 132, 086503 (2024). DOI: 10.1103/PhysRevLett.132.086503; arXiv:2306.08163; https://arxiv.org/abs/2306.08163
  24. J. T. Schneider, A. Ueda, Y. Liu, A. M. Läuchli, M. Oshikawa, and L. Tagliacozzo, 'Self-congruent point in critical matrix product states: an effective field theory for finite-entanglement scaling', SciPost Phys. 18, 4.142 (2025). DOI: 10.21468/SciPostPhys.18.4.142; arXiv:2411.03954; https://arxiv.org/abs/2411.03954
  25. P. Calabrese and J. Cardy, 'Entanglement entropy and quantum field theory', J. Stat. Mech. (2004) P06002. DOI: 10.1088/1742-5468/2004/06/P06002; arXiv:hep-th/0405152; https://arxiv.org/abs/hep-th/0405152
  26. U. Schollwöck, 'The density-matrix renormalization group in the age of matrix product states', Annals of Physics 326, 96-192 (2011). DOI: 10.1016/j.aop.2010.09.012; arXiv:1008.3477; https://arxiv.org/abs/1008.3477
  27. J. Jordan, R. Orús, G. Vidal, F. Verstraete, and J. I. Cirac, 'Classical simulation of infinite-size quantum lattice systems in two spatial dimensions', Phys. Rev. Lett. 101, 250602 (2008). DOI: 10.1103/PhysRevLett.101.250602; arXiv:cond-mat/0703788; https://arxiv.org/abs/cond-mat/0703788