Precise Control of Quantum Cat and Fock States in Driven-Dissipative Multi-Mode Kerr Cavities via Engineered Nonlinearity Ratios

Driven-dissipative Kerr cavities serve as a versatile platform for generating nonclassical states, such as quantum cat states and high-photon-number Fock states, which are essential for quantum information processing. Building on the quantized Lugiato-Lefever equation (QLLE), which predicts steady-state coherent state superpositions in multi-mode Kerr resonators, and the framework for engineering Kerr nonlinearity ratios to mitigate spectral crowding in coupled oscillators, we propose an integrated approach for precise quantum control. By engineering the Kerr ratio K1/K2 to approximate an incommensurate value using complex rational approximations, we eliminate parasitic degeneracies and enable selective addressing of transitions in the QLLE steady state. Employing a Magnus expansion, we derive an effective Hamiltonian that incorporates Stark-shift corrections, facilitating deterministic synthesis of entangled cat states and Fock states with fidelities exceeding 99.9%. Numerical simulations using QuTiP validate robustness against dissipation and thermal noise, demonstrating quantum hysteresis and Wigner negativity in multi-mode configurations. This synthesis bridges dissipative phase transitions with architectural control, offering a blueprint for scalable bosonic quantum processors in circuit QED. The Wigner function plays a pivotal role in characterizing these nonclassical states, providing a phase-space representation that reveals quantum interference effects. In our framework, the emergence of negative values in the Wigner function for the central mode confirms the non-Gaussian nature of the cat states, as dictated by Hudson’s theorem, which states that only Gaussian pure states have non-negative Wigner functions everywhere 23. This negativity arises from the interference term in the cat-state Wigner expression, quantifying macroscopic quantum coherence and serving as a witness for entanglement in multi-mode systems. For instance, our simulations at F = 2.0 yield a minimum Wigner value of approximately -0.006, highlighting the system’s departure from classical behavior and its utility for bosonic error correction, where cat codes leverage phase-space separation to suppress bit-flip errors exponentially with cat size. Furthermore, the relation to Wigner functions extends to dynamical properties: the steady-state Wigner negativity is robust against thermal noise due to the dissipative stabilization mechanism, aligning with results from driven Bose-Hubbard models where two-mode cat states maintain negativity under loss 50. In the QLLE context, the multi-mode extension introduces spatial localization, where the Wigner function’s negative regions correspond to soliton-like structures, as evidenced by photon population distributions 0.33, 0.90, 0.33. This negativity is not merely a feature but a consequence of fundamental limits on classical simulability, relating to the Gottesman-Knill theorem’s extension to continuous variables, where non-Gaussian operations are necessary for universal quantum computation. Quantum hysteresis, observed in our hysteresis plots, is intimately linked to the Wigner representation through the bifurcation dynamics. Near the codimension-2 point, the system’s bistability manifests as distinct paths in phase space, with Wigner functions switching between Gaussian-like (low-photon) and cat-like (high-photon) forms, supported by spectral theory of Liouvillians that guarantees metastable states in dissipative phase transitions 32. The forward and backward sweeps in our simulations demonstrate this path dependence, with the Wigner negativity peaking in the high-drive regime, underscoring the role of Kerr nonlinearity in engineering quantum resources. Overall, the Wigner function’s negativity serves as a diagnostic tool for the efficacy of our engineered ratios, directly tying to experimental verifiability in circuit QED setups, where homodyne measurements can reconstruct Wigner distributions to confirm cat-state fidelity. This work not only advances control protocols but also highlights the profound interplay between phase-space quasi-probabilities and dissipative quantum engineering, paving the way for fault-tolerant bosonic qubits.