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

книги из ГПНТБ / Волновые и флуктуационные процессы в лазерах

..pdf
Скачиваний:
25
Добавлен:
22.10.2023
Размер:
14.41 Mб
Скачать

310

ВЗАИ М О Д ЕЙ СТВИ Е ДВУХ МОД

[ГЛ. XVI

генерации обеих мод на одной частоте. Область нелинейного

захвата частот мод имеет порядок

При увели-

чении разности |Qjv— Qp | захватный режим сменяется одномо­ довым режимом генерации при сильном пространственном пе­ рекрытии мод (см. (16.28)). В случае слабого пространственного перекрытия мод (выполнено условие (16.27)) при увеличении разности частот |fiw — Qp | наступает режим нестационарных двухмодовых биений.

3.Устойчивость режима генерации зависит как от конку­

ренции, взаимной синхронизации,

так

и от

деформации

(iinnnp/^n Ф 0)

(см. § 4). При малой

разности

частот мод

| Qjv— Пр| имеет

место захват моды

Р на частоту

генерирую­

щейся моды cojv- Область деформационного захвата имеет поря­ док (см. (16.45))

При |Qjv— Qp| < | Q/v — Пр|Гр предпороговая мода Р остается в захвате при любой накачке и генерируется совместно с мо­ дой N на частоте cojv- При | QNQP| > ] Qn £2р |Гр мода Р при достаточной накачке выходит из захвата и при слабой конку­ ренции с модой N устанавливается двухмодовый режим гене­ рации.

Ч А С Т Ь 3

ЕСТЕСТВЕННЫЕ ФЛУКТУАЦИИ ИЗЛУЧЕНИЯ В ЛАЗЕРАХ

Г Л А В А X V II

ЕСТЕСТВЕННЫЕ ФЛУКТУАЦИИ В ГАЗОВОМ ЛАЗЕРЕ, РАБОТАЮЩЕМ В РЕЖИМЕ БЕГУЩЕЙ ВОЛНЫ

§ 1. Введение

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

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

Естественные флуктуации обусловлены молекулярной приро­ дой рабочего вещества и стенок резонатора и поэтому принци­ пиально неустранимы. Они определяют предельную стабильность частоты генерации квантового оптического генератора, предель­ ную чувствительность лазерного гироскопа и т. д. Так, например, стабильность частоты генератора определяется отношением Асо/мо, где Дш — естественная ширина линии излучения лазера. Предельная чувствительность лазерного (оптического) гироскопа определяется шириной линии сигнала разностной частоты и временем наблюдения.

Естественные флуктуации имеют более широкий спектр, чем технические флуктуации (ширина спектра технических флуктуа­ ций порядка 103—104 гц). Это дает возможность выделить сла­ бые естественные флуктуации на фоне более сильных техниче­ ских флуктуаций. Расчет естественных флуктуаций представляет, конечно, не только практический, но и научный интерес, так как исследование флуктуаций позволяет получить значительную до­ полнительную информацию о процессах в лазерах.

312

ФЛУКТУАЦИИ В РЕЖ И М Е БЕГУЩ ЕЙ ВОЛНЫ

[ГЛ. XVII

Можно указать два источника естественных флуктуаций. Это, во-первых, тепловые флуктуации в резонаторе, не связанные с переходами между рабочими уровнями. Эти флуктуации можно считать равновесными и рассчитывать по формуле Каллена — Вельтона. Вторым источником являются флуктуации поляриза­ ции активной среды. Они обусловлены атомной структурой активной среды и спонтанными переходами между рабочими уровнями. Эти флуктуации являются, естественно, неравновес­ ными. Расчет их и составляет одну из основных задач теории естественных флуктуаций в лазерах.

Чтобы не затенять сущность явления, начнем рассмотрение естественных флуктуаций для лазера, работающего в режиме бегущей волны (кольцевой лазер при условии подавления одной из встречных волн).

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

икольцевого лазеров.

Внастоящее время имеется значительное число работ, по­ священных экспериментальному и теоретическому изучению

естественных флуктуаций излучения лазеров. При последова­ тельном квантовомеханическом описании флуктуаций излучения лазера возможны два эквивалентных подхода. В качестве ис­ ходных можно использовать уравнения для матрицы плотности всех переменных атомов и поля — квантовый аналог классиче­ ского уравнения Лиувилля для функции распределения всех пе­ ременных системы. Такой подход развит, например, в работах Лэмба и Скалли [15, 16], Казанцева и Сурдутовича [17]. При другом подходе в качестве исходной используется система опе­ раторных уравнений для операторной матрицы плотности ато­ мов и уравнения для операторов электромагнитного поля (Лэкс

[10], Хакен [26] и др.).

Однако, поскольку число фотонов в генерируемой моде даже у самого порога генерации велико (порядка 104), для описания электромагнитного поля можно использовать классические урав­ нения. Квантовый характер излучения атомов рабочей среды учитывается в уравнениях поля посредством введения слу­ чайных источников. При таком подходе одна из основных задач теории естественных флуктуаций состоит в расчете статистиче­ ски неравновесных характеристик этих случайных источ­ ников.

Проведение расчетов на основе такой полуклассической тео­ рии значительно проще и дает возможность производить рас­ четы естественных флуктуаций в более сложных случаях. При-

ИСХОДНЫ Е У РА ВН ЕН И Я

313

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

§ 2. Исходные уравнения

Расчет флуктуаций будем проводить на основе системы урав­ нений для элементов матрицы плотности ра (д), рb(v), pab{v), Pba(v) и уравнения поля:

( 1

+ ' ж ) р .= т А

А „ -

Р.А.) в - V. (Р. -

Р?).

(17-1)

( ж +

» ж ) р. “ ~ т А

а . -

P- А я ) Е ~ v. (р. -

р?) ■

(17-2)

Здесь ЕЮ— источник тепловых флуктуаций.

Вектор поляризации активных атомов Р связан с элементами

матрицы плотности соотношением

 

(17.6)

Уравнения

(17.1) —(17.6) отличаются от использованных ра­

нее уравнений

(2.10) —(2.13), (2.20) тем, что теперь ра, рь, раь,

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

Во введении уже отмечалось, что существуют два источника естественных флуктуаций в лазерах. Первый — это тепловые флуктуации, т. е. флуктуации, не связанные с переходами между рабочими уровнями а, Ь. Найдем выражение для спектральной плотности тепловых флуктуаций.

Заметим лишь прежде, что в режиме бегущей волны поле Е

(г, t) можно задать в виде

 

 

E{r,t) = \ [ 8 ( t ) e lk°r + к. с.],

g{t) = eEe-l ^+w. (17.7)

Здесь

Е, ср — медленно меняющиеся

случайные амплитуда и

фаза,

е^- единичный вектор вдоль вектора Е.

314

ФЛУКТУАЦИИ В РЕЖ И М Е БЕГУЩ ЕЙ ВОЛНЫ

(ГЛ. XVII

Из уравнения (17.5), аналогично тому как это сделано в § 4 гл. II, получим уравнение для комплексной функции $ (t)

,

©о d<5 . о ее

 

 

dF- +

"Q '^ T + a)^ =

 

 

 

= - - f J

eP (r, t) e",v dr +

(t), (17.8)

где Eb}(t) — компонента

Фурье поля E iT)( r , i ) .

 

§ 3. Спектральная плотность источника тепловых флуктуаций

В уравнение поля (17.5) введен дополнительный член оаоЕ(т)— источник тепловых флуктуаций. Чтобы определить спектральную плотность этого источника, поступим следующим образом.

Разложим функцию £<т) (г, t) в ряд Фурье

EM{ r , t ) = ^ьl Ek{t)e^.

Средняя энергия тепловых флуктуаций с учетом, этого раз­ ложения может быть представлена в виде

 

Е(т)!

V:

14Т)

V.

(17.9)

 

 

 

 

 

Здесь

 

 

 

 

 

4 Т)

V = 2hck (п +

, й =

[ех р -^ ----l]

 

 

 

 

 

 

— средняя тепловая энергия двух взаимно перпендикулярных осцилляторов поперечного поля. Отсюда получаем (в расчете на одно колебание с k — k0) выражение

| еЕ(к} |2 V — 4яйсо0 (й + -j) • ®о = ck0,

(17.10)

которое определяет одновременную корреляцию пространствен­

ных компонент Фурье eE{ul{t).

Чтобы найти соответствующую спектральную плотность ис­ точника, стоящего в правой части уравнения (17.8), т. е. функ­

цию (еЕ*?)*» используем известный прием.

Поскольку нас интересует источник тепловых флуктуаций, т. е. флуктуаций, не связанных с флуктуациями поляризации Р,

рассмотрим уравнение (17.5) при 6Р =

0. Из него следует урав­

нение для компонент б£Й1(оо)

 

 

6 £ * 0 (со) (соо — со2 i

=

©оеЕ£1 (со).

5 3]

ИСТОЧНИКИ ТЕПЛОВЫХ ФЛУКТУАЦИИ

315

Отсюда получаем выражение для спектральной плотности

(17.11)

K - c o 2)2 + ( ^ f ) 2-

Получим выражение для функции (б£|„)ш иным путем. Запи­ шем уравнение (17.8) при 6Р = 0 и при отсутствии источника £<т) в виде системы уравнений первого порядка

d b E ——лI I

___________

' w- All —1—fjj J

dt

= бп,

^ -

+ ^ 6 П + ш26Е==0.

 

 

 

Умножим оба уравнения на бE(t') и усредним. Учитывая, что

бЕ (t) бЕ (/') = (бЕ ЬЕ)Х (т = / — t'),

получим

-^(6Е 6Е )Т= (6П6Е)Т,

(17.12)

 

dx (6П б£)т + ^ (6П бЕ)х+ cOq(6£ бЕ)х= 0.

Эту систему решаем при начальных условиях

(6Е б£)т„0 = | Е%|2 = -у- Пщ(я +

т ) . (бП б£)г_0 = 0. (17.13)

Первое из них определяет равновесную плотность энергии при k = k0 в расчете на одну плоскополяризованную волну, вто­ рое есть следствие равенства нулю среднего значения произведе­ ния обобщенных координат и скорости в один и тот же момент

времени.

Используем для решения уравнений (17.12) одностороннее

преобразование Фурье по т, например,

 

(б£ 6E)„ = | (6Е 6£)t eiax dx.

(17.14)

Из уравнений (17.12)

при начальных условиях (17.13) находим

(6Е6ЕГ

ШСОо

(17.15)

 

 

СОа — со

Из формулы (17.14) следует, что спектральная плотность равна

(6 £ 6 £ )a = 2 R e ( 6 £ 6 £ ) + .

(17.16)

316

ФЛУКТУАЦИИ В РЕЖ И М Е БЕГУЩ ЕЙ ВОЛНЫ

(ГЛ . XVII

Из (17.15), (17.16) находим искомое выражение

(17.17)

Сравнивая формулы (17.11), (17.17), находим выражение для спектральной плотности теплового источника флуктуаций поля

(17.18)

Из формул (17.18) видно, что спектральная плотность тепловых флуктуаций пропорциональна ширине полосы резонатора, кото­ рая определяется потерями, не связанными с поляризацией активной части рабочей среды.

Напомним, что величина щ/Q связана с эффективной прово­ димостью а (см. (2.31)). Величину а можно представить в виде суммы частей, характеризующих различные потери: за счет по­ глощения в зеркалах резонатора, выхода излучения, поглощения в неактивной среде лазера. В соответствии с этим спектральная плотность тепловых шумов может быть представлена в виде суммы соответствующих членов.

Естественно, что из-за различия температур зеркал и окру­ жающих тел функции п = ^ехр-^— lj , входящие, напри­

мер, в выражения для спектральных плотностей, обусловленных потерями в зеркалах и выходом излучения, будут различны. Это различие, однако, не является существенным в силу того, что для Не—Ne-лазеров отношение йсо/й7'>1. По этой причине мы будем использовать формулу (17.18), понимая в ней под Т некоторую среднюю температуру.

§

4. Уравнения для амплитуды и фазы излучения

в

гелий-неоновом лазере

Представим в уравнении поля (17.5) вектор поляризации

активных молекул в виде суммы двух частей

 

Р (г, t) = Р(инд) + 6Р(СП).

(17.19)

Здесь рсипд)(£) — индуцированная часть поляризации — отклик системы на поле Е = ( Е ) -j-6Е; 6Р<СП>— флуктуации поляриза­ ции, обусловленные атомной структурой активной среды. Эта часть вектора поляризации связана со спонтанным излучением (см. § 9 гл. XX). Конкретное выражение для спектральной плот­ ности поляризационного шума зависит от режима генерации.

§ 4]

У РА ВН ЕН И Я Д Л Я АМПЛИТУДЫ И ФАЗЫ

317

В настоящей главе будет рассмотрен лишь режим бегущей

волны.

Характер флуктуаций амплитуды и фазы излучения суще­ ственно зависит также от соотношения временных параметров Уа, уъ, уаь, Дсор. В Не—Ne-лазере обычно выполняются соотно­ шения

Уа> Уь> Yab>A“ p-

(17.20)

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

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

Р(инд) = (*' + Ы " ) Е .

(17.21)

Для режима бегущей волны действительная и мнимая части по­ ляризуемости активной среды определяются формулами (3.27).

Представим выражение (17.7) для поля в виде

 

E(r, t) = £ cos (со,/ — k0r + ф)

(17.22)

и вместо уравнения (17.8) запишем систему уравнений для мед­ ленно меняющихся амплитуды Е и фазы <р. С учетом равенств (17.19), (17.21) получим следующую систему уравнений:

+ х ( ^

+ 4ях'/) £ = (0‘> ^ ) ,

(17.23)

-g- + ^

- W = -f-g +(0.

(17.24)

Здесь введены обозначения для источников амплитудного и фа­ зового шума

=

у - J [ 4 я ( е 6 Р(сп)) +

( е £ (т)) ] sin ( с о ,/ —

kQr +

q > ) d r в g w - f | W

 

 

 

 

(17.25)

| ф =

— -p-1 [4я (е6Р(сп>) +

(еД(т))] cos (со0t -

k0r +

q>)dr^ | (фп) + $».

 

 

 

 

(17.26)

Из формул (17.25), (17.26) следует, что спектральные плот­ ности источников теплового шума определяются выражением

т . -

( « = } и и и +(•*$% +„].

(17.27)

 

318 ФЛУКТУАЦИИ В РЕЖИМЕ БЕГУЩЕЙ ВОЛНЫ [ГЛ. XVII

Правая часть

(17.27) в силу неравенства Дюр <С соо слабо зависит

от со, так как со <С соо. Поэтому при расчете флуктуаций ампли­

туды и фазы

достаточно знать спектральную плотность источ­

ника теплового шума на нулевой частоте (со =

0).

 

Из формул (17.18), (17.27) следует

 

 

йга=(isa-lw

t - f 4*'(s+?)•

<i7-28>

Из формул (17.25), (17.26) находим

 

 

(laTUT))o= (|фТ)|а })о= 0.

(17.29)

Это означает, что источники |1т),

статистически независимы.

Найдем соответствующие выражения для

источников

поля­

ризационного шума. Напомним, что вектор поляризации мы представили в виде суммы двух частей: индуцированной и спон­ танной (см. (17.19)). Проведем аналогичное разделение для элементов матрицы плотности р0, рь, раь

р = р(и«Д) +

6р(сп).

(17.30)

Из формулы (17.6) находим

 

 

р(инд) (Г) t) = П J

+ й аЪpj” rt) dr.

( 1 7 .3 1 )

Подставим выражения (17.30)

для элементов матрицы плот­

ности в уравнения (17.1) —(17.4). Уравнения для

индуциро­

ванных частей элементов матрицы плотности совпадают по форме с уравнениями (17.1) — (17.4) и соответственно с систе­ мой уравнений (2.10) —(2.13) гл. II.

Уравнения для элементов матрицы 6р(сп) при заданной на­

качке,

т.

е. без учета флуктуаций функций р^,

имеют вид

 

^■gf" 4

®

бра =

-g-(daj6pja

bPab^ba) Е

Yo^Pa»

(1 7 .3 2 )

( ж

+

V Ж

) брй =

г

бРЪа — ЬраьЛьа) Е — Yi 6рЬ,

( 1 7 .3 3 )

I

д

 

д

i®ab +

\

id .

 

( 1 7 .3 4 )

(^gj + V -g^T +

Yab) 6pa6 =

E (6pь — 6pa),

 

 

 

 

 

 

6p6o = 6p;6.

 

 

( 1 7 .3 5 )

(Здесь и ниже опускаем для сокращения записи индекс «сп».) Выражение для спонтанной части вектора поляризации свя­

зано с элементами матрицы 6р соотношением

бР (г, 0 = Я J (Аьа&РаЬ 4* брbadab) dv.

(17.36)

5 4]

УРАВН ЕН И Я Д Л Я АМ ПЛИТУДЫ И ФАЗЫ

319

Из формул (17.25), (17.26), (17.36) следуют выражения для

'-(п) 6(П)

источников поляризационного шума £W. ё<л> а *

 

X J {ed-ba6раЬ + tdab SPia) Sin (и</ ^0T+ ф) dr dv,

(17.37)

S(n):

4я . .

 

б*

— « X

 

 

X J (edba6pab+ edab 6pba) cos (co4 — k0r + qp) dr dv.

 

Введем медленно меняющиеся функции 6pab, 6pba:

 

 

6pa6= SPa6e‘ W - A”r+<P). fiP6a = 6P*a6-

07.38)

Подставляя выражения (17.38) в формулы (17.37) и выполняя интегрирование по г, получим выражения

£ (ап) (0 = 2 nni J (edba bpabedab bpba) dv,

 

£фп)(0 = — 2пп j (edba 6pab+ edab6pba) dv.

(17.39)

Отсюда следуют выражения для спектральных плотностей источников gw, g£>

a n . =

J

[ I « ь . p

»p«). -

 

 

 

 

— Иаб)2(6РбабРа&)ш +

К- ^]dv dv',

(17.40)

(&фП)’)а> =

4л2«2 J

[ I edba Г(6Pa6 6Pa>)a, +

 

 

 

 

+

(в<Га»)*(вРьавР«Д» +

K- C-] dv dv'-

(I7-41)

Таким образом, задача определения спектральных плот­

ностей

источников поляризационного шума gw, g^

сводится

к задаче нахождения спектральных плотностей (6pob6pab)a,

(foPbtfiPab)e>'

Из уравнений (17.32) —(17.35) для флуктуации бр следует,

что при однородном уширении линии ширина

спектра (брбр)<»

порядка уаь, а при неоднородном уширении,

когда ku уаЬ,

она порядка ku. Ширина спектра флуктуаций поля в Не—Ne-ла- зере порядка Дшр, т. е. много меньше (см. неравенства (17.20))