Stabilization of cat-state manifolds using nonlinear reservoir engineering

We introduce a novel reservoir engineering approach for stabilizing multi-component Schrödinger’s cat manifolds. The fundamental principle of the method lies in the destructive interference at crossings of gain and loss Hamiltonian terms in the coupling of an oscillator to a zero-temperature auxiliary system, which are nonlinear with respect to the oscillator’s energy. The nature of these gain and loss terms is found to determine the rotational symmetry, energy distributions, and degeneracy of the resulting stabilized manifolds. Considering these systems as bosonic error-correction codes, we analyze their properties with respect to a variety of errors, including both autonomous and passive error correction, where we find that our formalism gives straightforward insights into the nature of the correction. We give example implementations using the anharmonic laser-ion coupling of a trapped ion outside the Lamb-Dicke regime as well as nonlinear superconducting circuits. Beyond the dissipative stabilization of standard cat manifolds and novel rotation symmetric codes, we demonstrate that our formalism allows for the stabilization of bosonic codes linked to cat states through unitary transformations, such as quadrature-squeezed cats. Our work establishes a design approach for creating and utilizing codes using nonlinearity, providing access to novel quantum states and processes across a range of physical systems.

💡 Research Summary

The paper introduces a new paradigm for dissipative quantum state engineering called Nonlinear Reservoir Engineering (NLRE). Traditional reservoir‑engineering approaches rely on low‑order bosonic processes (e.g., single‑photon loss, two‑photon loss) to create dark states that encode information in bosonic modes. While effective for simple cat codes or GKP states, these methods become impractical when one wishes to stabilize higher‑dimensional rotation‑symmetric manifolds because the required nonlinear processes scale unfavourably with the code order.

NLRE overcomes this limitation by coupling a harmonic oscillator to a zero‑temperature auxiliary system through a Hamiltonian of the form
 H = Ω(K A† + h.c.),
where the engineered operator K contains a raising term a†ʳ multiplied by a function f( n̂ ) and a lowering term aˡ multiplied by g( n̂ ):
 K = a†ʳ f( n̂ ) − g( n̂ ) aˡ.
The matrix elements ˜f(k)=⟨k+r|a†ʳ f( n̂ )|k⟩ and ˜g(k)=⟨k|g( n̂ ) aˡ|k+l⟩ are called “Rabi frequencies”. Their ratio G(k)=|˜f(k)|/|˜g(k+r)| determines the direction of energy flow: G>1 drives the oscillator to higher photon numbers, G<1 drives it down. By engineering the functional forms of f and g such that G crosses unity at a chosen photon number k* (i.e., |˜f(k*)|=|˜g(k*+r)|≡h*), the system becomes dynamically attracted to the crossing point. Because the raising and lowering processes are part of the same Lindblad jump operator, they interfere destructively at the crossing, creating a dark subspace of pure states localized near k*.

The dimension of this stabilized manifold is d = r + l. In phase space the steady states exhibit d‑fold rotational symmetry, forming a set of “nonlinear coherent states” spaced by 2π/d. This framework therefore generalizes the familiar two‑component cat (r=0,l=2) and the d‑legged cat (K = a^d − α^d) to a whole family of codes where the same manifold can be obtained from many different (r,l) pairs. Crucially, the required physical processes need not be of order d; for example, a four‑legged cat (d=4) can be stabilized with (r,l) = (1,3) or (2,2), which involve only third‑order or second‑order nonlinearities, dramatically easing experimental constraints.

The authors develop an analytical approximation by linearizing ˜f and ˜g around k*. This yields closed‑form expressions for the mean photon number ⟨n̂⟩, its variance, and the effective confinement rate (the linewidth of the engineered dissipation). These quantities directly determine the bias between phase‑flip and loss errors and thus the error‑correction performance of the code.

Error‑correction analysis shows that the commutation relation between K and the rotation operator R = e^{2πi n̂/d} dictates which error channels are autonomously corrected. When (r,l) = (0,d) or (d,0) the code autonomously suppresses dephasing (phase‑flip) errors, reproducing the well‑known cat‑code bias. When (r,l) = (1,1) the code protects a specific quadrature (either q or p) because the stabilized states lie along that axis; this property underlies the robustness of quadrature‑squeezed cat states against photon‑loss and gain. By reducing the photon‑number variance through appropriate choice of the slopes of ˜f and ˜g at the crossing, the authors demonstrate that the noise bias can be enhanced by several orders of magnitude compared with conventional cat codes.

Implementation proposals are given for two leading quantum‑hardware platforms.

1. Trapped ions – By driving motional sideband transitions beyond the Lamb‑Dicke regime, the effective Rabi frequencies become Bessel‑function‑shaped functions of the phonon number. Adjusting laser detuning, intensity, and phase allows precise control over the crossing point k* and the curvature of the functions, thereby realizing arbitrary (r,l) pairs. This opens the possibility of stabilizing rotation‑symmetric codes that have not yet been demonstrated with ions, and highlights the richness of ion‑laser physics outside the usual low‑excitation approximation.

2. Superconducting circuit QED – The authors reinterpret recent two‑photon cat‑stabilization experiments in terms of NLRE, and then propose a concrete circuit consisting of a microwave cavity coupled capacitively to an asymmetrically threaded SQUID (ATS) biased with a DC voltage. The nonlinear inductance of the ATS yields photon‑number‑dependent coupling rates that realize the desired f( n̂ ) and g( n̂ ) functions. Because the same nonlinearities appear in the ion implementation, the circuit can generate the same family of manifolds, including d > 2 rotation‑symmetric codes, without requiring high‑order Josephson processes.

Finally, the paper extends the basic model to include multiple raising and lowering terms, showing that a Bogoliubov transformation of the ladder operators leads to yet richer manifolds and error‑correction properties. This demonstrates that NLRE is not limited to a single K operator but can be generalized to a whole class of engineered dissipators.

In summary, the work provides a systematic design methodology for dissipatively stabilizing multi‑component cat manifolds using nonlinear gain‑loss interference. By treating the crossing of photon‑number‑dependent rates as a tunable resource, NLRE enables the creation of high‑dimensional, rotation‑symmetric bosonic codes with flexible error‑bias properties, using experimentally realistic nonlinearities in both trapped‑ion and superconducting‑circuit platforms. This advances the toolbox for autonomous quantum error correction and paves the way toward more scalable bosonic quantum memories and processors.