Integrability and multisoliton solutions of the massive easy-axis O(3) sigma model, Lorentz-breaking anisotropy effects on the soliton sector, and nonperturbative soliton energies in the one-dimensional antiferromagnet

Haldane's 1983 nonlinear field theory of the large-spin one-dimensional easy-axis Heisenberg antiferromagnet derives, as the continuum limit, the O(3) nonlinear sigma model with an easy-axis mass term, constructs its single soliton (a domain wall between the two degenerate Néel states), and quantizes it semiclassically. The paper explicitly leaves open: (i) whether the classical equations of motion remain integrable for nonzero mass, allowing multisoliton generalizations of the single-soliton solution; and (iii) nonperturbatively renormalized soliton rest energies and precise critical thresholds beyond the leading semiclassical and parametric results; it also introduces (ii) a Lorentz-breaking anisotropy term of strength $\gamma_1=(\mu+\lambda)/2$ but sets it aside ("will not be considered further here"), leaving its consequences for the soliton sector undetermined. A fourth question raised by the same paper, the required nonstandard quantization for half-integer spin, is resolved by the topological $\theta=2\pi S$ term and the Lieb-Schultz-Mattis theorem and is not part of this problem. The three questions (i), (ii), and (iii) — the completion of the soliton theory of this model — remain open: no Lax pair, multisoliton construction, or non-integrability proof is known for the massive model; no dedicated quantitative treatment of the $\gamma_1$ term on the soliton sector is known; and the renormalized soliton energies are controlled only at the exactly solved $S=1/2$ and studied $S=1$ endpoints, with no controlled computation in between.

Background

A one-dimensional antiferromagnetic spin- $S$ chain with easy-axis anisotropy is described at the lattice level by $H=|J|\sum_n[\mathbf{S}_n\cdot\mathbf{S}_{n+1}+\lambda S_n^z S_{n+1}^z+\mu(S_n^z)^2]$ with $\lambda,\mu>0$ (the $XXZ$ chain with single-ion anisotropy). Its classical ground state is the Néel state, in which neighboring spins alternate between two opposite orientations; "easy-axis" means the anisotropy favors alignment along the $z$ axis, leaving two degenerate Néel states related by a global spin flip. In the large- $S$ (semiclassical) limit, coarse-graining over many lattice spacings $a$ yields a continuum field theory for the unit-length Néel order parameter $\Omega(x,t)\in S^2$ : the O(3) nonlinear sigma model (NLSM) — the universal field theory of slowly varying order-parameter dynamics on the sphere — augmented here by an easy-axis mass term, $\mathcal{L}=(1/2g)[c^{-1}|\partial_t\Omega|^2-c|\partial_x\Omega|^2-(\omega_0^2/c)\Omega_\perp^2]$ , where $g=2/\hbar S$ is the coupling, $c$ the spin-wave velocity, $\omega_0$ the anisotropy frequency, and $\Omega_\perp$ the components of $\Omega$ transverse to the easy axis. Because the two Néel states are degenerate, the model admits a soliton: a localized domain wall in which $\Omega$ rotates between the two easy-axis orientations, $\cos\theta=\pm\tanh[(s-s_0)/R]$ , where $\theta$ is the polar angle of $\Omega$ , $s$ the spatial coordinate in the soliton's comoving frame (related to $(x,t)$ by a Lorentz boost of velocity $v$ ), $s_0$ the soliton position, $R$ the width, and $\omega$ the soliton's internal azimuthal precession frequency, which fixes the width through $c^2/R^2=\omega_0^2-\omega^2$ (the static, non-precessing soliton has $\omega=0$ and $R=c/\omega_0$ ; after semiclassical quantization $\omega=m\omega_0/(m^2+S^2)^{1/2}$ ). The soliton's semiclassical dispersion is $E_m(P)=[(m^2+S^2)(\hbar\omega_0)^2+c^2P^2]^{1/2}$ , where the quantized soliton spin $S^z=m\hbar$ takes integer (half-integer) $m$ for integer (half-integer) $S$ by a time-reversal argument. A soliton is the field-theory avatar of the domain wall that flips the Néel order; in spin chains a deconfined domain-wall fragment carrying spin $1/2$ is called a spinon. Quantum tunneling between topologically distinct sectors is generated by instantons (finite-action Euclidean trajectories connecting the vacua), giving a tunneling gap $\varepsilon_0$ exponentially small in the semiclassical parameter. The magnon is the elementary collective excitation carrying $S^z=\pm1$ , obtained by semiclassical quantization of the spin waves, with dispersion $\varepsilon(P)=[(\hbar\omega_0)^2+c^2P^2]^{1/2}$ ; near the isotropic limit the nonlinear zero-point fluctuations ("instantons") renormalize the coupling $g$ to strong coupling and give the collective magnon mode a finite rest energy $\varepsilon_0=\xi^{-1}cg^{-1}\exp(-2\pi/g)$ , with $\xi\sim a$ the ultraviolet cutoff length scale. Semiclassically the soliton energy gap is much larger than the magnon gap, while near the isotropic limit the renormalized soliton rest energy becomes lower than the renormalized magnon rest energy; for integer $S$ the $m=0$ soliton gap vanishes at a finite critical anisotropy $(\hbar\omega_0)_c\sim\varepsilon_0$ , signaling a doublet-singlet transition of the Onsager/Ising-field-theory type to the singlet ground state of the nearly isotropic model exhibited by the strong-coupling rotator version of the model (for half-integer $S$ the $m=\pm1/2$ soliton gap vanishes only at the isotropic point $\omega_0=0$ , a direct transition of the Kosterlitz-Thouless/sine-Gordon type to the gapless easy-plane state). The sine-Gordon model is the integrable scalar field theory of a periodic potential whose kinks describe domain walls; when the azimuthal angle $\phi$ of $\Omega$ is held fixed, the easy-axis sigma model reduces statically to a sine-Gordon-type problem, and in the easy-plane limit the full chain maps to sine-Gordon itself. The Kosterlitz-Thouless (KT) transition is the topological vortex-unbinding transition that governs the easy-plane chain at the isotropic point. Two further anchors of the modern literature: the Lieb-Schultz-Mattis (LSM) theorem, originally for $S=1/2$ , and its generalizations prove that a half-integer-spin chain with one spin per unit cell cannot have a unique gapped ground state; and the topological theta term, a total-derivative term in the NLSM action with coefficient $\theta=2\pi S$ (Haldane's companion paper), encodes the distinction between integer and half-integer spin. At $\theta=\pi$ the O(3) NLSM is massless (flowing at long distances to the $SU(2)_1$ Wess-Zumino-Witten conformal field theory, the exactly solvable scale-invariant theory describing the critical half-integer-spin chain), while away from it the spectrum is massive. The $S=1/2$ $XXZ$ chain is exactly solvable by the Bethe ansatz, the algebraic method yielding exact spectra in the massive regime; its elementary excitations are kinks (spinons). With these definitions, the open questions addressed here concern the internal consistency and completion of the soliton picture: the integrability and multisoliton content of the classical massive model, the role of the Lorentz-breaking anisotropy $\gamma_1=(\mu+\lambda)/2$ that couples to the soliton sector through $H'=(\gamma_1 c/2)\int dx[gL^2+g^{-1}(\nabla\cos\theta)^2]$ (with $L$ the azimuthal spin density), and the true nonperturbative soliton energies beyond the leading semiclassical order.

Problem Statement

Complete the soliton theory of the one-dimensional easy-axis antiferromagnetic continuum model — the O(3) nonlinear sigma model with easy-axis mass term, with Lagrangian density

\mathcal{L}=\frac{1}{2g}\left[c^{-1}|\partial_t\Omega|^2-c|\partial_x\Omega|^2-\frac{\omega_0^2}{c}\Omega_\perp^2\right]

for the unit Néel vector $\Omega(x,t)\in S^2$ (Haldane, Phys. Rev. Lett. 50, 1153 (1983), Eqs. (8)-(10), DOI: 10.1103/PhysRevLett.50.1153; $g=2/\hbar S$ is the coupling for spin quantum number $S$ , $c$ the spin-wave velocity, $\omega_0\neq 0$ the anisotropy frequency, and $\Omega_\perp$ the components of $\Omega$ transverse to the easy axis) — by settling three linked questions the source paper left open. (a) Classical integrability and multisolitons: determine whether the equations of motion for $\omega_0\neq 0$ are integrable (the $\omega_0=0$ case is integrable by Pohlmeyer's reduction): either establish integrability and construct multisoliton solutions — exact multi-soliton scattering configurations and bound states generalizing the single soliton $\cos\theta=\pm\tanh[(s-s_0)/R]$ with $c^2/R^2=\omega_0^2-\omega^2$ (where $\theta$ is the polar angle of $\Omega$ , $s$ the comoving spatial coordinate, $s_0$ the soliton position, $R$ the width, and $\omega$ the soliton's internal azimuthal precession frequency), together with their scattering data — or establish non-integrability of the equations of motion for $\omega_0\neq 0$ . (b) Lorentz-breaking anisotropy: determine the consequences, for the soliton sector of the easy-axis model — in particular the resulting changes in the soliton rest energy, width, and dispersion — of the $\gamma_1=(\mu+\lambda)/2$ perturbation (with $\mu,\lambda$ the single-ion and exchange anisotropy parameters of the underlying lattice model) $H'=(\gamma_1 c/2)\int dx[gL^2+g^{-1}(\nabla\cos\theta)^2]$ , where $L$ is the azimuthal spin density, which the source paper introduces and then sets aside ("will not be considered further here"); the source-identified facts to build on are that this term breaks Lorentz invariance — the property the paper uses to obtain the moving finite- $P$ soliton by a boost of the static $P=0$ soliton — and that in the easy-plane limit it changes the collective-mode light velocity from $c$ to $(1+\gamma_1)^{1/2}c$ , with stability requiring $\gamma_1>-1$ . (c) Nonperturbative soliton spectroscopy: compute the renormalized soliton rest energies beyond leading semiclassical order and the precise critical value of the anisotropy at the integer- $S$ soliton-doublet to singlet transition (sharpening the parametric condition $(\hbar\omega_0)_c\sim\varepsilon_0$ ). In (a), "multisoliton" means the general $N$ -soliton solution for arbitrary soliton number $N$ ; in (b), the quantitative outcome is the change in the soliton rest energy, width, and dispersion to leading nontrivial order in $\gamma_1$ ; and (c) is posed for general spin $S$ . The fourth question raised in the same passage of the source — the required nonstandard quantization for half-integer $S$ — is resolved (topological $\theta=2\pi S$ term, Haldane's companion paper; Lieb-Schultz-Mattis theorem and its generalizations) and is not part of this problem.

References

F. D. M. Haldane, "Nonlinear Field Theory of Large-Spin Heisenberg Antiferromagnets: Semiclassically Quantized Solitons of the One-Dimensional Easy-Axis Néel State", Phys. Rev. Lett. 50, 1153-1156 (1983). DOI: 10.1103/PhysRevLett.50.1153. https://doi.org/10.1103/PhysRevLett.50.1153 (source paper)
F. D. M. Haldane, "Continuum dynamics of the 1-D Heisenberg antiferromagnet: Identification with the O(3) nonlinear sigma model", Phys. Lett. A 93, 464-468 (1983). DOI: 10.1016/0375-9601(83)90631-X. https://doi.org/10.1016/0375-9601(83)90631-X (companion paper introducing the theta=2piS topological term)
K. Pohlmeyer, "Integrable Hamiltonian systems and interactions through quadratic constraints", Commun. Math. Phys. 46, 207-221 (1976). DOI: 10.1007/BF01609119. https://doi.org/10.1007/BF01609119 (integrability of the massless O(3) sigma model, the omega_0=0 case)
J. D. Johnson, S. Krinsky, B. M. McCoy, "Vertical-Arrow Correlation Length in the Eight-Vertex Model and the Low-Lying Excitations of the X-Y-Z Hamiltonian", Phys. Rev. A 8, 2526-2547 (1973). DOI: 10.1103/PhysRevA.8.2526. https://doi.org/10.1103/PhysRevA.8.2526 (Bethe-ansatz solution including the massive regime of the S=1/2 XXZ chain)
Z.-P. Shi, G. Huang, R. Tao, "Envelope-kink excitations in a quantum anisotropic antiferromagnetic Heisenberg chain", Phys. Rev. B 43, 8583-8586 (1991). DOI: 10.1103/PhysRevB.43.8583. https://doi.org/10.1103/PhysRevB.43.8583
G. Gómez-Santos, "Variational approach to the XXZ spin-1 linear chain: Elementary excitations and Haldane conjecture", Phys. Rev. Lett. 63, 790-793 (1989). DOI: 10.1103/PhysRevLett.63.790. https://doi.org/10.1103/PhysRevLett.63.790
G. Gómez-Santos, "Haldane conjecture and elementary excitations in the XXZ spin-1 linear chain: A novel variational approach", J. Appl. Phys. 67, 5610-5612 (1990). DOI: 10.1063/1.345899. https://doi.org/10.1063/1.345899 (extended companion of the PRL above, on the 2D-Ising universality of the S=1 transition)
M. Lagos, G. G. Cabrera, "Coherent states as solutions of the anisotropic Heisenberg antiferromagnetic chain", Phys. Rev. B 38, 659-665 (1988). DOI: 10.1103/PhysRevB.38.659. https://doi.org/10.1103/PhysRevB.38.659 (a ground-state study, adjacent to but not resolving the soliton questions)
R. Balakrishnan, R. Blumenfeld, "On the twist excitations in a classical anisotropic antiferromagnetic chain", Phys. Lett. A 237, 69-72 (1997). DOI: 10.1016/S0375-9601(97)00685-3. https://doi.org/10.1016/S0375-9601(97)00685-3
D. C. Cabra, P. Pujol, C. von Reichenbach, "Non-Abelian bosonization and Haldane's conjecture", Phys. Rev. B 58, 65-68 (1998). DOI: 10.1103/PhysRevB.58.65. https://doi.org/10.1103/PhysRevB.58.65 (field-theory treatment for all S: integer gapped, half-integer gapless)
M. Oshikawa, "Commensurability, Excitation Gap, and Topology in Quantum Many-Particle Systems on a Periodic Lattice", Phys. Rev. Lett. 84, 1535-1538 (2000). DOI: 10.1103/PhysRevLett.84.1535. https://doi.org/10.1103/PhysRevLett.84.1535 (Lieb-Schultz-Mattis generalization)
M. B. Hastings, "Lieb-Schultz-Mattis in higher dimensions", Phys. Rev. B 69, 104431 (2004). DOI: 10.1103/PhysRevB.69.104431. https://doi.org/10.1103/PhysRevB.69.104431 (Lieb-Schultz-Mattis generalization)
D. Controzzi, G. Mussardo, "Mass Spectrum of the Two-Dimensional O(3) Sigma Model with a theta Term", Phys. Rev. Lett. 92, 021601 (2004). DOI: 10.1103/PhysRevLett.92.021601. https://doi.org/10.1103/PhysRevLett.92.021601 (form-factor analysis near theta=pi; double sine-Gordon effective theory)
V. Azcoiti, G. Di Carlo, A. Galante, "Critical Behavior of CP^1 at theta=pi, Haldane's Conjecture, and the Relevant Universality Class", Phys. Rev. Lett. 98, 257203 (2007). DOI: 10.1103/PhysRevLett.98.257203. https://doi.org/10.1103/PhysRevLett.98.257203 (numerical evidence for criticality at theta=pi)
B. Allés, A. Papa, "Mass gap in the 2D O(3) nonlinear sigma model with a theta=pi term", Phys. Rev. D 77, 056008 (2008). DOI: 10.1103/PhysRevD.77.056008. https://doi.org/10.1103/PhysRevD.77.056008 (masslessness at theta=3.10(5), consistent with pi)
G. Y. Cho, C.-T. Hsieh, S. Ryu, "Anomaly manifestation of Lieb-Schultz-Mattis theorem and topological phases", Phys. Rev. B 96, 195105 (2017). DOI: 10.1103/PhysRevB.96.195105. https://doi.org/10.1103/PhysRevB.96.195105 (modern 't Hooft-anomaly framing of the half-integer constraint)
S. B. Rutkevich, "Spinon confinement in the gapped antiferromagnetic XXZ spin-1/2 chain", Phys. Rev. B 106, 134405 (2022). DOI: 10.1103/PhysRevB.106.134405. https://doi.org/10.1103/PhysRevB.106.134405 (exploits the exact integrability of the S=1/2 chain in its gapped regime)
R. Blumenfeld, "Strange dynamics of domain walls and periodic stripes along classical antiferromagnetic chains", arXiv:cond-mat/0108470 (2001). DOI: 10.48550/arXiv.cond-mat/0108470. https://arxiv.org/abs/cond-mat/0108470 (unpublished preprint, no journal reference)
S. K. Kim, O. Tchernyshyov, Y. Tserkovnyak, "Thermophoresis of an antiferromagnetic soliton", Phys. Rev. B 92, 020402(R) (2015). DOI: 10.1103/PhysRevB.92.020402. https://doi.org/10.1103/PhysRevB.92.020402 (modern use of the single-soliton solution of the easy-axis antiferromagnet)
M. Daniel, A. R. Bishop, "Nonlinear excitations in the classical continuum antiferromagnetic Heisenberg spin chain", Phys. Lett. A 162, 162-166 (1992). DOI: 10.1016/0375-9601(92)90994-W. https://doi.org/10.1016/0375-9601(92)90994-W