Bosonic Kitaev Chain (BKC)

• Bosonic Kitaev Chain (BKC) is a quadratic bosonic lattice model featuring nearest-neighbor hopping and pairing that exhibits non-Hermitian dynamics and topological edge modes.
• The model leverages coherent hopping and parametric pairing to generate phenomena such as the non-Hermitian skin effect, exceptional points, and directional amplification.
• Experimental platforms including superconducting circuits, nano-optomechanical arrays, and cold-atom lattices validate its potential for quantum amplification and non-reciprocal transport.

The Bosonic Kitaev Chain (BKC) is a paradigmatic quadratic bosonic lattice model characterized by the interplay of nearest-neighbor hopping and onsite or nearest-neighbor two-boson pairing, realized in a variety of physical platforms ranging from nanomechanical arrays and superconducting circuits to cold-atom optical lattices. While formally analogous to the fermionic Kitaev chain—well known for its topological Majorana zero modes—the BKC exhibits distinct phases and dynamical phenomena rooted in the bosonic, non-number-conserving, and non-Hermitian effective dynamics induced by pairing. Its study provides a foundation for the investigation of non-Hermitian topological physics, non-reciprocal transport, skin effects, criticality, and the emergence of novel edge states in Hermitian systems of bosonic modes.

1. Model Definition and Canonical Hamiltonians

The canonical BKC consists of an array of bosonic modes with both coherent single-particle hopping and parametric (two-photon or two-phonon) pairing. In its minimal single-band, homogeneous form, the second-quantized Hamiltonian is

$$\hat H = \sum_{j=1}^{N-1} \left( i t\, \hat a_{j+1}^\dagger \hat a_j + i\Delta\, \hat a_{j+1}^\dagger \hat a_j^\dagger + \text{h.c.} \right),$$

where $\hat a_j$ denotes the annihilation operator at site $j$, $t$ is the hopping amplitude, and $\Delta$ is the (nearest-neighbor) parametric pairing amplitude. Various generalizations incorporate on-site energies, alternating couplings, multi-sublattice (SSH-type) unit cells, and time-periodic (Floquet) drives (McDonald et al., 2018, Slim et al., 2023, Busnaina et al., 2023, Fortin et al., 2024, Bomantara et al., 21 May 2025, Bomantara, 19 Nov 2025).

The model can equally be recast in quadrature variables: $$\hat x_j = \frac{\hat a_j + \hat a_j^\dagger}{2}, \qquad \hat p_j = \frac{\hat a_j - \hat a_j^\dagger}{2i}.$$ In this basis, for imaginary couplings, the Hamiltonian factorizes as

$$\hat H = \sum_{j=1}^{N-1}\left[ (t + \Delta)\, \hat x_j \hat p_{j+1} - (t - \Delta)\, \hat x_{j+1} \hat p_j \right].$$

This representation underlies the decoupling of the chain's dynamics into two quadrature sectors, which is essential for understanding the model's non-Hermitian properties (McDonald et al., 2018, Fortin et al., 2024).

2. Dynamical Matrix, Pseudo-Hermiticity, and Non-Hermitian Phenomena

The BKC's equations of motion for the Nambu spinor $$(\hat a_1, ..., \hat a_N, \hat a_1^\dagger, ..., \hat a_N^\dagger)^T$$ produce a Bogoliubov–de Gennes (BdG) matrix with an inherent pseudo-Hermitian structure: $\frac{d}{dt} \Psi = D \Psi,$ where $D$ is non-Hermitian due to the presence of pairing terms. Hermiticity of the underlying Hamiltonian ensures $J D J^{-1} = D^\dagger$ with $J = \mathrm{diag}(I_N, -I_N)$, meaning eigenvalues of $D$ come in $(\omega, -\omega^*)$ pairs (Wang et al., 23 Sep 2025). When $\Delta \ne 0$ and $t \cos \theta < \Delta$\,, bulk eigenvalues become complex, and exceptional points (EPs) emerge—these correspond to parameter values where eigenvalues and eigenvectors coalesce (Wang et al., 11 May 2025, Wang et al., 23 Sep 2025).

Non-Hermitian topology manifests in the BKC via point-gap winding invariants evaluated for the quasi-energy or dynamical matrix under periodic boundary conditions (PBC). The number of times the complex-energy band $E(k)$ encircles the origin defines a winding number $\nu$, which signals distinct topological phases and predicts the presence of skin modes and non-reciprocal amplification (McDonald et al., 2018, Fortin et al., 2024, Slim et al., 2023, Bomantara et al., 21 May 2025).

3. Non-Hermitian Skin Effect, Topological Invariants, and Edge Modes

The BKC generically displays the non-Hermitian skin effect: under open boundary conditions (OBC), all bulk eigenmodes localize exponentially at one edge. For the minimal chain, the localization length is given by $$\xi = \left[\ln\left(\frac{t+\Delta}{t-\Delta}\right)\right]^{-1}$$; the sign and location of localization depend on the quadrature and coupling asymmetry (Fortin et al., 2024).

The topological classification is governed by the point-gap winding: for Bloch Hamiltonian $E(k)$, the invariant is

$$\nu = \frac{1}{2\pi i} \int_{-\pi}^{\pi} \frac{d}{dk} \ln(E(k) - E_b) dk.$$

A nonzero $\nu$ under PBC implies, via bulk-boundary correspondence, the skin effect and OBC spectrum comprising edge-localized modes. Certain multi-sublattice or Floquet BKC generalizations map exactly to non-Hermitian SSH chains, where topological zero (and $\pi$) quasienergy modes are protected by winding numbers computed from analytic or Floquet block-decompositions (Bomantara et al., 21 May 2025, Bomantara, 19 Nov 2025).

The balanced pairing limit $t = \Delta$ realizes extreme non-reciprocity: propagation occurs strictly in one direction for $x$ and the opposite for $p$ quadratures, completely decoupling and isolating “Majorana-like” end-site quadratures that are exact zero modes under OBC (Slim et al., 2023).

4. Experimental Realizations and Observable Phenomena

Distinct physical signatures of the BKC include (i) phase-sensitive, chiral, directional amplification of coherent drives, (ii) exponential gain scaling $G^{N-1}$ with system size in the topological regime, and (iii) sharp sensitivity to boundary conditions (drastic spectrum changes between PBC/OBC, skin effect) (Slim et al., 2023, Busnaina et al., 2023, Fortin et al., 2024).

Implementation platforms include:

• Cavity/circuit QED and superconducting circuits: Frequency-modes of a parametric superconducting cavity realize the BKC in synthetic dimensions. Tunable parametric pumps at sum/difference frequencies engineer the hopping and pairing terms. Chiral transport, skin localization, and boundary condition sensitivity have been demonstrated experimentally (Busnaina et al., 2023).
• Nano-optomechanical arrays: Parametrically driven nanomechanical modes coupled via beamsplitter and two-mode squeezing interactions realize the BKC, enabling direct measurement of exponential gain and phase-dependent signal routing (Slim et al., 2023).
• Photonic and optical cavity arrays: Coupled resonators in nonlinear media, driven with site- and bond-modulated parametric processes, realize the model with tunable disorder and frequency detuning. Transmission/reflection spectra and local density-of-state measurements directly probe skin localization and zero modes (Bomantara et al., 21 May 2025).
• Ultracold atom systems: Zig-zag optical lattices with reservoir-induced pairing enable mapping to the BKC, with the possibility of accessing the full phase diagram including $\mathbb{Z}_2$ ordered (“Ising”) phases and Majorana end modes (Vishveshwara et al., 2020).

5. Extensions: Floquet, Staggered, and Interacting BKCs

Periodically driven (Floquet) bosonic Kitaev chains reveal an expanded topological landscape, supporting an arbitrary number of zero and $\pi$ edge modes depending on the Floquet winding numbers. Nontrivial combinations of intra- and intercell hoppings/pairings per drive half-period produce Floquet superoperators whose non-Hermitian skin effect and edge modes (with tunable multiplicity) persist up to moderate onsite frequency detuning or are regenerated by disorder (Bomantara, 19 Nov 2025).

Staggered (bi-sublattice) BKCs admit analytic criteria for exceptional points and localization-delocalization transitions of collective many-body states via the spectrum of a $4\times4$ Bloch core matrix. The position of EPs coincides exactly with sharp Fock-space transitions characterized using the layer-resolved inverse participation ratio, providing a practical tool for phase diagram mapping (Wang et al., 11 May 2025).

Introducing local interactions (two-/three-body repulsion) or constraining boson occupancy modifies the phase structure drastically. In the hard-core limit, the free BKC maps exactly to a transverse-field spin-$\frac12$ chain (additional Dzyaloshinskii–Moriya and compass terms for complex hopping/pairing phases) and, via Jordan–Wigner, to the fermionic Kitaev chain. In this interacting context, the non-Hermitian skin effect vanishes, and traditional many-body quantum criticality and symmetry-breaking phases dominate, with Ising universality and exact string-order parameters characterizing the transitions (Vishveshwara et al., 2020, Wang et al., 2022).

6. Effects of Dissipation, Disorder, and Topological Phase Transitions

The topological and amplification properties of the BKC exhibit sharp dependence on dissipative and disorder perturbations. Uniform onsite losses $\gamma$ enforce a topological transition at $\gamma_c = 4\Delta$; beyond this, exponential gain collapses and the point gap unwinds ($\nu \to 0$). More intricate loss patterns can enhance the system's resilience: chains of even length with loss on every other site preserve the topological phase for arbitrary loss values, while two-sublattice structures allow for compensation of large loss via a small counterpart on the alternate sublattice (Fortin et al., 2024). This tunability extends to disorder, where, in Floquet and modified BKCs, sufficient disorder restores partial skin localization or maintains edge modes even at nonzero on-site frequency (Bomantara et al., 21 May 2025, Bomantara, 19 Nov 2025).

The onset of complex-conjugate spectrum at EPs systematically matches the delocalization transitions in the many-body Hilbert space, providing a rigorous analytic framework for identifying critical points even in highly interacting and disordered regimes (Wang et al., 11 May 2025).

7. Relation to Fermionic Kitaev Chain and Broader Significance

Despite formal similarities, the BKC departs fundamentally from the fermionic Kitaev chain due to the absence of Pauli exclusion and the ensuing non-number-conserving (pairing) dynamics. In the BKC, pairing generically induces instability in periodic chains (complex eigenvalues, exceptional points) and the skin effect, while in the fermionic case the spectrum remains real and supports topologically protected Majorana modes with exponentially suppressed overlap (Wang et al., 23 Sep 2025, Vishveshwara et al., 2020). Strong interactions or hard-core constraints in the bosonic model, however, restore real spectra and enable a mapping to the fermionic topological paradigm.

The BKC presents a fertile setting for the physical realization and exploration of non-Hermitian topological phases, skin effect, edge zero-modes, and criticality in bosonic systems, supporting quantum amplification, enhanced sensing, and the study of complex band-structure topology in Hermitian parent models (McDonald et al., 2018, Slim et al., 2023, Fortin et al., 2024, Bomantara et al., 21 May 2025, Bomantara, 19 Nov 2025). The integration of parametric driving, dissipation engineering, periodic modulation, and strong interactions situates BKC physics at the intersection of quantum optics, condensed matter, and quantum information science, with ongoing experimental progress across multiple platforms.