A Rudimentary Resevoir Problem^{[2]}

Consider a complete system which consists of a single simple harmonic oscillator (the system component of "interest") coupled to aof many simple harmonic oscillators.reservoirIn particular, the complete interacting system has a simple Hamiltonian of the form

[ VII-1 ] where

ais the system variable ofinterest(the "state of the system") and the 's are theothersystem variables (the "reservoir states"). System operators commute with reservoir operators at a given time, we may easily generate the following set of equations of motion:

[ VII-2a ]

[ VII-2b ] Formally integrating Equation [ VII-2b ], we obtain [3]

[ VII-3 ] The first term in this equation describes the

evolution of the reservoir states in the absence of any interaction with the system and the second gives the modification of this free evolution as a consequence of the coupling to the system. Inserting Equation [ VII-3 ] into Equation [ VII-2a ] we obtainfree

[ VII-4 ] Here the second summation describes a source of fluctuations arising from the free evolution of the reservoir states and the first is a feedback through the reservoir which might be described as a

To remove the rapid variation in the system variable we transform to the Heisenberg interaction picture --radiation reaction.i.e.we take [4]

[ VII-5 ] so that

[ VII-6 ] where

[ VII-7 ] is the so called

noise operator.The first term on the RHS of Equation [ VII-6 ] may be reordered as

. [ VII-8a ] By assumption, the frequencies of the reservoir oscillators are closely spaced and widely distribution so that the function is very sharply peaked at and has a width of the order of the inverse of the reservoir's frequency bandwidth. Thus, may taken out of the integration and the limits of the integration may be extended to infinite --

i.e.

. [ VII-8b ] Again following arguments explicated by Heitler,[5] we see that

[ VII-9 ] If we neglect the imaginary part of Equation [ VII-9 ], which leads to a Lamb-type frequency shift, Equation [ VII-6 ] leads directly to a quantum mechanical version of the classical Langevin equation [6] --

viz.

[ VII-10 ] where

. [ VII-11 ] is the quantum mechanical

Langevindrift coefficient.In optical applications the averageintensityornumber operator -- i.e.-- is the experimental quantity of greatest interest. A useful equation of motion for that operator may obtained by making use the Langevin equation --i.e.Equation [ VII-10 ].To that end we first write

[ VII-12 ] To obtain and we make use of the identity

[ VII-13 ] so that

[ VII-14 ] The first term on the RHS vanishes

sinceThat is, cannot be correlated with a future value of the fluctuating force. Similarly, in the second term vanishes except in the infinitesimal small inteval where .causecannot precedeeffect.

Therefore

. [ VII-15a ] If the random force is stationary in time

[ VII-15b ] and, finally, if is long compared to the random force correlation time

. [ VII-15c ] Therefore Equation [ VII-12 ], the equation of motion for the average number operator, becomes

[ VII-16 ] which is one form of the

famous fluctuation-dissipation theorem.If the reservoir is in thermal equilibrium, we may directly evaluate the noise operator correlation to wit.

[ VII-17 ] Since the operators of the reservoir commute we can write

[ VII-18 ] and Equation [ VII-17 ] becomes

[ VII-19 ] where is of the form of a

in the classical Langevin theory. Two-time correlation functions of this form are characteristic ofdiffusion coefficientMarkoffianrandom processes and represent fluctuations in a reservoir with essentially zero memory. Substitution of this Markoffian correlation function into Equation [ VII-16 ] yields anEinstein relationof the form

. [ VII-20 ]

General Resevoir Problem

As the previous section demonstrates, the effects of random fluctuations in an assemblage of harmonic oscillators may be accounted for in precise detail. The problem can be analyzed completely since the first term on the RHS of Equation [ VII-6 ] -- the-- contains only the system variable. In general, this feedback term will included reservoir variables as well and the simplified model breaks down. The general case may, however, be analyzed if we take the simpler case as a guide. In particular, let us suppose that a set of systems operators are coupled to the reservoir and satisfy a set Langevin equations of motion asdamping term

[ VII-21 ] where is the drift term and is the Markoffian noise operator appropriate to the system operator . From the earlier results, we again assume that the two-time correlation of the random fluctuating force is of the form

[ VII-22 ] Again using the identity

[ VII-23 ] we find

[ VII-24 ] As discussed in connection with Equation [ VII-14 ], the first and second term on the RHS in Equation [ VII-24 ] vanishes

since. Thereforecausecannot precedeeffect

and by Equation [ VII-22 ]

[ VII-25 ]

[ VII-26 ] We may then form the following equation of motion

[ VII-27 ] which leads to a

generalized Einstein relationship

[ VII-28 ] With this equation it is possible to calculate the diffusion coefficients from the drift coefficients, provided one can independently calculate the equation of motion for .

We can also use Equation [ VII-21 ] to calculate the expectation values of two-time correlation functions. In paryicular, if we multiply it by where and average we find

[ VII-29 ] Since in the Markoff approximation, the present system operator cannot know about the future noise source, vanishes and

[ VII-30 ] which is the so called quantum regression theorm which states that the two-time correlation function satisfies the same equation of motion as the single-time .

Semiconductor Langevin EquationsAs discussed earlier, the effective Hamiltonian for optical interactions in a semiconductor may be written as

[ VII-29 ] Thus, the equation of motion motion for the effective dipole operator is given by

[ VII-30 ] where is the damping, , and is the k-dependent noise operator. We transform the various operators to a frame rotating at the electromagnetic field frequency --

viz.

[ VII-31a ]

[ VII-31b ]

[ VII-31c ] and find the transformed equation of motion in the rotating frame --

i.e.the Heisenberg interaction picture -- as

[ VII-32 ] Applying the generalized Einstein relation --

i.e.Equation [ VII-28 ] -- to the operatorwe obtain

[ VII-33a ] Similarly,

[ VII-33b ] where .

The equation of motion for the laser field operator in a frame rotating at the laser frequency may be written as

[ VII-34 ] where is the cold cavity resonant frequency and where we have replace the damping parameter in equation [ VII-10 ] by the experimental parameter . If we consider the Equation [ VII-32 ] in the

quasi-equilbriumlimit, we find

[ VII-35 ] where is the complex Lorentzian demoninator.[7] Thus in the quasi-equilbrium limit, Equation [ VII-34 ] becomes

where

[ VII-36 ]

. [ VII-37 ] Following arguments similar to those used in deriving Equation [ VII-16 ], we easily derive an equation of motion for the laser photon number operator --

viz.

[ VII-38a ] or

[ VII-38b ] where

. [ VII-39a ] and

. [ VII-39b ] Assuming that the cavity damping and carrier scattering reservoirs are uncorrelated, the two-time correlation function may be written as

[ VII-40a ] and if the noise operators for different are uncorrelated

. [ VII-40b ] Finally, using Equations [ VII-19 ] and [ VII-33a ] we obtain

. [ VII-40c ] and, thus, the "" diffusion coefficient for the laser field is

. [ VII-41a ] By a similar argument, the "" diffusion coefficient is

. [ VII-41b ] Using Equation [ VII-41a ], Equation [ VII-38b ] becomes

[ VII-42a ] Similarly

[ VII-42b ] As check on the consistency of these result, we see that the form of Equations [ VII-42a ] and [ VII-42b ] guarantees that the commutator is independent of time --

viz.

[ VII-43 ] to incorporate saturation effects into the noise spectrum we need to write a Langevin equation for the carrier-density operator --

viz.

[ VII-44 ] whereis the pumping due injected carriers, is the radiative recombination (spontaneous emission) term, is the nonradiative recombination term, is the relaxation due carrier-carrier scattering, and is the carrier noise operator. In the quasi-equilibrium approximation --

i.e.from Equation [ VII-35 ] -- we may sum Equation III-44 ] to obtain a Langevin equation for the total carrier-density operator

[ VII-45a] where

[ VII-45b ]

NOISE SPECTRA:To calculate noise spectra we make use of semiclassical transformations --

viz.

[ VII-46a ]

[ VII-46b ]

[ VII-46c ]

[ VII-46d ] where is the electric field per photon. Using the expression

[ VII-47 ] we see that

[ VII-48a ]

[ VII-48a ]

[1] Much of what follows draws heavily on material
in the ** Arizona Books **--

1. M. Sargent III, M. O. Scully and W. E. Lamb, Jr.,Laser Physics, Addison-Wesley (1974)

2. P. Meystre and M. Sargent III,Elements of Quantum Optics, Springer-Verlag (1992)

3. Weng W. Chow, Stephan W. Koch and Murray Sargent III,Semiconductor-Laser Physics, Springer-Verlag (1994)

[2] This model is particular valuable in treating
electromagnetic intermode coupling problems, but it also provides guidence
in treating more general problems which incorporate atomic states.

[3] As is pointed out in Section 14-3 of P. Meystre
and M. Sargent III, *Elements of Quantum Optics*, Springer-Verlag
(1992), the statement that "System operators commute with reservoir operators
at given time..." is a very critical point. From Equation [ II-3 ], we
see that as time evolves, the reservoir operators acquire some of the characteristics
of the system operator. Thus, it is crucial in this development to maintain
a strict and consistent ordering of operators at different times!

[4] The transformation does not alter the system's communtation relationship, so, of necessity, .

[5] In Chapter II, Section 8 of W. Heitler, *The
Quantum Theory of Radiation *(3rd edition)

[6] In the classical Langevin theory of Browian motion, the equation of motion for drift velocity of a particle suspended in a liquid is given by

where
is a damping constant and
is a random, rapidly fluctuation force with zero mean value. In this classical
problem the famous fluctuation-dissipation theorem is expressed as

where i s the so called Fokker-Planck diffusion coefficient.

[7] It may be recalled that in the notes entitled
*Lasers:
Models and Theories* for convenience we indroduced the Lorentzian functions
-- *viz.*

thecomplex Lorentzian denominator

thedimensionless Lorentzian function.

This page was prepared and is maintained by R. Victor Jones,

Last updated May 8, 2000