close

Вход

Забыли?

вход по аккаунту

?

Многослойная оптика

код для вставкиСкачать
Спецкурс
Многослойная оптика
Профессор А.В.Тихонравов
1. Введение
1.1 Общие сведения
Общий вид
многослойного покрытия
тонкие слои
d
Подложка:
d – от
подложка
долей мм до нескольких см
Поперечные размеры от долей мм до метров
Слои: толщина от нескольких мм до нескольких мкм
1.2 Области применения
-
Просветляющие покрытия для оптики
Зеркала (R до 99.9995% - лазерные гироскопы)
Фильтры (полосовые, отрезающие, ... )
Поляризаторы, светоделители, ...
Элементы WDM технологий
Защитные покрытия (например, Euro)
Декоративные покрытия (Swarowski)
Архитектурные покрытия, стекла автомобилей, ...
Спектральные диапазоны
1 мкм = 103 нм = 104 Å
200
400
700
10.6
мкм
ВУФ
ДУФ
УФ
Видимый
Ближний
ИК
Дальний
ИК
1.3 История развития
Первое покрытие (просветление оптики) – Германия 30-ые годы
Первая монография – А.Власов, Просветление оптики, 1946
Начало бурного развития – после появления лазеров
С начала 90-х годов – широчайшее использование, от новейших
технологий до бытовых приложений
Литература:
1. М.Борн, Э.Вольф Основы оптики, М., 1973
2. A.Thelen, Optical Interference Coatings, N.Y., 1989
3. Sh.Furman, A.Tikhonravov, Basicas of Optics of Multilayer Systems,
Editions Frontiers, 1992
1.4 Немного о технологии
Ранние процессы: sol-gel, химическая обработка поверхности
С середины 40-х годов: напыление в вакууме
До начала 80-х годов: различные модификации PVD
Последнее десятилетие: десятки новых высокоэнергетических
процессов
Основные тонкопленочные материалы: оксиды, фториды,
металлы
1.5 Экспериментальные методы
исследования тонких пленок и покрытий
Спектральная фотометрия (начиная с ВУФ)
Спектральная эллипсометрия (видимая область и ближний ИК)
Атомная силовая микроскопия
Электронная микроскопия
......
1.6 Основные задачи
Расчет и исследование спектральных коэффициентов
заданного покрытия.
Синтез (проектирование) покрытий с заданными
спектральными свойствами.
Исследование параметров тонких слоев по
спектральным данным.
Исследование параметров многослойных покрытий по
спектральным данным.
2. Спектральные характеристики
сплошной среды
2.1. Вывод основных уравнений
Рис.1. Модель слоистой среды
Уравнения для поля в слоистой среде
E 1 H ,
c t
E , H
E
4
H E
c t
c
- электрический и магнитный векторы.
Рассмотрим монохроматическую волну с частотой :
E E exp( i t ),
H H exp( i t )
Уравнения для поля примут вид:
4
E i H , H i E E
c
c
c
Введем обозначения
2 c
k 2
- длина волны в вакууме
c
~ i 4
c
- волновое число
- комплексная диэлектрическая проницаемость
Уравнения для поля запишутся в виде:
E ik H ,
~
H ik ( z ) E
Разложение на S- и P-волны.
Рис.1.2. Ориентация электрического и магнитного векторов в S- и P-случаях.
Покоординатная форма уравнения поля (S-случай)
E x
Hx 0 ,
H z
y
z
H y
z
ikH y ,
ik ( z ) E x ,
E x
y
H z
x
ikH z ,
0 ,
H y
x
0
Из 2-го, 3-го и 4-го уравнения следует
2Ex
y 2
2Ex
z 2
k 2 ~ ( z ) E x 0
Ищем E x в виде (метод разделения переменных)
E x ( y , z ) u ( z ) g( y )
Получим
1
2
d g
g( y ) d y 2
d 2u
2~
2 k ( z ) u ( z ) const
u( z) dz
1
Постоянную в м.р.п. положим равной k 2 2. Тогда
d 2g
k 2 2 g ( y ) 0
dy 2
Решение этого уравнения возьмем в виде
g( y ) exp( ik y )
При этом E x запишется так:
E x ( y , z ) u ( z ) exp( ik y )
Из 3-го уравнения в покоординатной записи уравнений поля следует что
H Z ( y , z ) u ( z ) exp( ik y )
Из 2-го уравнения в покоординатной записи следует что
такую же зависимость от y как и E x . Положим
Hy
имеет
H y v ( z ) exp( ik y )
где ( z ) - неизвестная функция (знак минус взят для удобства).
В целом имеем следующее представление для поля
E u ( z ), 0 , 0 exp( ik y ) , H 0 , v ( z ), u ( z ) exp( ik y )
При этом 1-ое, 3-е. 5-е и 6-е уравнение в покоординатной записи уже
выполняются.
Из 2-го и 4-го уравнений получаем
du ikv , dv ik ~ ( z ) 2 u
dz
dz
Это и есть основные уравнения, связывающие поля
E
и H в S–случае.
Для определения постоянной a рассмотри поле падающей волны
E E A exp( i k r i t )
Волновой вектор в однородной среде равен
k k l , l 0, sin a , cos a a
Здесь l - единичный направленный вектор.
Падающее поле в координатном виде
E E a exp ik a ( z cos a y sin a ) i t С другой стороны зависимость поля от
y есть exp( ikay )
откуда
a sin a
Поскольку зависимость поля от y должна быть одинакова во всех
средах, то
a sin a s sin s
Это соотношение есть закон Снеуллиуса
Покоординатный вид уравнений поля
E y
E z E y
E z
ikH x ,
0,
0,
y
z
x
x
H x
H x
Ex 0 ,
ik ~ ( z ) E y ,
ik ~ ( z ) E z
z
y
Ищем H x
и E x в виде
H x ( y , z ) v ( z ) exp( ik y )
,
E y ( y , z ) u ( z ) exp( ik y )
Из 5-го уравнения следует, что
E z ( y , z ) ~ v ( z ) exp( ik y )
(z)
Это означает, что
E 0 , u ( z ), ~
v ( z ) exp( ik y ) ,
(z)
H v ( z ), 0 , 0 exp( ik y )
Подставляя в уравнения для поля, получаем основные уравнения в P-случае
du ik 1 2 v ,
~(z) dz
dv ik ~ ( z ) u
dz
2.2 Амплитудные коэффициенты отражения и пропускания
S-случай
Из основных уравнений получаем, что во внешней среде
d 2 u k 2 ( 2 ) u 0
a
dz 2
Положим
qa =
2
ea - a =
e a cos g a
(a =
e a sin g a )
Во внешней среде есть два решения:
u ( z ) ~ ex p( ikq a z ) ( E ) x ~ ex p( ikq a z i t ) - п адаю щ ая волн а
u ( z ) ~ ex p( ikq a z ) ( E ) x ~ ex p( ikq a z i t ) - отраж ен н ая волн а
Из основных уравнений в случае падающей волны
u ( za ) E A ,
v( za ) qa E A
и в случае отраженной
u ( za ) E R ,
EA , ER
( z a ) qE R
- амплитуды тангенциальных компонент электрического
вектора падающей и отраженной волн на внешней границе
На границе с подложкой
u (0 ) E T ,
где
v (0 ) q s E T
q s s 2 s cos s
Амплитудные коэффициенты пропускания и отражения задаются
соотношениями
t
ET
,
EA
r
ER
EA
Положим E T 1 . Отсюда получим граничные условия для основных
уравнений поля:
u ( 0, k ) 1 ,
Из непрерывности u , v
v ( 0, k ) q s
на внешней границе следует, что
u ( za , k ) E A E R ,
v( za , k ) qa ( E A E R )
Отсюда
ER qau ( za , k ) v( za , k )
,
2qa
EA qau ( za , k ) v( za , k )
2qa
Окончательно для амплитудных коэффициентов получаются следующие
выражения
t (k ) 2qa
,
qau ( za , k ) v( za , k )
r (k ) qau ( za , k ) v( za , k )
qau ( za , k ) v( za , k )
P-случай
Из основных уравнений поля следует, что для падающей волны
u (za ) E
,
A
v(za ) qa E
A
,
а для отраженной
u(za ) E R ,
где
qa a
a 2
v( za ) qa E R
,
a
cos a
Аналогично для прошедшей волны
u (0 ) ET ,
где
qs S
2
v (0 ) q S ET
cos S
S
S
t,r
Как и в S-случае определяем
компонент электрического вектора.
E A , E R , ET
через отношения тангенциальных
- тангенциальные компоненты электрического вектора
падающей, отраженной и прошедшей волн.
Положим E T 1 . Отсюда следуют граничные условия для основных
уравнений поля:
u (0 , k ) 1 ,
v (0 , k ) q S
На внешней границе
u ( za , k ) E A E R ,
v( za , k ) qa ( E A E R )
Откуда получаем
t (k ) 2qa
,
qau ( za , k ) v( za , k )
r (k ) qau ( za , k ) v( za , k )
qau ( za , k ) v( za , k )
Замечание: Формулы для и внешне одинаковы в S и P-случае.
dv i k ~ ( z ) 2 u
dz
При этом в S-случае
du ikv ,
dz
и в P-случае
du ik 1 2 v ,
dz
~ ( z ) dv ik ~ ( z ) u
dz
Граничные условия для этих уравнений также внешне одинаковы
u (0 , k ) 1 ,
но
qa qa a 2
a
a 2
v (0 , k ) q S ,
, qs , qs s 2
s
s в S-случае
,
2
,
в P-случае
2.3. Формулы Френеля
Введем показатель преломления n~ ( z ) равенством
~ ( z ) n~ 2 ( z )
Замечание: В общем случае n~
- комплексная величина: n~ n i .
Закон Снеллиуса, выраженный через показатели преломления
n a sin a n s sin s
Параметры q в расчетных формулах
q a n a cos a , q s n s cos s ,
qa na
ns
, qs ,
cos a
cos s
( S - случай )
( P - случай )
В случае наличия только одной границы раздела
t
2qa
q qs
, r a
qa qs
qa qs
Отсюда следуют формулы Френеля:
2 n a cos a
n cos a n s cos s
t
, r a
( S - случай)
n a cos a n s cos s
n a cos a n s cos s
t
2 n a cos s
n cos s n s cos a
, r a
n a cos s n s cos a
n a cos s n s cos a
( P - случай )
2.4. Энергетические коэффициенты отражения и пропускания
c
S
Re E , H * 8
- вектор Пойнтинга
Поток энергии вдоль нормали к слоистой среде определяется величиной S z .
В S– и P-случаях
S z c Re( uv * )
8
Поэтому падающий и отраженный потоки энергии пропорциональны
соответственно Re q a | E A |2 и Re q a | E R |2 .
Отсюда следует, что энергетический коэффициент отражения равен
E
R R
EA
2
r2
Проходящий поток энергии пропорционален
Re q s | E T |2
Отсюда следует, что энергетический коэффициент пропускания равен
T R eq s E T
R eq a E A
2
R eq s R eq a t
2
Из закона сохранения энергии поглощение в среде равно
A 1 R T
В случае поляризованной волны, определяемой углом поляризации R R s sin 2 R p cos 2 , T T s sin 2 T p cos 2 Для неполяризованного света при наклонном падении
R 1 R s R p ,
2
T 1 T s T p 2
2.3. Матричный метод расчета
Рис.3. Многослойное покрытие состоящее из m однородных слоев
Вывод расчетных формул в S-случае
Уравнения для поля в j–ом слое имеют вид
du
d
2
2
ikv ,
ik n~ j u
dz
dz
Откуда следует уравнение
2
d u
dz
2
k
2
n~
2
j
2
u 0
Определим угол распространения волны в j–ом слое ~
n j sin j
n a sin j
равенством
a
При этом уравнение для u запишется в виде
2
d u
dz
2
k~
n j cos j
2
u0
Это уравнение имеет два линейно независимых решения
u 1 ( z ) cos k ~
n j cos j
z z j1 , u 2 ( z ) sin k ~n j cos j z z j1 Общее решение в слое записывается в виде
u ( z ) c 1 cos k ~
n cos z z j 1
c 2 cos k ~n cos z z j1 Полагая z z j , находим константы c 1 и c 2 :
c1 u ( z j 1 ), c 2 i ( z j 1 )
n~ j cos j
Дифференцируя выражение для u ( z ) находим v ( z ) . Далее,
полагая z z j , получаем
u ( z j ) u ( z j 1 ) cos j i v ( z j 1 )sin j ,
qj
v ( z j ) i q j u ( z j 1 )sin j v ( z j 1 ) cos j ,
где
j k n j cos j d j , d j z j z j 1 , q j n j cos j
.
Введём в рассмотрение матрицу
M
j
co s j
i q sin j
i
sin qj
co s j
j
j
Эта матрица называется характеристической матрицей слоя.
С её помощью можно записать
u
v z z
M
j
j
u
v zz
Переходя от слоя к слою, получим, что
u
v
где
M
z z
M
a
u
v
z 0
M M m M m 1 ... M 1
- характеристическая матрица слоя.
j 1
Используя начальные условия
u (0, k ) 1 ,
v (0, k ) q s
получим, что
u ( z a , k ) m11 m12 q s , v ( z a , k ) m 21 m 22 q s .
Подставляя эти соотношения в общее выражение для r и
t (k ) r (k ) t , получим
2qa
q a m1 1 q s m 2 2 q a q s m1 2 m 2 1
q a m11 q s m 22 q a q s m12 m 21
q a m11 q s m 22 q a q s m12 m 21
Таким образом, матричный метод расчета дает эффективный рекуррентный
алгоритм расчета амплитудных коэффициентов отражения и пропускания.
При нормальном падении матрица j–го слоя равна
M
j
cos j
in j sin i
sin j nj
cos j j
где j kn j d j - фазовая толщина j–го слоя.
При наклонном падении n j заменяется на эффективный показатель
преломления
q j = n j cos qj =
nj -
2
2
- в S-случае
2
nj
cos j
j
nj
2
nj
- в P-случае
2
При этом фазовая толщина слоя в S- и P–случае одинакова и равна
j = kn j cos j d
j
.
2.4 Адмитансный метод расчета
Отношение тангенциальных компонент магнитного и электрического
векторов называется адмитансом:
A( z, k ) ( z, k )
u ( z, k )
Импеданс вводится соотношением
1
Z ( z, k ) A( z, k )
u ( z, k )
( z, k )
Амплитудный коэффициент отражения выражается через
адмитанс на внешней границе покрытия:
r(k ) =
q a - A ( za , k )
q a + A ( za , k )
q a - эффективный показатель преломления внешней среды.
Дифференциальные уравнения для адмитанса
В S-случае
dA
dz
2
2
= ik { ~ ( z ) - - A }
В P-случае
ik ( z ) dz
dA
2
2 1 A (z) Начальное условие задается на границе с подложкой
A ( 0, k ) q s
,
где q s - эффективный показатель преломления подложки
( q s различается в S- и P–случаях).
Рекуррентные формулы адмитансного метода расчета.
Aj iq j sin j A j 1 co s co s j i
qj
A j 1 sin j
, A0 q s
j
A j - адмитанс на правой границе j–го слоя,
qj
j
- эффективный показатель преломления j–го слоя,
2 kn~ j cos j d
j
- фазовая толщина j–го слоя.
Амплитудный коэффициент отражения выражается через адмитанс
на внешней границе последнего (m-го) слоя
r q a Am
q a Am
Фазовая плоскость адмитанса
Адмитанс – комплексная величина, которую можно изобразить точкой
на комплексной плоскости (фазовой плоскости адмитанса).
Рекуррентную формулу для адмитанса можно преобразовать к виду
Aj q j iA j 1
i где
ctg A j 1
qj
j
Когда фазовая толщина слоя j изменяется от 0 до , пробегает
значения от до . При этом A j движется по окружности в фазовой
плоскости.
2.5. Однослойное просветляющее покрытие
(случай нормального падения волны)
Фазовая толщина слоя
knd 2
nd
Применим для расчета r адмитансный метод:
A0 n s , A1 n in s ctg ictg ns
,r n a A1
n a A1
.
n
Коэффициент отражения равен нулю, если
A1 na
(фазовая траектория адмитанса начинается с n s и заканчивается n a )
Это условие выполняется если
n
2
nans , 2
nd 4
Второе условие связывает толщину слоя с длиной волны, на которой
достигается полное просветление.
Просветляющие покрытия (общие замечания)
Как правило не существует материалов с
Типичные значения: n s 1.52 , 1.46 , ...
n
nans
.
n 1.38 , 1.46 , 1.65 , 2 .05 , ...
Поэтому однослойное покрытие не может обеспечить просветление даже
на одной длине волны.
Двухслойные покрытия могут обеспечить просветление на одной длине волны
(достаточно для лазеров).
Для большинства других приложений необходимо просветление в широком
диапазоне длин волн. Этого можно достичь лишь с помощью многослойных
покрытий.
2.7. Формальная запись многослойного оптического покрытия.
Замечание 1: Во всех формулах мы имеем только d
.
Поэтому при пропорциональном изменении всех толщин слоев спектральные
характеристики лишь сдвигаются по спектру.
Замечание 2: В случае нормального падения физические толщины слоев входят
в формулы лишь в сочетании nd .
Величина nd
называется оптической толщиной.
Выберем и фиксируем некоторую длину волны 0 (основная длина волны).
Измерим оптическую толщину слоя в долях 0 . Соответствующая единица
4
измерения называется QWOT.
Пример: 0 600 nm , 1 Q W O T
соответствует оптической толщине 150 nm.
Рассмотрим покрытие состоящее из слоев с двумя чередующимися показателями
преломления n L и n H .
Оптическая толщина 1 Q W O T для слоев с
для слоев с
nH
- H
n L обозначается L ,
.
Пример: Двухслойное покрытие
n s 1.52
Пусть
Тогда
n H 2 .0
n L 1.38
n a 1.00
d H 150 nm
d L 200 nm
0 600 nm
n H d H 300 nm 2 Q W O T
n L d L 276 nm 1.84 Q W O T
Формульная запись этого покрытия:
или
2H
или 1.84 L
2 H 1.84 L
Замечание : Если материалов больше чем два, по L и H приписываются тем,
у которых наименьший и наибольший показатель преломления,
остальным приписываются буквы M , A , B , ...
Пример: 4 H 0.5 LM H
слоев 4, 0.5 , 1
- четырехслойное покрытие с оптическими толщинами
и 1 Q W O T с показателями преломления n H , n L , n M , n H
Замечание : Повторяющиеся группы слоев объединяются в круглые скобки и число
повторений дается показателем степени этих скобок.
Примеры:
4
H LH LH LH LH ( H L ) H
3
H LH LH L 4 H LH LH LH ( H L ) 4 H ( LH )
3
.
2.8. Четвертьволновые диэлектрические зеркала.
Структура: чередующиеся показатели преломления
nH и
nL ,
все слои имеют одинаковую оптическую толщину 1 Q W O T
(при этом основную длину волны 0
центральной длиной волны зеркала).
называют также
Характерные спектральные свойства:
Зона высокого отражения в окрестности 0
Более узкие зоны высокого отражения в окрестности 0 3 , 0 5, ...
Осциллирующий коэффициент пропускания между этими зонами.
Приближенные формулы для четверть волновых зеркал
Первый от подложки слой имеет высокий показатель преломления и
число слоев m нечетно
R (0 ) 1 4
число слоев
m четно
R (0 ) 1 4
Пример:
nans nL nLnH nH m
ns nL nL nH m
n a 1 .0 0, n s 1 .5 2, n H 2 .3 5, n L 1 .4 5
m 1 1 R 0 .8 8 % , m 1 5 R 0 .1 3 %
Ширина основной зоны высокого отражения
upper
low er
1 2 arcsin 1 2 arcsin , где
Замечание : ширина растет с увеличением
nH nL
nH nL
nH nL
Документ
Категория
Презентации
Просмотров
35
Размер файла
1 090 Кб
Теги
1/--страниц
Пожаловаться на содержимое документа