We arrive at transition probabilities for a laserlike system by neglecting wave properties of matter and radiation, and treating the problem, in an analogy to chemical systems, entirely in a number representation. Thus, for example, stimulated emission is idealized as a collision between a photon and an excited atom and the transition rate is proportional to the amount of each species present in the cavity. In all, transition probabilities for stimulated emission, absorption of radiation, cavity loss, atomic deexcitation (other than stimulated emission), and pumping of atoms are included in the model. The resulting stochastic master equation analogous to that employed for the discussion of fluctuations in chemical systems is then used to obtain the deterministic behavior of atoms and field. This agrees with that of the quantum mechanical laser model of Casagrande and Lugiato; thus the rate parameters of the stochastic model may be related to the parameters of their model. Returning to the master equation, we project out the atomic variable and use adiabatic elimination to arrive at a closed, single-variable equation for P(n), the probability of having n photons in the cavity. The stationary solution of this equation is obtained and graphs of (n) and log[(n2)/(n)2] versus the pumping parameter show close agreement with the microscopic results of the Casa-grande-Lugiato model. A similar discussion of the fluctuations can also be given for a laser with saturable absorption. It shows a first-order phase transition. A comparison with the microscopic results of Lugiato et al. is outlined. Again a remarkable agreement of the stochastic approach via the master equation and the microscopic theory is obtained.
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 fully quantum-mechanical model of a single-mode laser is examined. A previously introduced method of solution where both the pump and loss reservoirs are treated simultaneously as a weak perturbation of the atom-light system is employed. This method is similar to perturbation theory for degenerate systems in that the reservoirs remove the nonuniqueness (degeneracy) of the steady state of the unperturbed atom-light system. In the limit of weak reservoir coupling, the exact steady-state solution of the master equation is derived. This solution still contains the atomic variables explicitly and therefore remains valid even in situations where the conventional adiabatic elimination of the atomic variables is no longer possible. A compact discrete representation of the master equation in terms of atomic occupation numbers and photon numbers is developed which allows one to perform all calculations in an elementary manner, requiring only appropriate summations over photon and occupation numbers. General exact relations are obtained for the steady-state probability distribution function
p(n,m) to find
n photons in the light field and
m atoms in the excited state. Our solution constitutes the first steady-state laser distribution function that is valid for an arbitrary number of laser atoms and in the complete region from far below to far above threshold. It provides a quantitative criterion for the range of validity of approximate methods in laser theory.
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.
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 analogy of a laser with an autocatalytic chemical reaction is used to write down a macroscopic stochastic master equation of the birth-and-death type for the photon probability function in a semiclassical laser formalism. This equation may be solved exactly in the steady state using a generating function. This leads to a prescription for calculating all the higher-order moments of the photon number as well as providing a unique truncation of the hierarchy of moment equations at any given order. Further, the variance is given by a Poisson distribution, in agreement with the strong-signal Scully-Lamb theory. This master equation is also shown to lead to an exact Fokker-Planck equation by using a modified Kramers-Moyal expansion. A comparison between our results and the (microscopic) Scully-Lamb theory shows that the two approaches give the same results at large photon number.
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.