POSITIVE-P REPRESENTATION FOR 2ND-HARMONIC GENERATION - ANALYTIC AND COMPUTATIONAL RESULTS

The problem of the formulation and simulation of quantized optical systems using the positive-P representation is considered. We show that feasible computational strategies are determined mainly by the eigenstructure of the system of stochastic differential equations that correspond to the positive-P Fokker-Planck equation. We demonstrate analytically that the stability and accuracy of numerical solutions are generally much improved by adopting time-centered implicit schemes. These findings are then related to the specific problem of second-harmonic generation. It is shown that the eigenstructure of the stochastic system possesses interrelationships that can lead to severe difficulties in numerical simulations. In particular we show how the problem of spiking encountered at high levels of quantum noise can be related to the numerical stiffness of the problem. We also demonstrate the advantages of time-centered implicit methods by highlighting the distortions that can occur in the Euler-type simulation of limit cycle solutions. Of physical interest is the fact that the phase-space morphology of the quantized second-harmonic generator admits families of limit cycles. This is in marked contrast to the purely classical picture in which, for a given parameter set, there exists only a single limit cycle. Finally we demonstrate by direct analytic solution that, under the conditions of adiabatic elimination, the positive-P representation for the second-harmonic generator with zero driving field fails, in the sense that it produces nonphysical results. This failure is in all respects analogous to the breakdown encountered in the problem of nonlinear optical damping. © 1990 The American Physical Society.