Negative Spin XXX Chain Model

• The negative-spin XXX chain is a noncompact quantum model with s=-1 that captures reggeized gluon dynamics in high-energy QCD.
• It is exactly solved using the Bethe ansatz, yielding a Fermi sea of real rapidities and fermionic topological solitons (lipatons) unique to the model.
• The model maps to a quantum lattice NLS chain and exhibits Luttinger liquid and conformal field theory behavior, highlighting its critical thermodynamics.

The negative-spin XXX chain is an integrable quantum spin chain where each site carries a spin representation $s=-1$ of $SU(2)$. It arises as an effective lattice model in high-energy quantum chromodynamics (QCD), specifically in the large $N_c$ limit, where reggeized gluon dynamics reduce to nearest-neighbor interactions of non-compact $SU(2)$ spins. This model is mathematically equivalent to a quantum lattice version of the repulsive nonlinear Schrödinger (NLS) equation, describing a chain of interacting bosonic harmonic oscillators. Uniquely, its elementary excitations are fermionic topological solitons (“lipatons”), and its vacuum and thermodynamics are distinct from conventional positive-spin ($s>0$) XXX chains. The model admits an exact solution by the Bethe ansatz, a well-defined thermodynamic Bethe ansatz (TBA), a conformal field theory (CFT) low-energy limit, and displays Luttinger liquid behavior with parameters and scaling dissimilar to its positive-spin and Lieb–Liniger counterparts (Hao et al., 2019, Zhong et al., 3 Feb 2026).

1. Hamiltonian Formulation, R-Matrix, and Mapping to Lattice NLS

The local Hamiltonian for the negative-spin XXX chain derives from the $SU(2)$ R-matrix for arbitrary spin $s$: $$R_{jk}^{(s,s)}(\lambda) = f(s, \lambda) \frac{\Gamma(i\lambda - 2s)\, \Gamma(i\lambda + 2s + 1)}{\Gamma(i\lambda - J_{jk})\, \Gamma(i\lambda + J_{jk} + 1)},$$ where $$J_{jk}(J_{jk} + 1) = 2\vec{S}_j \cdot \vec{S}_k + 2s(s+1)$$ and $f(s, \lambda)$ is a normalization factor. The local Hamiltonian density is given by

$$H_{jk} = \left.\frac{1}{i}\frac{d}{d\lambda}\ln R_{jk}^{(s,s)}(\lambda)\right|_{\lambda=0}.$$

For the holomorphic noncompact representation relevant to QCD, $s=-1$ so $s(s+1)=0$ and $J_{jk}(J_{jk}+1) = 2\vec{S}_j \cdot \vec{S}_k$, leading to

$$H_{jk}^{(s=-1)} = \frac{1}{i} \frac{d}{d\lambda}\ln R_{jk}^{(-1,-1)}(\lambda)\Big|_{\lambda=0}.$$

The total Hamiltonian is $H^{(s=-1)} = \sum_{j=1}^L H_{j,j+1}^{(s=-1)}$.

A direct mapping exists to a quantum lattice NLS chain with bosonic creation/annihilation operators $\Psi_j, \Psi_j^\dagger$ satisfying $[\Psi_j, \Psi_k^\dagger] = \delta_{jk}$. The Hamiltonian is

$$H_{\rm NLS} = -\sum_{j=1}^L (\Psi_j^\dagger\Psi_{j+1} + \Psi_{j+1}^\dagger\Psi_j) + c\sum_{j=1}^L \Psi_j^\dagger\Psi_j^\dagger\Psi_j\Psi_j - \mu\sum_{j=1}^L \Psi_j^\dagger\Psi_j,$$

with $c>0$ and $\mu$ the chemical potential. In the continuum limit, this yields the Lieb–Liniger model. The equivalence is made precise by identifying $s = -2/(\kappa \Delta)$, where $\kappa$ and $\Delta$ are coupling and lattice-spacing parameters, with $s=-1$ corresponding to $\kappa=1$, $\Delta=2$ (Hao et al., 2019, Zhong et al., 3 Feb 2026).

2. Bethe Ansatz: Spectrum, Equations, and Thermodynamic Limit

The spectrum is determined using the algebraic Bethe ansatz. For a chain of length $L$ and $N$ reversed spins (“reggeized gluons”), Bethe rapidities $\{\lambda_k\}$ satisfy periodic Bethe equations: $$\left(\frac{\lambda_k - i}{\lambda_k + i}\right)^L = \prod_{j \neq k} \frac{\lambda_k - \lambda_j + i}{\lambda_k - \lambda_j - i}, \qquad k = 1, \ldots, N.$$ The bare energy and momentum of each rapidity are

$$\varepsilon_0(\lambda) = -\frac{2}{\lambda^2 + 1}, \qquad p_0(\lambda) = i\ln\left(\frac{i+\lambda}{i-\lambda}\right).$$

All Bethe roots are real for $s=-1$, so no bound-state (string) solutions occur—this property contrasts with positive-$s$ chains (Hao et al., 2019, Zhong et al., 3 Feb 2026).

In the thermodynamic limit ($L,N \to \infty$ with $N/L = n$ fixed), root and hole densities $\rho(\lambda)$ and $\rho_h(\lambda)$ satisfy

$$\rho(\lambda) + \rho_h(\lambda) = \frac{1}{\pi(1+\lambda^2)} + \frac{1}{2\pi}\int_{-\infty}^\infty \frac{2}{1+(\lambda-\mu)^2}\rho(\mu)\,d\mu,$$

with ground-state support $[-\lambda_F, \lambda_F]$ determined by the filling $n$.

3. Thermodynamic Bethe Ansatz and Quantum Criticality

At finite temperature $T$ and chemical potential $h$, the thermodynamic Bethe ansatz (TBA) yields the dressed energy $\varepsilon(\lambda)$ as the solution to

$$\varepsilon(\lambda) = -\frac{2}{\lambda^2+1} - h - \frac{T}{2\pi} \int_{-\infty}^{+\infty} \frac{2}{1 + (\lambda - \mu)^2} \ln\bigl[1 + e^{-\varepsilon(\mu)/T}\bigr] d\mu.$$

The pressure per unit length is

$$p(T, h) = \frac{T}{2\pi} \int_{-\infty}^{+\infty} \frac{2}{1+\lambda^2} \ln[1 + e^{-\varepsilon(\lambda)/T}] d\lambda,$$

and all thermodynamic observables follow by differentiation: $$n = \partial_h p,\quad s = \partial_T p,\quad \kappa = \partial_h n,\quad c_V = T\partial_T s.$$ At $T=0$, a quantum critical point is located at $h_c = -2$:

• For $h < h_c$: vacuum ($n=0$),
• For $h > h_c$: gapless Luttinger-liquid phase.

Near $h\approx h_c$, scaling laws take the form $n(h,T) \sim \sqrt{T} \mathcal{F}_n[(h-h_c)/T]$, exhibiting critical exponents $z=2$, $\nu=1/2$ (Zhong et al., 3 Feb 2026).

4. Structure of Excitations and Soliton-Fermion Duality

In the ground state, the many-body spectrum consists of a Fermi sea of real rapidities. Removing (adding) a rapidity corresponds to creating a hole (particle) excitation. Single particle–hole excitations are described by a shift function $F(\lambda|\lambda_p,\lambda_h)$ which solves

$$F(\lambda|\lambda_p,\lambda_h) = \frac{1}{2\pi}\int_{-\lambda_F}^{\lambda_F} \frac{2}{1+(\lambda-\mu)^2} F(\mu|\lambda_p,\lambda_h) d\mu + \frac{\theta(\lambda-\lambda_p) - \theta(\lambda-\lambda_h)}{2\pi}.$$

Energy and momentum follows from

$$\Delta E = \epsilon(\lambda_p) - \epsilon(\lambda_h) - \int_{-\lambda_F}^{\lambda_F} \epsilon'(\lambda) F(\lambda|\lambda_p,\lambda_h) d\lambda,$$

$$\Delta P = \theta(\lambda_p) - \theta(\lambda_h) - \int_{-\lambda_F}^{\lambda_F} \theta'(\lambda) F(\lambda|\lambda_p,\lambda_h) d\lambda,$$

with $\epsilon(\lambda) = -2/(1 + \lambda^2)$, $\theta(\lambda) = 2\arctan\lambda$.

The elementary excitations—"lipatons"—are fermions, forming Slater-determinant–like states and showing $Z_2$ soliton statistics in the bosonic oscillator description. The absence of bound states is enforced by the real Bethe roots, in contrast to the positive-spin chain (Hao et al., 2019, Zhong et al., 3 Feb 2026).

5. Continuum and Lattice NLS Correspondence

There is an exact mapping at the level of the Bethe equations and thermodynamics between the $s=-1$ spin chain and the quantum lattice NLS model with repulsive interactions. At each site, one identifies bosonic operators and a realization of $su(2)$: $$[\Psi_j, \Psi_k^\dagger] = \delta_{jk}, \quad \varrho_j = \sqrt{1+\frac{\kappa\Delta}{4}\Psi_j^\dagger\Psi_j}, \quad s = -\frac{2}{\kappa\Delta}.$$

$$S_j^x = \frac{i}{\sqrt{\kappa\Delta}}(\Psi_j^\dagger\varrho_j + \varrho_j\Psi_j),\quad S_j^y = \frac{1}{\sqrt{\kappa\Delta}}(\varrho_j\Psi_j - \Psi_j^\dagger\varrho_j), \quad S_j^z = -\frac{2}{\kappa\Delta}(1 + \frac{\kappa\Delta}{2}\Psi_j^\dagger\Psi_j).$$

At $\kappa=1, \Delta=2$, this yields $s=-1$, matching the Bethe equations and excitation spectrum between the two models. In the continuum ($\Delta\to0$), the chain reduces to the Lieb–Liniger field theory of repulsive bosons, but the thermodynamics are not continuously connected to the positive-spin XXX case (Hao et al., 2019, Zhong et al., 3 Feb 2026).

6. Low-Energy Conformal Field Theory and Luttinger Liquid Regime

In the $T\to 0$, $h > h_c$ regime, the negative-spin XXX chain exhibits a CFT with central charge $c=1$. The dressed energy $\varepsilon(\lambda)\sim v_s(|\lambda|-\lambda_F)$ near the Fermi points, yielding a linear spectrum and identifying a sound velocity $v_s$. The low-energy effective Hamiltonian is a Luttinger liquid: $$H_{\rm LL} = \frac{v_s}{2\pi} \int dx \left[ K(\partial_x\Theta)^2 + \frac{1}{K}(\partial_x\Phi)^2 \right],$$ with the Luttinger parameter $K=\pi v_s\kappa$, where $\kappa$ is the zero-temperature compressibility. Scaling dimensions are

$$\Delta_{n,m} = \frac{n^2}{4K} + \frac{m^2 K}{4}, \quad n, m\in\mathbb{Z},$$

and entanglement entropy of an interval $y$ at $T=0$ follows $S(y)\simeq \frac{c}{3} \ln y$ (Hao et al., 2019, Zhong et al., 3 Feb 2026).

7. Distinctions from Positive-Spin XXX and Related Models

The negative-spin chain differs fundamentally from the positive-spin XXX model:

Feature	$s=-1$ XXX chain	$s>0$ XXX chain
Bethe roots	All real; no strings	Complex (string) solutions exist
Excitations	Lipatons (fermionic solitons)	Magnons, spinons (bosonic)
Hilbert space	Non-compact, infinite-dim.	Finite-dim., compact
Thermodynamic regime	No analytic continuation from $s>0$	Adiabatically connected
CFT exponents	$c=1$, distinct $\mathcal{Z}$	$c=1$, different $\mathcal{Z}$

Lipatons, absent bound states, and the non-compact nature of the representation reflect the model's unique physical content and its emergence as an effective QCD theory. The low-$T$ Luttinger-liquid regime and quantum phase transition at $h_c=-2$ further distinguish the negative-spin chain from both conventional XXX and Lieb–Liniger models (Hao et al., 2019, Zhong et al., 3 Feb 2026).