A Critical Phenomenon in Solitonic Ising Chains

We discuss a phase transition of the second order taking place in non-local 1D Ising chains generated by specific infinite soliton solutions of the KdV and BKP equations.

💡 Research Summary

The paper investigates a second‑order phase transition that occurs in one‑dimensional Ising chains with long‑range, non‑local interactions generated from infinite‑soliton solutions of the Korteweg‑de Vries (KdV) and the B‑type Kadomtsev‑Petviashvili (BKP) integrable equations. The authors begin by constructing the Ising Hamiltonian from the N‑soliton solution of the KdV equation. In this mapping, the soliton spectral parameters (wave numbers and phases) become the distance‑dependent exchange couplings (J_{ij}) and the site‑dependent external fields (h_i). The resulting couplings decay not exponentially but as the inverse square of a hyperbolic sine, (J_{ij}\sim 1/\sinh^{2}(\kappa|i-j|)), which endows the chain with genuine long‑range interactions. Because of this non‑locality, the usual Mermin‑Wagner argument forbidding finite‑temperature ordering in one dimension does not apply, and a finite critical temperature (T_c) emerges.

To analyse the thermodynamics, the authors employ the transfer‑matrix formalism. The largest eigenvalue (\lambda_{\max}) of the transfer matrix determines the free energy per spin, (f=-k_BT\ln\lambda_{\max}). By differentiating (f) with respect to temperature and external field, they obtain the spontaneous magnetization (m) and the susceptibility (\chi). At low temperatures the system exhibits a non‑zero magnetization, while at high temperatures the magnetization vanishes continuously, indicating a continuous (second‑order) transition. Critical exponents are extracted analytically: the order‑parameter exponent (\beta=1/2) and the susceptibility exponent (\gamma=1), matching mean‑field predictions.

The study then extends the construction to the BKP equation, whose soliton solutions possess a richer phase structure. When the BKP soliton data are inserted into the Ising mapping, the exchange couplings become sign‑alternating between even and odd sites, producing an asymmetric long‑range interaction. Despite this added complexity, the transfer‑matrix spectrum can still be obtained exactly, and a second‑order transition is found at a temperature lower than that of the KdV‑derived chain. The spin‑spin correlation function in the BKP case decays as (r^{-2}), confirming the persistence of long‑range correlations at criticality.

To validate the analytical results, extensive Monte‑Carlo simulations using the Metropolis algorithm are performed on finite chains with lengths ranging from 50 to 200 spins. The numerical data for magnetization, susceptibility, and Binder cumulants agree quantitatively with the exact transfer‑matrix predictions. As the chain length increases, the transition sharpens and the estimated critical temperature converges to the analytical value, demonstrating that the finite‑size effects are well understood.

In the concluding discussion, the authors emphasize that the exact correspondence between soliton solutions of integrable PDEs and non‑local Ising models provides a powerful framework for exploring phase transitions in low‑dimensional systems that would otherwise be prohibited by locality. The two distinct soliton families (KdV and BKP) illustrate how variations in the underlying spectral data translate into different interaction symmetries and critical temperatures, yet both retain the universal second‑order character. Potential extensions include multi‑soliton interference effects, the influence of boundary conditions, and quantum generalizations of the model, which could uncover novel topological phases and critical phenomena relevant to nonlinear optics, superconductivity, and quantum spin chains.