Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review

Quantum Theory of a Laser Model

C. R. Willis

Published Jul 8, 1966• 22 citations
Publisher

Abstract

We derive the kinetic equations for the coupled single-particle density matrix $\ensuremath{\rho}$ and the electromagnetic density matrix $R$ to lowest order in the dimensionless coupling constant ${\ensuremath{\beta}}^{2}\ensuremath{\equiv}{(\frac{{\ensuremath{\omega}}_{L}}{{\ensuremath{\omega}}_{D}})}^{2}$. The laser frequency ${\ensuremath{\omega}}_{L}$ is ${(4\ensuremath{\pi})}^{\ensuremath{-}\frac{1}{2}}{(\mathfrak{N}{r}_{0}{\ensuremath{\lambda}}^{2})}^{\frac{1}{2}}{\ensuremath{\omega}}_{0}$ where $\mathfrak{N}$ is the number of two-level systems per unit volume, ${r}_{0}$ is the classical electron radius, $\ensuremath{\lambda}$ is the wavelength of the radiation, and $k{\ensuremath{\omega}}_{0}$ is the two-level energy difference. The Doppler frequency ${\ensuremath{\omega}}_{D}$ characterizes the center-of-mass motion. For gas lasers ${\ensuremath{\beta}}^{2}$ is much less than 1 and, consequently, we generalize and use the Bogoliubov derivation of kinetic equations for weak interactions. We find solutions when the average field vanishes and which include spontaneous emission correctly. The single-particle density matrix and the radiation density matrix are coupled through their second moments. When we substitute the solution of the second-moment equations into the density-matrix equations, we find that each density matrix satisfies an uncoupled linear equation with known time-dependent coefficients. We introduce and discuss dissipation from the density-matrix point of view. With the use of the density-matrix formalism we indicate that the correct expansion parameter for higher order kinetic equations is ${\ensuremath{\beta}}^{2}$.

DOI
10.1103/physrev.147.406

Related papers

Quantum Theory of a Gas Laser

C. R. Willis

We derive the equation of motion for the quantum-mechanical radiation density matrix of a gas laser to lowest order in the dimensionless coupling constant. Our derivation is fully quantum mechanical and we can calculate the coherence properties from the radiation density matrix. Our model consists of $N$ two-level systems interacting with radiation in a cavity in the presence of dissipation, pumping, and collisions. The method we use is a generalization of the Bogoliubov derivation of the kinetic equation for a small parameter. Our derivation holds for any physically realizable pump power, and near threshold reduces to Lamb's near-threshold theory. With the equation of motion for the radiation density matrix, we obtain solutions both for when the average field is nonzero and for when it is zero. The steady-state electromagnetic density is the same in both cases except for a small spontaneous-emission term. We show that the reason gas lasers do not satisfy rate equations is the existence of zeroth-order correlations between the internal atomic variables and atomic center-of-mass variables. It is these same zeroth-order correlations which are responsible for the Lamb dip. Our derivation includes collisions and reduces the calculations of their effect to quadrature.

Physical Review•Apr 10, 1967•14 citations

Quantum Theory of a Gas Laser. II

C. R. Willis

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.

Physical Review•Jan 10, 1968•14 citations

Mean-Field Kinetic Equations for a Laser.

Richard H. Picard

A set of kinetic equations is derived for a single-mode laser based on the assumption that the mean optical field of the laser model need not vanish. They constitute a coupled system for the field and single-atom density operators and are obtained from the quantum-Liouville equation by a generalization of the method of Bogolyubov. The equations are derived to second order in the matter-field coupling constant, assuming an asymptotic state with no matter-field correlations. (Author)

•Jan 1, 1968•0 citations

A Model of Interacting Radiation and Matter. II

C. R. Willis

We extend the investigation of the long time behavior of a model consisting of N two-level atoms interacting with radiation to include dissipation, pumping, and center-of-mass motion. We show that when the effective expansion parameter γ̃Nγ̃ is small, the self-consistent field approximation remains the solution to lowest order in γ̃Nγ̃. The inclusion of the center-of-mass motion introduces another dimensionless parameter β into the theory. We show that when this parameter is small, the electromagnetic field amplitude varies slowly on the time scale of the center-of-mass motion. We solve the equations of motion in this slowly varying limit by an extension of the Bogoliuboff-Kryloff theory of quasilinearity to the problem of time-dependent integral kernals. We find the unique stable stationary state and show that in the slowly varying limit the stationary state is approached independent of initial conditions. We calculate the frequency shift to second order in β. The first-order frequency shift is the same as that calculated by Lamb. We compare our steady-state solution with recent experiments with lasers. We include the effect of collisions in the steady state.

Journal of Mathematical Physics•Dec 1, 1965•10 citations

Master-Equation Approach to Spontaneous Emission. III. Many-Body Aspects of Emission from Two-Level Atoms and the Effect of Inhomogeneous Broadening

G. S. Agarwal

A general theory of spontaneous emission was developed in Papers I and II of this series. In the present paper, the mastre equation describing spontaneous emission from $N$ identical two-level atoms is investigated further. It is shown that the equation of motion for the $p(<N)$-atom density matrix is coupled to the ($p+1$)-atom density matrix. We first consider the initial excitation of the system to state $|{\ensuremath{\theta}}_{0},{\ensuremath{\phi}}_{0}〉$ (${\ensuremath{\theta}}_{0}<\ensuremath{\pi}$) and solve the hierarchy of equations so obtained in the "uncorrelated approximation." It is found that the presence of other atoms is equivalent to an external field which is determined from a self-consistent analysis. It is also found that in this approximation the state of the atom at time $t$ can be described by a single parameter $\ensuremath{\theta}(t)$ ($\ensuremath{\theta}<\ensuremath{\pi}$), with ${sin}^{2}\frac{1}{2}\ensuremath{\theta}$ giving the probability that the atom is to be found in the excited state. Next, some improvements over the uncorrelated approximation are given. This is done by decoupling the equation of motion for higher-order mean values. It is found that the fluctuation $〈({S}_{i}^{Z}\ensuremath{-}〈{S}_{i}^{Z}〉)({S}_{j}^{Z}\ensuremath{-}〈{S}_{j}^{Z}〉)〉$ for any pair of atoms is of the order of $\frac{1}{N}$. The case when ${\ensuremath{\theta}}_{0}=\ensuremath{\pi}$ is considered separately and we obtain an expression for the radiation rate by making a "Hartree-Fock" type of approximation on the two-particle mean values. In this case, the behavior of the radiation rate is found to differ markedly from the "sech" behavior. It is then shown that the entire dynamics of two-level atoms emitting spontaneously can be described by a set of $2N$ coupled first-order equations which clearly exhibit the type of nonlinearity (which is the analog of the van der Pol type of nonlinearity) for this problem and provide a better understanding of spontaneous emission. Finally in Sec. V, the theory is extended to include the effects of inhomogeneous broadening, and the functional dependence of the radiation rate on the atomic-line-shape factor is obtained.

Physical Review A•Nov 1, 1971•62 citations

Quantum Theory of an Optical Maser. VI. Transient Behavior

Y. K. Wang, Willis E. Lamb

The transient behavior of a laser is discussed using the quantum theory as did Scully and Lamb. The formal solution of the density-matrix equation is expressed in terms of exponentially decaying eigenmodes. Some of the lower decay constants are obtained numerically. The equations for the moments of the density matrix are then derived and solved by a truncation method. The equations of motion are integrated numerically for the case where the average number of photons in a laser cavity has the realistically large value 1.3 \ifmmode\times\else\texttimes\fi{} ${10}^{5}$. An alternative Fokker-Planck-equation approach is discussed.

Physical Review A•Aug 1, 1973•21 citations
We derive the kinetic equations for the coupled single-particle density matrix \ensuremathρ\ensuremath{\rho}\ensuremathρ and the electromagnetic density matrix RRR to lowest order in the dimensionless coupling constant \ensuremathβ2\ensuremath≡(\ensuremathωL\ensuremathωD)2{\ensuremath{\beta}}^{2}\ensuremath{\equiv}{(\frac{{\ensuremath{\omega}}_{L}}{{\ensuremath{\omega}}_{D}})}^{2}\ensuremathβ2\ensuremath≡(\ensuremathωD​\ensuremathωL​​)2. The laser frequency \ensuremathωL{\ensuremath{\omega}}_{L}\ensuremathωL​ is (4\ensuremathπ)\ensuremath−12(Nr0\ensuremathλ2)12\ensuremathω0{(4\ensuremath{\pi})}^{\ensuremath{-}\frac{1}{2}}{(\mathfrak{N}{r}_{0}{\ensuremath{\lambda}}^{2})}^{\frac{1}{2}}{\ensuremath{\omega}}_{0}(4\ensuremathπ)\ensuremath−21​(Nr0​\ensuremathλ2)21​\ensuremathω0​ where N\mathfrak{N}N is the number of two-level systems per unit volume, r0{r}_{0}r0​ is the classical electron radius, \ensuremathλ\ensuremath{\lambda}\ensuremathλ is the wavelength of the radiation, and k\ensuremathω0k{\ensuremath{\omega}}_{0}k\ensuremathω0​ is the two-level energy difference. The Doppler frequency \ensuremathωD{\ensuremath{\omega}}_{D}\ensuremathωD​ characterizes the center-of-mass motion. For gas lasers \ensuremathβ2{\ensuremath{\beta}}^{2}\ensuremathβ2 is much less than 1 and, consequently, we generalize and use the Bogoliubov derivation of kinetic equations for weak interactions. We find solutions when the average field vanishes and which include spontaneous emission correctly. The single-particle density matrix and the radiation density matrix are coupled through their second moments. When we substitute the solution of the second-moment equations into the density-matrix equations, we find that each density matrix satisfies an uncoupled linear equation with known time-dependent coefficients. We introduce and discuss dissipation from the density-matrix point of view. With the use of the density-matrix formalism we indicate that the correct expansion parameter for higher order kinetic equations is \ensuremathβ2{\ensuremath{\beta}}^{2}\ensuremathβ2.
We derive the equation of motion for the quantum-mechanical radiation density matrix of a gas laser to lowest order in the dimensionless coupling constant. Our derivation is fully quantum mechanical and we can calculate the coherence properties from the radiation density matrix. Our model consists of NNN two-level systems interacting with radiation in a cavity in the presence of dissipation, pumping, and collisions. The method we use is a generalization of the Bogoliubov derivation of the kinetic equation for a small parameter. Our derivation holds for any physically realizable pump power, and near threshold reduces to Lamb's near-threshold theory. With the equation of motion for the radiation density matrix, we obtain solutions both for when the average field is nonzero and for when it is zero. The steady-state electromagnetic density is the same in both cases except for a small spontaneous-emission term. We show that the reason gas lasers do not satisfy rate equations is the existence of zeroth-order correlations between the internal atomic variables and atomic center-of-mass variables. It is these same zeroth-order correlations which are responsible for the Lamb dip. Our derivation includes collisions and reduces the calculations of their effect to quadrature.
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 RRR 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 RRR 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 RRR to answer time-dependent questions about quantities, such as the line shape, that depend on phase.
We extend the investigation of the long time behavior of a model consisting of N two-level atoms interacting with radiation to include dissipation, pumping, and center-of-mass motion. We show that when the effective expansion parameter γ̃Nγ̃ is small, the self-consistent field approximation remains the solution to lowest order in γ̃Nγ̃. The inclusion of the center-of-mass motion introduces another dimensionless parameter β into the theory. We show that when this parameter is small, the electromagnetic field amplitude varies slowly on the time scale of the center-of-mass motion. We solve the equations of motion in this slowly varying limit by an extension of the Bogoliuboff-Kryloff theory of quasilinearity to the problem of time-dependent integral kernals. We find the unique stable stationary state and show that in the slowly varying limit the stationary state is approached independent of initial conditions. We calculate the frequency shift to second order in β. The first-order frequency shift is the same as that calculated by Lamb. We compare our steady-state solution with recent experiments with lasers. We include the effect of collisions in the steady state.
A general theory of spontaneous emission was developed in Papers I and II of this series. In the present paper, the mastre equation describing spontaneous emission from NNN identical two-level atoms is investigated further. It is shown that the equation of motion for the p(<N)p(<N)p(<N)-atom density matrix is coupled to the (p+1p+1p+1)-atom density matrix. We first consider the initial excitation of the system to state $|{\ensuremath{\theta}}_{0},{\ensuremath{\phi}}_{0}〉$ (${\ensuremath{\theta}}_{0}<\ensuremath{\pi}$) and solve the hierarchy of equations so obtained in the "uncorrelated approximation." It is found that the presence of other atoms is equivalent to an external field which is determined from a self-consistent analysis. It is also found that in this approximation the state of the atom at time ttt can be described by a single parameter $\ensuremath{\theta}(t)$ ($\ensuremath{\theta}<\ensuremath{\pi}$), with ${sin}^{2}\frac{1}{2}\ensuremath{\theta}$ giving the probability that the atom is to be found in the excited state. Next, some improvements over the uncorrelated approximation are given. This is done by decoupling the equation of motion for higher-order mean values. It is found that the fluctuation $〈({S}_{i}^{Z}\ensuremath{-}〈{S}_{i}^{Z}〉)({S}_{j}^{Z}\ensuremath{-}〈{S}_{j}^{Z}〉)〉$ for any pair of atoms is of the order of 1N\frac{1}{N}N1​. The case when ${\ensuremath{\theta}}_{0}=\ensuremath{\pi}$ is considered separately and we obtain an expression for the radiation rate by making a "Hartree-Fock" type of approximation on the two-particle mean values. In this case, the behavior of the radiation rate is found to differ markedly from the "sech" behavior. It is then shown that the entire dynamics of two-level atoms emitting spontaneously can be described by a set of 2N2N2N coupled first-order equations which clearly exhibit the type of nonlinearity (which is the analog of the van der Pol type of nonlinearity) for this problem and provide a better understanding of spontaneous emission. Finally in Sec. V, the theory is extended to include the effects of inhomogeneous broadening, and the functional dependence of the radiation rate on the atomic-line-shape factor is obtained.
The transient behavior of a laser is discussed using the quantum theory as did Scully and Lamb. The formal solution of the density-matrix equation is expressed in terms of exponentially decaying eigenmodes. Some of the lower decay constants are obtained numerically. The equations for the moments of the density matrix are then derived and solved by a truncation method. The equations of motion are integrated numerically for the case where the average number of photons in a laser cavity has the realistically large value 1.3 \ifmmode\times\else\texttimes\fi{} 105{10}^{5}105. An alternative Fokker-Planck-equation approach is discussed.