We extend our previous quantum mechanical nonlinear treatment of laser noise to the following problem: We consider a set of atoms each with three levels, which support laser action of one or several modes. The laser action can take place either between the upper or the lower two levels. The atomic line is assumed to be homogeneously broadened. The broadening can be caused by the decay into the nonlasing modes, by the pumping process, lattice vibrations and other, non specified sources. The fluctuations of the atomic variables (or operators) are taken into account in a quantum mechanically consistent way using results of previous papers byHaken andWeidlich as well asSchmid andRisken. The laser modes are coupled to the thermal resonator noise usingSenitzky's method. In the first part of the present paper, we treat quite generally multimode laser action. It is shown, that each light mode chooses a specificcollective atomic “mode” to interact with. We introduce a set of suitable collective atomic “modes”, which leads to a simplification of the equations of motion for theHeisenberg operators of the light field and the atomic operators. From the new equations we can eliminate all atomic operators. We are then left with a set of coupled nonlinear, integro-differential equations for the light field operators alone. These equations, which are completely exact and valid both for running and standing waves, represent a considerable simplification of the original problem. In the second part of this paper, these equations are specialized to single mode operation, which is studied above laser threshold. In the vicinity of the threshold the laser equation can be simplified to an operator-equation, whose classical analogue is vander-Pol's equation with a noisy driving force. With increasing inversion, the full equation must be treated, however. Using the method of our previous paper, we decompose the light amplitude into a phase-factor and a real amplitude, which is expanded around its stable value. We determine the Fourier-transform of the intensity correlation function and the total intensity of the fluctuating part of the amplitude. Somewhat above threshold this intensity drops down with the inverse of the photon output power,P, while the inherent relaxation frequency increases withP. The noise intensity stems in this region from the off-diagonal elements of the noise operators and not from the diagonal elements, which are responsible for the shot noise. This result is insofar remarkable, as a rate equation treatment would include only the latter ones. Under certain conditions the intensity fluctuations can show resonances with increasing output power,P. At high inversion the vacuum fluctuations of the light field are dominant, while the other noise sources give rise to contributions which vanish with the inverse of the output power. As a by-product our treatment yields the following formula for the linewidth (half width at half power) which is caused by phase fluctuations: $$\Delta \nu = \frac{{\gamma _{3 2}^2 \kappa ^2 }}{{(\kappa + \gamma _{3 2} )^2 }}\frac{{\hbar \omega }}{P}\left( {\frac{1}{2}\frac{{(N_3 + N_2 )}}{{N_3 - N_2 }} + n_{Th} + \frac{1}{2}} \right)$$ κ is the cavity half-width, γ32 the atomic half-width for the transition (3→2). N3 and N2 are the steady state occupation numbers of the atoms, nTh is the number of thermal quanta. This formula is valid for running and standing waves (with finite wavelength) and a two, three or four level system. It contains previous results of the present author, as well as those ofLamb,Sauermann andLax, either as special cases, or generalizes the range of applicability, as discussed in the text.
Read full abstract