Документ взят из кэша поисковой машины. Адрес оригинального документа : http://hbar.phys.msu.ru/hbar/optmicro/lec10.pdf Дата изменения: Wed Sep 26 00:00:00 2007 Дата индексирования: Mon Oct 1 20:42:42 2012 Кодировка: Windows-1251

Спецкурс Оптические микрорезонаторы. Лекция 10. Связь с оптическими резонаторами с модами типа шепчущей галереи.
М.Л.Городецкий 26 сентября 2007 г.

(

Перевести все в СИ
i

)

С точки зрения простейшей геометрической оптики, если внутри сферического резонатора луч циркулирует так, что падает на его поверхность под углом большим чем угол полного внутреннего отражения, то такой луч никогда не сможет покинуть резонатор. Из соображений взаимности, если луч, циркулирующий внутри, не может выйти из резонатора, то и никакой луч, вошедший в резонатор, не может преворатиться в такую бесконечно долго циркулирующую волну. На одном из заданий в Вечерней Физической Школе МГУ, когда я еще был школьником, была задана следующая простая задачка (Рис. 1). Требовалось определить, как надо направить луч на стеклянный шар, чтобы после преломления луч упал на поверхность под углом большим полного внутреннего отражения. В рамках школьной программы задача решения не имеет. Легко показать, что все преломленные и отраженные лучи будут лежать в одной плоскости, проходящей через падающий луч и центр шара, и если если перерисовать картинку в этой плоскости, отрицательный ответ становится очевиден. Из соображений симметрии относительно оси

OC

угол падения

i

равен углу преломления

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

P

, и поскольку

@n

IAa)

от его поверхности, сравнимом с радиусом резона-

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

1


Рис. 1: Школьная задача

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

` a IHH (a ` a IHHH (a

9 IS мкм при a H:TQ мкм) получаем 9 P ? IH IP , 9 IRH мкм) получаем 9 T ? IH ISV ! Для того чтобы

тем не менее связаться с резонатором надо каким-то образом уменьшить потенциальный барьер, создать в нем что-то типа окошка, через которое с резонатором можно было бы связаться. Исторически впервые с кварцевым шариком удалось связаться потому, что на поверхность в области моды попала светорассеивающая пылинка [1]. Естественно, на пылинку всякий раз рассчитывать не приходится и надо изыскать более воспроизводимый способ. Вспомним квазиклассическое приближение и эквивалентный потенциал, образующий барьер для фотонов. Если мы поднесем в какой-то точке вблизи поверхности резонатора в области моды диэлектри с показателем преломления большим чем показатель преломления резонатора, мы можем тем самым создать то самое требуемое "окошко", пропускание которого можно легко регулировать, меняя расстояние до поверхности резонатора и толщину локального потенциального барьера. Эффективное возбуждение мод шепчущей галереи может быть осуществлено с помощью призменного элемента ввода, известного в интегральной оптике и служащего для связи плоских пучков с модами планарных световодов. Рассмотрим возбуждение сферического резонатора с модами типа шепчущей галереи пучком, падающим на внутреннюю поверхность призмы под углом, большим угла полного внутреннего отражения (Hис.2). Резонатор при этом помещается на некотором малом расстоянии

d

$

от призмы

так, что внешнее поле мод ШГ существенно проникает в полупространство призмы. Априори очевидно, что для эффективного взаимодействия с модами ШГ возбуждающий пучок должен быть сфокусирован в области "касания". Приведенный ниже расчет имеет целью выяснить оптимальные параметры возбуждающего пучка (ориентация, угловой спектр и т.п.). С точки зрения простейшей геометрической оптики, если угол падения луча, отражающегося от внутренней поверхности резонатора равен Фундаментальные моды типа

T E``q и T M``q

i

, то

могут эффективно возбуж-

даться с помощью гауссового пучка (Рис.2), сфокусированного на внутреннюю поверхность призменного элемента связи (с углом преломления большим чем показатель преломления резонатора

nc nm ) в области вблизи резонатора под оптимальным углом Е a rsin@nm =nc A. Для каждой моды можно рассчитать добротность связи (нагружение) Qc , которая характеризует потери энергии моды вследствие связи с излучательными модами полупространства призмы. Далее будем называть такие моды пространственными (угловыми) спектрами переизлучения мод шепчущей галереи с призменным элементом связи. Оптимальная связь достигается если соб-

2


Рис. 2: Схема возбуждения мод типа шепчущей галереи с помощью призмы

ственная добротность резонатора больше добротности связи и если угловой спектр возбуждающего пучка совпадает с угловым спектром переизлучения. Основной результат работы [2] состоял в том, что пространственный спектр переизлучения моды

Y``

представляет собой гауссов луч, совпада-

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

Y``

, наклоненная по отно-

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

второго выходного луча
Y`m

Y`m

?a

(Рис.3). Другое очевидное следствие

Ѓ ros@m=`A, возбудят появление

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

, симметричного входному лучу относительно экваториальной плоско-

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

h@IA @kH rA, `

но при больших порядках

`

может быть ап-

3


проксимировано экспонентой p резонатора и m

r? a =@P nP

ром является гауссовым

вания ближнего поля моды в дальней зоне. Для больших значений рованы следующим образом:

exp@ d=r? A, где d IA. Фактически, yP CzP для окном exp PRr?

расстояние от поверхности область связи с резонатоспектрального преобразо-

` m

,

присоединенные функции Лежандра в этом окне могут быть аппроксими-

P`m @os A 9 B`m os@@ =PA C @` mA=PA; P a `@` C IA mP ; ( @`Cm IA33@` m IA33 ; при четных ` m; @` mA3 B`m a @`CmA33@` mA33 ; при нечетных ` m: : @` mA3

(1)

Как известно из теории дифракции, поле в дальней зоне (угловой спектр) можно найти из распределения поля на поверхности экрана, в данном случае на поверхности призмы с помощью двумерного преобразования Фурье ([3]):

Z

I

I

p xP PP eikx PP dx a Pe k =P
ZZ

(2)

q I PPP F @kx a nP kH ky kz ; ky ; kz A a U @x a H; y; z Aeiky yCikz z dy dz c @PAP P P PP ? Uma` G e i`y=a z `=@Pa A @y Cz A=@Par A ; PP ? UmTa` G e imy=a @y Cz A=@Par A os@z C @` mA=PA (3) Для определенности взято распределение поля для моды TE и лишь основная компонента поля. Таким образом, выражая компоненты волнового вектора через углы (Рис. 3), получаем

U @?; ЕA G U @y; z Ae inc k sin Е os ?z inc k sin Е sin ?x dy dz a P? Par? e @m nc ka sin Е os ?A r =@PaA ? P? P? e @ nc ka sin Е sin ?A r =@PaA C @ IA` m e @Cnc ka sin Е sin ?A r =@PaA
Расчет интеграла (при

Z

(4) ,

` > I и kx ; ky < k

что справедливо при

nc > nm

при интегрировании можно перейти к бесконечным пределам) дает два гауссовых луча с одинаковыми фазами для четных нечетных

` m.

Поляризация лежит при

излучения для мод

TM

` m

и перпендикулярно ей

ный учет второй компоненты поля

E

` m и в противофазе для ( ` примерно в плоскости для мод T E . Дополнитель-

даст дополнительный, малый при

4


Рис. 3: Система координат в микросфере и элементе связи

Рис. 4: Связь с модами типа шепчущей галереи в диэлектрическом сфероиде. Входной луч связывается с копланарной фундаментальной циркулярной модой

Y``

. Прецессируя вокруг оси симметрии сфероида (совпадающей с

Y`m

осью ножки), фундаментальная мода

Y``

1) формирует стационарную моду

и 2) образует два гауссовых луча переизлучения в призму.

5


` m

m

поворот поляризации. Углы и угловые апертуры определяются

следующими выражениями (все внутри призмы):

@? ?H AP @Е ЕH AP U @?; ЕA G e P??P P?ЕP p sin ?H a Ѓ p 9 Ѓ I C Pp ` I=P ; `@` C IA " # p `@` C IA nm q ` P=Q sin ЕH a 9 n I P P n c kH a c p p P I P nm Pa ?? 9 nm I `nm nP kH a sinP ЕH cp p P I n nP Inm m ; ?ЕP a P m P 9 nc kH a os ЕH `@nP nP A c m если ` a m; то # " nm q ` P=Q ?H a H; sin ЕH 9 I nc PP p p m m P a nm C nP I 9 nm C nP I : ?? `nm nP ka sinP ЕH c

(5)

Использованные приближения для полиномов Лежандра справедливы если

> `,

p

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

ния по углу



. Однако, прямое численное интегрирование без использова-

ния приближения 1 показывает, что количественные оценки 5 существенно не меняются и двухлучевые спектры излучения характерны для всех мод вплоть до зависит от

` m a I q и ` (и в

и вплоть до малых

меньшей степени от

` $ SH. m при

Поскольку волновое число

k q

малой эллиптичности), углы

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

`

)I
P

зависимости от

pa`

m

и особенно

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

i

моды и мощности

R P P M jEj dv ! i kH n m R Qc a a P nc jEjP ds C

излучаемой в призму:

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

cn I Pa cH P @PAP

cn F P @kx ; ky Adkx dky a c H P
6

ZZ

jE jP dxdy a

(7)


Используя для аппроксимации формулу Стирлинга, окончательно получаем [4] (для простоты в окончательном выражении положено

n

nc

9 nm a

):

Qc
0.1

9

nP I Pa Q=P PkdpnP I e P n

(

p ; p IC nP I p P ` m;

q



при при

` a m; `>m

(8)

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

В работе [4] пространственные спектры переизлучения кварцевых микросфер в призму исследовались экспериментально. Наблюдения проводились в видимом диапазоне, на длине волны гелий-неонового лазера 633 нм. Поскольку частотный диапазон возможной перестройки частоты He-Ne лазера достаточно мал (менее 1ГГц; резонаторы диаметром около

?!=! < P ? IH T ), 250 мкм (` $ IVHH),

использовались большие что обеспечивало доста-

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

I Ѓ H:QA7).

Схема эксперимента очевидна из (Рис.2,4). Лазерный луч фокусировался с помощью широкоапертурного микрообъектива (фокусное расстояние 5 мм, числовая апертура 0.8) в призме (стекло К8,

a TQQ
? Е

nc a I:SIS

для

нм). Фокальное пятно помещалось в точке полного внутреннего от-

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

d

$ H:I)

поверхности призмы. Ориентацию входного луча (углы

) можно было изменять смещением лазера в плоскости перпендику-

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

Qc

(при

этом ширина резонансной кривой превышала диапазон перестройки и не

7


Рис. 5: Оптическая установка для наблюдения "чистых"спектров переизлучения мод типа шепчущей галереи и непосредственного наблюдения распределения поля в зоне контакта резонатора с призмой.

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

QH

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

q a I и @Q ? TA ? IHU ` m

для более глубоких мод с

QH

изменялась между

q

$ S.

@T ? VA ? IHV

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

$ IHH

Qc

$ P ? IHS

для мод с

` m a H и Qc

$ P ? IHT

в согласии с формулой 6. После этого свипирование

частоты лазера отключалось. Меняя углы

?

и

Е

можно было наблюдать

эффективную связь и переизлучение различных групп мод в соответствии с результатами предыдущего раздела. В нагруженном состоянии резонансы соседних мод в частотном области заметно перекрывались. При этом на заданной частоте лазера и при накрытии конусом отраженного луча лазером нескольких угловых пучков переизлучения от разных мод можно было одновременно возбуждать и несколько мод одновременно. На экране при этом наблюдались различные интерференционные картины, обусловленные перекрытием спектров переизлучения от разных мод и от входного пучка. После этого связь уменьшалась, резонатор отодвигался от поверхности призмы, и частота призмы подстраивалась под одну из мод с сузившимся резонансом. В этом случае мы наблюдали одномодовые угловые спектры переизлучения добавленные к отраженному лучу. (Рис.6.a). При включении медленного свипа частоты, наблюдались попеременные вспышки на экране пятен переизлучения различных мод. Для наблюдения чистых спектров переизлучения от мод мы использовали дополнительную призму, подносимую с другой стороны резонатора (схема эксперимента показана на (Рис.5). Выходной пучок попадал в апертуру входного объектива и выходил назад в сторону лазера, формируя два параллельных луча, соответствующих двум угловым максимумам

8


Рис. 6: Угловые спектры переизлучения с одной и двумя призмами.

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

` m a R.

Показанная на Рис.5 схема позволила получить и действительное изображение конфигурации поля моды области связи между призмой и резонатором. Это достигалось небольшим отводом входного объектива от призмы и формированием при этом оптической схемы типа микроскопа. Действительное изображение поля моды при этом получалось с увеличением около 1000 (размер изображения на экране около 1 см). Соответствующие изображения для разных мод показаны на Рис.7. Они воспроизводят фрагменты соответствующего распределения ченный для моды с при больших получающееся в области контакта призмы с резонатором. Рис.7.c, полу-

P`m @os A

видимый через гауссово окно,

`

m 9 UH

иллюстрирует корректность косинусного

приближения (1) для широтного распределения поля вблизи области связи

` m.

Кроме того, были непосредственно сфотографированы зоны локализации светового поля вблизи поверхности резонатора для разных мод. В частично нагруженном резонаторе (

Q

$ I ? IHU

) при остановленном сканиро-

вании частоты можно было удерживать моду в резонансе в течение достаточно долгого времени. При мощности лазера

$I

мВт и при времени экспо-

зиции около 1 мин оказалось возможным через микроскоп сфотографиро-

9


Рис. 7: Картины распределения поля мод типа шепчущей галереи в ближней зоне. (a)

` maH

(b)

` maP

(c)

` m

9 UH

вать картины спеклов, получающиеся в результате слабого поверхностного рассеяния на остаточных неоднородностях молекулярного размера энергии моды в области ее локализации (Рис.8). Идентификация мод в наших экспериментах осуществлялась посредством сравнения наблюдаемых спектров переизлучения с расчетными. Представленные в этой главе результаты иллюстрируют плодотворность прецессионного подхода для предсказания свойств высокодобротных мод шепчущей галереи в реальных твердотельных микрорезонаторах с малой асферичностью. Хотя понятие прецессии можно применить к описанию стационарных полей лишь условно, оно дает ясную физическое интерпретацию природы наблюдаемых явлений: 1) снятие частотного и пространственного вырождения мод шепчущей галереи в реальных сфероидальных резонатора; 2) возбуждение

T E @T M A`mq

мод наклонными гауссовыми пучками; 3)

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

Qc a QH

для выбранной моды.

10


Рис. 8: Микрофотографии спекл-картин, получающихся вследствие рассеяние на остаточных поверхностных неоднородностях для различных мод типа шепчущей галереи, полученные для того же резонатора 4. a) b)

` m a P,

c)

` m

9 IS,

d)

` m

9 UH

` m a H,

11


Рис. 9: Схема возбуждения микросферы с помощью призмы

IH T
1

? IH IH

с.

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

Моды типа шепчущей галереи можно возбудить в сферическом резонаторе несколькими способами. Для связи со свободными оптическими пучками наиболее естественным представляется использование, широко применяемого в интегральной оптике для возбуждения мод в планарных волноводах, призменного элемента связи [2, 4], основанного на эффекте нарушенного полного внутреннего отражения. Этот элемент подробно обсуждался в предыдущей главе. Для возбуждения резонатора луч лазера с помощью объектива вводится в призму и фокусируется на внутренней поверхности на которой происходит полное внутреннее отражение, а резонатор подносится к области фокального пятна снаружи на расстояние около

H:I

мкм.

Весьма интересен недавно предложенный и опробованный метод [5] при котором резонатор изготавливается из специального германиевого стекла, способного изменять в небольших пределах и запоминать изменения показателя преломления под действием ультрафиолетового облучения (подробнее см. [6]). На часть боковой поверхности готового резонатора ультрафиолетовым лучом с помощью маски наносится периодическая структура c пространственным периодом

? a =@n

sin A (

- требуемый угол падения

возбуждающего сфокусированного луча), также обеспечивающая согласование c свободным лазерным пучком. Такой элемент связи аналогичен также применяемому в интегральной оптике возбуждению планарных волноводов через гофрированный участок. Для практических применений весьма актуальна разработка элементов связи сферических микрорезонаторов с волоконными световодами. К настоящему времени продемонстрированы три способа: 1) Связь через сполированный сбоку до сердцевины участок волокна (аналогично тому как это делается в волоконных ответвителях) [7, 8, 9].

12


Рис. 10: Схема возбуждения микросферы с помощью сполированного волокна

Рис. 11: Схема возбуждения микросферы с помощью растянутого волокна

2) Связь через вытянутый на конус участок волокна [10], при этом моды сердцевины переходят в моды оболочки, проникающие в окружающее пространство, и в некоторой области конуса с оптимальной толщиной обеспечивается синхронизм с модами резонатора [10]; 3) Связь через торец косо срезанного под оптимальным углом волокна ("поросячий хвостик") [11]. Такой срез в общем-то аналогичен призме полного внутреннего отражения, и два волокна со срезами обеспечивают эффективный ввод и вывод. Некоторой комбинацией призменного и последнего волоконного элемента является связь через срезанный под углом в области фокуса градан (градиентная линза, для свободных пучков, представляющая собой стержень с градиентным по сечению распределением показателя преломления, фокусное расстояние такой линзы зависит просто от длины стержня). Такой элемент связи, как бы совмещающий в себе и призму и объектив недавно был успешно применен в схеме стабилизации полупроводникового лазера с помощью микрорезонатора [12]. Все эти системы ввода излучения в резонатор основаны на эффективном обмене энергии между направляемой циркулирующей в резонаторе направ-

Рис. 12: Схема возбуждения микросферы с помощью косо срезанного волокна

13


Рис. 13: Схема возбуждения микросферы с помощью градиентной линзы

ляемой волной полного внутреннего отражения (модой ШГ) и выпадающим полем моды волновода или пятна полного внутреннего отражения элементе связи. Априори очевидно, что эффективную связь можно ожидать при выполнении двух важных условий: 1) фазовый синхронизм и 2) существенное перекрытие полей двух волн, моделирующих, соответственно, моду ШГ и моду элемента связи. Хотя разумная связь и была продемонстрирована с разными элементами связи (до нескольких десятков процентов отбора мощности входного пучка в резонансе), их разработка и оптимизация сдерживалась отсутствием систематического теоретического подхода, позволившего бы количественно характеризовать работу разных устройств. Оставалось также неясным - возможно ли вообще обеспечить полный обмен энергией между бегущей волной в элементе связи и данной модой шепчущей галереи в сферическом высокодобротном микрорезонаторе. Ответы на эти вопросы являются критически важными для многих предложенных применений микросфер и особенно для экспериментов в области квантовой электродинамики резонатора. Ниже представлен общий подход для описания связи в ближнем поле высокодобротных мод шепчущей галереи с модами, распространяющимися в призме, планарной или волноводной структуре. Представляемые теоретические результаты дают полное описание и дают рецепты для получения оптимальной связи в элементах разного типа. Следует обратить внимание на важность для достижения оптимальной связи введенного параметра

Qc

нагруженной добротности и его соотношения с собственной добротностью

QH

.

1.1

Общие соображения

Рассмотрим возбуждение отдельной высокодобротной моды типа шепчущей галереи с высокой добротностью с помощью ных мод (

N

бегущих мод в элементе свя-

зи. При этом элемент связи может иметь бесконечное число пространствен-

Na

I

), как, например, призма [2, 4] или планарный волновод,

или только одну моду (

N a I,

одномодовая связь), связанную с резонато-

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

14


Рис. 14: Схема возбуждения мод шепчущей галереи в высокодобротном сферическом микрорезонаторе

Пусть

ражения в резонаторе (см. Рис. 14), моделирующая моду шепчущей галереи. Пусть мощность накачки распределена между ортонормированными модами в элементе связи так, что )и P поля, так что (

AH @tA N

амплитуда циркулирующей моды полного внутреннего от-

I

k

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

Bk @tA представляет амплитуду моды k B in @tA - медленно меняющаяся амплитуда возбуждающего in jBk @tAjP a jB in @tAjP определяет полную мощность накачки

Pa

резонатора, так что можно ввести локальные действительные

амплитудные коэффициенты передачи

Tk

, описывающие связь резонатора

со всеми модами элемента связи (как направляемыми, так и вытекающими) и коэффициент внутреннего отражения

R

. Будем обозначать массивы

коэффициентов передачи и амплитуд как векторы

Если добротность моды резонатора достаточно велика, то каждый оборот моды дает только малый вклад в установление моды и поэтому

TиB

, соответственно.

I

I

R (

. В этом случае (пренебрегая для простоты поглощением и рассеянием в

элементе связи, то есть, считая, что выполняется закон сохранения энергии

15


P P Rk a I Tk

), получаем

Ra

Y

Rk a I



X

P Tk =P a I T P =P

(9)

Уравнение для амплитуды моды резонатора и амплитуд в элементе связи запишется в следующем виде: X in AH @tA a i Tk Bk @tA C RAH @t H AeiPnm L= L=P out in Bk @tA a Rk Bk @tA C iTk AH @tA: (10) где H a nm L=c время одного оборота бегущей волны в сфере, а L 9 Pa примерно равна длине окружности сферы, длина волны, nm показатель преломления резонатора, c скорость света и коэффициент затухания мощности в резонаторе, обусловленный потерями на рассеяние, поглощение и излучение в резонаторе. В изложенном представлении микросфера эквивалентна кольцевому резонатору, образованному зеркалами с пропусканием глощающей средой с коэффициентом ослабления Перо длиной



Tk

, заполненному по-

, а в случае одномодовой

связи, как это было справедливо отмечено в работе [8], резонатору Фабри-

L=P

с глухим задним зеркалом. , где

Pc=H , nm L a mH AH @tA H dAH =dt из

Если потери на один оборот малы, то вблизи резонансной частоты:

m

целое число, раскладывая в ряд

(10) получим:

!H a AH @t H A a
(11)

dAH T C @c C H C i?!AAH a i B in ; dt H
где

Здесь введен

c I R TP H a Y c a a : Pnm RH PH коэффициент I:
a

(12)

TB
TB

in in

;

(13)

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

Pnm =, а c

H

обусловлено собственной добротностью резонатора

описывает нагружение, т.е. утечку энергии моды в моды эле-

QH a

мента связи. Понятно, что связь с резонатором является взаимной энергия не только поступает в моду резонатора из мод накачки, но и утекает в них. Именно балансом собственных потерь, коэффициента связи и нагружения описывается установившийся режим в резонаторе. Ниже все величины, связанные с элементом связи (coupler) будут обозначаться индексом ` ', а величины связанные со сферическим микрорезонатором `

m

c

'.

Получившееся уравнение (11) является классическим уравнением для медленно меняющейся амплитуды осциллятора под действием гармонической внешней силы.

16


Как будет показано ниже, коэффициенты

Tk

могут быть вычислены как

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

Tk

не являются фиксированными

параметрами, а сильно зависят от геометрии связи, и в частности, экспоненциально зависят от величины зазора между резонатором и элементом связи, и, таким образом, находятся в руках экспериментатора. Как уже подчеркивалось в [2], именно управляемое соотношение между и

H

c

определяет эффективность связи при заданной конфигурации, при

этом это соотношение учитывает одновременно перекрытие мод, синхронизм и оптимальное нагружение, обеспечивающее обмен энергией между резонатором и элементом связи. Стационарное решение (11) имеет типичный для теории колебаний вид: r Pc i B in iPc B in a AH a H C c C i?! T H C c C i?! H Амплитуда поля в резонаторе будет максимальна при

(14)

ственная добротность равна добротности связи, или полная "нагруженная"добротность в два раза меньше собственной). Выходная стационарная амплитуда при этом равна

c a H

(соб-

B
висимости:

out

a Bin

B in C PcC i?! T ; T H c

(15)

и полная выходная интенсивность имеет Лоренцов профиль частотной за P jBout jP a jB in jP I @ CRHAcC @?!AP ; H cP

(16)

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

a I),

получающийся при

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

B

in

=B in a T=T

, когда доля

c a H

, выходная интенсивность обращается в ноль, т.е. вся

K

на резо-

нансной частоте в частотной передаточной характеристике. Из выражения

Ka

RQ P RQH Qc P a @QH C Qc AP QH C Qc
17


Рис. 15: Моделирование связи с микрорезонатором направляемым ответвителем

I PH Pc I I a C a C Q ! ! QH Qc
При критической связи

(17)

K a IHH7

. Однако критическая связь может на-

блюдаться и в случае неидеального соответствия до тех пор пока

P P > I

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

Q,

когда

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

c a H =@P P IA.

Если

c

) H
.

(сильное нагружение), то



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

1.2

Моделирование связи волноводным ответвителем

Задача теперь состоит в том, чтобы определить параметры системы элемент связи резонатор из электродинамики. В недавней статье Д.К.Роуланда и Дж.Д.Лава [13] (последний известен основательной книгой по теории оптических волноводов, написанной им в соавторстве с А.В.Снайдером [14]), проблема связи с модами типа шепчущей галереи рассматривается на основе модели распределенной связи между бегущей поверхностной цилиндрического резонатора и заданной моды планарного волновода. В этом подходе вместо локальных коэффициентов передачи рассматривается распределенная связь, аналогично тому как анализируются распределенные волноводные ответвители. В наших обозначениях этот подход

18


ведет к следующей системе дифференциальных уравнений:

Коэффициенты сел циенты связи

H

и

k

мод резонатора и элемента связи) и распределенные коэффи-

dAH a i?H @z AAH C iCk @z Aei@k H Az Bk dz dBk a i?k @z ABk C iCk @z Ae i@k H Az AH ; dz ?H @z A и ?k @z A (описывающие возмущения

(18) волновых чи-

Ck

могут быть вычислены в явном виде через интегралы

перекрытия полей (см. [13] и ссылки в этой работе): Z !@nP IA c jeH jP dsY V C Z !@nP IA m jek jP dsY ?k @z A a V M Z Z !P Pa P IAe? eH ds @nP IAe? ek ds Ck @nc k m H TRP C M

?H a

(19)

Здесь

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

eH и ek

представляют эквивалентные направляемые волноводные мо-

MиC

означают что интегрирование проводится внутри микрорезо-

P Ck

были равны, и это часто просто по-

Tk

могут быть найдены

из уравнений (18). Однако аналитические выражения выходных полей кроме нескольких простейших случаев не могут быть найдены из уравнений (18). Возможно именно поэтому авторы [13] предложили лишь численное решение для своего частного случая. Более того, в общем случае передаточная матрица является комплексной матрицей 2x2 и не может быть описана единственным действительным параметром. Ситуация к счастью более благоприятна для микросфер с высокой добротностью связи следует, что диаметром

Qc a !=Pc , когда Tk ( I. И действительно, T P 9 Pm=Q. Если взять резонаторы из плавленого IRHm (m 9 IHHH) с сильным нагружением Q $ T ? IHU
19

из (12) кварца (в дей-

ствительности добротность в таких резонаторах может быть порядка

IHIH

)


то

T

9 I7

. Практически

да поля в резонаторе

зи меняется незначительно и поэтому может быть положена константой во втором уравнении (18), а также то, что стационарная амплитуда Поэтому оказывается возможным получение приближенного решения:

AH

T

имеет порядок

IH Q

. Это означает, что амплиту-

при однократном прохождении через область свя-

AH

) Bk

.

Aout a RAin e H H out in Bk a Bk C
где

iЕ

in C iTk Bk iTk AH ;
Z

(20)

Tk a

Z

I

I

Ck e

i@H k Az

dz Y

Еa

I

I

?H dz

(21)

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

!H



Е !@nP IA c H a !H a H VH

Z

C

jeH jP dv
 Z    

c
1.3

T P !P @nP IAP c a ka PH IPVP H

C

e? eH ei@H k A dv k

     

P
(22)

Вариационный подход

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

E m @r
где

; tA a e i!t

X

j

e Aj @tA^j @rA;

(23)

не имеющего собственных потерь.

^j e

ортонормированные собственные моды невозмущенного резонатора, Z

I R
(

"m^j I^?P dv a j I;j P e ej

(24)

j I;j P

здесь символ Кронекера).

20

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

Qrad

[15], однако в нашем рассмотрении мы можем легко из-

бежать эти проблемы, считая интересующие нас собственные частоты

!j

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

Aj

Qrad =

, но много

являются медленно

меняющимися во времени и отличаются от ранее введенных медленно ме-

Aj

циркулирующих волн не только

заменой переменной, но и нормировкой Легко видеть, что Z

j^j e jej

jP a c ^ h? ds a I : ^ jP R ej j j
M

(25)

Уравнение для поля в сфере в присутствии элемента связи запишется в следующем виде: "m @rA "c @rA I PH "m @rA @ P Em " @ rA I @ P E c Em C cP C cP C i ! cP a m P ; @ tP c @ tP H

r?r?

(26)

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

nPjc r

внутри и 1, соответственно, внутри и сферического резонатора и эле-

"rjc @rA

равны

мента связи. Подставляя (23) в (26) и затем умножая это уравнение на лучаем.

после интегрирования по всему объему, пренебрегая малыми членами, по-

^? eH

d AH i!@nP IA i!t m C AH @H C i?!H A a e dt V

Z

где

H ?!H a !H

!

M

Ec^? eH

dv;

(27)

и

H !H a !H



! V

Z

C

@nP c

IAj^H jP dv: e

(28)

является новой резонансной частотой сдвинувшейся в присутствии элемента связи, в полном согласии с (22). Поле в элементе связи запишем в виде разложения по модам, бегущим в направлении

z

(см. Рис. (14)):

Ec @r

; tA a e i!t

Z

B @z ; tAe @rAeiz d ; B
как

(29)

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

B a

k

~ Bk @ k A C B
21

(30)


Моды элемента связи нормируем так, чтобы

c R
(здесь

Z

eI ; h? P z ds a @I 

P A
zиt

(31)

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

@I P A z

дельта функция и

h

амплитуда магнитного поля, соот) описывают

. Амплитуды

B

(медленно меняющиеся по

распределение волны накачки в элементе связи. Подставляя (29) в волновое " @ rA I @ P E m "c @rA "m @rA I @ P Ec a c P ; Ec C cP C cP @ tP c @ tP

r?r?
получим

(32)

Z

@B  @z



i!P @"m IA B e ei PcP i@" IA!P a cP Em Pc

z

d
(33)

ei!t

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

h

?



, и проинтегрировав по поперечному сечения,

I P @nP IA Z !m e ; h? z ds; ? @z A a  Vc M

B a

in B ei?z

i!P C Vc

Zz

e

i!t i zH Ci? @z zH

A

Z

C

@nP c

IAEHm ; hH? z dsH dzH 
(34)

Подставляя (29) в (27) и используя (34) пренебрегая

?

, окончательно

получаем следующее выражение для амплитуды поля в резонаторе:

d AH i!@nP IA c C @H C c C i?!H AAH a dt V
и

Z

in B

Z

C

^H e? e

e iz dvd ;

(35)

M I C  P  Z Z !P @nP IAP  c ? e iz dv d ;   9 ^H e e   IPVP   C
22

c a !Q TRP c

Z Z Zz Z

@nP m

IA@nP IA ei@C?A@z zH A ^H ; hH? @e ^? AdsH dzH dvd c eH  z  eH 
(36)


в естественном согласии с (22). Полное согласие с (11-15) станет очевидным, если положить

C Для высокодобротных мод H 9 m=a, и так как поле спадает вне резона r ( P 9 kP @nP IA), зависимость eH от z может тора приблизительно как e m zP =Pa . Если быть аппроксимировано следующим образом: e 9 e@z a HAe элемент связи является протяженным в направлении z (как большинство
продемонстрированных до сих пор элементов связи), получаем:

T a

!@nP IA c V

Z

eH e? ei@H A 

z

dv

(37)

T a

r !@nP IA Pa c V

Z

C

P e @a mA =P a eH e? ds 

(38)

Список литературы
[1] . . . В. Б. Брагинский, Доклады академии наук

293

, 1358 (1987).

[2] С.П. Вятчанин, М.Л. Городецкий и В.С. Ильченко, Журнал прикладной спектроскопии.

56

, 274 (1992).

[3] М.Б.Виноградова, О.В.Руденко, А.П.Сухоруков, ука, 1979. [4] M. L. Goro detsky and V. S. Ilchenko, Opt. Comm.

Теория волн,
113

М.: На-

, 133 (1994).

[5] M. G. L. M. J. F. V.S. Ilchenko, D.S. Staro dub ov, Conference on Lasers and ElectroOptics, Baltimore, May 23-28, Technical Digest , 67 (1999). [6] J. F. B. K. S. J. D.S. Staro dub ov, V. Grubsky, Opt. Lett [7] S. A. A. Serp enguzel and G. Griel, Opt. Lett.

22

, 1086 (1997).

20

, 654 (1994).

[8] D. T. A. S. J. C. G. Griel, S. Arnold and N. Morris, Opt. Lett. (1995).

21

, 695

[9] D. L. V. S. J. H. V. L.-S. J. M. R. N. Dubreuil, J. C. Knight and S. Haro che, Opt. Lett.

20

, 1515 (1995).

[10] F. J. J. C. Knight, G. Cheung and T. A. Birks, Opt. Lett. [11] X. S. Y. V. S. Ilchenko and L. Maleki, Opt. Lett. [12] S. M. I. V. V. Vassiliev and V. L. Velichansky, (2003).

22

, 1129 (1997).

24

, 723 (1999).

Appl. Phys. B

76

, 521

[13] D. R. Rowland and J. D. Love, IEE Pro ceedings, J. Opto electronics 177 (1993).

140

,

23


[14] А. Снайдер, Дж. Лав, связь, 1987.

Теория оптических волноводов

,

М.: Радио и

[15] H.M. Lai, P.T. Leung, K. Young, P.W. Barb er, and S.C. Hill, Phys. Rev. A

41

, 5187 (1990).

24