REVIEW 5 minor 36 references
TeNeS-v2: Enhancement for Real-Time and Finite Temperature Simulations of Quantum Many-Body Systems
T0 review · 0 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read TeNeS-v2 extends a tensor-network solver to real-time and finite-temperature simulations of two-dimensional quantum lattices.
desk verdict A solid, honest software paper: the real-time and finite-temperature modes work on the demonstrated examples, and the QMC benchmarks give the finite-T claim real teeth. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the iTPS/iPEPS tensor network and its mixed-state generalization, the iTPO (infinite Tensor Product Operator), whose bond dimension $D$ controls accuracy. For finite temperature, the density matrix is evolved as $\rho(2\tau) = e^{-\tau H}\rho(0)e^{-\tau H}$ from $\beta = 0$, with simple-update truncation, and expectation values are contracted as a single-layer square-lattice tensor network via CTMRG with environment bond dimension $\chi$. This single-layer structure is what makes the thermal calculation computationally lighter than the double-layer contraction used for pure-state expectation values.
What would settle it
Run the finite-temperature mode on a small 2D spin model whose exact thermal density matrix is known, and check whether the computed energy density ever falls below the exact ground-state energy or whether a positive observable such as $\langle (S^z)^2\rangle$ becomes negative. If such unphysical values persist as the CTMRG bond dimension $\chi$ is increased without increasing the iTPO bond dimension $D$, then the single-layer iTPO representation itself, not just the contraction, is failing.
Extended reading notes
Core claim
The central claim is that TeNeS-v2 reliably performs real-time and finite-temperature simulations of infinite two-dimensional quantum lattice systems using iTPS/iPEPS. Real-time evolution is obtained by replacing the imaginary-time step $\tau$ with $-it/\hbar$ in the same Suzuki-Trotter update procedure used for ground states, then measuring observables at chosen time intervals. Finite temperature is obtained by imaginary-time evolving a single-layer infinite Tensor Product Operator (iTPO) starting from the infinite-temperature identity state, with CTMRG supplying the environment; because the traced iTPO is a single-layer network rather than a double layer, the thermal calculation is cheaper than a pure-state calculation at the same bond dimension. The paper demonstrates finite-temperature energy density, specific heat, and magnetization of the transverse-field Ising model against quantum Monte Carlo, and the real-time magnetization of a quenched Ising model, while noting that growing entanglement limits real-time accuracy at fixed bond dimension.
Load-bearing premise
The finite-temperature mode assumes that a single-layer iTPO, truncated by simple update and contracted by CTMRG, can represent the thermal density matrix closely enough that computed expectation values stay physical; the paper itself notes positive semidefiniteness is not guaranteed and shows an unphysical singularity in a $D=4$ triangular-lattice calculation.
Editorial extensions
If this is right
- Users can run ground-state, real-time, and finite-temperature calculations on the same 2D lattice models by changing one mode parameter, so dynamical and thermal studies no longer require separate codes.
- For the transverse-field Ising model, thermal energy density, specific heat, and magnetization computed with $D=6$ or $D=10$ match quantum Monte Carlo results, so the method is a viable alternative where QMC sign problems appear.
- Multi-site observable support makes quantities like scalar chirality directly accessible, enabling tensor-network study of non-coplanar magnetic phases on frustrated lattices.
- Real-time simulations at small bond dimension can show artifacts such as energy jumps, but increasing $D$ removes them, giving users a concrete convergence check for dynamics.
- Real-time accuracy at fixed bond dimension is limited to short times because entanglement grows during evolution, matching known behavior of iTPS/iPEPS dynamics.
Reading between the lines
- The single-layer iTPO approach could plausibly be applied to frustrated magnets or other 2D models where quantum Monte Carlo faces sign problems, but the positive-semidefiniteness caveat means users should validate results with energy bounds or independent methods.
- A natural extension would be to add a purified, double-layer iTPO mode alongside the supported single-layer mode, since the paper notes the purified representation avoids unphysical densities at higher computational cost.
- The observed $D=4$ scalar-chirality singularity in the triangular-lattice XXZ model suggests that the simple-update environment, rather than the iTPO representation alone, may be the limiting factor; testing against exact diagonalization on finite clusters would isolate the cause.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. TeNeS-v2 is a software paper for the open-source tensor-network solver TeNeS, which targets infinite two-dimensional quantum lattice systems using iTPS/iPEPS. The new version adds real-time evolution and finite-temperature simulation modes to the previously implemented ground-state solver. The real-time mode applies Suzuki-Trotter decomposition with simple and full updates and evaluates local observables along the trajectory. The finite-temperature mode represents the thermal density matrix as a single-layer iTPO, evolves it in imaginary time from the infinite-temperature identity state using the simple update, and contracts expectation values with CTMRG. The paper also introduces multi-site observables, exemplified by scalar chirality on the triangular lattice. Validation consists of finite-temperature energy, specific heat, and magnetization for the transverse-field Ising model compared with DSQSS quantum Monte Carlo, plus a real-time quench example checked through bond-dimension convergence and energy conservation. The manuscript explicitly discloses known limitations: iTPO positive semidefiniteness is not guaranteed, full update is not yet supported for iTPO, real-time simulations are limited by entanglement growth, and the triangular-lattice example shows an unphysical singularity at D=4 attributed to the simple update.
Significance. If the claims hold, TeNeS-v2 is a useful community resource: it provides an accessible, open-source implementation of real-time and finite-temperature iTPS/iPEPS simulations, with reproducible example data and scripts in an ISSP repository. The finite-temperature mode is benchmarked against an independent method, DSQSS QMC, and the agreement at D=6 and D=10 for the transverse-field Ising model is genuine evidence that the implementation works for that model. The complexity analysis (single-layer contraction with O(D^6) cost when chi is proportional to D) is also helpful for users. The central software-capability claim is therefore supported for the demonstrated cases. The main weakness is the breadth of the validation: only one spin model is checked at finite temperature, no Bose-model benchmark is presented, and the real-time example relies on internal consistency checks rather than an external baseline. These caveats are, however, disclosed in the paper itself, including the explicit warning in Sec. 3.2 that iTPO may violate positive semidefiniteness. I do not see a hidden technical flaw, but the abstract and conclusion should be aligned with the narrower demonstrated scope.
minor comments (5)
- [Abstract and Sec. 6] The phrase 'demonstrating TeNeS-v2's applicability to various quantum spin and Bose models' is broader than the evidence presented: the finite-temperature benchmark covers only the transverse-field Ising model and no Bose model is used in any new-feature example. Please either add an additional benchmark (for example, a hardcore-boson or Heisenberg model at finite temperature) or qualify the scope in the abstract and conclusion, for instance by saying 'spin and boson models with short-range interactions' while noting that the demonstrated examples are spin models.
- [Sec. 5.2, footnote 2] The footnote reporting that for hx=0 a large deviation occurs when D increases above D=2 is potentially confusing and is stated without explanation. Since the main text recommends checking bond-dimension convergence, this example shows that larger D does not always improve accuracy. Please add a short explanation or a forward reference, and give concrete advice on how a user should choose D for a new model.
- [Sec. 5.1, Figs. 4 and 5] The real-time demonstration uses bond-dimension dependence and energy conservation as accuracy diagnostics, but no quantitative tolerance is given. Reporting the maximum energy drift or the norm error after truncation would make the example more reproducible and would clarify what 'the discontinuity disappears' means quantitatively.
- [Sec. 4.2, Eq. (9)] The statement that the step size tau in the imaginary-time evolution operator corresponds to 2tau in the actual inverse-temperature step is correct but could be made more explicit. I suggest writing the cumulative relation beta = 2 N tau for N evolution steps and stating clearly in the FT output-file description that the first column is beta, not the number of steps.
- [Program Summary and references] There are several minor presentation issues: 'di fficult' appears in the Program Summary; the rendering 'T eNeS' is inconsistent in the abstract and body; reference [21] contains an incorrectly rendered exponent 'm = 1 2'; and in Fig. 4 the panel labels (a) and (b) are placed at inconsistent positions. These do not affect the technical content.
Circularity Check
No significant circularity; the new-mode software claims are supported by an independent QMC benchmark and the key approximation limitations are disclosed.
full rationale
The paper's central claim is a software-capability claim: TeNeS-v2 implements real-time evolution and finite-temperature iTPO/iTPS simulations. The algorithmic derivation in Sec. 3 is standard Suzuki-Trotter time evolution with simple/full updates and CTMRG contraction; no parameter is fitted to the results that are then reported as validations. The finite-temperature benchmark in Sec. 5.2 compares TeNeS-v2 output against QMC results from DSQSS. Although DSQSS is cited from overlapping authors, it is an independent stochastic method with its own published implementation, so the comparison is an external check rather than an input to the algorithm. The paper also explicitly discloses the known limitation that iTPO does not guarantee positive semidefiniteness of the density matrix (Sec. 3.2) and shows an unphysical chirality singularity at D=4 (Sec. 5.3), attributing it to the simple-update approximation; this is a disclosed limitation, not a circular derivation. The abstract's phrase 'various quantum spin and Bose models' is broader than the demonstrated benchmarks: only one spin model is compared against QMC and no Bose-model example is shown. That is an overstatement of scope, but it is not circularity. No load-bearing step reduces by definition to its own input, and no fitted input is renamed as a prediction.
Assumptions & free parameters
free parameters (3)
- Bond dimension D of iTPS/iTPO =
D = 3, 6, 10 in examples
- CTMRG bond dimension chi =
chi = 16 in Sec. 5.2
- Trotter time step tau =
e.g., tau = [0.01, 0.005, 0.05] schedules
assumptions (4)
- standard math Suzuki-Trotter decomposition is valid for the small time steps used in real- and imaginary-time evolution.
- domain assumption The iTPS/iTPO tensor network with finite bond dimension D can represent the states of interest well enough for the targeted accuracy.
- domain assumption CTMRG and simple-update contractions give sufficiently accurate expectation values for the single-layer iTPO network.
- domain assumption The simulated models (transverse-field Ising, triangular-lattice XXZ) are representative of the claimed range of applicability.
Cite this review
Pith. "Pith review of TeNeS-v2: Enhancement for Real-Time and Finite Temperature Simulations of Quantum Many-Body Systems." pith.science (2026). https://pith.science/paper/7U6ZFYAD
@misc{pith2026250107777,
author = {Pith},
title = {Pith review of: TeNeS-v2: Enhancement for Real-Time and Finite Temperature Simulations of Quantum Many-Body Systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/7U6ZFYAD}},
note = {Machine review of arXiv:2501.07777}
}
read the original abstract
Quantum many-body systems are challenging targets for computational physics due to their large degrees of freedom. The tensor networks, particularly Tensor Product States (TPS) and Projected Entangled Pair States (PEPS), effectively represent these systems on two-dimensional lattices. However, the technical complexity of TPS/PEPS-based coding is often too much for researchers to handle. To reduce this problem, we developed TeNeS (Tensor Network Solver). This paper introduces TeNeS-v2, which extends TeNeS with real-time and finite temperature simulations, providing deeper insights into quantum many-body systems. We detail the new algorithms, input/output design, and application examples, demonstrating TeNeS-v2's applicability to various quantum spin and Bose models on two-dimensional lattices.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
R. Or ´us, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics 349 (2014) 117–158. doi:https://doi.org/10. 1016/j.aop.2014.06.013
work page 2014
-
[2]
Or ´us, Tensor networks for complex quantum systems, Na- ture Reviews Physics 1 (9) (2019) 538–550
R. Or ´us, Tensor networks for complex quantum systems, Na- ture Reviews Physics 1 (9) (2019) 538–550. doi:10.1038/ s42254-019-0086-7
work page 2019
-
[3]
J. Eisert, M. Cramer, M. B. Plenio, Colloquium: Area laws for the entanglement entropy, Rev. Mod. Phys. 82 (2010) 277–306. doi:10.1103/RevModPhys.82.277
-
[4]
T. Nishino, Y . Hieida, K. Okunishi, N. Maeshima, Y . Akutsu, A. Gendiar, Two-Dimensional Tensor Product Variational For- mulation, Progress of Theoretical Physics 105 (3) (2001) 409–
work page 2001
-
[5]
F. Verstraete, J. I. Cirac, Renormalization algorithms for quantum-many body systems in two and higher dimensions (2004). arXiv:cond-mat/0407066
arXiv 2004
-
[6]
P. Corboz, F. Mila, Tensor network study of the Shastry- Sutherland model in zero magnetic field, Phys. Rev. B 87 (2013) 115144. doi:10.1103/PhysRevB.87.115144
-
[7]
T. Okubo, K. Shinjo, Y . Yamaji, N. Kawashima, S. Sota, T. To- hyama, M. Imada, Ground-state properties of na 2iro3 deter- mined from an ab initio Hamiltonian and its extensions contain- ing Kitaev and extended Heisenberg interactions, Phys. Rev. B 96 (2017) 054434. doi:10.1103/PhysRevB.96.054434
-
[8]
H. J. Liao, Z. Y . Xie, J. Chen, Z. Y . Liu, H. D. Xie, R. Z. Huang, B. Normand, T. Xiang, Gapless spin-liquid ground state in the s = 1/2 kagome antiferromagnet, Phys. Rev. Lett. 118 (2017) 137202. doi:10.1103/PhysRevLett.118.137202
Show all 36 references
-
[9]
H.-Y . Lee, R. Kaneko, T. Okubo, N. Kawashima, Gapless kitaev spin liquid to classical string gas through tensor networks, Phys. Rev. Lett. 123 (2019) 087203. doi:10.1103/PhysRevLett. 123.087203
2019 doi
-
[10]
H.-Y . Lee, N. Kawashima, Y . B. Kim, Tensor network wave function of s =1 kitaev spin liquids, Physical Review Research 2 (3) (Aug. 2020). doi:10.1103/physrevresearch.2. 033318
2020 doi
-
[11]
Mashiko, T
T. Mashiko, T. Okubo, Quantum phase transition between spin liquid and spin nematics in spin-1 kitaev honeycomb model, Phys. Rev. Res. 6 (2024) 033110. doi:10.1103/ PhysRevResearch.6.033110
2024
-
[12]
Corboz, T
P. Corboz, T. M. Rice, M. Troyer, Competing states in the t- j model: Uniform d-wave state versus stripe state, Phys. Rev. Lett. 113 (2014) 046402. doi:10.1103/PhysRevLett.113. 046402
2014 doi
-
[13]
Corboz, F
P. Corboz, F. Mila, Crystals of bound states in the magnetization plateaus of the shastry-sutherland model, Phys. Rev. Lett. 112 (2014) 147203. doi:10.1103/PhysRevLett.112.147203
2014 doi
-
[14]
Czarnik, J
P. Czarnik, J. Dziarmaga, P. Corboz, Time evolution of an infi- nite projected entangled pair state: An efficient algorithm, Phys- ical Review B 99 (3) (Jan. 2019). doi:10.1103/physrevb. 99.035115
2019 doi
-
[15]
Hubig, A
C. Hubig, A. Bohrdt, M. Knap, F. Grusdt, J. I. Cirac, Eval- uation of time-dependent correlators after a local quench in iPEPS: hole motion in the t-J model, SciPost Phys. 8 (2020)
2020
-
[16]
Kaneko, I
R. Kaneko, I. Danshita, Tensor-network study of correlation- spreading dynamics in the two-dimensional bose-hubbard model, Communications Physics 5 (1) (2022) 65. doi:10. 1038/s42005-022-00848-9
2022
-
[17]
R. T. Ponnaganti, M. Mambrini, D. Poilblanc, Tensor network variational optimizations for real-time dynamics: Application to 9 the time-evolution of spin liquids, SciPost Phys. 15 (2023) 158. doi:10.21468/SciPostPhys.15.4.158
2023 doi
-
[18]
Ponsioen, P
B. Ponsioen, P. Corboz, Excitations with projected entangled pair states using the corner transfer matrix method, Phys. Rev. B 101 (19) (2020) 195109. arXiv:2001.02645, doi:10.1103/ PhysRevB.101.195109
2020 arXiv
-
[19]
W.-L. Tu, L. Vanderstraeten, N. Schuch, H.-Y . Lee, N. Kawashima, J.-Y . Chen, Generating function for projected entangled-pair states, PRX Quantum 5 (2024) 010335. doi: 10.1103/PRXQuantum.5.010335
2024 doi
-
[20]
Kshetrimayum, M
A. Kshetrimayum, M. Rizzi, J. Eisert, R. Or ´us, Tensor net- work annealing algorithm for two-dimensional thermal states, Physical Review Letters 122 (7) (Feb. 2019). doi:10.1103/ physrevlett.122.070502
2019
-
[21]
doi:10.21468/SciPostPhys.8.2.021
-
[22]
J. L. Jim ´enez, S. P. G. Crone, E. Fogh, M. E. Zayed, R. Lortz, E. Pomjakushina, K. Conder, A. M. L¨auchli, L. Weber, S. Wes- sel, A. Honecker, B. Normand, C. R ¨uegg, P. Corboz, H. M. Rønnow, F. Mila, A quantum magnetic analogue to the criti- cal point of water, Nature 592 (...
2021 doi
-
[23]
Czarnik, M
P. Czarnik, M. M. Rams, P. Corboz, J. Dziarmaga, Tensor net- work study of the m = 1 2 magnetization plateau in the shastry- sutherland model at finite temperature, Phys. Rev. B 103 (2021) 075113. doi:10.1103/PhysRevB.103.075113
2021 doi
-
[24]
https://www.pasums.issp.u-tokyo.ac.jp/tenes/en
-
[25]
https://github.com/issp-center-dev/TeNeS
-
[26]
Or ´us, G
R. Or ´us, G. Vidal, Simulation of two-dimensional quantum systems on an infinite lattice revisited: Corner transfer matrix for tensor contraction, Physical Review B 80 (9) (Sep. 2009). doi:10.1103/physrevb.80.094403
2009 doi
-
[27]
Motoyama, T
Y . Motoyama, T. Okubo, K. Yoshimi, S. Morita, T. Kato, N. Kawashima, TeNeS: Tensor network solver for quantum lat- tice systems, Computer Physics Communications 279 (2022) 108437. doi:10.1016/j.cpc.2022.108437
2022
-
[28]
Jordan, R
J. Jordan, R. Or ´us, G. Vidal, F. Verstraete, J. I. Cirac, Classical simulation of infinite-size quantum lattice systems in two spa- tial dimensions, Physical Review Letters 101 (25) (Dec. 2008). doi:10.1103/physrevlett.101.250602
2008 doi
-
[29]
H. C. Jiang, Z. Y . Weng, T. Xiang, Accurate determination of tensor network state of quantum lattice models in two dimen- sions, Physical Review Letters 101 (9) (Aug. 2008). doi: 10.1103/physrevlett.101.090603
2008 doi
-
[30]
https://isspns-gitlab.issp.u-tokyo.ac.jp/ tenes-dev/tenes-gallery
-
[31]
H. N. Phien, J. A. Bengua, H. D. Tuan, P. Corboz, R. Or ´us, Infinite projected entangled pair states algorithm improved: Fast full update and gauge fixing, Physical Review B 92 (3) (Jul. 2015). doi:10.1103/physrevb.92.035142
2015 doi
-
[32]
Yamamoto, G
D. Yamamoto, G. Marmorini, I. Danshita, Quantum phase dia- gram of the triangular-lattice XXZ model in a magnetic field, Phys. Rev. Lett. 112 (12) (2014) 127203. doi:10.1103/ PhysRevLett.112.127203
2014
-
[33]
Motoyama, K
Y . Motoyama, K. Yoshimi, A. Masaki-Kato, T. Kato, N. Kawashima, Dsqss: Discrete space quantum systems solver, Computer Physics Communications 264 (2021) 107944. doi: https://doi.org/10.1016/j.cpc.2021.107944
2021
-
[34]
Motoyama, K
Y . Motoyama, K. Yoshimi, T. Kato, S. Todo, Materiapps live! and materiapps installer: Environment for starting and scaling up materials science simulations, SoftwareX 20 (2022) 101210. doi:https://doi.org/10.1016/j.softx.2022.101210. 10
2022
-
[35]
Sellmann, X.-F
D. Sellmann, X.-F. Zhang, S. Eggert, Phase diagram of the antiferromagnetic XXZ model on the triangular lattice, Phys. Rev. B Condens. Matter 91 (8) (2015) 081104. doi:10.1103/ PhysRevB.91.081104
2015
-
[417]
doi:10.1143/PTP.105.409
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.