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

arXiv:quant-ph/0108128 · Europhysics Letters (EPL)

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

L. I Plimak, M. K. Olsen, M Fleischhauer, M. J Collett

Published Nov 1, 2001• 36 citations
arXivHTMLPDF

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.

DOI
10.1209/epl/i2001-00529-8

Related papers

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

Lev 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. 1 I.

•Jan 1, 2001•23 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

Semiclassical Wigner distribution for a two-mode entangled state generated by an optical parametric oscillator

K. Dechoum, Marcelo Dumas Hahn, Raúl O. Vallejos +1 more

We derive the steady-state solution of the Fokker-Planck equation that describes the dynamics of the nondegenerate optical parametric oscillator in the truncated Wigner representation of the density operator. We assume that the pump mode is strongly damped, which permits its adiabatic elimination. When the elimination is correctly executed, the resulting stochastic equations contain multiplicative noise terms and do not admit a potential solution. However, we develop a heuristic scheme leading to a satisfactory steady-state solution. This provides a clear view of the intracavity two-mode entangled state valid in all operating regimes of the optical parametric oscillator. A non-Gaussian distribution is obtained for the above threshold solution.

arXiv:1107.2667•Physical Review A•Apr 27, 2010•12 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.
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.
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.
We derive the steady-state solution of the Fokker-Planck equation that describes the dynamics of the nondegenerate optical parametric oscillator in the truncated Wigner representation of the density operator. We assume that the pump mode is strongly damped, which permits its adiabatic elimination. When the elimination is correctly executed, the resulting stochastic equations contain multiplicative noise terms and do not admit a potential solution. However, we develop a heuristic scheme leading to a satisfactory steady-state solution. This provides a clear view of the intracavity two-mode entangled state valid in all operating regimes of the optical parametric oscillator. A non-Gaussian distribution is obtained for the above threshold solution.
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.