A fully quantum-mechanical treatment is given of field and population fluctuations in a homogeneously broadened laser. The field and population reservoirs have been replaced by dissipation coefficients (or transition probabilities) and noncommuting noise sources. The atomic polarization variables are adiabatically eliminated. The mean motion of an arbitrary operator function $M(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, {N}_{1}, {N}_{2})$ of the field variables
b, ${b}^{\ifmmode\dagger\else\textdagger\fi{}}$ and the population variables
N1 and
N2 is calculated. The equation of motion for the density matrix $\ensuremath{\rho}(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, {N}_{1}, {N}_{2}, t)$ is deduced, as well as the Fokker-Planck equation for the associated classical distribution function $P(\ensuremath{\beta}, {\ensuremath{\beta}}^{*}, {N}_{1}, {N}_{2}, t)$. The diffusion coefficients in the latter equation demonstrate the shot-noise origin of the fluctuations. The Fokker-Planck equation for
P is replaced by Langevin equations for $\ensuremath{\beta}$, ${\ensuremath{\beta}}^{*}$,
N1, and
N2 to facilitate adiabatic elimination of the population variables (when valid). The resulting Langevin equations for $\ensuremath{\beta}$, ${\ensuremath{\beta}}^{*}$ are converted to a Fokker-Planck equation for $P(\ensuremath{\beta}, {\ensuremath{\beta}}^{*}, t)$, to an operator equation for $\ensuremath{\rho}(b, {b}^{\ifmmode\dagger\else\textdagger\fi{}}, t)$, and to an equation for ${\ensuremath{\rho}}_{\mathrm{mn}}$, the field density matrix in the photon-number representation. The transformation $\ensuremath{\beta}={I}^{\frac{1}{2}}{e}^{\ensuremath{-}i\ensuremath{\phi}}$ to an ($I, \ensuremath{\phi}$) representation leads (a) to the steady-state solution
P(I) valid at all operating levels, and (b) to the phase linewidth. The transformation to new variables
N=21(m+n), $k=m\ensuremath{-}n$, permits the equation for ${\ensuremath{\rho}}_{\mathrm{mn}}$ to be replaced by a Fokker-Planck equation for $\ensuremath{\rho}(N, k, t)$. The latter permits a verification of the previously obtained phase linewidth and a solution for the steady-state $\ensuremath{\rho}(N)$.
The equation of motion for the slowly varying field amplitude of the laser containing an intracavity saturable absorber (LSA) is where a is the net gain (pump parameter) and b is a saturation parameter of the laser. In our experiment the absorber saturates faster than the amplifier and the parameter b is positive. q(t) is a complex Gaussian white noise random process representing spontaneous emission noise. From Eq. (1) we obtain for the steady-state probability density of light intensity I = | E |2. From this equation various moments of light intensity can be calculated and compared with experimentally measured quantities. For b > 0 the LSA exhibits bistability and a first order phase transition occurs at a = -3b2/16. For negative values of b there is a continuous transition at a = 0. Experiments were carried out on a He—Ne laser with an intracavity Ne absorber cell. The laser was operated close to the center of the absorption line. Slow drifts of light intensity were controlled by an electronic feedback loop. The parameter b was varied by changing the current in the absorber tube. Pump parameter a was then varied by changing passive losses in the cavity by means of an adjustable knife edge. The main beam from the LSA was allowed to fall on a detector/am- plifier combination. The voltage output from the detector was digitized by a multichannel analyzer and a histogram of voltage amplitudes was built up. Since the voltage output from the detector is directly proportional to the light intensity incident on the detector, voltage histogram is a direct representative of the intensity probability distribution. From the measured histogram moments of light intensity were calculated. These results are compared with the theoretical predictions.
The optical bistability (OB) behavior of single-mode Fabry- Perot lasers containing saturable absorber (LSA) is presented on the basis of the rate equations with allowance for spontaneous emission and spatial hole-burning. Special attention is paid to the case of dominant inhomogeneous broadening in both Lorentzian and Gaussian models. The OB occurrence conditions are clearly displayed in phase diagrams. The influence of the main LSA parameters on OB characteristics as well as on OB curves' stability are investigated in detail. Some practical suggestions for optimizing the bistable operation in LSA are also given.