---
title: Hermitian Bosonic Kitaev Chains
url: https://www.emergentmind.com/topics/hermitian-bosonic-kitaev-chains
type: topic
---

# Hermitian Bosonic Kitaev Chains

Hermitian bosonic Kitaev chains are one-dimensional pairing lattices whose second-quantized Hamiltonians are self-adjoint, yet whose bosonic Bogoliubov or dynamical matrices are generically non-Hermitian. In this sense, they are bosonic analogues of the fermionic Kitaev chain, but with a decisive structural difference: the relevant linear generator for bosonic Heisenberg dynamics can exhibit complex spectra, exceptional points, skin effects, and boundary-sensitive stability even when the many-body Hamiltonian satisfies \(\hat H^\dagger=\hat H\) [2607.08638]. A second, distinct usage of the term appears in interacting spin and fermion models that map exactly onto Hermitian Kitaev-chain forms while supporting emergent bosonic edge excitations, such as bosonic zero modes that are even under total fermion parity [2105.04326].

## 1. Canonical definitions and model classes

A canonical bosonic Kitaev chain is the nearest-neighbor quadratic pairing Hamiltonian
\[
\hat{H}_{\rm BKC} = \frac{1}{2}\sum_{j} \left( i w\,\hat{a}_{j+1}^{\dag}\hat{a}_{j} + i\Delta\,\hat{a}_{j+1}^{\dag}\hat{a}_{j}^{\dag} + {\rm H.c.} \right),
\]
with bosonic modes \(\hat a_j\), real \(w\), and real \(\Delta\) [2607.08638]. A more general translationally invariant family allows normal and pairing couplings at arbitrary range,
\[
\hat{H} = \frac{1}{2} \sum_{d\ge 0} \sum_{j} \Big( (g_d + i w_d)\,\hat{a}_{j+d}^{\dag}\hat{a}_j + (i\Delta_d + \eta_d)\,\hat{a}_{j+d}^{\dag}\hat{a}_j^{\dag} + {\rm H.c.} \Big),
\]
with real \(g_d,w_d,\Delta_d,\eta_d\); the bare BKC corresponds to \(w_d = w\,\delta_{d,1}\), \(\Delta_d = \Delta\,\delta_{d,1}\), and \(g_d=\eta_d=0\) [2607.08638].

The same label also covers several closely related Hermitian constructions. A minimal dimer version is
\[
H = \left(a^{\dag}+a\right)\left(b^{\dag}+b\right) + \mu\left(a^{\dag}a + b^{\dag}b\right),
\]
which already displays exceptional points in its effective bosonic core matrices [2502.19179]. A staggered chain introduces two sublattices and alternating hopping, pairing, and on-site potentials, leading to a \(4\times 4\) Bloch core matrix in Nambu space [2505.07017]. A modified two-sublattice bosonic Kitaev chain with intracell and intercell couplings yields an excitation Hamiltonian exactly mappable to a non-Hermitian SSH model [2505.15908]. In a different interacting direction, a spin-\(\tfrac12\) fermion chain with \(XY\) antiferromagnetic interaction maps onto a Hermitian Kitaev chain at half filling and supports a bosonic zero mode protected by spin-up parity [2105.04326].

This multiplicity of constructions suggests that “Hermitian bosonic Kitaev chain” is best understood as a family resemblance rather than a single model class. The common thread is a Hermitian many-body Hamiltonian with hopping-plus-pairing structure whose effective low-energy or dynamical description reproduces hallmark Kitaev phenomena in a bosonic setting.

## 2. Hermiticity of the Hamiltonian and non-Hermiticity of the dynamics

In this literature, “Hermitian” refers to the many-body Hamiltonian, not to the bosonic dynamical matrix. For quadratic bosonic systems written in Nambu form,
\[
\hat{H} =
\begin{pmatrix}
\hat{\mathbf a}^\dagger{}^T & \hat{\mathbf a}^T
\end{pmatrix}
H
\begin{pmatrix}
\hat{\mathbf a} \\
\hat{\mathbf a}^\dagger
\end{pmatrix},
\]
the dynamical matrix is
\[
M = i Z H,\qquad Z=\bigoplus_{j=1}^L \sigma_z,
\]
and the Heisenberg equations read \(i\,d\mathbf v/dt=M\mathbf v\) [2607.08638]. Even if \(H\) is Hermitian, \(M\) is generically non-Hermitian because of the bosonic symplectic structure. A recurring misconception is that Hermiticity of \(\hat H\) should force Hermiticity of the BdG generator; the bosonic case is precisely the counterexample [2607.08638, 2003.03405].

For the bare BKC, the momentum-space dynamical matrix is
\[
M_{\rm BKC}(k)= w \sin k\,\mathbb{I} + i \Delta \cos k\,\sigma_x,
\]
with eigenvalues
\[
E_\pm(k)= w\sin k \pm i\Delta\cos k,
\]
so periodic-boundary spectra are generically complex when \(\Delta\neq 0\) [2607.08638]. In the quadrature basis,
\[
\hat q_j=\frac{\hat a_j+\hat a_j^\dagger}{\sqrt2},\qquad
\hat p_j=-i\frac{\hat a_j-\hat a_j^\dagger}{\sqrt2},
\]
the BKC becomes block diagonal:
\[
M_{\rm BKC}^{\rm qp}(k)=
\begin{pmatrix}
w\sin k + i\Delta\cos k & 0\\
0 & w\sin k - i\Delta\cos k
\end{pmatrix},
\]
so the \(q\) and \(p\) quadratures evolve independently [2607.08638]. This decoupling is equivalent to an effective particle-hole symmetry of the dynamical matrix,
\[
M_{\rm BKC}(-k)^*=-M_{\rm BKC}(k),
\]
and each quadrature sector is a single-band non-Hermitian Hatano–Nelson problem [2607.08638].

The same structure underlies the phase-dependent chiral transport emphasized in the bosonic Kitaev-Majorana chain. In the quadrature representation, the equations of motion decouple and acquire asymmetric nearest-neighbor couplings, so one quadrature propagates preferentially in one direction and the other in the opposite direction; as \(t\to\Delta\), the transport becomes strictly one-way in the idealized limit [1805.12557]. The superconducting-circuit realization later observed this behavior directly as chiral transport and quadrature-dependent localization [2309.06178].

## 3. Symmetry structure, topology, and skin effects

Two effective symmetries organize much of the Hermitian bosonic Kitaev-chain phenomenology. The first is sublattice symmetry, defined at the second-quantized level by
\[
\hat{H}\big(\{\hat{a}_j\}\big) = -\,\hat{H}\big(\{(-1)^j \hat{a}_j\}\big).
\]
For translationally invariant bosonic pairing Hamiltonians, this symmetry implies an effective transpose-type time-reversal symmetry of the dynamical matrix. In the BKC, that relation can be written as
\[
\sigma_y\,M(k)\,\sigma_y = M(\pi-k),
\]
or, in the folded-zone formulation, \(T \mathcal M[k]^T T^{-1}=\mathcal M[-k]\) with \(TT^*=-\mathbb I\) [2607.08638]. This effective time-reversal symmetry is not physical time-reversal symmetry, but it supports a \(\mathbb Z_2\) invariant \(\nu_{\rm sl}(E)\) and a symmetry-protected skin effect [2607.08638].

The second organizing symmetry is quadrature particle-hole symmetry, which is equivalent to decoupling of \(q\) and \(p\) dynamics. When it holds, the problem reduces to two single-band non-Hermitian sectors with winding numbers
\[
W_{q/p}(E) = \frac{1}{2\pi i} \int_0^{2\pi} dk\,\partial_k \log\big(M_{q/p}(k)-E\big),
\]
and the BKC non-Hermitian skin effect can be understood as two decoupled Hatano–Nelson chains, one for each quadrature [2607.08638].

These symmetry principles clarify why Hermitian bosonic chains display non-Hermitian bulk-boundary phenomena. In the bare BKC, periodic-boundary eigenmodes are extended plane waves, whereas open-boundary eigenmodes become exponentially localized at one edge; at the same time, periodic boundaries can be dynamically unstable while open boundaries are completely stable over a broad parameter regime [2607.08638]. This boundary-sensitive stability was already identified in the bosonic Kitaev-Majorana chain, where the boundary-less system has delocalized unstable modes while a finite open chain can be dynamically stable and described by localized modes [1805.12557].

The modified bosonic Kitaev chain sharpens the topological content. For \(\omega=0\), its excitation Hamiltonian is exactly similar to a non-Hermitian SSH model with effective couplings
\[
\tilde\Delta_1=\sqrt{\Delta_1^2-J_1^2},\qquad
\tilde\Delta_2=\sqrt{\Delta_2^2-J_2^2},
\]
and the topologically nontrivial regime is
\[
|\tilde\Delta_2|>|\tilde\Delta_1|.
\]
In that regime, zero-energy edge modes coexist with a skin effect in the excitation Hamiltonian; finite onsite potential \(\omega\) rapidly destroys the skin effect and some, but not all, edge modes [2505.15908]. Disorder does not simply wash out this structure: at \(\omega=0\) the topological zero modes remain robust, and disorder in \(\omega\) can partially recover skin-effect-like localization at nonzero onsite potential [2505.15908].

## 4. Exceptional points, stability transitions, and dynamical braiding

Exceptional-point physics is central to Hermitian bosonic Kitaev chains because the non-Hermitian object is the dynamical matrix, not the many-body Hamiltonian. The minimal dimer provides the simplest analytic example. Its Nambu core matrices
\[
h^\pm = \frac{\mu\pm 1}{2}\sigma_z \pm \frac{i}{2}\sigma_y
\]
have eigenvalues
\[
\lambda_{1,2}^{\pm} = \pm \frac{1}{2}\sqrt{\mu^2 \pm 2\mu},
\]
with exceptional points at \(\mu=0,\pm2\) [2502.19179]. These EPs partition parameter space into four regions: two regions with two harmonic-oscillator sectors, and two mixed regions with one harmonic oscillator and one inverted harmonic oscillator. The nonequilibrium quantum phase transition is then detected through the second-order intensity correlation \(g^{(2)}(t_1,t_2)\); the derivative of its long-time average with respect to \(\mu\) develops pronounced valleys near \(\mu=\pm2\), which deepen with averaging time [2502.19179].

The extended bosonic Kitaev chain shows a related hidden-EP mechanism. For the periodic model
\[
H=\sum_{j=1}^{N}\Big[\,it\,b_{j+1}^{\dag }b_{j} +i\Delta\, b_{j+1}^{\dag }b_{j}^{\dag }+\text{H.c.} +\mu\left( 2b_{j}^{\dag }b_{j}+1\right)\Big],
\]
the Nambu core matrix
\[
h_k=
\begin{pmatrix}
\mu & i\Delta_k\\
i\Delta_k & -\mu
\end{pmatrix}
+T_k\,\mathbf 1,\qquad
T_k=t\sin k,\quad \Delta_k=\Delta\cos k,
\]
becomes non-diagonalizable at
\[
\mu=\Delta_k=\Delta\cos k.
\]
That exceptional point coincides with a localization-delocalization transition in an equivalent single-particle problem formulated in Fock space, obtained by projecting onto a BCS-like pairing basis [2410.15967]. In the localized regime \(\mu>|\Delta_k|\), a Bogoliubov vacuum exists and the inverse participation ratio stays finite; for \(\mu<|\Delta_k|\), the standard Bogoliubov vacuum ceases to exist and the equivalent Fock-space eigenstates delocalize [2410.15967]. The same hidden-EP scenario appears in the Dicke-model reduction to a two-site bosonic Kitaev model, where the average photon number after a quench from the empty state detects the transition [2410.15967].

A staggered bosonic Kitaev chain extends this logic to a \(4\times4\) non-Hermitian Bloch core matrix,
\[
h_{k}= \begin{pmatrix} 2g_{2} & i\Lambda_{k} & -iT_{k} & 0 \\
i\Lambda_{-k} & -2g_{1} & 0 & iT_{-k} \\
iT_{-k} & 0 & 2g_{1} & i\Lambda_{-k} \\
0 & -iT_{k} & i\Lambda_{k} & -2g_{2}
\end{pmatrix},
\]
with \(T_k=t_1+t_2 e^{ik}\) and \(\Lambda_k=\Delta_1+\Delta_2 e^{ik}\) [2505.07017]. In analytically tractable regimes, the EP conditions reduce to simple algebraic relations such as
\[
(g_1+g_2)^2=(\Delta_1\pm\Delta_2)^2
\]
or
\[
(t_1\pm t_2)^2=(\Delta_1\pm\Delta_2)^2,
\]
and these EP loci coincide with sharp localization-delocalization transitions of collective eigenstates in effective Fock-space networks, diagnosed by layer-resolved block inverse participation ratios [2505.07017].

A broader stability theory is supplied by pseudo-Hermiticity, generalized \(\mathcal{PT}\) symmetry, and Krein stability theory. For any quadratic bosonic Hamiltonian,
\[
G=\tau_3 H,\qquad G^\dagger=\tau_3 G\tau_3,
\]
so the effective BdG generator is pseudo-Hermitian and therefore generalized-\(\mathcal{PT}\)-symmetric [2003.03405]. Dynamical stability corresponds to the unbroken generalized-\(\mathcal{PT}\) phase, where \(G\) is diagonalizable and all eigenvalues are real. Instability appears either through exceptional points or through Krein collisions, where degenerate real eigenvalues split into complex-conjugate pairs while the matrix remains diagonalizable [2003.03405]. The Krein phase rigidity introduced there vanishes at both kinds of transition and extends standard non-Hermitian phase rigidity to bosonic indefinite-metric problems [2003.03405].

These non-Hermitian spectral features can be organized topologically. In multiband Hermitian bosonic Kitaev chains, complex dynamical eigenvalues can braid in the complex-frequency plane when loops in parameter space encircle exceptional points. Explicit two-strand and three-strand braids arise from square-root and cubic-root branch structures near second- and third-order EPs, and symmetry constrains the braids by pairing \(\lambda\) with \(-\lambda^*\) at fixed momentum and relating \(K\) to \(-K\) by complex conjugation [2509.18879].

## 5. Interacting, mapped, and spin-chain realizations

One major interacting realization begins from a spin-\(\tfrac12\) fermion chain with nearest-neighbor hopping and \(XY\) antiferromagnetic interaction. After a Majorana representation, a Jordan–Wigner transformation on a doubled chain, a Mattis–Nam transformation, and a second Jordan–Wigner transformation, the \(XY\) sector reduces exactly to
\[
H_0=-i\frac{J}{4}\sum_{j=1}^{L-1} (1-\gamma)a_j b_{j+1} -(1+\gamma) b_j a_{j+1},
\]
which is a Hermitian Kitaev chain at half filling and \(\mu=0\) [2105.04326]. The corresponding zero mode
\[
\eta_0=\frac{a_1+i b_L}{2}
\]
is “bosonic” in the sense that it is built from an even number of microscopic fermions and is even under total fermion parity \(P\), while remaining odd under spin-up parity \(P_\uparrow\), the protecting \(\mathbb Z_2\) symmetry of the topological phase [2105.04326]. In that model, the bosonic zero mode is accompanied by spontaneous breaking of \(P_\uparrow\), long-range end-to-end order in \(\langle S_1^x S_L^x\rangle\), and at least twofold entanglement-spectrum degeneracy in the ordered phase [2105.04326]. When the hopping term is turned on, the system undergoes a continuous transition to a trivial phase; at \(\gamma=0.4\), finite-size scaling gives \(\beta=\tfrac12\) and \(\nu=1\), consistent with the 2D classical Ising universality class [2105.04326].

A related but explicitly bosonic interacting construction is the paired Bose–Hubbard chain
\[
H = -w \sum_{\langle ij\rangle}(b_i^\dagger b_j + b_j^\dagger b_i)
+\sum_i\left[\frac{U}{2}\hat n_i(\hat n_i-1)-\mu \hat n_i\right]
-\Delta \sum_{\langle ij\rangle}(b_i^\dagger b_j^\dagger + b_i b_j),
\]
which supports a gapped \(\mathbb Z_2\) Ising phase with number fluctuation but no off-diagonal long-range order [2012.02380]. In the strongly interacting limit, the low-energy Hilbert space truncates to two occupations per site, yielding an anisotropic XY spin chain and, by Jordan–Wigner, the fermionic Kitaev chain. The bosonic and fermionic systems then share identical energy spectra, including the ground-state doublet, while their wavefunctions remain markedly different [2012.02380]. The bosonic doubly degenerate phase is described there as a gap-protected macroscopic qubit [2012.02380].

The bosonic Kitaev-Hubbard chain adds explicit on-site interactions and, in some cases, a three-body constraint. In the hard-core limit it maps to a spin-\(\tfrac12\) XY chain with Dzyaloshinskii–Moriya interaction,
\[
H_S = -\frac{1}{2}\sum_{r=1}^{L} \Big[(t\cos\theta + \Delta)\,\sigma_r^x \sigma_{r+1}^x + (t\cos\theta - \Delta)\,\sigma_r^y \sigma_{r+1}^y \Big]
+ \frac{1}{2}\sum_{r=1}^{L} t\sin\theta (\sigma_r^x \sigma_{r+1}^y - \sigma_r^y \sigma_{r+1}^x)
-\frac{1}{2} \sum_{r=1}^{L} \mu(1 - \sigma_r^z),
\]
and the non-Hermitian skin effect of the free-boson BdG description disappears exactly in this hard-core regime [2206.11827]. For chains with a three-body constraint and for unconstrained soft-core bosons, DMRG finds direct transitions between bond-pairing insulators and trivial insulators with correlation-length exponent \(\nu\approx 1\), again consistent with Ising criticality [2206.11827].

Spin-chain generalizations broaden the class of Hermitian bosonic Kitaev analogues. The spin-1 Kitaev-AKLT chain
\[
\hat{H}_{\theta} = \sum_{j} \Big[ K \big( \hat{S}_{2j}^{x}\hat{S}_{2j+1}^{x} + \hat{S}_{2j+1}^{y}\hat{S}_{2j+2}^{y} \big) + Q \big( (\hat{S}_{2j}^{x}\hat{S}_{2j+1}^{x})^{2} + (\hat{S}_{2j+1}^{y}\hat{S}_{2j+2}^{y})^{2} \big) \Big],
\]
with \(K=\cos\theta\) and \(Q=\sin\theta\), is Hermitian and exactly solvable at \(\theta=\pi/4\), where it becomes a sum of projectors onto maximal spin components along bond directions [2510.12880]. At that point it ունի an exponential ground-state degeneracy \(2^N+1\), admits a fractionalized spin-\(\tfrac12\) description, and each ground state can be written as a matrix product state [2510.12880]. This suggests that Hermitian bosonic Kitaev-chain physics extends naturally from quadratic bosons to constrained spin systems with projector structures, fractionalization, and exact tensor-network ground states.

## 6. Experimental implementations and broader significance

The most direct experimental realization to date is a multimode superconducting parametric cavity implementing the bosonic Kitaev chain in synthetic dimensions. There, lattice sites are encoded in cavity frequency modes, while complex hopping and pairing are generated in situ by parametric pumping at mode-difference and mode-sum frequencies, respectively [2309.06178]. The experiment demonstrated chiral transport, quadrature wavefunction localization, and sensitivity to boundary conditions, all identified as precursors of nontrivial topology and the non-Hermitian skin effect in a Hermitian bosonic chain [2309.06178].

Several other platforms are explicitly discussed in the literature surveyed here. The Hermitian bosonic dimer can be realized by two coupled cavities or mechanical resonators with both beam-splitter and two-mode-squeezing interactions, with \(\mu\) controlled by detunings or onsite frequencies; in that setting, second-order coherence \(g^{(2)}\) serves as an EP witness [2502.19179]. The hidden-EP chain analysis maps naturally onto the Dicke model, where the time evolution of the average photon number from the empty state detects the localization-delocalization transition associated with the exceptional point [2410.15967]. The modified bosonic Kitaev chain admits an optical-cavity implementation using alternating cavities embedded in a \(\chi^{(2)}\) nonlinear medium with parametrically modulated frequencies, giving tunable intracell and intercell hopping and pairing amplitudes [2505.15908]. More generally, optomechanical systems, circuit QED architectures, and cold-atom settings recur as natural platforms for realizing bosonic Kitaev chains and their generalizations [2502.19179, 2505.07017].

Two broad conclusions emerge from this body of work. First, Hermitian bosonic Kitaev chains provide controlled access to non-Hermitian topology without abandoning Hermitian many-body dynamics: skin effects, exceptional points, spectral winding, symmetry-protected skin effects, and dynamical braiding all arise from the bosonic dynamical matrix rather than from explicit gain-loss terms in the Hamiltonian [2607.08638, 2509.18879]. Second, interacting and mapped realizations show that “bosonic Kitaev physics” is not confined to quadratic bosons: it also appears in spinful fermion chains, paired Bose–Hubbard models, hard-core bosonic Kitaev-Hubbard systems, and spin-1 projector chains, where the protected low-energy object may be a bosonic zero mode, a bond-pairing Ising phase, or a fractionalized spin-sector edge degree of freedom rather than a literal bosonic Majorana excitation [2105.04326, 2012.02380, 2206.11827, 2510.12880].

In that broader sense, Hermitian bosonic Kitaev chains define a research program rather than a single model: they are Hermitian one-dimensional pairing systems whose bosonic structure turns the usual distinction between Hermitian and non-Hermitian physics into a problem of representation, symmetry, and stability.

Source: https://www.emergentmind.com/topics/hermitian-bosonic-kitaev-chains