Analytical diagonalization of the open-boundary bosonic Kitaev chain: An asymmetric plane-wave ansatz approach

arXiv:2608.29684 · quant-ph, math-ph, math.MP · Submitted 2026-08-30 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: Today's paper: "Analytical diagonalization of the open-boundary bosonic Kitaev chain".

Mira: The open-boundary bosonic Kitaev chain (OBKC) is a model of interest due to its realization in driven-dissipative systems and its non-Hermitian boundary physics, despite being governed by a Hermitian Hamiltonian.

Kai: First, who's behind it and why it matters.

Title and authors: Kai: So we're looking at a paper called "Analytical diagonalization of the open-boundary bosonic Kitaev chain: An asymmetric plane-wave ansatz approach." It seems like they're tackling a model that has been around for a while because it shows up in driven-dissipative systems, even though the underlying Hamiltonian is Hermitian.

Mira: That title tells us immediately that they aren't just doing another numerical diagonalization; they are proposing an entirely different algebraic way to solve this specific problem, which is quite intriguing.

Lev: From a hardware standpoint, if this method works exactly for any finite chain length N, that takes the pressure off needing massive computational resources just to get spectral information.

Kai: Exactly. The authors are proposing an asymmetric plane-wave ansatz with unequal left- and right-moving momenta to solve the 2N times 2N non-Hermitian associated matrix, which is a big step toward finding analytical solutions for these boundary conditions.

Mira: They are aiming for a purely algebraic solution using standard bosonic Bogoliubov transformation, which is appealing because it bypasses some of the geometric intuition involved in other methods.

Lev: If they can provide exact closed-form expressions for the eigenvalues and eigenvectors for any finite N, that’s what makes this relevant for quantum error correction researchers like myself; we need to know exactly what the system supports before we try to map it onto physical qubits.

Kai: It sounds like they are trying to provide a self-contained solution that doesn't rely on transforming into the position-momentum representation, which is something I find really useful for understanding how these systems behave in real experiments.

Mira: And they’re going ahead and discussing the implications right away, suggesting this new approach addresses some of the non-uniqueness issues seen in other methods.

The paper's summary: Kai: So, to summarize what they’ve done with this paper, they’re presenting an alternative method for solving the open-boundary bosonic Kitaev chain that relies only on the standard bosonic Bogoliubov transformation. They use an asymmetric plane-wave ansatz with unequal left- and right-moving momenta to tackle the 2N times 2N non-Hermitian associated matrix.

Mira: The core of their summary is showing how this ansatz leads to N distinct eigenvalues, each twofold degenerate, and then they show a rigorous algebraic way using bosonic commutation relations to select only the single physical Bogoliubov mode from those two eigenvectors.

Lev: That degeneracy resolution part is crucial because in real hardware, we can't just pick an eigenvector randomly; we need a concrete way to construct the actual quasiparticle operators.

Kai: They then explicitly construct these N physical Bogoliubov quasiparticle operators, and they show that the diagonal form of the Hamiltonian they get is equivalent to the original one.

Mira: The paper also points out that this construction reveals non-uniquenesses that in special cases perfectly mirror the freedom we see in local squeezing transformations, which is a very important point for theorists looking at these systems.

Lev: So, it’s not just finding the modes; it's proving that there is only one physical way to build them up from the mathematical solutions they found. That kind of rigorous proof is what we need for reliable implementations.

The paper's improvements: Kai: The authors highlight several improvements in their approach, focusing on how this method stands apart from previous work, especially when compared to the local squeezing transformation method they reviewed.

Mira: They specifically state that their proposed asymmetric plane-wave ansatz and degeneracy-resolution technique are not restricted just to the OBKC model, but can be generalized to other bosonic pairing systems, including those with inhomogeneous hopping or pairing terms.

Lev: Generalizing beyond the specific chain structure is huge; if this algebraic technique works for more complex models, it means we don't have to re-derive everything from scratch every time we move to a slightly different physical setup.

Kai: They also emphasize that their method works directly with creation and annihilation operators without having to transform to the position-momentum representation or invoking the generalized Brillouin zone formalism, which I think is much more practical for experimentalists.

Mira: And they provide exact closed-form expressions for eigenvalues and eigenvectors for any finite chain length N, rather than just focusing on the continuous spectrum in the thermodynamic limit that some other theories address.

Lev: Having those exact expressions for finite N is incredibly useful because it gives us concrete numbers we can compare against simulations or experimental measurements on finite samples.

Conclusion: Kai: So, wrapping up the paper "Analytical diagonalization of the open-boundary bosonic Kitaev chain: An asymmetric plane-wave ansatz approach," they’ve shown a purely algebraic path to solving this model that yields exact closed-form expressions for eigenvalues and eigenvectors for any finite chain length N.

Mira: The main implication is that they provide a self-contained solution to the problem, explicitly resolving the degeneracy issue in a mathematically transparent manner using bosonic commutation relations.

Lev: For error correction, this means we have a precise tool for constructing the physical modes without relying on ambiguous transformations like local squeezing.

Kai: They also showed that the resulting N Bogoliubov transformations contain N/two

(N − one)/two: real free parameters, which they link back to the arbitrary constant lambda in the original squeezing transformation method as a special case.

Mira: That non-uniqueness is shown to reflect that the system doesn't have a well-defined ground state, which is a critical physical caveat for theorists considering these models.

Lev: I just want to reiterate that the method’s ability to handle inhomogeneous systems and provide exact finite N results makes it a very robust tool for testing the limits of our current theoretical frameworks.

Center for Quantum Technology Research · Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements (MOE) · Center for the Cross-disciplinary Research of Space Science and Quantum-technologies (CROSS-Q)

quant-ph, math-ph, math.MP

Submitted: 2026-08-30

Updated: 2026-09-28

Comments: 13 pages, 1 figure

Journal ref: J. Phys. A: Math. Theor. 59 395201 (2026)

DOI: 10.1088/1751-8121/aea940

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 79/100

The gist: The open-boundary bosonic Kitaev chain (OBKC) is a model of interest due to its realization in driven-dissipative systems and its non-Hermitian boundary physics, despite being governed by a Hermitian

Key concepts

Open-Boundary Bosonic Kitaev Chain (OBKC)
This is a model of interest in driven-dissipative systems that has a Hermitian Hamiltonian but exhibits non-Hermitian boundary physics. It is relevant because it appears in physical systems and is studied for its boundary conditions.
Asymmetric Plane-Wave Ansatz
The authors propose using an ansatz with unequal left- and right-moving momenta to solve the 2N times 2N non-Hermitian associated matrix. This method aims to provide an algebraic solution without needing transformations into the position-momentum representation.
Bogoliubov Transformation
This is a standard bosonic transformation used by the authors to select a single physical Bogoliubov mode from the two degenerate eigenvectors found by their ansatz. This step is crucial for constructing concrete quasiparticle operators.
Degeneracy Resolution
The method shows how to rigorously select one physical Bogoliubov mode from the N distinct eigenvalues, each of which is twofold degenerate. This selection is done using bosonic commutation relations to ensure a concrete way to build the actual physical operators.

Terminology

Summary

The open-boundary bosonic Kitaev chain (OBKC) is a model of interest due to its realization in driven-dissipative systems and its non-Hermitian boundary physics, despite being governed by a Hermitian Hamiltonian. The paper presents an alternative, purely algebraic solution that relies entirely on the standard bosonic Bogoliubov transformation. For an N-site chain, the authors propose an asymmetric plane-wave ansatz with unequal left- and right-moving momenta to analytically solve the associated 2N × 2N non-Hermitian “associated matrix.” The left eigenvalue problem yields N distinct eigenvalues, each of which is twofold degenerate. By carefully resolving these degeneracies using the bosonic commutation relations, they construct the N physical Bogoliubov quasiparticle operators. The construction reveals non-uniquenesses that in special cases exactly mirror the freedom in the local squeezing transformations of the original approach. The diagonal form of the Hamiltonian is obtained explicitly and is shown to be equivalent to the original Hamiltonian. The proposed asymmetric plane-wave ansatz and degeneracy-resolution technique are not limited to this model and can be generalized to other bosonic pairing systems, including those with inhomogeneous pairing or hopping terms.

The paper contrasts its method with the original solution, which employs a local squeezing transformation in the position-momentum representation:

"Under the squeezing transformation, the Hamiltonian is converted to HB = N/X−1 j=1 p δ+δ− x˜jp˜j+1 − δ− δ− x˜j+1p˜j = (PN−1 j=1 p δ+δ−i(˜a† j+1a˜j − a˜† j a˜j+1), δ− > 0, PN-1 j=1 p −δ+δ−i(˜a† j+1ㆠj − ãjãj+1), δ− < 0."

The authors note that the original approach is indirect: "While elegant, this approach relies on the (x, p)-representation and the geometric intuition of “squeezing”, which, from a purely algebraic standpoint, is somewhat indirect and tricky. The method seems simpler since some elaborate manipulations have been hidden in the squeezing, gauge, and Fourier transformations."

The proposed new approach consists of five steps:

  1. We cast the Hamiltonian into the standard bosonic quadratic form and perform a conventional Bogoliubov transformation, leading to a non-Hermitian 2N × 2N associated matrix M [8].

  2. "To solve the left eigenvalue problem of M, ϕkM = Λkϕk, we propose an asymmetric plane-wave ansatz with unequal left- and right-moving momenta, which is a natural generalization of the symmetric plane waves used for Hermitian tridiagonal systems [9–12]."

  3. "Substituting the ansatz into the bulk equations and boundary conditions yields explicit quantization conditions. Solving these conditions explicitly gives the discrete set of allowed momenta and the closed-form eigenvalues."

  4. "A key technical subtlety emerges: each eigenvalue is twofold degenerate in the left-eigenspace, with two linearly independent eigenvectors. We demonstrate, through a rigorous algebraic proof, that the bosonic commutation relations allow only a single physical Bogoliubov mode to be constructed from the two degenerate eigenvectors. This resolves the degeneracy in a systematic and mathematically transparent manner."

  5. "The resulting N Bogoliubov transformations contain N/2 [(N − 1)/2] real free parameters for even (odd) N. This non-uniqueness is shown to be equivalent to the arbitrary constant λ in the squeezing transformation method as a special case, and reflects the fact that the system has no well-defined ground state."

The method distinguishes itself by working directly with creation and annihilation operators without transforming to the (x, p)-representation or invoking the generalized Brillouin zone formalism. Compared to previous works, it provides exact closed-form expressions for eigenvalues and eigenvectors for any finite chain length N, rather than the continuous spectrum in the thermodynamic limit. Furthermore, it explicitly resolves the degeneracy and discusses non-uniqueness.

The analysis of F1 (related to one set of modes) shows that the allowed momenta are determined by Eq. (49): αm = −π/2 + mπ/N + 1, m = 1, 2,..., N. The resulting physical Bogoliubov mode is given by Eq. (50): Vj (αm) = 2(−i) j−1 sin πmj/N + 1 e −βj. Similarly, for F2, the eigenvector is given by Eq. (51): "Wj (αm) = 2(−i) j−1 sin πmj/N + 1 e βj.

Improvements for AI systems

Based on the provided scientific paper, here are specific improvements for AI systems, categorized by the capabilities they could gain:


)Improved AI System Capabilities:

  1. Generation of Novel Topological/Non-Hermitian Quantum Architectures (Quantum Simulation & Materials Science):

  2. Solving Complex Non-Hermitian Boundary Value Problems (Analytical/Numerical Solver):

  3. Modeling Driven-Dissipative Open Systems (Dynamics & Transport Modeling):

  4. Developing Robust, Non-Bloch Band Theory Frameworks:

)Specific Improvements for AI Systems:

  1. AI can design and predict the behavior of synthetic quantum materials (like bosonic Kitaev chains) with engineered non-Hermitian boundary conditions (inhomogeneous hopping/pairing).

  2. The AI can analytically determine the exact spectral properties (eigenvalues and eigenvectors) of these complex systems, bypassing computationally expensive numerical diagonalization for finite chains.

  3. The system can model and predict the transport characteristics (e.g., chiral transport, non-Hermitian skin effects) in driven-dissipative quantum devices with high precision.

  4. The AI can develop generalized mathematical frameworks (like the asymmetric plane-wave ansatz) to solve quadratic Hamiltonians beyond simple, periodic systems, specifically those with inhomogeneous boundary conditions or non-Bloch features.

)Detailed Functionality:

  1. AI will be able to generate optimal lattice structures for quantum simulators where edge effects (boundary physics) are critical, by analytically solving the corresponding non-Hermitian associated matrices.

  2. It can perform exact spectral analysis of open-boundary systems, providing closed-form expressions for energy bands and quasiparticle modes, which is far superior to traditional numerical methods that rely on truncating the system size or relying on continuous approximations (like the non-Bloch band theory).

  3. The AI will be able to characterize the stability and vacuum states of these open systems by precisely identifying whether negative single-particle energies exist, allowing it to distinguish between dynamically unstable and well-defined ground states based on system parameters.

  4. It can automatically resolve inherent mathematical degeneracies in the diagonalization process (the twofold degeneracy mentioned in Section III.A) by applying rigorous algebraic proofs derived from bosonic commutation relations, ensuring only physically constructible quasiparticle operators are generated, thereby eliminating ambiguity often found in alternative transformation methods (like the local squeezing approach).

  5. The AI can generate customized Bogoliubov transformations that explicitly map the original Hamiltonian to a diagonal form using non-unique parameters (the free real parameters derived from the solution), providing a richer set of physically equivalent descriptions for experimental realization.

Related papers