Stabilization of finite-energy grid states of a quantum harmonic oscillator by reservoir engineering with two dissipation channels

1. Problem Statement

Quantum Error Correction (QEC) using bosonic codes, specifically the Gottesman-Kitaev-Preskill (GKP) code, offers a pathway to reduce hardware requirements by encoding logical qubits into the infinite-dimensional Hilbert space of a single harmonic oscillator. While discrete-time control protocols (using ancilla qubits) have successfully generated and stabilized GKP states, they are limited by error propagation from the ancilla to the encoded qubit.

Continuous-time autonomous stabilization via reservoir engineering (coupling the system to a dissipative environment) promises significantly longer logical lifetimes. However, previous theoretical proposals (e.g., Ref. [11]) required four distinct dissipation channels and stringent hardware constraints (such as specific impedance values or high control complexity), making experimental realization difficult.

The Core Challenge: Can the protocol be simplified to use fewer dissipation channels and relaxed parameters while maintaining effective stabilization of finite-energy GKP states for both error correction and quantum metrology?

2. Methodology

The authors propose a simplified Lindblad master equation that replaces the four-dissipator dynamics of previous works with a two-dissipator model.

• The Proposed Dynamics:
  The evolution of the density matrix $\rho$ is governed by:
  $$\frac{d}{dt}\rho = \mathcal{L}(\rho) = \mathcal{D}[M_1](\rho) + \mathcal{D}[M_2](\rho)$$
  where $\mathcal{D}[M](\rho) = M\rho M^\dagger - \frac{1}{2}\{M^\dagger M, \rho\}$ is the standard Lindblad dissipator. The jump operators are:
  $$M_1 = \sin(\eta q) + i\epsilon \cos(\eta q)p$$
  $$M_2 = \sin(\eta p) - i\epsilon \cos(\eta p)q$$
  Here, $q$ and $p$ are position and momentum operators, and $\epsilon$ is a regularization parameter ensuring finite energy.

• Key Simplifications:

  1. Reduction of Channels: The model uses only $M_1$ and $M_2$, discarding the $M_3$ and $M_4$ terms found in the four-dissipator proposal.
  2. Parameter Adjustment: The lattice constant $\eta$ is halved compared to the standard square GKP lattice ($\eta = \sqrt{\pi}$ for a qubit, $\eta = \sqrt{\pi/2}$ for a "qunaught"). The authors argue intuitively that the constraints imposed by $\sin(\eta q)=0$ and $\sin(2\eta q)=0$ are effectively equivalent in the context of the stabilizer group, allowing the reduction of operators.

• Analytical Approach:

  • Energy Bounds: The authors derive a priori estimates on the photon number expectation value $\langle N \rangle$ to prove that the energy of the system remains bounded during evolution.
  • Convergence Analysis: They analyze the convergence to the codespace by studying the spectral properties of the adjoint Lindblad operator acting on periodic observables. This involves mapping the problem to a singular differential operator on periodic functions and establishing a spectral gap using weighted Poincaré and Hardy inequalities.
  • Noise Robustness: Numerical simulations are performed to study the effect of photon loss (modeled by an additional dissipator $\kappa \mathcal{D}[a]$) on logical coherence and steady-state fidelity.

3. Key Contributions

1. Simplified Protocol: Introduction of a two-dissipator Lindblad dynamics that autonomously stabilizes finite-energy GKP states, significantly reducing the experimental complexity compared to the four-dissipator standard.
2. Rigorous Energy Bounds: Proof that the solutions to the proposed Lindblad equation have bounded energy, ensuring physical realizability (finite energy states).
3. Convergence Rate Estimates: Derivation of explicit convergence rates for logical observables. The authors relate the convergence to the spectral gap of a specific differential operator, proving exponential convergence to the codespace.
4. Metrological Application: Demonstration that by adjusting parameters (specifically $\eta$), the same dynamics can stabilize a single GKP "qunaught" state (1D subspace), which is a valuable resource for quantum metrology (simultaneous measurement of modular position and momentum).
5. Implementation Feasibility: Analysis of how this reduction in dissipation channels alleviates hardware constraints in circuit QED (impedance requirements) and trapped ions (Lamb-Dicke parameter tuning).

4. Results

• Stabilization Success: Numerical simulations starting from the vacuum state show convergence to the GKP qubit state (fidelity > 92% for specific initial conditions) and the GKP qunaught state.
• Logical Decoherence under Noise:
  • When photon loss ($\kappa$) is introduced, the system remains confined to the codespace, but logical coherences decay.
  • Scaling Law: The logical decoherence rate $\Gamma$ scales as a power law: $\Gamma \propto \kappa^n / \epsilon^r$ (with fitted $n \approx 0.88, r \approx 0.57$).
  • Comparison: This scaling is qualitatively worse than the exponential suppression ($\Gamma \propto e^{-1/\sigma}$) observed in the four-dissipator model. This indicates a trade-off: the two-dissipator model is easier to implement but offers less robustness against photon loss.
• Metrology: Simulations confirm that the fine periodic structure of the GKP qunaught is preserved in the steady state even with photon loss, making it viable for metrological applications.

5. Significance

• Experimental Viability: The primary significance is the drastic reduction in experimental overhead. Implementing four distinct, high-fidelity dissipation channels is a major hurdle. Reducing this to two channels (or potentially implementing them via stroboscopic control with fewer Trotter steps) makes autonomous GKP stabilization much more accessible for current quantum hardware (superconducting circuits and trapped ions).
• Trade-off Analysis: The paper provides a clear trade-off analysis. While the four-dissipator model offers superior error suppression (exponential scaling), the two-dissipator model offers a "good enough" stabilization with significantly lower hardware complexity. This makes it a compelling candidate for near-term experimental validation and a stepping stone toward more complex codes.
• Dual Utility: The framework is versatile, capable of stabilizing both logical qubits (for QEC) and qunaughts (for metrology) simply by tuning the lattice parameter $\eta$.
• Theoretical Foundation: The mathematical derivation of energy bounds and spectral gaps for this specific class of singular differential operators provides a rigorous foundation for future analysis of autonomous bosonic codes.

In conclusion, this work proposes a pragmatic simplification of autonomous GKP stabilization. It sacrifices some theoretical robustness against noise for a substantial gain in experimental feasibility, potentially accelerating the path toward practical bosonic quantum error correction and metrology.