Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Beyond the Fokker-Planck equation: Stochastic Simulation of Complete Wigner representation for the Optical Parametric Oscillator

Lev I. Plimak, M. K. Olsen, M. Fleischhauer, Matthew J. Collett

Published Jan 1, 2001• 23 citations
OpenAlex

Abstract

We demonstrate a method which allows the stochastic modelling of quantum systems for which the generalised Fokker-Planck equation in the phase space contains derivatives of higher than second order. This generalises quantum stochastics far beyond the quantum-optical paradigm of three and four-wave mixing problems to which these techniques have so far only been applicable. To verify our method, we model a full Wigner representation for the optical parametric oscillator, a system where the correct results are well known and can be obtained by other methods. 1 I.

Related papers

Beyond the Fokker-Planck equation: Stochastic simulation of complete Wigner representation for the optical parametric oscillator

L. I Plimak, M. K. Olsen, M Fleischhauer +1 more

We demonstrate a method which allows the stochastic modelling of quantum systems for which the generalised Fokker-Planck equation in the phase space contains derivatives of higher than second order. This generalises quantum stochastics far beyond the quantum-optical paradigm of three- and four-wave mixing problems to which these techniques have so far only been applicable. To verify our method, we model a full Wigner representation for the optical parametric oscillator, a system where the correct results are well known and can be obtained by other methods.

arXiv:quant-ph/0108128•Europhysics Letters (EPL)•Nov 1, 2001•36 citations

Numerical solution methods for quantum stochastic processes

Eli Pollak

The study of quantum stochastic processes presents severe difficulties, both on the theory level as well as on technical grounds. The numerically exact solution remains prohibitive even today. In this paper we review and present new results for three different methods used for the modelling of quantum stochastic processes. These include a mixed quantum classical approach, semiclassical initial value representations of the quantum propagator and the reduced density matrix approach as typified by the quantum Wigner-Fokker-Planck equation. A new semiclassical initial value representation that does away with cumbersome prefactors which depend on the monodromy matrix elements but is exact for a harmonic oscillator is presented and its properties analysed. A recently proposed systematic method for improving semiclassical initial value representations is reviewed. The generalization of the Wigner-Fokker-Planck equation to stochastic processes with memory is obtained by using a novel integral equation representation.

Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIE•May 8, 2003•1 citation

Langevin equation for the squeezing of light by means of a parametric oscillator

Trevor W. Marshall, Emilio Santos

We show that the Langevin equation for a nonlinear-optical system may be obtained directly from the Heisenberg equation of motion for the annihilation operators, provided a certain linearization procedure is valid. We apply the technique to the parametric oscillator used to generate squeezed light and compare our results to those obtained from Fokker-Planck-type equations. We argue that, only when the Wigner, as opposed to the P or Q, representation of quantum optics is used, do we get a correct description of the underlying stochastic process. We show how the linearization procedure may be carried out to describe the operation of the parametric oscillator both below threshold, where a squeezed vacuum state results, and above threshold, where we find a squeezed coherent state. In the region of the threshold a heuristic extension of the method leads to a possible description of the system by means of a nonlinear Langevin equation.

Physical Review A•Feb 1, 1990•29 citations

Quantum-theoretical treatments of three-photon processes

M. K. Olsen, Lev I. Plimak, Michael Fleischhauer

We perform and compare different analyses of triply degenerate four-wave mixing in the regime where three fields of the same frequency interact via a nonlinear medium with a field at three times the frequency. As the generalized Fokker-Planck equation (GFPE) for the positive-$P$ function of this system contains third-order derivatives, there is no mapping onto genuine stochastic differential equations. Using techniques of quantum field theory, we are able to write stochastic difference equations that we may integrate numerically. We compare the results of this method with those obtained by the use of approximations based on semiclassical equations, and on truncation of the GFPE leading to stochastic differential equations. In the region where the difference equations converge, the stochastic methods agree for the field intensities, but give different predictions for the quantum statistics.

Physical Review A•Apr 22, 2002•12 citations

Quantum theory of nonlinear fiber optics: Phase-space representations

S. J. Carter

In this paper the equations for quantum optical pulse propagation in nonlinear and dispersive single-mode fibers are presented in terms of two phase-space formulations based on the positive-P and the Wigner distributions. Included are the effects due to the coupling of the electromagnetic modes to the vibrational states of a vitreous silica fiber. By making use of the well-known equivalence of Fokker-Planck and Ito stochastic equations, we demonstrate two alternative methods for formulating the equations of motion as coupled stochastic c-number equations for the propagating field. The first method involves a representation of the density operator in terms of the positive P distribution function. This leads to exact stochastic equations of motion. The second method makes use of the Wigner distribution function. This method, which requires truncation of third-order derivative terms in the corresponding Fokker-Planck equation, is necessarily approximate. However, we discuss certain advantages to the Wigner approach that have made it the preferable method for exploratory work.

Physical Review A•Apr 1, 1995•50 citations

Simulation approach to Fokker-Planck modeling of lasers and optical amplifiers

Giacomo Mauro D’Ariano, Chiara Macchiavello

A novel numerical method to solve Fokker-Planck equations in quantum optics is presented, based on a Monte Carlo simulation of the probability diffusion process. The method is especially useful for multimode analysis and hence for studying realistic models of nonlinear optical systems. Two simple examples are given: the first is a 1D Fokker-Planck equation in the number representation, which describes a simple model of optical amplifier; the second is a 2D equation in the P-function representation-the quasi-probability for normal-ordered averages-which corresponds to the customary Van der Pol model of the laser threshold. In this case also the field correlation function and spectrum are numerically simulated.

Proceedings of SPIE, the International Society for Optical Engineering/Proceedings of SPIE•May 3, 1994•0 citations
We demonstrate a method which allows the stochastic modelling of quantum systems for which the generalised Fokker-Planck equation in the phase space contains derivatives of higher than second order. This generalises quantum stochastics far beyond the quantum-optical paradigm of three and four-wave mixing problems to which these techniques have so far only been applicable. To verify our method, we model a full Wigner representation for the optical parametric oscillator, a system where the correct results are well known and can be obtained by other methods. 1 I.
We demonstrate a method which allows the stochastic modelling of quantum systems for which the generalised Fokker-Planck equation in the phase space contains derivatives of higher than second order. This generalises quantum stochastics far beyond the quantum-optical paradigm of three- and four-wave mixing problems to which these techniques have so far only been applicable. To verify our method, we model a full Wigner representation for the optical parametric oscillator, a system where the correct results are well known and can be obtained by other methods.
The study of quantum stochastic processes presents severe difficulties, both on the theory level as well as on technical grounds. The numerically exact solution remains prohibitive even today. In this paper we review and present new results for three different methods used for the modelling of quantum stochastic processes. These include a mixed quantum classical approach, semiclassical initial value representations of the quantum propagator and the reduced density matrix approach as typified by the quantum Wigner-Fokker-Planck equation. A new semiclassical initial value representation that does away with cumbersome prefactors which depend on the monodromy matrix elements but is exact for a harmonic oscillator is presented and its properties analysed. A recently proposed systematic method for improving semiclassical initial value representations is reviewed. The generalization of the Wigner-Fokker-Planck equation to stochastic processes with memory is obtained by using a novel integral equation representation.
We show that the Langevin equation for a nonlinear-optical system may be obtained directly from the Heisenberg equation of motion for the annihilation operators, provided a certain linearization procedure is valid. We apply the technique to the parametric oscillator used to generate squeezed light and compare our results to those obtained from Fokker-Planck-type equations. We argue that, only when the Wigner, as opposed to the P or Q, representation of quantum optics is used, do we get a correct description of the underlying stochastic process. We show how the linearization procedure may be carried out to describe the operation of the parametric oscillator both below threshold, where a squeezed vacuum state results, and above threshold, where we find a squeezed coherent state. In the region of the threshold a heuristic extension of the method leads to a possible description of the system by means of a nonlinear Langevin equation.
We perform and compare different analyses of triply degenerate four-wave mixing in the regime where three fields of the same frequency interact via a nonlinear medium with a field at three times the frequency. As the generalized Fokker-Planck equation (GFPE) for the positive-PPP function of this system contains third-order derivatives, there is no mapping onto genuine stochastic differential equations. Using techniques of quantum field theory, we are able to write stochastic difference equations that we may integrate numerically. We compare the results of this method with those obtained by the use of approximations based on semiclassical equations, and on truncation of the GFPE leading to stochastic differential equations. In the region where the difference equations converge, the stochastic methods agree for the field intensities, but give different predictions for the quantum statistics.
In this paper the equations for quantum optical pulse propagation in nonlinear and dispersive single-mode fibers are presented in terms of two phase-space formulations based on the positive-P and the Wigner distributions. Included are the effects due to the coupling of the electromagnetic modes to the vibrational states of a vitreous silica fiber. By making use of the well-known equivalence of Fokker-Planck and Ito stochastic equations, we demonstrate two alternative methods for formulating the equations of motion as coupled stochastic c-number equations for the propagating field. The first method involves a representation of the density operator in terms of the positive P distribution function. This leads to exact stochastic equations of motion. The second method makes use of the Wigner distribution function. This method, which requires truncation of third-order derivative terms in the corresponding Fokker-Planck equation, is necessarily approximate. However, we discuss certain advantages to the Wigner approach that have made it the preferable method for exploratory work.
A novel numerical method to solve Fokker-Planck equations in quantum optics is presented, based on a Monte Carlo simulation of the probability diffusion process. The method is especially useful for multimode analysis and hence for studying realistic models of nonlinear optical systems. Two simple examples are given: the first is a 1D Fokker-Planck equation in the number representation, which describes a simple model of optical amplifier; the second is a 2D equation in the P-function representation-the quasi-probability for normal-ordered averages-which corresponds to the customary Van der Pol model of the laser threshold. In this case also the field correlation function and spectrum are numerically simulated.