ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2008, том 18, № 1, c. 3-15
= ИССЛЕДОВАНИЯ, ПРИБОРЫ, МОДЕЛИ -
И МЕТОДЫ АНАЛИЗА
УДК535.5.511: 531.7 © А. И. Семененко, И. А. Семененко
О НОВЫХ ВОЗМОЖНОСТЯХ МЕТОДА ЭЛЛИПСОМЕТРИИ, ОБУСЛОВЛЕННЫХ "НУЛЕВОЙ" ОПТИЧЕСКОЙ СХЕМОЙ. ЭЛЛИПСОМЕТРИЯ РЕАЛЬНЫХ ПОВЕРХНОСТНЫХ СТРУКТУР.
11. ПРОБЛЕМА ПОВЫШЕНИЯ ТОЧНОСТИ В "НУЛЕВОЙ" ЭЛЛИПСОМЕТРИИ. ОПРЕДЕЛЯЮЩАЯ РОЛЬ ИЗМЕРИТЕЛЬНЫХ
КОНФИГУРАЦИЙ ПРИБОРА
Работа посвящена изучению измерительных конфигураций прибора, связанных с положением "быстрой" оси фазового компенсатора эллипсометра. Показано, что они играют огромную роль в повышении точности экспериментального определения поляризационных углов Д и Рассмотрены частные случаи (Д = 0, п, п/2, 3п/2), допускающие полное аналитическое рассмотрение. Для этих случаев получены и исследованы выражения, определяющие вторые производные от интенсивности светового пучка на выходе прибора по угловым положениям поляризатора и анализатора. В результате выявлены измерительные конфигурации, которым отвечают максимальные значения указанных производных, обеспечивающие достаточную выраженность минимума интенсивности, а значит, и необходимую точность в экспериментальном определении углов Д и Сделан вывод, что процесс измерения углов Д и ¥ в определенных ситуациях, связанных с малыми значениями угла может быть разделен. В то же время отмечено, что во многих случаях необходимо стремиться к выбору некоторой общей оптимальной измерительной конфигурации, обеспечивающей достаточное разрешение по обоим поляризационным углам. На основе результатов, относящихся к частным случаям, проанализирована измерительная ситуация в окрестности угла Брюстера. Обсужден общий случай произвольных значений поляризационных углов Д и Кроме того, в работе обсуждены вопросы, требующие продолжения анализа измерительных конфигураций прибора. В частности, сделан вывод о необходимости детального изучения роли измерительных конфигураций в "нулевой" эллипсометрии анизотропных сред.
ВВЕДЕНИЕ. ПОСТАНОВКА ЗАДАЧИ
Настоящая работа является непосредственным продолжением работы [1], посвященной анализу зонных соотношений "нулевой" эллипсометрии для случая произвольной ориентации "быстрой" оси компенсатора относительно плоскости падения. Фиксируя в эллипсометрическом эксперименте то или иное положение "быстрой" оси, мы определяем тем самым конкретную измерительную конфигурацию прибора. Как известно, существуют такие области значений ^ и А, которым соответствует слабая выраженность (по одному или двум параметрам гашения) минимума интенсивности светового пучка на выходе анализатора. В этом случае может существенно понизиться точность измерения положений гашения, а значит, и точность определения ^ и А . Особенно сильно это проявляется в окрестности угла Брюстера, которой отвечают малые значения поляризационного угла ^ . Есть все основания считать, что можно существенно усилить выраженность минимума интенсивности светового пучка на выходе анализатора, подбирая необходимую измерительную
конфигурацию. Данной проблеме и посвящена настоящая работа. Для полного решения задачи надо исходить из выражения для интенсивности световой волны на выходе анализатора. В соответствии с этим изложение ведется по следующему плану:
• выражение для интенсивности светового пучка на выходе прибора и общая схема его исследования;
• определение положений гашения поляризатора и анализатора через поляризационные углы ^ и А, возможности моделирования процесса прохождения светового пучка через оптический тракт прибора;
• выбор измерительных конфигураций прибора, обеспечивающих оптимальную точность в экспериментальном определении положений гашения оптических элементов; изучение частных случаев, допускающих полное аналитическое рассмотрение; общая схема исследования для произвольного случая;
• обсуждение вопросов, требующих продолжения анализа измерительных конфигураций прибора.
1. ВЫРАЖЕНИЕ ДЛЯ ИНТЕНСИВНОСТИ СВЕТОВОГО ПУЧКА НА ВЫХОДЕ ПРИБОРА.
ОБЩАЯ СХЕМА ИССЛЕДОВАНИЯ ХАРАКТЕРА МИНИМУМА ИНТЕНСИВНОСТИ
Прежде чем переходить к выражению для интенсивности светового пучка на выходе анализатора, сделаем одно существенное уточнение. Оно касается переменных, определяющих положение оптических элементов в пределах ]-й измерительной зоны. Для рассматриваемой здесь оптической системы РК^А (поляризатор—компенсатор— образец—анализатор) измерительная зона определяется (см. [2-4]) сочетанием типов ориентации анализатора (^^) и компенсатора (ц.) относительно плоскости падения
(Z j, h j),
(1)
где Zj и hJ — единичные знакопеременные параметры, связанные с номером зоны j следующими соотношениями:
Z, = (-1)j+1, h, =
\ 1, j = 1,2; [-1, j = 3,4.
ляется переопределением. В такой ситуации их надо обозначать как уа, у к и ур . В то же время положения гашения оптических элементов по-прежнему будем обозначать через у(аз), у kj) и ур). Поскольку "быстрая" ось компенсатора в
пределах j-й зоны обычно фиксируется (в том числе и в настоящей работе), то ее положение, независимо от рассматриваемого случая, будем определять через у 'к) или же через угол в, (см. [1, 5])
0 j = у j - Р.
(3)
(2)
Знак каждого из этих параметров определяет положительный (+) или отрицательный (-) тип ориентации соответствующего оптического элемента. Положительно (отрицательно) ориентированный элемент (точнее, направление пропускания анализатора, или "быстрая" ось компенсатора) отклоняется от плоскости падения на угол, не превышающий 90° и принадлежащий независимо от типа ориентации интервалу (0, 90°). Во всех предыдущих работах обозначение этого угла содержало верхний индекс (j), указывающий на принадлежность j-й измерительной зоне. Для анализатора — это угол у (аз) , а для компенсатора — укз) . Что касается поляризатора, то он не участвует в определении измерительных зон системы РК8А, и его положение относительно плоскости падения определяется углом, который не ограничен никаким интервалом и имеет знак, связанный с выбранным направлением отсчета положительных угловых значений, одинаковым для всех оптических элементов. Этот угол обозначался как ур),
т. е. опять-таки с верхним индексом. Но в этих работах речь шла в основном о положениях гашения оптических элементов в пределах j-й измерительной зоны, поэтому наличие соответствующего верхнего индекса было оправдано. В общем же случае, когда при выбранной измерительной зоне эти углы свободно варьируются (для анализатора и компенсатора в пределах интервала (0, 90°)), наличие верхнего индекса уже неоправданно и яв-
Это совершенно естественный подход, т. к. положения гашения анализатора и поляризатора, как и фиксированное положение "быстрой" оси, зависят от номера зоны.
Выражение для интенсивности светового пучка на выходе оптической системы PKSA определим, используя полную комплексную амплитуду волны E0A) на выходе анализатора в пределах j-й зоны. Эта амплитуда определяется с использованием аппарата матриц Джонса [2, 3] и должна содержать свободные углы для анализатора и поляризатора и фиксированный угол для компенсатора. В силу изложенных выше соображений это должны быть углы ya и gp (без верхнего индекса) и угол y k) (0j). Воспроизведем здесь приведенное в работе
[1] выражение для полной комплексной амплитуды на выходе анализатора. При этом перейдем к новым обозначениям E(01) ® E0A) и E0O) ® E0P), используем общую (исходную относительно ya) форму записи [2, 3], а также учтем переход к указанным свободным углам:
E0A) ~ E0P) (j jRp cosy a + b2 jR sinya), (4)
Ь j = [(1 + P) - (1 - P) sin 20j ] cos gp + + hj (1 - P) cos 20j sin gp,
b2 j = [(1 + p) + (1 - p) sin 20j ]sin gp + h j (1 - p )cos20j cos gp,
(5)
(6)
где Е(0Р) — полная комплексная амплитуда волны на выходе поляризатора; Яр и — комплексные
амплитудные коэффициенты отражения Френеля для р- и 8-волн. В формулы (4)-(6) входят параметры Zj и , сочетание которых (см. (1)) определяет j-ю измерительную зону. Как и в работе [1], используется идеальный фазовый компенсатор с комплексным параметром р = /ехр(-/5) . Для дальнейшего в выражении (4) целесообразно вы-
делить поляризационные углы ^ и А :
^ ~ (С А , 1ё ^ 008 уа + Ь2; 8шу„) . (7)
С23 = (1 -/2)81п2^..
т(3 )
(18)
Функция 1А (У , уа) достигает минимума (ну-
Используя (7), (5) и (6), найдем выражение для интенсивности светового пучка на выходе оптической системы РК^А:
г( 3)
7( А)
левого значения) в точке (ур', уа3)), образованной положениями гашения поляризатора и анализатора (при фиксированном положении в3 "быстрой" оси),
А Р)
|Д| ^ 0082уа 182^ + 0,81П2уа
1А )(ур.), у а3)) = 0.
+2С. 8пуа 008 уа ^ Ч(И. 008 А + д 81п А)], (8) В этой точке выполняются условия
или
А 3)
7( Р )
1
1 + Ъ У а
х[ FJ182 ^ + 0.182 у а +
+ 1С 3 Хё у а Хё 008 А + д. 81п А)],
где
= Яе( 3.), д. = ТшАД;)
3 2 3
^ = Ы2 = 3*3,
3 2 3^'
03 = |Ь2\ 2 = Ь23Ь2*3 .
(9)
(10)
Величины (10) определены в работе [1], однако для рассматриваемой задачи имеет смысл использовать и другую форму записи этих величин:
= V3 (а03 + а23) 81п2 Ур + 2а1381п 7Р 008 7Р +
3 (-а03 + а23 ) 00 82 Ур,
дз = 3 81п2 Ур + М381п Ур 008 Ур +
+Т1 3 0082 Ур,
^ = С0 3 81п2 Ур + 2Т13 (-а0 3 + а2 3 ) 81п Ур 008 Ур +
+(С3 - 2С3 ) 00 82 Ур ,
О3 = (С3 + 2с2 3 )81п2 Ур +
3 (а03 + а23 ) 81п Ур 008 Ур + С03 0082 Ур , в которой
а0. = (1 + /2 - 2/ 008 5) 81п 2в3 008 2в3,
3
а1 = (1 + /2) 0082 2в 3 + 2 / 008 5 81п2 2в 3,
а2 3 = (1 - /2 )008 2в3, 3 = 2 / 81п 5 008 2в3 й13 = 2 / 81п 5 81п 2в3,
С3 = (1 + /2 - 2/0085)0082 2в3, С3 = (1 + /2 )(1 + 81п2 2в3) + 2/ 008 5 0082 2в3
(11) (12)
(13)
(14)
(15)
(16) (17)
81(3) А = 0,
8У р
81(3) А = 0,
8у а
(19)
(20)
которые, в принципе, позволяют получить зонные соотношения, определяющие поляризационные углы А и ^ через положения гашения. Однако это не является задачей настоящей работы, зонные соотношения получены в работе [1] из уравнений гашения. Целью данной работы является изучение характера минимума функции IА3(Ур, уа) и установление тех измерительных конфигураций, связанных с положением "быстрой" оси компенсатора, которым отвечает максимальная выраженность минимума этой функции и, следовательно, максимальная точность в измерении положений гашения оптических элементов.
Равенство нулю первых производных, всегда выполняющееся в точке экстремума, само по себе, еще не обеспечивает существования экстремума в соответствующей точке. Это необходимые условия существования экстремума. Достаточным же условием является выполнение неравенства, которое запишем применительно к нашему случаю:
821А3) 821^
8Ур2 8 у а
(3)
(
821
(3) ^
8У р 8у а
> 0.
(21)
Из (21) следует, что вторые производные по каждой переменной в точке экстремума имеют одинаковый знак, плюс для точки минимума и минус для точки максимума. Это означает, что в нашем случае в точке (У(р), у(а3 ) выполняются неравенства
8 21 (3) ^ > 0, 8Ур
2 А 3)
821
8у а
> 0.
(22)
Чем сильнее неравенства (22), тем острее минимум, тем больше его выраженность. Поэтому выбор измерительных конфигураций прибора должен осуществляться путем исследования неравенств (22). Это исследование должно быть проведено для случая произвольного отражающего об-
разца (8), т. е. для случая произвольного набора поляризационных углов ¥ и А. Поэтому в следующем разделе будут получены выражения, определяющие положения гашения поляризатора
(g(p¡)) и анализатора (yj) через углы Y и А .
(j) -
2. ОПРЕДЕЛЕНИЕ ПОЛОЖЕНИИ ГАШЕНИЯ ПОЛЯРИЗАТОРА И АНАЛИЗАТОРА ЧЕРЕЗ ПОЛЯРИЗАЦИОННЫЕ УГЛЫ ¥ И А
Как хорошо известно, положение гашения поляризатора обеспечивает линейную поляризацию отраженной от образца световой волны (на входе анализатора). Из этого и будем исходить. Для этого воспользуемся процедурой, связанной с использованием матриц Джонса. Последовательно действуя матрицами Джонса поляризатора (МР), компенсатора (Мк ) и образца (М5) на амплитудный вектор-столбец Q0 с элементами, образованными р- и 8-составляющими полной амплитуды Е(00) на входе поляризатора, получим вектор-столбец QS после отражения от образца (на входе анализатора):
Qo =
Qs =
f E(0) ^ no p
E (0>
V^ 0s
f E (s) ^
n0 p E (s)
V 0 s 0
= MsMKMpQo = E0 p )
f tg Ye1 А b1
2 j
(23)
где b1. и b2. определяются формулами (5) и (6), а
180°. В этом случае отношение
E
(s)
b
2 j
E (s) П0 p
tg Ye1 Abj j
(25)
является действительной величиной и, следовательно, угол ур, обеспечивающий линейную поляризацию, можно найти из условия равенства нулю мнимой части отношения (25). Учитывая (10), разделим в (25) действительную и мнимую части
7( S )
2 j
e"1А j j
E0SJ tg Ye1 A bj; tg Ybj*.bj
hj cos A + qj sin A + i(qf cos A - hj sin A) = tg Y-Fj
и запишем соответствующее уравнение: (qf cos A- h. sin A) = 0 .
(27)
Преобразуем это уравнение, используя (11), (12) и (15)-(18):
Lj tg2 gp - 2h M0tg gp - N0= 0,
if 0j
0 j
где
L0j = d0 j cos A + (a0 . + a2 .) sin A
0 j 1 2j>
M0j = d1 j cos A - a1 j sin A, N0j = d0j cos A + (a0f - a2f) sin A.
Из (28) находим:
h m 0 j ±
(28)
(29)
tg gP =■
L
(30)
'0 j
E0P ) = Efp cos Yp + E00 'sin Yp . (24)
Линейную поляризацию световой волны после отражения от образца, очевидно, обеспечивает то значение угла gp, при котором разность фаз элементов E(0Sp и E0S-1 вектора QS составляет 0 или
-0 p
1Г1 j
R0 = (sin S cos A - cos S sin 20j sin A)2 + +cos2 20j sin2 A. (31)
Углы gp , определяемые формулой (30), фиксируют те положения поляризатора (положения гашения), которые обеспечивают линейную поляризацию световой волны после отражения от образца и, следовательно, возможность гашения светового пучка на выходе анализатора. Как следует из формулы (30), угол gp зависит от типа ориентации "быстрой" оси компенсатора, определяемого параметром hj, и положения этой оси относительно
плоскости падения (от угла 0. ). Двойной знак в (30) означает, что для каждого типа ориентации "быстрой" оси существуют (при заданном 0 j ) два
значения gp , обеспечивающие указанную линейную поляризацию. Можно также сказать, что этим значениям gp отвечают два вектора линейной поляризации, имеющие разный тип ориентации относительно плоскости падения. Это означает, что отношение (25), являющееся в рассматриваемом случае действительной величиной, меняет знак при переходе от одного значения gp к другому.
Этот вывод непосредственно вытекает из общих свойств измерительных зон прибора. Необходимо конкретизировать двойной знак в (30), связав его с типами ориентации вектора линейной поляризации, а значит, и с типами ориентации анализатора, находящегося в положениях гашения. Имеется
в виду, что вектор линеинои поляризации и анализатор в соответствующем положении гашения имеют противоположные типы ориентации. В случае такоИ конкретизации формула (30) определит положения гашения ур) для каждой из 4 измерительных зон. Покажем, что окончательное выражение для положений гашения ур) имеет следующий вид:
*ТР" = Ц М'0- 2 ^, (у = 1.2.3.4). (32)
L
0 j
(h. cosD + q, sinD).
(33)
Zj =
1
cos2 gp) j
(h, cos D + qf sin D). (34)
где
L, = -d0. sin D + (a0. + a2. )cos D M1. = d1. sin D + a1. cos D,
(36)
N1. = d0. sin D + (-a0. + a2. )cos D.
Величина Z,, очевидно, удовлетворяет неравенствам
Z13 <0, Z24 >0,
откуда, в частности, следует
22 - 21 > 0, 24 - 2Ъ > 0 .
(37)
(38)
2^./) в (32) приводит к изменению символов неравенств (37), а значит, и неравенств (38) на противоположные. Из сказанного ясно также, что для доказательства можно ограничиться рассмотрением упрощенного варианта, когда во всех зонах
в, = в, (j = 1,2,3,4).
(39)
Знак отношения (25) определяется, как следует из (26), величиной
В этом случае в формулах (15) и (16) надо сделать замену (39) и опустить индекс у в определяемых величинах, т. е. перейти к величинам
а0, а1, а2 и ё0, ё1. (40)
Соответственно опускаем индекс у и переходим к величинам
Lo,Mо, No, R и Ц,Mx, N1,
(41)
Поскольку модуль этой величины роли не играет, то ее можно умножить на любую положительную функцию, зависящую, в частности, и от номера зоны. В нашем случае наиболее удобной для использования является величина
которые по-прежнему определяются выражениями (29), (31) и (36), но с учетом формулы (39) и перехода к величинам (40). Теперь для упрощенного варианта (39) можем записать:
(.) ц М0 - 2£ /Ж 18У(р0) = 0 0 ^ , (0 = 1,2,3,4), (42)
Lo
Преобразуем 2у , снабдив, в соответствии с (32), угол ур в выражениях (11) и (12) для И. и индексом у :
2. = ЛА. 182 УР" + М18 Ур'} + , (35)
Zj = л.Ц tg2 gp) + 2M1 tg gp) + hN. (43)
Наконец, используя (42) и (43) и принимая во внимание характер параметров и Ц., находим
22 - 21 = 24 - 23 = 8/^2°(АМ0 + МЛ). (44)
А)
Подставив в (44) выражения для А0, М0 и А, М1, придем к конечному результату:
Z 2 Z1 = Z 4 Z3 =
= 16 f2 sin 8 cos2e
Lo
X [(1 + sin 2в) + f 2(1 - sin 2в)] > 0.
(45)
Неравенства (37) и вытекающие из них неравенства (38) имеют место при любых в. (у = 1,2,3,4), в
том числе и тех, которые в совокупности образуют набор независимых значений. Поэтому для проверки правильности формулы (32) необходимо показать, что она обеспечивает выполнение неравенств (38) при любых в.. Этого достаточно, т. к. изменение знака перед корнем (перед множителем
Правая часть равенства (45) положительна при любом значении угла в , что находится в соответствии с (37) и (38) и подтверждает справедливость формулы (32).
Проведенное доказательство позволяет легко увидеть, что изменение знака (с минуса на плюс) перед 2^./ в (32) переводит, как уже указывалось, символы неравенств (38) в противоположные, т. е. не является разрешенным. Таким образом, формула (32) полностью определяет положения гашения поляризатора для каждой измерительной зоны. Появляющаяся здесь неопределен -
ность для ур) типа ±р не играет принципиальной
роли, она не проявляется ни в зонных соотношениях для А и У, ни в выражении для интенсивности светового пучка на выходе анализатора. Что касается формулы для положений гашения у(аз) анализатора, то она, очевидно, непосредственно связана с проведенным здесь рассмотрением и включает в себя положения гашения ур). Она
легко получается из зонного соотношения для угла У [1], но так же легко ее получить из выражения (25). Учитывая характер относительного расположения вектора линейной поляризации и анализатора в положении гашения, находим
ч(0)
шуУ =
Ь(0) 2 3
■18 У =
>3
18
(46)
Ур = УР ):
(47)
а анализатор при заданном типе ориентации сво-
боден. При выполнении условия (47) формула (9) для интенсивности, с учетом выражений (10), определяющих величины и О^, запишется:
1
ЗУ?, Уа) * ЕР1 1^1,,
1 + 182 У а
у + \Ь2У\ У а + +2£3. ¡8Уа ¡8 У(Л|0) 008 А + д(0) 8т А)],
(48)
где верхний индекс (0) указывает на соответствие условию (47). Затем используем соотношения, определяющие поляризационный угол А (см. [1]):
008 А =
И Ь20)|
8Ш А =
-с1
1Ь(0) Ь20)
13 2 3
(49)
где верхний индекс (0) указывает на то, что соответствующие величины определены для положений гашения ур).
Формулы (32) и (46), определяющие положения гашения ур) и у(а 3) также и в зависимости от поляризационных углов А и У отражающего объекта, позволяют осуществить моделирование процесса прохождения светового пучка через оптическую систему РК^А. Такое моделирование должно включать в себя и исследование выраженности минимума интенсивности светового пучка на выходе системы РК8А, что непосредственно связано с точностью измерения положений гашения оптических элементов, а значит, и поляризационных углов А и У .
3. ВЫБОР ИЗМЕРИТЕЛЬНЫХ КОНФИГУРАЦИЙ ПРИБОРА, ОБЕСПЕЧИВАЮЩИХ МАКСИМАЛЬНУЮ ТОЧНОСТЬ В ОПРЕДЕЛЕНИИ ПОЛОЖЕНИЙ ГАШЕНИЯ
Для выбора оптимальных измерительных конфигураций прибора необходимо исследовать вторые производные от интенсивности 1р \ур, уа) по
каждому из углов ур и уа в точке минимума
(Ур), У а3 ), удовлетворяющие неравенствам (22).
Интенсивность светового пучка определяется формулами (8) и (9). Выражение (9) предпочтительнее, оно легко преобразуется к удобному виду в случае, когда поляризатор зафиксирован в положении гашения
В результате получим следующее выражение для интенсивности в точке (у(3), Уа):
1А'\у р\У а)
7( Р)
№
—^ (IЬ20)|18Уа -К0118^)2. (50)
1 + ¡82 У^ ] ^ 1 3 1 ;
Выражение (50) допускает дальнейшее преобразование, основанное на использовании формулы (46):
^(Г РУ а ) *
Е0Р )12 |я|2
1 +182 У аГ 23
Ь20)| (¡8 У а - 18Уа0))2 =
= Е0Р)|2 №2(1 + ^а 18Уа0))2 |Ь^р х
1 + ¡8 У а х ¡82(Уа -Уа0)) »
|ЕР )|2| Л.1 | Ь2°)|2(У а -У а-))2. (51)
I N 1 + ¡82 Уа I 3 I
Окончательная приближенная запись в формуле (51) подразумевает относительно слабые отклонения анализатора от положения гашения.
Используя (51), запишем вторую производную от интенсивности по Уа в точке минимума
(ур), У а3)):
^213 ^ 2
V ~ т а 00
дУ а
Е0Р Т №2 21 Ь20)Г(1+¡82 У а(3)) =
= \Е,
( р )|
№ 2(Ь
(0)1 23
кя
¡82 У) =
= Е0Р )|2 \я,\2 2(О (0) + е,(0)г82 У),
(52)
где нижний нулевой индекс в обозначении производной указывает на то, что производная определена в точке минимума.
Запишем вторую производную в точке минимума и по параметру Ур , но при этом, очевидно, надо использовать выражение для интенсивности, в котором свободны оба угловых параметра (ур и
уа). Выбрав выражение (9), найдем:
'д 2iAj ) 1
V дГв 2
V 'p 0о
' д2 F, ^
E0P)l2\r\2 х
Vdgр 0
tg2 Y
'д2G 1
i+tg2 y a.)
,( j)
ду '
V lP 0
2„,( j)
tg2 y
i+tg2 y aj)
+2Zj tgY T+gb)
''д2 hJ 1
VVdgP 0o
i + 4 y
cos D
'д2 q 1 V^р 0o
sin D
0
(53)
= (i + f 2)sin2 20. + 2 f cos 8 cos2 20.,
' д2 h} 1
дУ 2 V 'р 0o
'д2 q 1
V^p 0o
= 4 [-a j sin 2gpj) + h a j cos 2gp) ],
= 4[-4j sin 2gpj) - hjd0j cos2gPj)].
(57)
(58)
(59)
Производные (55), (56) и (58), (59) можно легко выразить и через tg 0. . Исследуем вторые производные для некоторых частных случаев, допускающих полное аналитическое рассмотрение.
3.1. Частный случай D = 0
Для этого случая sin D = 0, cos D = i и формула (32), определяющая положения гашения поляризатора, существенно упрощается. С учетом выражений (i6) для d0j и dij она запишется
tg g p) =
h A j - 2Z jf sin 8 h j sin20. - Z.
•(j) _ 'J V
cos20..
(60)
Числитель и знаменатель каждой дроби этого выражения умножим на 0{р и перейдем от 1ёу(а° к 18 ¥, согласно формуле (46). В результате получим:
"0 j -----J
С помощью элементарных преобразований, связанных с переходом от sin 20. и cos 20. к tg 0., находим:
' д2 I(AJ) 1
V дув2
V 'p 00
E0 I2 R2 х
tg2 Y
G (.0) + F(0) tg2 Y
'д2F} ^
Vдур 00
G( 0)
'д 2G.. ^
"дУ"2
V 'р 00
FJ
(0)
+2Z
''д 2h 1
ду
VV 'p 00
cos D-
'д 2q 1 • ^
-- sin D
. ду 2
V 'p 00 0
. (54)
-ZJ + hJ tg o.. = tg(-zJ ^+tg(h°J J t + zhjtgqj " i+tg(zjp4)tg(h)
tg уp) =
p
= 18(-С. 4+). (61)
Из формулы (61) легко определяется явная связь положений гашения у) с углами в.:
Производные, входящие в выражения (53) и (54), определим, используя формулы (ii)—(i4):
уPJj) = -ZJ p+hOj ± p,
(62)
или
'д2FJ ^
Vdj2
V 'р 00
'д 2G 1
= 4[-hj (-00 J + a- j )™2уР
+ (-сз. + с-. )cos2TpJ')],
(j).
ду
V 'p 00
= 4[-h j (00 J + o- j )sin2уp
•( j)
+ (C3 j + c- j )cos2уPJ)],
C3j = -(cij - C0j ) =
(55)
(56)
p
ypk) = ур) - h o. =-Z j-± p.
(63)
Для этого случая рассмотрим сначала вторую производную по уа. Опустив в выражении (52) несущественный общий множитель, запишем эту производную в наиболее простой форме:
'д2 iA J) 1
^ a
(G(0) ■
F,(0)tg2 Y).
(64)
Подставив в (i3) и (i4) уи) из (62), найдем выра-
жения для Fj ) и Gj
(0) .
F = [(1+/2) - Z h, (i - /2)](1+Z h,81п20,), (65)
Gj0) = [(i+/2) - Z h, (i - /2)](i - Z h,sin20, )• (66)
Эти выражения можно получить и более простым способом, используя соответствующие формулы для Fj0> и Gj0) из работы [1], которые имеют другую структуру и включают в себя величину y{j, определенную соотношением (63). Данные выражения отличаются только знаком перед величиной Zjhj sin 29j в простых скобках. С учетом (65) и (66) вторая производная (64) запишется
Га2 iA j) л
ау
2
a 00
[(1 + /2) - Z h, (1 - /2)] X
х[(1 + tg2 Y) - Zh, sin 20, (1 - tg2 Y)].
Га 2 ia л л
аУ
2
a 00
® p, т.е. Cjhj sin20j ®1,
j ® 2[(1 + /2) -Zh,(1 -/2)], G
(0) j
Га 2 i a j л
аУ
2
a 00
В другом же:
Р, т.е. Cjhjsin20j ®-1,
F;U) ® 0, Gj
(0)
' 2[(1 + /2) - Z h, (1 - / 2)];
Га2 iA j) л
ау
2
a 00
2[(1 + /2) - Z h, (1 - /2)].
Из приведенных формул ясно, что для каждой измерительной зоны в одном из этих предельных случаев "быстрая" ось перпендикулярна плоскости падения, а в другом — параллельна. При одном положении "быстрой" оси вторая производная имеет минимальное значение, а при другом — максимальное. При переходе от одного положения к другому производная меняется. Уменьшается она или увеличивается — это зависит от номера измерительной зоны, а также от знака величины
1 - tg2 Y .
(75)
(67)
Исследуем выражение (67) для второй производной от интенсивности по Уа . Величина в первой квадратной скобке в (67) может принимать только два значения — 2 и 2/2 » 2, т. е. она не играет принципиальной роли. Основную роль играет вторая скобка. Для классического варианта, когда "быстрая" ось компенсатора располагается под углом 45° к плоскости падения, а в этом случае в3 = 0 , данная производная обеспечивает достаточную выраженность минимума по Уа :
Если поляризационный угол Y мал (tg2 Y • 1), то для одного из этих положений вторая производная близка к нулю, а для другого достигает значения, в два раза превышающего ее значение для классического варианта. Но надо иметь в виду, и это отмечено в работе [1], что для предельных положений "быстрой" оси измерительная процедура "нулевой" эллипсометрии теряет смысл. Поэтому можно говорить лишь о приближении к предельным положениям "быстрой" оси.
Затем рассмотрим выражение (54) для второй производной от интенсивности по gp . Как и в
случае второй производной по ya, опустим в нем несущественный общий множитель и учтем, что в рассматриваемом частном случае sin Д = 0, cos Д = 1. В результате, получим
• [(1 + /2) - Z h , (1 - /2)](1 + tg2 Y). (68)
г а2i{j) Л dg p2
tg2 Y
G
(0)
-f:
(0)
tg2 Y
При произвольном же расположении "быстрой" оси (в3 Ф 0) надо различать два предельных случая. В одном из них:
г а2 F, ^
2
' Р 00
(69)
(70)
dg ' i
+2Z
Gf>
га2g л
dg
2
Р 00
F,
(0)
Га2 к
аg
2
Р 00
(76)
2[(1 + /2) - Z h, (1 - /2)]tg2 Y . (71)
Здесь нам понадобятся вторые производные в точке минимума от функций , О3 и Н по ур . Найдем их, используя (55), (56) и (58), а также соотношение (62) для у |Р
/ j): р •
(72)
(73)
(74)
Га2 F л
дg 2
V ' р 00
Га 2g, л
аg2
V ' Р 00
Га2h л
аg2
V ' Р 00
= -4Zh,[(1 + /2)sin20, - (1 -/2)], (77)
= 4Z h, [(1 + / 2)sin20, + (1 - /2)], (78)
= 4Z,(1 + /2)cos20,. (79)
x
Преобразуем (76), используя выражения (65), (66), определяющие ^(0) и Gj0), а также производные (77)-(79):
Га/^) ^
¿7 р 2
8[(1 + /2) + СЛ1 (1 - /2)]tg2 ^ (1 + 2У) - С Л, вш 20; (1 - tg2 У)
. (80)
Для классического варианта (0j = 0) производная (80) запишется
Г а21У) ^ ¿7 р2
8[(1 + /2) + С Л, (1 - /2)]
tg2 У
(1 + tg2 У)
, (81)
Г д21к) ^
. д7Р
а для случая (72):
4[(1 + /2) + С Л, (1 - /2)]:
(82)
Г д21у) ^
д7 р2
• 4[(1 + /2) + С Л, (1 - /2)]tg2 У . (83)
Как видим, в первом предельном случае производная не зависит от угла У. При движении в сторону классического варианта, а затем — второго предельного случая зависимость от угла У возникает и усиливается.
Необходимо дать важное разъяснение, относящееся к рассмотренным предельным случаям (69) и (72). Частично оно уже дано и относится к положению "быстрой" оси. Здесь мы приведем более конкретное описание ситуации. Для первого предельного случая формула (62) определяет, если
учесть (69), стремление 7(1) к нулю
7°) ® 0.
(84)
измерительных зон, то
0■ ® ^
к
к
, у^) =И + 0 .
1 4 Гк 4 1 2
(85)
зон, то
0, 4,
j 4
к
У Г = 4 + 01
0.
(86)
где величина в квадратных скобках, как и соответствующая величина в выражениях для второй производной по уа (см. (71) и (74)), отличающаяся знаком перед , принимает те же два близких значения — 2 и 2 /2 » 2 .
Теперь исследуем производную (80) для тех же предельных случаев. Для предельного случая (69):
Затем учтем, что направление пропускания поляризатора, находящегося в положении гашения, в рассматриваемом частном случае (А = 0) параллельно одной из главных осей идеального компенсатора. Это обеспечивает сохранение линейной поляризации светового пучка после прохождения компенсатора и, следовательно, сохранение ее после отражения от образца [3]. В результате на основе формул (84)-(86) можно сказать следующее. В 1-й и 4-й измерительных зонах "быстрая" ось стремится к 8-направлению, перпендикулярному к плоскости падения, а направление пропускания поляризатора, находящегося в положении гашения, должно быть перпендикулярно "быстрой" оси, т. к. только в этом случае оно стремится, согласно формуле (84), к р-направлению, параллельному плоскости падения. Во 2-й и 3-й зонах "быстрая" ось стремится к р-направлению, и, следовательно, направление пропускания поляризатора должно совпадать с "быстрой" осью, стремясь вместе с ней к р-направлению, что опять-таки согласуется с (84). Из всего этого наиболее существенным является то, что в первом предельном случае вектор линейной поляризации после отражения от образца независимо от значения угла У стремится к плоскости падения, а это означает, что
у а1)
к ®—. 2
(87)
Во втором предельном случае формулы (62) и (72) дают следующие результаты для положений
гашения поляризатора (7р ) и положений "быстрой" оси компенсатора (0■):
7
(■) =-С р Ь 1 „ :
1 =1,4 =+1, 01 ^ укт=К4+01 ®o,
1 = 2,3 сл, =-1,
Если С Л, =+1, а это имеет место для 1-й и 4-й
_ к (■) к к
0, , у к) = —+ 0, ® —.
1 4 к 4 1 2
(88)
(89)
(90)
Если же С Л =-1, что имеет место для 2-й и 3-й
Не будем повторять рассуждения, аналогичные предыдущим. Отметим только следующее. Поскольку в формулу (88) входит параметр С1 , то
при анализе соотношений (89)-(90) надо учитывать тип ориентации "быстрой" оси. Кроме того, надо иметь в виду, что тип ориентации вектора линейной поляризации световой волны в данном
частном случае (А = 0) не меняется при отражении от образца [3]. Вообще-то, все это надо было учитывать и в предыдущем случае, но там это не столь критично, т.к. формула (84) имеет очень простой вид. Основная особенность второго предельного случая состоит в том, что вектор линейной поляризации после отражения от образца независимо от значения угла ^ стремится к 8-направлению, а это означает, что
y aj) ® 0.
D» 0,
Y • i.
p
Для классического варианта (y{¿) = —, 0. = 0)
4 J
уPJ) = Z . p+h 0J ± p.
ляющие для произвольных значений 0. величины
F(0), G
'д2 h.
(0)
ду
и
V ' р 00
'д2 iA л i 2
д21(J) 1
'д2 F. >
а также
a0
ду
'д2G 1
V ' р 00
ду
V ' р 00
ду p
, получаются из формул (65)-
0
(9i)
(67) и (77)-(80) путем перехода
Zj ®-Zj.
(94)
Рассмотренные предельные случаи не являются равноправными. Это особенно становится ясным, если рассмотреть малые значения угла ^ (малые значения величины 18 ^). Похожая ситуация возникает при измерении поляризационных углов в окрестности угла Брюстера. Справа от угла Брю-стера
(92)
Таким же способом, исходя из (68) и (81), определяются вторые производные по уа и ур для
классического варианта (в. = 0). Характер этих
производных не меняется. В предельных же случаях положение изменяется. Кроме перехода (94), здесь происходит еще обмен соответствующими выражениями между предельными случаями. Рассмотрим эту ситуацию подробнее. В первом предельном случае ((69)):
p
вторая производная (68) по ya имеет значение, обеспечивающее достаточную выраженность минимума интенсивности по ya. В то же время, вторая производная (8i) по ур в классическом варианте имеет второй порядок малости по Y , поэтому в традиционном эксперименте при определении поляризационного угла D возникают большие затруднения, связанные со слабой выраженностью минимума интенсивности по ур . Если же обратиться к предельным случаям, то из формул (7i), (74) и (82), (83) следует, что для точного определения поляризационных углов для угла Y надо использовать измерительную ситуацию, близкую ко второму, а для угла D - к первому предельному случаю. Это означает, что процесс измерения углов D и Y в рассматриваемом частном случае (D = 0) и близком к нему случае (92) должен быть разделен. В принципе, при не слишком малых углах Y возможен выбор некоторой оптимальной измерительной конфигурации, обеспечивающей достаточное разрешение как по ур , так и по ya.
Обратимся ко второму частному случаю.
3.2. Частный случай D = p
Для этого случая sin D = 0, cos D = -i и формула (32) с помощью тех же преобразований легко определяет положения гашения поляризатора:
'д2 iA J) 1
у p) ® Z,
2[(i+f2)+Z л j (i - f2)],
(95)
(96)
a /0
'д21A j ) 1
ду p
4[(i + f2) - Zh j (i - f)]tg2Y . (97)
Во втором же предельном случае ((72)):
у(j) ® 0, ¡ p ■>
'д2 iA J) 1 2
д 21(J) 1
^ a
дуp2
2[(i + f2) + Zлj(i- f 2)]tg2 Y .
4[(i + f2) - Z л j (i - f2)].
(98)
(99)
(i00)
0
(93)
Для данного частного случая формулы, опреде-
Чтобы проследить в предельных случаях за процессом прохождения светового луча через оптическую систему РК^А, необходимо учитывать соотношения (95) и (98), а также следить за типами ориентации "быстрой" оси компенсатора и вектора линейной поляризации светового луча перед образцом и после отражения от образца. При этом надо иметь в виду, что для данного частного случая (А = р) тип ориентации вектора линейной поляризации в результате отражения от образца меняется на противоположный [3].
Как и в предыдущем частном случае, большой интерес представляют малые значения угла ^, т. к. сходное положение возникает при измерении
поляризационных углов слева от угла Брюстера, где
D»p, Y • i.
(i0i)
' Я - T7 1 'д-G ^
д2 F,
V^p 00
ду
Л
= 0,
p 0
д 2 q ду
=8Zj . (i06)
p 0
Общий вывод здесь такой же. Процесс измерения углов D и Y в рассматриваемом частном случае и близком к нему случае (101) должен быть разделен. Как следует из формул (96), (97) и (99), (100), для точного определения поляризационных углов для угла Y надо использовать измерительную ситуацию, близкую к первому, а для угла D — ко второму предельному случаю. Для угла Y, как и в первом частном случае, можно использовать и классический вариант.
Перейдем к следующему частному случаю.
Вторые производные от интенсивности найдем на основе формул (105) и (106), используя выражения (52) и (54), в которых опущены те же самые несущественные общие множители, что и в предыдущих частных случаях. В этом случае вторая производная по уа определится выражением (64), которое после подстановки (105) даст следующий результат:
'д2 iA■ 1
C)yc
• 2 (i + tg2 Y).
(i07)
a0
p
3.3. Частный случай D = —
Для этого случая sinD = i, cosD = 0 , и формула (32), определяющая положения гашения поляризатора ур), упрощается. Однако, чтобы довести
аналитическое рассмотрение задачи до конца, воспользуемся идеальными значениями параметров компенсатора
8 = 90°, f = i, т. е. cos 8 = 0, sin 8 = i. (i02) Тогда формула (32) запишется: (.) -h ,■ cos20. - Z ■
tg ур) - lj J =
sin 20j hj sm|2(p4 -)] + Z■ cos[2(p4 -O.)]
(i03)
А дальше выражение (103) преобразуется так же, как и формула (60). В результате, пренебрегая несущественной неопределенностью типа ±р, находим:
ур) =-(Z ■ + h j ) f+h jqj,
p
(i04)
ур) = ур) - h 0J =-(Z j + h J
Подставив в (i3), (i4) и (55), (56), (59) выраже-
ние (i04) для угла у , найдем величины F
(0)
G
(0)
'д2 F, >
ду
'c2g 1
V ' р 00
ду
V ' р 00
'д2 q 1
ду : 1р 0
F(0) = G(0) = 2.
j j
(i05)
Если учесть значение угла А (А = "2), то вторая производная от интенсивности по ур запишется в виде
'д2 I[f) 1
ду р
tg2 Y
0
G (0) + F,
д2 F,
(0)
tg2 Y
2
1 P 00
G
(0)
'д 2G 1
ду
+2Z j
2
V • p 00
ду p
Fj
(0)
(0)
'д 2 q 1
ду
2
P 00.
(i08)
откуда с учетом (i05) и (i06) найдем
' д21A j ) 1
ду p
i6
tg2 Y
0
i
tg2 Y
(i09)
Таким образом, вторые производные от интенсивности по уа и ур в данном частном случае не зависят от положения "быстрой" оси компенсатора (от угла в. ). Обращает на себя внимание, что выражения (107) и (109) для этих производных при выполнении условий (102) полностью совпадают с результатами, которые наблюдаются в предыдущих частных случаях для классического варианта (в. = 0). Этот факт можно объяснить следующим
образом. После прохождения компенсатора (перед образцом) световая волна должна быть поляризована по кругу, ибо только в этом случае (при
А = ) после отражения от образца возникает линейная поляризация. Это, очевидно, наблюдается при любом положении "быстрой" оси (при любом в.) Но при круговой поляризации р- и 8- направ-
и
ления эквивалентны для световой волны. То же самое можно сказать и о линейно поляризованной волне на выходе компенсатора в предыдущих частных случаях, если рассматривать классический вариант. Но это означает, что в обоих случаях положение вектора линейной поляризации после отражения от образца будет одинаковым, что и объясняет отмеченное выше совпадение.
Аналогично рассматривается частный случай 3
А = — к . Он ничем не отличается от только что 2
рассмотренного. Результаты описываются теми же формулами (107) и (109). Используя полученные результаты, перейдем к оценке общего случая.
3.4. Переход к общему случаю произвольных значений углов А и У
В многоугловой эллипсометрии особый интерес представляет рассмотрение углов падения, попадающих в окрестность угла Брюстера. В непосредственной окрестности угла Брюстера проявляются особенности поляризационных углов А и У, использование которых имеет большое значение для успешного исследования поверхности. Но именно в этой окрестности значительно снижается точность измерения поляризационного угла А , что обусловлено малыми значениями угла У . Поведение угла А здесь характеризуется размытой ступенькой. Если идти слева направо, переходя через угол Брюстера, то А изменяется от значения А»к до значения А»0(2к), т. е. переходит через точку А = к/2 (3к/2). Как уже отмечалось выше, слева и справа от угла Брюстера можно кардинально повысить точность измерения
положений гашения поляризатора 7р), а значит, и
угла А путем выбора соответствующих измерительных конфигураций прибора. С приближением (слева) к углу Брюстера, т. е. к значению А = к/2(3к/2), ситуация ухудшается и, наконец, при достижении этого значения проявляется, причем при любых значениях угла 01 , неблагоприятный результат, характерный для классического варианта. При последующем отходе (вправо) от угла Брюстера ситуация снова резко улучшается. Таким образом, значительные трудности возникают в сравнительно небольшой окрестности угла Брю-стера. Вне этой окрестности характер ступеньки, определяемый верхним и нижним размытыми участками, проявляется полностью, и именно в этой области может быть обеспечена достаточная точность измерений. Но надо еще раз отметить, что процесс измерения углов А и У в данной области ради достижения максимальной точности в измерении обоих углов может быть разделен. Понятно также, что небольшая окрестность угла Брюстера,
в которой могут быть допущены большие ошибки в измерении угла А , должна быть исключена из процесса измерений. Это не приведет к упущению каких-то существенных дополнительных возможностей, связанных с описанием ступеньки.
В общем случае произвольных значений поляризационных углов А и У используется расчетная процедура, основанная на применении соответствующей математической программы. С помощью формулы (32) для любой пары (А, У) определяются для набора значений угла 0 1 положения гашения поляризатора 7р), по которым рассчитываются вторые производные от интенсивности световой волны по переменным уа и 7р . Эта
расчетная процедура позволяет выбрать те положения "быстрой" оси (те измерительные конфигурации), которым отвечают максимальные значения указанных производных. Если эти положения сильно различаются, то процесс измерения углов А и У можно разделить. В определенных ситуациях можно выбрать некоторую общую оптимальную измерительную конфигурацию, обеспечивающую достаточное разрешение как по 7р , так и по уа . Что касается привязки к конкретному эксперименту, то здесь есть различные возможности. Например, можно использовать предварительные приближенные экспериментальные данные о характере зависимостей углов А и У от угла падения светового луча, которые могут быть получены на приборе с классической измерительной конфигурацией (01 = 0). Если же процесс измерений автоматизирован по "нулевой" схеме, то соответствующая расчетная процедура может быть включена в единый измерительный процесс.
ЗАКЛЮЧЕНИЕ
Полученные в настоящей работе результаты указывают на огромную роль выбора измерительных конфигураций прибора для повышения точности экспериментального определения поляризационных углов. Анализ измерительных конфигураций прибора будет продолжен. Во-первых, необходимо конкретизировать результаты численного моделирования, проиллюстрировав их соответствующими графиками. Во-вторых, не исчерпаны полностью возможности аналитического рассмотрения. В-третьих, необходимо подключить к проблеме повышения точности "нулевую" эллипсо-метрию анизотропных сред. В методе обобщенных измерительных зон отступление от классической измерительной конфигурации является обязательным [5]. Кроме того, этот метод очень чувствителен к точности определения положений гашения.
При этом в каждой простой измерительной зоне процесс определения положений гашения нельзя разделять, относя его к разным измерительным конфигурациям. А это означает, что необходимо искать единые для каждой зоны оптимальные конфигурации. Наконец, проблема повышения точности в "нулевой" эллипсометрии тесно связана с усовершенствованием методов прецизионного определения параметров компенсатора. Все эти вопросы будут рассмотрены в следующей работе.
СПИСОК ЛИТЕРАТУРЫ
1. Семененко А.И., Семененко И.А. // Научное приборостроение. 2007. Т. 17, № 4. С. 42-54.
2. Семененко А.И. // Оптика и спектроскопия. 1978. Т.45. С. 199-201.
3. Ржанов А.В., Свиташев К.К., Семененко А.И.
и др. Основы эллипсометрии. Новосибирск: Наука, 1979. 422 с.
4. Семененко А.И., Семененко И.А. // Научное приборостроение. 2005. Т. 15, № 4. С. 74-82.
5. Семененко А.И., Семененко И.А. // Научное приборостроение. 2007. Т. 17, № 2. С. 20-34.
Институт прикладной физики НАН Украины, г. Сумы (Семененко А.И.)
Институт аналитического приборостроения РАН, Санкт-Петербург (Семененко И.А.)
Материал поступил в редакцию 15.01.2008.
ON THE NEW POTENTIALS OF ELLIPSOMETRY ARISING FROM THE NULL OPTICAL CIRCUIT.
ELLIPSOMETRY OF REAL SURFACE STRUCTURES.
11. PROBLEM OF IMPROVING ACCURACY IN THE NULL ELLIPSOMETRY MEASUREMENTS.
DEFINITION OF THE ROLE OF THE INSTRUMENT METERING
CONFIGURATION
1 2 A. I. Semenenko , I. A. Semenenko
1 Institute of Applied Physics NAS, Ukraine, Sumy 2Institute for Analytical Instrumentation RAS, Saint-Petersburg
The paper considers the instrument metering configurations dependent on the position of the ellipsometer phase compensator fast axis. We have shown that they play a very important role in increasing the accuracy of experimental estimation of polarization angles A and Y. Particular cases (A = 0, p, p/2, 3p/2) permitting complete analytical consideration have been studied. For those cases, relations for second derivatives of the light beam intensity at the instrument output by the polarizer and analyzer angular positions have been constructed. As a result, metering configurations providing maximum derivatives of the mentioned type, which ensure sufficient pronouncement of the intensity minimum and, therefore, the necessary accuracy of experimental estimation of the A and Y angles, have been revealed. A conclusion was made that, in certain situations when angle Y is small, the process of measuring angles A and Y can be subdivided. At the same time, we have noticed that in many cases it is necessary to make efforts to select some general optimal metering configuration ensuring sufficient resolution by both polarization angles. Based on the results obtained in particular cases, the metering situation in the vicinity of the Brewster angle has been analyzed. The paper considers a general case of arbitrary polarization angles A and Y. In addition, the paper discusses questions that need further analysis of the instrument metering configurations. Particularly, a conclusion has been made that it is necessary to thoroughly study the role of metering configurations in null ellipsometry of anisotropic media.