We develop a systematic projection-operator technique for constructing Gaussian approximations and their perturbative corrections in bosonic nonlinear models. As a case study, we apply it to the driven dissipative Kerr oscillator. In the absence of external driving, the model can be solved exactly within a low-dimensional Fock subspace, leading to strongly non-Gaussian states. Nevertheless, we demonstrate that the evolution of first- and second-order moments is captured by our Gaussian scheme with high accuracy even in this regime, providing a natural benchmark. For the general case with external driving, our approach reduces the equations of motion to a closed system for means and covariances and allows one to compute systematic corrections beyond the Gaussian level in closed form. We also calculate the dynamics of linear and quadratic combinations of creation and annihilation operators in both weak- and strong-drive regimes.
We develop a systematic projection-operator technique for constructing Gaussian approximations and their perturbative corrections in bosonic nonlinear models. As a case study, we apply it to the driven dissipative Kerr oscillator. In the absence of external driving, the model can be solved exactly within a low-dimensional Fock subspace, leading to strongly non-Gaussian states. Nevertheless, we demonstrate that the evolution of first- and second-order moments is captured by our Gaussian scheme with high accuracy even in this regime, providing a natural benchmark. For the general case with external driving, our approach reduces the equations of motion to a closed system for means and covariances and allows one to compute systematic corrections beyond the Gaussian level in closed form. We also calculate the dynamics of linear and quadratic combinations of creation and annihilation operators in both weak- and strong-drive regimes.
Abstract This paper studies a class of single-mode bosonic models of open quantum systems in which the dynamics of the creation and annihilation operators remains linear despite the strongly non-Gaussian nature of the induced dissipation. It is shown that, in contrast to isolated bosonic systems—where linear dynamics uniquely corresponds to quadratic Hamiltonians and Gaussian states—linearity of the operator dynamics in open systems does not guarantee Gaussianity of the state evolution. We construct quantum channels and Gorini–Kossakowski–Sudarshan–Lindblad generators for which the dynamics of the creation and annihilation operators is linear and closed, while Gaussian states are driven beyond the convex hull of Gaussian states during the evolution. In particular, the emergence of negativity of the Wigner function at a finite evolution time is demonstrated.
Many recent advancements in quantum computing leverage strong drives on nonlinear systems for state preparation, signal amplification, or gate operation. However, the interplay within such strongly driven system introduces multi-scale dynamics that affects the long-time behavior of the system in non-trivial ways that are very difficult to model. Therefore, the analysis of these systems often relies on effective Hamiltonian models that introduce additional nonlinear processes which approximate the long-time dynamics so that highly oscillatory terms may be ignored. However, the removal of such high frequency transitions can only be performed rigorously within a systematic framework of time-coarse graining, which is a fundamentally irreversible operation. This implies that standard approaches with unitary effective models cannot accurately capture the long-time behavior of strongly driven nonlinear quantum systems in general. We introduce a systematic perturbation theory for obtaining the complete non-unitary effective model of the time-coarse grained (TCG) dynamics of a driven quantum system to any order in the coupling strengths. We derive a closed-form analytical formula for both unitary and non-unitary contributions, in the form of an effective Hamiltonian and non-unitary (pseudo-)dissipators. Remarkably, even though the effective theory presumes unitary time evolution at the microscopic level, the time-coarse grained dynamics is found to follow a non-unitary time evolution in general. This occurs even when there is no open heat reservoir for the system to become entangled with or dissipate into. We demonstrate the effectiveness of the new method using several typical models of driven nonlinear systems in superconducting circuits, and show that it generalizes and improves on existing methods by providing more accurate results and explaining phenomena that have not been accounted for.
We construct Feynman rules for Lindblad master equations in the operator language of quantum optics: a fixed algebra of elementary operator moves and closed-form propagators for perturbation theory about an exactly solvable, in most cases Gaussian, generator. The expansion separates into a kinematic part (the shift algebra on the Wick-ordered eigenoperators Omn=:(c†)mcn:, from which the vertex table for a specified interaction is assembled) and a dynamical part that enters only through the resolvent of the unperturbed generator, diagonal and additive for any stable, diagonalizable Gaussian L0, including thermal and squeezed reservoirs and correlated decay, and equally available for a two-level emitter, driven or undriven. In the driven Kerr cavity the loop expansion is the semiclassical expansion and the series is asymptotic. In a strongly driven two-level atom, dressed lines replace a series of finite radius by one accurate deep into saturation, one line per Mollow component. For Kerr rings and chains the cost is polynomial in the number of sites at fixed perturbative order. On a three-site ring the diagrams recover the non-Gaussian part of the three-photon correlation that the Gaussian cumulant closure lacks; carried to tenth order they give the smallest g(3) error of the methods compared for U≤0.25κ at a cost below that of a fourth-order closure, and third- and fourth-order closures are more accurate at moderate coupling. The same rules, extended to the counting-field-tilted generator, give the photon-counting cumulants of a disordered eight-site chain beyond direct diagonalization; quantum-jump trajectories at three couplings agree with them to within 1.6 standard errors, within one at the strongest coupling, where they place the third cumulant six standard errors above the Gaussian value.
Driven nonlinear quantum oscillators are a central platform for quantum technologies, yet their dissipative dynamics are typically described using Lindblad or Caldeira-Leggett master equations derived under assumptions that exclude nonlinearities and driving. Here, we derive a generalized Caldeira-Leggett master equation for driven nonlinear oscillators by retaining the full nonlinear and time-dependent system dynamics in the construction of the dissipator. For position- and momentum-dependent system-bath coupling, the dissipator itself becomes dynamically dressed, generating nonlinear and drive-dependent dissipative channels beyond conventional fixed-dissipator approaches. This produces nonlinear damping together with dissipation-induced corrections to the effective drive. The resulting dissipative dynamics suppress large-amplitude excitations and reduce phase-space fluctuations. For a driven Kerr oscillator, this leads to the suppression of bistability, asymmetric resonance responses, and strongly modified fluctuation distributions. More broadly, our results establish a microscopic framework in which nonlinear dynamics and driving directly reshape the dissipative sector of driven open quantum systems.
Driven-dissipative systems have recently attracted great attention due to the existence of novel physical phenomena with no analog in the equilibrium case. The Keldysh path-integral theory is a powerful tool to investigate these systems. However, it has still been a challenge to study strong nonlinear effects implemented by recent experiments, since in this case the photon number is few and quantum fluctuations play a crucial role in the dynamics of the system. Here we develop an approach for deriving exact steady states of driven-dissipative systems by introducing the Keldysh partition function in the Fock-state basis and then mapping the standard saddle-point equations into Keldysh-Heisenberg equations. We take the strong Kerr nonlinear resonators with and without the nonlinear driving as two examples to illustrate our method. It is found that, in the absence of the nonlinear driving, the exact steady state obtained does not exhibit bistability and agrees well with the complex P-representation solution. While in the presence of the nonlinear driving, the multiphoton resonance effects are revealed and are consistent with the qualitative analysis. Our method provides an intuitive way to explore a variety of driven-dissipative systems especially with strong correlations.
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 [21]. 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 [46]. 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 [30]. 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.