pith. sign in

arxiv: 2606.01918 · v1 · pith:GRSXZFF4new · submitted 2026-06-01 · 🧮 math-ph · math.MP

Exact spectrum of the XX spin chain with constrained non-diagonal boundary fields

Changqing Liu , Xiaotian Xu , Xin Zhang This is my paper

Pith reviewed 2026-06-28 12:25 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords XX spin chainBethe Ansatznon-diagonal boundary fieldsthermodynamic limitground state energyexact spectrumelementary excitationsintegrable models
0
0 comments X

The pith

The XX spin chain with constrained non-diagonal boundaries has an analytical ground state energy in the thermodynamic limit.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper solves the Bethe Ansatz equations for the XX spin chain under constrained non-diagonal boundary fields. The equations require that the number of Bethe roots has definite parity and that every root lies at a zero of a unary function. Numerical inspection of finite chains identifies the root patterns that realize the ground state and the first excited state. From these patterns the authors extract a closed analytical expression for the ground-state energy once the chain length tends to infinity. A reader would care because the result supplies an exact benchmark for open-boundary integrable magnets that can be checked without diagonalizing ever-larger matrices.

Core claim

The spectrum is obtained exactly from the Bethe Ansatz equations whose roots must satisfy a definite parity condition and lie among the zeros of a unary function. Numerical evidence shows that the ground state and first excited state correspond to particular root configurations, and that elementary excitations arise from the simultaneous displacement of a pair of roots. These configurations yield an explicit closed-form expression for the ground-state energy in the thermodynamic limit.

What carries the argument

Bethe Ansatz equations whose roots are constrained to the zeros of a unary function and obey a fixed parity on their number.

If this is right

  • The ground-state energy is given directly by the closed analytical formula without solving the full Bethe equations for large but finite chains.
  • Elementary excitations always involve the cooperative shift of exactly two Bethe roots.
  • The parity constraint on the total number of roots partitions the spectrum into allowed sectors.
  • Any eigenenergy is determined once the positions of the roots at the unary zeros are known.
  • The thermodynamic-limit energy expression is independent of the particular finite-size root pattern beyond the ground-state configuration.

Where Pith is reading between the lines

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

  • The paired-root excitation rule may simplify the computation of dynamical correlation functions or quench dynamics in the same model.
  • The unary-zero condition could be exploited to obtain exact expressions for boundary observables such as local magnetizations.
  • Similar parity and unary-function structures may appear in related open XXZ chains, allowing analogous thermodynamic energies to be derived.
  • The analytical energy formula supplies a precise target for testing approximate methods such as variational tensor networks on open spin chains.
  • Large-scale numerics on chains of length several hundred would directly confirm or refute the extrapolation from the observed finite-size root patterns.

Load-bearing premise

The Bethe-root configurations identified numerically for the ground and first excited states remain valid when deriving the closed-form thermodynamic-limit energy expression.

What would settle it

Compute the ground-state energy by exact diagonalization or DMRG on chains of several hundred sites, extrapolate to infinite length, and test whether the result matches the proposed analytical formula within the expected 1/L corrections.

Figures

Figures reproduced from arXiv: 2606.01918 by Changqing Liu, Xiaotian Xu, Xin Zhang.

Figure 1
Figure 1. Figure 1: Configuration of zeros of the auxiliary function f(u) in different regions of the boundary parameters w±. The zero configurations in regions B1(D1) and B2(D2) are quite similar, yet no direct relation exists between them. 3.1. Finite N case Let us start with finite-size systems. First, the properties of f(u) at some specific points can be analyzed f(±1) = 0, f(0) = −1, (22) f(w±) = −(w 2 ±g + 1)(1 − w 2 ±)… view at source ↗
Figure 2
Figure 2. Figure 2: Distributions of {u1, . . . , uN+1} for different boundary parameters. The blue circles denote {u1, . . . , uN+1}, the gray dashed line corresponds to Im(u) = 0, and the yellow dashed line represents the semicircle of unit radius. Panels (a) and (b) correspond to region A. The real zeros undergo qualitative changes between N = 13 and N = 14. Panels (c) and (d) correspond to region B. The real zeros in the … view at source ↗
Figure 3
Figure 3. Figure 3: Variation of the real zero with the system size N for the boundary parameters g = 0.04, w+ = 1.75, and w− = 0.36. The red circles denote the numerical values of the zeros, while the blue solid line indicates their evolution trend as N increases. The gray auxiliary line corresponds to uj = w −1 + . Real and imaginary quasi-momentum p have different bare energies, and the imaginary p plays a critical role wh… view at source ↗
Figure 4
Figure 4. Figure 4: hσ z ni versus n for certain eigenstates. Here N = 10, M = 1, β = 0.50, {α−, α+, θ−} = {0.73, 1.00, 0.78}. The blue, red, green, gray, and purple dots represent the eigenstates corresponding to p = {0.0926i}, p = {0.3495}, p = {0.6267}, p = {2.8491}, and p = {1.4501}, respectively. 4 Distributions of Bethe roots in the ground state and first excited state 4.1. Ground state Since all Bethe roots are non-int… view at source ↗
Figure 5
Figure 5. Figure 5: (a): Distribution of the Bethe roots for N = 100, w+ = −1.50, w− = −1.22, g = 0.04. Here, p − 2 = 1.5270, p − 1 = 1.5582, p + 1 = 1.5893, p + 2 = 1.6205, the single quasi-particle energy is ǫ(p − 2 ) = 0.1752, ǫ(p − 1 ) = 0.0506, ǫ(p + 1 ) = −0.0741, ǫ(p + 2 ) = −0.1987. When M is odd, the ground state energy is E(pg) = −129.0713. The elementary excitation has three probabilities (42): (1): replace p + 1 i… view at source ↗
Figure 6
Figure 6. Figure 6: “Phase” diagram of the distribution of the M Bethe roots corresponding to the ground and first excited state for different lattice sizes. In the four panels, regions a1, a2, b1, and b2 each exhibit distinct root distribution patterns for the ground and first excited states. The patterns vary across regions but remain consistent within each region. For a given region, the primary difference for different la… view at source ↗
Figure 7
Figure 7. Figure 7: (a): Bethe roots distributions for the parameters N = 76, w+ = 2.00, w− = 0.50, g = 0.04. Here, p − 2 = 1.5097, p − 1 = 1.5504, p + 1 = 1.5912, p + 2 = 1.6319, the energy corresponding to a single root is ǫ(p − 2 ) = 0.2444, ǫ(p − 1 ) = 0.0815, ǫ(p + 1 ) = −0.0815, ǫ(p + 2 ) = −0.2444. When M = N1 + N2 + 2k + 1, the ground state energy is E(pg,1) = E(pg,2) = −98.4305 and the first excited state energy is E… view at source ↗
Figure 8
Figure 8. Figure 8: The curves of p ± 1 versus the system size N, with fixed boundary parameters w+ = 1.74, w− = 0.45, and g = 0.04. Here the blue circles and red squares denote the numerical data of p − 1 and p + 1 , respectively, while the yellow and purple curves correspond to the fitting functions p − 1 (N) = π/2 − 1.41N−1 + 0.95N−2 and p + 1 (N) = π/2 + 1.60N−1 − 1.07N−2 , respectively. The gray dashed line in the middle… view at source ↗
read the original abstract

We study the exact spectrum of the XX spin chain with constrained non-diagonal boundary fields, which can be analyzed by solving the associated Bethe Ansatz equations. In these equations, the number of Bethe roots has a definite parity, and all Bethe roots are located at the zeros of a unary function. We investigate the possible positions of the Bethe roots. Based on numerical observations, we analyze the Bethe root configurations for the ground state and the first excited state. Our results show that elementary excitations are characterized by the cooperative change of a pair of Bethe roots. Furthermore, we obtain an analytical expression for the ground state energy in the thermodynamic limit.

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

1 major / 1 minor

Summary. The manuscript analyzes the Bethe Ansatz equations for the XX spin chain with constrained non-diagonal boundary fields. It reports that the number of Bethe roots has definite parity and all roots lie at zeros of a unary function. Based on numerical observations, the authors identify root configurations for the ground state and first excited state, observe that elementary excitations involve cooperative shifts of a pair of roots, and derive a closed-form analytical expression for the ground-state energy in the thermodynamic limit.

Significance. An exact closed-form expression for the thermodynamic-limit ground-state energy of this boundary-value problem would be a useful addition to the literature on integrable spin chains with non-diagonal boundaries, provided the underlying root configurations are rigorously justified.

major comments (1)
  1. [Abstract] Abstract and the derivation of the thermodynamic-limit energy: the closed-form expression is obtained by taking the continuum limit of a sum over Bethe roots whose configuration (all roots real, fixed parity, located at zeros of the unary function) is identified numerically for finite N. No analytic argument is supplied showing that this pattern is the global energy minimizer for every N, that no admissible alternative set of roots (different cardinality or complex roots) yields lower energy, or that the pattern survives without rearrangement as N→∞. This assumption is load-bearing for the central claim.
minor comments (1)
  1. [Abstract] The phrase 'unary function' in the abstract is nonstandard; clarify whether 'univariate' is intended.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for their careful reading and for identifying the central assumption underlying the thermodynamic-limit result. We respond to the major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract and the derivation of the thermodynamic-limit energy: the closed-form expression is obtained by taking the continuum limit of a sum over Bethe roots whose configuration (all roots real, fixed parity, located at zeros of the unary function) is identified numerically for finite N. No analytic argument is supplied showing that this pattern is the global energy minimizer for every N, that no admissible alternative set of roots (different cardinality or complex roots) yields lower energy, or that the pattern survives without rearrangement as N→∞. This assumption is load-bearing for the central claim.

    Authors: We agree that the closed-form expression rests on the root configuration identified numerically. The manuscript already states that the configurations are determined from numerical observations. Extensive checks for system sizes up to several hundred sites confirm that the ground-state roots are real, obey the fixed parity, and occupy zeros of the unary function, while other admissible sets (complex roots or altered cardinality) produce higher energies when compared with exact diagonalization for small N. Nevertheless, the manuscript supplies no analytic demonstration that the observed pattern is the unique global minimizer for every N or that it persists without rearrangement as N→∞. Such a proof lies outside the scope of the present analysis. revision: no

standing simulated objections not resolved
  • Absence of an analytic proof that the numerically observed Bethe-root configuration is the global energy minimizer for all N and survives unchanged to the thermodynamic limit.

Circularity Check

0 steps flagged

No circularity: thermodynamic-limit expression derived from conjectured root patterns without reduction to fitted inputs or self-citations.

full rationale

The paper identifies Bethe-root configurations via numerical observation and then derives a closed-form thermodynamic-limit ground-state energy from the resulting sum over those roots. This is a standard conjecture-then-derive procedure rather than any of the enumerated circular patterns: no parameter is fitted and renamed as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and the final expression is not definitionally equivalent to the numerical input. The dependence on the observed pattern is an assumption about which configuration is minimal, but the derivation itself does not collapse to that assumption by construction. Hence the central claim remains non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on the standard Bethe-Ansatz framework for the XX chain; the abstract supplies no explicit free parameters, invented entities, or ad-hoc axioms beyond those of the method itself.

axioms (1)
  • domain assumption The transfer-matrix or algebraic Bethe-Ansatz construction remains valid under the stated boundary constraints.
    Invoked implicitly when the authors state that the spectrum is obtained by solving the associated Bethe Ansatz equations.

pith-pipeline@v0.9.1-grok · 5634 in / 1167 out tokens · 25903 ms · 2026-06-28T12:25:47.734896+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

41 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    Heisenberg W 1928 Zeitschrift f¨ ur Physik49 619–636

    1928

  2. [2]

    Mourigal M, Enderle M, Kl¨ opperpieper A, Caux J S, Stunault A and Rønnow H M 2013 Nature Physics 9 435–441

    2013

  3. [3]

    Mikeska H J and Kolezhuk A K 2008 One-dimensional magnetism Quantum magnetism (Springer) pp 1–83

    2008

  4. [4]

    Liu W Y, Gong S S, Li Y B, Poilblanc D, Chen W Q and Gu Z C 2022 Science bulletin 67 1034–1041

    2022

  5. [5]

    Broholm C, Cava R J, Kivelson S, Nocera D, Norman M and Senthil T 2 020 Science 367 eaay0668

  6. [6]

    Jiang H C 2021 npj Quantum Materials 6 71

    2021

  7. [7]

    Baxter R J 1982 Exactly Solved Models in Statistical Mechanics (Academic Press)

    1982

  8. [8]

    Korepin V E, Korepin V E, Bogoliubov N and Izergin A 1997 Quantum Inverse Scattering Method and Correlation Functions (Cambridge University Press)

    1997

  9. [9]

    Wang Y, Yang W L, Cao J and Shi K 2015 Off-diagonal Bethe ansatz for exactly solvable models (Springer)

    2015

  10. [10]

    Bethe H 1931 Z. Phys. 71 205

    1931

  11. [11]

    Baxter R J 2002 J. Statist. Phys. 108 1–48

    2002

  12. [12]

    Nepomechie R I 2003 J. Phys. A 37 433–440

    2003

  13. [13]

    Nepomechie R I and Ravanini F 2003 J. Phys. A 36 11391

    2003

  14. [14]

    Yang W L and Zhang Y Z 2006 Nucl. Phys. B 744 312–329

    2006

  15. [15]

    Takhtadzhan L and Faddeev L D 1979 Russ. Math. Surv. 34 11–68

    1979

  16. [16]

    Cao J, Lin H Q, Shi K J and Wang Y 2003 Nucl. Phys. B 663 487–519

    2003

  17. [17]

    Belliard S and Cramp´ e N 2013 SIGMA 9 072

    2013

  18. [18]

    Takahashi M 1999 Thermodynamics of one-dimensional solvable models (Cambridge University Press, Cam- bridge)

    1999

  19. [19]

    Cao J, Yang W L, Shi K and Wang Y 2013 Phys. Rev. Lett. 111 137201

    2013

  20. [20]

    Cao J, Yang W L, Shi K and Wang Y 2013 Nucl. Phys. B 875 152–165

    2013

  21. [21]

    Cao J, Yang W L, Shi K and Wang Y 2013 Nucl. Phys. B 877 152–175

    2013

  22. [22]

    Niccoli G 2012 J. Stat. Mech. 2012 P10025

    2012

  23. [23]

    Niccoli G 2013 Nucl. Phys. B 870 397–420

    2013

  24. [24]

    Faldella S, Kitanine N and Niccoli G 2014 J. Stat. Mech. 1401 P01011

    2014

  25. [25]

    Zhang X, Kl¨ umper A and Popkov V 2021 Phys. Rev. B 103 115435 21

    2021

  26. [26]

    Zhang X, Kl¨ umper A and Popkov V 2021 Phys. Rev. B 104 195409

    2021

  27. [27]

    Sklyanin E K 1988 J. Phys. A 21 2375–2389

    1988

  28. [28]

    Cherednik I V 1984 Theor. Math. Phys 61 977–983

    1984

  29. [29]

    Li Y Y, Cao J, Yang W L, Shi K and Wang Y 2014 Nucl. Phys. B 884 17–27

    2014

  30. [30]

    Sun P, Xin Z R, Qiao Y, Hao K, Cao L, Cao J, Yang T and Yang W L 2019 J. Phys. A 52 265201

    2019

  31. [31]

    Qiao Y, Cao J, Yang W L, Shi K and Wang Y 2021 Phys. Rev. B 103 220401

    2021

  32. [32]

    Dong J S, Lu P, Sun P, Qiao Y, Cao J, Hao K and Yang W L 2023 Chinese Physics B 32 017501

    2023

  33. [33]

    Kozlowski K and Pozsgay B 2012 Journal of Statistical Mechanics: Theory and Experiment 2012 P05021

    2012

  34. [34]

    Pozsgay B and R´ akos O 2018 Journal of Statistical Mechanics: Theory and Experiment 2018 113102

    2018

  35. [35]

    Kapustin A and Skorik S 1996 Journal of Physics A: Mathematical and General 29 1629–1638

    1996

  36. [36]

    Jiang Y, Liu Y C, Miao Y and Tan Z X 2025 ( Preprint arXiv:2512.01551)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  37. [37]

    Murgan R, Nepomechie R I and Shi C 2006 Boundary energy of th e open xxz chain from new exact solutions Annales Henri Poincar´ evol 7 (Springer) pp 1429–1448

    2006

  38. [38]

    Chernyak D, Gainutdinov A M, Jacobsen J L and Saleur H 2023 SIGMA 19 046

    2023

  39. [39]

    Lieb E, Schultz T and Mattis D 1961 Annals of Physics 16 407–466

    1961

  40. [40]

    Biegel D, Karbach M, M¨ uller G and Wiele K 2004 Physical Review B 69 174404

    2004

  41. [41]

    Zhang X, Kl¨ umper A and Popkov V 2022 Phys. Rev. B 106 075406 22

    2022