Emergent $D_8^{(1)}$ spectrum and topological soliton excitation in CoNb$_2$O$_6$

Quantum integrability emerging near a quantum critical point (QCP) is manifested by exotic excitation spectrum that is organized by the associated algebraic structure. A well known example is the emergent $E_8$ integrability near the QCP of a transverse field Ising chain (TFIC), which was long predicted theoretically and initially proposed to be realized in the quasi-one-dimensional (q1D) quantum magnet CoNb$_2$O$_6$. However, later measurements on the spin excitation spectrum of this material revealed a series of satellite peaks that cannot be described by the $E_8$ Lie algebra. Motivated by these experimental progresses, we hereby revisit the spin excitations of CoNb$2$O$6$ by combining numerical calculation and analytical analysis. We show that, as effects of strong interchain fluctuations, the spectrum of the system near the 1D QCP is characterized by the $D{8}^{(1)}$ Lie algebra with robust topological soliton excitation. We further show that the $D{8}^{(1)}$ spectrum can be realized in a broad class of interacting quantum systems. Our results advance the exploration of integrability and manipulation of topological excitations in quantum critical systems.

💡 Research Summary

This paper revisits the long‑standing puzzle of the excitation spectrum in the quasi‑one‑dimensional Ising magnet CoNb₂O₆. While early experiments suggested that the emergent E₈ integrability of the transverse‑field Ising chain (TFIC) could be observed near the one‑dimensional quantum critical point (QCP), later high‑resolution THz and neutron scattering measurements revealed numerous satellite peaks that cannot be accommodated within the E₈ mass hierarchy.

The authors construct a realistic minimal model that captures the essential intra‑chain XXZ Ising physics together with the frustrated inter‑chain Ising couplings characteristic of CoNb₂O₆’s isosceles triangular lattice. By treating two neighboring chains exactly (cluster mean‑field) and the rest at a mean‑field level, the Hamiltonian reduces to a two‑leg ladder with a longitudinal effective field ˜h originating from the surrounding chains. In the limit where the inter‑chain coupling J_i dominates over ˜h (J_i/˜h → ∞), the low‑energy field theory becomes two coupled critical Ising conformal field theories (c = 1/2 each) perturbed by a bilinear σ₁σ₂ term. This perturbation is known to be integrable and is described by the affine Lie algebra D₈⁽¹⁾.

The D₈⁽¹⁾ spectrum consists of a single soliton (and its antisoliton) of mass m_s and six breathers B_n with masses m_n = 2 m_s sin(nπ/14) for n = 1…6. Importantly, only the even‑parity breathers couple to the ground state, while odd‑parity breathers are symmetry‑forbidden, a prediction confirmed by the numerical data.

Using infinite‑time‑evolving block decimation (iTEBD), the authors compute the dynamical structure factor of the ladder at the critical transverse field (H₁D^c) with J_i = 0.36 J and ˜h = 0. The resulting k = 0 spectrum displays sharp peaks precisely at the D₈⁽¹⁾ masses, including a prominent low‑energy peak around 0.74 meV identified as a topological soliton—an excitation corresponding to a single domain wall that is normally forbidden but becomes robust due to the inter‑chain confinement. Multiparticle continua and higher‑breather thresholds are also reproduced, matching the complex satellite structure observed experimentally.

A direct comparison with THz absorption data shows that most experimental peaks up to ~2 meV align with the D₈⁽¹⁾ mass predictions. Two residual peaks not captured by the pure ladder model are explained by (i) a zone‑folding effect caused by next‑nearest‑neighbor antiferromagnetic coupling J_AF, and (ii) additional bound states induced by the domain‑wall interaction λ_dw. These features are model‑specific and do not alter the underlying D₈⁽¹⁾ algebraic organization.

Field‑dependence analysis reveals that the soliton‑to‑second‑breather mass ratio (m₁ + m_s)/m₂ ≈ 1.665 remains nearly constant over a finite field range, reflecting the topological protection of the soliton. In contrast, other mass ratios vary more strongly with field, consistent with their non‑topological nature. The scaling of masses follows m_n ∝ J^{4/7}, as expected for the D₈⁽¹⁾ integrable theory, distinct from the E₈ scaling m_n ∝ J^{8/15}.

The authors argue that the dominance of strong, frustrated inter‑chain fluctuations near a three‑dimensional QCP naturally drives the system into the D₈⁽¹⁾ universality class, a scenario that should be realizable in a broad family of quasi‑1D quantum magnets with similar geometry. The presence of a robust topological soliton and symmetry‑restricted breathers offers potential avenues for quantum information applications, as these excitations are protected against local perturbations.

In conclusion, the paper provides compelling theoretical and numerical evidence that CoNb₂O₆ does not realize the E₈ spectrum but instead exhibits an emergent D₈⁽¹⁾ integrable structure with a distinctive set of soliton and breather excitations. This work expands the landscape of emergent integrability in quantum critical systems and opens new directions for exploring topological quasiparticles in low‑dimensional magnets.