Quantum Noise. X. Density-Matrix Treatment of Field and Population-Difference Fluctuations

Abstract

Starting with the quantum-noise-source model of a maser developed in the fourth paper in this series and specializing to the case in which all atomic parameters except the population difference vary rapidly compared with the field, an equation is developed for the density operator $\ensuremath{\rho}(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, D, t)$ that depends on the electromagnetic field variables $b$, ${b}^{\ifmmode\dagger\else\textdagger\fi{}}$ and the population difference $D$. The corresponding equation obeyed by the density matrix ${\ensuremath{\rho}}_{\mathrm{mn}}(D, t)$ in the "$n$ representation" is obtained. The antinormal ordering correspondence $\ensuremath{\rho}(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, D, t)=\mathcal{a}{\overline{\ensuremath{\rho}}}^{(a)}(\ensuremath{\beta}, {\ensuremath{\beta}}^{*}, D, t)$ [in which $\ensuremath{\beta}$ is replaced by $b$ and ${\ensuremath{\beta}}^{*}$ by ${b}^{\ifmmode\dagger\else\textdagger\fi{}}$ with all ${b}^{\ifmmode\dagger\else\textdagger\fi{}}'\mathrm{s}$ to the right of all $b'\mathrm{s}$] defines an associated classical function ${\overline{\ensuremath{\rho}}}^{(a)}$. The equation for $\frac{\ensuremath{\partial}\ensuremath{\rho}}{\ensuremath{\partial}t}$ corresponds to a Fokker-Planck equation for $\frac{\ensuremath{\partial}{\overline{\ensuremath{\rho}}}^{(a)}}{\ensuremath{\partial}t}$. Under the Markoffian assumption, an associated classical random process is uniquely defined. The quantum-mechanical average of normal-ordered time-ordered operators is given by simple integrals over the corresponding classical distribution functions. For one-time operators this result is well known. For two-time operators, a proof is given using the quantum-regression theorem of the second paper in this series. Amplitude and phase fluctuations are discussed by rewriting the Fokker-Planck equation in the variables $I$, $\ensuremath{\varphi}$, $D$, where $\ensuremath{\beta}={I}^{\frac{1}{2}}\mathrm{exp}(\ensuremath{-}i\ensuremath{\varphi})$. The Langevin description of the classical random process associated with ${\overline{\ensuremath{\rho}}}^{(a)}(\ensuremath{\beta}, {\ensuremath{\beta}}^{*}, D, t)$ is used to make a classical adiabatic elimination of $D$. The resulting equations for ${\overline{\ensuremath{\rho}}}^{(a)}(\ensuremath{\beta}, {\ensuremath{\beta}}^{*}, t)$ and $\ensuremath{\rho}(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, t)$ are shown to agree with the corresponding results of the ninth paper in this series, obtained by a more difficult quantum adiabatic elimination of $D$. The equation for $\ensuremath{\rho}(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, t)$, rewritten in the $n$ representation agrees with corresponding results of Scully and Lamb.