удк 563.75:(535.3;535.182)
МОДЕЛИРОВАНИЕ СТРУКТУР И ХАОСА В ГЕНЕРАТОРЕ С ПРОТЯЖЁННЫМ КРИСТАЛЛОМ, ОБЛАДАЮЩИМ ТЕРМОЗАВИСИМОЙ ОПТИЧЕСКОЙ АКТИВНОСТЬЮ
Т.В. Голякевич, И.В. Измайлов*, Б.Н. Пойзнер*, В.М. Трухан, Д.А. Шергин*
Объединённый Институт физики твёрдого тела и полупроводников, НАНБ, г. Минск *Томский государственный университет E-mail: [email protected]
Построена математическая модель нелинейных процессов - в приближении протяжённой среды с термозависимой оптической активностью - для кольцевого резонатора. Построены карты динамических режимов. Сделаны оценки параметров резонатора и фосфида кадмия, благоприятных для осуществления устройства конфиденциальной связи оптического диапазона.
В кольцевых резонаторах, содержащих нелинейный элемент (насыщающийся поглотитель [1] либо материал, в котором имеет место высокочастотный эффект Керра [2-4]), открыты явления самоорганизации и детерминированного хаоса. Универсальность этих феноменов подтверждается тем, что нелинейные свойства элемента способны обеспечить многообразные физические эффекты [4]. Случай электрооптических явлений в немагнитном диэлектрике активно изучается с 1990-х гг. (например, [2-5]). Случай же оптической активности исследован явно недостаточно.
Поэтому целесообразно выяснить, насколько пригодны кристаллы, поворачивающие плоскость поляризации света, для реализации нелинейного элемента в кольцевом резонаторе. Это позволит определить границы аналогии между электрическими и магнитными явлениями применительно к процессам в нелинейном кольцевом интерферометре. Пример использования эффекта Фарадея для генерации хаотических колебаний можно обнаружить в [4].
В связи с этим возникают следующие задачи: оптико-физическое описание динамики оптического поля в нелинейном кольцевом интерферометре; построение соответствующей математической модели, которая содержит ту или иную нелинейную (например, оптически активную) среду, а также проведение моделирования.
Рис. 1. Схематическое изображение нелинейного кольцевого интерферометра
Предварительно целесообразно рассмотреть классический случай, когда кольцевой интерферометр (рис. 1) содержит керровскую среду. Ход лучей
показан для случая поворота пучка (элементом 6) в плоскости хОу на 120°. Евх и Евых - напряжённость входного и выходного оптического поля; НС - нелинейная среда длиной I вдоль направления OZ; 6 - линейный элемент, производящий крупномасштабное преобразование поля; зеркала М1, М2 обладают коэффициентом отражения по интенсивности Я, зеркала М3, М4 - «глухие».
Модель, учитывающая диффузию поляризованных
молекул керровской среды
Опишем процессы в нелинейном кольцевом интерферометре (НКИ) для случая, когда поле на его входе состоит из двух компонентов с круговой поляризацией (рис. 2):
Е(г,/) = е(0(г,/)) Л(г,0 ео8(ю/+ф(г,/))+ + е(0(г,/)+п/2) Дг,/) 8т(ю/+ф(г,0), (1) где а - основная частота светового поля, 0(г,/)=у(г,/)+О/ - угол между вектором в(0), задающим направление поляризации и осью Ох, лежащим в плоскости хОу поперечного сечения пучка (О может быть сравнима с а), О - частота синхронного вращения векторов поляризации е, лежащих в плоскости, называемой нами плоскостью поляризации; А(г,/), В(г,/), ф(г,/), у(г,/) - незначительно меняющиеся за время Т=2п/а амплитуды, фаза, положение плоскости поляризации светового поля.
На рис. 2, а, толстые линии соответствуют мгновенному состоянию векторов напряжённости при а/+ф(г, /)=30°, 0(г, /)=45°, А/В=2/3, пунктирные линии изображают возможные состояния этих векторов при произвольных а/+ф(г, /), а направление их вращения показано для 0>0. На рис. 2, б, в, толстые линии соответствуют мгновенному состоянию векторов напряжённости при а/+ф(г, 0=30°, 0(г, 0=90°, А/В=1/5, направления вращения (противоположные (б), когда а>О и одинаковые (в), когда а<О) показаны для а>0, 0>0. Пунктирные линии соответствуют компоненту поля с частотой а+О и амплитудой а(г, /), а сплошные толстые - компоненту с а-О и Ь(г, /).
Если ввести обозначения а(г,/)ЦЛ(г,/)+В(г,/))/2, Ь(г,/)ЦЛ(г,/)-В(г,/))/2, то (1) выразится через проекции Е(г,/) в форме:
В кп.1(м/+ф)
У
ю>0
'У
I 1 1
1 1 / Iх
(□>0) I = а ехр ЮСсм-О^+ф+у]}"
а б в
Рис. 2. Структура бихроматического оптического излучения Е(г, а - для записи в форме (1), б, в - для записи в форме (2)
Ex(r,t)=a(r,t)cos[(a+D.)t+q>(r,t)+^y{r,t)] +
+Ь(г,0ео8[(ю-0)/+ф(г,0-у(г,0], Еу(г,0=а(г,/)8т[(ю+0)/+ф(г,0+у(г,0]-
-Ь(г,/)8т[(ю-0)/+ф(г,0-у(г,0]. (2) Тогда моделью процессов в НКИ (с учетом многих проходов поля через контур обратной связи) является уравнение:
т„(г)Щг,0М =А(г)ди(г,о- и(г,о+Дг,о, Лг,0=п2(г)/кап<Е2нс(г,0 > г=апп2(г)/к[а2нс(г,0+Ьнс(г,0], Енс(г, /)=(1-Л)05Е(г,/)+[т(г ',/)/ст/2]ЕнС(г ',-). (3) Здесь тп и Д - время релаксации и эквивалентный коэффициент диффузии поляризованных молекул НС; Ц(г,/) - нелинейный фазовый набег в НС; п2 и / - коэффициент нелинейной рефракции и длина НС; к=а/с=2п/Х; т=т(г',/)=4(г ',0+и(г',-е(г\Ц)/ы -полное время запаздывания поля в НКИ; 7(г ',/)=2Лк(г ',/)Сп(г ) - удвоенный коэффициент передачи (по амплитуде ап) поля за один проход НКИ (к и Сп характеризуют элемент О и НС); а -коэффициент сжатия поля в элементе О; Т - период оптических колебаний.
Эта модель учитывает диффузию вещества (поляризованных молекул НС), а член /(г,/)-Ц(г,/) играет роль интенсивности источников и поглотителей вещества. В такой модели для широких интервалов значений её параметров возможны и регулярные процессы, и режим детерминированного хаоса [5]. Иначе говоря, НКИ интересен как генератор хаоса - компонент систем конфиденциальной связи.
Модель процессов в кольцевом интерферометре, содержащем протяжённую среду с изменяющимися (не)оптическими параметрами
Выберем модель оптического поля. В рамках обсуждаемой задачи вектор Е(г,/) можно описывать формулами (1) или (2) при 0=0 (рис. 3).
Рис. 3. Иллюстрация описанию поля в форме (1) или (2) при 0=0. А и В - величины полуосей эллипса, угол у задаёт его ориентацию в плоскости хОу
Дадим операторное описание динамики в НКИ с учётом взаимодействия оптических и неоптических свойств среды. Для этого необходимо найти поле на выходе контура обратной связи (КОС) Евых КОС(г,0 и нелинейной среды Евых НС:
Евых кос(г,0 ^КОС(Евхкос(г,0,Ркос),
Евыхнс(Г,0=^с(Еи нс(Г,0,Рнс).
Здесь оператор ЬКОС описывает преобразование функции Евх КОС(г,г)=ЕВых нс(г,0 в функцию Евыхкос(г,^)=Евх нс(г,0-Еи нки(г,0, т. е. преобразование поля в КОС; оператор описывает преобразование функции Евх тс(г,0 в функцию Евых нс(г,0, т. е. преобразование поля в НС; параметры операторов рКОС и рНС описывают оптические свойства КОС и НС.
Тогда, очевидно, модель процессов в НКИ есть выражение вида
Евых НС(г,0 =
=^С(ЕКОС(Евых НС (г,0,Ркос(г,0)+Евх НКИ (г,/), Рнс(г,/)) либо
Евх НС(г,/)=
=Евх нки(г,0+Ькос№с(Евх нс(г,/),Рнс(г,/)),Ркос(г,/)).
Причём, в свою очередь,
Рнс(г,/)=Р(Евх нс(г,/), Рнс(г,/), Рп(г,/)), а Рп(г,/)= =Рп(Евх нс(г,/), Рнс(г,/), Рп(г,/)), где рп - неоптические свойства среды; Рп и Р - операторы, задающие рп, рНС, соответственно. В общем случае под обозначениями рНС(г,?), рп(г,/), Евх НС(г,/), Евых нс(г,0, Евх нки(г,0 в аргументах операторов следует понимать функции, определённые на интервале времени (-да, /].
Описание распространения поля через НС
и контур обратной связи
Процесс распространения поля - с комплексной амплитудой Ес в прозрачном веществе со средним («фоновым») значением показателя преломления <п> и диэлектрической проницаемостью е=е0+5е (е0=<п>2) - можно описать параболическим уравнением
21к< п>дЕс(х,у,г)/дг=Д1Ес(х,у,г) +
к2(е(х,у,1)-<п>2)Ес(х,у,1), (4)
где Д1=д2/дх2+д2/ду2. Для поля в форме (1) при 0=0 следует рассматривать две комплексные амплитуды:
Ес+=а(г, Оехр[ /ф(г, 0+/у(г, 0] и Ес+=¿(г, /)ехр[ /ф(г, ¡)-1ц/(г, Щ.
Для простоты пренебрежём дифракцией поля в КОС, возможностью изменения направления и вида поляризации. Тогда преобразование поля в КОС НКИ описывается соотношением
Еых кос(г,0=[Як(г',0/^]Евх кос(г',/-о(г',0).
Здесь /0(г ',/) - время запаздывания поля в КОС НКИ. В соответствии с выбранной моделью поля оно полностью описывается величинами ф(г,/), у(г,0, а(г,0, ¿(г,/) или Л(г,t), В(г,/). Поэтому в сделанных приближениях описание распространения поля в НС сводится к системе соотношений между этими характеристиками на входе и выходе НС. Вид соотношений зависит от типа НС.
Так, при отсутствии дихроизма и в пренебрежении дифракцией поля в НС для а(г,/), ¿(г,/) (Л(г, /), В(г,/)) имеем
Озых Нс(г,0=С„(г)аи нс(г,0, ¿вых нс(г,0=Сп(г)Ьи нс(г,0.
Здесь Сп(г ) - коэффициент передачи (по амплитуде) поля за один проход НС НКИ, и учтено, что
аи нс(г, О-авх нс(г, /„(г, -/„(г, /))), Ь„х Нс(г,/)~ЬВХ Нс(г,/-/„(г,/-/„(г,/))), где I„ - время распространения поля в среде. Тогда, следуя [2], можно ввести параметр 7(г ',/)=2Лк(г ',/)С„(г ) - удвоенный коэффициент передачи (по амплитуде) поля за один проход НКИ. В дальнейшем ограничимся сюжетом, когда НКИ содержит вещество, обладающее естественной либо искусственной оптической активностью.
Случай термооптических эффектов
в средах с оптической активностью
Как известно, в протяжённой среде длиной /, обладающей естественной оптической активностью, угол поворота в плоскости поляризации света равен в=а/=п(п--п+)//А, где а - вращательная способность, п_, п+ - показатели преломления для двух лево- и правополяризованных круговых волн. Аналогично, для эффекта Фарадея имеем в=ЦЩ/, где V - постоянная Верде, Н - напряжённость постоянного магнитного поля.
Пусть неоптическое свойство, влияющее на оптические свойства среды, есть температура, причём влияние это правомерно описать линейной зависимостью: а=а0+а1Т и V=V0+V1T.
Как известно, пространственное распределение температуры подчиняется уравнению теплопроводности
д(рс,Т)/д1 =У(АУТ)+/.
Здесь р(х,у,1,Т) - плотность вещества, с,(х,у,1,Т) - теплоёмкость среды при постоянном объёме, Т(х,у,1) - температура, Л(х,у,1,Т) - коэффициент теплопроводности, / - плотность тепловых источников. Очевидно, что плотность тепловых источников / пропорциональна интенсивности света (<Е2нс(г,/))=а2нс(г,/)+Ь2нс(г,/)) и коэффициенту его поглощения 1-Сп2(г), т. е.
Дг,0=(1-СпШа2Ж0+Ь2нс(г,0]//. Не только температура, но и световое поле, воздействуя на среду, влияют на её оптические свойства. Предположим, что эти воздействия непостоянны во времени, р0(г) - невозмущенное значение вектора параметров, 5р11ш(г,/) и 5р(г,/) - равновесный и текущий отклики на воздействие. Тогда установление равновесного р11п(г,/)=р0(г)+5р11п(г,/) состояния вектора параметров р(г,/)=р0(г)+5р(г,/) среды происходит асимптотически. Пусть асимптотический процесс установления каждого из параметров происходит по экспоненциальному закону с постоянной времени т. Тогда для переменных воздействий справедливо уравнение
дд(г,0/д*=(Ат,(г,0-Д(г,0)/т,
Если параметры р(г,/) НС изменяются из-за диффузии с коэффициентом Д(х,у,г) её молекул (например, поляризованных или ориентированных) либо носителей заряда, то необходимо использовать уравнение диффузии
дSp (г, f)/дt=V{D^8pi{г,f))+Fi.
Здесь 7<=(Лп г(г,0-л(г,0)Л=(Фшп (г,/)-8р(г,!))/т - интенсивность источников и поглотителей молекул или носителей заряда, определяющих величину параметра р. Например, в случае естественной оптической активности вращательная способность а(г,/)=а0+5а(г,/), поэтому угол поворота плоскости поляризации в=а/=в0+5в(г,/), где в0=а0/ и 5в(г,/)=/5а(г,/). Следовательно, д5в(г,/)/д/=У[^аУ5в(г,/)]+[а1^7(г,/)-5в(г,/)]/ту. (5) Для эффекта Фарадея V(г,t)=V0+8V(г,t) зависит от температуры, при этом текущий отклик на воздействие SV(г,t) удовлетворяет уравнению дSV(г,t)/дt =V[DvVSV(г,t)]+[V1T(г,t)-SV(г,t)]/тv. (6) Динамика 8V(г,t), в свою очередь, определяет поведение угла поворота в плоскости поляризации при наличии второго воздействия - магнитного поля Н
д8в(г,()/д1=
^[ВвЧЗвгМ+М+ЗЦгМт-Звт/тв. (7) Здесь по аналогии с моделью в [1] введём эффективный коэффициент диффузии Бе(х,у,1) - коэффициент, нормированный к времени релаксации НС и к квадрату некоторого характерного размера г0 пучка: Д(х,у,г)=т„(х,у,.г)Дх,у,.г)/г02. Тогда, вновь обозначая (х,у,г):=(х,у,г)/г0, получим соответствующие уравнения (5)-(7). тадSв(г,t)/дt='V[DаeУSв(г,t)]+[а1/T(г,t)-Sв(г,t)], тVдSвV(г,t)/дt=V[DVeVSV(г,t)]+[V1 Т(г,/)^(Г,Щ, твд8вг,1)/д= =ЩВаУ8в(г,Щ+{М+8ЦгМН\/-8в(г,{)}. Согласно [6. С. 106-109], п-=п0+ас/а, п+=п0-ас/а, где п0 - показатель преломления при распространении поля вдоль оптической оси вещества. Аналогично, для эффекта Фарадея п=п0+ЦН\с/т, n+=n0-V\H\c/ю [6. С. 113]. Очевидно, что для обоих случаев п±=п0-±в/(кГ),
п±=п0-±в/(к1). Тогда моделями процессов для случаев естественной и искусственной оптической активности служат следующие формулы
д(рсДО/=У(ЛУ 7)+(1-Сп(г)2)[аУг,0+^с(г,0]Л
Евх нс(г,0=(1-Юевх(г,0+Мг;0/°/2]х
х[Евх нс+(г;—+)+ЕВх нс-(г;/-Т-)], (8)
где угол в для случая естественной оптической активности подчиняется дифференциальному уравнению
тдбв^Л/д^У^Убв^Щ+аТг^-бв^^)], (9) а в случае эффекта Фарадея - уравнению
гвд5в(г, 1)/д=У[ВвеУ5в(г, /)]+ {[П+5К(г,0]|ЯМв(г,0}. (10)
Здесь /0(г ',/) - время запаздывания поля в КОС НКИ; 4(г ',0-/0(г ',0+П)(г ',/-/0(г ',/)) - эквивалентное время запаздывания в НКИ (доля времени распространения светового поля, обусловленная наличием линейной части п0 показателя преломления п(г,0 НС и наличием времени /0; т±=/е(г ',/)-±/в(г ',/-/0(г ',/)) - полное время запаздывания для право- и левополяризованного компонента поля; 4(г,/)=0(г,/)/(2ю) - отличие во времени распространения в НС кругополяризованных волн относительно времени распространения вдоль оптической оси; р(х,у,1,Т) - плотность; су(х,у,1,Т) -теплоёмкость при постоянном объёме; Л(х,у,£, Т) -коэффициент теплопроводности НС.
Модель (8)—(10) описывает нелинейные процессы в НКИ, содержащем среду с термозависимой оптической активностью - естественной либо искусственной. Общий вид модели затрудняет её численный анализ. Попытаемся её упростить.
Модель в виде дискретного отображения. Результаты моделирования
Для преобразования соотношений (8)—(10) к модели в виде дискретного отображения, сделаем ряд обоснованных приближений:
• благодаря линейной зависимости угла вращения от температуры и толщины элементарного слоя среды правомерен переход к <Т(г,/)>г;
• допустимо пренебречь диффузией молекул, т. е. результатом применения оператора набла в выражениях У [ЛГе6У У(г, 0], У [ЛаеУ6в(г, 0];
• допустимо пренебречь инерционностью изменения (магнито)оптических свойств среды по сравнению с характерными скоростями теплопередачи;
• коэффициенты теплопередачи и диффузии постоянны в пространстве;
• благодаря малому поперечному градиенту интенсивности поля и организации эффективного теплоотвода преимущественно в направлении Oz, поперечный теплоперенос невелик, поэтому оператор ЛAT»Лд2T/дz2;
• так как молекулы с возбуждёнными электронными оболочками в кристаллах слабо подвержены диффузии, то допустимо пренебречь членом Д^ДЩг,/);
• поле на входе НКИ постоянно, и Сп(г)=Сп. Тогда модель в виде дискретного отображения
для случая естественной оптической активности имеет вид
вх НС ¡-1 а ' > вх НС ¡
-Arg(Ac¡,¡, Asa),
ьвхнс ¡ч=(Лcí¡2+Лsí¡2)0•5, Фьвх нс ,=А^(Ась, Л$ь), (11) где амплитуды авх НС и Ьвх НС нормированы на авх и Ьвх,
Аса=(1-^)ео8(Фа вх ¡)+
+(7/о/2)авх НС ¡-1е08(Фа вх НС ¡-1+вм-ад,),
А5а,=(1-^)81п(Фа вх ¡) + +( У/ ^/2)авх НС ¡-151п(Фа вх НС ¡-1 + в-1 - щ),
Ась,=(1-^)ео8(Фь вх ¡)+ +(/// °/2)»вх НС ¡-1ео8(Фь вх НС ¡-1- в-1- щ),
А5ь=(1-^)81п(Фь вх ¡) +
+(у/и/2)Ьвх НС ¡-151П(Фь вх НС м-вч-щ),
в,=-/+а1Щ //4+а1Т(0)//4, /=Ка а2ж/(1-Щ+Кь ь2нс,/(1-^),
Ка=0аК, Кь=(1-0а)К, 0а-авх2/(авх2+ьвх2), КЦ1-Д)а1(1-Сп2)(авх2+ьвх2)/2/(24Л), Фа(г,^Щг,Ъ+у(г,Ъ, Фь(г,1)=Щг,0-у(г,0. Здесь Клогично назвать (суммарным) коэффициентом нелинейности, Ка и КЪ- коэффициентами нелинейности по компоненту а и ь, 0а - долей интенсивности компонента а. Обратимся к результатам моделирования.
у у
0.8- 0.8-
0.7- 0.7-
0.6- о.е-
0.5- 0.5-
0.4- 0.4-
0.3- 0.3-
0.2- 0.2-
0.1- ▼ | ЧгГО 0.1-
0 10 20 30 а о" * 2 4 6 8 10 12 К
Рис. 4. Карты ляпуновских характеристических показателей Лc(K, у) для модели НКИ с термозависимой оптической активностью (11) при щ0=0, 0=1 (а) и для модели НКИ с керровской средой (б)
Одним из наиболее эффективных методов численного исследования свойств динамических систем является построение карт ляпуновских характеристических показателей Лс на плоскости параметров, например К и у Карты Лс(К, у) приведены на рис. 4 и свидетельствуют о возможности как регулярной (серые районы), так и хаотической динамики в модели (11).
Теперь попытаемся разработать проект экспериментальной установки и оценить параметры, обеспечивающие хаотический режим.
Вариант экспериментальной установки и требования к её элементам
Поскольку данный эксперимент предполагает наблюдение динамики оптического поля на выходе НКИ, удобно использовать непрерывный лазер видимого диапазона. Очевидно, для минимизации размеров НС (и, следовательно, НКИ) требуется возможно более высокая мощность лазерного излучения. Ориентируясь на данные табл. 1, можно ожидать, что в качестве НС пригодны кристаллы СёР2 и ZnP2. Из них предпочтительнее применять в данном эксперименте фосфид кадмия, т. к. для него температурная зависимость вращательной способности сильнее (значение коэффициента а1 выше).
Опираясь на значения параметров выбранного кристалла СёР2, оценим его длину, необходимую для обеспечения режима детерминированного хаоса в НКИ, в предположении, что диаметр лазерного пучка 2 мм. Для этого построим карту динамических режимов в модели (11) - рис. 5. Из рис. 5 следует, что оптимальная длина кристалла близка к 3 см, при этом хаотический режим возникает, если мощность лазерного излучения превышает 4,5 Вт.
Таблица 1. Параметры кристаллов ZnP-l и СдР2
Параметры сар2 7пР2
1/(1=632,8 нм), мин/(Э.см) 0,55 0,279
а (7=298 К), градус/мм 676 (1=612,3 нм) 503 (1=576,9 нм)
а, градус/(см.К) 5,74 2,3
1п( Сп)/1, см-1 1,25...2 (1=694,1 нм), 0,9...1,9 (1=1106 нм) 1,3.2 (1=694,1 нм), 1,19 (1=1106 нм)
/., Вт/(м.К) 9,65 -
Возможны различные варианты реализации НКИ с нелинейной средой: на основе двух призм (типа «крыша»), на основе четырёх либо трёх зеркал и т. п. Но, исходя из критерия минимального числа оптических элементов и юстировочных уз-
лов, а также стремления упростить процесс юстировки, целесообразно использовать трёхзеркаль-ную схему (рис. 6). Считая, что удвоенная длина кристалла не должна превышать длины НКИ, а последняя - длины когерентности излучения лазера (около 30 см), получаем ограничения на геометрическую длину резонатора: не менее 6 см и не более 30 см. Поскольку предполагается использование пучка лазерного излучения диаметром порядка 2 мм, поперечный размер кристалла, очевидно, должен быть не менее этой величины.
4 -2
О —I—I—I—I-1-1-1-
0 1 2 3 4 5 6 7 I, см
Рис. 5. Карта динамических режимов в модели (11) в координатах: длина кристалла - мощность лазерного излучения на входе нелинейной средыI. Тёмным цветом обозначены неустойчивые режимы
Оценим число, типы и точность требуемых юстировочных узлов. Пусть при 1Я/10<П1 влиянием /т на динамику в НКИ можно пренебречь, где /т -интенсивность некоторой составляющей поля после зеркала М1, т - количество обходов НКИ рассматриваемой составляющей поля. Тогда можно ввести эффективное (и нецелочисленное) количество проходов поля в НКИ
тэ=0,51п(п/)/1п(0,5/).
Рис. 6. Схема трёхзеркального НКИ с элементом ввода лазерного излучения в интерферометр (а) и изображения степеней свободы/ юстировочных элементов (б): п - нормаль к отражающей поверхности зеркала, V, Ь - вертикальная и горизонтальная оси. Наклон зеркал производится поворотом относительно оси Ь, ортогональной к п и V
Кроме того, можно ввести параметр п - максимально допустимое линейное смещение центра пучка после тэ обходов НКИ.
Пусть П/=5 %, п=5 %, диаметр пучка й=2 мм, база интерферометра Х=30 см. Тогда для выбранной схемы требования точности настройки указаны в табл. 2. В качестве элемента О можно использовать призму Дове. Тогда для неё необходим поворотный механизм (точность около градуса), например, механизм с осью вращения в центре отверстия с призмой.
Таблица 2. Требования к точности юстировки для трёхзеркаль-ного НКИ при п=5 %, Пс1=5 %, d=2 мм, L=30 см
Настроечные степени свободы Точность настройки
7=0,5 (тэ=1,08) Y=0,8 (тэ=1,64)
3 наклона 27" 7"
1 поворот 10' 1'
1 сдвиг вдоль п 50 мкм
Следовательно, несмотря на достаточно жёсткие требования к оптической схеме, она принципиально реализуема, т. к. коммерчески доступны оптико-механические узлы с приемлемыми параметрами. В качестве источника излучения может использоваться Кг+-лазер (диапазон длин волн 400...800 нм, мощность до 20 Вт) [7. С. 97].
Заключение
Предложено обобщённое операторное описание сложной динамики в модели нелинейного кольцевого интерферометра с учётом трёхстороннего взаимодействия: оптического поля, оптических и иных свойств среды.
СПИСОК ЛИТЕРАТУРЫ
1. Garcia-Ojalvo J., Roy R. Spatiotemporal communication with Synchronized Optical Chaos [Электронный ресурс]: Статья. -2000. - Режим доступа: http://xxx.lanl.gov/abs/nlin.CD/0011012, свободный.
2. Новые физические принципы оптической обработки информации / Под ред. С.А. Ахманова и М.А. Воронцова. - М.: Наука, 1990. - С. 263-326.
3. Izmailov I.V., Lyachin A.V., Poizner B.N., Shergin D.A. Discrete Maps as a Model of Spatial Deterministic Chaos // Nonlinear Phenomena in Complex Systems. - 2006. - V. 9. - № 1. - P. 48-67.
4. Розанов Н.Н. Оптическая бистабильность и гистерезис в распределённых нелинейных системах. - М.: Наука, 1997. - 336 с.
В приближении отсутствия дифракции построена модель нелинейного кольцевого интерферометра, содержащего среду с диффузией поляризованных молекул и оптической активностью (естественной либо искусственной, обусловленной эффектом Фарадея). Учитываются как термооптические эффекты, так и обратное («оптотермическое») влияние светового поля на температуру среды, что делает её нелинейной.
С помощью рассчитанной и построенной карты динамических режимов в модели нелинейного кольцевого интерферометра оценены значения параметров кольцевого резонатора и нелинейного кристалла, а также уровня мощности лазерного излучения, при которых в системе реализуется режим детерминированного хаоса. Проведён анализ вариантов схемы эксперимента, обоснованы требования к числу, типу и точности юстировочных элементов.
Согласно выполненным оценкам, для Кг+-ла-зерного пучка диаметром 2 мм оптимальная длина кристалла CdP2 близка к 3 см, а хаотический режим возникает, если мощность лазерного излучения превышает 4,5 Вт.
Предлагаемая модель расширяет перспективы разработки генераторов детерминированного хаоса оптического и СВЧ диапазонов.
Исследование поддержано грантом Президента РФ МК-4701.2006.9.
Работа доложена на VIII Международной конференции «Atomic and Molecular Pulsed lasers», Tomsk, 10-14 September, 2007.
5. Измайлов И.В., Лячин А.В., Пойзнер Б.Н. Детерминированный хаос в моделях нелинейного кольцевого интерферометра. - Томск: Изд-во Том. ун-та, 2007. - 256 с.
6. Ярив А., Юх П. Оптические волны в кристаллах. - М.: Мир, 1987. - 468 с.
7. Справочник по лазерной технике: Пер. с нем. / Под ред. А.П. Напартовича. - М.: Энергоатомиздат, 1991. - 544 с.
Поступила 07.12.2007г.