Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Квантовые и оптические процессы в твердых телах. Теория и практика. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
802 Кб
Скачать

ek - единичный вектор в направлении распространения вол-

ны;

eE , eH - единичные векторы в направлении E и H , r - ра- диус-вектор из начала координат в точку, где рассматриваются колебания векторов E и H ;

(K × r ) - скалярное произведение волнового вектора K и ра-

диус-вектора r .

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

r r r

+

D

, (I)

r r

B

, (II)

(56)

[ÑH ]= J

[ÑE]= -

 

 

t

 

 

t

 

 

(ÑD)= ρ , (III)

(ÑB)= 0

(IV)

(57)

дополненных двумя уравнениями связи для однородной изотропной среды:

 

 

 

D = ε ×ε0E ,

 

B = μ × μ0 H ,

(58)

и законом Ома в дифференциальной форме:

 

 

 

 

 

 

 

J = σ × E .

 

(59)

 

r

r

r

 

r

 

 

Здесь

Ñ = i

 

 

+ j

 

+ k

 

 

-векторный дифференциаль-

 

x

y

z

 

 

 

 

 

 

 

 

 

ный оператор Гамильтона “набла”;

i , j , k

- единичные векто-

ры в направлении осей x, y, z; D - вектор электрического смещения или, что то же самое, вектор электростатической индукции,

измеряемый в Кл/м2; B - вектор магнитной индукции ([B] = Тл =

=H/(А×м) = Вб/м2); J - вектор плотности тока проводимости ([J] = А/м2); σ - удельная электропроводность среды, имеющая ед. измерения Ом-1×м-1; ρ - объемная плотность заряда ([ρ]=Кл/м3); ε0

=8,85×10-12 Ф/м - электрическая постоянная; μ0 = = 4π×10-7 Гн/м -

магнитная постоянная; ε, μ - относительные диэлектрическая и магнитная проницаемости среды.

101

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

εε0

E =

μμ0

H .

(60)

В самом деле, подставив (54) и (55) в (56) с учетом (57) и (58), получим (60) (более подробно см. приложение 1). Учтем, что

для диэлектриков μ = 1, тогда из (60) следует:

 

 

 

H0 = kS × n1 × E0 p ,

H1 = kS × n1 × E1p ,

H2 = kS × n2 × E2 p ,

(61)

 

 

 

kS

=

 

 

 

 

 

 

 

 

где коэффициент равен:

 

 

ε0 μ0

и имеет единицу измере-

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ния (Ф/Гн)

 

=Ом-1,

n1 =

 

ε1μ1 и n2

=

 

ε2μ2

- абсолютные по-

2

казатели преломления сред 1 и 2.

 

 

 

 

 

 

 

 

Разделив (52) на cosϕ1 , получим:

 

 

 

 

 

 

 

 

 

 

E0 p - E1p

 

= E2 p (cosϕ2 / cosϕ1) .

(62)

Подставив (61) в (53) и разделив на kS ×n1 , имеем

 

 

 

 

 

 

E0 p + E1p

= E2 p (n2 / n1) .

 

 

(63)

Исключив из уравнений (62) и (63) E2 p

и учитывая (46),

находим амплитудный коэффициент отражения Френеля Rp :

R

p

= (E

/ E

) =

n2 cosϕ1 - n1 cosϕ2

.

 

(64)

 

 

 

1p

0 p

 

 

 

n1 cosϕ2

+ n2 cosϕ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Складывая уравнения (62) и (63) и учитывая (47), получаем

амплитудный коэффициент пропускания Френеля:

 

Tp

= (E2 p

E0 p )=

 

 

 

 

2 × n1

× cosϕ1

 

.

(65)

n1

× cosϕ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ n2 ×cosϕ1

 

Пользуясь законом преломления света:

 

 

 

 

 

 

 

 

n1 ×sinϕ1 = n2 ×sinϕ2

 

 

(66)

можно выразить входящий в (64) и (65) косинус угла преломления:

102

(64/)

cosϕ2 =

 

.

 

1- (n1 n2 )2 sin2 ϕ1

(67)

Таким образом, коэффициенты Френеля Rp для отражения света от границы раздела сред 1, 2 и пропускания Тр через эту границу зависят только от абсолютных показателей преломления сред и угла падения ϕ1 света на образец.

Отметим, что с использованием (66) формула (64) может

быть написана в более компактном виде:

Rp = tg((ϕ1 − ϕ2 )) . tg ϕ1 + ϕ2

При ϕ1 + ϕ2 = 90о имеем tg(ϕ1 + ϕ2 )= ∞ и Rp → 0, то есть амплитуда отраженного луча равна нулю. Соответствующий угол

падения

ϕ1 называется

углом Брюстера ϕБр. Так как

 

sinϕБр

 

)=

n2

, то tgϕБр

=

n2

 

sin(90o − ϕ

 

 

 

.

 

Бр

n1

n

 

 

 

 

 

1

 

Б. Вектор E перпендикулярен плоскости падения света p, то есть лежит в плоскости s (рис. 52, б).

Рассмотрев проекцию вектора E на ось у и вектора H на ось х (рис. 52, б) и рассуждая аналогично п. А, получим:

E0s + E1s = E2s , H0 cosϕ1 - H1 cosϕ1 = H2 cosϕ2 .

После проведения выкладок, аналогичных выполненным в п. А, получаются выражения для амплитудных коэффициентов отражения Rs и пропускания Ts Френеля для s- компоненты:

R =

n1 cosϕ1

- n2 cosϕ

2

,

 

 

(68)

n1 cosϕ1

+ n2 cosϕ

 

 

 

s

 

 

 

 

 

 

 

 

2

 

 

 

 

Ts

=

 

2 × n1 ×cosϕ1

 

 

 

.

(69)

n1

× cosϕ1 + n2 ×cosϕ

2

 

 

 

 

Выражения (64), (65), (68) и (69) называются формулами Френеля. Из (64) и (68) следует, что амплитуды р- и s- компонент

103

световой волны после отражения от диэлектрика уменьшаются по сравнению с амплитудами падающей волны.

В. Падающий свет поляризован в плоскости, составляющей угол 45о с плоскостью p.

В этом случае вектор E1 отраженной волны совершает колебания в плоскости, наклон которой определяется отношением р- и s- компонент E1 (рис. 53).

Рис. 53.

Поворот

плоскости

поляризации плоскополя-

ризованного света, отраженного от диэлектрика. E0 - вектор напряжен-

ности электрического поля падающей волны, E1 - вектор напряжен-

ности электрического поля отраженной волны, ψ - угол наклона плоскости поляризации отраженной волны относительно плоскости S

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

tgψ = E1p E1s .

(70)

104

Так как в падающей волне Еор = Еos (см. рис. 53), то, разде-

лив (46) на (47), получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rp

 

Rs = E1p E1s = tgψ .

 

 

 

 

 

 

 

 

(71)

 

Подставив в (71) выражения для Rp

(64) и Rs (68), с уче-

том (67), находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

 

 

 

 

 

 

öæ

 

 

 

 

 

 

 

 

 

 

 

 

 

ö

 

 

 

 

æ

 

n1

ö

2

 

 

 

 

 

 

 

 

æ

n1

ö

2

 

 

 

 

 

ç

 

 

 

2

 

 

֍

 

 

 

 

 

 

 

2

 

÷

 

 

 

 

ç

 

÷

 

 

 

 

 

 

 

 

ç

÷

 

 

 

 

 

 

 

çn2cosϕ1

- n1

 

 

sin

 

ϕ1

 

÷çn1cosϕ1

+ n2

 

 

 

 

 

sin

 

ϕ1

÷

 

 

 

 

1 - ç

 

 

 

÷

 

 

1 - ç

 

 

÷

 

 

 

 

 

ç

 

 

è n2

ø

 

 

 

 

֍

 

 

 

 

 

è n2

ø

 

 

 

 

÷

 

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

øè

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

= tgψ .

æ

 

 

 

 

 

 

 

 

 

 

 

 

 

öæ

 

 

 

 

 

 

 

 

 

 

 

 

 

ö

 

 

æ

 

n1

ö

2

 

 

 

 

 

 

 

 

 

æ

n1

ö

2

 

 

 

 

 

 

ç

 

 

 

2

 

 

֍

 

 

 

 

 

 

 

2

 

÷

 

 

 

 

 

ç

÷

 

 

 

 

 

 

 

 

ç

÷

 

 

 

 

 

 

 

çn2cosϕ1

+ n1

 

 

sin

 

ϕ1

 

÷çn1cosϕ1

- n2

 

 

 

 

 

sin

 

ϕ1

÷

 

 

 

 

1 - ç

 

 

 

÷

 

 

1 - ç

 

 

÷

 

 

 

 

 

ç

 

 

è n2

ø

 

 

 

 

֍

 

 

 

 

 

è n2

ø

 

 

 

 

÷

 

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

øè

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

Полученное уравнение может быть решено аналитически

относительно n2:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= n sinϕ

 

æ

1m tgψ

tgϕ

ö

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

2

 

 

1

ç

÷

+1 .

 

 

 

 

 

 

 

(72)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

ç

1 ± tgψ

 

 

1 ÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь выбираются верхние знаки,

если ϕ1

<

ϕ Бр,

то есть

угол падения не превышает угол Брюстера, и нижние - если

ϕ1

> ϕ Бр, где угол Брюстера ϕ Бр находится из условия

 

 

 

 

 

ϕ Бр =

arctgæ n2 ö .

ççè n1 ÷÷ø

Чтобы исключить неоднозначность в определении n2 по формуле (72), угол ϕ1 выбирают заведомо большим, чем ϕ Бр для исследуемой группы материалов (например, стекол). Выбрав угол падения ϕ1 , помещают исследуемый диэлектрик в среду с из-

вестным показателем преломления n1 и измеряют угол ψ. Затем по формуле (72) рассчитывают показатель преломления диэлектрика n2.

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

105

Рис. 54. Упрощенная оптическая схема эллипсометра компенсационного типа. П - поляризатор, О - образец, А - анализатор, s - след плоскости образца, S - след плоскости, перпендикулярной плоскости падения и проходящей через отраженный луч

Он состоит из поляризатора П, то есть поляроида, плоскость пропускания которого составляет угол 45о с плоскостью падения луча Р, и другого поляроида, называемого анализатором А. Перед измерениями плоскость пропускания анализатора А устанавливается параллельно плоскости Р (положение А0 на рис. 55). Вышедшая из поляроида П плоскополяризованная волна с Eop = Eos , отражаясь от образца О, изменяет плоскость поляри-

зации. После этого она попадает в анализатор А. Вращая анализатор А, добиваются, чтобы интенсивность прошедшего через него света стала равна нулю. В соответствии с законом Малюса это происходит, если плоскость пропускания поляризатора А установлена перпендикулярно плоскости колебаний Е1 отраженного луча (положение А1 на рис. 55). Как видно из рис. 55, угол между плоскостями А0 и А1, то есть угол поворота анализатора, равен ψ. На практике применяют более сложные системы для измерения эллипсометрических параметров.

Другой способ определения угла ψ состоит в измерении интенсивностей I1 и I2 р- и s- компонент отраженного луча (рис.56).

106

Рис. 55. Измерение угла ψ компенсационным эл-

липсометром. Е0 и Е1 - положение плоскости

колебаний в падающем и отраженном лучах, Р - плоскость падения луча, Ао и А1 - исходное и конечное положения плоскости пропускания анализатора при измерении

Рис. 56. Упрощенная оптическая схема эллипсометра с двумя анализаторами А1 и А2 с взаимно перпендикулярными плоскостями пропускания

Поставив на пути отраженного луча поочередно анализаторы А1 с плоскостью пропускания P, совпадающей с плоскостью падения луча р, и А2 с плоскостью пропускания S, перпендикулярной плоскости Р (плоскости Р и S на рис. 53) измеряют интенсивности прошедшего света I1 и I2, а tgψ находят из соотноше-

ния:

107

tgψ =

E1p

=

 

I1

 

.

(73)

 

 

 

 

 

E1s

I2

 

Последнее равенство учитывает, что интенсивность света I пропорциональна квадрату амплитуды напряженности электрического поля в электромагнитной волне (I ~ E2).

3.2. Отражение света от проводящих сред (металлы, полупроводники)

Если выражения (54) и (55), описывающие плоскую электромагнитную волну, подставить в уравнения Максвелла (56) (I и II), то можно получить следующие соотношения, связываю-

щие векторы E и H с волновым вектором плоской электромаг-

нитной волны (см. приложение 1):

 

[K × E]= ωμμ0H , [K × H ]= -ωεε 0E .

(74)

При выводе (64) считалось, что среда диэлектрическая, то есть σ = 0 и, следовательно, j = 0.

Из (74) следует выражение (60), которое вместе с выражениями (50), (51) положено в основу вывода уравнений Френеля

(64), (65), (68), (69).

Если среда проводящая (σ ¹ 0), то уравнения (74) приобре-

тают несколько иной вид:

 

 

 

 

[K × E]= ωμμ0H , [K × H ]= -ωεε 0E + iσE .

(75)

Систему (74) можно свести к системе (75), если e заменить

на комплексную величину

σ

 

 

 

εc = ε - i

 

,

 

ωε

0

 

 

 

 

где i = -1 - мнимая единица.

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

~ = ε n c .

108

Формулы Френеля, являющиеся следствием (74), (60), (50), (51), остаются справедливыми и для случая отражения от границы раздела с проводящей средой. При этом показатель преломления проводящей cреды n2 и зависящий от него cosϕ2 становятся

комплексными величинами. Из (64) и (68) следует, что тогда и коэффициенты Френеля Rp, Rs будут комплексными величинами.

Выражения, описывающие колебания p- и s- компонент напряженности электрического поля для удобства расчетов записывают в комплексном виде:

Ep (t) = Ep × ei(ωt p ) ,

(76)

Es (t) = Es ×ei(ωts ) .

(77)

Чтобы от (76) и (77) перейти к записи уравнений колебаний в действительном виде (43) и (44) достаточно воспользоваться формулой Эйлера:

eiz = cos z + i ×sin z

(78)

иотбросить не имеющие физического смысла мнимые части (76)

и(77).

Перепишем (76) и (77) в виде:

~

Величины Ep , ми Ep и Es:

~

Ep (t)

= ~ × iωt

Ep (t) Ep e ,

= ~ × iωt

Es (t) Es e .

~

Es называются комплексными амплитуда-

= Ep ×e

iδ p

~

(t) = Es

i

δ s .

(79)

 

, Es

× e

Комплексные коэффициенты Френеля для отражения света от проводящей среды по аналогии с (46), (47) можно представить в виде

~

 

~

 

E1p ×e

iδ1 p

~

~

 

E

×eiδ1s

 

 

E1p

 

 

 

 

E

 

 

Rp

=

~

=

 

 

 

 

;

Rs =

~1s

=

 

1s

 

 

.

(80)

E

 

i

 

E

 

i

 

 

 

E

0 p

 

×e

δ0 p

 

E

 

0s

×e

δ0s

 

 

 

 

 

0 p

 

 

 

 

 

0s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

109

 

 

 

~

 

~

 

 

 

 

Если в падающей волне E0 p = E0s , то

 

 

~

 

~

 

~

 

~

 

 

Rp

E1p

E

E1p

 

 

~

=

~

 

 

~1s

=

~

.

(81)

Rs

E0 p

E0s

E1s

 

 

Подставив в (81) выражение (79), получим основное урав-

нение эллипсометрии

 

 

 

 

 

 

 

 

~

 

E

 

 

 

 

 

R

 

 

 

 

 

 

 

~p

=

1p

×ei1 p −δ1s ) ,

 

 

Rs

 

E1s

 

 

 

 

 

которое с учетом (70) обычно представляют в виде:

 

~

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

~p

= tgψ × ei ,

 

 

(82)

 

Rs

 

 

 

 

 

 

 

 

где D = δ1p -δ1s характеризует разность фаз между р-

и s- со-

ставляющими E в отраженной волне, а

E1p E1s = tgψ характе-

ризует отношение их амплитуд.

Волна, отраженная от проводящей поверхности, имеет не равную нулю или 180о разность фаз между р- и s- компонентами

E (D ¹ 0, 180о). Такая волна является эллиптически поляризо-

ванной, то есть конец вектора E с течением времени описывает в пространстве эллипс.

Практически величину D можно определить, если ввести в

схему эллипсометра два дополнительных анализатора А3 и А4 с плоскостями пропускания, составляющими угол 45о с плоскостью S, и на пути отраженного луча, идущего через анализатор А4, поставить компенсатор С, который вносит сдвиг фаз в 90о между р- и s- компонентами при прохождении через него световой волны

(рис. 57).

Необходимость использования компенсатора при определении угла D можно понять из следующих рассуждений. Пусть при отражении от диэлектрического образца |Rp| = |Rs|, то есть то есть

110

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]