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

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

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

180

УРАВНЕНИЯ МНОГОМОДОВОЙ ГЕНЕРАЦИИ

[ГЛ. XI

§ 9. Влияние поперечного распределения поля на нелинейную поляризацию активной среды

Рассчитаем поляризацию инверсно-населенного перехода с учетом поперечного распределения поля. В отличие от преды­ дущих параграфов номер моды будем обозначать прописными латинскими буквами (N,K,...). Каждая мода N характери­ зуется совокупностью трех индексов — продольным qN и двумя поперечными mN и nN. Вклад поляризации первого порядка по полю Pjv(0 в моду N в силу ортонормированности собственных функций определяется нолем самой N-u моды:

2ясо

АсОрА^о

 

е P{N](t)

2Wпор MNENe~‘(*«*+*«) + к. с.

(11.70)

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

,3) _2я®

jVoAcop

;

i+

(t))

 

г »

е

2Nnop

е

 

 

 

*

 

 

 

X

^

o,EjEKE Q\i.JKQNAjKQe ^Jkqn _)_ к. с.,

(11.70а)

где

(Pjkqn (со, +

 

— со0 — Мд,)/ + Фу + Фх — qp0 — q>N.

Коэф­

фициент \.i j k g n

характеризует перекрытие собственных функций

резонатора Еа(г)

 

 

 

 

 

 

 

V j k g n

=

j

E,(r)EK(r)Eb(r)E*N(r)dV,

(11.71)

 

 

 

 

 

Ко

 

 

 

где

Ко — объем трубки

с

активной средой; коэффициент A j k g

не зависит от распределений полей мод, а определяется только

их частотами со/, сок,

(см. (11.65), (11.67)):

Ajkg

WjKQ, волны бегут в одном направлении,

HJKa,

волна 7 бежит

навстречу волнам К. и G.

В формулы (11.70а) входят коэффициенты

WN N N

N>

w N K K : :Enk

(КФЩг HnKK^L N K -

*) При расчете поляризации движущихся атомов нужно учитывать из­

менения поля Еа(г), действующего на атом, со временем вследствие движе­ ния атома г' = г + vt. Мы будем пренебрегать изменением поперечной коор­

динаты атома, что справедливо, если расстояние иу-1, проходимое атомом за

время жизни у”1, мало по сравнению с масштабом рх, у изменения поля в по­

перечном направлении.

% 10] УРАВНЕНИЯ МНОГОМОДОВОЙ ГЕНЕРАЦИИ 181

В формулах для поляризации среды (11.70), (11.70а) коэффи­ циенты MN, RN, Fnk и Wjkg при однородном уширении контура усиления даются формулами (11.60), в которых нужно учесть, что скорость атома v = 0. В случае газовой среды с широким доплеровским контуром усиления коэффициенты M n , R n , F n k , L n k , W j k g и Hjkg даются формулами (11.66).

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

§10. Уравнения многомодовой генерации в газовом кольцевом лазере

Для того чтобы получить уравнения генерации газового ла­ зера, нужно в уравнения поля (11.43) подставить выражение для поляризации активной среды на N-м типе колебаний (11.70) и (11.70а). В приближении медленно меняющихся амплитуд и

фаз поля 8 Ns(t) = ENs(t) e~lWNt

1 dENs

■“C CONs,

d<pNs

■C ®Ns

(11.72)

' Ns

dt

dt

 

 

 

 

 

получим укороченные уравнения генерации для комплексных

амплитуд бегущих

волн

ENs(t) =

ENs(t)e~i<fNsii) (N — индекс

моды,

s — индекс направления волны,

s = l , 2)

dENs

+

[ 4

*

 

 

 

 

 

P+ N

dt

■ г'(со ^

^ л г ) j E p js

g

 

 

 

 

2

^ ns p E p s 4 "

 

4" 2 S ^ n p ^ ps' "b i

Д©N

Nn

 

{

 

 

2

Nnop

 

 

 

 

 

 

 

 

 

 

 

+

S

V-NP&NPaE2ps' +

S

tXNPF NPaE2ps ] E Ns +

 

 

P

 

 

 

P * N

 

 

J

 

4 *

2

 

\l p M ] N ^ p /M a E p s E J s / E M s,e

ss 4 *

 

 

P < J, M

 

 

 

 

 

 

 

 

 

 

 

4 -

2

 

 

 

 

( n -7 3 >

 

 

 

 

P < J ,M

 

1,2;

 

 

1

 

 

 

 

 

(s, s '=

s ф s'),

где

 

— ширина линии резонатора для N-й моды,

 

 

 

 

Ф и ' =

(fflps ”f- COys'

COMs' — COjVs) tt

 

 

 

 

=

(C0p^

COys

C°Ms

’ <%s) E

182

УРАВНЕНИЯ МНОГОМОДОВОЙ

ГЕНЕРАЦИИ

[ГЛ. XI

В сумме 2

отсутствуют члены,

у которых одновременно

р < / , м

равны индексы P = N и ] = М. Коэффициент нелинейного про­ странственного перекрытия V>PMjN= V 1pmjn==V-pmjn (п -71) опре­ деляется перекрытием собственных функций резонатора (см.

(14.22) — (14.45):

M-aip = n p p n J IFn (r) |2! EP (r) |2 dV,

Комплексные коэффициенты усиления, насыщения и нелиней­

ного взаимодействия

бегущих волн через активную среду MN,

R n , S ’n p , F n p , H p j m

и W p j m в случае доплеровского уширения

линии усиления определяются формулами (11.66), (11.67). В случае однородного уширения линии усиления нелинейное взаимодействие встречных волн такое же, как и взаимодейст­

вие

однонаправленных. Поэтому

в этом

случае коэффициенты

S ’n p

и Fnp, H p j m и W p j m попарно совпадают. Выражения для

коэффициентов даны

формулами

(11.60),

в которых нужно по­

ложить скорость v =

0.

 

 

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

(11.74)

где

В приближении плоских волн p,jv = 1.

Если взаимодействуют лишь две встречные волны одной моды, то в приближении плоских волн уравнение (11.73) совпа­ дает с уравнениями (4.4), (4.5). При взаимодействии четырех и более волн в уравнениях (11.73) имеется не равный нулю член

(при М = N и Р =

/), описывающий взаимодействие пары од­

нонаправленных волн с парой встречных

волн

(см.

(11.74)),

Этот член осциллирует с частотой

 

 

 

d<t> ,

® N s ~ 2 ((Op

<В;у)

2 (Qp

Qyy)>

^ == ®Р$ ”Ь

где йр — Qn — разность частот двух продольных мод резонатора

(Qp — Qn =

2nc/L).Такие члены вызывают осцилляции амплитуд

с той же

частотой, т. е. EN =

+ £ $ cos [2 (Qp — QN)t + ф].

§ ю]

УРАВНЕНИЯ МНОГОМОДОВОЙ ГЕНЕРАЦИИ

183

Оценим глубину модуляции Ем/Е$( и з амплитудного урав­ нения (11.73). Из уравнения нулевого приближения получаем, предполагая, что стационарные амплитуды мод равны:

(е % Т =(е№ =

где

55:5kM, N М,

%n m >

Рлгм>

kN kM.

Тогда для амплитуды осцилляций Е$ получим из уравнения первого приближения следующую оценку:

*9

1а ( Е (0)) 2Н рма

 

 

 

 

 

 

2i (<ор -

<в„)

 

 

 

 

 

 

(V nop) A(J\v

 

 

 

 

 

 

 

 

 

 

 

HpMG^N

 

 

 

 

4/a# + У] Ь N M

2i (С0р - <DW

 

 

 

Ил,—

 

 

 

 

 

м

 

(V nop) Д®.

Из формулы для коэффициентов Нрма (11.66), (11.67)

видно,

что \ НРма\^.2а. Отсюда в силу условия

 

 

Д с о N n

%

I ®р

млг I

(11.75)

 

2

~fj

получим

*

' v пор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ИЛ1___________

 

^

(Упор) % A(0N <

1.

 

a p - < » N

V

| _ а

 

4 ( ® P “ m Af)

 

 

(N0/Nuop)A“J

%

 

 

(11.75а)

 

 

 

 

 

 

 

Таким образом, в силу условия (11.75), означающего, что ха­ рактерное время изменения поля 1/(AtojvT)) гораздо больше пе­ риода осцилляции с разностной частотой TPN = 2n/\Qp— Qjv|, амплитуда этих осцилляций мала по сравнению со средней амплитудой поля:

^

ДЦдгИлг (NolNnop) < J

Е®

| £2Р —&N|

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

184 УРАВНЕНИЯ МНОГОМОДОВОЙ ГЕНЕРАЦИИ ГГЛ. XI

систему

 

амплитудных

и фазовых

уравнений

для

усредненных

за время

ГР* = ■

2к

 

1

величин

Ем, со*

Qp ~

q n

<

11 Дсо.

dt

^

f- f -

- i (CD* -

Q *)JENs=

/

/ nop

X

 

d E Ns

 

 

Aco ,r

 

 

,1

 

 

Д о .,

N n

 

 

 

 

Г АММ

 

 

 

 

 

 

 

 

 

 

 

X |

["л*

a Na ^ 2Ns

2

^ N P a ^ P s '

^

^ N P a ^ P s

Ens

 

 

l L

 

 

P

 

 

 

P ^ N

 

 

 

2

 

 

a (P p j m co s Фp j m n +

I p j m s in Фрумлг) E ps E js ^ M s ^P JM N | •

p <if, M

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

(11.76)

rf(PNs

+ <%s

&Ns

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

4" “-

J . „ + i A

2

Nnop

+

(IpjM C0S ФРУМЛС

P<J, m

+ S > A + S V ^ +

РФ N

P/VM З^Пфруд!*) - PS JS MS >PJMN "Ns

Г Л А В А

X I I

 

 

ВЗАИМОДЕЙСТВИЕ ПРОДОЛЬНЫХ

МОД

 

§ 1. Введение

 

 

Оптический резонатор имеет

бесконечный набор частот про­

дольных мод. В соответствии

с этим уравнения

генерации

(11.76)

представляют собой бесконечную зацепляющуюся си­

стему уравнений. Однако обычно условия генерации

(эффектив­

ное усиление больше потерь) выполнены только для небольшого числа п-мод.

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

Таким образом, задача об устойчивости n-модового режима стационарной генерации распадается на две: исследование внут­ ренней и внешней устойчивости. На основе системы Ап времен­ ных связанных уравнений для амплитуд и фаз бегущих встреч­ ных волн генерируемых мод решается задача о внутренней устойчивости 2п-волнового стационарного режима генерации относительно гашения одной или нескольких волн или возникно­ вения автомодуляционного (пичкового) режима. На основе временных уравнений для бегущих волн одной предпороговой моды решается задача о внешней устойчивости n-модового ре­ жима по отношению к возникновению новых генерируемых мод. Так как взаимодействие мод носит характер конкуренции, то может возникнуть генерация только такой предпороговой моды, у которой линейное усиление выше потерь.

Учет комбинационного взаимодействия существенно меняет дело. Если волны с волновыми векторами k и km бегут в одном направлении, то волновой вектор комбинационного тона 2k km равен й_т — волновому вектору волны, бегущей в том же

186

ВЗАИМОДЕЙСТВИЕ ПРОДОЛЬНЫХ МОЛ

[ГЛ. ХИ

направлении

и принадлежащей моде —т с частотой

о)~т =

— 2оз т . Моды т и —т, расположенные симметрично отно­ сительно моды с частотой со (сот — со — со — ш_т ), оказываются связанными друг с другом. На языке квантовой теории излуче­ ния комбинационное взаимодействие означает, что в активной среде при поглощении двух квантов частоты со возможно излу­ чение двух квантов с частотами <вто и co-m. Возможен и обратный процесс.

Этот вывод не распространяется на встречные волны. В этом

случае

волновой вектор комбинационного

тона

на частоте

со-m =

2со — сот равен 2k + km = 2km+ fe-m,

т. e.

амплитуда

комбинационного тона осциллирует вдоль оси лазера с про­ странственной частотой 2km. При усреднении по длине активной среды I вклад такого комбинационного члена в моду —т имеет исчезающе малый порядок Х/(21).

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

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

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

§ 2. Стационарная генерация сильной бегущей волны

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

§ 2] ГЕНЕРАЦИЯ СИЛЬНОЙ БЕГУЩЕЙ ВОЛНЫ 187

интенсивность генерации, нужно заново решить уравнения, оп­ ределяющие поляризацию активной среды.

При условии равенства ширин уровней генерации \ а = уь уравнения для матрицы плотности (11.46) могут быть сведе­ ны к дифференциальному уравнению второго порядка для

поляризации Р =

с?(раь + р&а)

и уравнению

первого

порядка

для разности заселенностей N =

раа рьь-

 

 

 

Полная система уравнений имеет вид

 

 

 

d2E (г, О +

д д Е (г , 0

„2 д2Е (z, t)

•4я

d2U

 

dt2

dt

dz2

^ d

t 2

 

д2Р

Ъ л ^ + А Р - - ■ ~ ^ N E ( Z , t ) ,

( 12. 1)

dt2

dt

 

 

 

 

3N

4a(N-

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

Электрическое поле должно удовлетворять условию периодич

ности

( 12.2)

E(L, t ) = E(0, t ) ,

где L — периметр кольцевого резонатора. При выводе уравнений

(12.1) предполагается, что активная среда

целиком заполняет

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

Система (12.1) пригодна для описания генерации как при однородном, так и неоднородном уширении линии усиления. В случае неоднородного уширения частоты переходов отдель­ ных атомных систем сааь смещены относительно центра линии:

®аЬ == ®о + х.

(12.3)

Будем считать, что функция

распределения частот соаь имеет

гауссовскую форму с дисперсией Г

 

f(x) =

^ i — e-wn*.

 

' w

VnT

 

Для однородного уширения f(x) — 8(х). В общем случае поля­ ризация Р(х) и разность населенностей N являются характери­ стиками атомов, обладающих резонансной частотой (о0ь, в то время как в уравнение поля входит средняя поляризация среды (средний дипольный момент единицы объема)

00

 

U = | P{x)f(x)dx.

(12.4)

Наше рассмотрение не зависит от причины, вызывающей не­ однородное уширение. Для твердотельных лазеров такой

188

ВЗАИМОДЕЙСТВИЕ ПРОДОЛЬНЫХ МОД

[ГЛ. XII

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

Запишем решения системы (12.1) в виде

Е (г,

t) =

Re {Е (z, t) e~l

 

P(z, t, x) =

Re{— iP(z, t,

 

U (z, t) =

Re{— iU(z,

t)e~l(at~kz)},

(12.5)

 

 

oo

 

 

0(Z,

t ) =

J P(z, t, x) f (x) dx.

 

Здесь E(z, t), P(z,

t, x) и U (z, t)—

медленно

меняющиеся

комплексные амплитуды поля и поляризации.

Тогда в обычном

оптическом приближении (заключающемся в

том, что ширина

спектра частот и релаксационные константы

рассматриваемой

задачи много меньше самой частоты) система укороченных уравнений для медленных переменных, описывающая взаимо­ действие поля с доплеровским ансамблем атомов, имеет вид

дЕ

дЕ

Д(0п

 

Ж + С~Ш +

• I (со — Q J Е — 2n<oU,

 

 

 

 

Щ' + [Yab— I(о — ©Об)]Р — 1X

( 12.6)

d N + y a ( N - N 0) = : - ^ r Ф Ё ' + E h ,

 

dt

 

 

где Qq = ck — резонаторная частота.

Волновой

вектор

k = -j- q (q — большое

целое число),

и

медленно меняющееся поле E(z,t)

удовлетворяет условию

пе­

риодичности

(12.2).

 

 

 

 

 

Если амплитуды поля и поляризации постоянны во времени

и пространстве, то решения (12.5)

описывают стационарный ре­

жим генерации плоской бегущей

 

волны. Такие стационарные

величины определяются из уравнений (12.6):

 

 

 

М - *0 0 +*/«)•

,

— rf!Wo

£ст

 

 

°

Ьу“Ь ! + / + /£*’

 

 

 

 

 

tfoO + fo*)

 

 

 

(12.7)

 

 

 

со — co0

 

ЛГст

,2

*

 

 

 

 

fo =

vаь

 

У+ I + fo

 

 

 

 

§ 2] ГЕНЕРАЦИЯ СИЛЬНОЙ БЕГУЩЕЙ ВОЛНЫ 189

где / =

а£ст==^2—

---- безразмерная интенсивность

генерации.

Здесь

 

^аЬ^а

для определения средней поляризации

и в дальнейшем

среды

U по формуле

(12.4) нам понадобятся интегралы вида

 

 

1

Г

e~xW dx

2i

 

 

 

-L

a0+ (fo+ xf

Zi (Si» У ,

 

 

 

h + t2

( 12.8)

 

 

 

 

 

 

 

 

 

1

Г

(fo + x) e~x V

^ 2 (Si>У .

 

 

 

V n

 

a20 +

(f0 + x f

 

 

 

 

 

 

где ti =

ia0— f0,

t2 = ia0+ f0. При

условии Re {(a0 — if0) b) > 0

функции

Zj и Z2

следующим образом выражаются

через спе­

циальную плазменную функцию Z (£):

Zi(h,

h ) = i j - [ z m + z m ] ,

z 2(Si,

y = 4 t Z { b W - Z [ b l2)l

Напомним, что

 

 

il

 

y ir

Z (|) = 2i J e~^+^dt = 2ie-V

2

—oo

 

 

(12.8a)

e-t'dt .

Если a0, fo и b — вещественные величины, то g2 = — и, исполь­ зуя свойство плазменной функции Z*(£) = —Z(— |*), получим

Zi( y

- I I ) =

Im Z(6i1),

Z2(g„

- Q =

( 12.86)

Re Z(6Ej).

Стационарная поляризация среды f/CT находится по форму­ лам (12.4), (12.7), (12.8):

Ост =

[17f 17 Z, (£„ У + iZ2 (£„ у ] ;

(12.9)

функции Z] и Z2 определены формулами (12.86), где

12= - t ; = * Y T + T + f 0, 6 = Yaft/r .

В соответствии с (12.86) функции Zi и Z2 вещественны. Подставляя (12.9) в первое уравнение (12.6), получим