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

Обоснование и практическое приложение методов решения математически некорректной обратной задачи эллипсометрии. 1. О роли угла Брюстера в решении обратной задачи. Общие свойства и выбор комплекса углов падения светового пучка Текст научной статьи по специальности «Физика»

CC BY
45
12
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Научное приборостроение
ВАК
RSCI
Область наук
Ключевые слова
ЭЛЛИПСОМЕТРИЯ / ПОЛЯРИЗАЦИОННЫЕ УГЛЫ / МАТЕМАТИЧЕСКИ НЕКОРРЕКТНАЯ ОБРАТНАЯ ЗАДАЧА / КРИТЕРИЙ / ОПТИМАЛЬНОЕ РЕШЕНИЕ / ЧИСЛЕННЫЙ ЭКСПЕРИМЕНТ / СВЕРХТОНКАЯ ПЛЕНКА / ПОДЛОЖКА / ОПТИЧЕСКИЕ ПОСТОЯННЫЕ / ELLIPSOMETRY / POLARIZATION ANGLES / MATHEMATICALLY INCORRECT INVERSE PROBLEM / CRITERION / OPTIMUM SOLUTION / NUMERICAL EXPERIMENT / SUPER- THIN FILM / GROUND / OPTICAL CONSTANTS

Аннотация научной статьи по физике, автор научной работы — Семененко Альберт Иванович, Семененко И. А.

В работе проведен анализ влияния наборов углов падения светового пучка прибора на характер решения математически некорректной обратной задачи эллипсометрии. Показано, что существование только одного набора углов падения, на котором реализуется правильный предельный переход к идеальной ситуации, и заметный разброс оптимальных значений d opt и n opt параметров пленки по угловым комплексам некоторого множества этих комплексов не является ограничивающим возможности метода фактором. Более того подобные общие свойства угловых комплексов приводят к дополнительным возможностям метода решения математически некорректной обратной задачи. Рассмотрена также роль угла Брюстера в формировании используемых для решения обратной задачи комплексов углов падения.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Семененко Альберт Иванович, Семененко И. А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

FEASIBILITY AND PRACTICAL SUPPLEMENT OF SOLUTION METHODS FOR MATHEMATICALLY INCORRECT INVERSE PROBLEM OF ELLIPSOMETRY. 1. TOWARDS AN UNDERSTANDING OF BREWSTER ANGLE FOR SOLUTION OF INVERSE PROBLEM. GENERAL PROPERTIES AND COMPLEX SELECTION OF THE ANGLES OF INCIDENCE OF LIGHT BEAM

In this work the impact analysis of kits for the angles of incidence of mathematically incorrect inverse problem of ellipsometry. It’s shown that the existence only one kit of the angles of incidence, on which the correct passage to the limit towards the ideal situation is realized, and the noticeable scatter of optimal values d opt and n opt of the film parameters in accordance with the angular complexes of some multitude of these complexes isn’t the factor limiting the method’s possibilities. Besides, the similar general properties of angular complexes lead to the additional possibilities of solution method for the mathematically incorrect inverse problem. The role of Brewster angle in formation used complexes of the angles of incidence for solution of the inverse problem also is considered

Текст научной работы на тему «Обоснование и практическое приложение методов решения математически некорректной обратной задачи эллипсометрии. 1. О роли угла Брюстера в решении обратной задачи. Общие свойства и выбор комплекса углов падения светового пучка»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2012, том 22, № 4, c. 30-37

== ПОВЕРХНОСТЬ =

УДК 535.5.511:531.7 © А. И. Семененко, И. А. Семененко

ОБОСНОВАНИЕ И ПРАКТИЧЕСКОЕ ПРИЛОЖЕНИЕ МЕТОДОВ РЕШЕНИЯ МАТЕМАТИЧЕСКИ НЕКОРРЕКТНОЙ ОБРАТНОЙ

ЗАДАЧИ ЭЛЛИПСОМЕТРИИ.

1. О РОЛИ УГЛА БРЮСТЕРА В РЕШЕНИИ ОБРАТНОЙ ЗАДАЧИ.

ОБЩИЕ СВОЙСТВА И ВЫБОР КОМПЛЕКСА УГЛОВ ПАДЕНИЯ

СВЕТОВОГО ПУЧКА

В работе проведен анализ влияния наборов углов падения светового пучка прибора на характер решения математически некорректной обратной задачи эллипсометрии. Показано, что существование только одного набора углов падения, на котором реализуется правильный предельный переход к идеальной ситуации, и заметный разброс оптимальных значений dopt и пор{ параметров пленки по угловым комплексам некоторого множества этих комплексов не является ограничивающим возможности метода фактором. Более того, подобные общие свойства угловых комплексов приводят к дополнительным возможностям метода решения математически некорректной обратной задачи. Рассмотрена также роль угла Брюстера в формировании используемых для решения обратной задачи комплексов углов падения.

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

ВВЕДЕНИЕ. ПОСТАНОВКА ЗАДАЧИ

Исследование поверхностной структуры твердых тел и жидкостей относится к числу сложных и практически важных научных задач. Такая структура проявляется в самых различных формах. Это — сверхтонкие пленки естественного окисла на полупроводниках и металлах. Это — нарушенные слои на поверхности твердых тел, возникающие в результате той или иной обработки. Кроме того, это также — естественная структура поверхности жидкостей, связанная с их общими свойствами и изменяющаяся при добавках различных веществ. Во всех этих случаях обратная задача эллипсометрии является математически некорректной. Математическая некорректность проявляется в том, что даже относительно слабые экспериментальные ошибки заметно искажают истинные значения параметров, характеризующих поверхностную структуру. Очень часто эти искажения приобретают абсурдный характер. Однако надо иметь в виду, что такое проявление математической некорректности непосредственно связано с классическим подходом к решению обратной задачи. В этом случае в качестве решения выбираются параметры, определяющие точку абсолютного минимума функционала обратной задачи. При наличии экспериментальных ошибок причи-

ны, приводящие к математической некорректности, могут быть самыми разными. Для естественных окислов выраженная некорректность обратной задачи обусловлена их сверхмалой (порядка нескольких нанометров) толщиной. Для нарушенных слоев толщина не играет особой роли, и основной причиной заметного проявления некорректности является слабое отличие оптических параметров таких слоев от их объемных значений. Нарушенные слои имеют место всегда, и это обстоятельство надо учитывать при исследовании сверхтонких пленок, имеющих другую природу. Но и сам нарушенный слой, как правило, имеет более выраженную по оптическим параметрам верхнюю часть с относительно небольшой, в том числе и сверхмалой, толщиной. Подобную структуру, по-видимому, имеет и поверхность жидкостей. Таким образом, задача по исследованию поверхностной структуры твердых тел и жидкостей является довольно сложной. Она не сводится только к определению параметров некоторой сверхтонкой поверхностной пленки при заданных параметрах подложки.

В работе [1] предложен новый подход к решению математически некорректной обратной задачи эллипсометрии для однослойной системы. Этот подход основан на использовании измерений поляризационных углов, отвечающих набору углов

падения светового пучка. В работе [1] приведено общее описание правил отбора оптимального решения. Использование данных правил позволяет с высокой точностью найти решение обратной задачи относительно толщины и показателя преломления пленки при известных значениях оптических постоянных подложки. Это убедительно продемонстрировал численный эксперимент, основанный на моделировании экспериментальных ошибок. Однако для реальных исследуемых объектов оптические постоянные подложки не могут быть заданы точно. Во-первых, это обусловлено сложной структурой поверхности любого материала, а значит, и подложки. Во-вторых, между пленкой и подложкой всегда возникает некоторый переходный слой. Данные факторы являются причиной того, что подложка характеризуется уже некоторыми эффективными значениями оптических параметров. Эти эффективные значения могут заметно отличаться от оптических параметров, определяющих объемные свойства подложки. В работе [1] показано, что математическая некорректность обратной задачи значительно усиливается из-за неточного задания оптических параметров подложки. Это наблюдается даже тогда, когда неточность в задании показателя преломления подложки проявляется в четвертом или даже пятом знаке после запятой. Поэтому речь может идти только об одновременном определении всех параметров однослойной системы — как поверхностной пленки, так и подложки.

В последующих работах для решения математически некорректной обратной задачи эллипсо-метрии в ее общей постановке предложен метод последовательных приближений. Основу этого метода составляют правила отбора оптимального решения для параметров сверхтонкой пленки или нарушенного слоя при фиксированных на каждом этапе решения оптических параметрах подложки. При этом рассматривается процесс по сближению двух точек. Одна из этих точек определена параметрами пленки в точке абсолютного минимума функционала

выполнялось условие

п).

(1)

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

о ор1

).

(2)

значения параметров «°°), к"°°) подложки. Выбор

(0)

„(°)

-(0)

значений «° , к° носит относительный характер. Здесь важно лишь, чтобы на первом шаге при любом типе экспериментальных ошибок обязательно

«тш < ПоРг,

орг '

(3)

Условие (3) достигается легко. Это связано с тем, что точка (1) очень подвижна даже при слабом изменении показателя преломления подложки. В то же время точка оптимальных значений (2) обладает ограниченной подвижностью. Это обеспечивает переход точки минимума через точку оптимальных значений с последующим расхождением данных точек

«тш > Пор1

< d

орг '

(4)

Очевидно, обеим точкам соответствуют одни и те же фиксированные для каждого шага метода

Основную роль в перемещении точки минимума (1) играет показатель преломления подложки «°, т. е. его изменение вдоль некоторой возрастающей последовательности значений. Вследствие этого переход точки минимума (1) через точку оптимальных значений (2) наблюдается при любой величине коэффициента поглощения подложки к°. Задача состоит в том, чтобы из множества переходов точки минимума через точку оптимальных значений выбрать тот переход, которому соответствует оптимальное решение, наиболее приближенное к точному решению обратной задачи. И поскольку каждому такому переходу соответствует определенное значение к°, то правильный выбор перехода означает и выбор величины к° как одного из параметров оптимального решения обратной задачи. Данная процедура подробно описана в работе [2]. При этом особое внимание уделено установлению признака, по которому определяется оптимальное значение коэффициента поглощения подложки.

Метод последовательных приближений отработан на численном эксперименте. Данный эксперимент включает в себя моделирование экспериментальных ошибок в задании поляризационных углов ^ и А для набора углов падения светового пучка. Использование набора углов падения светового пучка на образец является необходимым условием для успешного применения метода. В то же время метод демонстрирует надежность и хорошие возможности при исследовании также и реальных отражающих систем типа "полупроводниковая подложка—прозрачная сверхтонкая пленка". Он успешно опробован при исследовании кремния и арсенида галлия с естественными сверхтонкими окисными пленками. При этом речь идет об определении всех 4 параметров таких систем.

Следует отметить, что для успешного использования предложенного метода необходима высокая точность вычислительного процесса. В связи с этим существенно модифицирован хорошо из-

вестный комплексный метод Бокса [3], положенный в основу решения обратной задачи.

Еще раз отметим, что новый подход к решению математически некорректной обратной задачи эл-липсометрии основан на использовании измерений поляризационных углов, отвечающих некоторому набору углов падения светового пучка на образец. Это является необходимым условием для успешного применения метода. В то же время, как показывает анализ, выбор комплекса углов падения светового пучка не является простой задачей. И дело не только в количестве углов в таком наборе, что само по себе также играет большую роль. Есть еще один крайне важный момент, связанный с предельным переходом к идеальному случаю, когда отсутствуют экспериментальные ошибки в используемых поляризационных углах У и А, а отражающая система является строго однослойной. Такой предельный переход при использовании методов решения математически некорректной задачи реализуется не для любого набора углов падения светового пучка. Без учета данного обстоятельства возможны дополнительные ошибки в определении параметров отражающей системы. И когда выше указывается на успешное применение метода при исследовании реальных образцов, то прежде всего имеется в виду подтверждение основных принципов, положенных в основу нового подхода к решению математически некорректной обратной задачи. Возможные дополнительные ошибки, связанные с выбором комплекса углов падения, не искажают и тем более не отменяют эти принципы.

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

ОБЩИЕ СВОЙСТВА КОМПЛЕКСА УГЛОВ ПАДЕНИЯ СВЕТОВОГО ПУЧКА.

О ВЛИЯНИИ УГЛА БРЮСТЕРА И ЕГО ОКРЕСТНОСТИ

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

Для пленки принимаем значения параметров

d = 2.50 нм, п = 1.50, (5)

а для подложки

п0 = 3.865, к0 = 0.023. (6)

Длина световой волны определяется значением Л = 632.8 нм. Что касается значений углов падения р0, то они устанавливаются по результатам исследования свойств некоторого множества угловых комплексов.

Как и ранее, в качестве функционала £0 обратной задачи будем рассматривать среднеквадратичную невязку следующего вида:

1 Ы

^ = ™£[(АИ - А.И)2 + (У,М -У.М)2]. (7)

2М 1=1

Суммирование в выражении (7) производится по набору углов падения ( 0 светового пучка, т. е. величины А, [1], У, [1] и А. [1], У. [1] из данного выражения — это соответственно теоретические и экспериментальные значения поляризационных углов, отвечающие 1-му углу падения. Набор углов падения должен удовлетворять определенным требованиям. В предыдущих работах предполагалось, что данному набору должно соответствовать тщательное экспериментальное описание тех участков кривых У(р0) и А(р0), которым соответствуют достаточно широкий интервал вокруг острия для угла У и вся размытая ступенька для угла А . Это условие связывалось с тем, что указанные участки наиболее чувствительны к изменению параметров поверхностных структур. Вроде бы все правильно, и угол Брюстера, как и его ближайшая окрестность, в той или иной мере должны быть представлены в используемом для решения обратной задачи комплексе углов падения. Однако ситуация гораздо сложнее, чем это кажется на первый взгляд. Именно существенное изменение (на наборе углов падения) чувствительности поляризационных углов к изменению параметров пленки является причиной того, что значение параметра-критерия для всех точек траектории, ведущей (при фиксированных оптических параметрах подложки) к точке абсолютного минимума функционала £0 (см. [1]), существенно возрастает. При этом правила отбора, основанные на использовании свойств параметра-критерия [1], дают оптимальное решение, которое может заметно отличаться от точного решения даже в идеальном случае. Причем, это проявляется тем сильнее, чем ближе подходит угловой комплекс к углу Брюстера и его ближайшей окрестности. Это и означает нарушение предельного перехода к идеальной ситуации. Однако такое нарушение не является общим правилом. Всегда существует такой угловой комплекс, который восстанавливает указанный предельный переход. Ниже мы продемонстрируем это

Зависимость оптимальных значений d и n

opt rlopt

параметров пленки от положения перемещающегося углового комплекса

а (град) (град) <Р0[ Nf ] (град) d t opt (нм) nopt

- 1.500 46.500 72.500 2.440 1.531

- 0.500 47.500 73.500 2.457 1.522

0.000 48.000 74.000 2.475 1.512

0.600 48.600 74.600 2.500 1.500

1.000 49.000 75.000 2.511 1.495

1.425 49.425 75.425 2.523 1.490

1.494 49.494 VBr 2.507 1.496

1.700 49.700 75.700 2.560 1.473

1.850 49.850 75.850 2.562 1.472

2.000 50.000 76.000 2.561 1.473

2.500 50.500 76.500 2.550 1.478

3.000 51.000 77.000 2.544 1.480

3.500 51.500 77.500 2.562 1.472

4.000 52.000 78.000 2.575 1.466

4.500 52.500 78.500 2.576 1.465

5.000 53.000 79.000 2.565 1.471

5.500 53.500 79.500 2.584 1.463

6.000 54.000 80.000 2.592 1.460

на ряде примеров, а затем опишем способ выбора такого углового комплекса.

Основные свойства угловых комплексов установим, принимая для оптических параметров подложки их точные значения (6). Обратимся сначала к угловому комплексу, для которого углы падения определяются выражением

Р°И = Р°ъ + 0' -Р, 1 = 1,...,Ы/ , (8)

где ср°ъ — начальный угол комплекса; h — шаг, с которым изменяются углы; Ы/ — число углов

в комплексе.

Введем также угловой параметр а, определяющий движение комплекса (8) в целом к углу Брюстера или же от него. С учетом этого параметра выражение для углов перемещающегося комплекса запишется

Р°И = р°ъ + (1 -+ а, 1 = 1,...,Ы/ . (9)

При монотонном увеличении параметра а от некоторого заданного значения рассматриваемый угловой комплекс, очевидно, перемещается в сторону угла Брюстера. Определим величины р°ъ, h и Ы, значениями

р°ъ = 48.°°°, h = 2.°°°, = 14. (1°)

При этих значениях конечный угол комплекса (8) запишется:

р° / = 74.°°°, (11)

Что касается начального и конечного углов перемещающегося комплекса (9), то они, очевидно, определятся соотношениями

Р°[1] = Р°ъ + а, р°[ Ы/ ] = р°/ + а. (12)

В целях удобства для углового комплекса (9) введем обозначение

F (а), или F (р°[ Ы/ ]). (13)

Укажем также значение угла Брюстера для подложки (6)

рВг = 75.494° . (14)

Меняя значения параметра а , рассмотрим различные положения углового комплекса (9) и определим отвечающие им оптимальные значения dopt

и порг. При этом необходимо учитывать, что в

идеальном случае при точно заданных параметрах подложки точка абсолютного минимума функционала определяет значения dmln и «тт, совпадающие, причем для любого углового комплекса, с точными значениями (5) параметров пленки:

dmln = 2.5° нм, «т1п = 1.5°. (15)

Соотношения (15) уже дают точные значения параметров пленки, однако интерес представляет предельный переход, реализуемый с помощью метода решения математически некорректной задачи, т. е. путем нахождения оптимального решения и сравнения его с точным решением. Это связано с интересами решения обратной задачи для реального случая. В следующем разделе при рассмотрении предельного перехода относительно всех 4 параметров мы учтем соотношения (15), имеющие место только для идеального случая и только при точно заданных параметрах подложки.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Оптимальные значения dopt и пор(, отвечающие

различным значениям параметра а , т. е. различным положениям углового комплекса (9), приведены в таблице. Величины Ро[1] и Ро[Ы/], приведенные в этой таблице, представляют значения начального и конечного углов перемещающегося комплекса (9) (см. (12)). Для наглядности целесообразно следить за зависимостью величин dopt

и « от конечного угла р°[Ы/ ], который легко

1.53

1.50-

1.46-

80 0 РсЕ ^ ]

Рис. 1. Схематическое изображение зависимости пор1 (р0[]) от положения конечного угла перемещающегося углового комплекса

2.50'

2.44'

1 _1_ 1 1 1 1 1 1 1 1 _1_1

72.5

Рвг

80.0

РсЕ К; ]

Рис. 2. Схематическое изображение зависимости dopt (р0[ ]) от положения конечного

угла перемещающегося углового комплекса

сравнивать с углом Брюстера. Из таблицы видно, что величины йор1 и пор1 при изменении конечного угла в пределах интервала

р0[ ] е [72.50; 80.00] (16)

могут значительно отличаться от точных значений параметров d и п, совпадающих со значениями dmin и птш этих параметров в точке абсолютного

минимума функционала S0. С приближением конечного угла р0[ ] к углу Брюстера зависимость величин dopt и по0 приобретает характер осцил-

ляций, что наглядно представлено на рис. 1 и 2. В этом случае оптимальные значения dopt и по0

особенно сильно отличаются от точных значений параметров пленки. К тому же происходит это на относительно небольшом интервале углов падения светового пучка, включающем в себя и угол Брю-стера. А это не выгодно в любом случае, в частности по причине того, что на угле Брюстера, как и в его окрестности, точность измерения поляризационных углов У и А (особенно угла А) резко снижается. В результате проявляется только один угловой комплекс F (р0[ ]), на котором реализуется правильный предельный переход к идеальному случаю

d , ^ d

орt т

порг ^ пт:п.

Это происходит на угловом комплексе с конечным углом

Данный угол достаточно удален от угла Брю-стера, чтобы не вызывать затруднений в экспериментальных измерениях.

В реальной ситуации нецелесообразно определять комплекс углов падения, на котором реализуется предельный переход (17), рассматривая некоторое множество положений перемещающегося углового комплекса (9). Это связано с очень большим объемом экспериментальной работы по измерению поляризационных углов. Есть другой, гораздо более простой способ. Суть его следующая. Перемещается только конечный угол, а остальные углы фиксируются. Для примера рассмотрим конкретный угловой набор

( + (1 -1)к, 1 = 1,...,(К7 -2); Ро[1 ] = | (N -2] + к,, 1 = (-1); (19)

(0; + 1 = ,

где величины р0ь, р0;, к, N и а имеют тот же смысл, что и для предыдущего варианта, а параметр к — это шаг от угла р0[- 2] к углу

Р0[N -1]. Определим эти величины значениями

р0Ь = 50.50°, к = 2.00°, 1.00°, = 14, р0/ = 74.00°.

(20)

(17)

Р0[N] = 74.60° .

(18)

Параметр а будем определять, начиная от нуля. Здесь мы не будем изучать угловой комплекс (19), рассматривая различные положения конечного угла р0[] относительно угла Брюстера. Отметим только, что и в этом случае ситуация усложняется с приближением конечного угла к углу

Брюстера. Ограничимся рассмотрением углового комплекса (19) при следующих значениях параметра а :

а = 0.00; 0.10; 0.20; 0.30; 0.35; ( )

0.40; 0.50; 0.60; 0.70. ( )

Меняя параметр а согласно распределению (21), убеждаемся, что и в случае углового комплекса (19) существует только один набор углов падения, отвечающий значению а = 0.35° (значению конечного угла р0[Nf ] = 74.35°), на котором

реализуется правильный предельный переход (17) к идеальной ситуации. При остальных значениях величины а оптимальные значения dopt и nopt

параметров пленки заметно отличаются от точных значений этих параметров. Это особенно проявляется на крайних значениях величины а из распределения (21). Так, при нулевом значении а

dopt = 2.480, nopt = 1.510, (22)

а при максимальном значении а из (21)

dopt = 2.520, nopt = 1.491. (23)

Существование только одного набора углов падения, на котором реализуется правильный предельный переход к идеальной ситуации, и заметный разброс оптимальных значений dop t и nopt по

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

О ВЫЯВЛЕНИИ И ИСПОЛЬЗОВАНИИ ОСНОВНОГО КОМПЛЕКСА УГЛОВ ПАДЕНИЯ

СВЕТОВОГО ПУЧКА ПРИ ИССЛЕДОВАНИИ РЕАЛЬНЫХ ОТРАЖАЮЩИХ ОБЪЕКТОВ

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

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

Будем исходить из того, что в идеальной ситуации неизвестны все параметры однослойной системы. Рассмотрим некоторое множество угловых комплексов F(a) (см. (13)), придавая параметру а различные значения. Для всех углов падения каждого комплекса измерены, причем без ошибок, поляризационные углы ^ и А. Будем исходить из того, что в идеальном случае при точно заданных параметрах подложки точка абсолютного минимума функционала определяет значения dmln и пт1п, совпадающие, причем для

любого углового комплекса рассматриваемого множества таких комплексов, с точными значениями параметров пленки. Такое свойство позволяет сформулировать правило, с помощью которого сначала находятся оптические параметры подложки, а затем и основной угловой комплекс. Указанное свойство имеет место только для идеального случая и только при точно заданных параметрах подложки. А это означает, что разброс величин dmln и пт1п по угловым комплексам, отличаясь

при произвольных параметрах подложки от нуля, стремится к нулю только в том случае, если параметры подложки стремятся к своим точным значениям (по)1ги и (Ко)1ги

П0 ^ (П0 ^¡ги , К0 ^ (К0)1ги . (24)

Поскольку задача решается в условиях, когда углы ^ и А измерены, т. е. считаются известными, то нет никакой необходимости знать параметры пленки d и п. В этом случае обратная задача решается относительно оптических параметров подложки п0 и к0. Решением будут те значения этих параметров, которые обеспечивают нулевой разброс величин dmln и пт1п по угловым комплексам множества. Данные значения оптических параметров подложки в рассматриваемом идеальном случае совпадут с их точными значениями. После этого способом, указанным в предыдущем разделе, из множества угловых комплексов выделяется основной комплекс. При этом определяются и параметры пленки. Для этого используются не величины dmln и пт1п, что возможно в идеальной ситуации, а метод решения математически некорректной задачи. Это связано с интересами перехода к реальным отражающим объектам.

Теперь можно сформулировать и способ выявления основного углового комплекса для случая реальных отражающих объектов. Он представляет собой обобщение описанного выше способа, соответствующего идеальной ситуации. Прежде всего решается обратная задача относительно оптических параметров подложки п0 и к0. Значения этих параметров, обеспечивающие минимизацию разброса величин dmln и пШ1П по угловым комплексам множества, представляют собой решение данной обратной задачи. При этом минимизация в общем случае достигается при отличном от нуля разбросе величин dmln и пш1п . По значениям оптических параметров подложки, обеспечивающим указанную минимизацию, для каждого углового комплекса находятся оптимальные значения d

opt

и n

opt

параметров пленки. Основным будет тот

угловой комплекс, для которого выполняются соотношения

d d

opt m

opt

(25)

Если соотношения (25) не выполняются ни для одного углового комплекса множества, то данное множество необходимо расширить. Найденный описанным способом основной угловой комплекс одновременно определяет и параметры d и п , совпадающие с указанными в (25) величинами.

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

СПИСОК ЛИТЕРАТУРЫ

1. Семененко А.И., Семененко И.А. // Научное приборостроение. 2010. Т. 20, № 4. С. 132-142.

2. Семененко А.И., Семененко И.А. // Научное приборостроение. 2012. Т. 22, № 1. С. 44-51.

3. Box M.J. A new method of constrained optimization and a comparison with other methods // Comp. Journ. 1965. V. 8. P. 42-51.

Украина, г. Сумы (Семененко А.И.)

Первый Московский государственный медицинский университет им. И.М. Сеченова (Семененко И.А.)

Контакты: Семененко Альберт Иванович, sem199@mall.ru

n

mill

ЗАКЛЮЧЕНИЕ Материал поступил в редакцию 2.07.2012

Проведенный в настоящей работе анализ ком-

FEASIBILITY AND PRACTICAL SUPPLEMENT OF SOLUTION METHODS FOR MATHEMATICALLY INCORRECT INVERSE

PROBLEM OF ELLIPSOMETRY. 1. TOWARDS AN UNDERSTANDING OF BREWSTER ANGLE FOR SOLUTION OF INVERSE PROBLEM. GENERAL PROPERTIES AND COMPLEX SELECTION OF THE ANGLES

OF INCIDENCE OF LIGHT BEAM

1 2 A. I. Semenenko , I. A. Semenenko

lSumy, Ukraine

2I.M. Sechenov First Moskow State Medical University

In this work the impact analysis of kits for the angles of incidence of mathematically incorrect inverse problem of ellipsometry. It's shown that the existence only one kit of the angles of incidence, on which the correct passage to the limit towards the ideal situation is realized, and the noticeable scatter of optimal values d opt and n opt of the film parameters in accordance with the angular complexes of some multitude of these complexes isn't the factor limiting the method's possibilities. Besides, the similar general properties of angu-

lar complexes lead to the additional possibilities of solution method for the mathematically incorrect inverse problem. The role of Brewster angle in formation used complexes of the angles of incidence for solution of the inverse problem also is considered.

Keywords: ellipsometry, polarization angles, mathematically incorrect inverse problem, criterion, optimum solution, numerical experiment, super- thin film, ground, optical constants

i Надоели баннеры? Вы всегда можете отключить рекламу.