Добавил:
kiopkiopkiop18@yandex.ru t.me/Prokururor I Вовсе не секретарь, но почту проверяю Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Ординатура / Офтальмология / Английские материалы / Optics Learning by Computing with Examples using MATLAB_Dieter Moller_2007

.pdf
Скачиваний:
0
Добавлен:
28.03.2026
Размер:
11.34 Mб
Скачать

7.4.

BANDWIDTH

289

(w)

 

 

w

 

 

Application 7.6.

1.For the Lorentzian line shape, change the lifetime τ to more realistic values in the range of 103 to 108 and find a corresponding frequency range.

2.Change the lifetime and find the bandwidth. Calculate Q as ω0/ ω and compare with the data from a graph.

7.4.3Natural Emission Line Width, Quantum Mechanical Model

In Section 7.4.2 we saw how the width of the Lorentzian line shape depends on the lifetime τ of the electron in the upper state. The uncertainty principle of quantum mechanics relates the lifetime τ to the width E of the corresponding energy state, in which the electron is placed. The uncertainty principle may be expressed as

Eτ h/2π.

(7.34)

In Figure 7.12 we indicate three possible transitions from the higher to the lower state. One corresponds to the transition from center-to-center of the two bands, and the others to the transition between two different energies of the bands, corresponding to lower and higher frequencies. Quantum mechanics shows that the center-to-center transition has a higher transition probability than the other possible transitions. The shape of the emission line shows these differences. The Einstein coefficient of spontaneous emission A21 describes a transition probability, and is related to the lifetime of the emission process as

A21 1/τ.

(7.35)

7.4.4 Doppler Broadening (Inhomogeneous)

The classical Doppler effect may be observed when a train passes with a blowing horn. The signal seems to have a higher frequency when the train is traveling towards us and a lower frequency when traveling away from us. The atoms in a

290

7. BLACKBODY RADIATION, ATOMIC EMISSION, AND LASERS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FIGURE 7.12 The upper energy level has the width E2, the lower E1. Transitions of higher and lower frequencies are indicated, with respect to the peak frequency of the Lorentzian line shape.

gas also move while emitting light. The frequency shift, related to their particular velocity v, is given as

ν ν0 ν0(v/c).

(7.36)

The distribution of the velocities depends on the temperature and follows Maxwell’s velocity distribution law. One gets

I (ν) Iνo exp(Ek /kT ) Iνo exp(mv2/2kT )

 

Iνo exp{(mc2/2kT )((ν ν0)2/(ν0)2)},

(7.37)

where I0 is a constant, m the mass of the oscillator, c the speed of light, k the Boltzmann constant, and T the absolute temperature. As we did in Section 7.4.2 for the Lorentzian line shape, we define the Doppler line shape as I (ν) I0gd(ν) and get

gd(ν) [(

 

 

 

π(ln 2)(2/π ν)] exp{−(ν ν0)2(ln 2)/( ν/2)2}.

(7.38)

The line shape of Eq. (7.38) has a Gaussian profile which is different from the Lorentzian, schematically shown in Figure 7.13a.

The halfwidth at half height is calculated to be

ν 2ν0 [(2kT /mc2)(ln 2)].

(7.39)

The Gaussian line shape may be looked at as the envelope of a large number of Lorentzian line shapes. Each oscillator moves with a different velocity in a different direction and the peak frequencies of the emitted light have a Gauss distribution, (shown in schematically Figure 7.13b). When discussing an exam-

7.5. LASERS

291

∆ν

FIGURE 7.13 (a) Comparison of Lorentzian and Gaussian line shapes of approximately same line width at half height; (b) Gaussian line shape is the envelope of Lorentzian line shapes emitted statistically at different velocities in different directions.

ple of the Ne–He laser below, we give numerical values for the halfwidth of Eq. (7.39).

7.5LASERS

7.5.1Introduction

We discussed blackbody radiation in Section 7.2, atomic emission in Section 7.3, the Fabry–Perot in Chapter 2, and modes in Chapter 6. We now discuss the following simplified model for a two-level laser. We consider a Fabry–Perot filled with an active medium of atomic oscillators. First we need population inversion. The oscillators will be excited into the upper state 2 by using the light of frequency νP from the outside. We need the excited states to remain excited for sometime, for example, 103sec. Then we need stimulated emission to transfer the energy of the excited oscillators to the modes of the Fabry–Perot. To have laser light emitted from the Fabry–Perot, one must couple out part of the light of the modes. This is done by making one of the mirrors of the Fabry–Perot not totally reflective, such as when there is a small hole in the mirror.

292

7. BLACKBODY RADIATION, ATOMIC EMISSION, AND LASERS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FIGURE 7.14 Schematic of a three-level laser. S23 is the probability of a radiationless transition, Wik and Wki are probabilities of induced absorption and emission, and Aik is the probability of spontaneous emission.

7.5.2 Population Inversion

7.5.2.1 Two-Level System with Stimulated and Spontaneous Transitions

When discussing the blackbody radiation law, we used for the occupation of the different energy levels in thermal equilibrium the Boltzmann distribution NeE/kT , where k is Boltzmann’s constant. This distribution law tells us that, at temperature T , the lower states are more occupied than the higher states. For laser action we need just the reverse, which is more electrons in a higher energy state. This is called population inversion and is in contrast to the population of the electrons in thermal equilibrium.

We assume that we have to deal with oscillators having just two levels, as shown in Figure 7.14. We call the upper level E2 and the lower level E1. The change in time of the number of photons N1 of E1 and N2 of E2, the rate equations, is

dN1

/dt N1B12u(ν) + N2B21u(ν) + N2A21

(7.40)

dN2

/dt +N1B12u(ν) N2B21u(ν) N2A21,

(7.41)

where u(ν) is the radiation density. The coefficient B12 is the Einstein coefficient of stimulated absorption. The coefficient B21 is the Einstein coefficient of stimulated emission and A21 is the coefficient of spontaneous emission. We have used these coefficients when deriving Planck’s radiation law, (see Section 7.2), and used the relation

A21/B12 8πh(ν/c)3

(7.42)

as well as B12 B21. This relation has to be slightly changed for atomic oscillators. We have to take into account that the atomic energy levels may be degenerate, which means there are several energy levels having the same energy. For example, if there is a threefold degeneracy, we have to use the weight g 3

7.5. LASERS

293

for the transition. Therefore we have to use

 

g1B12 g2B21.

(7.43)

7.5.2.2 Changes in the Upper Level Considering Stimulated Transitions

A necessary condition for laser action is population inversion which as stated before, means there must be more photons in the upper state than in the lower state. Using only the stimulated emission and absorption processes, we have for the number of photons in the upperstate.

N1B12u(ν),

(7.44)

and the number in the lowerstate

 

N2B21u(ν).

(7.45)

The change in the number of photons in time is then

 

dnp/dt N2B21u(ν) N1B12u(ν).

(7.46)

With Eq. (7.43) we write

 

dnp/dt [N2 N1(g2/g1)]B21u(ν).

(7.47)

To have laser action one needs

 

 

g2

 

N2 N1

 

> 0

(7.48)

g1

and one has the condition for laser action, N2/g2 > N1/g1.

7.5.3Stimulated Emission, Spontaneous Emission, and the Amplification Factor

To have the condition N2/g2 > N1/g1 fulfilled, one has to make the stimulated emission larger than the spontaneous emission. The stimulated emission transfers the photons from the upper energy state E2 into the modes of the Fabry–Perot. The bandwidth for the emission process was discussed above and we assume that we have a Lorentzian line shape. The modes of the Fabry–Perot are much narrower and close to delta functions (see Figure 7.15). The probability of stimulated emissions per unit time is calculated by considering Eq. (7.46), written for a number of n photons as

dn/dt [N2 (g2/g1)N1]B21u(ν).

(7.49)

The radiation density is u(ν) hνn. Taking into account the line shape of the emission process, (see Eq. (7.30)), we have to multiply the coefficient of stimulated emission B21 by u(ν)gl(ν), that is, B21u(ν)gl(ν), and get

du(ν)/dt [N2 (g2/g1)N1]hνgl(ν)B21u(ν).

(7.50)

294

7. BLACKBODY RADIATION, ATOMIC EMISSION, AND LASERS

ν−1 ν0 ν1

FIGURE 7.15 The width of the modes of the cavity at frequencies ν1, ν0, and ν1 are almost delta functions compared to the bandwidth of the Lorentzian-shaped emission line.

We can express the change of u(ν) with respect to the time dt by considering the length dz in which the light travels in the time dt This gives us dt dz/c . Here we use the speed of light c , modified in the medium of the laser cavity. The change of the radiation density over the interval dz is then

du(ν)/dz [N2 (g2/g1)N1](hν/c )gl(ν)B21u(ν, z).

(7.51)

Using A21 1, and A21/B21 8πh(ν/c )3, we have

 

du(ν)/dz {(c 2/8πν2τ )[N2 (g2/g1)N1]gl(ν)}u(ν, z).

(7.52)

The expression in the curly braces is called the amplification factor ε(ν):

ε(ν) {(c 2/8πν2τ )[N2 (g2/g1)N1]gl(ν)}.

(7.53)

The gain of the beam ε(ν) is “per length” and an example for the Ruby laser is calculated in FileFig 7.7.

FileFig 7.7

(L7RUBYS)

 

 

 

´

Gain calculation for the example of the Ruby laser. (See also Kneub’uhl and Sigmst, 1988.)

L7RUBYS is only on the CD.

7.5.4The Fabry–Perot Cavity, Losses, and Threshold Condition

The Fabry–Perot cavity has two plane parallel mirrors with reflectivities R1 and R2, and the modes are standing waves. These standing waves may be considered as two traveling waves, moving in opposite directions. If R1 and R2 have high reflectivity, the traveling waves will pass forward and backward through the volume filled with the oscillators in the excited state E2. By stimulated emission,

7.5. LASERS

295

energy will be picked up and the intensity of each mode will increase. This process will come to an end when the increase in energy is set off by the losses. Large distances between the two mirrors of the Fabry–Perot correspond to a high order of longitudinal modes and these high-order modes correspond to a very narrow bandwidth, schematically shown in Figure 7.15. The intensity depending on the length z of the light traveling in the active medium of the Fabry–Perot is related to the radiation density as I (z) c u(ν, z). From Eq. (7.50) and (7.51) we have

dI (z)/dz ε(ν)I (ν, z).

(7.54)

Integration gives us

 

I (z) I (0)eε(ν)z.

(7.55)

We see the exponential increase of the light traversing the active medium; energy is transferred to the modes of the Fabry–Perot.

In order to use some of the energy one has to couple it out of the cavity. This may be achieved by using a mirror with a small hole for one of the two mirrors of the Fabry–Perot. The length of one round trip is z 2L and the losses are taken into account by the factor (α(ν)) in the exponent of Eq. (7.55). We multiply Eq. (7.55) with the reflectivities R1R2 to account for the reflection losses of one round trip and get

I (2z) I (0)e[ε(ν)α(ν)]2L(R1R2).

(7.56)

The thereshold condition is then obtained from Eq. (7.56) when

 

e[ε(ν)α(ν)]2LR1R2 1

(7.57)

or [ε(ν) α(ν)] (1/2L)ln(1/R1R2). The threshold gain is

 

εT (ν) α(ν) + (1/2L)ln(1/R1R2)

(7.58)

We rewrite Eq. (7.58) by calling σ [N2 N1(g2/g1)] and insert εT (ν) in Eq. (7.53) and get

σT [α(ν) + (1/2L)ln(1/R1R2)](8πν2τ/c 2)(1/gl(ν)).

(7.59)

´

This is the famous threshold condition by Schawlow and Townes (Kneub’uhl and Sigrist, 1988, p.38). In FileFig 7.8 we present a calculation with numerical values for the He–Ne laser.

FileFig 7.8 (L8HENES)

Gain calculations for the example of the Ne–He laser. We have used for m in mc2 the mass of the proton times 20 for Ne. The Lorentzian line width has been

296 7. BLACKBODY RADIATION, ATOMIC EMISSION, AND LASERS

replaced by the Doppler line width. The line shape gd(ν) is approximately 1 divided by (2ν0 (2kT /mc2) ln 2) (see Eq. (7.39).)

L8HENES is only on the CD.

7.5.5 Simplified Example of a Three-Level Laser

We consider an energy schematic of the oscillators shown in Figure 7.14 with an upper state 3, a lower state 2, and the ground state 1. We have for transitions between 3 and 1 induced absorption, induced emission, and spontaneous emission, as discussed in Section 7.5.2. The transition probabilities are called W31, W13, and A31, respectively. A similar description holds for the transitions from 2 to 1. The transitions from 3 to 2 are now special. They are radiationless transitions and their probability is called S32. The occupation number of the oscillators in the states 1, 2, and 3 are called N1, N2, N3 and the total number N0 is assumed to be a constant.

N0 N1 + N2 + N3.

 

(7.60)

The change in time of the number of oscillators in state 3 is

 

dN3/dt W13N1 (W31

+ A31 + S32)N3

(7.61)

and in state 2

 

 

dN2/dt W12N1 (W21

+ A21)N2 + S32N3.

(7.62)

For the steady state, the time derivatives are zero. We assume that A31 W31, which gives us for the ratio N2/N1 of the oscillators

N2/N1 {(W13S32)/(W13 + S32) + W12}/(A21 + W21).

(7.63)

From the discussion of blackbody radiation, W13 W31 and W12 W21, and we assume that the radiationless transition probability S32 is much larger than the transition probability W13 (the probability to empty, state 3). We then have from Eq. (7.63)

N2/N1 {(W13 + W12)/(A21 + W12)}.

(7.64)

Using Eqs. (7.60) as N0/N1 and (7.61) to (7.63) we may get after some calculations

(N2 N1)/N0 {(W13 A21)/(W13 + A21 + 2W12)}.

(7.65)

The threshold condition for operation is obtained when, in the steady state, there are as many oscillators in state 1 as in state 2, N1 N2. From Eq. (7.65) one has W13 A21. The minimum power for operation is calculated from the condition that the power corresponding to induced absorption of level 3 is equal to the power corresponding to the spontaneous emission of state 2. In FileFig 7.9 we calculate P NA21for a process involving a metastable state and spontaneous emission.

7.6. CONFOCAL CAVITY, GAUSSIAN BEAM, AND MODES

297

FileFig 7.9 (L9MINPOWS)

Model calculation of minimum power for operation of a laser, P NA21:

1.for a metastable state, and

2.for spontaneous emission.

L9MINPOWS is only on the CD.

Application 7.9. How would the numbers change if S32 must be taken into account?

Take S32 10W13 or S32 1000W13.

7.6CONFOCAL CAVITY, GAUSSIAN BEAM, AND MODES

7.6.1 Paraxial Wave Equation and Beam Parameters

Laser light travels forwards and backwards between the two mirrors of a resonance cavity. The first Ne–He laser used a Fabry–Perot cavity of a length of one meter between the two plane parallel mirrors and was extremely difficult to align. Later it was found that the confocal cavity was much easier to align, and was also very efficient from the point of view of diffraction losses. In Chapter 1, we discussed several types of resonators from the perspective of geometrical optics, and in Chapter 2 the Fabry–Perot as a resonance interferometer.

In this section we discuss the confocal cavity with radii of curvature Rm d, where d is the distance between the two mirrors placed at z ±d/2. Here the origin of the coordinate system is at z 0 in the middle between the mirrors. The laser beam in the cavity is approximated by a wave traveling in the z direction and having a bell-shaped profile in the transversal (x, y) direction (Figure 7.16). For the mathematical presentation of these propagating modes, we use Cartesian coordinates and rectangular mirrors. We follow H. Kogelnick and

FIGURE 7.16 Waist of a Gaussian beam in a confocal cavity depending on the z coordinate.

298

7. BLACKBODY RADIATION, ATOMIC EMISSION, AND LASERS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FIGURE 7.17 (a) Coordinates for discussion of the radius of curvature of the wavefront of a mode in a confocal cavity; (b) radius of curvature of the wavefront is indicated for different values of the z coordinate.

T. Li and start with the scalar wave equation

 

δ2u/δx2 + δ2u/δy2 + δ2u/δz2 + k2u 0,

(7.66)

where k 2π/λ and try to find a solution of the wave equation of the form

u ψ(x, y, z) exp(ikz).

(7.67)

Inserting Eq. (7.67) into Eq. (7.66) and assuming for the calculation that the second derivative δ2ψ/δz2 may be neglected, we further consider the paraxial wave equation

δ2ψ/δx2 + δ2ψ/δy2 2ikδψ/δz 0.

(7.68)

Solutions of Eq. (7.68) describe a Gaussian beam profile in the transverse direction, depending on the distance r from the axis, where r2 x2 + y2 (Figure 7.17). A solution of Eq. (7.68) is

ψ exp{−i[P (z) + k(x2 + y2)/2q(z)]}

(7.69)

with

δq(z)/δz 1

(7.70)

and

δP (z)/δz i/q(z),

(7.71)

where the two functions q(z) and P (z) are not independent of each other. We are mainly interested in the solution for q(z) and can solve for P (z) if we need it, which is a phase factor.

We write the solution of Eq. (7.70) as

q(z) izR + z,

(7.72)