Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Страховский Г.М. Основы квантовой электроники учеб. пособие

.pdf
Скачиваний:
60
Добавлен:
25.10.2023
Размер:
18.19 Mб
Скачать

Эта функция удовлетворяет уравнению Шредингера с гамильто­

нианом (1.35):

 

 

 

 

 

 

 

 

 

mJ± = mp.

 

(1.39)

 

 

 

 

dt

 

 

Теперь сформулируем задачу. Необходимо вычислить

коэффициен­

ты Эйнштейна Атп

и Втп.

Для этого можно воспользоваться форму­

лой (1.3а),

связывающей

три величины: вероятность поглощения в

единицу времени

Wnm,

спектральную плотность

энергии

электромаг­

нитного поля р ѵ

и коэффициент Эйнштейна Впт.

Если вычислить ве­

личину Wnm

для известного

значения р ѵ , то их отношение дает коэффи­

циент Эйнштейна Впт.

Зная же коэффициент Впт

~ Втп,

из формулы

(1.13) можно определить коэффициент Эйнштейна

Атп.

 

Величина Wnm выражается через коэффициенты an(t) волновой функции (1.38) возмущенного состояния. Вычислим эти коэффициенты. Для этого подставим функцию вида (1.38) в уравнение (1.39):

т ( S - ^Г Ч>по + 2 а п ~ Г )=Н*Ъап

У по + Н'2ап

$ п 0 .

Последний член левой части и первый член правой части этого ра­ венства взаимно приводятся вследствие соотношения (1.36). В резуль­ тате получаем:

і й 2 - ^ „ о = Я в 2 а , Ж о -

(1-40)

Умножая правую и левую части этого равенства на \j)|0 " инте­ грируя по объему V с учетом условия ортогональности волновых функ­ ций, имеем:

m^JL ^^aJ^lofr^dV.

(1.41)

dt п V

В интеграл в правой части уравнения подставляем явный вид вол­ новых функций а|з|0 и г|5п 0 . Тогда

S ^ * o Ä > n 0 d V

= exp ^^~^t^u%HsundV.

 

(1.42)

Если ввести обозначения

 

 

 

^ ^ г ^ .

 

(1-43)

 

h

 

 

 

HL=~=\t*kHaundV,

 

(1.43а)

 

V

 

 

то уравнение (1.41) примет окончательный вид

 

 

ІП

dt = 2аn пHinexp(icohn

t).

(1.44)

20

Найдем его решение, считая гамильтониан взаимодействия малым возмущением.

Предположим, что при отсутствии возмущения система находилась только в одном из стационарных состояний, например в состоянии т. Это значит, что в момент, когда было включено возмущение (t = 0), коэффициенты a{k0) были равны:

4 ° Ч °

k ¥ =

m '

(1-45)

(I

k =

m.

 

Это нулевое приближение; индекс (0) при коэффициенте ah

и обо­

значает нулевое приближение. Для определения первого приближе­

ния

найдем ah

в виде ak = û40 ) -\- a(k)

+

т. е. для состояния

k =

m

a m

= 1 + йт)

+

а для остальных

состояний ak •— 0 ^ ' +

... .

Пра­

вая часть уравнения (1.44) содержит малые члены Hin (малое возму­

щение), поэтому в нее следует подставить выражение

ап =

1 (нулевой

порядок) для п = т. Все остальные коэффициенты

а(п0)

при пфт

равны нулю, поэтому ьместо суммы остается только один член. В пер­ вом приближении вместо уравнения (1.44) будем иметь

 

dt

іѢ

# L e x p (i(ùkmt)

(1.46)

[индекс (1) при коэффициенте ak

означает, что речь идет о вычислениях

в первом приближении].

 

 

 

Уравнение (1.46) легко проинтегрировать. В результате получим:

ар

=-L\Hlnexv(i<ùhmt')dt'.

(1.47)

 

 

о

 

 

Далее необходимо вычислить

матричные элементы

гамильтониана

взаимодействия Н\т.

Для этого

нужно знать вид гамильтониана вза­

имодействия.

 

 

 

 

Рассмотрим наиболее простую модель частицы (атома). Пусть атом состоит из бесконечно тяжелого ядра с эффективным потенциалом Ѵа и электрона с массой m и зарядом е. Тогда для электрона, движущего­ ся в электромагнитном поле с векторным потенциалом Л, скалярным потенциалом <рц и в стационарном поле с потенциалом Ѵа, полный га­ мильтониан имеет вид

где р — оператор импульса

21

Раскроем круглую скобку в выражении для гамильтониана и учтем, что коммутатор операторов р и А имеет вид рА— Ар = — ihàivA. Тогда гамильтониан H приобретает вид

H = ~ J

r

Уп - — Ah-

!г- <ЗіѵЙ +

2тс2

Ä2

+ е>Фп. (1.48а)

 

тс

2тс

 

 

С помощью градиентного преобразования можно всегда выбрать потенциал электромагнитного поля в отсутствие зарядов таким обра­ зом, что ф п = 0 и div А = 0. Кроме того, обычно член, пропорцио­ нальный А2, значительно меньше члена, содержащего векторный по­ тенциал в первой степени, поэтому им можно пренебречь. Учитывая это, получаем выражение для полного гамильтониана системы части­ ца плюс поле:

H = U

L + VB)

e~Äp.

(1.49)

\

I

тс

 

Члены, заключенные в круглую скобку, представляют собой не­ возмущенный гамильтониан Н°, а член --Ар — часть полного

гамильтониана, определяющую взаимодействие. Таким образом, га­ мильтониан взаимодействия имеет вид

Нв =

—Ар.

(1.50)

 

тс

 

Пусть электромагнитное поле представляет собой монохроматиче­

скую волну с волновым вектором k, т. е.

 

 

 

A (rt) = Л exp [i (kr(i>t)\ -f- Ä\ exp [—i (k~r — cat)].

(1.51)

Тогда гамильтониан взаимодействия Нв

имеет вид

 

Нв =

—{Ai [рехр (/ kr) ] exp ( — Ш) + Ä* [р ехр (— i k г) ] exp

.

 

тс

 

(1.52)

 

 

 

Вычисление матричных элементов Н\т

сводится к вычислению мат­

ричных элементов выражений, стоящих в квадратных скобках этого выражения, а именно:

Н\т —

{Aï exp (ik~r)]km

ехр

- f

 

тс

 

 

+

А\ [рехр (— ikr)]hmехр

(Ш)}.

( 1.53)

22

 

 

 

Подставляя матричные элементы (1.53) в уравнение (1.47) и инте­ грируя его, получаем:

ai1» (t) = ~

At

[рехр( ^r)lf t ,n^lS^=^Ih±

+

/ясл

 

й ) д т — (о

 

+ - V Л?[ p e x p ( - ï 7 ë r ) J A m ^ ^ - ^ i b i .

(1.54)

Из этого выражения видно, что первый член в нем велик, если ча­ стота электромагнитной волны со близка к частоте перехода coA m , т. е.

если выполняется

условие « ж w h m

--= ———™, или

W f t ^ W m + f i « .

Таким

образом, первый член

связан

с переходом в состояние k, кото­

рое выше

состояния m на

величину энергии На. При таком

пе­

реходе

внутренняя

энергия

частицы

увеличивается

за счет энергии

электромагнитного

поля, т. е. первый член описывает процесс резо­

нансного

поглощения

(или

просто

поглощения). Второй же член

выражения

(1.54)

велик, если выполняется условие

со = — cof t m

=

Wm~~K W~'

т '

е ' ^ h ~

^ m

Следовательно, в этом случае состо­

яние k ниже состояния

m на

величину энергии Йю. Тогда при пе­

реходе в состояние k внутренняя энергия частицы уменьшается на ве­

личину Гш, энергия отдается электромагнитному полю. Таким образом,

второй член в выражении (1.54) описывает переход за счет индуциро­

ванного излучения.

 

 

 

 

 

 

Поскольку нашей задачей является вычисление вероятности по­

глощения,

входящей

в

формулу (1.3а), то в коэффициенте ai1 '(О

м ы

оставим лишь первый член, описывающий процесс поглощения.

Известно,

что вероятность обнаружить систему в энергетическом

состоянии k в момент времени t равна | ah(t)

| 2 , т. е. | a{kl)(l)

| 2 в рамках

данной задачи

 

 

 

I о I 2 - -

i V

I ^ I2 1 ^ е х р

]*» I 2 L [ e x p ' / m

f e m ~ t ~ 1 1 1 2 - ( L 5 5 )

 

т2сЧ2

 

(Щгп — <ù)2

Здесь вместо скалярного произведения Аф введена проекция вектора р на направление вектора Ах, а именно рА. Тогда Агр = АгрА- Приведем формулу для квадрата модуля матричного элемента, входящего в вы­ ражение (1.55), в дипольном приближении без вычислений:

\ І Р а

exp(ikr)]km\'=

/ Г Г Ю km I - [2

(1.56)

— ^ - \ r k m \ \

где r h m — матричный

элемент оператора смещения.

 

Нетрудно видеть, что квадрат модуля в выражении (1.55) можно

преобразовать к виду

 

 

Uexpi'Km — (ù)t 1]|2

=

 

= 2 [ 1 - c o s K m - и) t] = 4 sin2

[ ^ыь-ѵѴ j .

(1.57)

23

Подставляя выражения (1.56) и (1.57) в равенство (1.55), получаем, что вероятность нахождения системы в состоянии к в момент времени / равна:

 

 

 

 

2 Г Л 2

 

4 sin2

(

 

 

 

J t

 

 

 

 

 

 

 

W | F f

c J '

t

K

 

\ t

'

(1-58)

 

 

 

 

ЗА2

с 2

 

 

 

 

 

( м

Д т с о ) 2

 

Рассмотрим это выражение, т. е. изучим поведение множителя

 

 

 

 

 

 

4 sin2

 

 

/

 

 

 

 

 

 

при достаточно больших значениях

/.

v>hm

 

 

 

 

 

 

 

Во-первых, если выбрать величину

 

— и

столь

малой,

что

(ОЪт — СО

1

 

СОьт

— СО , . .

 

 

 

 

 

 

 

 

Л

2

<^ - ,

или

2

г <

1,

т. е.

рассмотреть

множитель

Qx

вблизи точки

со =

km,

то sin ( — ^ ^ 2 — м о ж н о

разложить в ряд,

ограничившись первым членом разложения. Тогда

 

 

 

 

 

 

 

 

 

Qi =

=

t.

 

 

 

 

 

 

 

 

 

 

((ühm

— (ü)2t

 

 

 

со = ыІіт

 

 

Таким образом, множитель Qt

вблизи

точки

возрастает

с ростом t и при / - >

оо становится бесконечным.

 

 

 

 

 

Во-вторых,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

Q,ci(Aco') = j ^ - ( 4 ^ - d ( Ä c o ' )

=

K,

 

(1.59)

где

введено

обозначение

Асо'

= [ ù

f e m ~ g >

. Но

оба

эти

условия озна­

чают, что множитель Q1 удовлетворяет всем требованиям, предъявля­ емым к ô-функции. Поэтому при больших значениях t его можно аппрок­

симировать

б-функцией,

а именно

записать Q1 =

по ^ Ю й т 2 ~ т

^ . Мно-

4 5 І П 2

^ m f t m - t O j ,

 

 

 

 

 

 

житель 7

гт;

, очевидно,

можно

тогда

представить

в

виде

. . о і

" f e m CO , ,

tQ1 = ntb

(

) = 2 ^ Ô K m - t o ) ,

(1.60)

4 sin2

t

(Cùftm — C O ) 2

 

V

2

 

 

 

где мы воспользовались следующим свойством ô-функции: ô(ax) =

— Êi^L (х — аргумент, а — некоторый множитель).

24

Подставляя формулу (1.60) в выражение (1.58), получаем:

2™ 2

 

Hm I а ^ Ч О I я = ~ ^ г | Л Г Ы 2 2 я / О К ; п - а > ) .

(1.61)

Разделив это выражение на t, получим вероятность перехода-в е д и ­ ницу времени*:

'

t

= - ^ І Л і Н ^ т Г 2 я о K m - ö , ) .

(1.62)

Наличие ô-функции в формулах (1.61) и (1.62) выражает закон со­ хранения энергии. Именно, вероятность процесса поглощения стано­ вится бесконечной при резонансе (со = cofem ) и обращается в нуль,

если со ф(дкт.

Спектральная линия перехода в рассмотренном приближении для системы с дискретными уровнями энергии бесконечно узка. Следует отметить, что при учете конечной ширины энергетических уровней, или учете ширины линии, связанной с индуцированным излучением, вероятность поглощения, хотя и максимальна при выполнении резо­ нансного условия (со = o)7 i m ), но конечна, а спектральная линия пе­ рехода имеет конечную ширину. В конце этой главы дается строгое решение задачи о двухуровневой системе в электромагнитном поле. Приведенное там решение может служить иллюстрацией этого поло­ жения.

Теперь определим входящую в формулы величину | ^ x | 2 через аб­ солютную величину вектора потока энергии S:

 

 

 

 

 

 

 

дА 12

 

 

 

 

 

 

 

 

 

An

4лс

dt

 

 

 

 

 

Подставляя

сюда выражение (1.51) и считая, что комплексный век­

тор Аг может

быть представлен в виде

У4х

-= е |

| ехр

(гЧр), где ф —

некоторая

фаза, а е — единичный

вектор в направлении вектора

Аи

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\s\

= — sm2(îir~(ùt~(f)\A1\2.

 

 

 

 

(1.63)

 

 

 

 

ПС

 

 

 

 

 

 

 

Усредним величину потока энергии

по периоду

колебаний

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

|S |

= — IAA2.

Отсюда

 

 

 

 

 

 

 

 

 

 

 

2пс

 

 

 

 

 

 

 

 

 

 

 

 

 

Л Г =

Мер -

 

 

 

(1.64)

*

При

поставленных

в задаче начальных

условиях | a

k (t) f

есть не только

вероятность

нахождения

частицы на уровне

k, но и вероятность перехода

на

этот уровень. Знак l i m в

(1.61) и дальше означает «при

очень

больших tu

 

25

Подставляя выражение (1.64) в формулу (1.62), получаем:

hm

I а * ] 4 0 I я е * < п Г

,» 4п»с

(1.65)

__!_ ^ _____ I r f t m I . IS | с р 6 ( и к т — с о ) .

При выводе этой формулы мы считали электромагнитное поле стро­ го монохроматическим. Между тем оно обычно обладает некоторой ко­ нечной шириной спектра. Поэтому в общем случае вместо средней по периоду колебаний величины потока энергии следует ввести спектраль­ ную плотность интенсивности излучения S (со) и для спектрального интервала ширины do> заменить | 5 | с р на произведение S((a)d(ù. Фор­ мула (1.65) в этом случае будет определять вероятность перехода в еди­ ницу времени под действием излучения в полосе частот doi и примет вид

 

 

-

I

 

 

s м

«

к , - » ) . - ( 1 .ее)

 

 

I

ЗТі2с^

 

со2

 

 

 

 

Проинтегрируем это выражение по узкому спектру частот в

области

со «

coftTO. Результат дает

полную

вероятность

перехода

с поглощением кванта,

т. е.

Wkm:

 

 

 

 

 

 

 

^ m = ^ S K J k f t m | 2 -

 

 

(1-67)

Наконец,

выражая

спектральную

плотность

интенсивности излу­

чения

через

спектральную

плотность

энергии

по формуле

S (со) =

SM

с

 

 

 

 

 

 

 

 

ѵ / •р , получаем

 

 

 

 

 

 

 

 

 

^ m

= | f Ы

Ч -

 

 

(1.68)

Для получения коэффициента Эйнштейна Bhm в дипольном при­ ближении сравним формулу (1.68) с формулой (1.3а). Имеем

Bnm-=^fY'\m\2-Bmh-

 

 

 

 

(1-69)

Для того чтобы получить в явном виде второй

 

коэффициент

Эйнштейна A m h , воспользуемся установленным

в § 1.1

[см. форму­

лу (1.13)] соотношением между коэффициентами Вктп

и

Akm:

8nhvlm

б 4 л 4 е 2

і - is

4co|m е 2

 

| 2

Hhm — "hm

Т 7 . Т

I rhm

I —

$сз

\rhm\ >

 

3hXk

 

 

 

 

 

где были приведены два выражения для коэффициента Эйнштейна

Ahm' один — через частоту перехода a k m

(как и коэффициент

Bkm),

другой — через длину

волны излучения,

соответствующую переходу

k <-» m.

 

 

 

Полезно провести

некоторые оценки. Пусть величина е \ r h m

| (ди-

польный момент перехода) равна 5 • Ю - 1 8

СГС (5 дебаев). Это соответ-

26

ствует

разрешенному

электрическому дипольному переходу.

Пусть

Ккт --

1 см ( Ы О 4 мкм) (радиодиапазон).

Тогда

из формулы

(1.70)

следует, что Ahm=0,8

Ю - 5 сек _ 1 . Если же величина e | r h m | порядка

JQ-2Ö

с г с

(магнитный

дипольный

переход), то

A h m Ю - 1 1

сек.-1.

Величина

т с п — -^—

определяет

время

жизни

частиц на уровне k

 

 

Акт

 

 

 

 

 

за счет спонтанных переходов на уровень т. Проведенные оценки по­ казывают, что в радиодиапазоне это время очень велико: ж 105 сек (для электрических дипольных переходов) и » 10й сек (для магнитных дипольиых переходов), и, следовательно, спонтанное излучение в ра­ диодиапазоне не играет обычно большой роли. В оптическом диапазоне роль спонтанного излучения существенно возрастает. Действительно, если Xkm = 0,4 мкм = 4 - Ю - 5 см (фиолетовая граница видимой об­ ласти спектра), то, как следует из формулы (1.70), коэффициент Эйн­ штейна Akm возрастает примерно в 1015 раз (остальные величины в формуле берем такими же, как в проведенных выше оценках для радио­ диапазона).

§ 1.5. Форма и ширина спектральной линии

До сих пор мы говорили об энергетических уровнях, предполагая, что они бесконечно узки. Между тем энергетические уровни имеют некоторую ширину, т. е. энергия каждого атомного состояния не яв­ ляется строго фиксированной, а несколько размыта. Следовательно, существует некоторое распределение интенсивности поглощения (из­

лучения)

по

частоте,

или,

говорят, линия поглощения (излучения)

имеет некоторую форму.

Количе­

ственно

эта "форма

характеризует­

ся функцией

g((ù),

которая

носит

название

 

ф о р м - ф а к т о р а

спектральной

линии,

или

 

просто

формы линии. Функция g(a)

нор­

мируется

таким образом, что

 

 

ou

 

1.

 

 

(1.71)

 

Jg(<ö)rf«Ö:

 

 

Парис . 1.3 показана типичная форма линии излучения.

Важнейшей количественной ха­ рактеристикой линии является ее

Рис. 1.3. Типичная форма линии из­ лучения

ширина.

Шириной линии называется

интервал

частот А со

между

точками,

для которых

интенсивность

излучения

(поглощения)

падает

в два раза

по сравнению

с максимальной.

Иногда оперируют понятием не

ширины, а полуширины спектральной линии.

Некоторая размытость энергетических уровней и связанное с этим появление у спектральной линии конечной ширины может быть обус-

27

ловлено различными причинами. В следующей главе перечисляются причины уширения линий в газовой радиоспектроскопии. Сейчас же мы рассмотрим лишь одну из этих причин—конечную ширину энерге­ тических уровней, связанную с конечностью времени жизни частицы на энергетическом уровне (естественная ширина линии).

Дело в том, что если время жизни частицы на энергетическом уров­

не равно At, то согласно соотношению

неопределенностей энергия

h

 

уровня имеет неопределенность A F Ä ; ^ ,

Т. е. уровень «размазан».

Отметим, что естественная ширина линии определяет тот предел, ниже которого ширина спектральной линии быть не может.

Для вычисления формы спектральной линии, связанной с конеч­ ностью времени жизни частицы на энергетическом уровне, можно ис­ пользовать результаты, полученные в предыдущем параграфе. Конечность времени жизни частицы на уровне проще всего учесть,

введя в волновые

функции невозмущенного состояния

затухание.

В этом случае

 

 

yno

= unexp(-jrWnt-^y

(1.72)

т. е. в отличие от формулы (1.37) в показателе экспоненты

содержится

еще дополнительное

слагаемое — "Ц^-, характеризующее

конечность

времени жизни частицы на энергетическом уровне. Действительно, вероятность нахождения частицы в состоянии с энергией Wn пропор­ циональна величине 'фпо'Фпо. т - е - пропорциональна ехр (—ynf), дру­ гими словами, как функция времени затухает с константой затухания у п . Время жизни частицы At на уровне, таким образом, равно At =

~~Уп'

Повторим все выкладки предыдущего параграфа, используя в ка­ честве волновых функций невозмущенного состояния функции (1.72). Проще всего это сделать следующим образом. Подставим волновые функции (1.72) в уравнение (1.41). Тогда интеграл (1.42) примет вид

J %о Hipпо dV = ехр

Уи + Уп , • « V

•Wn

 

j UkH°tihdV.

V

 

 

Введя обозначения

 

 

Wk — Wn

Vk + Уп

У h - f r V n

ш * » = — *

5 г - = ш » » - - ^ г " '

mn

= $ti%HBundV,

(1.73)

 

V

 

приходим сразу к уравнению (1.44), в котором вместо (ofe n стоит (ù'kn. Поэтому все последующие выкладки остаются теми же, но надо всюду поменять частоты hn [см. формулу (1.43)] на комплексные частоты

28

(i>'kn [см. формулу (1.73)1. В результате вместо выражения (1.54) полу­ чаем

" Y f t,

ехр

Ay [pexp{ikr)]hn

 

l (Cùftm— Cû)-

Y m +

Yft

 

 

 

 

 

ехр

 

- Yft.

 

« (О>Й/Л + Ю ) < -

 

 

(1.74)

me, Л*[рехр( — / £ r ) ] f t m

 

Y m +

Yft

Действительно, обратимся, например, к показателю экспоненты первого члена уравнения (1.54), т. е. к i(ahm — со)/. В соответствии с указанным выше его следует заменить на

І (©*„, — ©) t = i[ кт Щ^ — <Й И = j K , „ -

Рассмотрим опять, как и в предыдущем параграфе, процесс погло­ щения, и поэтому оставим в решении только первый член. Кроме того,

пусть t ^> — \ — . Это условие означает, что нас интересует решение

Уh + Уm

в достаточно далекие моменты времени (установившийся процесс). Поскольку t > , т. е. ^ + ^"1 / ^> і( члены, определяю­

щие процесс установления, затухают по экспоненте и ими можно пре­ небречь. Тогда

<A»(t) =

- т г Л [pexp(ikr) Iftm

1

 

/

 

Y m +

 

 

 

 

 

 

Yft

 

 

 

 

 

« (»fem С О )

 

 

а вероятность

перехода

 

 

 

 

 

е

1 А |2

I [рлехр {ik r)]h

12

 

 

 

(1.75)

i 2 c 2 Â 2

 

Jftm I

 

( ( 0 A m —Cù)2 +j

Y m +

Yft

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

где как и в формуле (1.55), введена

величина Р А —проекция век­

тора р на направление вектора А1. Наконец, подставляя в (1.75) выра­ жение (1.56) для квадрата модуля матричного элемента | [exp(ikr)pA]hm |2

и выражение (1.64) для квадрата модуля вектора Аи

получаем:

^ 1 ) | 2 = _2яе2_ | 5

(1.76)

ср Irkra \

кт ш)2 + Y m +

Yft

29

Соседние файлы в папке книги из ГПНТБ