---
title: Extended Kitaev Chain Models
url: https://www.emergentmind.com/topics/extended-kitaev-chain
type: topic
---

# Extended Kitaev Chain Models

The extended Kitaev chain is a family of one-dimensional Kitaev-chain models in which the standard spinless \(p\)-wave superconducting chain is generalized beyond the nearest-neighbor, uniform, Hermitian setting. In the literature cited here, the extension appears in several distinct but related forms: longer-range hopping and pairing, next-nearest-neighbor terms, dimerized or trimerized unit cells, spatially varying bond phases, quantum-dot or Andreev-bound-state inhomogeneities, and non-Hermitian deformations. The common consequence is a richer bulk-edge structure than in the basic nearest-neighbor chain, including multiple topological sectors, multicriticality, massive edge modes, non-endpoint Majorana bound states, and, in broader spin-chain usages, Jordan–Wigner-related Majorana descriptions of bond-dependent spin models [1703.10086] [2307.11996] [2310.05006].

## 1. Canonical extended-range formulation

A canonical fermionic definition is the spinless \(p\)-wave superconducting chain with couplings beyond nearest neighbors,
\[
H= -\sum_{j=1}^{L}\mu\left(a_j^\dagger a_j-\frac12\right) +\sum_{\ell=1}^{r}\sum_{j=1}^{L-\ell} \left( -w_\ell e^{i\varphi_\ell} a_j^\dagger a_{j+\ell} +\Delta_\ell a_j a_{j+\ell} +\text{h.c.} \right),
\]
where \(r\) is the maximum number of coupled neighbors, \(w_\ell\) and \(\Delta_\ell\) are hopping and pairing amplitudes, and \(\varphi_\ell\) can break time-reversal symmetry [1703.10086]. A closely related truncated-range formulation writes
\[
H = - \mu \sum_{j = 1}^N c_j^\dagger c_j^{\,} - \sum_{j = 1}^{N - 1} \sum_{\ell =1}^{q} \big( t_\ell^{} c_j^\dagger c_{j+\ell}^{\,} - \Delta_\ell^{} c_{j}^\dagger c_{j+\ell}^\dagger + \mbox{H.c.} \big),
\]
with algebraic couplings
\[
\Delta_\ell = \Delta \left(\frac{1}{\ell}\right)^\alpha, \qquad t_\ell = t \left(\frac{1}{\ell}\right)^\beta,
\]
so that \(\alpha\) and \(\beta\) control the decay of pairing and hopping [2606.10043].

Within this range-extended setting, the truncated-range scenario has as many distinct topological phases as the number of coupled neighboring sites [2606.10043]. A particularly important special case is the extended Kitaev chain with next-nearest-neighbor hopping and pairing, where the real-space Hamiltonian includes nearest-neighbor amplitudes \(\lambda_1\), next-nearest-neighbor amplitudes \(\lambda_2\), and, in the non-Hermitian deformation, a complex chemical potential \(\mu+i\gamma\) [2307.11996]. In that model the next-nearest-neighbor term enlarges the topological structure from the usual two-phase picture of the basic Kitaev chain to three phases with winding numbers
\[
W=0,\;1,\;2
\]
[2307.11996].

This family already shows that “extended” is not merely quantitative. Adding further-neighbor couplings changes the allowed winding sectors, the number of phase boundaries, and the multiplicity of edge modes. In the short-range limit \(\alpha,\beta\to\infty\), these formulations reduce to the standard nearest-neighbor Kitaev chain [1703.10086] [2606.10043].

## 2. Bulk topology and phase structure

For the finite-range time-reversal-symmetric chain, the Bloch/BdG Hamiltonian can be written as
\[
\mathcal{H}(k)= h_y(k)\sigma_y+h_z(k)\sigma_z,
\]
with
\[
h_z(k)= -\frac{\mu}{2}-\sum_\ell w_\ell \cos(\varphi_\ell)\cos(k\ell),\qquad
h_y(k)= \sum_\ell \Delta_\ell\sin(k\ell),
\]
and the corresponding BDI invariant is the winding number
\[
{\sf w} = \frac{1}{2\pi}\oint d\theta_k
\]
[1703.10086]. When time-reversal symmetry is broken by \(\varphi_\ell\neq 0\), the symmetry class reduces to class D and the relevant invariant becomes
\[
\upsilon=\operatorname{sign}\!\bigl(h_z(0)\,h_z(\pi)\bigr)
\]
[1703.10086].

A central finite-range result is that, with time-reversal symmetry preserved and with both hopping and pairing extended to \(r\) neighbors, one can have as many Majorana zero modes per edge as the number of coupled neighbors \(r\) [1703.10086]. By contrast, if only hopping or only pairing is extended, one still gets at most one Majorana zero mode per edge [1703.10086]. In the specific uniform-coupling case analyzed there, the winding number satisfies
\[
{\sf w}=1 \quad \text{for} \quad -2r<\mu/w_0\le 0,\qquad
{\sf w}=r \quad \text{for} \quad 0<\mu/w_0<2
\]
[1703.10086].

Dimerized extensions reorganize the same topological content in a different way. For the SSH-Kitaev chain at \(\mu=0\), an exact real-space extended Bogoliubov transformation decomposes a \(2N\)-site chain into two independent \(N\)-site Kitaev chains, with the full winding number given by
\[
\mathcal N=\mathcal N_A+\mathcal N_B,
\]
so that the full chain supports phases with \(\mathcal N=0,1,2\) [2310.05006]. The phase boundaries follow from the standard subchain condition \(\mu_\sigma/J_\sigma=\pm1\), which in SSH-Kitaev parameters gives four critical curves:
\[
\Delta-2\lambda+1=0,\quad 2\lambda \Delta-\Delta+1=0,\quad
\Delta+2\lambda-1=0,\quad 2\lambda \Delta-\Delta-1=0
\]
[2310.05006].

Long-range models add further distinctions. In the infinite-range cases, standard bulk invariants do not always fully distinguish between massless Majorana zero modes and massive edge modes, so exact diagonalization and finite-size scaling are additionally required [1703.10086]. In the time-reversal-symmetric long-range hopping-and-pairing case, the critical lines are
\[
\mu_{c1}=-2w_0\zeta(\alpha), \qquad \mu_{c2}=2w_0\eta(\alpha),
\]
while for sufficiently long-ranged couplings the formal winding number can become \(\pm \tfrac12\) because of nonanalyticity at \(k=0\) [1703.10086].

## 3. Edge modes, finite-size diagnostics, and spatial distribution

The extended Kitaev chain is defined as much by its edge phenomenology as by its bulk winding. In the generalized transfer-matrix treatment for arbitrary neighbor range \(r\), zero-energy wavefunctions take the form
\[
\phi_{o,j}=\sum_{s=1}^{2r} c_s^\phi \lambda_s^j,\qquad
\psi_{o,j}=\sum_{s=1}^{2r} c_s^\psi \lambda_s^{-j},
\]
and the number of independent zero-mode solutions per edge is
\[
N=\max(n^<,n^>)-r \le r,
\]
where \(n^<\) and \(n^>\) count roots inside and outside the unit circle [1703.10086]. This real-space result matches the bulk statement that the maximum number of Majorana zero modes per edge is controlled by the interaction range [1703.10086].

Extended chains can also support massive edge modes. In the infinite-range models studied in the same work, these are edge-localized states with finite energy that remain separated from the bulk by a finite gap, and they can coexist with Majorana zero modes in crossover regions [1703.10086]. This sharpens the distinction between “edge-localized” and “zero-energy Majorana” in long-range superconducting chains.

A more recent finite-chain analysis introduces two additional diagnostics: the Majorana average position and the occupation of edge-to-edge non-local fermion states, defined as an effective parity [2606.10043]. That work reports a direct correlation between the ground-state fermion parity and the edge occupation numbers, which are translated into localization and delocalization of the Majorana average position [2606.10043]. In that formulation, the distribution of Majorana modes is not reducible to the bulk winding number alone, but is tied to how non-local fermion occupation reorganizes edge localization in finite systems [2606.10043].

Operational diagnostics can also be formulated directly in terms of end operators. For the nearest-neighbor open Kitaev chain, the end-to-end observable
\[
\hat Q=2(\hat c_1\hat c_N^\dagger+\hat c_N\hat c_1^\dagger)
\]
defines the thermodynamic indicator
\[
Z=\lim_{N\to\infty} |\langle \hat Q\rangle|,
\]
with the empirical expression
\[
Z=\max\left[ \frac{4|w\Delta|}{(|\Delta|+|w|)^2} \left(1-\left(\frac{\mu}{2w}\right)^2\right),\,0 \right]
\]
[1804.10584]. That work explicitly notes that this formula is not derived for extended chains and should not be transferred unchanged to longer-range hopping or pairing without separate derivation [1804.10584].

## 4. Dimerized, trimerized, phase-textured, and inhomogeneous variants

One major branch of the subject replaces the uniform lattice by an enlarged unit cell. In the SSH-Kitaev chain, the exact decomposition into two subchains gives a direct explanation of why the full system can host \(0\), \(1\), or \(2\) Majorana zero modes per edge [2310.05006]. The same work also introduces two BCS-pair order parameters, \(O_{A,\mathrm g}\) and \(O_{B,\mathrm g}\), associated with the two decomposed sectors, so that the full ground state is a product of two BCS condensates [2310.05006].

A trimerized extension goes further. The three-sublattice model obtained by hybridizing a modified SSH chain with trimerized unit cells and the standard Kitaev chain has a hexamer structure in the Majorana basis and supports not only Majorana zero modes but also nonzero-energy edge-localized modes that are still Majorana in nature [2509.04420]. Its BdG Hamiltonian has a \(6\times 6\) structure, remains in class BDI, and uses a winding number
\[
\mathcal{W} = \frac{1}{2 \pi i} \int_{-\pi}^{\pi} \text{Tr} \left( h^{-1}(k)\, \frac{\partial h(k)}{\partial k} \right) dk
\]
for zero-energy topology, while open-boundary spectra reveal the additional \(\pm E\) edge-mode pairs [2509.04420].

A different extension is generated not by unit-cell enlargement but by a spatial phase texture. In the linearly phase-modulated chain
\[
H_{C}=-\mu\sum_{n=1}^{N}c_{n}^{\dag}c_{n}-t\sum_{n=1}^{N-1}e^{i\phi_n}(c_{n}^{\dag}c_{n+1}+h.c.) +\Delta \sum_{n=1}^{N-1}e^{i \theta_n}(c_{n}c_{n+1}+h.c.),
\]
a hopping phase difference \(\delta\phi\) generates zero-energy bound states localized at non-endpoint sites [2403.02023]. In the finite-size calculations of that work, twofold-degenerate zero-energy states first emerge at \(\delta\phi=0.25\pi\); for \(\delta\phi=\pi\) there are two non-endpoint Majorana bound states, for \(\delta\phi=2\pi\) fourfold-degenerate Majorana bound states appear, and for \(\delta\phi=3\pi\) sixfold-degenerate Majorana bound states appear [2403.02023]. In ring geometries, the same phase control can move interior Majoranas and implement braiding without T/Y junctions [2403.02023].

An inhomogeneous end segment provides yet another extension. For a finite Kitaev chain coupled to a quantum dot through a step-like potential \(V(x)\) and position-dependent pairing \(\Delta(x)\),
\[
\tilde{H}_{BdG}= -\left(\partial_x^2+\mu-V\left(x\right)\right)\tau_z-i\Delta\left(x\right)\partial_x\tau_y,
\]
a robust Andreev bound state localized in the quantum-dot region appears as the generic lowest-energy solution in the topologically trivial phase [1808.02495]. That state does not exist in the bare uniform chain and is attributed to a partial decoupling of the component Majorana bound states over the length of the dot potential [1808.02495]. The same work stresses that the signatures of this trivial ABS in local tunneling are identical to the signatures of topologically protected Majorana zero modes [1808.02495].

## 5. Non-Hermitian and fractional generalizations

Non-Hermitian extensions have produced a distinct body of “extended Kitaev chain” results. One route is asymmetric pair creation and annihilation,
\[
\mathcal{H} =-t\sum_{j=1}^{N}(c_{j}^{\dag }c_{j+1}+\text{H.c.}) -\mu \sum_{j=1}^{N}\left( 1-2n_{j}\right) -\sum_{j=1}^{N}(\Delta _{\alpha }c_{j}^{\dag }c_{j+1}^{\dag } +\Delta _{\beta }c_{j+1}c_{j}),
\]
with \(\Delta_\alpha\neq \Delta_\beta\) [1707.04718]. In the unbroken time-reversal-symmetric region the spectrum is real, the topology is characterized by a biorthogonal extended Zak phase,
\[
\mathcal Z_\pm= \begin{cases} -\pi\,\mathrm{sgn}\!\left(\dfrac{\Delta_\alpha+\Delta_\beta}{t}\right), & |\mu|<1,\\ 0, & |\mu|>1, \end{cases}
\]
and Majorana edge zero modes survive under open boundary conditions [1707.04718]. The topological phase is destroyed in the broken region, where eigenstates coalesce at exceptional points [1707.04718].

A more explicitly range-extended non-Hermitian chain adds next-nearest-neighbor hopping and pairing together with a complex chemical potential \(\mu+i\gamma\) [2307.11996]. There the next-nearest-neighbor term produces the three phases \(W=0,1,2\), while the non-Hermitian factor qualitatively reorganizes the critical structure. The gap-closing momentum moves continuously along the critical line, multiple gapless points can appear, one Hermitian critical line collapses into an isolated critical point, and the number of multicritical points is reduced [2307.11996]. The paper also reports unconventional critical exponents,
\[
z \simeq 0.5,\qquad \nu \simeq 1,\qquad \gamma_{CE}\simeq 1,
\]
in sharp contrast to the Hermitian extended Kitaev chain, where different multicritical points can have \(z=1\) or \(z=2\) [2307.11996].

A more algebraic extension replaces ordinary Pauli matrices in the BdG blocks by rational powers, producing a centrally extended Clifford algebra and a pseudo-metallic regime with rational-valued winding number [2204.13104]. In that construction the critical twist is
\[
\gamma^*=\frac12,
\]
the bulk gap closes at
\[
\nu=\frac12,
\]
and the topological Majorana boundary mode of the untwisted chain evolves into an extended midgap state in the pseudo-metallic phase [2204.13104]. The paper interprets this as a projective or fractional extension of Kitaev-chain topology rather than a longer-range coupling extension [2204.13104].

## 6. Artificial realizations, transport, and engineered chain links

Artificial implementations have made the extended Kitaev chain a device-level concept. In an alternating array of normal quantum dots and proximitized hybrid segments hosting Andreev bound states, the low-energy spin-polarized dot orbital and the spin-polarized ABS quasiparticle become the two sites of an effective two-site Kitaev chain,
\[
H_{K2} = \varepsilon_1 f^\dagger_1 f_1 + \varepsilon_2 f^\dagger_2 f_2 + t_{12} f^\dagger_2 f_1 + \Delta_{12} f^\dagger_2 f^\dagger_1 + h.c.,
\]
with
\[
t_{\mathrm{eff}}=-tu,\qquad \Delta_{\mathrm{eff}}=-t_{so}v
\]
[2309.15777]. When the system is scaled to three sites, second-order virtual tunneling generates next-nearest-neighbor couplings,
\[
H^{(2)}_{\text{DAD,eff}} = t_{DD} d^\dagger_{L\downarrow} d_{R\downarrow} + \Delta_{DD} d_{L\downarrow} d_{R\downarrow} + h.c.,
\]
so the longer device is explicitly an extended Kitaev chain rather than an ideal nearest-neighbor one [2309.15777]. The same work reports that the effective coupling between sites at distance \(k\) scales as
\[
\Gamma_{k} \sim \frac{t^k_0}{(2E_Z)^{k-1}},
\]
so the longer-range terms are exponentially short-ranged [2309.15777].

A complementary implementation is the flux-controlled two-site Kitaev chain, where two spin-polarized quantum dots are coupled through extended ABSs in a flux-tunable Josephson junction [2501.15912]. The phase difference
\[
\Delta \phi = \phi_2 - \phi_1
\]
controls the ABS energy and coherence factors, thereby tuning the ECT-to-CAR ratio [2501.15912]. The poor man’s Majorana sweet spot is identified by
\[
\Gamma_{\mathrm{e}}=\Gamma_{\mathrm{o}},
\]
and the paper reports a continuous sweet-spot line in \((V_\mathrm{ABS},\Delta\phi)\) space rather than isolated points [2501.15912]. An additional middle probe detects a zero-bias conductance peak from the hybrid region, indicating that the two poor man’s Majorana wavefunctions \(\gamma_1\) and \(\gamma_2\) both reside partially in the ABS region [2501.15912].

Transport calculations through longer-range chains connected to metallic leads show that range extension and time-reversal-symmetry breaking are not interchangeable perturbations. In the N-TS-N junction built from a chain with
\[
\gamma_l=\gamma_0 l^{-\alpha},\qquad \Delta_l=\Delta_0 l^{-\beta},\qquad \phi_l=\phi_0 l,
\]
the short-range proxy \(\alpha=\beta=10\) and the long-range case \(\alpha=\beta=0.5\) respond differently as \(\phi_0\) varies [2606.30483]. In the short-range chain a quantized zero-bias peak of height \(2e^2/\hbar\) survives for \(\phi_0=0\) and \(\phi_0=\pi/4\), weakens at \(\phi_0=\pi/3\), and disappears when the bulk gap closes at \(\phi_0=\pi/2\) [2606.30483]. In the long-range chain there is no quantized zero-bias peak even at \(\phi_0=0\), because long-range interactions convert zero-energy Majoranas into massive subgap Dirac modes, and for nonzero \(\phi_0\) the subgap states disappear from the gap region [2606.30483].

## 7. Related spin-chain usages and broader generalizations

In a broader but well-established usage, “extended Kitaev chain” also refers to bond-dependent spin chains whose Jordan–Wigner images are Majorana or BdG chains. The spin-\(\tfrac12\) Kitaev-XX-\(\Gamma\) chain is one such example. After Jordan–Wigner transformation it becomes a quadratic spinless-fermion Hamiltonian with a two-site unit cell, and the exact solution yields six phases: four gapped ordered phases and two gapless phases [2509.21901]. The gapless phases contain two branches of helical Majorana fermions, and the transition lines include deconfined quantum critical lines with \(z=1\) and quadratic critical lines with \(z=2\) [2509.21901].

A more materials-oriented example is the \(JK\Gamma\) chain proposed for CoNb\(_2\)O\(_6\), where the dominant microscopic exchange is ferromagnetic Kitaev and a sizable antiferromagnetic \(\Gamma\) term fixes the easy axis [2403.14754]. In that work the chain is described as a one-dimensional Kitaev chain in disguise as an Ising chain, with local-frame Hamiltonian
\[
H_{ij}^0 = J \mathbf s_i\cdot \mathbf s_j + K s_i^\gamma s_j^\gamma + \Gamma(s_i^\alpha s_j^Z+s_i^Z s_j^\alpha)
\]
and fitted parameters
\[
J=-0.8~{\rm meV},\qquad K=-1.1~{\rm meV},\qquad \Gamma=0.56~{\rm meV}
\]
[2403.14754].

Coupled-chain extensions of the spin Kitaev-\(\Gamma\) chain show a different kind of enlargement. In the anisotropic honeycomb model with \(d=0\) corresponding to decoupled Kitaev-\(\Gamma\) chains, the emergent \(SU(2)_1\) Tomonaga–Luttinger liquid persists for finite interchain coupling and develops into an extended quantum spin liquid with spinon-like excitations [2212.11000]. That phase differs from the conventional Kitaev spin liquid and is explicitly compared with sliding Luttinger liquids [2212.11000].

Other spin-chain descendants probe transport and entanglement rather than bulk topology. The Heisenberg–Kitaev chain interpolates between a one-dimensional spin Kitaev chain and an XXZ Heisenberg chain, and its high-temperature energy transport shows ballistic response at the integrable points, a topological gap near the Kitaev point, and quantum-chaotic regions away from integrability [1312.4954]. The spin-\(\tfrac12\) twisted Kitaev chain and generalized bond-dependent XY chain can be rewritten in terms of Majorana degrees of freedom and driven by local pulse sequences into rainbow-like states of long-distance entangled Majorana pairs, with maximal bipartite entanglement entropy and mirror-symmetric Majorana pairing [2301.10777].

These spin-chain usages are not identical to the standard fermionic extended Kitaev chain with longer-range \(p\)-wave hopping and pairing. They nevertheless belong to the same conceptual field because they preserve the defining combination of one-dimensionality, bond-dependent or BdG-like structure, and Majorana-based descriptions of edge, bulk, or entanglement phenomena [2509.21901] [2403.14754] [2212.11000] [1312.4954] [2301.10777].

Source: https://www.emergentmind.com/topics/extended-kitaev-chain