CC BY
110
НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ май-июнь 2015 Том 15 № 3 ISSN 2226-1494 http://ntv.ifmo.ru/
SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS
May-June 2015
Vol. 15 No 3 ISSN 2226-1494
http://ntv.ifmo.ru/en
УДК 535.51:666.192
ПРИМЕНЕНИЕ МЕТОДА ЭЛЛИПСОМЕТРИИ В ОПТИКЕ НЕ ОДНОРОДНЫХ СРЕД А.Н. Горляк", И.А. Храмцовский\ В.М. Солонухаb
а Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В.И. Ульянова (Ленина),
Санкт-Петербург, 197376, Российская Федерация
ь Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация
Адрес для переписки: [email protected]
Информация о статье
Поступила в редакцию 20.11.14, принята к печати 10.02.15
^:10.17586/2226-1494-2015-15-3-378-386
Язык статьи - русский
Ссылка для цитирования: Горляк А.Н., Храмцовский И.А., Солонуха В.М. Применение метода эллипсометрии в оптике неоднородных сред // Научно-технический вестник информационных технологий, механики и оптики. 2015. Т. 15. № 3. С. 378-386.
Аннотация
На основе обобщения теории отражения поляризованного света в приближении Борна от границы раздела неоднородных сред усовершенствованы физико-математические методы моделирования и измерений поляризационно-оптических характеристик отражающей системы «неоднородный поверхностный слой - неоднородная среда». Дано теоретическое и экспериментальное обоснование одновременного применения методов многоугловой и сканирующей эллипсометрии для элементов градиентной оптики при наличии неоднородного слоя на поверхности оптической детали. Измерение основных эллипсометрических параметров - азимута линейной восстановленной поляризации ¥ и разности фаз Д между взаимно ортогональными p- и ^-компонентами отраженного поляризованного светового пучка -осуществлялось на двухканальной установке, где в первом канале используется схема измерений поляризационных параметров ¥ и Д по методу компенсационной (нулевой) эллипсометрии; во втором канале - по фотометрической схеме измерений параметров ¥ и Д с применением метода Фурье-детектирования выходного сигнала. Такое разделение каналов позволяет измерять поляризационные параметры отраженного светового пучка раздельно и тем самым повысить точность эллипсометрических измерений по отношению к исследуемой физической величине объекта измерений. На основе измерений поляризационных параметров ¥ и Д при угле падения светового пучка ф=45° в рамках метода сканирующей эллипсометрии и модели геометрически плоской границы раздела однородных сред были рассчитаны значения показателя преломления п(г) и показателя поглощения Мт) по «аксиальному значению» градана. Полученные результаты указали на некорректность модели геометрически плоской границы раздела сред для поляризационных измерений распределения показателей преломления п(г) и поглощения к(г) в элементах градиентной оптотехни-ки. Для отражающей системы «неоднородный поверхностный слой - неоднородная среда» изменение показателя преломления п(г) по «аксиальному сечению» элемента г пропорционально разности фаз Д(т) между двумя взаимно ортогональными компонентами отраженного светового пучка. При этом азимут линейной восстановленной поляризации ¥(т) по сечению градиентного оптического элемента изменяется незначительно.
На основе уравнения эллипсометрии, полученного в рамках теории отражения поляризованного света Борна, приведены основные соотношения для расчета показателя преломления пс и толщины 4с поверхностного слоя в каждой локальной точке г градана. Полученные результаты экспериментальных исследований позволяют дать количественную оценку качества обработки поверхности элементов оптотехники. Разработанные эллипсометрические методики рекомендованы в технологическом производстве оптических узлов с элементами градиентной и волоконной оптики. Ключевые слова
эллипсометрия, градиентный элемент, поверхностный слой, показатель преломления, толщина слоя.
ELLIPSOMETRY METHOD APPLICATION IN OPTICS OF INHOMOGENEOUS
MEDIA
A.N. Gorlyaka, I.A. Khramtsovkyb, V.M. Solonukhab
a Saint Petersburg State Electrotechnical University, Saint Petersburg, 197376, Russian Federation b ITMO University, Saint Petersburg, 197101, Russian Federation Corresponding author: [email protected] Article info
Received 20.11.14, accepted 10.02.15 doi: 10.17586/2226-1494-2015-15-3-378-386 Article in Russian
For citation: Gorlyak A.N., Khramtsovky I.A., Solonukha V.M.. Ellipsometry method application in optics of inhomogeneous media.
Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2015, vol.15, no. 3, pp. 378-386.
Abstract
Physico-mathematical methods of modeling and measurements of polarization optical characteristics of "inhomogeneous surface layer - inhomogeneous medium" reflective system have been improved based on generalization of the polarized light reflection theory in the Born approximation from inhomogeneous media interface. Theoretical and experimental justification for simultaneous application of multi-angular and scanning ellipsometry methods is given for gradient optics elements in inhomogeneous layer presence on the surface of the optical detail.
The measurement of the main ellipsometric parameters - linear reduced polarization azimuth ¥ and phase difference Д between mutually orthogonal p- and s- components of a reflected polarized light beam - was carried out on a dual-channel installation, in which the scheme of measurements of polarization parameters ¥ and Д by balanced (null) ellipsometry method is used in the first channel; the photometric scheme of measurements of parameters ¥ and Д by Fourier method of output detection is used in the second channel. Such separation of channels gives the possibility for measuring of polarization parameters of reflected light beam separately and, thereby, improving the accuracy of ellipsometric measurements in relation to observable physical quantity of measurements object.
The values of refractive index n(r) and absorption index k(r) were calculated by axial section of selfoc on the basis of measurements of polarization parameters ¥ and Д with angle of incidence of the light beam ф=45о within the limits of scanning ellipsometry method and a model of geometrically flat interface of homogeneous media. Empirical data indicated the model incorrectness of the geometrically flat interface for polarization measurements of refractive index distribution n(r) and absorption index k(r) in gradient optical engineering elements. For the reflective system "inhomogeneous surface layer -inhomogeneous medium" the refractive index n(r) change by axial section of an element r is proportional to phase difference Д(г) between mutually orthogonal components of a reflected light beam. At the same time the linear reduced polarization azimuth ¥(r) in the section gradient optical element varies slightly.
The main proportions for calculation of refractive index nc and thickness dc of superficial layer in each local point r of selfoc are adduced on the basis of the ellipsometry equation, obtained within the framework of the polarized light reflection Born theory. Empirical data of experimental researches make it possible to give quantitative assessment of surface processing quality for optical engineering elements. Elaborated ellipsometric techniques are recommended in technological production of optical base assemblies with gradient and fiber optics elements. Keywords
ellipsometry, gradient element, surface layer, refractive index, layer thickness.
Введение
В настоящее время разработка и усовершенствование оптико-электронных систем различного функционального назначения (техники связи и систем обработки информации, автоматизированного управления технологических процессов оптического производства, эндоскопов технического и медицинского назначения и т.п.) непосредственно связаны с использованием элементов градиентной и волоконной оптики [1, 2].
Основными физико-техническими характеристиками градиентных элементов являются распределение показателя преломления n(r) и угловой пространственной частоты передачи изображения Q(r) по «аксиальному сечению» r элемента, апертура градана NA, аберрационные коэффициенты hik трансляции изображения [1-7].
Преимущество элементов градиентной и волоконной оптики, используемых в системах трансляции изображения и обработки информации, по сравнению с обычными линзами проявляется при условии реализации распределений показателя преломления по «аксиальному сечению» градана n(r), обеспечивающих минимальные оптические аберрации [3-5]. Вследствие целого ряда причин, возникающих при технологическом процессе изготовления элементов градиентной оптики (selfoc imaging lens, plane-convex, Rod Lens), реализуемое на практике распределение показателя преломления n(r) отличается от идеального (теоретически обоснованного) и подлежит технологическому контролю [4-7].
Для измерений оптических характеристик элементов градиентной оптотехники, полученных различными технологическими способами [8-10], используются методы интерферометрии [4, 5] и эллипсо-метрии [6, 7]. Однако метод интерферометрии «тонкого среза» [4] имеет существенные недостаток: он является разрушающим и трудоемким. Точность этого метода (10-25%) определяется точностью локального измерения доли интерференционной полосы и недостаточна для граданов с дифракционно-ограниченными параметрами и блоков граданов. Увеличение толщины «среза» образца d^. по сравнению с оптимальной толщиной d^. приводит к возрастанию погрешности измерений из-за искривления траектории лучей в неоднородной среде [4].
Метод иммерсионной интерферометрии [5] основан на сравнении показателя преломления n(r) в исследуемой точке r градиентного образца с показателем преломления по' оптически однородного эталона. Точность метода определяется, главным образом, той точностью, с которой известен показатель преломления эталона по, а также величиной ошибки в измерении дробной части интерференционной полосы при определении ее смещения, связанной с возможными изменениями показателя преломления иммерсии и ее температурной нестабильности [5].
Методы отражательной эллипсометрии в основном разработаны для анализа поляризационно-оптических характеристик неоднородного поверхностного слоя (ПС), образующегося при различной технологической обработке поверхности детали, изготовленной из оптически однородного по своим
свойствам материала [11—14]. Для определения оптических постоянных (показателя преломления n и показателя поглощения k) материала при наличии на поверхности образца шероховатости или неоднородного слоя используются методы иммерсионной [15-17] и многоугловой эллипсометрии [18-20]. Однако предложенные методы эллипсометрии не снимают в полной мере научно-технических проблем, связанных с определением физико-технических характеристик элементов градиентной оптики по поля-ризационно-оптическим параметрам отраженного светового пучка, учитывающих влияние поляризаци-онно-оптических свойств неоднородного ПС.
Цель настоящей работы состояла в усовершенствовании метода отражательной эллипсометрии, предназначенного для измерения физико-технических характеристик градиентных элементов при наличии на поверхности оптической детали неоднородного ПС.
Использование новых методологических подходов в эллипсометрии неоднородных оптических сред для отражающей системы «неоднородный поверхностный слой - неоднородная подложка» позволит усовершенствовать поляризационно-оптические методы технологического контроля при изготовлении элементов градиентной оптотехники и выйти на качественно более высокий уровень измерения физико-технических характеристик оптических узлов с элементами градиентной и волоконной оптики.
Особенности методов измерения и расчета основных эллипсометрических параметров отраженного светового пучка
Одна из основных задач, которую необходимо решить для практического применения метода эллипсометрии при анализе поляризационно-оптических свойств отражающих систем, состоит в установлении функциональной связи между параметрами эллипса поляризации отраженного светового пучка и оптическими характеристиками исследуемой оптической системы. Для описания эллипса поляризации отраженной световой волны используются следующие параметры [11]: азимут большой полуоси эллипса X [0; п]; азимут восстановленной линейной поляризации ¥ [0; п/2]; эллиптичность у = ± arctg(b/a), где b, а
- соответственно малая и большая полуоси эллипса [- п/4; + п/4]; разность фаз А между p- и ^-компонентами поляризованного светового пучка. Поскольку основные эллипсометрические параметры А, ¥, X, Y связаны между собой тригонометрическими соотношениями
cos2¥ = - cos2y • cos2x, tgA = tg2y • sin2x, (1)
sin2y = sin 2¥ • sin 2A, tg2x = -tg2¥ • cos A, то для полного описания состояния поляризации отраженного светового пучка необходимо измерить только два параметра x, Y или ¥, А отраженного светового пучка.
Поскольку в исследуемых объектах изменение показателя преломления происходит по двум независимым направлениям - по сечению градиентного элемента n(r) и по глубине неоднородного ПС n(z), то используются два метода измерения поляризационных характеристик отраженного светового пучка -метод сканирующей и многоугловой эллипсометрии.
На рис. 1 приведена обобщенная оптическая схема измерения эллипсометрических параметров от -раженного светового пучка. В предложенной схеме измерений используются два канала оптико-электронной измерительной системы. По первому каналу применяется схема измерений поляризационных параметров ¥ и А по методу компенсационной (нулевой) эллипсометрии PSCA2; по второму каналу
- по фотометрической схеме PSA\ измерений поляризационных параметров х и y. Такое разделение каналов в оптико-электронной схеме прибора позволяет проводить измерение поляризационных параметров отраженного светового пучка ¥, А и y, X раздельно по каждому каналу и, тем самым, повысить точность эллипсометрических измерений.
Рис. 1. Оптическая схема экспериментальной установки для измерения основных эллипсометрических параметров по методу компенсационной (нулевой) эллипсометрии РЭОЛ2 и фотометрическому методу РБЛь Р - поляризатор; 8 - объект измерения; С - компенсатор; А1, А2 - анализатор; Д - делитель светового пучка; 01, 02 - фотодетектор
Отметим, что под ^-компонентой светового пучка подразумевается компонента электромагнитного поля Е, когда электрический вектор Е^ ориентирован перпендикулярно плоскости падения светового пучка (плоскость ХОТ?). В другой составляющей вектора электромагнитного поля Е, т.е. /»-компоненте, электрический вектор Ер ориентирован параллельно плоскости падения. Ось Т направлена по нормали к границе раздела сред (плоскостьХОТ), угол падения светового пучка ф (рис. 1).
В фотометрической схеме эллипсометрии РБА1 может использоваться метод «стандартных положений» поляризующих элементов и метод «вращающего анализатора» [21, 22]. В первом методе измеряемыми параметрами являются интенсивности отраженного светового пучка /0, /45, 190 при азимуте падающего линейно поляризованного света tg¥п и азимутах анализатора А1 = 0°, А2 = 45°, А3 = 90° [21]. Поляризационные параметры отраженного светового пучка определяются по формулам
У=tgтп-,//07745 , ^д = (2/45 -¡0 -/^^./т^т;;. (2)
Практическая реализация фотометрических измерений интенсивностей 7(А) при использовании метода «вращающего анализатора» с цифровым Фурье-детектированием представлена в [22]. В этом случае зависимость интенсивности светового пучка от азимута анализатора /(А) можно записать как
7(А) = а0 + а2 - ^2А + Ь2 - sin2A , (3)
■ (4)
1 NT 2 N 4л . , 2 NT
= — > I , a2 = — > I ■ cos — i , b2 = — > I
\T ¿—i ' 2 \T i—1 1 N 2 AT 1
N£f 1 2 NtT
Nt1
• sin
n'
Поляризационные параметры светового пучка определяются по формулам
г- —.1/2
J О - О 1 / О 1 <
л 1 b2 X = ---arctg— , tgy = 4 2 a2
a0 - (a22 + b22)1/2
a0 + (a22 + b22)1/2
... (5)
a0 I (a2 1 b2 ) _
Если в методе фотометрической эллипсометрии обычно определяются два параметра эллипса поляризации x и у, то для измерения двух других параметров эллипса поляризации ¥ и Д широко используется метод компенсационной (нулевой) эллипсометрии PSCA2 (рис. 1). Все эти параметры связаны соотношениями (1)-(5).
В компенсационной (нулевой) схеме элллипсометрии PSCA2 на отражающую поверхность объекта S подает линейно поляризованный свет с такой ориентацией плоскости поляризации относительно плоскости падения, при которой после отражения от объекта S свет эллиптически поляризуется так, чтобы, пройдя компенсатор С, он имел бы линейную поляризацию. В этом случае отраженный световой пучок может быть погашен с помощью поворота анализатора на азимутальный угол А. Линейное состояние поляризации на входе оптической системы Р и азимут анализатора А в положении «гашения» светового пучка определяются из следующих выражений [11]:
tgP = Р, •[Р, - С • tg(АС)] • [ • tgC + tg(A - С)]-1, (6)
р, = tg¥, • 'Д, = R(pVR(,:i, Рс = tg¥, • exp'5c = T{c p]/T^, (7)
где р,, рс - эллипсометрические соотношения и поляризационные характеристики (Д, ¥) объекта измерений (S) и компенсатора (C) соответственно. Приборная реализация метода компенсационной эллипсометрии представлена в [23].
Таким образом, применение двухканальной схемы измерений поляризационных параметров отраженного светового пучка позволяет объединить преимущества методов фотометрической и компенсационной (нулевой) эллипсометрии и устранить недостатки, присущие каждому из этих методов.
Теория метода эллипсометрии неоднородных оптических сред
В эллипсометрии однородных сред амплитудные коэффициенты отражения Ro(p,,) для р- и ,-компоненты поляризованного света определяются формулами Френеля [24], которые для наклонного падения светового пучка взаимосвязаны и определяются соотношениями (8)-(10):
R(Р) = R0,). R) - cos 2ф) • (1 - R0,) • cos 2ф)-1, (8)
Ф = R(s) • (р + cos 2ф) • (1 + р • cos 2ф)-1, (9)
р = R(p) /R(,) = tgT • exp(iД), (10)
где ¥ - азимут линейной восстановленной поляризации; Д - разность фаз между р- и ,-компонентой поляризованного света; Ф - вспомогательный параметр для вычисления амплитудных коэффициентов отражения взаимно ортогональных поляризаций;
R0^), r(p^ - коэффициенты отражения для границы раздела однородных сред и неоднородной подложки; ф - угол падения светового пучка. Для однородной подложки значение параметра |Ф| ~ 1 (соотношение (9)). При этом показатель преломления n и показатель поглощения к для геометрически плоской границы раздела однородных сред определяются по следующей формуле [11]:
N = П - ik = Пв • БШф ■
1+1т-Р I2 •
(11)
где «В - показатель преломления внешней среды.
Для неоднородной подложки, т.е. оптически неоднородной среды, при наличии на поверхности элемента ПС соотношение (8) для двух взаимно ортогональных компонент поляризованного излучения не выполняется. При этом значения оптических параметров « и k зависят не только от угла падения светового пучка ф, но и способа их вычисления по аналогам формул Френеля, приведенным в [11, С. 31].
В общем случае рассмотрение поляризационно-оптических свойств отражающей системы «внешняя среда - поверхностный слой - неоднородная подложка» проводится по двум взаимно ортогональным направлениям:
1. по «аксиальному сечению» градана п(г), т.е. на границе раздела «ПС - неоднородная подложка», на глубине ПС г^-да;
2. по глубине ПС значение «(г) в каждой локальной точке г градана, т.е. на границе раздела «внешняя среда - ПС».
При физико-математическом моделировании оптического профиля неоднородной отражающей системы удобно использовать следующее математическое выражение для показателя преломления ПС по его глубине:
«(г, г) = по (г) + (п(0, г) - по (г)) • ехр(-г / й), (12)
где «(г,г) - распределение показателя преломления по глубине г ПС в локальной точке г градиентного элемента; по(г) - значение показателя преломления ПС в локальной точке градана г по его «аксиальному сечению», т.е. перпендикулярно нормали к поверхности (ось 2 рис. 1) при г^-да; п(0,г) - значение показателя преломления ПС в локальной точке градана г при г = 0; й - толщина ПС.
В глубине ПС, т.е. когда г^да, изменение показателя преломления «(г) по радиусу г градана описывается следующими зависимостями [4, 5]:
п(г) = по (0) • (1 - [1 - («2/«о (0))2] • (г • а)а )1/2, (13)
О (г ) =
л ■ 1
д п
-д п
(14)
Д« = [По2 (0) - «2 ] • п0-2 (0), НА = [а /(а + 2)] • [«2 (0) - «2 Г , (15)
где по(0) - показатель преломления на оптической оси ОО/ градана, т.е. при г = 0; «2 - показатель прелом-
ления на краю градана, т.е. при г = а, где а - эффективный радиус, отсчитываемый от оптической оси ОО градана; показатель степени «а» может принимать значения 1 < а < да: для параболического закона изменения «(г) величина а = 2, для ступенчатого оптического профиля «(г) значение а^-да; НА - апертура градиентного элемента; 0.(г) - распределение угловой пространственной частоты переноса изображения по «аксиальному сечению» градана г.
Определяя значения показателей преломления в локальных точках градана «(^), «(0), «2, по формулам (13)-(15) несложно определить значение параметра «а» профиля показателя преломления градана, распределение угловой пространственной частоты передачи изображения 0.(г) по «аксиальному сечению» градана г, значение апертуры градиентного элемента НА, а по формулам, приведенным в [4], -аберрационные коэффициенты трансляции изображения к^. Тем самым можно определить все физико-технические характеристики градиентного элемента [1-7].
При анализе поляризационно-оптических характеристик отражающей системы «неоднородный поверхностный слой - неоднородная подложка» используется уравнение эллипсометрии в приближении Борна [11, 12]. Для профиля показателя преломления неоднородного ПС, описываемого формулой (12), оптические параметры слоя определяются по формулам, приведенным в работах [11, 12, 25]: Е = Яе(р/Ро) -1, В = Яе (р/Ро), й = Е/(В • ио • к),
5« = 2В • (1 + ио2 • к2 • й2) • (А • к • й)-1, 5« = «(0, г) - «о(г),
(16)
р = гяТ • ехр(/Д), Ро = гяТо • ехр(/До), ио = [По2(г) - эш2 ф]1/2, к = 2п / X, А = -(1 + tgT)3 • [4 • (1 - tgT) • tgT • эЬ2ф • tg2ф]-1. где Е и В - вспомогательные параметры, введенные в решение обратной задачи эллипсометрии; й - толщина слоя; к - волновое число; ио - адмиттанс однородной среды локальной точки градана радиуса г; А - параметр, зависящий только от условия проведения многоугловых измерений; 5« - отклонение показателя преломления на поверхности градана «(0,г) при г=0, от показателя преломления в объеме неоднородной среды «о(г), р, ро - эллипсометрическое отношение; ¥о - азимут линейной восстановленной
поляризации; Д, До - разность фаз между ортогональными компонентами поляризованного излучения для неоднородной и однородной подложки в локальной точке градана r соответственно.
Теоретический анализ неоднородных отражающих систем показал [25], что наличие на поверхности градана неоднородного слоя позволяет в методе эллипсометрии перейти от зависимости ¥ = ¥[по(г)] к зависимости Д = Д[по(г)] и тем самым повысить чувствительность метода и дать более детальный анализ аберраций угловой пространственной частоты передачи изображений ^(r). При этом изменение показателя преломления 5no(r) по радиальному сечению градана прямо пропорционально отклонению фазового сдвига Д от фазового сдвига однородной подложки До.
S«(r)9. const = [n(r) - n0(r0)] * (5Д)г= var = [A(r) - До],
Sn(z)r=const = [n(z) -По(Го)] * (5Д)2_ = [А(ф) - До]2. (17)
Отклонение 5n(z) показателя преломления неоднородного слоя n(z) от объемного значения no в локальной точки градана r¡ пропорционально квадрату отклонения фазового сдвига Д от фазового сдвига однородной подложки До (формула (17)). Если измерение эллипсометрических параметров ¥(r) и Д(г) проводится по сечению градана методом локального зондирования поляризованным световым пучком поверхности, имеющей неоднородный ПС, для которого по сечению градана (в пределах погрешности измерения поляризационных углов Sv и SA) средние значения толщины dcp и показателя преломления ncp неоднородного ПС для различных видов профиля n(z) можно считать близкими, то из уравнения эллипсометрии в приближении Борна несложно получить следующие соотношения [25]:
5(Д) = Д(г. ) - Д(г.+i) = Aу • [n(r ) - n(rj+i)],
(18)
Aj = -(1 + tg¥j )3 • [4 • (1 - tg¥j) • tg¥j • sin2ф • tg2ф]-1.
При углах падения светового пучка ф, где параметр ¥/ принимает малые значения и, согласно (18), коэффициент пропорциональности Aj резко возрастает, т.е. увеличивается чувствительность измерений по параметру Д, но возрастает также среднеквадратическая погрешность его измерения [7]. В связи с этим измерение разности фаз Д(г,) по сечению градиентного элемента должно проводиться при углах падения ф, при которых сохраняется высокая чувствительность метода эллипсометрии к измерению параметра Д(г,) и незначительна его погрешность измерения Sj,A, которая определяется по формулам, приведенным в работе [11].
Для определения оптического профиля градана следует сопоставить теоретические зависимости профиля показателя преломления п(т)(г), рассчитанные по формуле (6) при различных значениях а, с экспериментальной зависимостью п(э)(г), полученной из результатов эллипсометрических измерений поляризационных параметров Д и ¥. Сопоставление проводится по значениям функции ошибок F(a) принятия решений о зависимости n(r) по методу максимального правдоподобия, описанному в работах [20-26]:
F(а) = — £[П(3)(r-)-П2(т)(г)]2 , (19)
' ' Mtí (Sn)2
где Sn - погрешность измерения показателя преломления; М - количество измерительных ситуаций, определяемых количеством локальных значений п(э)(г,). Вычисления по формуле (19) для различных видов распределений показателя преломления п(э)(г,} позволяют принять за истинное распределение теоретическую зависимость показателя преломления п(т) (r), при которой значение функции F(a) < 1. В этом случае вероятность ошибки второго рода минимальна при наименьшем количестве проводимых эллипсометри-ческих измерений.
Эксперимент и обсуждение результатов
В качестве объектов исследования были использованы элементы градиентной оптики, полученные методом ионообменной диффузии [8, 9], имеющих эффективный радиус а = 6-12 мм. Полирование образцов осуществлялось по стандартной технологии [11]. Измерение эллипсометрических параметров ¥ и Д проводились на приборе ЛЭФ-3М (X = 0,6328 мкм) по двухзонной методике эллипсометрических измерений [11] при углах падения светового пучка ф = 45-60° и при диаметре светового пучка D < 1 мм. Погрешность измерений этих параметров не превышала Sv = 1 и SA = 2-3.
По измеренным значениям основных эллипсометрических параметров Д и ¥ при угле падения светового пучка ф = 45° в различных локальных точках градиентного элемента X по формуле (11) были рассчитаны эффективные значения показателя преломления n* и показателя поглощения к* (рис. 2). Из представленных результатов видно, что показатель поглощения к* изменяется закономерно от значений к, соответствующих прозрачным диэлектрикам к<10-3, до значений, соответствующих полупроводникам (в частности, для кремния к ~ 0,023).
Это указывает на несоответствие модели геометрически плоской границы раздела для элементов градиентной оптики. Исходя из этого, анализ поляризационно-оптических свойств элементов градиентной оптики должен проводиться в рамках модели «неоднородный поверхностный слой -неоднородная подложка» [25].
к-103 20
15 10
5 1
/
Г
* >
/ * 2 2
I
ш
/о
!
.К*
X ♦ • • •
\ I
» * —
Ч..? х -*
1,5240
1,5220
0
4
12
1-1
16 X, мм
Рис. 2. Изменение значений показателя преломления п (кривая 1) и показателя поглощения к (кривая 2)
по сечению градана, рассчитанных по формуле (11)
На рис. 3 показано изменение поляризационных параметров ¥(г) (кривая 2) и |5Д| (кривая 1), где |5Д| = |Д(г)- До|, Д(г) - экспериментально измеренное значение разности фаз для взаимно
ортогональных поляризаций; До - разность фаз для геометрически плоской границы раздела сред, значение которой при ф = 60° равно До = 0°.
|8Д|, угл. мин 100
75
ОС А
50
25
-6
Дп-104 75
50
25 ¥
з
5°06' 5°00' 4°56'
4 г, мм
Рис. 3. Изменение поляризационных параметров отраженного светового пучка и показателя преломления по «аксиальному сечению» градана. Кривая 1 - | 6Д|; кривая 2 - Ф(г); кривая 3 - бп(г); кривая 4 - Фо(г);
за ось симметрии градана принято значение г = 0, диаметр градана составлял 12 мм; различные физические параметры в локальных точках градана «г» помечены различными символами: изменение показателя преломления бп(г), - «х», а изменение фазового сдвига | 6Д | - «о»
При малом изменении параметра ¥(г) по «аксиальному сечению» градана изменение параметра 5Д(г) носит регулярный характер, а изменение показателя преломления 5«(г) (кривая 3) по радиальному сечению градана пропорционально изменению разности фаз 5Д(г). Расчет показал, что в локальной точке градана при г1 =а значение показателя преломления равно «ог=а = 1,5256 при погрешности измерения 8п = ± 2 10-4.
По показателю преломления «о(г) материала градиентного элемента рассчитывались истинные значения азимута линейной восстановленной поляризации ¥о для геометрически резкой границы раздела в локальных точках «внешняя среда - неоднородная среда», т.е. без учета оптических свойств ПС (кривая 4). Затем по формулам (16) определялись оптические характеристики ПС: показатель преломления «(0,г) на поверхности градана и толщина слоя й (таблица).
Погрешность в измерении показателя преломления «(0,г) составила 8« < ± 2-10 4, а толщины слоя -< ± 10 нм. Расчет параметров ПС при локальном значении показателя преломления градана «(0,г) = «о(г=а) + 5«(г), т.е. при учете распределения показателя преломления «(г) по радиальному сечению градана, дает хорошую сходимость результатов, что подтверждает правомерность используемой ме-
«
1
тодики определения оптического профиля градана и топографии поверхности, оцениваемой по параметрам ПС.
r, мм d, нм n(0,r) na(r)
- 6 143 1,5278 1,5254
- 4 406 1,5376 1,5278
0 145 1,5434 1,5299
4 492 1,5383 1,5284
6 140 1,5290 1,5258
Таблица. Оптические параметры поверхностного слоя по «аксиальному сечению» градана
Для данного распределения показателя преломления градана по формулам (13)-(15) были получены следующие физико-технические параметры градана: а = 1,76; ß(r=0) = 0,07678 мм1, NA = 0,04655.
Заключение
На основе разработанного методического подхода к эллипсометрическому анализу оптических характеристик отражающей системы «поверхностный слой - неоднородная подложка», проводимому на основе результатов исследований поляризационно-оптических характеристик элементов градиентной оптики, можно сделать следующие выводы.
Модель геометрически плоской границы раздела однородных сред в локальной точке градиентного элемента не соответствует физическому содержанию исследуемой отражающей системы.
Показано, что при наличии на поверхности градиентного элемента неоднородного слоя азимут линейной восстановленной поляризации практически не зависит от изменений показателя преломления в глубине поверхностного слоя, а разность фаз между взаимно ортогональными компонентами поляризованного светового пучка полностью определяет распределение показателя преломления по «аксиальному сечению» градиентного элемента.
Разработанный метод эллипсометрии неоднородных сред позволяет определять распределение показателя преломления по «аксиальному сечению» градиентного элемента при наличии на поверхности детали неоднородного слоя.
Предложена и обоснована эллипсометрическая методика определения показателя преломления и толщины поверхностного слоя при изменении показателя преломления по «аксиальному сечению» градиентного элемента в объеме материала.
Литература
1. Васильева В.И., Петровский Г.Т., Черных В.Д., Ягмуров В.Х. Гибкие регулярные элементы для эндоскопии // Оптический журнал. 1994. Т. 61. № 12. С. 55-57.
2. Архипова Л.Н., Ивашевский С.Н., Карапетян Г.О. и др. Градиентная оптика для медицинских эндоскопов // Оптический журнал. 1994. Т. 61. № 12. С. 51-54.
3. Ильин В.Г., Меланьина Т.М., Ремизов Н.В. Оценка распределения показателя преломления при измерении угловой пространственной частоты передачи изображения в граданах // Оптика и спектроскопия. 1984. Т. 57. № 5. С. 933-936.
4. Ильин В.Г., Ремизов Н.В. Интерференционный метод измерения распределения показателя преломления в передающих изображение градана // Письма в ЖТФ. 1984. Т. 10. № 2. С. 105-110.
5. Ильин В.Г., Карапетян Г.О., Полянский М.Н. Измерение локальных значений показателя преломления неоднородных сред // Журнал прикладной спектроскопии. 1978. Т. 28. № 1. С. 160-163.
6. Демидов И.В., Лисицын Ю.В., Михайлов Ю.В., Храмцовский И.А. Определение пространственной частоты передачи изображений градиентными оптическими элементами методами эллипсометрии // Сборник тезисов V Международной конференции "Прикладная оптика-2002". СПб., 2002. Т. 1. С. 156.
7. Землянский В.С., Горляк А.Н., Степанчук А.А., Храмцовский И.А. Эллипсометрический метод технологического контроля элементов лазерной техники и градиентной оптики // Научно-технический вестник СПбГУ ИТМО. 2007. № 9 (43). С. 81-87.
8. Ильин В.Г., Карапетян Г.О., Ремизов Н.В. Стекла для градиентной оптики // Сборник тезисов VII Всесоюзного симпозиума "Оптические и спектральные свойства стекол". Л., 1989. С. 234-235.
9. Стеклообразное состояние / Под ред. Е.А. Порай-Кошица. Л.: Наука, 1988. 170 с.
10. Серков М.М., Кондратьев Ю.Н. Градиентные среды на основе кварцоидов // Физика и химия стекла. 1984. Т. 10. № 1. С. 160-163.
11. Пшеницын В.И., Абаев М.И., Лызлов Н.Ю. Эллипсометрия в физико-химических исследованиях. Л.: Химия, 1986. 152 с.
12. Храмцовский И.А., Пшеницын В.И., Каданер Г.И., Кислов А.В. Учет оптических характеристик поверхностного слоя при определении коэффициентов отражения и пропускания прозрачных диэлектриков // Журнал прикладной спектроскопии. 1987. Т. 46. № 2. С. 272-279.
13. Землянский В.С., Храмцовский И.А., Горляк А.Н., Степанчук А.А. Методы эллипсометрического анализа поляризационно-оптических свойств неоднородных поверхностных слоев элементов оптоэлек-троники // Оптика и спектроскопия. 2008. Т. 105. № 2. С. 346-351.
14. Прокопенко В.Т., Храмцовский И.А., Землянский В.С., Лисицын Ю.В., Секарин К.Г. Эллипсометрия оптических соединений элементов оптоэлектроники // Изв. вузов. Приборостроение. 2008. Т. 51. № 10. С. 57-65.
15. Любинская Р.И., Мардежов А.С., Швец В.А. Исследование неоднородных структур с использованием иммерсионных эллипсометрических измерений // Сб. Эллипсометрия: теория, методы, приложение / Под ред. А.В. Ржанова и Л.А. Ильина. Новосибирск: Наука, 1987. С. 59-67.
16. Гайнер А.В., Сурдутович Г.И. О критериях шероховатости и неоднородности поверхностного слоя образцов при эллипсометрических измерениях // Сб. Эллипсометрия: теория, методы, приложение / Под ред. А.В. Ржанова и Л.А. Ильина. Новосибирск: Наука, 1987. С. 50-53.
17. Мардежов А.С., Швец В.А. К вопросу интерпретации результатов эллипсометрических измерений // Сб. Эллипсометрия: теория, методы, приложение / Под ред. А.В. Ржанова и Л.А. Ильина. Новосибирск: Наука, 1987. С. 83-86.
18. Golyak A.N., Cramtsovsky I.A. The ellipcometry of the roug suface on an inhomogeneous substrate // Proc. of 4th Int. Conf. on Spectroscopic Ellipsometry (ICSE 4). 2007. Stockholm, Sweden, 2007. Р. 74.
19. Gorlyak A.N., Khramtsovskij I.A. Diagnostics of the physicochemical state of the surface of optoelecronics elements by ellipsometry method // Proc of. 5th Workshop Ellipsometry. Zweibruecken, Germany, 2009. P. 63.
20. Gorlyak A.N., Khramtsovsky I.A. Research of the optical characteristics of surface layers by methods of an ellipsometry and of the resonance excitation of waveguide modes // Proc. 6th Int. Conf. on Spectroscopic Ellisometry (ICSE VI). Kyoto, Japan, 2013. P. 195.
21. Маслов В.П., Одарич В. А. Эллипсометрические исследования механически полированных образцов некоторых оптических стекол // Оптико-механическая промышленность. 1983. № 3. С. 60-61.
22. Алексеев С.А., Горляк А.Н., Прокопенко В.Т., Сухорукова М.В., Устинов С.Н., Фэн Ли Шуан. Автоматический спектральный эллипсометр // Приборы и техника эксперимента. 1996. № 3. С. 166.
23. Алгазин Ю.Б., Иощенко Н.Н., Леоненко А.Ф., Панькин В.Г., Рыхлитский С.В., Свиташев К.К. Лазерный фотоэлектрический эллипсометр ЛЭФ-3М-1 // Приборы и техника эксперимента. 1987. № 6. С. 204.
24. Azzam R.M.A. Direct relation between Fresnnel's interface reflection coefficients for the parallel and perpendicular polarizations // Journal of the Optical Society of America. 1979. V. 69. N 7. P. 1007-1016.
25. Пшеницын В.И., Храмцовский И.А., Крылова Н.А., Подсекаев А.В., Туркбоев А. Эллипсометрия отражающей системы «неоднородный слой - неоднородная подложка» // Сборник «Методы прикладной математики в транспортных системах» / Под ред. Ю.М. Кулибанова. СПб.: ГУВК, 2000. Вып. 3. С. 3237.
26. Горляк А.Н., Солонуха В.М., Храмцовский И.А. Применение метода секционирования в эллипсомет-рии неоднородных отражающих систем // Научно-технический вестник СПбГУ ИТМО. 2014. № 3 (91). С. 24-29.
Горляк Андрей Николаевич
Храмцовский Игорь Анатольевич
Солонуха Владимир Михайлович
кандидат технических наук, ведущий научный сотрудник, Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В.И.Ульянова (Ленина), Санкт-Петербург, 197376, Российская Федерация, [email protected]
кандидат технических наук, ведущий инженер, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
аспирант, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
Andrey N. Gorlyak
Igor A. Khramtsovsky Vladimir M. Solonukha
PhD, leading scientific researcher, Saint Petersburg State Electrotechnical University, Saint Petersburg, 197376, Russian Federation, [email protected]
PhD, leading engineer, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
postgraduate, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
