Reservoir-Engineered Mechanical Cat States with a Driven Qubit

TL;DR

The paper addresses deterministic generation of macroscopic quantum superpositions in a mechanical resonator by engineering the reservoir through a single driven qubit with both longitudinal and transverse couplings. It derives an effective master equation for the mechanics, featuring parity-preserving two-phonon dissipation $\mathcal{D}[a^{2}]$ and a coherent two-phonon squeezing term $\chi a^{2} + \chi^{*} a^{\dagger 2}$, enabled by driving near the two-phonon resonance $\omega_q \approx 2\omega_m$ and $\omega_d \approx \omega_q$, and adiabatic elimination of the dissipative qubit. The authors validate the mechanism with Wigner-function simulations showing the vacuum state bifurcating into a two-peak cat state with interference fringes, and they provide design rules for achieving robust parity-protected states, including the role of Kerr nonlinearity $\delta_k$ and temperature via $n_{th}$. The approach is scalable, cavity-free, and applicable to circuit QED and related platforms, offering a practical route to mechanical bosonic encodings and quantum-enhanced sensing, with possible extensions to multi-resonator networks for entangled cat states.

Abstract

Macroscopic quantum superpositions, such as mechanical Schrödinger cat states, are central to emerging quantum technologies in sensing and bosonic error-correcting codes. We propose a scheme to generate such states by coupling a nanomechanical resonator to a coherently driven two-level system via both transverse and longitudinal interactions. Driving the qubit at twice the oscillator frequency activates resonant two-phonon exchange processes, enabling coherent conversion of drive energy into phonon pairs and their dissipative stabilization. Starting from the full time-dependent Hamiltonian, we derive an effective master equation for the mechanical mode by perturbative elimination of the lossy qubit. The reduced dynamics feature engineered two-phonon loss and a coherent squeezing term, which together drive the resonator into a deterministic Schrödinger-cat state. Our approach requires only a single driven qubit and no auxiliary cavity, offering a scalable and experimentally accessible route to macroscopic quantum superpositions in circuit-QED and related platforms.

Figures (5)

• Figure 1: Time evolution of the mechanical resonator's state under the effective master equation (\ref{['eq:EffME']}), visualized via the Wigner quasi-probability distribution. The panels from left to right represent the Wigner function of the mechanical resonator at $\Gamma^{-}_{2} t = 0, 0.3, 0.6, 1$. The resonator starts in the vacuum state, appearing as a single peak at the origin in phase space at $t=0$. Subsequently, the Wigner distribution bifurcates into two lobes, indicating the formation of a superposition of two coherent-state components. By the steady state (final frame), an even Schrödinger-cat state is realized, featuring two distinct phase-space peaks and a set of interference fringes between them. Parameters: $g_z/2\pi = 6~\text{MHz}$, $g_x = 0.1 g_z$, $\omega_m/2\pi = 100~\text{MHz}$, $\kappa/2\pi = 100~\text{kHz}$, $\gamma/2\pi = 15~\text{Hz}$, $\varepsilon = 4g$. We set the qubit and drive frequencies resonantly at $\omega_q = \omega_d = 2\omega_m$ to activate the two-phonon dynamics and assume a zero-temperature environment for the mechanical bath. The mechanical Hilbert space is truncated to 60.
• Figure 2: Left: Wigner snapshots of the mechanical state from the full driven qubit–mechanics model (top row) and from the effective master equation (bottom row) at $\kappa t=0,\,1.31,\,19.66,\,39.00$. Both descriptions exhibit the bifurcation into two lobes and the development of parity-protected interference fringes, converging to a steady cat state. Right: (top) Uhlmann fidelity $F(\rho_{\mathrm{full}}(t),\rho_{\mathrm{eff}}(t))$ between the full model and the reduced effective model; (bottom) mean phonon number $\langle n\rangle(t)$ from the full (solid black) and effective (red dashed) dynamics. The effective description closely tracks the full evolution and attains $F\gtrsim 0.98$ at late times. Parameters: $g_z/2\pi=6~\mathrm{MHz}$, $g_x=0.1\,g_z$, $\omega_m/2\pi=100~\mathrm{MHz}$, $\kappa/2\pi=100~\mathrm{kHz}$, $\gamma/2\pi=15~\mathrm{Hz}$, $\varepsilon = 2g$, $\omega_q=\omega_d=2\omega_m$, and $n_{\mathrm{th}}=0$.
• Figure 3: Fidelity $\mathcal{F}$ between the reduced oscillator state obtained from the full driven qubit--oscillator master equation \ref{['eq:master']} and the oscillator state obtained from the effective oscillator-only master equation \ref{['eq:EffME']}, evaluated at a fixed comparison time $\kappa t^\ast = 20$. Each panel shows $\mathcal{F}$ as a function of a single control parameter: (a) drive amplitude $\Omega$, (b) two-phonon detuning $\Delta_{2-} = \omega_q - 2\omega_m$, (c) transverse coupling $g_x$, (d) longitudinal coupling $g_z$, (e) mechanical frequency $\omega_m$, and (f) mechanical damping rate $\gamma$. The vertical axis is plotted on a linear scale and zoomed near $\mathcal{F}=1$ to highlight small deviations between the two descriptions. In each panel, only the indicated parameter is varied while all others are held fixed. (a) $(\Delta_{2-}, g_x, g_z, \omega_m, \gamma) = (0,\ 1.0,\ 1.0,\ 120,\ 3\times 10^{-3})$, (b) $(\Omega, g_x, g_z, \omega_m, \gamma) = (0.20,\ 1.0,\ 1.0,\ 120,\ 3\times 10^{-3})$, (c) $(\Omega, \Delta_{2-}, g_z, \omega_m, \gamma) = (0.20,\ 0,\ 1.0,\ 120,\ 3\times 10^{-3})$. (d) $(\Omega, \Delta_{2-}, g_x, \omega_m, \gamma) = (0.20,\ 0,\ 1.0,\ 120,\ 3\times 10^{-3})$, (e) $(\Omega, \Delta_{2-}, g_x, g_z, \gamma) = (0.20,\ 0,\ 1.0,\ 1.0,\ 3\times 10^{-3})$, and (f) $(\Omega, \Delta_{2-}, g_x, g_z, \omega_m) = (0.20,\ 0,\ 1.0,\ 1.0,\ 120)$. All frequencies, couplings, and rates are expressed in units of $\kappa$, and the mechanical bath is assumed to be cold, $n_\mathrm{th} = 0$.
• Figure 4: Wigner distributions $W(x,p)$ of the mechanical resonator for different thermal occupations of its phonon bath. The panels, from left to right, correspond to $n_{\mathrm{th}} = 0$, $0.5$, $1$, and $2$, respectively, and show the Wigner function at a fixed evolution time $\Gamma_{2}^{-} t = 3$. For $n_{\mathrm{th}} = 0$, the steady state is a nearly ideal even Schrödinger-cat, characterized by two well-separated phase-space lobes and pronounced interference fringes. As $n_{\mathrm{th}}$ increases, thermal phonons progressively smear out the interference pattern and reduce the negativity of $W(x,p)$, illustrating the gradual degradation of cat-like nonclassicality under realistic mechanical decoherence. The system parameters are the same as in Fig. \ref{['fig:fig1']}, except $\varepsilon = 10 g$. The evolution is governed by the effective master equation (\ref{['eq:EffME']}) and the mechanical Hilbert space truncated at $N=60$.
• Figure 5: Steady-state mean phonon number $\bar{n}_{\mathrm{ss}}$ as a function of the longitudinal coupling $g_z/2\pi$. The blue solid line with circles shows the result obtained from the effective master equation \ref{['eq:EffME']}, while the orange dashed line with squares corresponds to the cat-state approximation given in Eq. \ref{['eq:A12']}. Parameters: $g_x/g_z = 0.1$, $\omega_m/2\pi = 100~\mathrm{MHz}$, $\omega_q = 2\omega_m$, $\kappa/2\pi = 100~\mathrm{kHz}$, $\gamma/2\pi = 15~\mathrm{Hz}$, $n_{\rm th} = 5$, $\Omega/2\pi = 0.10~\mathrm{MHz}$, and $N = 80$.