Математическая и программная обработка данных спектрофотометрических измерений Текст научной статьи по специальности «Математика»

УДК 535.8

МАТЕМАТИЧЕСКАЯ И ПРОГРАММНАЯ ОБРАБОТКА ДАННЫХ СПЕКТРОФОТОМЕТРИЧЕСКИХ ИЗМЕРЕНИЙ Е.В. Альтшулер, Э.С. Путилин

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

Введение

При изготовлении интерференционных покрытий в вакууме одной из основных задач является точный контроль толщины слоев в процессе осаждения. В настоящее время широко распространены два метода контроля толщины: резонансный по массе и спек-трофотометрический по R или T. Оба метода - косвенные и имеют свои недостатки: первый - пористость покрытий (контрольная масса), второй - зависимость показателя преломления n от толщины d и технологических факторов, что влияет на значение оптической толщины nd. [1].

Поскольку возрастает интерес к металлодиэлектрическим покрытиям, то при их изготовлении необходимо не только точно контролировать толщину слоя, но также и отслеживать значения n и к, так как они оказывают влияние на фазовый сдвиг при отражении. Приведем формулы расчета R и T для поглощающего слоя, т.е. прямой задачи [2]. Обозначим ф = ф - iф' - фазовая толщина слоя, где

ф = 2nn d/Л, ф' =2пк d/Л , (1)

Л - длина волны падающего излучения. При таком представлении матричные элементы можно записать в следующем виде: m11 = cos ф еЬф', m22 = cos ф еЬф', m'n = si^ shф', m'22 = si^ shф', n si^ Лф' + к cosф shф' n2 + к2 '

к si^ Лф' - n cosф shф' . (2)

m12 =

m

12 _ 2»2 ' п + к 2

т21 = п sinф еЬф' - к cosф shф',

ш'21 = -к мпф chф' - п cosф shф'.

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

я = 4+4, Т = , (3)

Хг2 . ЛТ2 ' Л/-2 I т/2 4 '

где

V = n0mn + m2' - nt (n0m'2 + т22 ) + к1 {n0mn - m'22 ) Z = nom'' - т2' + щ (Пот'2 - m'22 ) + к1 fem^ + т22 ) X = nmii - m2' - nt (nom'2 - т22 ) + к1 ^т^ + m'22 )

Y = nom'i + т21 + ni (nomi2 + m22 ) + к1 (nomi2 - т22 )

(4)

Постановка задачи

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

Поэтому целесообразно проводить измерения Я и Т на двух подложках с разным п1 , что к тому же позволит снизить машинное время, так как п1 входит только в конечные формулы (4), а X - во все элементы матрицы интерференции. Следует заметить, что время расчета - весьма немаловажный фактор, так как слои металлов осаждаются с большой скоростью.

Методы исследований

В процессе осаждения можно измерить величины Т1, Т2, Я1, Я 2 . Поскольку система из четырех уравнений (3) может иметь несколько решений, вводим дополнительные условия, накладывающие определенные рамки на диапазон решений: во-первых, диапазон, в котором ожидается нахождение п и к, а во-вторых, условие, что d возрастает от измерения к измерению, начиная с нуля и далее с некоторым конечным приращением.

Для решения обратной задачи можно выделить 3 варианта.

Вычисление данных методом последовательного поиска. Задаются начальное значение для толщины слоя d, а также граничные условия для параметров п (от 1.8 до 2.2) и к (от 2.7 до 3.3). Задаются допустимые отклонения расчетных параметров Т1, Т2, Я1, Я 2 от реальных. При совпадении всех четырех условий значения п, к, d, для которых это совпадение произошло, заносятся в буфер. Если при том же значении параметра d , но уже других п и к, опять происходит выполнение условий, то текущие значения п и к снова заносятся в буфер. Затем все занесенные в буфер значения п и к усредняются. Алгоритм работает достаточно точно, но его применение нецелесообразно в силу низкой скорости выполнения. Как правило, нужно контролировать значения параметров п, к, d во время осаждения.

Метод поиска решения из сформированного массива данных. Этот метод отличается от предыдущего тем, что решение прямой задачи и поиск обратного решения разделены во времени. Предварительно проводится решение прямой задачи для всего диапазона данных и инициализируется массив. При этом адрес каждого элемента массива равен значению Т1, Т2, Я1, Я 2 , а сам элемент содержит значения показателей п и к и толщину d. После того как массив создан, все значения п, к, d усреднять не имеет смысла в силу огромных вычислительных затрат. Далее, в процессе осаждения из получающихся значений Т1, Т2, Я1, Я 2 формируется адрес массива, считываются его элементы и затем усредняются. Скорость получения результатов при таком способе огромна.

В представленных формулах выделить обратную зависимость

п = 8 1 С^ Т 2 , Я ^ Я 2 )

< к = 8 2 (X, Т 2, Я1, Я 2 ) d = 8 3 (Т , Т 2, Я 1, Я 2 )

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

X - к, й) = 0,

Кг - /(п,к,й)= 0,

<ТХ - / (п, к, й) = 0,

т2 - /4 (п, к, й) = 0

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

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

Остановимся подробнее на численном варианте решения обратной задачи.

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

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

Рассмотрим нелинейную систему уравнений: 'А(( ^к , хп)= °

/2 ((

'2''

,)= 0,

(5)

Ап (х1, Х2, ■■■ , Хп )= 0,

или, в векторной форме,

/ (Х )= 0,

где

/ =

" /1" Х1

/2 — Х 2

, Х =

_ /п _ _ Хп _

(6)

(7)

Для решения системы (5) будем пользоваться методом последовательных приближений. Предположим, известно к-е приближение,

Хч - IX

(к) = (Х (к) Х (к)

Х(к))

= (хп

Х1 , Х 2 -

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

(8)

векторного уравнения (6). Тогд а

2

х = х(к } + А х(к) , (9)

где Ах(к) = (ах(к\ А х(к, Ах(к)) - поправка (погрешность корня). Подставляя выражение (9) в (6), будем иметь

7 (х(к } + А х(к })= 0. (10)

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

7 (х(к) + А х(к ))= 7 (х(к))+ 7'(х(к ))А х(к ) = 0, (11)

или, в развернутом виде,

7(к ) + А х<к) + А х<к))= ),■.., х„( )}+А + К + А = а

дх1

дх„

(12)

7п + А х/к) ,к, х^ к + А х}-к >) = 7п ( ,К, х^к ))+А х1(к + К + А х^к = 0.

Из формул (11) и (12) вытекает, что под производной 7 (х) следует понимать матрицу Якоби системы функций 7\, 72, ■ ■ ■, 7 п относительно переменных х1, х2, ■.., х п т. е.

7 (х)= W(х) =

7 7 71

д х1 дх 2 дхп

д/2 72

д х1 дх 2 дхп

7 7

7

дх дх.,

дх„

или, в краткой записи,

7 (х)= W(х) =

7

дх,

(г, ] = 1,2, ■ , п).

Поэтому формула (11) может быть записана в следующем виде: 7 (х(к))+ W (х(к ))А х(к ) = 0.

7"

Если W (х )= ёй

дх

* 0, то Ах(к) = - W -1 (х(к))7(х(к)) Отсюда

видно, что ме-

тод Ньютона решения системы (5) состоит в построении итерационной последовательности

х (к+1) = х (к)- W -1 (х(к))7(х(к)) к = 0,1,2, К .

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

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

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

Теперь вернемся к нашей обратной задаче. Исходная система уравнений, вообще говоря, не является алгебраической. Модифицируем уравнения (2), введя замену переменных:

cos ср = Zj, ch <р' = Z1, sin ср = Z2, sh ср' = Z 2. (13)

Тогда уравнения (2) примут вид:

т11 = Z1Z1, m22 = Z1Z1, т11 = Z2 Z2, m22 = Z2Z 2,

m = nZ2Z1 + kZ1Z2 = kZ2Z1 - nZ1Z2 (14)

12 _ 2 ,2 ' OT12 " 2 ,2 3

n2 + k2 n2 + к2

m21 = nZ2 Z1 - kZ1Z^, m21 = -kZ2 Z1 - nZ1Z2. а система (4) будет выглядеть следующим образом (в предположении, что k¡ = 0 ):

' ' ' kZ Z '-nZ Z '

V = n0Z1Z1 - kZ2Z1 - nZ1Z2 - nino -—2-" - niZ1Z1 =

n + k

(no - n)Z1Z1 - k( 1 + Jf^ ZZ/-ní 1 - -Jn+kT IZ1Z2

, , , nZ Z '+kZ Z ,

Z = noZ2Z2 - nZ2Z1 + kZ1Z2 + nino 2 T2 niZ2Z2 =

n + k

(n0 - ni )Z2Z2 + k+

-^r 1^1 - Z Z1

n + k I V n + k I

, , , kZ Z ,- nZ Z , ,

X = n0Z1Z1 + kZ2Z1 + nZ1Z2 - nln0—^-+ nlZ1Z1 =

n + k

(n0 + ni )Z1Z1+ k t1 - nn^ )Z ■ Z1,+n(1+nTF hZ 2

F = n0Z2Z2 + nZ2Z1 kZ1Z2 + nin0 2 21 21 2 + niZ2Z2 =

n + k

(0 + ni )Z2Z2 - k(1 - )Z1Z2'+njj + nnf^Z2Z1'

Если ввести обозначения

а = n0 - ni d = n0 + ni

b = k

С nn ^ ^

1 + I e = k

n2 + k2 I

2 , 12

V n +k I V

1-

n0 ni

nn

с = -n

n0 ni

1 —Ffr I / = n

V

n + k

1 + -

n0 ni

n2 + k2

V

то уравнения (15) - (18) примут вид: V = aZ1Z1 - bZ 2 Z1 + cZ1Z 2

Z = aZ 2 Z 2 + bZ1Z 2 + cZ 2 Z1

X = dZ1Z1 + eZ2 Z1 + fZ1Z 2 '

F = dZ2 Z22 - eZ1Z2 + 2 Z1

(15)

(16)

(17)

(18)

n2 + k2 J , (19)

(20)

V2 = a2+ Ь2Z22Z12 + с2Z12Z22 - 2abZ 1Z, Z2 + 2acZ 12Z1Z2 - 2bcZ 1Z^2Z2 Z2 = a2 Z 22 Z 22 + Ь 2 Z12 Z 22 + c2 Z 22 Z12 + 2abZ 2 Z 22 ^ + 2acZ 22 Z 2 7 + 2bcZ 1 ZlZ 2 Z 2 X2 = d2Z12+ е2Z22Z1,2 + 72Z12Z22 + 2deZ 1 Zl2Z2 + 2d/Z2ZlZ2 + 2e/Z 1 ZlZ2Z2 . (21) Y2 = d2Z22Z22 + е2Z12Z22 + 72Z22Zl2 - 2deZ2Z22Z1 + 2d/Z22ZlZ2 - 2е^ 1 ZlZ2Z2 Перепишем формулы (3) следующим образом:

|ЯХ2 +ЯУ2 -V2 -72 =0 \ТХ2 + ТУ2 -4п0щ = 0 .

(22)

Тогда

(хd2 - а2 )(т^2+ 722722)+ (яе2 - Ь2))22+ 72722 )+ (я72 - с2^722 + 722)-

2(Xed + аЬ)((22 - 2 + 2((7d - ас) + 722))2 = 0,

Td 2 (2 + 7 22 7 22 )+ Те 2 ((22 + 7 22 )+ Т 7 2 (2 7 22 + 7 22 )-2Ted((22 - ^^2 + 2Т7d((2 + 722))2 - 4п0щ = 0,

Я = {Я;, Я2 }, Т = {Т;, Т2 }.

Решая эту систему, мы получим значения неизвестных 21,71^2,72. Далее остается решить систему уравнений (13) относительно неизвестных п, к, d:

„ , 2яnd Л , (2лМЛ „ . (2япй Л . , (2яkd А = С08\ —:— | 7, = сп\ -:— I 72 = 81и| -;- | 7 = ЯП\ ■

Л ) 1 \ Л Г2 л Г2 "Л л ). (23)

Так как в системе (23) на 4 уравнения приходится всего 3 неизвестных, то решений будет несколько. Для выбора достоверного ответа предполагается руководствоваться физическими соображениями.

Заключение

Рассмотренные методы решения обратной задачи позволяют определять оптические постоянные и толщину слоя с погрешностью <4% при измерении коэффициента пропускания (Я) и отражения (Т) с погрешностью 0.1%.

Метод выборки из массива позволяет обрабатывать измеряемые данные со скоростью, превышающей быстродействие (более 60 Гц) матрицы ПЗС.

Для повышения быстродействия и точности расчетов ведется работа по созданию алгоритма решения обратной задачи методом Ньютона.

Литература

1. Андреев С.В., Карасев Н.Н., Путилин Э.С., Шакин А.О. Автоматизация фотометрического контроля толщины осаждаемых слоев // Известия вузов. Электроника. -2003. - № 6. - С. 85-90.

2. Борн М., Вольф Э. Основы оптики. - М.: Наука, 1970. - 721 с.

Альтшулер Евгений Владимирович Путилин Эдуард Степанович

ООО «Промкоммуникации», программист, xrocky@rambler.ru

Санкт-Петербургский государственный университет информационных технологий, механики и оптики, доктор технических наук, профессор, заведующий кафедрой оптических технологий, eputilin@yandex.ru