Квантовые и оптические процессы в твердых телах. Теория и практика. Учебное пособие
.pdfek - единичный вектор в направлении распространения вол-
ны;
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
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(ωt+δs ) . |
(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 |
×ei(δ1 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
