---
title: 'Soliton Hopping: Discrete Nonlinear Transport'
url: https://www.emergentmind.com/topics/soliton-hopping
type: topic
---

# Soliton Hopping: Discrete Nonlinear Transport

Searching arXiv for the cited papers and closely related "soliton hopping" work.
Soliton hopping denotes a class of processes in which a localized nonlinear or topological excitation undergoes a discrete change of position, host, or dynamical state rather than a purely continuous drift. Across recent literature, the term has been used for adiabatic shuttling of topological mid-gap modes in a bosonic Su-Schrieffer-Heeger chain implemented in a Josephson-junction array, for electric-pulse-driven three-dimensional hops of Hopf solitons in chiral liquid crystals, for integer or fractional unit-cell transport of nonlinear lattice solitons under Thouless pumping, and for periodic or chaotic transfer phenomena involving dissipative Kerr solitons in coupled or injection-locked microresonators [2312.03456] [2112.14260] [2502.06131] [2508.09921] [2509.09108]. In this broad usage, the common element is discrete, protocol-dependent transport of a localized object together with the retention of a defining invariant such as a mid-gap branch, a Hopf index, a Chern number of a self-consistent effective Hamiltonian, or a periodic-orbit structure.

## 1. Scope and taxonomy

In the literature considered here, “soliton hopping” is not a single universal mechanism. In one-dimensional bosonic SSH systems, it refers to controlled motion of a topological zero-mode between neighboring dimers by local electrostatic gates. In chiral liquid crystals, it refers to a net three-dimensional displacement acquired after a full electric-field on/off cycle because the heliknoton \(\to\) hopfion and hopfion \(\to\) heliknoton transformations are non-reciprocal. In nonlinear AAH lattices, it denotes quantized adiabatic transport over one or several pump periods. In microresonator photonics, it can denote either periodic transfer between coupled resonators or random transitions between synchronized repetition-rate states [2312.03456] [2112.14260] [2502.06131] [2508.09921] [2509.09108].

| Platform | Localized object | Hopping signature |
|---|---|---|
| Josephson metamaterial | topological mid-gap soliton states | final fidelities \(F_{s_1},F_{s_2}\) exceed \(0.998\) |
| Chiral liquid crystal | heliknoton or hopfion | \(\Delta x\approx1\,\mu{\rm m}\), \(\Delta z\approx0.5\,p_0\) per cycle |
| Nonlinear off-diagonal AAH chain | soliton | average shift \(1\), \(1/2\), \(1/3\), or \(1/4\) unit cells per cycle |
| Coupled or driven microresonators | dissipative Kerr soliton | periodic transfer between resonators or random repetition-rate jumps |

A recurrent misconception is to equate hopping with ordinary translational motion of a continuum soliton. The surveyed works do not support that reduction. In several cases the process is instead branch-following in an instantaneous spectrum, a topology-preserving reconversion between distinct embedded textures, or a transition between coexisting attractors. This suggests that “hopping” functions as a phenomenological label for discrete transport of localized nonlinear states across markedly different dynamical settings.

## 2. Bosonic SSH soliton shuttling in a Josephson metamaterial

A particularly explicit formulation appears in the bosonic SSH realization derived from a Josephson-junction array. In units \(\hbar=2e=1\), the circuit Hamiltonian is
\[
H_{\rm JJA}=
\frac12\sum_{i,j}\tilde C^{-1}_{ij}(n_i-n_{g,i})(n_j-n_{g,j})
-\sum_{\langle i,j\rangle}J_{\langle i,j\rangle}\cos(\varphi_i-\varphi_j),
\]
with \(\tilde C_{ij}\) set by gate and coupling capacitances, \(n_{g,i}=-C_{g,i}V_{g,i}\), and \([\varphi_i,n_j]=i\delta_{ij}\). In the charging-dominated regime, linearizing Cooper-pair hopping yields a non-interacting bosonic tight-binding model,
\[
H_{\rm CPH}=-\sum_{i,j}t_{ij}b_i^\dagger b_j+\sum_iV_i n_i,
\]
which becomes an SSH chain under alternating couplings \(t_{2j-1,2j}=t_1\) and \(t_{2j,2j+1}=t_2\) [2312.03456].

The time-dependent driven Hamiltonian is
\[
H(t)=-\sum_{j=1}^N\Bigl[t_1\,b_{2j-1}^\dagger b_{2j}
+t_2\,b_{2j}^\dagger b_{2j+1}+\mathrm{h.c.}\Bigr]
+\sum_{i=1}^{2N}V_i(t)\,n_i,
\]
with only two adjacent islands pulsed to move the soliton. For the static chain with \(V_i=0\), the bulk bands obey
\[
\varepsilon_\pm(k)=\pm\sqrt{t_1^2+t_2^2+2t_1t_2\cos k},
\]
so the band gap is \(\Delta=|t_2-t_1|\). A domain wall or a strong detuning defect produces mid-gap states pinned near \(\varepsilon=0\). In the semi-infinite case, the zero-mode has support
\[
u_{2j-1}\propto\Bigl(-\frac{t_1}{t_2}\Bigr)^{j-1},\qquad
u_{2j}=0,\qquad \varepsilon=0,
\]
while on a finite ring with two detuned sites there is a pair of nearly degenerate mid-gap wavefunctions localized with length \(\xi_{\rm SSH}=\ln(t_2/t_1)\) [2312.03456].

The shuttling protocol moves the soliton from dimer \(j_0\to j_0+1\) by a three-segment sequence: ramp down \(V_{2j_0-1},V_{2j_0}\) from \(V_0\to0\), hold near zero for a buffer time, then ramp up \(V_{2j_0+1},V_{2j_0+2}\) from \(0\to V_0\). A convenient choice uses linear ramps of duration \(\tau\) and a relative shift \(\delta\approx\tau\), with amplitudes \(1.1V_0\) and \(0.9V_0\) on the two sites of a dimer so that the soliton never sees two strong gates at once. The evolution is tracked in the instantaneous basis \(|u^{(\alpha)}(t)\rangle\) and propagated by \(U(t)=T\,e^{-i\int^tH(t')dt'}\).

Adiabaticity is controlled by
\[
\max_t\frac{|\langle u^{(m)}(t)|\dot H(t)|u^{(s)}(t)\rangle|}
{[\varepsilon_m(t)-\varepsilon_s(t)]^2}\ll1,
\]
with the rough requirement \(\|\dot H\|\sim V_0/\tau\ll\Delta^2\). The fidelity of a mid-gap branch,
\[
F_\alpha(t)=\bigl|\langle u^{(\alpha)}(t)|U(t)|u^{(\alpha)}(0)\rangle\bigr|^2,
\]
satisfies
\[
F(T)=1-\mathcal O\!\bigl((V_0/(\Delta^2\tau))^2\bigr),
\]
and a Landau-Zener estimate gives \(F\approx1-\exp[-\pi\Delta^2/(2\,V_0/\tau)]\). Numerical benchmarks used a ring of \(N=21\) dimers, \(t_1=1\), \(t_2=2\), \(V_0/t_1=100\), small sub-dimer splitting \(\pm10\%\), and effective temperature \(T/t_1=5\times10^{-4}\). With total single-hop time \(t_{\rm prot}=3.17\times10^3\,t_1^{-1}\), protocol \(P1\) keeps the mid-gap states away from bulk crossings and yields final fidelities exceeding \(0.998\), whereas protocol \(P2\) without the proper time-shift fails with \(F\to0\) [2312.03456].

## 3. Three-dimensional hopping of Hopf solitons

In chiral liquid crystals, soliton hopping is formulated for three-dimensional topological textures characterized by the Hopf invariant. The director field \(\mathbf n(\mathbf r)\), with \(|\mathbf n|=1\), is governed by a Frank-Oseen free energy
\[
F[\mathbf n]=\int d^3r\,\bigl(f_{\rm Elastic}[\mathbf n]+f_{\rm Dielectric}[\mathbf n]\bigr),
\]
where
\[
f_{\rm Elastic}=K_{11}(\nabla\!\cdot\!\mathbf n)^2
+K_{22}\bigl(\mathbf n\!\cdot\!(\nabla\times\mathbf n)+q_0\bigr)^2
+K_{33}|\mathbf n\times(\nabla\times\mathbf n)|^2
\]
and
\[
f_{\rm Dielectric}=-\frac12\epsilon_0\Delta\epsilon(\mathbf E\!\cdot\!\mathbf n)^2.
\]
The Hopf charge is the linking number of preimage loops and may be written as
\[
Q=H=\frac{1}{32\pi^2}\int \epsilon_{ijk}A_iF_{jk}\,d^3r,
\qquad
F_{ij}=\mathbf n\cdot(\partial_i\mathbf n\times\partial_j\mathbf n)
\]
[2112.14260].

Two embedded textures are emphasized. In the uniform background with \(\mathbf n\to\mathbf e_z\), a localized Hopf texture is termed a hopfion. In a helical background winding around an axis \(X_0\) with period \(p_0\), the corresponding embedded object is a heliknoton, with a background ansatz
\[
\mathbf n_{\rm helix}(z)=
\bigl(\sin\theta_0\cos(2\pi z/p),\,
\sin\theta_0\sin(2\pi z/p),\,
\cos\theta_0\bigr).
\]
Stability depends strongly on elastic anisotropy. By mixing a bent-core liquid crystal with a rod-like liquid crystal, the ratio \(k=K_{33}/K_{22}\) can be tuned from \(\approx0.4\) to \(\approx1.5\), and numerical minimization shows that for \(k\le1.2\) and zero field heliknotons are energetically metastable, and in certain \(k\)-ranges even thermodynamically favored, without any applied field or cell confinement [2112.14260].

The hopping mechanism is driven by electric pulses. With \(E\) on, the unwound background is favored and the soliton deforms into a hopfion in a uniform surround. When \(E\) is switched off, the helical background is restored and the texture becomes a heliknoton again. The critical point is that the director evolution is non-reciprocal: the path from helix to uniform to helix does not retrace itself. Each full on/off cycle therefore produces a net displacement while the Hopf index remains rigorously conserved at every moment. The evolution is modeled by overdamped relaxation,
\[
\frac{\partial\mathbf n}{\partial t}=-\Gamma\,\frac{\delta F}{\delta\mathbf n},
\]
with the unit-length constraint and projected torques [2112.14260].

The reported transport is discrete and three-dimensional. Under periodic on/off pulses with period \(T\lesssim3\,{\rm s}\), the soliton drifts by \(\Delta x\approx1\,\mu{\rm m}\), rotates its long axis by \(\Delta\phi\approx10^\circ\) per cycle, and exhibits a vertical hop \(\Delta z\approx0.5\,p_0\) each cycle. The full heliknoton \(\to\) hopfion \(\to\) heliknoton cycle takes \(O(10\text{–}50\,{\rm s})\), and the two half-cycles are asymmetric, with hopfion \(\to\) heliknoton taking \(\approx3\times\) longer than the reverse. Simulations typically used a box of \(4p_0\times4p_0\) in \(x,y\), height \(d=1.7\text{–}4p_0\), a mesh of \(24\) grid points per \(p_0\), periodic boundaries in \(x,y\), strong Dirichlet boundaries in \(z\), and material parameters \(K_{11}=6.4\,{\rm pN}\), \(K_{22}=3\,{\rm pN}\), \(K_{33}=1.2\text{–}1.5\,K_{22}\), \(\Delta\epsilon=13.8\), \(p_0\approx2.33\,\mu{\rm m}\) [2112.14260].

## 4. Fractional Thouless pumping of nonlinear lattice solitons

A different use of soliton hopping arises in the nonlinear off-diagonal Aubry-André-Harper model, where transport is induced by adiabatic modulation of nearest-neighbor hoppings. The linear Hamiltonian is
\[
H^{\rm lin}(\phi)=\sum_n\Bigl[-t_{n,n+1}(\phi)|n\rangle\langle n+1|
-t_{n-1,n}(\phi)|n\rangle\langle n-1|
+m_n|n\rangle\langle n|\Bigr],
\]
with
\[
t_{n,n+1}(\phi)=J_0+J_a\cos\!\Bigl(\phi-\frac{2\pi q}{p}(n-1)\Bigr),
\qquad \phi\in[0,2\pi].
\]
After adding on-site Kerr nonlinearity, the amplitudes obey a discrete nonlinear Schrödinger equation,
\[
i\frac{d\psi_n}{dt}
=\sum_{n'}H_{n,n'}^{\rm lin}(\phi(t))\psi_{n'}
+g_n|\psi_n|^2\psi_n,
\]
with conserved norm \(N=\sum_n|\psi_n|^2\) [2502.06131].

The central mechanism is self-induced topology. In the instantaneous approximation, the stationary soliton profile \(\psi_n(\phi)\) satisfies a nonlinear eigenproblem with eigenvalue \(\mu(\phi)\). Its density produces an effective on-site shift
\[
\widetilde m_n(\phi)=g_n|\psi_n(\phi)|^2,
\]
so that the modified linear Hamiltonian
\[
H^{\rm sc}(\phi)=H^{\rm lin}(\phi)+\sum_n\widetilde m_n(\phi)|n\rangle\langle n|
\]
can carry nontrivial topology even though the bands of \(H^{\rm lin}\) are topologically trivial. The pump parameter is ramped as
\[
\phi(t)=\frac{2\pi}{T}t,\qquad T\gg1,
\]
ensuring adiabatic evolution [2502.06131].

In a supercell of length \(L\), one defines a Chern number
\[
C=\frac{1}{2\pi}\int_0^{2\pi}d\phi
\int_{\rm BZ}dk^{\rm sc}\,\Omega_{k^{\rm sc},\phi}.
\]
If after \(m\) cycles the soliton has moved by exactly \(n\) unit cells, the average displacement per cycle is
\[
\Delta x_{\rm avg}=\frac{n}{m}\ \text{unit cells}.
\]
The reported cases realize one-cell transport in one, two, three, or four periods, corresponding to average shifts of \(1\), \(1/2\), \(1/3\), or \(1/4\) unit cells per cycle. The soliton center is tracked by
\[
x_c=\sum_n n\,|\psi_n|^2/N,
\]
and exhibits perfectly quantized plateaux of \(\Delta x=1/k\) per cycle for \(k=1,2,3,4\) [2502.06131].

Representative parameter sets are explicitly given. For the integer pump, \(p=3\), \(q=1\), \(N=2.3\), \(J_a=0.05\), \(m=(-0.2,-0.2,0.04)\), and \(g=(-0.8,-0.8,-1)\), and the soliton moves exactly one unit cell as \(\phi\) goes from \(0\) to \(2\pi\). For the half-cell pump, \(p=5\), \(q=2\), \(N=2.2\), \(J_a=0.05\), \(m=(-0.2,-0.2,0.2,-0.02,0.2)\), \(g=(-0.67,-0.67,-1.33,-1.18,-1.33)\), and one unit cell is traversed in two cycles. Analogous one-third-cell and one-quarter-cell pumps are reported for \(p=7,q=3\) and \(p=9,q=4\), respectively. The proposed implementation is a one-dimensional array of evanescently coupled nonlinear waveguides, with typical parameters \(J_0\approx0.15\,{\rm mm}^{-1}\), modulation period \(T\approx(2\text{–}8)\times10^3\), and physical sample lengths of order \(10\text{–}50\,{\rm m}\) per pump cycle [2502.06131].

## 5. Dissipative Kerr-soliton hopping in microresonators

In coupled Kerr-nonlinear ring resonators, soliton hopping appears as a periodic orbit of coupled Lugiato-Lefever equations rather than as adiabatic motion in a static spectral gap. For a dimer or trimer of identical resonators, the normalized mean-field equations are
\[
\partial_t\Psi_j=-(1+i\zeta)\Psi_j
+i\,d_2\,\partial_\phi^2\Psi_j
+i|\Psi_j|^2\Psi_j
-i\sum_{k=1}^N J_{jk}\Psi_k+f_j.
\]
A stationary supermode soliton can lose stability when a pair of complex-conjugate eigenvalues crosses the imaginary axis, producing a periodic orbit in which the soliton periodically hops between resonators. For the dimer with \(J=10\) and \(f^2=64\), the stationary branch \(SS^{(2)}\) is stable in a window of \(\zeta\) and undergoes a subcritical Hopf at \(\zeta_H^{(2)}\approx2.77\), with \(\omega_H\approx33.3\). The resulting branch \(SH^{(2)}\) initially emerges unstable and becomes stable only after saddle-node bifurcations of the orbit. For the trimer at the same \(J\) and \(f^2\), the stationary branch \(SS^{(3)}\) is already unstable to a real eigenvalue, and a supercritical Hopf at \(\zeta_H^{(3)}\approx5.79\) creates an unstable \(SH^{(3)}\) that stabilizes after a subsequent fold. These bifurcation differences explain why trimer hopping can be experimentally accessible at lower pump power than dimer hopping, and why the dimer exhibits hysteresis, multistability, and strong path dependence in laser scans [2508.09921].

A distinct microresonator realization concerns chaotic group-velocity hopping of a single dissipative Kerr soliton under Kerr-induced synchronization to an externally injected weak reference laser. The underlying single-ring dynamics follow the Lugiato-Lefever equation
\[
\frac{\partial a(\theta,t)}{\partial t}
=-\Bigl(\frac{\kappa}{2}-i\delta\omega_0\Bigr)a
+i\mathcal D(a)-i\gamma|a|^2a+iF_0.
\]
With a reference laser near comb mode \(\mu_s\), synchronization fixes
\[
\dot\varphi=\omega_{\rm ref}-(\mu_s\omega_{\rm rep}+\omega_0)=0,
\]
and phase modulation \(\phi_{\rm ref}(t)\to\phi_{\rm ref}(t)+A\cos(\omega_{\rm rf}t)\) allows locking to the carrier or sidebands. The reduced dynamics are described by a second-order Adler equation,
\[
\beta\frac{d^2\varphi}{d\tau^2}+\frac{d\varphi}{d\tau}
=\alpha+\sin\varphi,
\]
and, under modulation, by the autonomous three-dimensional system
\[
\frac{d\varphi}{d\tau}=\Omega,\qquad
\frac{d\psi}{d\tau}=\Omega_{\rm ext},\qquad
\beta\frac{d\Omega}{d\tau}=-\Omega+\alpha+\sin(\varphi+A\cos\psi).
\]
Because the second-order term supplies effective inertia, the modulated system supports multiple coexisting attractors and strange attractors, yielding chaotic switching between repetition-rate states [2509.09108].

The experiment used a \(23\,\mu{\rm m}\) radius \({\rm Si_3N_4}\) microring pumped in TE mode at \(282.3\,{\rm THz}\) with \(\sim150\,{\rm mW}\) on chip, a counterpropagating TM cooler laser at \(307.6\,{\rm THz}\), and a weak TE reference at \(192.5\,{\rm THz}\) phase-modulated up to \(850\,{\rm MHz}\). Detection employed EO-comb bridging of \(56\) sidebands at \(\sim17.84\,{\rm GHz}\). Near \(\alpha\approx0.58\) for \(\Omega_{\rm ext}=850\,{\rm MHz}\) and \(A\approx1\), the theory yields a positive Lyapunov exponent \(\lambda\approx2\pi\times5.3\times10^{-2}\) and the \(0\text{–}1\) chaos test gives \(K\approx1\). Experimentally the RF spectrum shows two sharp tones separated by
\[
\Delta\omega_{\rm rep}=\frac{\omega_{\rm rf}}{\mu_s}\approx2\pi\times9.44\,{\rm MHz},
\]
and time-resolved spectrograms reveal random jumps between them over a \(100\,{\rm MHz}\)-wide detuning range. Binary tests on the hopping sequence yield \(p\)-values all \(>0.8\) [2509.09108].

## 6. Related discrete transport problems and conceptual distinctions

A related but not identical framework is provided by one-dimensional mechanical linkage models. There the chain consists of rigid rotors with angles \(\theta_i(t)\), arm lengths \(r_i\), pivot spacings \(a_i\), and springlike constraints encoded by the discrete Lagrangian
\[
L=\sum_{i=1}^N\frac12Mr_i^2\dot\theta_i^2-\sum_{i=1}^{N-1}V_i(\theta_i,\theta_{i+1}),
\]
with
\[
V_i(\theta_i,\theta_{i+1})
=\kappa_i\bigl[l_i^2(\theta_i,\theta_{i+1})-\bar l_i^2\bigr]^2.
\]
Topological zero-energy solitons are defined as trajectories on the configuration-space manifold \(C=\{\theta\in\mathbb T^N:l_i(\theta_i,\theta_{i+1})=\bar l_i\ \forall i\}\), and scattering at a junction between chains with different parameters yields reflection and transmission amplitudes
\[
R=\frac{Z_2-Z_1}{Z_1+Z_2},\qquad
T=\frac{2Z_2}{Z_1+Z_2},
\]
with an effective impedance \(Z_i\). If both sides have the same topological index \(\nu_i=\operatorname{sign}[1-(r_i/a_i)\sin^2\bar\theta_i]\), the soliton is fully transmitted; if \(\nu_1\neq\nu_2\), the zero-energy cycle does not exist across the junction and the soliton is totally reflected [1804.06007].

This scattering problem clarifies an important distinction. Not every discrete change in a soliton’s position is a driven hop in the sense of the later SSH, Hopf-soliton, or microresonator literature. In the linkage model, the central issue is transmission or reflection at an interface; in the Josephson SSH model it is adiabatic branch-preserving shuttling; in chiral liquid crystals it is non-reciprocal cycle-to-cycle displacement with strictly conserved Hopf index; in the nonlinear AAH problem it is quantized topological pumping mediated by a self-induced potential; and in photonic microresonators it is either periodic resonator transfer or chaotic switching between synchronized states [1804.06007].

Taken together, these works indicate that the most stable cross-disciplinary meaning of soliton hopping is discrete transport of a localized nonlinear object under a protocol or bifurcation structure that suppresses ordinary dispersive delocalization. What is preserved differs from system to system: a mid-gap identity in bosonic SSH chains, a linking invariant \(Q\) in Hopf textures, a Chern-number-controlled displacement law in nonlinear pumping, or an invariant periodic orbit or synchronization branch in dissipative photonics. This suggests that the term is best understood not as a single mechanism but as a family resemblance across topological, nonlinear, and driven-dissipative transport phenomena.

Source: https://www.emergentmind.com/topics/soliton-hopping