Статья:

Анализ спектров эквивалентных структур на основе пленок ZnSe, As2Se3, BaF2, PbF2, GeTe

Конференция: XXX Международная научно-практическая конференция «Научный форум: технические и физико-математические науки»

Секция: Оптика

Выходные данные
Котликов Е.Н., Новикова Ю.А., Терещенко Г.В. Анализ спектров эквивалентных структур на основе пленок ZnSe, As2Se3, BaF2, PbF2, GeTe // Научный форум: Технические и физико-математические науки: сб. ст. по материалам XXX междунар. науч.-практ. конф. — № 1(30). — М., Изд. «МЦНО», 2020. — С. 35-50.
Конференция завершена
Мне нравится
на печатьскачать .pdfподелиться

Анализ спектров эквивалентных структур на основе пленок ZnSe, As2Se3, BaF2, PbF2, GeTe

Котликов Евгений Николаевич
д-р физ.-мат. наук, профессор, Санкт-Петербургский государственный университет аэрокосмического приборостроения, РФ, г. Санкт-Петербург
Новикова Юлиана Александровна
канд. физ.-мат. наук, доцент, Санкт-Петербургский государственный университет аэрокосмического приборостроения, РФ, г. Санкт-Петербург
Терещенко Георгий Викторович
канд. физ.-мат. наук, доцент, Санкт-Петербургский государственный университет аэрокосмического приборостроения, РФ, г. Санкт-Петербург

 

ANALYSIS OF EQUIVALENT STRUCTURE SPECTRA BASED ON ZnSe, As2Se3, BaF2, PbF2, GeTe FILMS

 

Evgeny Kotlikov

Doctor of Physical and Mathematical Sciences, Professor, St. Petersburg State University of Aerospace Instrumentation, Russia, St. Petersburg

Yuliana Novikova

Candidate of Physics and Mathematics, Associate Professor, St. Petersburg State University of Aerospace Instrumentation, Russia, St. Petersburg

Georgy Tereshchenko

Candidate of Physics and Mathematics, Associate Professor,  St. Petersburg State University of Aerospace Instrumentation, Russia, St. Petersburg

 

Аннотация. Проведен анализ спектров эквивалентных структур, имитирующих в широкой области спектра пленку с требуемым показателем преломления. В эквивалентную структуру входят чередующиеся пленки из двух веществ с высоким и низким показателями преломления. Рассчитаны эквивалентные показатели преломления и спектры эквивалентных структур на основе пленок ZnSe, As2Se3, BaF2, PbF2, GeTe.

Abstract. The spectra of equivalent structures that simulate a film with a required refractive index in a wide spectral region are analyzed. The equivalent structure includes alternating films of two substances with high and low refractive indices. The equivalent refractive indices and spectra of equivalent structures based on ZnSe, As2Se3, BaF2, PbF2, GeTe films were calculated.

 

Ключевые слова: пленки; интерференционные покрытия; эквивалентные структуры; показатели преломления.

Keywords: films; interference coatings; equivalent structures; refractive indices.

 

В оптическом приборостроении широко используются многослойные интерференционные покрытия (ИП). ИП используются в различных целях: просветление оптических элементов, создание оптических фильтров, зеркал, светотоделителей, интерференционных поляризаторов. Для изготовления ИП необходимы оптические материалы с различными показателями преломления. Набор таких материалов ограничен [1-3]. Для расширения технологических возможностей для создания ИП можно использовать многослойные эквивалентные структуры, имитирующие пленку с заданным показателем преломления [4-7]. Это позволяет упростить разработку и реализацию ИП с требуемыми характеристиками. Для замены слоя с определенным показателем преломления в ИП в большинстве случаев достаточно эквивалентной трехслойной структуры. Экспериментальный опыт по созданию ИП различных назначений показал, что условие идентичности показателя преломления лучше выполняется при использовании эквивалентного слоя со структурой LHL, где L – слой с низким показателем преломления nL, а H – слой с высоким показателем преломления nH. В настоящей работе рассмотрена методика формирования эквивалентных слоев и спектров на их основе. Эта методика используется в программе синтеза интерференционных покрытий Film Manager [8].

Возможность использования эквивалентных структур следует из матричного метода расчета спектра ИП [4,5]. В случае нормального падения света на ИП, решая граничную задачу для многослойного покрытия можно получить матричное уравнение связывающее полные тангенциальные составляющие векторов электрического  и магнитного полей  в j и j1 слоях при нормальном падении электромагнитной волны [4]:

                                                                                  (1)

где  – фазовая толщина,  – комплексный показатель преломления слоя, nj – показатель преломления, kj –показатель поглощения, lj – толщина j-го слоя ИП.

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

  .                                                                         (2)

Матрица Msys называется характеристической матрицей интерференционной системы. Пусть mij – матричные элементы характеристической матрицы многослойной системы, тогда

 .                                                                                                       (3)

Характеристические матрицы одиночного слоя и симметричной комбинации пленок имеют равные диагональные элементы. Перемножая две характеристические матрицы А и B слоев с различными показателями преломления nA и nB , и оптическими толщинами hА и hB, можно убедиться в том, что матрица произведения C имеет неравные диагональные элементы. Однако, если матрицу C умножить на матрицу A, то матрица ABA, также имеет равные диагональные элементы. Можно показать, что матрица Msys, являющаяся произведением любой симметричной комбинации ABCDEEDCBA имеет равные диагональные элементы. Матрице Msys можно придать вид матрицы одиночного слоя если ввести следующие величины: эквивалентную фазовую толщину

                                                                                                          (4)

и эквивалентный показатель преломления

                                                                                                            (5)

При этом матрица Msys принимает вид

 .                                                                                          (6)

Так как элементы матрицы Msys зависят от длины волны, то величины ФЕ и NE также являются функциями длины волны. Эти зависимости имеют сложный характер, причем спектральные интервалы в которых ФE и NE принимают действительные значения, чередуются с областями где эти величины не являются действительными. Из выражений (4) и (5) следует, что ФЕ и NE действительны при m11<1. При m11 >1 NE является чисто мнимой величиной, а ФE – комплексной.

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

Используя понятия эквивалентный показатель преломления и эквивалентная толщина, можно любую пленку Р для фиксированной длины волны заменить симметричной комбинацией пленок АВС...СВА. В данном случае эквивалентность пленки Р и комбинации АВС...СВА следует понимать в том смысле, что обе эти комбинации описываются одной и той же матрицей. На рис. 1, в качестве примера, представлены зависимости эквивалентного показателя преломления и коэффициента отражения от нормированной частоты для трехслойных симметричных покрытий типа HLH и LHLЗдесь Н и L – четвертьволновые слои с высоким nH и низким nL показателями преломления, соответственно.

 

Рисунок 1. Зависимости эквивалентного показателя преломления NE и энергетического коэффициента отражения R от нормированной частоты n симметричных комбинаций HLH(1) и LHL(2) на подложке с ns=2,4. (nL=1,35, nH=2,65, оптическая толщина всех слоев h=l0/4).

 

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

Рассмотрим комбинацию, состоящую из трех слоев, выполненных из материалов с показателями преломления nА и nв, фазовыми толщинами ФА и Фв и оптическими толщинами hA и hB. Матричные элементы такой системы имеют вид

                                                                                (7а)

                                                                           (7б)

                                                                                   (7в)

                                                                                 (7г)

Из соотношений (7) можно получить выражения для эквивалентного показателя преломления и фазовой толщины:

                                                                                          (8)

                                                                                       (9)

При конструировании покрытий обычно бывает необходимо по заданным значениям эквивалентных показателя преломления и оптической толщины hE, показателям преломления nА и nв вычислить оптические толщины симметричной комбинации hA и hB.

Из выражений (7) можно получить соотношение

                                                                                           (10)

Знак перед sinФв выбирают таким образом, чтобы все действительные значения NE были положительными. В этом случае

 где k=0,1,2,3, …..                                                               (11)

Таким образом величина Фв известна. Чтобы найти ФА надо решить уравнение (9). Введем обозначения:

                                                                                              (12)

Тогда соотношение (9) примет вид:

                                                                                                     (13)

Для решения уравнения (13) используется универсальная подстановка:

при 2FA=π+2πk, где k=1,2,3, …

Решением уравнения (13) является выражение

                                                                                      (14)

Здесь abc вычисляются по формулам (12).

Если (b2+c2)<a2 то уравнение (13) не имеет действительных решений. На практике более удобно использовать решение при k=0. В этом случае оптические толщины слоев FA и FB имеют наименьшие значения. При проектировании интерференционных покрытий удобнее пользоваться не фазовыми толщинами слоев ФА и ФB, а оптическими hA и hB, измеренными в четвертях масштабной длины волны. Если длина волны измеряется в относительных длинах волн lотн=l/l0, то из (11) - (14) получаем выражения для вычисления оптических толщин

,

,                                                                             (15)

где lE  длина волны, на которой выполняется условие эквивалентности.

В таблице 1 приведены рассчитанные нами значения оптических толщин слоев А и В конструкции АВА, вычисленные при заданных NE и hE=l0/4, по формулам (15). Для слоев А и В использованы три пары веществ с показателями преломления 1,35 и 2,3; 1,35 и 2,65; 2,4 и 4.

Таблица 1.

Рассчитанные оптические толщины слоев hA и hB эквивалентной конструкции АВА с hE=l0/4 при различных значениях эквивалентного показателя преломления nЕ.

nA=1,35; nB=2,3

nA=1,35; nB=2,65

nA=2,4; nB=4

NЕ

hA

hB

NЕ

hA

hB

NЕ

hA

hB

1,35

1,000

0,000

1,35

1,000

0,000

2,4

1,000

0,000

1,4

0,476

0,042

1,5

0,443

0,093

2,6

0,446

0,096

1,5

0,431

0,121

1,6

0,408

0,151

2,8

0,395

0,187

1,6

0,388

0,198

1,7

0,374

0,207

3,0

0,346

0,277

1,7

0,346

0,274

1,8

0,342

0,263

3,2

0,297

0,368

1,8

0,305

0,350

1,9

0,311

0,318

3,4

0,248

0,464

1,9

0,264

0,429

2,0

0,281

0,374

3,6

0,196

0,571

2,0

0,222

0,514

2,1

0,250

0,432

3,8

0,134

0,702

2,1

0,176

0,609

2,2

0,220

0,493

4

0

1

2,2

0,121

0,728

2,3

0,189

0,558

 

2,3

0,000

1

2,4

0,155

0,630

 

2,5

0,117

0,717

2,65

0,000

1,000

 

Анализ данных из таблицы 1 позволяет получить следующие приближенные соотношения для эквивалентной толщины и эквивалентного показателя преломления симметричной комбинации, состоящей из m слоев с показателями преломления ni и толщинами hi

 и                                                                                                   (16)

В процессе анализа и синтеза многослойных периодических систем часто используют трехслойные симметричные комбинации с фиксированным отношением оптических толщин слоев. В этих случаях выражения для NE и FЕ упрощаются. При выборе ФВ=2ФА выражения для эквивалентного показателя преломления NE и эквивалентной фазовой толщины FЕ приобретают вид:

 ,                                                                       (17)

                                                                                                (18)

Из выражения (17) следует, что NE принимает комплексные значения при выполнении неравенства:

                                                                                              (19)

Решением неравенства (19) являются значения фаз овой толщины слоя В удовлетворяющие условию

 где k=0,1,2…                                                                           (20)

Фазовая толщина Фв зависит от длины волны. Следовательно, диапазоны длин волн в которых эквивалентный показатель преломления принимает мнимые и действительные значения будут чередоваться. Центры полос заграждения комбинации 0.5АВ0.5А при 2hA=hB=l0/4 находятся на длинах волн, для которых выполняется условие

 где k=0,1,2,3…                                                                                               (21)

На рис. 2 изображены зависимости эквивалентных показателей преломления NE от относительной частоты n для комбинаций типа 0.5LH0.5L, 0.5HL0.5H, состоящих из пар веществ с показателями преломления n: 1,35 и 2,65; 2,4 и 4; 1,35 и 1,7. Этим парам соответствуют следующие отношения высокого и низкого значений показателей преломления nH/nL: 1,96, 1,67, 1,26. Из рис. 2 видно, что системы 0.5LН0.5L, 0.5НL0.5Н имеют области заграждения в окрестности точек n=1,3,5 ... , а также, чем больше отношение nH/nL, тем шире полосы заграждения, и меньше полосы пропускания.

Рассмотрим зависимости относительной ширины полос пропускания от нормированной частоты для эквивалентных слоев типа 0.5АВ0.5А (см. рис. 2). Отсчет по­лос будем вести от точки l=0. Для первой полосы − , для второй полосы Dl~3, для третей полосы − Dl~1,66, для четвертой Dl~1,4. Таким образом при увеличении частоты происходит уменьшение относительной ширины полос пропускания симметричных трехслойных комбинаций. Самой широкой является первая полоса пропускания. Поэтому для синтеза ИП можно будет пользоваться именно этой спектральной областью.

 

Рисунок 2а. Зависимость эквивалентного показателя преломления NE симметричной комбинаций 0.5HL0.5H от нормированной частоты nНа рис. (а) и (б) 1 − nL=1,35; nH=2,65; (nH/nL=l,96); 2 −  nL=2,4; nH=4; (nH/nL=l,66); 3 −  nL=1,35; nH=1,7 (nH/nL=1,26)

 

Рисунок 2б. Зависимость эквивалентного показателя преломления NE симметричной комбинаций 0.5LH0.5L (б) от нормированной частоты n

 

Рисунок. 3. Зависимости эквивалентных показателей преломления NE от относительной длины волны l для симметричных трехслойных конструкций:

1 – конструкция LHL (ZnSe и Ge), nL=2,4; hL=0,346; nH=4; hH=0,277;

2 – конструкция HLH (Ge и ZnSe), nH=4; hH=0,297; nL=2,4; hL=0,368;

 3 – конструкция LHL(BaF2 и As2Se3), nL=1,35; hL=0,281; nH=2,65; hH=0,374;

4 – конструкция HLH (As2Se3 и BaF2), nH=2,65; hH=0,346; nL=1,35; hL=0,257.

 

Из рис. 3 видно, что условие эквивалентности выполняется точно только на длине волны l=1. При λ>1 эквивалентные показатели преломления принимают значения близкие к тем, которые они принимают при 1. А на λ<1 условие эквивалентности выполняется только на очень узких спектральных участках, разделенных полосами заграждения. На рис. 3 штриховкой выделена полоса заграждения в области l=0,5. В спектрах эквивалентных слоев, которые имеют равные эквивалентные оптические толщины hE, полосы заграждения имеют одинаковое положение.

Поэтому поиск структур просветляющих покрытий (ПП) с фиксированными четвертьволновыми слоями на начальных этапах синтеза удобен тем, что полосы заграждения слоев, замененных эквивалентными, будут иметь одинаковое положение. Так как при толщинах слоев равных l0/4 длинноволновая граница полосы заграждения lstop находится на относительной длине волны близкой к 0,6, при поиске структуры ПП на начальных этапах мы использовали диапазоны длин волн, удовлетворяющих условию l>0,6. При замене слоя эквивалентным не во всех случаях использовалась lE=1. В некоторых случаях для получения ПП с более качественными спектральными характеристиками для каждого слоя подбиралась своя lЕ немного большая или меньшая 1. При замене слоев покрытия эквивалентными конструкциями происходит некоторое сужение области просветления из-за больших значений дисперсии nЕ при l<1. Поэтому на начальных этапах синтеза ПП предпочтительнее получать ширину области просветления несколько большую, чем это необходимо. Например, если после замены слоев на эквивалентные было нужно получить относительную ширину области просветления Dl равную 4, то на начальных этапах проводился поиск ПП с Dl=5.

На рисунке 4 приведены рассчитанные нами зависимости NE(l) для симметричных трехслойных конструкций LHL (PbF2/GeTe/PbF2) с эквивалентными толщинами hE=0,5: 1,0; 1,5 и NE=3 (n(GeTe)=3,54; n(PbF2)=1,65 в диапазоне l=4-10 мкм [3]). Толщины слоев комбинаций были вычислены по формулам (17) при lE=1. Из сравнения спектров следует, что длинноволновая граница полосы заграждения lstop прямо пропорциональна эквивалентной толщине hE и при lЕ=1 выполняется приближенное соотношение: lstop~0.6hE. Чем меньше эквивалентная толщина, тем ближе значения эквивалентного показателя преломления NE к заданному n в требуемом спектральном диапазоне.

Рисунок 4. Зависимости эквивалентных показателей преломления NE от относительной длины волны l для симметричных трехслойных конструкций типа LHL(nL=1,65; nH=3,5) с эквивалентными толщинами hE=0,5 (1); hE=l (2) и hE=1,5 (3).

 

Формулы (15) можно использовать для конструирования не только трехслойных эквивалентных комбинаций, но и для комбинаций с большим числом слоев. Например, для получения пятислойной конструкции слой с показателем преломления n и толщиной h, подлежащий замене эквивалентным слоем делится на два слоя с показателями преломления n и толщинами h/2. После этого каждый слой заменяется эквивалентной трехслойной комбинацией.

Толщины смежных слоев с одинаковыми показателями преломления складываются. В таб. 2. представлены результаты конструирования трех-, пяти- и семислойных эквивалентных конструкций с nЕ=3 и hE=l на основе комбинации GeTe/As2S3/GeTe (n(GeTe)=3,5; n(As2S3)=2,3 в диапазоне l от 5 до 10 мкм.

Таблица 2.

Конструкции эквивалентных структур с различным числом слоев на основе GeTe и As2S3.

 

ИП 3 слоя (HLH)

ИП 5 слоев (HLHLH)

ИП 7 слоев  (HLHLHLH)

№ слоя

n

h, l0/4

n

h, l0/4

n

h, l0/4

1

3,5

0,374

3,5

0,164

3,5

0,106

2

2,3

0,233

2,3

0,163

2,3

0,115

3

3,5

0,374

3,5

0,327

3,5

0,212

4

 

2,3

0,163

2,3

0,115

5

3,5

0,164

3,5

0,212

6

 

2,3

0,115

7

3,5

0,106

 

На рисунках 5 а, б, в, в качестве примера, приведены спектры отражения эквивалентных конструкций HLH, HLHLH, HLHLHLH и пленки с n=3 и h=l на подложке с n3=4. Сравнивая спектры можно сделать вывод, что чем больше количество слоев в эквивалентной конструкции, тем больше степень совпадения спектральных характеристик исходной и эквивалентной конструкций покрытия. Семислойной эквивалентной комбинации вполне достаточно, чтобы с высокой степенью точности имитировать спектральные характеристики однородной пленки. Однако, на практике число слоев покрытия ограничено. Поэтому несмотря на то, что имитационные свойства эквивалентных слоев улучшаются при увеличении количества слоев, авторами в основном используются трехслойные конструкции эквивалентных слоев.

 

Рисунок. 5а. Спектр отражения трехслойной эквивалентной комбинации HLH (1) в сравнении со спектром отражения пленки (4) с показателем преломления n=3 и h=l

 

Рис. 5б. Спектр отражения пятислойной эквивалентной комбинации HLHLH (2) в сравнении со спектром отражения пленки (4) с показателем преломления n=3 и h=l

 

Рисунок 5в. Спектры отражения семислойной эквивалентной комбинации HLHLHLH (3) в сравнении со спектром отражения пленки (4) с показателем преломления n=3 и h=l

 

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

 

Список литературы:
1. Котликов Е.Н., Кузнецов Ю.Н., Лавровская Н.П., Тропин А.Н. Оптические пленкообразующие материалы для инфракрасной области спектра // Научное приборостроение. 2008. №3. С. 32-36.
2. Котликов Е.Н., Терещенко Г.В. Исследование оптических констант плёнок, используемых для синтеза широкополосных просветляющих покрытий // Оптика и спектроскопия. 1997.Т. 82, № 4. С. 653–659.
3. Котликов Е.Н., Новикова Ю.А., Юрковец Е.В. Метод определения оптических констант плёнок на поглощающих подложках //Оптический журнал. 2018. Т. 85. № 10. С. 64-69.
4. Бернинг П.Х. Теория и методы расчета оптических свойств тонких пленок. в кн. Физика тонких пленок. Т.1. М.; Мир,1967. С.91-151.
5. Яковлев П.П., Мешков Б.Б. Проектирование интерференционных покрытий. М.; Машиностроение,1987. - 192с.
6. Котликов Е.Н., Терещенко Г.В. использование халькогенидных соединений для изготовления просветляющих покрытий в средней ИК-области спектра // Оптический журнал. 1997. Т. 64. № 3. С. 110.
7. Котликов Е.Н., Новикова Ю.А., Плехоткина Г.Л., Прилипко В.К.
Проектирование ахроматических просветляющих покрытий // Датчики и системы. 2016. № 1 (199). С. 54-57.
8. Котликов Е.Н., Новикова Ю.А., Коваленко И.И. Программа синтеза и анализа интерференционных покрытий Film Manager. // Научный журнал, информационно-управляющие системы. Политехника. 2015 г. – № 3(76). – С. 51-59.