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

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

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

270

НЕЛИНЕЙНОЕ РАСЩЕПЛЕНИЕ ЧАСТОТ

[ГЛ. XV

где

N — основная мода 00q, р — четно-четная мода

(р = 2т,

2п).

 

 

•' В этом случае нелинейное дифракционное расщепление ча­ стот генерации встречных волн много больше области их син­ хронизации и при отсутствии других невзаимных эффектов разность частот генерации a»i — «2 полностью определяется нелинейным расщеплением частот встречных волн. В других условиях нелинейное расщепление частот вносит вклад в раз­ ность частот генерации о>х— сог аддитивно с другими линейными невзаимными эффектами, приводящими к различию частот гене­ рации встречных волн.1

Для того чтобы уравнения (15.11) описывали одномодовый режим (15.1), будем считать, что условия генерации выполнены только для самой добротной моды с индексом N

А ».,

 

г\<

Дсо0 Доз,

(15.12)

0 < — — —

(рФЩ.

^ 2

N.пор

 

 

 

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

Из уравнения для предпороговой моды определим комплекс­ ную амплитуду моды Ера (р Ф N). Из уравнений для генерируе­ мой моды найдем частоты генерации и вещественные ампли­ туды ENs.

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

Будем решать уравнения (15.11) методом последовательных приближений, при условии малости амплитуд предпороговых мод

(P¥*N).

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

г«2

=

г*2 __ р2

%

t

 

— £о-

<Ч*дг (а + Р) ’

 

 

 

 

 

До,

Nn

Р + т Даз^о

®ЛМ— ®N2— — £2N.

Nпор

-а +

а + р

 

 

§ 2]

ОДНОМОДОВЫЙ РЕЖИМ ГЕНЕРАЦИИ

271

Из уравнений (15.11) для предпороговых мод (р ф N) и условия (15.12) видно, что при пренебрежении связью с генери­ руемой модой амплитуды предпороговых мод равны нулю.

Первое приближение для предпороговой моды

(о», -

Ш ! — [ ? % - +

^

V

+ И ) а й )

где

 

 

 

(15.13)

Qnp = Qn Qp, Асо^р = со^

Ир

{рфЫ’ ).

 

Уравнения (15.13) являются неоднородными. В правой части (15.13) стоит вынуждающая сила, пропорциональная амплитуде генерируемой моды.

Из уравнений (15.13) следует

Ёр1 =

RpsNEo,

(15.14)

4pN ц «

 

Nn

 

(W + H) aE20

*^NNNp

о

N

 

R p N = m

 

[bap

A«v \

N ^ 2

2 /

Малость амплитуд предпороговых волн \Ера\ при линейной и нелинейной связи мод определяется условиями:

Дсоп

{рФЫ, 5 = 1 , 2).

 

- % - 1{

По терминологии теории колебаний эти неравенства являются условиями малой связанности мод N и р для многомодовой колебательной системы с затуханием.

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

m'pN= ~ тТы= mNp

(15.15)

Согласно формуле (14.23)

^AWJVp = ^JVAfp = V-pNNN-

(15.15а)

Из выражения (15-14) в силу соотношений (15.15) и (15.15а) получим

R p N

Ё - Ы р ‘

(15. 156)

272

н е л и н е й н о е р а с щ е п л е н и е ча с то т

[ГЛ. XV

Первое приближение для генерируемой моды. Из системы (15.11) с учетом выражений нулевого и первого приближений получим следующие уравнения:

<> +

т г ~ V

е , (v e 'S. +

« а д

= -

Е

-

 

 

 

 

 

b-jbN

 

 

Л'пор

S { ^WVAW

2

^

+ W7 Re (Мадууд/Ям) +

 

 

b+ N

 

 

 

 

 

 

 

 

 

+

Я Re

 

 

(s, s' = 1,

2,

s Ф s').

 

Решая эти уравнения, получим выражение для сдвига частоты генерации волны ©^

 

тЕ

Р + Т

ДсОддЛ

 

(дО) =

а + Р

2 /

+

S

 

 

 

ЬФЫ

 

 

 

(15.16)

С помощью (15.15) и (15.15а) получим, что со',1’= ©у. Таким

образом, в первом

приближении по коэффициентам распреде­

ления RSpN частоты

генерации ©0 + ©(1) одинаковы для встреч­

ных волн.

 

* Смещение частоты генерации ©(1) состоит из трех членов: линейного а*1», пропорционального интенсивности ©у* и пропор­

ционального квадрату интенсивности ©^>.'

Первый член в а*,1’ связан со смещением резонаторной ча­

стоты Qn из-за дифракции, второй определяет смещение частоты генерации, пропорциональное изменению интенсивности, возник­ шему при изменении резонаторн ix потерь.» Сдвиг частоты Qjy и ширины резонанса A©jv при учете связи мод могут быть объяс-

§ 2]

ОДНОМОДОВЫЙ РЕЖИМ ГЕНЕРАЦИИ

273

йены на основе формул (14.60), в которых параметр связи т2 нужно считать равным

- т 1 т &Р-

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

ЛСдвиг

частоты <4" линейно

зависит

от

коэффициентов

связи ’ть%.

Если электрические

векторы полей

мод b и N па­

раллельны,

т. е. (е*влг) = 1 и

коэффициенты

отражения

равны

гх = гу = г,

то согласно (11.28)

mfN можно представить

в виде

mt% = 2/v0£ (бтьт„6nbnN- г М Ъ ) е - 1*™ Ы .

(15.17)

Здесь Mbh — интеграл перекрытия полей мод Ь и N на l-м зер­

кале (см. (11.32), (11.33)); (fsbN(zl) =

(kb- k N) z i—(fb(zl) + (pN(zl)—

разность фаз

мод для волн, бегущих в направлении s в точке zt

нахождения

/-го зеркала (ф^ (z) =

—q>lbN(г));

п — число зеркал

кольцевого

резонатора.

и (14.22),

можно показать,

Используя

выражения (15.17)

что суммирование в формуле (15.16) нужно

проводить лишь

по модам, у которых хотя бы один поперечный индекс отличен

от

индексов

моды N. Для таких мод коэффициент DbNss

=

Im

« V ' b N N N

можно представить в виде

'bN

и*

 

 

 

(15.17а)

г'~Т

где

Ф&лг (z) = (kb — kN) z — (ть — mN)arctg ( J j) —

— («» — «#)arctg

коэффициент h0 определен формулой (14.28); z0= (2, -f z2)/2 — центр трубки; l — (z2— 2,) — длина трубки. Величина и знак

коэффициента D% зависят от расстояния между центром трубки

274

НЕЛИНЕЙНОЕ РАСЩЕПЛЕНИЕ ЧАСТОТ .

[ГЛ. XV

и диафрагмой;

при

\ zQ— гг I -j

3L п ц\

В при-

-4- D&jv меняет знак.

ближении короткой трубки / <С

 

 

Dw

( ^ ™ - ) cos\щы{z0) - q bN (zt)} МУЪ.

(15.176)

bN ■

Знак смещения

частоты

определяется спектром

частот

резонатора QbN и, согласно (15.176), расстоянием между диа­ фрагмой и центром среды |z0 — Zi\.

' Сдвиг частоты Ц1* не зависит от дифракции и определяется нелинейной деформацией поля волны в активной среде.* Коэф-

фициент деформации

VbNNN

 

abNNN

hbNNN

, Г Д е

flbNNN

VN

\

“V

hn

 

 

трубки

определен формулой (14.26).

В

пределе

короткой

1Ti&jvjviv| = hb (см. (14.33)).

Тот факт, что продольные моды с одинаковыми поперечными индексами не вносят вклада в ©0), означает, что неполное отра­ жение зеркал, связывающее только продольные моды с коэффи­ циентом связи т |^ ~ ( 1— г), не приводит к смещению частоты

генерации в первом приближении по rj.

Перейдем к определению разности частот генерации встреч­ ных волн. Согласно формулам (15.7а) и (15.9) разность coi — м2 является величиной второго порядка малости. При выполнении условий

I m ( W — Я )

PcNNN

 

Re («7 - Н)

Im | S № cN

 

 

 

I^cNNNI ^

(15.18)

 

 

 

выражение для

разности частот (15.7а) можно

представить

в виде

 

 

 

ДсОрЯ^! w ■ НI

 

©! — (02

( W - H )

 

^CiV

2 Im

b, с

 

 

(15.18a)

 

 

 

Из выражения (15.7а)

следует также, что ©i — ю2 =

О при W =

= Н. Стремление разности частот ©i — со2 к нулю наблюдалось на экспериментах [2, 3] при генерации вблизи центра линии уси­ ления одноизотопного Не—Ne-лазера на X = 0,63 мкм, где

W w Н (см. рис. 15.1).

§3]

ВЛИЯНИЕ ПАРАМЕТРОВ РЕЗОНАТОРА И СРЕДЫ

275

§ 3. Зависимость расщепления частот генерации встречных волн от параметров кольцевого резонатора

иактивной среды

Впредыдущем параграфе выведена формула (15.18а) для разности частот генерации встречных волн. В этом параграфе мы займемся детальным анализом этой формулы.

Подставив в (15.18а) выражения (15.14) для коэффициентов

распределения амплитуд./?^ (s = 1, 2) и значение аЕо, получим

 

 

(р — т)2 + (а — Р)2

Дш

 

IV,

 

 

 

со, — со2 =

а 2 — р2

 

/■ дГ-^лгХ

 

 

 

 

^

/v пор

 

 

 

х £

 

 

 

 

Q N c +

 

 

 

1 VoQfjb^NcPbc +

Ь, сфЫ {Qllb+ ('

 

(

A °> N c ) 2

 

 

 

 

 

Дш,.

N„

 

 

\

2

 

 

 

+

 

 

 

 

 

 

Р +

Т ^

2v0 -

0

г л

гг

( п е л

к

 

 

v° ~ 2

Nnop ^N^NcZcb \ 2Q jv6

 

а +

р / ~~

/ Д ю „ N

\2

Г

 

 

 

 

Р + Т

— 2 Naop

^Ьс Y^Nb^Nc — 4А(ВMc^Nb а

+

+ Аюдйс,^-^!-)2J, (15.19)

где

 

 

Im {v-bcNNmbHmlcN)

 

 

 

be '

 

 

 

 

HiVv0

 

 

 

 

 

 

 

 

Z-Cb

 

t1MJVJV6,ncJv)

(15.20)

 

 

 

ll>0

 

 

 

 

 

 

 

(-’be

,m (V -b cN N V N K N b P -N N N c)

 

 

 

 

llW

 

 

 

 

 

 

Pn ^ P nnnn»

 

с

ДШд.

Nn

 

 

= —r—7=

т

—— Лм — коэффициент усиления

моды /V;

=

L v e

2

Nnop

 

Qn — Йь; Дсолгь = Дшл— Асоь. Выражение (15.19)

справедливо при выполнении (15.18). Первое из условий (15.18) существенно лишь для лазера на чистом изотопе, где оно приоб­ ретает вид

Таким образом, формула (15.19) не описывает область вблизи центра линии для чистого изотопа. Подставив выражение для ml'fj (15.17) и выражения для цш1 (14.22) в формулы (15.20),

276

НЕЛИНЕЙНОЕ РАСЩЕПЛЕНИЕ ЧАСТОТ

[ГЛ. XV

найдем выражения для коэффициентов РЬс, Zcb и Lbc-

р ьс= —

77 ^

S

(йт&'плАьг,л — riMbh) X

 

'

N ' 0

<=1

А=1

 

 

 

 

 

 

 

 

X (PmcmN^ncnN

rkMc^j P[cn\

Zcb— '

b c N N a N N N b

2

n

 

 

 

~j2

(PmcmN ^ncnN

r l^ c N )

ZcbNt

 

 

 

 

 

 

1=1

A A

 

 

i

 

 

 

 

 

__ [ abcNN°NNNbuNNNc \

^bcff

(15.20a)

 

 

hn

 

4

J

$

 

 

 

 

Коэффициенты поперечного перекрытия мод abcdp даны в § 2

гл. XIV (см. формулы (14.25), (14.27). Интегралы PbcN, Zb%, LbcN и h определяются формулами

 

*о+4

 

 

 

 

 

{f£

 

P&N =

г

Sin fo»JV( Z ) <P&JV( Zl) +

 

 

(g)

J

ФcN (Z) VcN (Zk)] ~~-{z) *

Z^bN=

z0+4J

*»+4 sinJ[фсЛГ(z) —ycN (zt) -f- q)jiV(z) — (fbN(z')] X

 

 

 

 

 

 

X

 

dz dz'

 

 

 

 

 

 

Px ( z ) p y ( z ) p x ( г ' ) Py ( г ' ) •

 

 

 

 

 

 

 

z,+4

Sin [<pftA(2) -

 

(20 + <PcAf (Z) -

<fcN (Z")\

,

 

- bcN

Ш "

 

 

 

P x ( z ) p y ( z ) Px ( z ' ) Py ( z ' ) p x ( z " ) p y ( z " )

a Z

»

 

 

 

 

/

 

 

 

 

 

 

 

^o+“2

 

 

 

 

 

 

- J

d z

 

(15.21)

 

 

Px ( z )

Py ( z )

 

 

 

 

 

 

 

 

 

 

Выражения для масштабов px(z)

и py(z) и приведенных ра­

диусов Ьх, Ьу для кольцевого резонатора со сферическими зер­ калами даны в § 1 гл. XIV (формулы (14.5) —(14.19)).

Интеграл /0 пропорционален коэффициенту нелинейного про­ дольного перекрытия мод Л0 (14.28), I0 = L2hQ. Оценка инте­ грала Л0 произведена в формулах (14.29) —(14.34).

8 3]

ВЛИЯНИЕ ПАРАМЕТРОВ РЕЗОНАТОРА И СРЕДЫ

277

Если в кольцевой резонатор помещена короткая трубка с активной средой (такая, что длина активной среды I bXlby), то можно получить приближенные выражения для интегралов

(15.21)

pH. к)

“bcN

L cbN

s i n ( k b N + k cN) l / 2

(kbN kcN)

X sin [(p&iV(z0) —q>bN (zt) -f ФсЛг (Zo) — фсм (Zft)], (15.22)

sin ( k cN + k bN) 1/2

s'n (ЬЬы42)

 

sin[qpcjv (z0) — tfeN(z,)]

(bCN + * b N ) 4*

ЬЬц112

при / <

bx,

by

и l < -J

 

 

 

 

 

 

 

 

 

cN

 

1

 

 

 

bcN

 

1

1

+

 

 

 

 

I2\b*x +4z$

2 j ^0l*ki,ffkeif (k^N “b kcft),

где

 

 

 

by + 4zq

 

 

 

 

 

 

 

 

bx/2

 

 

by! 2

kbN === k b

k ff

•(mb — tnN)

 

(nb — nN)

 

 

(byllY + 4,2 ‘

 

 

 

 

 

( b x/ 2) 2 + z 2Q

 

 

Для кольцевых резонаторов с большим радиусом сферических

зеркал Rt^_L приведенные радиусы bx

 

w b v растут пропорцио­

нально

V R-

В

этом случае

условие

I

bx, Ьу выполняется,

даже когда активная среда распределена по всей длине кольце­ вого резонатора I т L .

В формуле для разности частот (15.19) стоит двойная сумма по всем модам, отличным от генерируемой моды N . Рассмотрим суммирование по продольным модам, имеющим такое же попе­ речное распределение, как и генерируемая мода N . Потери всех продольных мод одинаковы, т. е. А&т = 0; разность фаз таких мод равна

 

ФNb(z) = sNb^ ,

(15.23)

где

SNb = QNЩь пробегает целочисленные значения

от — оо

д о

+оо. Соответственно спектр частот эквидистантен, т.

е.

 

&Nb 2jtsAfjVo.

(15.24)

Из выражений (15.20а), (15.21) и (15.23) следует, что коэффи­ циенты Р ьс, Z cb и L bc меняют знак при одновременной замене

SbN—►—sbN и scjv—*■—sCN- Например,

P bc (s bN> s cn) = P b c i s bN> s cn)-

Учет этого свойства коэффициентов P bc, Z cb для спектра (15.24) приводит к тому, что

и L bc и выражения двойная сумма по

ш НЕЛИНЕЙНОЕ РАСЩЕПЛЕНИЕ ЧАСТОТ [ГЛ. XV

продольным модам в формуле (15.19) равна нулю. Таким обра­ зом, суммирование в формуле для разности частот (15.19) прово­ дится лишь по таким модам с и Ь, у которых по крайней мере один поперечный индекс отличен от поперечного индекса моды N. При этом продольные индексы мод с, b и N могут как совпадать, так и отличаться друг от друга. Этот факт означает, что в одно­ мерной модели плоских волн невозможно получить эффект расщепления частот, необходима трехмерная модель с учетом неоднородности поперечного распределения поля.

Выражения для коэффициентов Рсь и Zcb для различных по­ перечных мод можно упростить:

Pcb

/

abcNN

\

 

г М№ г

Р<1-*>

V

aN

)

u

' t m bN к

cNr bcN >

 

i= l k=l

(15 .2 5 )

 

 

 

 

 

 

7 -

( abcNNaN N Nb

)

т Ь .

 

t'cb

<4

 

 

 

{

 

/

I q <=i

 

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

(Mbjv=0), то Р сь 0 и Z cb = 0. Таким образом, расщепление частот за счет этих коэффициентов не происходит.

Из выражения (15.19) видно, что расщепление частот яв­ ляется нелинейным эффектом, причем в (15.19) входят члены, отвечающие разным физическим механизмам, с различной сте­ пенью зависимости от цм~аЕ1. Разность частот (15.19) можно

представить в виде

ю , — со2 = Ай,11лг+

+ ASgTft.

(15.26)

Здесь Доп содержит члены, пропорциональные Рсь, т. е. члены второго порядка по дифракционным коэффициентам связи мод. Член Дацц является результатом дифракционного искажения полей встречных волн. Дй2 ~ Z cb зависит линейно от дифрак­ ционных коэффициентов связи и линейно от коэффициентов aNNNb, определяющих нелинейную деформацию поля в активной среде. Член Дй2г]2 возникает в силу интерференции дифракцион­ ного изменения и нелинейной деформации поля, разной для встречных волн. Член ДйзТ^, пропорциональный Lbc, не зависит

от дифракции и имеет чисто нелинейный характер. Каждый из членов AcOjrpv, Доу!2,, Дю3т]^ имеет наинизший для своего эф­

фекта ПОРЯДОК ПО T)jv.

При соблюдении условия

N .

aNNNb

(15.27)

= \ м b N »

----%

aN NNN

T V 8

N пор

 

§ ]

ВЛИЯНИЕ ПАРАМЕТРОВ РЕЗОНАТОРА И СРЕДЫ

279

3

 

 

выполняются неравенства

 

 

 

>

Дй2т]^ > Дй3г1^.

(15.28)

Коэффициенты Мъп и

определяются выражениями (11.32),

(11.33)

и (14.276). Коэффициенты c i n n n n даны в

табл. 14.1.

Условие (15.27) означает, что линейная связь мод за счет ди­ фракции много больше нелинейной связи мод за счет деформа­ ции поля в активной среде.

Рассмотрим зависимость параметров выражения для рас­ щепления частот (15.19) от расстройки частоты генерации ю == = (coi-f м2)/2 относительно центра линии усиления озаь- В вы­ ражение (15.19) входят три коэффициента, зависящих от рас­ стройки и — соа&:

р _ _ ( р — т ) г + ( а — Р )2

^2

_ Р + т

~jjT~ “П* (15.29)

1 ~

а 2 — Р2

— а + р *

Два из них, Fi

и (N0/NnOp)r\, — четные функции

со — шаь, коэф­

фициент F2— нечетная функция со соаь.

Для чистого изотопа с широким доплеровским контуром уси­

ления (уаь •с ku,

| 0) — СОаЬ | <

ku)

 

 

 

 

 

Fi (со — ®ab) =

(м-

(0аЬ)2 + уа2Ь

 

 

 

 

 

(® -

®ай)2 + 2уаЬ2 ’

(15.30)

 

 

 

 

(О - даЬ) уаЬ

 

F2((О— &ab) =

 

 

(аз -

<оа 6) 2 +

2Yaь

 

 

 

 

 

 

В соответствии

с (11.67)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(15.30а)

где

 

 

 

 

 

 

 

 

Для 50%-ной смеси изотопов (уаЬ<

к

» Дсонз, | 0 — и I<.ku)

Nn

 

 

 

 

 

 

 

т пор - М0

Nnop ^ = ( ^ о р - |2) ( 1 - 2 е 2).

 

 

k u

=>пор

k u

 

62

_ ( N\ + Nn

Л

 

 

(15.31)

 

П°Р

I

Л/пор

)

1'

■2ег0

Дш^з

где ю0 — центр

суммарного

контура

усиления,

е0 — ■2k u ~'

Дсоиз— ®аы — ааЬ2— изотопический сдвиг (Дюиз > 0). Выражения для коэффициентов р, т, р, а даны формулами (3.54).