Физико-математическое моделирование многоэтапных процессов формирование одномодовых интегрально-оптических волноводов в стекле K-8 методом ионного обмена

Физико-математическое моделирование многоэтапных процессов формирование одномодовых интегрально-оптических волноводов в стекле К-8 методом ионного обмена

М.М. Векшин, В.А. Никитин Кубанский государственный университет, Краснодар

Аннотация: Проведено физико-математическое моделирование многоэтапных процессов формирования волноводов в стекле К-8 на основе совместного решения двумерного уравнения диффузии, уравнения электростатики и волнового уравнения. Уравнение диффузии решалось путем использования неявной конечно-разностной двухслойной схемы Крэнка-Николсона. Рассмотрен общий случай трехэтапного технологического процесса, включающего в себя термический ионный обмен "серебро-натрий", электростимулированное заглубление волноводного канала и его отжиг. Определены условия формирования одномодовых волноводов, которые имеют размеры модового поля, подходящие для их эффективной торцевой стыковки с волоконными световодами на длине волны 1.55 мкм. Построена модель волновода с селективным (частичным) заглублением, который обладает асимметричным профилем показателя преломления и гибридными модами.

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

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

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

интегрально-оптических схем различного функционального назначения [3-4]. Цель настоящей работы построение двумерной модели процесса формирования оптических волноводов в силикатном стекле К-8 методом многоэтапного термического и электростимулированного ионного обмена.

Формирование пассивных одномодовых волноводных структур в стекле состоит из последовательности этапов: термического ионного обмена; заглубления волновода путем стимулированной миграции ионов; термического отжига волновода. Процесс формирования волновода производится путем ионного обмена Ag+•oNa+ в стекле, помещенном в расплав смеси солей AgNO3:NaNO3. Заглубление волноводного канала производится путем приложения внешнего напряжения к подложке, помещенной в расплав соли NN0^ Последний этап отжига применяется только для подстройки параметров сформированного волновода. На рис.1 показаны этапы технологического процесса изготовления волноводов [3-4].

Рис. 1. - Основные этапы технологии формирования волноводов в стекле методом ионного обмена: (а) термический ионный обмен, (Ь) - электростимулированное заглубление волновода, (с) отжиг

Описывающее любой этап данного технологического процесса уравнение диффузии имеет вид [3, 5]:

где САё - относительная концентрация внедряемого иона серебра , СА8=сА§/с0, сАё - концентрация ионов серебра, с0 - начальная концентрация ионов натрия в стекле; ДАё - собственный коэффициент диффузии ионов серебра, М -отношение собственных коэффициентов диффузии ионов серебра и натрия, М=ДА<ё1Дш; д,/, к, Т - заряд электрона, коэффициент корреляции, постоянная Больцмана и температура по Кельвину соответственно. ? - время диффузионного процесса. Е^х,^) - напряженность стимулирующего электрического поля в стеклянной подложке. Решение уравнения (1) -пространственное распределение концентрации внедренных ионов серебра определяет профиль показателя преломления формируемого волновода:

где п - показатель преломления стеклянной подложки, Лп - максимальное приращение показателя преломления волновода.

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

Нелинейное дифференциальное уравнение (1) в общем случае не имеет

и

{Х у ) = п8 +ЛпСАё {Х У ),

аналитического решения. Для решения уравнения диффузии использовался конечно-разностный метод Дюфорта-Франкеля [5]. Установлено [6], что стабильности у решения согласно этому методу нет, если М меньше, чем 0.1, и от плотности дискретизации по оси времени и по пространственным координатам стабильность зависит весьма незначительно. Значение величины М для ионного обмена Ag+•Na + зависит от состава стекла, и для обычного силикатного стекла она меньше, чем единица [7, 8], а иногда значительно меньше единицы [9].

Реконструкция параметров ионного обмена Ag+•Na+в стекле типа К-8, которая было проведено авторами, приводит к значению параметра М=0.02+0.03, что сходится с аналогичным значением М=0.01 [9] у стекла типа ВК-7, его зарубежного аналога.

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

В декартовых координатах определяющее уравнение (1) в случае двух пространственных измерений имеет вид:

А*

д 1 -(1 - м ^

д 2С д2С Л2. | +.

(1 - М) Грс ЛЯ ^

дх2 ду2 1 -(1 - М )С

Ля

уЧ V - у у

Расчетная область покрывается равномерной сеткой с шагом АИ = Ах = Ау. Производные от концентрации в уравнении

дх

Л

Ля

ду

Я

1кТ

( дС. дС. Л дх ду

:

аппроксимируются конечными разностями. Согласно методу Крэнка-Николсона, дискретизированное уравнение в матричной записи имеет вид:

г^Б+1 Г^Б

СА§ ~ СА§ = £

Л(

( с*+1 + с* ^

где С^ = ]с1 ] - вектор-столбец значений концентрации в узлах сетки общим

число пххпу на шаге *, пх и пу - число узлов сетки по обеим пространственным координатам; введена единая нумерация узлов сетки с индексами у,

к Ь

(с*+1 + с* ^

дискретизированная правая часть уравнения (2).

Уравнение (3) трансформируется в разностное уравнение с постоянными коэффициентами, если принять значение концентрации в

знаменателях коэффициентов в правой части уравнения (2) равным его значению С*. из предыдущего шага б [10]. Кроме того, в целях линеаризации

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

(дС Л ду

в узле (¡у) представляется следующим

образом:

Ag V 1+1, У 1-1, уМ 1+1, У 1-1, ]'

ду

2ЛИ

2Лк

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

(5С„„ Л

Ag

V дх ;

Применение метода Крэнка-

Николсона сводится к решению на каждом шаге по времени системы линейных уравнений с разреженной матрицей.

Решение уравнения (2) требует заданных условий нормализованной концентрации ионов Ag+ на всех границах расчетной области. Вычислительное окно выбирается таким образом, чтобы его боковые стороны и дно находились достаточно далеко от формируемого волновода и на этих границах можно использовать условие Дирихле 0^= 0. В случае

2

2

2

2

границы стеклянной подложки с расплавом граничные условия определяются геометрией маски и разновидностью технологического процесса. В месте щели в маске, где расплав соли соприкасается со стеклянной подложкой, нормализованная концентрация серебра С'Лё принимается равной 1 в случае термического ионного обмена и нулю в случае заглубления уже созданного волновода под действием внешнего электрического поля. В том месте, где металлическая маска закрывает доступ расплава к стеклу, применяется условие Неймана дСЛё/ду=0, что означает, что потока ионов в вертикальном направлении нет.

По расчетным значениям концентрации находится профиль показателя преломления волновода п(х, у ) и, далее, оптические характеристики волновода определяются путем его модового анализа. Расчет мод волновода проводился путем решения волнового уравнения методом конечных разностей [11]. Скалярное волновое уравнение для модового анализа имеет вид:

д2 Е д2 Е (, 2 ( у „2 ) -

+ + Кп(У) - Р у= 0,

дх ду

где Е=Е(х,у) - распределение электрического поля моды, в - константа распространения моды, к0=2п/к - волновой вектор вакуума, к - длина волны.

Реконструкция параметров волноводов Ап, 0Лё, М была проведена по измеренным значениям эффективных показателей преломления изготовленных многомодовых планарных Ag+-волноводов. Получены

15 2

следующие данные: коэффициент диффузии Д4Я=0.7-10~ м /с при температуре ионного обмена 380 0С, приращение показателя преломления на поверхности Ап=0.057 (для изготовления использовался расплав нитрата серебра и нитрата натрия в соотношении молярных долей 1:10), отношение коэффициентов диффузии ионов серебра и натрия М=0.03; /=0.5.

Результаты моделирования трехэтапного процесса создания Ag+-

волновода, эффективно стыкуемого с волокном на длине волны 1.55 мкм, представлены на рис.2. Двумерное распределение нормализованной концентрации внедренных ионов серебра после первого, второго и третьего этапа формирования волновода показаны на рис.2(a1-c1), а на рис^^^) показаны соответствующие распределения электрического поля фундаментальной моды ТЕ0.

Рис. 2. - Распределения концентрации ионов серебра в волноводе - (г1^1) и соответствующие распределения электрического поля фундаментальной моды - ^2^2); (а1^2) - после термического ионного обмена, (Ь1-Ь2) -после электростимулированного заглубления волновода, ^1^2) - после

отжига волновода

Численные результаты проведенных расчетов состоят в следующем. Волновод, полученный путем термического ионного обмена Ag+•^Na+ из расплава солей AgNO3:NaNO3 (1:10) через маску с шириной щели 2 мкм в течении 15 минут при температуре 380 °С, заглубленный за время 1500 с при

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

п = п + Ап exp

( (

2 2

ху

+

й2 й2

V V х у у у

со следующими параметрами: приращение показателя преломления волновода Ап=0.0057, эффективная ширина и высота 2йх=8.6 мкм, 2йу=8.0 мкм.

Нормированное на свое максимальное значение электрическое поле также было аппроксимировано двумерной гауссовой функцией:

Е = exp

( ( 2 2 ЛЛ

х у — + —

22 и и

V V х у у у

где значения параметров 2их и 2иу (эффективных размеров поля моды по горизонтальной и по вертикальной оси) составляют 10.6 мкм и 9.9 мкм, соответственно. Расчетное значение эффективности торцевой стыковки волновода, показанного на рис.2(с1-с2) со стандартным волоконным световодом БМР-28, применяемым в линиях оптической связи, составляет 99.2%.

Установлено, что различие процесса заглубления Ag+-волновода при различных напряжениях состоит в том, что форма распределения концентрации является более симметричной по вертикальной оси при низких напряжениях и больших временах процесса заглубления. Лучшее согласование волновода с волокном при таких стимулирующих напряжениях заглубления обусловлено большим влиянием диффузионной составляющей данного процесса по сравнению с дрейфовой составляющей.

Далее на основе разработанного способа решения уравнения диффузии

было проведено исследование более сложного варианта формирования волноводов с селективным (частичным) заглублением канала. Суть такой методики состоит в том, что для электростимулированного заглубления волновода на поверхность стекла наносится маскирующий слой, который располагается на поверхности подложки таким образом, что он закрывает часть канала, сформированного на первом этапе ионного обмена. В этом случае заглубляется только та часть волновода, которая не защищена маской. Если нанести маску клинообразным образом относительной продольной геометрической оси волновода, то это даст возможность создать плавный волноводный переходной канал, связывающий заглубленный и приповерхностный волновод [12], что может использоваться для создания эффективно стыкуемых с волокном волноводных микролазеров, различных интегрально-оптических датчиков и т.д. Также такая методика позволяет формировать волноводы с асимметричной формой поперечного сечения канала, на основе которых можно создавать пассивные преобразователи плоскости поляризации оптического излучения [13-14].

Проведен расчет влияния времени заглубления на форму частично заглубляемого канала, то есть на профиль показателя преломления волновода. На рис.3(г^) показано пространственное распределение показателя преломления в сечении ионообменного Ag+-волновода. На рис.3(Ь^) также показаны силовые линии электростатического стимулирующего поля, полученные путем решения уравнения электростатики. Маска закрывает левую половину волновода. Заглубляемый волновод изготовлен термическим обменом Ag+•oNa+ в течении 15 минут при указанной выше температуре, ширина щели в маске 2 мкм. На рис.3^) показано распределение показателя преломления исходного волновода, подвергаемого заглублению, на рис.3(Ь^) показаны распределения

показателя преломления после заглубления в течение 250 с, 900 с, и 1500 с, соответственно. Стимулирующее напряжение 350 В, температура - 380 С0.

5 10

15 .......... *..........................

20 -10 0 1 0 20

20 -10 0 10 20

-20 -10 0 10 20

20 -10 0 и 20

X, \1КМ

Рис. 3. - Контуры распределения показателя преломления волновода, формируемого путем селективного заглубления ионообменного Ag+-волновода, совмещенные с силовыми линиями электростатического поля

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

фундаментальными модами с соразмерными поперечными компонентами электрического и магнитного полей. При подборе определенной длины волновода его можно использовать в качестве четвертьволновой пластинки, преобразующей линейно-поляризованное ТЕ- или ТМ-излучение в излучение круговой поляризации. Аналогично, на основе такого волновода можно создать волноводный ТЕ-ТМ преобразователь поляризации излучения.

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

Благодарности

Работа выполнена при поддержке гранта РФФИ 16-48-230392 р_а "Новые методы физико-математического моделирования многоэтапных процессов формирования оптических микро- и наноструктур в стеклах".

Литература

1. Нуреев И.И. Сенсорные пассивные оптические сети и ключевые вопросы применения в них волоконных брэгговских решеток // Инженерный вестник Дона, 2016, №2 URL: ivdon.ru/ru/magazine/archive/n2y2016/3605.

2. Петров С.Н. Современное оптическое производство и некоторые тенденции его развития // Инженерный вестник Дона, 2009, №3 URL: ivdon.ru/magazine/archive/n3y2009/141.

3. Tervonen A., West B.R., Honkanen S. Ion-exchanged glass waveguide technology: a review. Optical Engineering. 2011. Vol. 50. №7. pp. 071107-1^ 0711071-15.

4. Никитин В.А., Яковенко Н.А. Электростимулированная миграция ионов в интегральной оптике. Краснодар: КубГУ, 2013. 245 с.

5. Najafi S. I. Introduction to Glass Integrated Optics. Artech House, 1992. 170 p.

6. Swillam M.A., Morshed A.H., Khalil D.A. 3D design of MMI splitter with ion exchange technology. SPIE Proceedings. 2005. Vol. 5970, pp. G1-G6.

7. Никаноров Н.В., Петровский Г. Т. Стекла для ионного обмена в интегральной оптике: современное состояние и тенденции дальнейшего развития (обзор) // Физика и химия стекла. 1999. № 1. C. 21-69.

8. Doremus R. H. Exchange and diffusion of ions in glass. J. Phys. Chem. 1964. Vol.68. № 8. pp. 2212-2218.

9. Cheng H., Ramaswamy R.V. Simulation of tapered transitions in ion-exchanged channel waveguides. Applied Optics. 1990. Vol. 29. № 8. pp. 11501156.

10. Guntau M., Brauer A., Karthe W., Pobner T. J. Numerica solution of ion exchange in glass for integrated optical components. Lightwave Technol. 1992. Vol. 10. № 3. pp. 312-315.

11. Xu S.L., Huang W.P. Stern M.S., Chaudhuri S.K. Full-vectorial mode calculations by finite difference method. IEE Proceedings (Optoelectrtronics). 1994. Vol. 141. № 5. pp. 281-286.

12. West B.R., Madasamy P., Peyghambarian N., Honkanen S. Modeling of ion-exchanged waveguide structures. Journal of Non-Crystalline Solids. 2004. Vol. 347. pp. 18-26.

13. Jordana E., Ghibaudoa E., Boucharda A., Blanc-Jamonb M.-F., Royerb F., Broquin J.-E. Design of a waveguide with optics axes tilted by 45° and

realized by ion exchange on glass. Proc. of SPIE. 2016. Vol. 9750. Paper 975009.

14. Векшин М.М., Кулиш О.А., Яковенко Н.А. Поляризационные элементы и устройства интегральной оптики. Краснодар: КубГУ, 2016. 240 с.

References

1. Nureev I.I. Inzenernyj vestnik Dona (Rus), 2016, №2 URL: ivdon.ru/ru/magazine/archive/n2y2016/3605.

2. Petrov S.N. Inzenernyj vestnik Dona (Rus), 2009, №3 URL: ivdon.ru/magazine/archive/n3y2009/141.

3. Tervonen A., West B.R., Honkanen S. Optical Engineering. 2011. Vol. 50. №7. pp. 071107-1- 0711071-15.

4. Nikitin V.A., Yakovenko N.A. Elektrostimulirovannaya migratsia ionov v integralnoi optike [Electric field-assisted ion migration in integrated optics]. Krasnodar: KubSU, 2013. 245 p.

5. Najafi S. I. Introduction to Glass Integrated Optics. Artech House, 1992. 170 p.

6. Swillam M.A., Morshed A.H., Khalil D.A. SPIE Proceedings. 2005. Vol. 5970, pp. G1-G6.

7. Nikanorov N.V., Petrovskii G.T. Phisica i chimiya stekla. 1999. № 1. pp. 21-69.

8. Doremus R. H. J. Phys. Chem. 1964. Vol.68. № 8. pp. 2212-2218.

9. Cheng H., Ramaswamy R.V. Applied Optics. 1990. Vol. 29. № 8. pp. 11501156.

10. Guntau M., Brauer A., Karthe W., Pobner T. J. Lightwave Technol. 1992. Vol. 10. № 3. pp. 312-315.

11. Xu S.L., Huang W.P. Stern M.S., Chaudhuri S.K. IEE Proceedings (Optoelectrtronics). 1994. Vol. 141. № 5. pp. 281-286.

12. West B.R., Madasamy P., Peyghambarian N., Honkanen S. Journal of Non-Crystalline Solids. 2004. Vol. 347. pp.18-26.

13. Jordana E., Ghibaudoa E., Boucharda A., Blanc-Jamonb M.-F., Royerb F., Broquin J.-E. Proc. of SPIE. 2016. Vol. 9750. Paper 975009.

14. Vekshin M.M., Culish O.A., Yakovenko N.A. Polarizatsionnye elementy i ustroistva integralnoi optiki [Integrated-optic polarization elements and devices]. Krasnodar: KubSU, 2016. 240 p.