Site Basis Excitation Ansatz for Matrix Product States
Pith reviewed 2026-05-18 17:32 UTC · model grok-4.3
The pith
Site basis excitation ansatz computes high-accuracy one-magnon dispersion for the S=1 Heisenberg chain from a single diagonalization step.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The author claims that an infinite matrix product state ground state, prepared with a simple finite-system DMRG procedure, combined with a site-basis excitation ansatz formed by a single multi-state diagonalization, produces the full one-magnon dispersion relation through repeated diagonalization of small non-orthogonal overlap and Hamiltonian matrices. Leaving the basis non-orthogonal proves essential for convergence, while imposing a left-orthonormal gauge degrades performance. The construction further permits exact recovery of all momentum modes from a single localized Wannier excitation translated across sites.
What carries the argument
The site basis excitation ansatz, which generates a compact non-orthogonal set of excitation tensors from one DMRG-like diagonalization to enable efficient momentum-space spectrum calculation.
If this is right
- The one-magnon dispersion of the S=1 Heisenberg chain is obtained efficiently and with high accuracy.
- Not imposing a gauge condition on the excitation basis is crucial for the method to converge well.
- Wannier excitations can be built that, when translated, exactly reconstruct the plane-wave magnon modes at every momentum.
- The simple finite DMRG approach serves as an effective alternative for obtaining the infinite MPS ground state.
Where Pith is reading between the lines
- If the approach generalizes, it could reduce the cost of computing excitations in other strongly correlated 1D systems.
- The non-orthogonal basis choice may inspire similar simplifications in higher-dimensional tensor network methods.
- Localized Wannier excitations open a path to studying real-space properties of magnons more directly.
Load-bearing premise
The infinite matrix product state must accurately represent the true ground state, and the single diagonalization must generate an excitation basis adequate for the desired accuracy without additional refinement.
What would settle it
A direct comparison showing that the SBEA-computed dispersion curve for the S=1 Heisenberg chain differs substantially from well-established numerical benchmarks at key momenta would indicate the method falls short.
Figures
read the original abstract
We introduce a simple and efficient variation of the tangent-space excitation ansatz used to compute elementary excitation spectra of one-dimensional quantum lattice systems using matrix product states (MPS). A small basis for the excitation tensors is formed based on a single diagonalization analogous to a single site DMRG step but for multiple states. Once overlap and Hamiltonian matrix elements are found, obtaining the excitation for any momentum only requires diagonalization of a tiny matrix, akin to a non-orthogonal band-theory diagonalization. The approach is based on an infinite MPS description of the ground state, and we introduce an extremely simple alternative to variational uniform matrix product states (VUMPS) based on finite system DMRG. For the $S=1$ Heisenberg chain, our method -- site basis excitation ansatz (SBEA) -- efficiently produces the one-magnon dispersion with high accuracy. We also examine the role of MPS gauge choices, finding that not imposing a gauge condition -- leaving the basis nonorthogonal -- is crucial for the approach, whereas imposing a left-orthonormal gauge (as in prior work) severely hampers convergence. We also show how one can construct Wannier excitations, analogous to the Wannier functions of band theory, where one Wannier excitation, translated to all sites, can reconstruct the single magnon modes exactly for all momenta.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces the Site Basis Excitation Ansatz (SBEA), a variation of the tangent-space excitation ansatz for elementary excitations in 1D quantum lattice systems with matrix product states. A small non-orthogonal basis for excitation tensors is constructed via a single diagonalization step analogous to a multi-state DMRG sweep; momentum-dependent excitations then follow from a small generalized eigenvalue problem. The ground state is represented by an infinite MPS obtained from a simple finite-system DMRG procedure presented as an alternative to VUMPS. For the S=1 Heisenberg antiferromagnet the method is shown to reproduce the one-magnon dispersion with high accuracy; the authors further demonstrate that omitting any left-orthonormal gauge condition is essential for convergence and that the resulting modes admit an exact Wannier reconstruction.
Significance. If the reported accuracy is confirmed by systematic benchmarks, SBEA would constitute a lightweight, easily implementable route to excitation spectra that avoids both full tangent-space optimization and the more involved VUMPS machinery. The explicit demonstration that a non-orthogonal basis is required, together with the Wannier-function construction, supplies concrete technical insight that could be reused in other MPS-based excitation calculations.
major comments (2)
- [Numerical results] Numerical results section: the central claim of “high accuracy” reproduction of the one-magnon dispersion is stated without quantitative error tables, maximum deviations from literature values, or convergence data versus bond dimension and basis size; such tables are required to substantiate the accuracy assertion and to rule out post-hoc parameter adjustment.
- [Ground-state construction] Ground-state construction paragraph: the finite-system DMRG procedure used to generate the iMPS is described only qualitatively; a direct comparison of its energy or correlation length against a converged VUMPS reference for the same S=1 Heisenberg chain is needed, because any systematic error in the ground state directly propagates into the excitation spectrum.
minor comments (2)
- [Abstract] The abstract and introduction should state the bond dimension and the dimension of the excitation basis explicitly rather than referring only to “small basis.”
- [Figures] Figure captions for the dispersion plots should include the precise literature reference curve against which the SBEA data are compared.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation and the specific suggestions that will help strengthen the manuscript. We address each major comment below.
read point-by-point responses
-
Referee: [Numerical results] Numerical results section: the central claim of “high accuracy” reproduction of the one-magnon dispersion is stated without quantitative error tables, maximum deviations from literature values, or convergence data versus bond dimension and basis size; such tables are required to substantiate the accuracy assertion and to rule out post-hoc parameter adjustment.
Authors: We agree that quantitative benchmarks are needed to support the accuracy claim. In the revised manuscript we will add tables reporting the maximum deviation of our one-magnon energies from established literature values over the Brillouin zone, together with convergence data for the excitation energies as functions of MPS bond dimension and site-basis size. These additions will allow readers to assess the precision systematically. revision: yes
-
Referee: [Ground-state construction] Ground-state construction paragraph: the finite-system DMRG procedure used to generate the iMPS is described only qualitatively; a direct comparison of its energy or correlation length against a converged VUMPS reference for the same S=1 Heisenberg chain is needed, because any systematic error in the ground state directly propagates into the excitation spectrum.
Authors: We acknowledge that a quantitative comparison would improve the presentation. We will compute and include a direct comparison of the ground-state energy per site and correlation length obtained from our finite-system DMRG procedure against a converged VUMPS reference for the S=1 Heisenberg chain. This will confirm that the ground state is sufficiently accurate for the subsequent excitation calculations. revision: yes
Circularity Check
No significant circularity; derivation self-contained
full rationale
The paper presents SBEA as a variation on the standard tangent-space excitation ansatz for infinite MPS, with the ground state obtained from a simple finite-system DMRG procedure explicitly positioned as an alternative to VUMPS and excitations generated by a single diagonalization to form a non-orthogonal basis followed by a small generalized eigenvalue problem. Accuracy for the S=1 Heisenberg one-magnon dispersion is reported by direct comparison to established external literature benchmarks rather than by fitting or self-referential prediction on the same data. Gauge-choice tests and Wannier reconstruction are performed internally with explicit convergence demonstrations. No load-bearing step reduces by construction to a fitted input, self-definition, or unverified self-citation chain; the method remains falsifiable against independent numerical or experimental results.
Axiom & Free-Parameter Ledger
free parameters (1)
- MPS bond dimension
axioms (1)
- domain assumption An infinite MPS ground state can be obtained via the paper's finite-system DMRG procedure as a practical alternative to VUMPS.
Forward citations
Cited by 2 Pith papers
-
Preparation and detection of quasiparticles for quantum simulations of scattering
A method using dressed creation operators from MLWFs enables selective preparation and detection of quasiparticles in lattice theories, tested via MPS on hardcore QCD ladders to separate known excitations from resonances.
-
Checkerboard Bose Hubbard Ladders using Transmon Arrays
Sublattice bias in the checkerboard Bose-Hubbard model makes the commensurate superfluid phase experimentally accessible via transmon arrays and supplies new probes for superfluid and insulating phases.
Reference graph
Works this paper leans on
-
[1]
S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett.69, 2863 (1992)
work page 1992
-
[2]
S. R. White, Density-matrix algorithms for quantum renormalization groups, Phys. Rev. B48, 10345 (1993)
work page 1993
-
[3]
E. Jeckelmann, Dynamical density-matrix renormalization-group method, Physical Review B 66, 045114 (2002)
work page 2002
-
[4]
Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Phys
G. Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Phys. Rev. Lett.93, 040502 (2004)
work page 2004
-
[5]
A. J. Daley, C. Kollath, U. Schollw¨ ock, and G. Vidal, Time-dependent density-matrix renormalization-group using adaptive effective hilbert spaces, J. Stat. Mech. , P04005 (2004)
work page 2004
-
[6]
S. R. White and A. E. Feiguin, Real-time evolution us- ing the density matrix renormalization group, Phys. Rev. Lett.93, 076401 (2004)
work page 2004
-
[7]
J. Haegeman, C. Lubich, I. Oseledets, B. Vandereycken, and F. Verstraete, Unifying time evolution and optimiza- tion with matrix product states, Phys. Rev. B94, 165116 (2016)
work page 2016
-
[8]
U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011)
work page 2011
-
[9]
S. Paeckel, T. K¨ ohler, A. Swoboda, S. R. Manmana, U. Schollw¨ ock, and C. Hubig, Time-evolution methods for matrix-product states, Annals of Physics411, 167998 (2019)
work page 2019
-
[10]
S. R. White and I. Affleck, Spectral function for the $s=1$heisenberg antiferromagnetic chain, Phys. Rev. B 77, 134437 (2008)
work page 2008
-
[11]
J. Haegeman, B. Pirvu, D. Weir, J. I. Cirac, T. J. Os- borne, H. Verschelde, and F. Verstraete, Variational ma- trix product ansatz for dispersion relations, Phys. Rev. B85, 100408 (2012)
work page 2012
-
[12]
L. Vanderstraeten, J. Haegeman, and F. Verstraete, An introduction to tensor networks for lattice field theory, SciPost Phys. Lect. Notes7, 1 (2019)
work page 2019
-
[13]
V. Zauner-Stauber, L. Vanderstraeten, M. T. Fishman, F. Verstraete, and J. Haegeman, Variational optimization algorithms for uniform matrix product states, Phys. Rev. B97, 045145 (2018)
work page 2018
-
[14]
F. D. M. Haldane, Nonlinear field theory of large- spin heisenberg antiferromagnets: Semiclassically quan- tized solitons of the one-dimensional easy-axis n´ eel state, Physics Letters A93, 464 (1983)
work page 1983
-
[15]
F. D. M. Haldane, Continuum dynamics of the 1d heisen- berg antiferromagnet: Identification with the o(3) non- linear sigma model, Physical Review Letters50, 1153 (1983)
work page 1983
-
[16]
I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Rig- orous results on valence-bond ground states in antiferro- magnets, Physical Review Letters59, 799 (1987)
work page 1987
-
[17]
Affleck, Quantum spin chains and the haldane gap, Journal of Physics: Condensed Matter1, 3047 (1989)
I. Affleck, Quantum spin chains and the haldane gap, Journal of Physics: Condensed Matter1, 3047 (1989)
work page 1989
-
[18]
S. R. White and D. A. Huse, Numerical renormalization- group study of the spin-1 antiferromagnetic heisenberg chain, Phys. Rev. B48, 3844 (1993). 8
work page 1993
-
[19]
I. P. McCulloch, Infinite size density matrix renormaliza- tion group, revisited (2008), arXiv:0804.2509
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[20]
M. Fishman, S. R. White, and E. M. Stoudenmire, The ITensor Software Library for Tensor Network Calcula- tions, SciPost Phys. Codebases , 4 (2022)
work page 2022
-
[21]
R. Or´ us and G. Vidal, Infinite time-evolving block deci- mation algorithm beyond unitary evolution, Physical Re- view B78, 155117 (2008)
work page 2008
-
[22]
V. Zauner-Stauber, L. Vanderstraeten, J. Haegeman, I. P. McCulloch, and F. Verstraete, Continuum of spinon and holon excitations in one-dimensional systems, Phys. Rev. B97, 235155 (2018)
work page 2018
-
[23]
J. J. Osborne and I. P. McCulloch, Efficient and sys- tematic calculation of arbitrary observables for the ma- trix product state excitation ansatz, Physical Review Re- search7, 023018 (2025)
work page 2025
-
[24]
J. Haegeman, J. I. Cirac, T. J. Osborne, I. Piˇ zorn, H. Ver- schelde, and F. Verstraete, Time-dependent variational principle for quantum lattices, Phys. Rev. Lett.107, 070601 (2011)
work page 2011
-
[25]
L. Vanderstraeten, F. Verstraete, and J. Haegeman, Scat- tering particles in quantum spin chains, Phys. Rev. B92, 125136 (2015)
work page 2015
-
[26]
N. Marzari and D. Vanderbilt, Maximally localized gen- eralized wannier functions for composite energy bands, Phys. Rev. B56, 12847 (1997)
work page 1997
-
[27]
N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Maximally localized wannier functions: Theory and applications, Rev. Mod. Phys.84, 1419 (2012)
work page 2012
-
[28]
M. Van Damme, L. Vanderstraeten, J. De Nardis, J. Haegeman, and F. Verstraete, Real-time scattering of interacting quasiparticles in quantum spin chains, Phys. Rev. Res.3, 013078 (2021)
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.