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

книги из ГПНТБ / Маркузе, Д. Оптические волноводы

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

340

Глава 7

Квадратичные среды обладают тем свойством, что

эффективный центр

мощности параксиальных пучков

движется в них по законам оптических лучей (см. разд. 3.6). Геометрооптический метод исследования таких структур физически оправдан. Он дает основные представления о свойствах распространения световых пучков.

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

До тех

пор

пока световые

лучи не достигают области

с /г < 1

или пока поле внутри этой области мало, можно

использовать

квадратичный

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

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

Среда с квадратичным законом изменения показателя преломления (7.1.1) может рассматриваться как первое приближение для среды с симметричным законом распре­ деления показателя преломления. Если разложить произ­ вольное распределение показателя преломления в ряд Тейлора, то можно получить члены порядка выше второго. Но во многих практических случаях ряды Тейлора схо­ дятся быстро, так что первых двух членов достаточно для первого приближения. Можно еще считать, что коэффи­ циенты при х2 и г/2 разные. Однако поскольку такое обобщение не дает больших преимуществ, то можно огра­ ничиться распределением показателя преломления вида (7.1.1) . Как было показано, среда с законом изменения показателя преломления 1/cli х является идеальной фоку­ сирующей структурой даже для непараксиальных лучей

Распространение света в квадратичных средах

341

[75, 118, 119]. Однако здесь мы ограничимся рассмотре­ нием только квадратичной среды и лучей в ней в пара­ ксиальном приближении.

7.2. ЛУЧЕВАЯ ОПТИКА КВАДРАТИЧНОЙ СРЕДЫ

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

rn+i

гн __ dr

 

(7.2.1)

L

dz

 

 

и

 

 

 

('n+i — rn)/Zi — ( r n r n _ i ) / L _

d Pr

(7.2.2)

L

 

dz2 '

 

 

Разностное уравнение (5.2.4) можно, таким образом, запи­ сать в виде

d2r __

1

(7.2.3)

dz2 ~ ~ ~ Lf

 

Допустим теперь, что /

оо, так что

 

1

Tli

(7.2.4)

Lf

п0

 

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

d2r

7*!

(7.2.5)

dz2

щ

 

для квадратичной среды. Дифференциальное уравнение (7.2.5), по-видимому, идентично параксиальному диффе­ ренциальному уравнению (3.5.36), в котором выражение (7.1.1) включено в правую часть. Таким образом, можно сделать вывод, что линзовый волновод становится квадра­ тичной средой в случае бесконечно близко расположенных линз. Решение для траектории луча линзового волновода может быть преобразовано подобным же образом к реше­ нию для случая квадратичной среды. Из формулы (5.2.6) получим cos 0 = 1 в пределе при L 0 и / -> оо, и,

342

Глава 7

 

кроме того,

 

 

si„ e = e = / f = / - a - i .

(7.2.6)

Используя выражение

 

 

 

z = пЬ,

(7.2.7)

получим из формулы

(6.4.26)

 

г = г °cos V w

z+ V i t sinV i t z

(7.2.8)

Очевидно, что (7.2.8) является решением уравнения пара­ ксиальных лучей (7.2.5). Положение луча и наклон1) г( выражены здесь соответственно через г0 и г', чтобы

показать,

что

они являются функциями

положения

и наклона

при

z = 0. Приведенный вывод

траектории

луча в квадратичной среде подчеркивает близкую связь между такой средой и линзовым волноводом. Разумеется, уравнение (7.2.5) и его решение (7.2.8) могут быть полу­ чены из уравнения параксиальных лучей (3.5.36) и фор­

мулы (7.1.1)

для показателя преломления без

ссылки

на линзовый

волновод. Итак, можно сделать

вывод,

что среда с квадратичным законом изменения показателя преломления ведет себя как световод.

Можно найти также траекторию луча в квадратичной среде с искривленной осью, используя результаты, полу­

ченные для

линзового

волновода.

Из

формул (5.4.8)

и (5.4.9) с помощью (7.2.6) и (7.2.7)

и с учетом того, что

dz = L при L -v 0,

получим

 

 

г (z)= r0 cos f / ^ L 2 + р Л К г’ Sin ] / ^ Z

+

 

__

2

___

 

+

/ V

1 W

si" /

 

<7-2-9)

 

 

0

 

 

 

Решение однородного уравнения [первые два члена (7.2.9)] записано в форме (7.2.8). Функция R (и) является радиусом кривизны изгиба оси волновода, образованного квадратичной средой, при z = и.)*

*) Величина г' является тангенсом угла наклона траектории,—

Прим. ред.

Распространение света а квадратичных средах

343

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

р (2)= Ро cos У ^ Z+ У

р; sin У ^ Z+

__

2

__

+ ] / V

I 5 (U) s in ] / -^ - (z -u )d u . (7.2.10)

 

о

 

Функция S (z) описывает отклонение оси квадратичной среды от прямой линии.

Эквивалентность двух решений (7.2.9) и (7.2.10) для траектории луча в квадратичной среде с изогнутой осью может быть легко продемонстрирована. Радиус кривизны изгиба оси квадратичной среды связан с линейным откло­ нением осп от идеальной прямой S (z) следующим при­ ближенным соотношением:

1 _ <г25

(г) — _ ; dz2 •

Это приближение справедливо, если первая производная dS/dz много меньше единицы. Отрицательный знак в (7.2.11) соответствует положительному радиусу кри­ визны на фиг. 5.4.2. Подстановка (7.2.11) в (7.2.9) при­ водит после интегрирования по частям к результату

Г (! ) = г , c o s ] / - Ь - * 4

- . 5 ( 1 ) 4

+ s < ° ) c o s / ^ + / l 4 4 f ) „ „ 5in / - £ г+

Ч

f 5 (“) s i n ] / - ^ - ( z —lt)

(7-2.12)

 

о

 

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

с формулой (5.4.10) имеем

соотношение

 

р (z) = г (z) + S (z).

(7.2.13)

Отсюда видно, что выражение (7.2.12) идентично (7.2.10), что показывает эквивалентность этих двух решений.

Среднеквадратичное отклонение луча а, вызванное случайными отклонениями оси квадратичной среды, по­

344

Глава 7

лучим из формулы (5.5.36). Это выражение является аппроксимацией более общего выражения (5.5.35). Однако в большинстве практических случаев это приближение хорошо оправдывается и удобно для использования. В пределе для очень малых промежутков между линзами выражение (5.5.36) принимает вид

оо__

a' M = T l t Z J Дс(и)С0 8 ] / T ^ udu- (7 -2Л4)

— ОО

Функция корреляции смещения оптической оси S (г) определяется как среднее по ансамблю произведение

*S (2) и S (z и)

R c (и) = (S (z) S (z - и)).

(7.2.15)

Выражение (7.2.14) показывает, что варпанс отклонения пучка пропорционален расстоянию вдоль волновода, на котором пучок еще наблюдается. Кроме того, оп пропор­ ционален преобразованию Фурье корреляционной функ­ ции. Преобразование Фурье корреляционной функции содержит только члены с косинусами, так как Rc (и) — четная функция.

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

7.3. МОДЫ КВАДРАТИЧНОЙ СРЕДЫ

Рассмотрим волноводные свойства квадратичной среды с точки зрения волновой оптики. Строго говоря, нужно найти решение уравнений Максвелла в среде с показателем преломления (7.1.1). Однако, как видно из разд. 1.3, скалярному волновому уравнению (1.3.6) должна удовле­ творять каждая составляющая электромагнитного поля. Для оптической среды с переменным показателем преломления волновое уравнение не является точным эквивален­ том уравнений Максвелла. Приближение, которое имеет место при использовании волнового уравнения вместо

Распространение свет.- в квадратичных средах

345

уравнений Максвелла, является приемлемым, если пока­ затель преломления (или диэлектрическая проницаемость) меняется незначительно на расстоянии порядка длины оптической волны. Это требование означает, что величина R (1.3.28) должна быть много меньше единицы. Такое требование с физической точки зрения также естественно, как и параксиальное приближение, которое было исполь­ зовано при решении задачи о распространении луча в квадратичной среде. Если нас удовлетворяют решения, которые вытекают из параксиального приближения, то для определения свойств распространения волны в квад­ ратичной среде вместо уравнений Максвелла можно использовать волновое уравнение. Таким образом, нужно рассмотреть распространение скалярных волн в квадра­ тичной среде [68, 69]. Скалярная волновая задача является хорошим приближением к задаче об электромагнитном поле, пока не рассматриваются поляризационные эффекты. Известно, что каждая составляющая электромагнитного поля должна быть приближенным решением скалярного

волнового уравнения. Скалярная волновая

задача в ква­

дратичной

среде

применима также к

проблеме распро­

странения

звука

в океане. В океане

имеются слои

с такими

изменениями плотности,

которые приводят

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

Предположив, что волны имеют гармонический во вре­ мени характер, можно использовать приведенное волновое

уравнение в форме

 

 

"

«2/^ф = 0

(7.3.1)

с постоянной распространения в свободном пространстве

7с0= с о / ^ Г о = ^

(7.3.2)

и относительной диэлектрической проницаемостью (ква­ драт показателя преломления)

n2 = nl — n0nl (x2-{-if).

(7.3.3)

Если считать исходным показатель преломления (7.1.1), то (7.3.3) является приближенным соотношением. С дру­ гой стороны, можно также считать выражение (7-3.3)

346

Глава 7

точным и рассматривать (7.1.1) как некоторое приближение. Такой подход имеет то преимущество, что п2, а ие п являет­ ся в данном случае основной величиной. Попытаемся решить приведенное волновое уравнение с помощью пробной функции

Ф = / (я) g (}/) e~ipz-

(7.3.4)

Здесь не известны не только функции / и g, но также постоянная распространения (3. Решение вида (7.3.4) есть мода квадратичной среды. Подставляя выражения

(7.3.3) и (7.3.4) в (7.3.1), получаем

( у - - S r + < К - Ра - Ц щ Щх 2 ) !

- | - ( y - | - - W ) - 0 .

(7.3.5)

Это уравнение должно быть справедливым для всех значений х и у, что возможно лишь тогда, когда зависимая от х и зависимая от у части порознь равны константе. Сумма этих двух постоянных должна равняться нулю. Таким образом, получаем два уравнения

+ (КК — к2 — Р'2 — к;п0щх2) / = 0

(7.3.6)

и

 

- g - + (x2- / ^ , z / 2)ir==0.

(7.3.7)

Чтобы преобразовать эти два уравнения к хорошо изве­ стной стандартной форме, введем следующие новые пере­ менные и параметры:

\ = V h ( n 0n{)Uk X,

(7.3.8)

г1=1А*о (п0щ)и!,у,

(7.3.9)

___

(7.3.10)

k0 У п0щ

 

__ к2

(7.3.11)

ко ~Vnoni

 

Дифференциальные уравнения запишем теперь в виде

^ + ( a - l 2) f = 0

(7.3.12)

Распространение света в квадратичных средах

347

И

 

^ r + (p - r)2)g = 0.

(7.3.13)

Дифференциальные уравнения (7.3.12) и (7.3.13) хорошо исследованы в теории гармонического осциллятора в квантовой механике [22]. Они определяют энергетическую собственную функцию и собственные значения, которые допустимы для гармонического осциллятора. В нашей теории, так же как и в теории гармонического осциллятора в квантовой механике, необходимо потребовать выполне­ ния граничных условий. Моды в квадратичной среде долж­ ны распространяться вблизи оси структуры. Это значит, что нужно потребовать, чтобы функции стремились к нулю при х оо и у оо. Такие же требования накладываются иа волновую функцию гармонического осциллятора. Они играют роль граничных условий. В работах по квантовой механике показано, что ограниченные решения такого типа могут существовать только при следующих условиях:

а =

+

1,

где

р = 0,

1,

2,

3, . . .,

(7.3.14)

о =

2q +

1,

где

q = 0,

1,

2,

3, . . ..

(7.3.15)

Эти условия вместе с (7.3.10) и (7.3.11) определяют воз­ можные значения постоянной распространения [5. Они, следовательно, являются собственными значениями урав­ нений мод квадратичной среды. Мы еще вернемся к этому важному понятию. Решения дифференциальных уравне­ ний (7.3.12) и (7.3.13) представляют собой хорошо изве­ стные функции Эрмита — Гаусса, которые уже встре­ чались нам неоднократно. Чтобы показать это, положим

и

/ = # ,( £ ) e-was*

 

(7.3.16)

g = H q(r])e-W№.

 

(7.3.17)

 

 

Подстановка в (7.3.12) и (7.3.13) дает

 

 

(PHг

 

ш

 

(7.3.18)

dl2

■2£

d\Р + 2 рНр (|) =

0

и

■2г)

9 '■2?Я,(т,) =

0.

(7.3.19)

<РНП

 

 

d,H

 

 

dip

 

dr)

 

 

348

Глава 7

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

l / 2 /я

Фру (я, У, 2):

~|/2Р+Чр\ q\ w

X e-№+v2Mw2e

с полушириной модового пучка

о

----

2

W~=

 

к0 V п0щ '

х

(7.3.20)

(7.3.21)

Модовая функция (7.3.20) здесь надлежащим образом нормирована.

Выражение для постоянной распространения $р<1 полу­ чается из формул (7.3.10), (7.3.11), (7.3.14) и (7.3.15):

Pp3 = K A " - * o /» o n 1[(2p+l)-l-(2g + l)]}1/2. (7.3.22)

Для сравнения с модами линзового волновода положим в формулах (6.4.47) и (6.4.48) L 0 и / ->- оо. Используя выражение (7.3.4), действительно получим полуширину мод квадратичной среды (7.3.21). При сравнении необхо­ димо помнить, что показатель преломления квадратичной среды на оси есть п0. Это означает, что в формулах (6.4.47) и (6.4.48) нужно использовать X /п0к0. Сравнение формул (7.3.20) и (6.3.20) показывает, что моды квадра­ тичной среды действительно тождественны модам линзового волновода c L - > 0 . Для полноты сравнения.покажем, что фазовый множитель в (6.3.20) переходит в (7.3.22) при отсутствии деления среды па линзы. Прежде всего по указанной выше причине необходимо в (6.3.20) к заме­

нить на и0к. Положим также

m = p i z n

= q n B пределе

при L —>■0 (z < L) получим

из (6.3.20)

с учетом (7.3.21)

рь=тг0/с0- 1 ] / l ± ( 2 p + i + 2q+l). (7.3.23)

Распространение света в квадратичных средах

349

Согласно формуле (6.6.10), радиус фазового фронта при d-i = L 0 становится бесконечным. Две постоянные распространения (7.3.22) иг (7.3.23), очевидно, совпадают только в первом приближении. Величину (7.3.23) можно получить из (7.3.22), если разложить в ряд квадратный корень и сохранить только первые два члена разложения. Из выражения (7.3.23) видно, что в первом приближении и при пренебрежении дисперсией среды моды квадратичной среды имеют замечательное свойство — их групповая скорость (др/<9со)-1 не зависит от порядка моды.

Таким образом, мы получили интересный результат: моды квадратичной среды совпадают с модами линзового волновода с точностью до фазового множителя в пределе при уменьшении до нуля расстояний между линзами. Фазовые множители совпадают только в первом прибли­ жении. Этот результат обусловлен тем, что моды квадра­ тичной среды (7.3.20) — (7.3.22) являются точными реше­ ниями волнового уравнения (7.3.1) с показателем прелом­ ления (7.3.3) . Ранее же мы отмечали, что моды Эрмита — Гаусса (6.3.20) являются только приближенным реше­ нием волнового уравнения. Приближение вносится тогда, когда вторая производная от и по z в (6.3.3) не учитывается. Таким образом, решение не является строгим, если оно получается путем перехода к пределу при стремлении расстояний между линзами к нулю. Однако в большинстве практических случаев величина

(?=

J A ot4

(7.3.24)

4

"5*0

 

обычно очень мала, так что величины (7.3.22) и (7.3.23) находятся в хорошем согласии, по крайней мере для малых модовых чисел. По и настоящая теория оказывается несправедливой, когда величина Q (7.3.24) достигает единицы. Скалярное волновое уравнение можно исполь­ зовать только в том случае, когда показатель преломле­ ния изменяется очень мало на расстоянии порядка длины волны. Это требование также приводит к Q 1. Чтобы убедиться в этом, образуем в соответствии с формулой (1.3.28) и с учетом (7.3.3) выражение

j q __пб п~“ ___ 1 nQiijX-

g 25^

‘о

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