Stabilization of cat-state manifolds using nonlinear reservoir engineering

1. Problem Statement

Quantum error correction (QEC) for bosonic systems (e.g., superconducting cavities, trapped ions) often relies on reservoir engineering, where dissipation is engineered to stabilize specific quantum states (manifolds) that encode logical information.

• Limitations of Current Methods: Standard approaches typically operate in the low-excitation regime (e.g., the Lamb-Dicke regime for ions or weak nonlinearity for circuits). In this regime, the stabilizing jump operators $\hat{K}$ are restricted to low-order boson processes (e.g., $\hat{K} = \hat{a}^d - \alpha^d$).
• The Challenge: Higher-order boson processes (required for robust codes like $d$-legged cat codes with $d>2$) are exponentially suppressed in the low-excitation regime, making them experimentally inaccessible. Furthermore, existing analytical frameworks struggle to describe systems where the coupling strength depends non-trivially on the photon number (strong coupling/nonlinear regimes).
• The Gap: There is a lack of a general, intuitive theoretical framework to design and analyze dissipative stabilization schemes that utilize strong nonlinearities to create complex, rotationally symmetric bosonic codes with enhanced error correction capabilities.

2. Methodology: Nonlinear Reservoir Engineering (NLRE)

The authors introduce Nonlinear Reservoir Engineering (NLRE), a paradigm shift from standard analytic treatments.

• Core Principle: The method relies on destructive interference between competing boson gain and loss processes.
  • The system is coupled to a zero-temperature auxiliary bath via a Hamiltonian $\hat{H} = \Omega(\hat{K}\hat{A}^\dagger + \text{h.c.})$, where $\hat{K}$ is a nonlinear operator acting on the target oscillator.
  • $\hat{K}$ is defined as a combination of a boson raising process of order $r$ and a lowering process of order $l$:
    $$\hat{K} = \hat{a}^{\dagger r} f(\hat{n}) - g(\hat{n})\hat{a}^l$$
  • The scalar functions $\tilde{f}(k)$ and $\tilde{g}(k)$ (Rabi frequencies) represent the strength of these processes for Fock state $|k\rangle$.
• Mechanism of Stabilization:
  • The system is driven toward higher energies where gain dominates ($\tilde{f} > \tilde{g}$) and toward lower energies where loss dominates ($\tilde{g} > \tilde{f}$).
  • Stabilization occurs at the crossing point $k^*$ where $\tilde{f}(k^*) \approx \tilde{g}(k^* + r)$.
  • At this crossing, the gain and loss terms destructively interfere, creating a "dark subspace" (manifold) where the population accumulates.
• Analytical Framework:
  • By linearizing $\tilde{f}$ and $\tilde{g}$ around the crossing point $k^*$, the authors derive that the steady-state boson distribution follows a Conway–Maxwell–Binomial (CMB) distribution (or its limits, Conway–Maxwell–Poisson).
  • This allows for the analytical calculation of the mean boson number, variance, and the effective confinement rate ($\kappa_{conf}$), which dictates how quickly the system returns to the manifold after a perturbation.

3. Key Contributions

A. Generalization of Cat Codes

The authors demonstrate that a $d$-dimensional cat manifold (where $d = r + l$) can be stabilized by various combinations of $(r, l)$, not just the standard $(0, d)$ or $(d, 0)$.

• New Degrees of Freedom: Schemes like $(1, 1)$, $(2, 2)$, or $(1, 3)$ can stabilize manifolds of dimension $d=2, 4, 8$, etc.
• Relaxed Experimental Requirements: This removes the need to engineer high-order processes of order $d$ directly. Instead, lower-order processes (e.g., 1st and 2nd order) can be combined to stabilize higher-dimensional codes.

B. Error Correction Capabilities

The paper analyzes the error correction properties based on the symmetry of the $(r, l)$ scheme:

• Dephasing Protection: Schemes with $(r, l) = (0, d)$ or $(d, 0)$ possess strong symmetry and autonomously correct dephasing errors (standard cat code behavior).
• Momentum/Quadrature Protection: The novel $(1, 1)$ scheme lacks strong symmetry for dephasing but exhibits autonomous correction against momentum errors (errors in the $\hat{p}$ quadrature). This is crucial for stabilizing squeezed cat states.
• Noise Bias Enhancement: By engineering the slopes of $\tilde{f}$ and $\tilde{g}$, the authors show that the Mandel Q parameter (variance) of the steady state can be reduced (number-squeezing). This significantly increases the confinement rate ($\kappa_{conf}$), leading to a massive enhancement in the noise bias (ratio of phase-flip to bit-flip rates) without increasing the mean photon number.

C. Generalization via Unitary Transformations

The framework is extended to General NLRE using isomorphic transformations (Bogoliubov transformations).

• By applying a unitary transformation $\hat{U}$ (e.g., squeezing), the standard NLRE model can be mapped to stabilize quadrature-squeezed cat states.
• This explains recent proposals for autonomous correction of squeezed cat codes and suggests that a minimum of four boson processes (two 1st-order and two 3rd-order) is required to autonomously correct both dephasing and loss errors simultaneously.

4. Experimental Results and Implementations

The paper proposes and simulates two physical implementations:

1. Trapped Ions (Beyond Lamb-Dicke Regime):

   • Mechanism: Utilizes the full nonlinear light-matter interaction of a trapped ion driven by bichromatic laser fields. The Rabi frequencies are approximated by Bessel functions ($J_r$).
   • Result: Simulations for a $^9$Be$^+$ ion in a Penning trap show that a 4-legged cat manifold ($(r,l)=(0,4)$) can be stabilized with fidelities $>90\%$ even in the presence of motional heating, dephasing, and photon recoil.
   • Significance: This validates the use of strong coupling (outside the Lamb-Dicke regime) as a resource rather than a limitation.

2. Superconducting Circuit QED:

   • Mechanism: Uses a high-Q LC oscillator coupled to a DC-biased Asymmetrically Threaded SQUID (ATS). The DC bias and flux drives allow for the engineering of specific nonlinear jump operators.
   • Result: The authors demonstrate that this setup can realize the toy model Hamiltonian, stabilizing 4-legged cat states with a confinement rate of $\sim 10$ MHz.
   • Significance: This provides a versatile route to realizing rotation-symmetric codes in superconducting circuits without the spurious terms (like AC Stark shifts) often found in standard Josephson junction implementations.

5. Significance and Impact

• New Design Paradigm: NLRE provides a systematic, intuitive "design recipe" for creating bosonic codes by engineering the crossing of gain and loss rates, moving away from restrictive low-order expansions.
• Enhanced Robustness: The ability to tune the boson distribution's variance (via the slopes of the crossing) allows for the optimization of error correction rates and noise bias, potentially reducing the overhead for fault-tolerant quantum computing.
• Access to New States: The method enables the stabilization of non-standard states, such as squeezed cat states and high-dimensional rotationally symmetric codes, which were previously difficult to access or analyze.
• Universality: The framework unifies existing dissipative stabilization schemes (like standard cat codes and squeezed vacuum stabilization) under a single theoretical umbrella and extends them to new regimes applicable to ions, circuits, and potentially optomechanical systems.

In conclusion, this work establishes a robust theoretical and experimental foundation for utilizing strong nonlinearities in quantum systems to engineer highly protected, complex quantum states, offering a promising pathway toward scalable, fault-tolerant bosonic quantum computing.