Systematic analytic finite-size corrections for the periodic spin-1/2 isotropic antiferromagnetic Heisenberg chain: convergent Bethe-ansatz rapidity asymptotics, higher-order RG energy series, and excited states

For the spin-1/2 isotropic antiferromagnetic Heisenberg chain, the ground-state energy obeys $E_L = L\varepsilon_\infty - \frac{\pi^2}{6L}\big(1 + \frac{C}{\ln^3 L} + \ldots\big)$ with $C = 3/8$, set by the marginally irrelevant current-current operator of the $SU(2)_1$ Wess-Zumino-Witten (WZW) continuum limit. Nickel's Lagrange-interpolation/Newton-Raphson solutions of the Bethe ansatz equations to $L = 2^{21}$ confirmed this amplitude numerically and resolved a decades-old discrepancy between analytic ($C\approx 0.3433$) and numerical ($C\approx 0.36516$) determinations, but left the analytic side open: the first-order rapidity corrections $\nu_n^{(1)}$ are only semi-quantitative, the higher-order RG structure was treated as unknown fitting freedom, and the $O(\ln L)$-variable method is unproven for excited states. Later literature supplies a key constraint the source paper did not use: in Zamolodchikov's scheme the renormalization-group (RG) flow near the $SU(2)_1$ point is exactly known and reduces on the $SU(2)$ line to $\beta(g) = -g^2/(1-g/2)$. The problem is to convert this into a closed, systematically improvable analytic theory of the finite-size corrections - an energy series beyond the known $\frac{3}{8}g^3$ term, convergent rapidity-level corrections, and an excited-state extension - validated to high precision against direct Bethe-ansatz solutions.

Background

Consider the periodic spin-1/2 isotropic (XXX) antiferromagnetic Heisenberg chain $H = \frac12\sum_{n=1}^{L}\boldsymbol{\sigma}_n\cdot\boldsymbol{\sigma}_{n+1}$. The ground-state energy density is Hulthén's $\varepsilon_\infty = \frac12 - 2\ln 2 \approx -0.886294$, and the ground state is the $M = L/2$ real-rapidity solution of the Bethe ansatz equations (BAE) $\big(\frac{\lambda_j + i/2}{\lambda_j - i/2}\big)^L = \prod_{k\neq j}\frac{\lambda_j - \lambda_k + i}{\lambda_j - \lambda_k - i}$. Continuum theory (Affleck-Gepner-Schulz-Ziman; Cardy) identifies the scaling limit as the $SU(2)_1$ Wess-Zumino-Witten conformal field theory perturbed by a marginally irrelevant current-current operator (a coupling that flows to zero only logarithmically slowly with distance, which is what generates the inverse powers of $\ln L$ below), giving the finite-size form $E_L = L\varepsilon_\infty - \frac{\pi^2}{6L}\big(1 + \frac{C}{\ln^3 L} + \ldots\big)$ with universal amplitude $C = 3/8$. Historically this amplitude was contentious: the analytic-continuation calculation of Woynarovich-Eckle gave $C\approx 0.3433$, while Nomura's careful numerics gave $C\approx 0.36516$, and direct high-precision BAE numerics at moderate $L$ could not discriminate cleanly against $3/8 = 0.375$. The exact low-energy theory is in fact stronger than the source paper used. Lukyanov's exact effective Hamiltonian for the XXZ chain, building on Zamolodchikov's exact renormalization-group (RG) scheme, fixes the flow near the $SU(2)_1$ point exactly: $r\,\mathrm{d}g_\parallel/\mathrm{d}r = -2g_\perp^2/(2-g_\parallel)$, $r\,\mathrm{d}g_\perp/\mathrm{d}r = -2g_\parallel g_\perp/(2-g_\parallel)$, which on the $SU(2)$-invariant line $g_\parallel = g_\perp = g$ reduces to the closed single-variable flow $\beta(g) = -\frac{g^2}{1-g/2} = -g^2 - \frac12 g^3 - \frac14 g^4 - \ldots$, with the running coupling fixed by $g^{-1} + \frac12\ln g = \ln L + \text{const}$ with exactly known constants. The finite-size ground-state energy then takes the RG-improved form $E_L = L\varepsilon_\infty - \frac{\pi^2}{6L}\big(1 + \frac{3}{8}g^3 + O(g^4)\big)$, and the excited $SU(2)_1$ towers carry their own exact $g$-series (in the spinon quantum numbers $s,\theta$ labeling the towers of the chain's fractionalized spin-$1/2$ excitations: coefficients $\{1 - \frac{g}{2} + \frac{g^2}{4} - \frac{7}{32}g^3\}$, $\{2g^2 - g^3\}$, $\{1 + \frac{g}{2} + \frac{3g^3}{32}\}$). Lukyanov-Terras clarified the scheme structure: odd coefficients of the running-coupling relation are scheme-invariant, even ones are convention. On the numerical side, Nickel solved the BAE for chains up to $L = 2^{21}$ to about 40 digits using Lagrange interpolation of the rapidities on a uniform grid ($M = \mathrm{round}(4\log_2 L + 8)$ basis rapidities) followed by Newton-Raphson iteration in the $O(\ln L)$ basis variables, confirmed the leading-log convergence and $C = 3/8$ with the Nomura data for $L\le 2^{14}$, and fitted the non-universal scale $L_0 = 0.48364(2)$ using the two-loop flow $\beta(g) = -g^2 - \frac12 g^3 - g^4\alpha(g)$ with unknown $\alpha(g)$. Nickel also derived an improved first-order correction $\nu_n^{(1)}$ to the Hulthén rapidities $\lambda_n^{(0)} = -\frac{2}{\pi}\ln\tan\frac{n\pi}{L}$ with half-integer indices $n = \frac12, \frac32, \dots, \frac{L-1}{2}$ (the correction $\nu_n^{(1)}$ is a gamma-function integral, eq. (33) of the source), which however is only semi-quantitatively accurate. Bortz-Karbach-Schneider-Eggert independently compared lattice and continuum perturbation theory, confirming the $-3/(8\ln^3 N)$ structure and providing systematic $1/N$ expansions, but not beyond-$O(g^3)$ energy coefficients nor any rapidity-level scheme.

Problem Statement

Goal: replace curve-fitting of finite-size data with a closed, systematically improvable analytic theory of the finite-size corrections of the periodic spin-1/2 isotropic antiferromagnetic Heisenberg chain $H = \frac12\sum_{n=1}^{L}\boldsymbol{\sigma}_n\cdot\boldsymbol{\sigma}_{n+1}$, anchored in the exactly known RG flow and validated to high precision against direct solutions of the Bethe ansatz equations (BAE). The energy scale is fixed by $E_L = L\varepsilon_\infty - \frac{\pi^2}{6L}\,c_{\mathrm{eff}}(L)$ with $\varepsilon_\infty = \frac12 - 2\ln 2$, and the exactly known RG input is the Zamolodchikov-scheme flow $\beta(g) = -g^2/(1-g/2)$ with the running coupling fixed by Lukyanov's exact relation, whose ground-state ($s = \theta = 0$) case is $g^{-1} + \frac12\ln g = \ln\big(2^{1/2}\pi^{-1/2}e^{1/4}L\big) + \gamma_E$ (excited sectors carry the additional $s,\theta$-dependent shift given in Lukyanov (1998)). Deliverables: (a) Energy series beyond the known order: compute the $O(g^4)$ and $O(g^5)$ coefficients $c_4, c_5$ in $c_{\mathrm{eff}} = 1 + \frac{3}{8}g^3 + c_4 g^4 + c_5 g^5 + \ldots$ (equivalently, the $1/\ln^4 L$ and $1/\ln^5 L$ terms including their $\ln\ln L$ content), stating which coefficients are scheme-invariant; and perform the matching between this scheme and the lattice Bethe-ansatz fit parametrization, including its non-universal constants and normalization conventions. (b) Convergent rapidity-level corrections: construct an analytic, iteratively improvable scheme for the ground-state rapidities $\lambda_n = \lambda_n^{(0)} + \nu_n$, where $\lambda_n^{(0)} = -\frac{2}{\pi}\ln\tan\frac{n\pi}{L}$ (half-integer indices $n = \frac12, \frac32, \dots, \frac{L-1}{2}$) is the Hulthén distribution and the known first-order $\nu_n^{(1)}$ is only semi-quantitative, with controlled successive terms (a stated asymptotic ordering, and ideally a convergence proof), such that the induced energies agree with directly computed high-precision BAE values (e.g. Nickel's data, exact to about 40 digits for $2^{10} \le L \le 2^{21}$) to the accuracy that the constructed asymptotic ordering permits; the scheme must resolve the analytic-continuation ambiguity of de Vega-Woynarovich-type finite-size formulas in a principled way rather than by ad hoc branch choices. (c) Excited states: extend the $O(\ln L)$-variable Lagrange-interpolation/Newton-Raphson method (or an equivalent reduced-variable method) to excited-state rapidity distributions that lack the ground state's smoothness, demonstrating stability and accuracy of the same class for at least the low-lying $SU(2)_1$ towers, and cross-validating the resulting excitation energies against the exact tower $g$-series ($s^2$, $|s|$, $\theta^2$ coefficients). (d) Numerical validation: with (a)-(c), extract the leading amplitude $C$ from independent fits at the best precision the constructed theory supports and check consistency with the exact continuum value $C = 3/8$; for context, Nickel's truncated two-loop fits could not determine $C$ better than about 1 part in $10^{5}$. A complete solution provides the closed-form or algorithmic results for (a)-(c) together with the numerical validation of (d); partial progress on any component is meaningful if the ordering/convergence structure is explicitly characterized.

References

1. B. Nickel, Scaling corrections to the ground state energy of the spin-1/2 isotropic anti-ferromagnetic Heisenberg chain, J. Phys. Commun. 1 (2017) 055021, doi:10.1088/2399-6528/aa9a15, https://doi.org/10.1088/2399-6528/aa9a15
2. B. Nickel, Observations on the Bethe ansatz solutions of the spin-1/2 isotropic anti-ferromagnetic Heisenberg chain, J. Phys. Commun. 3 (2019) 025007, doi:10.1088/2399-6528/ab030b, https://doi.org/10.1088/2399-6528/ab030b
3. F. Woynarovich, H.-P. Eckle, Finite-size corrections and numerical calculations for long spin 1/2 Heisenberg chains in the critical region, J. Phys. A: Math. Gen. 20 (1987) L97-L104, doi:10.1088/0305-4470/20/2/010, https://doi.org/10.1088/0305-4470/20/2/010
4. K. Nomura, Logarithmic corrections of the one-dimensional S=1/2 Heisenberg antiferromagnet, Phys. Rev. B 48 (1993) 16814-16817, doi:10.1103/PhysRevB.48.16814, https://doi.org/10.1103/PhysRevB.48.16814
5. I. Affleck, D. Gepner, H. J. Schulz, T. Ziman, Critical behaviour of spin-s Heisenberg antiferromagnetic chains: analytic and numerical results, J. Phys. A: Math. Gen. 22 (1989) 511-529, doi:10.1088/0305-4470/22/5/015, https://doi.org/10.1088/0305-4470/22/5/015
6. J. L. Cardy, Conformal invariance and universality in finite-size scaling, J. Phys. A: Math. Gen. 17 (1984) L385-L387, doi:10.1088/0305-4470/17/7/003, https://doi.org/10.1088/0305-4470/17/7/003
7. H. J. de Vega, F. Woynarovich, Method for calculating finite size corrections in Bethe ansatz systems: Heisenberg chain and six-vertex model, Nucl. Phys. B 251 (1985) 439-456, doi:10.1016/0550-3213(85)90271-8, https://doi.org/10.1016/0550-3213(85)90271-8
8. A. B. Zamolodchikov, Mass scale in the sine-Gordon model and its reductions, Int. J. Mod. Phys. A 10 (1995) 1125-1150, doi:10.1142/S0217751X9500053X, https://doi.org/10.1142/S0217751X9500053X
9. S. Lukyanov, Low energy effective Hamiltonian for the XXZ spin chain, Nucl. Phys. B 522 (1998) 533-549, arXiv:cond-mat/9712314, doi:10.1016/S0550-3213(98)00249-1, https://doi.org/10.1016/S0550-3213(98)00249-1
10. S. Lukyanov, V. Terras, Long-distance asymptotics of spin-spin correlation functions for the XXZ spin chain, Nucl. Phys. B 654 (2003) 323-356, arXiv:hep-th/0206093, doi:10.1016/S0550-3213(02)01141-0, https://doi.org/10.1016/S0550-3213(02)01141-0
11. M. Bortz, M. Karbach, I. Schneider, S. Eggert, Lattice versus continuum theory of the periodic Heisenberg chain, Phys. Rev. B 79 (2009) 245414, arXiv:0904.0001, doi:10.1103/PhysRevB.79.245414, https://doi.org/10.1103/PhysRevB.79.245414