Generating the negative Wigner functions where the corresponding Wigner states are nonclassical has been recognized as a powerful tool for successfully performing quantum information and computing protocols beyond the scope of classical computers. Here, we present the possibility to generate and engineer the negative Wigner function at a steady state using parametric (two-photon) driving and homodyne-based feedback in a quantum van der Pol (vdP) oscillator. Specifically, we employ a quantum master equation approach for calculating the Wigner function of the vdP oscillator field in phase space and, furthermore, quantifying its negativity content. We clearly show that the negative-value magnitudes, regions, and shapes of the Wigner function can be effectively tuned by the parametric driving phase and the parametric driving amplitude, as well as the feedback coefficient within a large range. We identify different contributions of these involved parameters to the Wigner-function negativity. In the present scheme, more complex quantum coherence and interference phenomena are introduced via the parametric driving and feedback, which stabilizes the phase of the vdP oscillator field and renders the capability to generate the negative Wigner function. Therefore, the enhanced Wigner-function negativity can be achieved under these optimized system parameters. Our in-depth study provides insight into the formation and in situ control of the desirable Wigner nonclassical states. The obtained results are not limited to the vdP oscillator systems and should be generally applicable to other coherent coupled systems within the reach of modern experimental facilities.
Generating the negative Wigner functions where the corresponding Wigner states are nonclassical has been recognized as a powerful tool for successfully performing quantum information and computing protocols beyond the scope of classical computers. Here, we present the possibility to generate and engineer the negative Wigner function at a steady state using parametric (two-photon) driving and homodyne-based feedback in a quantum van der Pol (vdP) oscillator. Specifically, we employ a quantum master equation approach for calculating the Wigner function of the vdP oscillator field in phase space and, furthermore, quantifying its negativity content. We clearly show that the negative-value magnitudes, regions, and shapes of the Wigner function can be effectively tuned by the parametric driving phase and the parametric driving amplitude, as well as the feedback coefficient within a large range. We identify different contributions of these involved parameters to the Wigner-function negativity. In the present scheme, more complex quantum coherence and interference phenomena are introduced via the parametric driving and feedback, which stabilizes the phase of the vdP oscillator field and renders the capability to generate the negative Wigner function. Therefore, the enhanced Wigner-function negativity can be achieved under these optimized system parameters. Our in-depth study provides insight into the formation and in situ control of the desirable Wigner nonclassical states. The obtained results are not limited to the vdP oscillator systems and should be generally applicable to other coherent coupled systems within the reach of modern experimental facilities.
We propose a method of producing steady-state coherent light with negative Wigner functions in nonlinear media combined with feedback control. While the nonlinearities are essential to produce the Wigner negativities, this alone is insufficient to stabilize steady-state light with negativities. Using feedback control to control the phase in the cavity, we find that this produces significant total negativities for reasonable experimental parameters. The negative Wigner function is produced continuously and does not appear to be restricted to low-amplitude light. The technique is applicable to systems such as exciton-polaritons, where strong natural nonlinearities are present.
Scalable methods for detecting and quantifying the nonclassical nature of a quantum state in noisy environments are challenging due to a complex relationship between noise and quantum coherence. In particular, identifying experimentally accessible signatures of synchronization in such regimes remains an open problem. By leveraging promising experimental implementation, we underpin what possible direct measures of nonclassicality are available. This work outlines accessing quantum synchronization (QS) in the steady state of a driven quantum van der Pol oscillator (vdPo) using two distinct figures of merit: (i) the nonclassical area δ and (ii) the second-order correlation function g(2)(0), both of which are viable in experimental architectures. The nonclassical area quantifier based on homodyne tomography allows us to assess the nonclassical nature of the vdPo state directly from the tomogram without requiring full state reconstruction or Wigner function negativity. Within a well-defined parameter regime of drive strength and detuning, both δ and g(2)(0) exhibit pronounced signatures of synchronization that complements the phase coherence between the drive and the vdPo. We derive an analytical expression for the steady state density matrix and the corresponding tomogram of the system, valid for arbitrary strengths of the harmonic drive. Analysis of the quantum tomogram uncovers clear phase locking behaviour, enabling the identification of the synchronization region (Arnold tongue) directly in terms of g(2)(0) and δ. Furthermore, the behaviour of g(2)(0) provides a statistical perspective that reinforces the tomographic signatures of QS. By analyzing the interplay between the aforementioned metrics, our findings indicate a scalable and experimentally relevant framework for characterizing QS in the driven vdPo.
In this work we demonstrate numerically that the nonlinearity provided by a continuously driven two-level system allows for the generation of Wigner-negative states of the electromagnetic field confined in one spatial dimension. Wigner-negative states, also known as Wigner nonclassical states, are desirable for quantum information protocols beyond the scope of classical computers. Focusing on the steady-state emission from the two-level system, we find the largest negativity at the drive strength where the coherent reflection vanishes.
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.
Abstract In this paper, we propose a measurement nonlinear feedback control scheme to generate Wigner-function negativity in an optical cavity having dynamics described as a linear quantum system. In general, linear optical quantum systems can be easily constructed with reliable devices; therefore, the idea of constructing the entire system with such an optical system and nonlinear feedback is reasonable for generating Wigner-function negativity. However, existing studies have insufficiently examined the realizability or actual implementation of feedback control, which essentially requires fast responses from the sensors and actuators. In order to solve this problem, we consider the realizable feedback control of the optical phase of a pumping beam supplied to a cavity by using electro-optical modulation, which can be utilized as a fast control actuator. Then, we introduce mathematical models of the feedback-controlled system and evaluate its effect on the generation of the Wigner-function negativity by using numerical simulation. Through various numerical simulations, we show that the proposed feedback control can effectively generate the negativity of the Wigner function.
We study quantum van der Pol oscillator in which single-photon pumping competes with a nonlinear two-photon dissipative channel to produce a quantum limit cycle, and the question of how two such oscillators with non-identical intrinsic frequencies synchronize in the sense of phase locking when coupled dissipatively to one another. For a single oscillator, starting from a Lindblad description of the dynamics, we derive the exact Wigner-function evolution equation, and show that its Fokker--Planck approximation admits in the appropriate regime an equivalent classical Itô stochastic description that reproduces remarkably well the full Lindblad dynamics. Notably, the stochastic approach provides an efficient means of capturing system dynamics in the semiclassical regime, where a large number of energy eigenstates contribute and direct numerical simulations become computationally demanding. Restricting to the three lowest Fock states and using a parity symmetry, we obtain in closed form the stationary-state density matrix and its Wigner-function representation that displays pronounced nonclassical features. For two non-identical oscillators, we study the full Lindblad dynamics numerically and characterize synchronization through the stationary-state order parameter as a function of frequency mismatch and dissipative coupling. We find that increasing the effective nonlinear dissipative strength suppresses the synchronized region. Restricting the two oscillators to the three-level qutrit subspace, a parity symmetry reduces the stationary-state problem to an exact 41-dimensional linear system, allowing to compute the synchronization order parameter. In this reduced description, results in the limit of large dissipation agree closely with those from the full dynamics, while exhibiting in the limit of small dissipation a double-peak synchronization structure accompanied by a phase reversal.