arxiv: 2605.05589 · v1 · submitted 2026-05-07 · ❄️ cond-mat.str-el · math-ph· math.MP

Recognition: unknown

Galois Solvability of Finite-Size Bethe Solutions in the Heisenberg Chain

Oliver R. Bellwood , William J. Munro

Authors on Pith no claims yet

Pith reviewed 2026-05-08 06:08 UTC · model grok-4.3

classification ❄️ cond-mat.str-el math-phmath.MP

keywords Heisenberg antiferromagnetic chainBethe ansatzGalois solvabilityfinite-size solutionsalgebraic complexityintegrable quantum modelsspin chains

0 comments

The pith

Bethe roots for the Heisenberg chain become Galois unsolvable at eight sites and above.

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

The paper examines the algebraic structure of exact ground states for the spin-1/2 Heisenberg antiferromagnetic chain obtained via the coordinate Bethe ansatz. Symbolic computations for chains up to ten sites show that the Bethe roots lose solvability by radicals starting at eight sites, while the wavefunction coefficients and energies follow at ten sites. This reveals that algebraic complexity increases rapidly with system size, preventing explicit closed-form expressions despite the model's formal integrability. A reader would care because it measures the practical boundary between exact solvability and usable analytic formulas in quantum many-body systems.

Core claim

Symbolic analysis demonstrates that the Bethe roots for chains of eight or more sites possess minimal polynomials with unsolvable Galois groups, so they cannot be expressed in radicals. The ground state wavefunction coefficients and energy become similarly unsolvable at ten or more sites. This establishes a lack of explicit analytic tractability arising from algebraic complexity in a canonical integrable model.

What carries the argument

The Galois group of the minimal polynomial satisfied by the Bethe roots and by the algebraic numbers appearing in the wavefunction coefficients, which decides whether the quantities admit expression by radicals.

If this is right

Explicit radical expressions for Bethe roots cease to exist for all chains of eight sites or larger.
Ground-state energies and wavefunctions lack closed-form algebraic expressions beyond nine sites.
Integrability alone does not guarantee that finite-size solutions remain symbolically tractable.
Numerical solution of the Bethe equations becomes the only practical route for moderate and large chain lengths.

Where Pith is reading between the lines

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

The same rapid growth of algebraic complexity may occur in other Bethe-ansatz integrable models when finite-size formulas are sought.
Methods that bypass explicit root extraction could restore tractability for larger systems.
The observed thresholds might shift under periodic versus open boundaries or for higher-spin variants, providing a direct test.

Load-bearing premise

The symbolic computations performed for sizes up to ten sites correctly detect the onset of Galois unsolvability and that this pattern holds for larger sizes in general.

What would settle it

An explicit computation of the Galois group for the Bethe roots of a nine-site chain that turns out to be solvable would contradict the reported onset of unsolvability.

Figures

Figures reproduced from arXiv: 2605.05589 by Oliver R. Bellwood, William J. Munro.

**Figure 1.** Figure 1: FIG. 1. Exact symbolic ground states of the two-, four-, six-, view at source ↗

**Figure 2.** Figure 2: FIG. 2. Algebraic complexity of the finite-size Bethe ansatz view at source ↗

read the original abstract

The spin-1/2 Heisenberg antiferromagnetic chain is the canonical example of an integrable quantum many-body model. Despite its exact solvability, explicit finite-size solutions are typically only accessible via numerical evaluation of the Bethe ansatz equations. Here, we analyse the algebraic structure of the exact, symbolic ground states for chains up to ten sites using the coordinate Bethe ansatz. We show that both the ground state wavefunction and the Bethe-roots rapidly develop algebraic complexity with respect to system size, but at different rates. The Bethe-roots appear to become Galois unsolvable for chains of eight or more sites, whereas the ground state wavefunction coefficients and energy appear to become unsolvable for ten or more sites. This demonstrates a lack of explicit analytic tractability in a quantum integrable model due to algebraic complexity.

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.

Desk Editor's Note private letter to a colleague

Symbolic computations show Bethe roots for the finite Heisenberg chain first acquire unsolvable Galois groups at N=8 and wavefunction coefficients at N=10, but the thresholds rest on unverified computer algebra.

read the letter

The main thing to know is that the paper turns the coordinate Bethe equations for the spin-1/2 Heisenberg antiferromagnet into polynomials over the rationals, extracts minimal polynomials for the roots and for the wavefunction coefficients, and then runs Galois-group calculations to test solvability by radicals. It reports that the roots become unsolvable starting at eight sites while the coefficients and energy do so at ten sites. This is a direct, concrete application of computational algebra to the exact finite-size solutions that people usually just solve numerically. It gives specific numbers for when explicit closed forms stop being feasible, which is a useful observation even if it is limited to small N. The work is honest about its scope and sticks to what the computations actually show. The soft spot is exactly what the stress test flags: everything depends on symbolic steps up to N=10 with no machine-checked certificates, no supplied code or raw polynomials, and no independent re-derivation. Galois-group algorithms for degree five and higher are reliable only when the input polynomial is correct, and a slip in clearing denominators or isolating the minimal polynomial would move the reported onsets. The paper uses “appear to,” which is the right cautious language, but it leaves the thresholds as observations rather than proven statements for all larger chains. There is no attempt to show the pattern persists or to find exceptions. This is for people who work on integrable models and care about the algebraic limits of the Bethe ansatz at finite size. A reader who wants to know why closed-form expressions dry up quickly will find it relevant; someone looking for new physics or large-N results will not. It is coherent on its own terms and shows clear engagement with the mathematics, so it deserves a serious referee who can check the computational pipeline and ask for the explicit polynomials or reproduction details.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes the algebraic structure of exact symbolic ground states for the spin-1/2 Heisenberg antiferromagnetic chain up to N=10 using the coordinate Bethe ansatz. It reports that the Bethe roots appear to acquire non-solvable Galois groups for chains of eight or more sites, while the ground-state wavefunction coefficients and energy appear to do so for ten or more sites, thereby demonstrating a lack of explicit analytic tractability arising from algebraic complexity.

Significance. If the reported thresholds hold under independent verification, the work supplies concrete evidence that algebraic complexity (measured via Galois solvability) sets in at small finite sizes even in an integrable model, quantifying the gap between formal Bethe-ansatz solvability and explicit closed-form expressions. This has implications for the practical limits of analytic solutions in quantum spin chains and similar integrable systems.

major comments (2)

[Results (Bethe-root analysis for N=8)] The central claim that Bethe roots become Galois-unsolvable at N=8 rests on symbolic conversion of the Bethe equations to polynomials over Q, extraction of minimal polynomials, and Galois-group computation. The manuscript provides no explicit minimal polynomials, no description of the clearing-denominators procedure, and no certificate or algorithm trace for the Galois-group determination (which for degree >4 is non-trivial). Without these, an error in root isolation or group computation could shift the reported onset, undermining the threshold statement.
[Results (wavefunction and energy analysis for N=10)] The parallel claim that wavefunction coefficients and energy become unsolvable at N=10 likewise depends on the same symbolic pipeline applied to the eigenvector components and the energy expression. No explicit polynomials, Galois-group output, or verification steps are supplied for these quantities either, leaving the N=10 threshold unsupported by reproducible evidence.

minor comments (2)

[Abstract and Introduction] The repeated use of 'appear to' in the abstract and conclusions is appropriately cautious given the finite-N scope, but the title and introductory framing could more explicitly qualify the results as computational observations up to N=10 rather than a general statement.
[Methods] Notation for the Bethe roots and the mapping from the coordinate Bethe equations to the minimal polynomials should be introduced with a short explicit example for a small N (e.g., N=4 or N=6) to aid readability before the larger-N claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and for highlighting the need for greater reproducibility in the algebraic computations. We agree that the absence of explicit polynomials and verification details weakens the presentation of the thresholds. We will revise the manuscript to include these elements as described below.

read point-by-point responses

Referee: [Results (Bethe-root analysis for N=8)] The central claim that Bethe roots become Galois-unsolvable at N=8 rests on symbolic conversion of the Bethe equations to polynomials over Q, extraction of minimal polynomials, and Galois-group computation. The manuscript provides no explicit minimal polynomials, no description of the clearing-denominators procedure, and no certificate or algorithm trace for the Galois-group determination (which for degree >4 is non-trivial). Without these, an error in root isolation or group computation could shift the reported onset, undermining the threshold statement.

Authors: We accept this criticism. The original manuscript omitted the explicit minimal polynomials for the N=8 Bethe roots, the clearing-denominators steps, and the Galois-group certificates. In the revised version we will add the minimal polynomials (obtained via resultant elimination over Q), a concise description of the denominator-clearing procedure, and the Galois-group output traces (computed with SageMath/Magma) that confirm the non-solvable groups. These additions will make the N=8 threshold independently verifiable. revision: yes
Referee: [Results (wavefunction and energy analysis for N=10)] The parallel claim that wavefunction coefficients and energy become unsolvable at N=10 likewise depends on the same symbolic pipeline applied to the eigenvector components and the energy expression. No explicit polynomials, Galois-group output, or verification steps are supplied for these quantities either, leaving the N=10 threshold unsupported by reproducible evidence.

Authors: We likewise accept the point. The revised manuscript will incorporate the minimal polynomials for the N=10 ground-state coefficients and energy, together with the corresponding Galois-group computations and verification steps. This will place the N=10 threshold on the same reproducible footing as the Bethe-root results. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation rests on independent symbolic computation

full rationale

The paper's central results follow from applying the coordinate Bethe ansatz to the Heisenberg chain, converting the resulting equations into polynomials over Q, extracting minimal polynomials for the Bethe roots and wavefunction coefficients, and computing their Galois groups via standard computer-algebra algorithms. These steps are direct, externally verifiable operations that do not reduce to any fitted parameter, self-referential definition, or load-bearing self-citation; the reported thresholds (N=8 for roots, N=10 for coefficients) are outputs of the computation rather than inputs. No ansatz is smuggled in, no known result is merely renamed, and the algebraic complexity claim is a straightforward consequence of the explicit calculations performed.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper rests on the standard Bethe ansatz framework and algebraic number theory without introducing new free parameters or entities.

axioms (2)

domain assumption The coordinate Bethe ansatz yields the exact eigenstates of the finite Heisenberg antiferromagnetic chain.
This is the foundational assumption for using Bethe ansatz to find symbolic ground states.
standard math Galois theory can be applied to determine if the minimal polynomials for Bethe roots and wavefunction coefficients are solvable by radicals.
Standard application of Galois group computations to algebraic numbers.

pith-pipeline@v0.9.0 · 5443 in / 1392 out tokens · 78294 ms · 2026-05-08T06:08:26.107407+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

35 extracted references

[1]

Bethe, Zur theorie der metalle: I

H. Bethe, Zur theorie der metalle: I. Eigenwerte und eigenfunktionen der linearen atomkette., Z. Phys.71, 205 (1931)

1931
[2]

R. J. Baxter, Exactly Solved Models in Statistical Me- chanics, inIntegrable Systems in Statistical Mechanics (World Scientific, 1985) pp. 5–63

1985
[3]

Takahashi,Thermodynamics of One-Dimensional Solvable Models(Cambridge University Press, 1999)

M. Takahashi,Thermodynamics of One-Dimensional Solvable Models(Cambridge University Press, 1999)

1999
[4]

Caux and J

J.-S. Caux and J. Mossel, Remarks on the notion of quan- tum integrability, Journal of Statistical Mechanics: The- ory and Experiment2011, P02023 (2011)

2011
[5]

Hulth´ en,¨Uber das austauschproblem eines kristalles, Ark

L. Hulth´ en,¨Uber das austauschproblem eines kristalles, Ark. Mat. Astron. Fys.26A, 1 (1938)

1938
[6]

Staudacher, Review of AdS/CFT Integrability, Chap- ter III.1: Bethe Ans¨ atze and the R-Matrix Formalism, Letters in Mathematical Physics99, 191 (2012)

M. Staudacher, Review of AdS/CFT Integrability, Chap- ter III.1: Bethe Ans¨ atze and the R-Matrix Formalism, Letters in Mathematical Physics99, 191 (2012)

2012
[7]

Larson, Integrability versus quantum thermalization, Journal of Physics B: Atomic, Molecular and Optical Physics46, 224016 (2013)

J. Larson, Integrability versus quantum thermalization, Journal of Physics B: Atomic, Molecular and Optical Physics46, 224016 (2013)

2013
[8]

M. T. Batchelor and H.-Q. Zhou, Integrability versus ex- act solvability in the quantum Rabi and Dicke models, Phys. Rev. A91, 053808 (2015)

2015
[9]

Links, Bethe Ansatz solution for a model of global- range interacting bosons on the square lattice, Journal of Statistical Mechanics: Theory and Experiment2026, 023103 (2026)

J. Links, Bethe Ansatz solution for a model of global- range interacting bosons on the square lattice, Journal of Statistical Mechanics: Theory and Experiment2026, 023103 (2026)

2026
[10]

C. N. Yang and C. P. Yang, One-Dimensional Chain of Anisotropic Spin-Spin Interactions. I. Proof of Bethe’s Hypothesis for Ground State in a Finite System, Phys. Rev.150, 321 (1966)

1966
[11]

Faddeev, Algebraic aspects of the Bethe ansatz, Inter- national Journal of Modern Physics A10, 1845 (1995)

L. Faddeev, Algebraic aspects of the Bethe ansatz, Inter- national Journal of Modern Physics A10, 1845 (1995)

1995
[12]

O. R. Bellwood and W. J. Munro, See Supplemental Ma- terial at [URL-to-be-inserted], which includes Refs. [30- 35]. (2026)

2026
[13]

Lieb and D

E. Lieb and D. Mattis, Ordering Energy Levels of Inter- acting Spin Systems, Journal of Mathematical Physics3, 749 (1962)

1962
[14]

Karbach and G

M. Karbach and G. M¨ uller, Introduction to the Bethe Ansatz I, Computer in Physics11, 36 (1997)

1997
[15]

Karbach, K

M. Karbach, K. Hu, and G. M¨ uller, Introduction to the Bethe Ansatz II, Computer in Physics12, 565 (1998)

1998
[16]

Rumer, E

G. Rumer, E. Teller, and H. Weyl, Nachrichten der Gesellschaft der Wissenschaften zu G¨ ottingen, Math.- Phys. Klasse , 409 (1932)

1932
[17]

H. N. V. Temperley and E. H. Lieb, Relations be- tween the ‘percolation’ and ‘colouring’ problem and other graph-theoretical problems associated with regular pla- nar lattices: some exact results for the ‘percolation’ prob- lem., Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences322, 251 (1971)

1971
[18]

Saito, A Proof of the Completeness of the Non Crossed Diagrams in Spin 1/2 Heisenberg Model, Journal of the Physical Society of Japan59, 482 (1990)

R. Saito, A Proof of the Completeness of the Non Crossed Diagrams in Spin 1/2 Heisenberg Model, Journal of the Physical Society of Japan59, 482 (1990)

1990
[19]

Beach and A

K. Beach and A. W. Sandvik, Some formal results for the valence bond basis, Nuclear Physics B750, 142 (2006)

2006
[20]

Liang, B

S. Liang, B. Doucot, and P. W. Anderson, Some New Variational Resonating-Valence-Bond-Type Wave Func- tions for the Spin-½Antiferromagnetic Heisenberg Model on a Square Lattice, Phys. Rev. Lett.61, 365 (1988). 5

1988
[21]

Oguchi and H

T. Oguchi and H. Kitatani, Theory of the Resonating Valence Bond in Quantum Spin System, Journal of the Physical Society of Japan58, 1403 (1989)

1989
[22]

K. Seki, T. Shirakawa, and S. Yunoki, Symmetry-adapted variational quantum eigensolver, Phys. Rev. A101, 052340 (2020)

2020
[23]

Poilblanc, Exact states and variational states of a pe- riodic two-dimensional cluster of the heisenberg model, Phys

D. Poilblanc, Exact states and variational states of a pe- riodic two-dimensional cluster of the heisenberg model, Phys. Rev. B39, 140 (1989)

1989
[24]

Oguchi, H

T. Oguchi, H. Nishimori, and Y. Taguchi, Ground State of Antiferromagnetic Quantum Spin Systems on the Tri- angular Lattice, Journal of the Physical Society of Japan 55, 323 (1986)

1986
[25]

Harary and E

F. Harary and E. M. Palmer,Graphical enumeration (Academic Press, 1973)

1973
[26]

D. W. Walkup, The number of plane trees, Mathematika 19, 200 (1972)

1972
[27]

R. P. Stauduhar, The Determination of Galois Groups, Mathematics of Computation27, 981 (1973)

1973
[28]

Fieker and J

C. Fieker and J. Kl¨ uners, Computation of Galois groups of rational polynomials, LMS Journal of Computation and Mathematics17, 141–158 (2014)

2014
[29]

Available at: https://magma.maths.usyd.edu.au/magma/
[30]

Kouzoudis, Exact analytical partition function and energy levels for a Heisenberg ring of N=6 spin 1/2 sites, Journal of Magnetism and Magnetic Materials189, 366 (1998)

D. Kouzoudis, Exact analytical partition function and energy levels for a Heisenberg ring of N=6 spin 1/2 sites, Journal of Magnetism and Magnetic Materials189, 366 (1998)

1998
[31]

B¨ arwinkel, H.-J

K. B¨ arwinkel, H.-J. Schmidt, and J. Schnack, Structure and relevant dimension of the Heisenberg model and ap- plications to spin rings, Journal of Magnetism and Mag- netic Materials212, 240 (2000)

2000
[32]

B. Buchberger, Bruno Buchberger’s PhD thesis 1965: An algorithm for finding the basis elements of the residue class ring of a zero dimensional polynomial ideal, Journal of Symbolic Computation41, 475 (2006)

1965
[33]

OEIS Foundation Inc., The On-Line Encyclopedia of In- teger Sequences, A006840 (2026)

2026
[34]

Fredricksen and I

H. Fredricksen and I. J. Kessler, An algorithm for gen- erating necklaces of beads in two colors, Discrete Math- ematics61, 181 (1986)

1986
[35]

OEIS Foundation Inc., The On-Line Encyclopedia of In- teger Sequences, A002995 (2026)

2026