We have solved a master equation to obtain the diagonal elements of the density matrix for laser light. The master equation represents a generalization of a master equation derived earlier for a single radiation mode interacting with
N stationary two-level atoms. The generalization takes into account the pumping scheme which characterizes a three-level laser. For positive temperatures and a lossless cavity, the equilibrium solution of the master equation gives the correct statistical-mechanical description of the atoms and radiation. For negative temperatures, the distribution function for the number of photons in the mode is either exponential or peaked, depending on whether the laser is operating below or above threshold. The calculated intensity fluctuations are in good agreement with semiclassical results for lasers operating slightly above threshold.
The steady-state intensity distribution for a laser operating somewhat above threshold is determined by solving a master equation. The master equation describes the interaction between a single radiation mode and
N two-level atoms, it contains the effects of pumping and radiation loss from the cavity, and it is applicable to both homogeneous and Doppler broadening. Forms of the master equation applicable to both three- and four-level lasers are considered. Solutions are obtained under the assumption that the distribution of atoms and radiation intensity is a bivariate Gaussian. Fluctuations in the output of giant-pulse lasers are also considered.
The interaction is considered between
N stationary two-level atoms and a radiation field described by a single cavity mode. The state vector for the complete system of atoms plus radiation is expressed as a linear superposition of states constructed from a product of photon states in the
n representation and products of Pauli spin eigenstates describing all combinations of atoms in the lower and upper energy levels. Equations for the corresponding probability amplitudes are derived by substituting this superposition into the Schr\"odinger equation. The resulting equations are combined into bilinear form and phenomenological damping contributions are added. After the neglect of certain of the bilinear quantities, a master equation is derived which governs the probability
Pmn of having
m atoms in the lower level and
n photons in the mode. This master equation takes account of multiple single-quantum absorption and emission processes but not of simultaneous multiple processes involving two or more atoms at a time. The equation which governs the expected number of photons
〈n〉 derived from the master equation bears a close resemblance to a rate equation. The effect of radiation loss from the cavity is incorporated into the master equation. Numerical calculations for a
Q-spoiled laser show that the statistics of the number of photons in the mode bear a qualitative resemblance to Poisson statistics.
The quantum theory of the three-energy-level two-mode laser in which all the levels are being pumped is investigated and the master equation of the laser operation is obtained. It contains two variables,
n1, the photon number of the first mode, and
n2, the photon number of the second mode. There appear four terms in the master equation which represent two-photon processes between the two modes, the absorption of a photon from one mode, and emission of a photon to the other. They are absent in the single-mode case and in the two-mode case with pumping only to the upper level. The master equation can be represented by a diagram of probability in two dimensions which can be extended to infinity, and each arrow of it represents a term on the right-hand side of the master equation. By summing each variable (
n1 or
n2) of the master equation, two equations can be obtained. Each equation contains only one variable, and can be expressed by a diagram of probability flow in one dimension. By considering the correspondence between macroscopic equilibrium and microscopic detailed balance the equations of motion in the steady state are obtained. By introducing a parameter
H, the equation of motion is simplified and a formal solution can be deduced which is dependent on
H. It appears that the parameter
H cannot be fully determined, but some properties of it can be deduced which prove that the introduction of
H is reasonable. The operation characteristics of the laser under different conditions are discussed on the basis of the formal solution: among them are the threshold condition, the condition for one-mode operation, and the change of the photon-statistical distribution. We especially discuss the laser output power curves under two-mode operation above threshold. The curve of one mode, which oscillates first, shows a bending-down phenomenon as the excitation increases when the other mode goes beyond its threshold as demonstrated experimentally by Otsuka. This phenomenon is explained qualitatively and is attributed to the two-photon processes between the two modes.
We generalize the method of deriving a kinetic equation for a gas laser developed in the first paper of this series to include radiation-matter correlations. As a consequence, we are able to show that the first Born approximation with asymptotic conditions which contain radiation-matter correlations is sufficient to explain the nonthermal photon distribution observed in photon-counting experiments. Our derivation includes the spatial and velocity effects of the motion of the gas atoms. After adiabatically eliminating the matter variables, we obtain a master equation for the radiation density matrix
R with nonlinear coefficients. Our radiation master equation for stationary atoms without spatial effects is the same as the generalized Fokker-Planck equation of the Langevin-noise-equation treatment of lasers. We show that the generalized Fokker-Planck equation is the result of the first-Born-approximation treatment of the j\ifmmode\cdot\else\textperiodcentered\fi{}A interaction and is not dependent upon the microscopic structure of the reservoirs. In the first Born approximation, the diagonal and off-diagonal matrix elements of
R do not interact with each other. Consequently, an initially diagonal radiation density matrix remains diagonal. However, we show that even though the stationary state is diagonal it is necessary to know the propagator of the off-diagonal part of
R to answer time-dependent questions about quantities, such as the line shape, that depend on phase.
The quantum theory of two-level atom single-mode laser system of Sargent, Scully and Lamb is extended to two-mode laser in a three-level atomic system with common lower level, and master equation is obtained for the system with all the three levels being pumped. There appear four terms in the master equation that are absent in the single-mode case and those interpreted as 2-photon processes play an especial role in the present case. The master equation is represented by a probability flow diagram in two-dimensions with the photon numbers n1 and n2 of the two-mode as variables. The probability flow diagram can be extended to infinity, and each arrow of it represents a term of the right hand side of master equation. By taking sum over one of the variables n1 or n2, the master equation is reduced to two equations, each of them can be represented by an one-dimensional probability flow diagram. The equation of motion under steady state is obtained by considering the correspondence between the macroscopic balance and the microscopic balance.A parameter H is introduced to reduce the equation of motion and to facilitate the mathematical procedure and a formal solution is thus obtained. Using the equation of motion and its formal solution various operation characteristics are discussed. Among them are: the threshold condition for each mode; the condition of only one mode operation; the variation of photon statistics between single-mode operation and two-mode operation and the variation of photon statistics with lower level being pumped and not pumped.