pith. sign in

arxiv: 2509.06241 · v2 · submitted 2025-09-07 · ❄️ cond-mat.str-el · quant-ph

Site Basis Excitation Ansatz for Matrix Product States

Pith reviewed 2026-05-18 17:32 UTC · model grok-4.3

classification ❄️ cond-mat.str-el quant-ph
keywords matrix product statesexcitation ansatzHeisenberg chainmagnon dispersionquantum spin modelsDMRGone-dimensional systems
0
0 comments X

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.

The paper presents the site basis excitation ansatz as a simplified way to calculate elementary excitations in one-dimensional quantum lattice models described by matrix product states. It forms a small basis of excitation tensors through one diagonalization step similar to a DMRG update on the infinite ground state. This basis then allows the momentum-dependent energies to be found by diagonalizing small matrices for each wavevector. The method is demonstrated on the spin-one Heisenberg chain where it reproduces the magnon dispersion curve accurately. Sympathetic readers would see this as a practical advance that reduces the computational effort needed to study spin-wave spectra compared to full tangent-space optimizations.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2509.06241 by Steven R. White.

Figure 2
Figure 2. Figure 2: FIG. 2. Spin profile of single site excitations [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. One-magnon dispersion of the spin-1 Heisenberg [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Spin profile of a Wannier excitation for the Heisen [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.”
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

1 free parameters · 1 axioms · 0 invented entities

The approach rests on standard MPS representational assumptions and the new SBEA construction; no new physical particles or forces are introduced.

free parameters (1)
  • MPS bond dimension
    Numerical accuracy depends on the chosen bond dimension, a standard tunable parameter in MPS calculations.
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.
    The method is built on an infinite MPS description of the ground state obtained this way.

pith-pipeline@v0.9.0 · 5761 in / 1407 out tokens · 48496 ms · 2026-05-18T17:32:42.732730+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Preparation and detection of quasiparticles for quantum simulations of scattering

    quant-ph 2026-04 unverdicted novelty 6.0

    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.

  2. Checkerboard Bose Hubbard Ladders using Transmon Arrays

    cond-mat.other 2026-05 unverdicted novelty 4.0

    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

28 extracted references · 28 canonical work pages · cited by 2 Pith papers · 1 internal anchor

  1. [1]

    S. R. White, Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett.69, 2863 (1992)

  2. [2]

    S. R. White, Density-matrix algorithms for quantum renormalization groups, Phys. Rev. B48, 10345 (1993)

  3. [3]

    Jeckelmann, Dynamical density-matrix renormalization-group method, Physical Review B 66, 045114 (2002)

    E. Jeckelmann, Dynamical density-matrix renormalization-group method, Physical Review B 66, 045114 (2002)

  4. [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)

  5. [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)

  6. [6]

    S. R. White and A. E. Feiguin, Real-time evolution us- ing the density matrix renormalization group, Phys. Rev. Lett.93, 076401 (2004)

  7. [7]

    Haegeman, C

    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)

  8. [8]

    Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011)

    U. Schollw¨ ock, The density-matrix renormalization group in the age of matrix product states, Annals of Physics 326, 96 (2011)

  9. [9]

    Paeckel, T

    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)

  10. [10]

    S. R. White and I. Affleck, Spectral function for the $s=1$heisenberg antiferromagnetic chain, Phys. Rev. B 77, 134437 (2008)

  11. [11]

    Haegeman, B

    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)

  12. [12]

    Vanderstraeten, J

    L. Vanderstraeten, J. Haegeman, and F. Verstraete, An introduction to tensor networks for lattice field theory, SciPost Phys. Lect. Notes7, 1 (2019)

  13. [13]

    Zauner-Stauber, L

    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)

  14. [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)

  15. [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)

  16. [16]

    Affleck, T

    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)

  17. [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)

  18. [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

  19. [19]

    I. P. McCulloch, Infinite size density matrix renormaliza- tion group, revisited (2008), arXiv:0804.2509

  20. [20]

    Fishman, S

    M. Fishman, S. R. White, and E. M. Stoudenmire, The ITensor Software Library for Tensor Network Calcula- tions, SciPost Phys. Codebases , 4 (2022)

  21. [21]

    Or´ us and G

    R. Or´ us and G. Vidal, Infinite time-evolving block deci- mation algorithm beyond unitary evolution, Physical Re- view B78, 155117 (2008)

  22. [22]

    Zauner-Stauber, L

    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)

  23. [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)

  24. [24]

    Haegeman, J

    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)

  25. [25]

    Vanderstraeten, F

    L. Vanderstraeten, F. Verstraete, and J. Haegeman, Scat- tering particles in quantum spin chains, Phys. Rev. B92, 125136 (2015)

  26. [26]

    Marzari and D

    N. Marzari and D. Vanderbilt, Maximally localized gen- eralized wannier functions for composite energy bands, Phys. Rev. B56, 12847 (1997)

  27. [27]

    Marzari, A

    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)

  28. [28]

    Van Damme, L

    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)