ЛАЗЕРНЫЙ ИНТЕРФЕРОМЕТР С ОСТРОСФОКУСИРОВАННЫМ ПУЧКОМ В СИСТЕМЕ ТРЕХ СВЯЗАННЫХ ИНТЕРФЕРОМЕТРОВ
Д.В. Лякин, В.В. Лычагов, В.П. Рябухо, М.И. Лобачев Институт проблем точной механики и управления РАН Саратовский государственный университет
Аннотация
Рассматривается система связанных интерферометров - измерительного и двух вспомогательных, построенных по схеме Майкельсона, для контроля пространственного положения (формы) объекта или параметров слоистых структур. Измерительный интерферометр функционирует на эффекте согласования интерферирующих полей и формирует интерференционный импульс совпадения фокуса зондирующего лазерного пучка с контролируемой поверхностью. Вспомогательный низкокогерентный интерферометр формирует опорный интерференционный импульс, относительно которого определяется положение измерительного импульса в шкале разности хода. Для этого служит вспомогательный лазерный интерферометр. Рассматривается теория интерферометров, схемные решения и алгоритмы компьютерной обработки сигналов. Приводятся экспериментальные результаты, анализируется точность измерений.
Введение
Лазерные системы широко используются для размерного контроля, пространственного положения поверхности объекта и границ раздела сред в слоистых структурах, макроформы и микроструктуры поверхности и т.п. [1-10]. Для решения этих измерительных задач используются различные методы и подходы, основанные, в основном, на высокой степени пространственной и временной когерентности лазерного излучения, и, как следствие, на возможности использования острой фокусировки лазерного пучка, эффектов интерференции и дифракции света.
Метод точной фокусировки лазерного пучка [6, 7, 10] нашел широкое практическое использование в силу относительной простоты реализации и достаточно высокой точности измерений. Этот метод относится к так называемым амплитудным методам оптических измерений и основан на регистрации сигнала совпадения отражающей поверхности с фокусом зондирующего лазерного пучка. Для этого используется точечная апертура фотоприемника, оптически сопряженная с фокусом зондирующего лазерного пучка. Поэтому данный метод чаще всего называют конфокальной микроскопией [11-14]. Точность метода определяется числовой апертурой используемой фокусирующей оптической системы и длиной волны излучения. Для увеличения точности измерений в конфокальной микроскопии используют когерентный опорный пучок и, соответственно, интерференционный эффект, что позволяет регистрировать положение отражающей поверхности с точностью до долей длины волны [11, 12, 15-17].
В конфокальной микроскопии использованием точечной апертуры фотоприемника исключается влияние фоновой засветки, что позволяет получать сигнал от объемной слоистой структуры рассеивающих объектов.
В определенном смысле альтернативным методом определения пространственного положения объекта и параметров слоистых структур является
метод оптической низкокогерентной интерферометрии [18-23]. В этом методе используется широкополосный источник света с малой длиной временной когерентности излучения, но с достаточно малой излучающей поверхностью и, как следствие, с высокой степенью пространственной когерентности излучения. Метод заключается в регистрации интерференционного сигнала совпадения длин оптических путей в опорном и объектном плечах интерферометра. Этот сигнал совпадения представляет собой интерференционный импульс временной когерентности в шкале разности хода с высокочастотной несущей, определяемой длиной волны, и огибающей, определяемой модулем функции временной когерентности излучения. Ширина огибающей этого импульса в шкале разности хода определяется длиной временной когерентности используемого излучения, а период осцилляций несущей - средней длиной волны излучения.
В настоящей работе рассматривается новый интерференционный метод точной фокусировки лазерного пучка, который сочетает в себе принципы конфокальной микроскопии и низкокогерентной интерферометрии. Метод основан на эффектах согласования и интерференции предметного и опорного волновых полей и формировании интерференционного импульса совпадения фокуса зондирующего пучка с поверхностью контролируемого объекта. Рассматривается практическая реализация этого метода с использованием системы трех связанных интерферометров с лазерными и низкокогерентным источниками света для определения пространственного положения поверхности объекта и параметров слоистых структур.
1. Система трех связанных интерферометров
Рассматриваемая в настоящей работе интерференционная система включает в себя три связанных интерферометра, построенных по схеме Майкельсо-на (рис.1).
Рис. 1. Система трех связанных интерферометров: Int.l - основной сканирующий интерферометр с остросфокусированным лазерным пучком: LD1 - лазерный диод; BS1 - делитель пучка; SP - сканирующая платформа;
RM - опорное зеркало; S - объект контроля; MO1 и MO2 - микрообъективы; PD1 - фотодетектор; u1 - выходной сигнал интерферометра; Int.2 - вспомогательный интерферометр с низкокогерентным источником света: LED - светодиод;
BS2 - делитель пучка; M1 - опорное зеркало; M2 - зеркало, жестко закрепленное на сканирующей платформе SP интерферометра Int.1; PD2 - фотодетектор; u2 - выходной сигнал интерферометра; Int.3 - вспомогательный интерферометр с параллельным лазерным пучком на базе оптических элементов низкокогерентного интерферометра Int.2: PD3 - фотодетектор; u3 - выходной сигнал интерферометра
Первый, основной, сканирующий интерферометр Int. 1 с остросфокусированным лазерным пучком используется для зондирования контролируемого объекта и формирования интерференционного импульса u1 совпадения фокуса лазерного пучка с поверхностью объекта или отражающей границей слоистой структуры.
Второй, вспомогательный, интерферометр Int.2 с низкокогерентным источником используется в качестве следящей системы для формирования интерференционного импульса u2 фиксации пространственного положения первого интерферометра. Интерференционные импульсы сканирующего (u1) и следящего (u2) интерферометров разнесены во времени и, соответственно, разнесены в шкале разности хода интерферирующих волн (рис. 1).
Третий, также вспомогательный, интерферометр Int.3 с квазипараллельным лазерным пучком используется в качестве средства измерения смещения сканирующего интерферометра - фактически для измерения расстояния между импульсами сканирующего и следящего интерферометров в шкале разности хода.
На выходе этого интерферометра формируется практически амплитудно немодулированный сигнал u3 , интерференционные осцилляции которого
используются в качестве шкалы отсчета.
Ниже рассмотрены принципы действия и функциональное назначение каждого интерферометра и всей системы в целом.
2. Интерферометр с остросфокусированным предметным пучком и интегрирующей апертурой фотоприемника Для формирования импульса совпадения контролируемой поверхности и фокуса зондирующего лазерного пучка в рассматриваемой системе используется интерферометр, функционирующий на принципах согласования объектного и опорного волновых полей на протяженной апертуре фотоприемника. На выходе такого интерферометра формируется импульсный интерференционный сигнал с максимальным коэффициентом модуляции интерференционных осцилляций при максимальном согласовании интерферирующих полей [24]. Эффект рассогласования интерферирующих полей при продольном смещении контролируемой поверхности относительно фокуса зондирующего лазерного пучка зависит от числовой апертуры пучка NA и длины волны X. Чем больше NA, тем более узкий интерференционный импульс в шкале разности хода формируется на выходе интерферометра.
На рис. 2 представлена оптическая схема сканирующего интерферометра с остросфокусированным предметным (зондирующим) лазерным пучком. Излучение лазерного диода LD1 разделяется светоделителем на опорный и предметный пучки, которые фокусируются, соответственно, на опорное зеркало RM и на поверхность контролируемого объекта S. Отраженные опорный и предметный пучки света с использованием поворотной призмы P направляются на фотодетектор PD1. Весь интерферометр за ис-
ключением объекта 8 располагается на сканирующей платформе 8Р, совершающей колебания вдоль оси Ъ с достаточно большой амплитудой 10 и 1 мм и низкой частотой /0 и 1 Гц. Таким образом изменяется разность хода волн в интерферометре и происходит расфокусировка лазерного пучка на поверхности объекта.
U (Az) ~ //1 ( p)d2 p = // [ IR ( p) + Is ( p, Az, Az,) + +2 Ar ( p) As ( p, Az, Az,) +
+ cosI [2Az + 2 Az, ] + ^2Az01 + A^RS (p, Az, Az,) l Я, Я, j
d2 P , (1)
Uj(Az)^|r
Рис. 2. Оптическая схема лазерного интерферометра с остросфокусированным предметным пучком и протяженной интегрирующей апертурой фотоприемника: LD1 - лазерный диод; BS1 - делитель пучка; MO1, MO2 - фокусирующие линзы (микрообъективы); RM - опорное зеркало; S - объект контроля; P - поворотная призма;
PD1 - фотодетектор; SP - сканирующая платформа, на которой расположен интерферометр; Az - смещение сканирующей платформы;
AD - апертурная диафрагма фотодетекора PD1; u, (Az) - выходной сигнал интерферометра
При совпадении фокусов опорного и предметного пучков света, соответственно, с поверхностью опорного зеркала RM и поверхностью объекта S оптические поля оказываются максимально согласованными на достаточно протяженной апертуре AD фотоприемника PD1. Это означает, что разность фаз интерферирующих полей A<SR (x, y) имеет минимальные изменения в плоскости XY апертуры фотоприемника, а амплитудные распределения этих полей Ar (x, y) и AS (x, y) в максимальной степени совпадают, AR (x, y)/AS (x, y) » const.
Отметим, что полному согласованию полей могут помешать аберрации оптических элементов в плечах интерферометра, несоосность фокусирующих объективов, а также шероховатость поверхности объекта, в результате которой возникает спекл-модуляция предметного волнового поля и комплексная амплитуда этого поля US (x, y) = AS (x, y)exp(S (x, y)) приобретает случайный пространственный характер.
Для выходного фотоэлектрического сигнала интерферометра u,(Az ) в предположении полной взаимной когерентности интерферирующих полей можно записать следующее выражение:
где Дг - смещение интерферометра относительно поверхности объекта, I (р) - распределение интенсивности света в интерференционной картине на апертуре Бр фотоприемника р(х,у); 1К(р), 15(р) -
распределения интенсивности в опорном и предметном полях; Дг1 - начальное расстояние между фокусом предметного пучка и поверхностью объекта (начальная расфокусировка), в предположении, что опорный пучок точно сфокусирован на опорном зеркале; 2Дг01 - оптическая разность хода волн в интерферометре при точной фокусировке пучков в предметном и опорном плечах интерферометра, возникающая, например, при неравенстве оптических путей в плечах интерферометра от светоделителя до фокусирующих свет объективов; ДфЯ5 (р, Дг, Дг1) - разность фаз опорного и предметного полей на апертуре фотоприемника, возникающая при смещении Дг + Дг1 фокуса предметного лазерного пучка относительно поверхности объекта;
- длина волны света лазерного диода ЬБ1.
Первое слагаемое в (1) Ц1К (р)d 2р определяет
постоянную составляющую сигнала фотоприемника, второе слагаемое Ц15 (р, Дг, Аг^2р - медленно
яр
меняющуюся с изменением Дг составляющую. Третье, быстро осциллирующее слагаемое, определяет переменную составляющую с интерференционными осцилляциями с периодом в шкале разности хода Д и огибающей, определяемой эффектом рассогласования интерферирующих полей по амплитудным, Лк (р) и Л5 (р, Дг, Дг1), и/или фазовым ДфЯ5 (р, Дг, Дг1) распределениям на апертуре
фотоприемника.
Влияние эффектов амплитудного и/или фазового рассогласований на огибающую интерференционного слагаемого зависит от положения апертуры фотоприемника относительно плоскости (рис. 3), оптически сопряженной с поверхностью опорного зеркала КМ. При совпадении апертуры фотодетектора ЯР с плоскостью - фокусом опорного пучка света, определяющую роль играет эффект амплитудного рассогласования полей (см., например, [25]). При значительном удалении апертуры фотоприемника от плоскости определяющую роль играет фазовое рассогласование полей при практическом отсутствии амплитудного рассогласования.
. Р(Х,У)
MOI М02
Рис. 3. К определению влияния фазового рассогласования опорного и предметного оптических полей
Однако оба этих эффекта рассогласования и их совместное действие приводят, что подтверждается теоретически и экспериментально, к одинаковым количественным результатам, определяющим огибающую интерференционного сигнала. Действительно, при амплитудном рассогласовании, когда апертура фотоприемника располагается в плоскости SR - в фокусе опорного пучка, который играет роль эффективной апертуры системы, подобно точечному отверстию в конфокальной микроскопии, смещение фокуса зондирующего лазерного пучка Az + Azj относительно поверхности предмета приводит к расфокусировке предметного пучка, увеличению его поперечного сечения и соответствующему уменьшению интенсивности этого пучка в пределах фокуса опорного пучка.
Фазовое рассогласование A<RS (р) доминирует при расположении апертуры фотоприемника вдали от плоскости фокусировки SR опорного пучка [24].
В предположении интерференции сферических опорной и предметной волн для разности фаз A<rs (р) на апертуре фотоприемника (см. рис. 3) может быть получено выражение
(
ДФ RS (р) = :—
А,
1 1
Л
Р
п Az' 2 2 Р '
s
А
(2)
1 Z R
где гк и - расстояния от центров кривизны опорной и предметной волн до плоскости апертуры фотоприемника, Дг ' = г'к - .
Интегрирование по круглой апертуре фотоприемника с диаметром Вр в предположении постоянства амплитуд Лк и опорного и предметного полей в пределах этой апертуры приводит к следующему выражению для третьего слагаемого в (1):
m,(2Az )'
sin(nD2 Az '/8A1zR2 ) ( p Az '/SAizR2 )
(
x cos
— [2Az + 2Az, + 2Az01 ] -А1
Dp
(3)
8А, zR
Az
Используя условие синусов Аббе [27] для хода лучей через фокусирующий микрообъектив в предметном плече интерферометра, можем записать
sin(6 ') « 0' = — sin(0) = — ,
M M
Az ' = (2Az + 2Az, )M 2,
(4)
(5)
где M - поперечное увеличение объектива в предметном плече, 0 и 0 ' - углы прохождения крайних лучей через апертуру объектива в пространстве предмета и его изображения, NA - числовая апертура фокусируемого на объект лазерного пучка, M2 -продольное увеличение.
Подставляя (4), (5) в (3) получаем выражение для переменной составляющей сигнала интерферометра с остросфокусированным зондирующим пучком в следующем виде:
sin
m,(2az )'
Xcos
п "211
2п
— NA2[ 2Az + 2Az, ] 2 Я,
-П- NA2[ 2Az + 2Az, ] 2 Я,
[ 2Az + 2Az, + 2Az01 ]-
(6)
NA2[ 2Az + 2Az, ]
Отметим, что согласно уравнению (6) начальная
разность хода 2Д?01 в данном интерферометре в отличие от низкокогерентного интерферометра не влияет на положение интерференционного импульса. Это определяет важное практическое значение рассматриваемого интерферометра по отношению к низкокогерентному интерферометру - не требуется точно выравнивать оптические пути в плечах. Положение интерференционного импульса в шкале оптической разности хода (или шкале относительного смещения объект-интерферометр) будет определяться значением начальной расфокусировки Д?1.
Из уравнения (6) легко сделать оценку полуширины ^с интерференционного импульса в шкале изменения разности хода 2Дг по спаданию амплитуды интерференционных осцилляций этого сигнала до половины максимальной амплитуды
1,2 ^
L
NA
(7)
Величина Ьс определяет продольное разрешение рассматриваемого интерферометра по аналогии с длиной временной когерентности 1с в низкокогерентной интерферометрии.
При использовании высокоапертурных микрообъективов ЫЛ и 1 в интерференционном импульсе будет проявляться всего несколько осцилляций. Таким образом, ширина импульса рассматриваемого интерферометра может составлять всего 1-3 мкм,
\ --1
V zR
что предопределяет заметно более высокое пространственное разрешение по сравнению с низкокогерентной интерферометрией.
3. Вспомогательный низкокогерентный интерферометр
Для регистрации пространственного положения сканирующего лазерного интерферометра используется интерферометр Майкельсона с низкокогерентным (широкополосным) источником света (например, светодиодом). Одно зеркало этого интерферометра жестко закреплено на сканирующей платформе 8Р лазерного интерферометра с остросфокусиро-ванным зондирующим пучком (рис.1). На выходе низкокогерентного интерферометра формируется сигнал с переменной составляющей и2 (2Дг), интерференционные осцилляции которого модулированы по амплитуде модулем функции временной когерентности излучения источника Г(2Дг):
~2 (2Дг) ~ |Г(2Д2 + 2Дг02 )| >
( 2п _ 2п - . Л х 008 — 2Дг +--2Дг0
ЧХ 2
X 2
(8)
где X 2 - средняя длина волны низкокогерентного источника, 2Дг02 - начальная разность хода волн в плечах интерферометра; 2Дг + 2Дг02 = Д2 - полная разность хода волн в этом интерферометре.
Для модуля Г(Д2) можно использовать приближение в виде гауссовой функции
|Г(Д2) и 10 ехр|-
^2
2
Ч 1с У
(9)
где 1С иХ2/ДХ2 - длина временной когерентности, ДХ 2 - ширина спектрального контура излучения. Таким образом, ширина интерференционного импульса низкокогерентного интерферометра определяется длинной когерентности IС.
Из (8) следует, что при продольном смещении Дг сканирующего лазерного интерферометра на выходе низкокогерентного интерферометра формируется импульс временной когерентности с максимальной амплитудой осцилляций при Дг = -Дг 02 . Этот импульс в рассматриваемой интерференционной системе используется в качестве опорного импульса начала отсчета пространственного положения платформы лазерного сканирующего интерферометра. Таким образом, низкокогерентный интерферометр, вернее, пространственное положение его опорного зеркала М1, служит опорой в определении пространственного положения поверхности контролируемого объекта.
Отметим, поскольку в рассматриваемом низкокогерентном интерферометре используются плоские зеркала, то возможно использование протяженного
источника света, излучение которого пространственно частично когерентно. Влияние малого радиуса поперечной пространственной когерентности излучения рс легко исключается соответствующим юстировочным наклоном неподвижного зеркала интерферометра.
В наших экспериментах в качестве низкокогерентных источников использовались красные (12 ~ 0,64 мкм) светодиоды с длиной временной когерентности 1С и 8 мкм.
4. Вспомогательный лазерный интерферометр для определения величины смещения платформы сканирующего интерферометра
Оптические элементы низкокогерентного интерферометра 1Ш-2 (рис.1) - делительная призма-куб и зеркала, использовались одновременно для вспомогательного лазерного интерферометра с параллельным пучком света. Иными словами, на входе одного интерферометра использовались сразу два источника - светодиод и лазерный диод, а на выходе - два фотоприемника, пространственное положение которых обеспечивало раздельную регистрацию сигналов этих интерферометров. Лазерный канал интерферометра использовался для получения практически немодулированного по амплитуде квазисинусоидального сигнала с известным периодом осцил-ляций, определяемым длиной волны лазера Х3 :
(
М3(2Дг) ~ 008
2п 2п
2Дг +— 2Дг, Х
ЧХ 3
03
3
(10)
где 2Дг03 - как и в (8), начальная разность хода волн
в интерферометре.
Поскольку в системе использовалось не равномерное, а колебательное движение платформы сканирующего интерферометра, то сигнал (10) имеет временную частотную модуляцию, поскольку уравнение смещения платформы можно записать в виде
МО = /0^(2/ + ф0 ), (11)
где 10 и /0, соответственно, амплитуда и частота колебаний платформы. Однако изменение частоты осцилляций сигнала и3 не влияет на результат измерений, поскольку величина Дг определяется числом этих осцилляций, каждая из которых связана с изменением Дг на X 3/2.
5. Алгоритм обработки сигналов интерферометров и экспериментальные результаты
Таким образом, на выходе рассматриваемой интерференционной системы формируются три интерференционных сигнала:
(1) - сигнал совпадения фокуса лазерного пучка с контролируемой поверхностью, ~ (Дг);
(2) - сигнал начала отсчета положения сканирующей платформы, и2 (Дг);
(3) - сигнал отсчета величины смещения платформы, щ3 (Дг).
При изменении положения поверхности объекта вдоль оптической оси объектива в предметном плече сканирующего интерферометра (изменении величины расфокусировки Дг1), изменяется положение измерительного импульсного сигнала щ(Дг) сканирующего интерферометра относительно импульса м2(Дг) низкокогерентного интерферометра, играющего роль опорного интерференционного импульса, в шкале смещения сканирующей платформы Дг. Величина смещения измерительного импульса относительно опорного импульса определяется с помощью так называемой «интерференционной линейки» - сигнала и3(Дт) интерферометра с параллельным лазерным пучком. Цена деления этой линейки Дгт1П, как величина смещения платформы, при которой максимум интерференционной осцилляции заменяется минимумом, равна Дгт1П = А3/4.
Процедура определения величины смещения поверхности объекта измерительного импульса сканирующего интерферометра сводится к следующему алгоритму:
(1) определяется расстояние Дгт1 между измерительным и опорным импульсами при исходном положении объекта;
(2) определяется расстояние Дгт2 между измерительным и опорным импульсами при смещенном положении объекта;
(3) находится величина смещения (или высоты рельефа) поверхности объекта ДА = Дгт2 - Дгт1.
Величина ДА может определять продольное смещение объекта (вдоль оптической оси сканирующего интерферометра), либо определять рельеф или наклон поверхности объекта при его поперечном смещении Дх относительно оптической оси.
Погрешность определения величины ДА определяется погрешностями определения положений опорного и измерительных интерференционных импульсов в шкале интерференционных осцилляций сигнала м3 (Дг) интерферометра с параллельным лазерным пучком. Нами экспериментально установлено, что величина средней квадратической ошибки среднего значения положения измерительного импульса относительно опорного импульса для объекта с гладкой поверхностью не превышает 0.5А3.
Ожидается, что для объектов с шероховатой поверхностью эта погрешность будет несколько больше и составит порядка (1-1,5) Х3.
В наших экспериментах для удобства оцифровки и дальнейшей компьютерной обработки выходных сигналов описанной системы трех связанных интерферометров сигнал щ2(Дг) низкокогерентного
интерферометра и сигнал щ (Дг) сканирующего ла-
зерного интерферометра со сфокусированным предметным пучком электрически смешивались и формировали сигнал положения поверхности объекта. Сигнал положения поверхности объекта и интерференционная линейка (сигнала и3 (Дг) интерферометра с параллельным лазерным пучком) синхронно оцифровывались и записывались в виде числовых файлов с помощью двухканальной платы аналого-цифрового преобразования (ЛА-н150-14, ЗАО «Руднев-Шиляев») и программы - регистрирующего осциллографа, поставляемой вместе с этой встраиваемой в компьютер платой АЦП.
На рис. 4 приведены цифровые осциллограммы сигнала положения поверхности объекта (рис. 4а) и интерференционной линейки (рис. 46). В качестве объекта служило зеркало. Цифрами обозначены: измерительный импульс (1), формирующийся в сигнале щ (Дг) лазерного сканирующего интерферометра, и опорный импульс (2), формирующийся в сигнале щ2 (Дг) низкокогерентного интерферометра.
Отсчеты
Рис. 4. Цифровые осциллограммы сигнала Щ (Дг) + И2 (Дг) положения поверхности объекта (а) и интерференционной линейки М3(Дг) (б): (1) - интерференционный импульс сканирующего лазерного интерферометра (измерительный);
(2) - интерференционный импульс низкокогерентного интерферометра (опорный)
Положение опорного импульса (2) в осциллограмме могло быть изменено путем увеличения или уменьшения оптической разности длин плеч низкокогерентного интерферометра при плавном смещении опорного зеркала (изменении величины 2Дг02 в формуле (8). Положение измерительного импульса (1) в осциллограмме менялось при продольном смещении поверхности объекта вдоль оптической оси фокусирующего свет объектива (при изменении расфокусировки ДгД
В силу неравномерности движения сканирующей платформы, как сигнал положения поверхности объекта, так и сигнал интерференционной линейки являлись частотно-модулированными в полосе часто от 0 до 38 кГц. Частота оцифровки этих сигналов с помощью АЦП составляла 1 МГц.
Алгоритм определения расстояния между опорным и измерительным импульсами в шкале смещения Дг сканирующей платформы (или изменения оптической разности хода 2Дг) заключался в следующем. Брался числовой файл осциллограммы интерференционной линейки в виде столбца чисел - значений сигнала, номер у -ой строки этого столбца совпадал с номером отсчета при оцифровке данного сигнала АЦП, и находились точки перехода этого осциллирующего сигнала через ноль, вернее точки, наиболее близко расположенные к нулевому значению. Для этого все значения перебирались таким образом, что каждый раз в цикле рассматривались и сравнивались по три следующих друг за другом значения сигнала линейки (5 3 у-1, 5 3 у
и 53+1). За «точку перехода через ноль» бралось то значение сигнала 53 у , соседние с которым значения 5 3-1 и 5 3+1, имеют различные знаки. Номер строки
(номер отсчета), содержащей это значение, запоминался. Таким образом формировался столбец номеров отсчета сигнала интерференционной линейки, когда этот сигнал переходит через ноль. Соседние значения этого столбца определяют половину интерференционной осцилляции, т.е. соответствуют смещению сканирующей платформы относительно объекта на величину, равную в шкале этого смещения Дг = Х3/4 (или в шкале оптической разности хода 2Дг = Х3/2 ). В результате становится возможным каждому значению сигнала интерференционной линейки (каждому номеру строки, каждому номеру отсчета) поставить в соответствие определенное значение смещения Дг или разности хода 2Дг, зная сколько точек сигнала находится между «точками перехода через ноль» и полагая скорость движения сканирующей платформы постоянной на каждом из участков, равных Х3 /4 (на каждой половине интерференционной осцилляции). Поскольку сигнал ^(Дг) + м2(Дг) положения поверхности объекта и сигнал и~3 (Дг) интерференционной линейки оцифровываются и записываются в компьютер синхронно и с одинаковой частотой, то значениям отсчетов сигнала положения поверхности ставятся в соответствие те же значения смещения Дг или разности хода 2Дг , что и в сигнале линейки. Таким образом, далее необходимо определить положение опорного и измерительного импульсов уже в шкале смещения Дг сканирующей платформы.
Для этого, как и в случае с сигналом интерференционной линейки и~3 (Дг) , в цифровой осциллограмме сигнала положения поверхности объекта и1(Дг) + м2(Дг) ищутся «точки перехода через ноль»
по описанной выше процедуре. Отличие этого сигнала от сигнала интерференционной линейки заключается в том, что интерференционные осцилляции существуют не на всем протяжении сигнала. Это обуславливает необходимость использования интерференционной линейки для перевода количества отсчетов АЦП в сигнале положения поверхности объекта в смещение Дг или разность хода 2Дг , задаваемых движением сканирующей платформы, через длину волны Х3 лазера, используемого в интерферометре с квазипараллельным пучком. При этом длины волн света Х1 и Х2 , используемые в
сканирующем лазерном интерферометре с острос-фокусированным предметным пучком и низкокогерентном интерферометре, соответственно, могут быть не равны друг другу и длине волны Х3 . В наших экспериментах для создания интерференционной линейки мы использовали излучение полупроводникового лазера с А3 = 0,65 мкм.
После нахождения «точек перехода через ноль» в сигнале положения поверхности объекта И] (Дг) + и2 (Дг), на каждом из интервалов между соседними точками «перехода через ноль» значения сигнала брались по модулю, и находилось максимальное значение модуля сигнала, его номер строки (номер отсчета) и соответствующее этому номеру значение смещения Дг , определенное по интерференционной линейке. То есть находились локальные максимумы в сигнале, взятом по модулю, на каждом полупериоде интерференционных осцилляций исходного сигнала. Таким образом, фактически, по вершинам осцилляций снималась огибающая интерференционного сигнала положения поверхности объекта. На рис. 5 приведена огибающая сигнала и1(Дг) + ~2(Дг), показанного на рис. 4а, построенная в шкале смещения сканирующей платформы Дг . Максимумы на этом графике, обозначенные цифрами (1) и (2), соответствуют измерительному и опорному интерференционным импульсам.
Рис. 5. Огибающая интерференционного сигнала положения поверхности объекта с найденным в шкале смещения Дг сканирующей платформы расстоянием между максимумами (1) и (2),
соответствующими измерительному и опорному интерференционным импульсам
Далее по графику этой огибающей находилось расстояние в шкале смещения Дг между максимумами, соответствующими измерительному и опор-
ному интерференционным импульсам. Для этого в матрице значений огибающей находилось максимальное значение среди ее элементов, его номер и соответствующее смещение Дг . Затем исходная матрица нормировалась на это максимальное значение и шкала значений смещения Дг сканирующей платформы смещалась таким образом, чтобы максимальному значению нормированной огибающей соответствовало нулевое значение смещения Дг (рис. 5). Далее осуществлялся поиск значения Дгт положения второго пика огибающей сигнала в шкале смещения Дг сканирующей платформы.
Этот поиск осуществлялся следующим образом. Программно, по процедуре, сходной с нахождением «точек перехода через ноль», находились значения смещения Дг и соответствующие им номера элементов матрицы значений огибающей, нормированной описанным выше способом, соответствующие значению 0,5 от максимального (равно 1). Таких значений получается два (8г1 < 0 и 5г2 > 0 ,
и ) - в обе стороны от максимального значения в шкале смещения Дг. Величина |5г1 + |8г2| приблизительно равняется 1с или Ьс, в зависимости от того, какой из интерференционных импульсов, опорный или измерительный был больше по амплитуде. Вообще, перед началом измерительных экспериментов опорный и измерительный импульсы разводятся в осциллограмме на расстояние, много превышающее сумму 1С + Ьс, во избежание наложения этих импульсов друг на друга. Добиться этого можно изменением длины одного из плеч низкокогерентного интерферометра (изменяя 2Дг02 ).
Далее в матрице значений нормированной огибающей находились номера элементов, для которых смещение Дг приблизительно равно 2Ъгг и 25г2. Этого достаточно, чтобы значения огибающей в этих точках существенно упали по сравнению с максимальным. После этого производился поиск максимальных значений нормированной огибающей среди элементов с номерами от одного, до номера, которому соответствует смещение и 2Ъгг, и среди элементов с номерами, начина с номера, соответствующего смещению и 25г2 до максимального номера. При этом фиксировались значения смещения Дг (Дг ' < 0 , Дг " > 0), соответствующие этим двум значениям нормированной огибающей. Затем эти значения нормированной огибающей сравнивались, выбиралось наибольшее, и величина абсолютного значения смещения |Дг '| или |Дг"|, соответствующего большему значению огибающей, принималась за искомое расстояние Дгт между опорным и измерительным максимумами.
По описанному выше алгоритму найдено расстояние между максимумами огибающей в шкале смещения сканирующей платформы, график которой приве-
ден на рис.5. Найденное положение второго максиму -ма отмечено на этом рисунке штриховой меткой.
В качестве примера измерительных возможностей системы приведены результаты измерения расстояния Дгт между опорным и измерительным импульсами в сигнале Щ1(Дг) + и2(Дг) в шкале смещения сканирующей платформы при смещении Дх зеркального объекта в направлении, поперечном оптической оси предметного пучка (рис. 6). Поперечное смещение Дх осуществлялось с помощью микрометрического столика, цена деления шкалы микрометрического винта которого составляла 10 мкм.
Кружками на рис.6 обозначены значения Дгт, усредненные по 10 измерениям. Вертикальные отрезки показывают доверительный интервал полученных средних значений Дгт (для доверительной вероятности Р = 0,99 и числа измерений к = 10 по распределению Стьюдента). Пунктирная прямая является результатом линейной регрессии полученных экспериментальных данных. Наклон этой прямой позволяет определить угол наклона поверхности плоского зеркала, служившего объектом, по отношению к плоскости, перпендикулярной оптической оси фокусирующего свет объектива в предметном плече сканирующего интерферометра. Этот угол наклона, определенный по графику на рис. 6 составляет 0,3° (18 угловых минут).
?0 100 150 200 2?0 J00
Дх, мкм
Рис. 6. Результаты измерений расстояния bzm между опорным и измерительным максимумами в сигнале положения поверхности объекта при поперечном смещении Дг плоского зеркального объекта
Рассматриваемая система, как отмечалось выше, может быть использована для контроля толщины слоев слоистых объектов [27]. В этом случае сигнал сканирующего интерферометра с остросфокусированным пучком будет содержать последовательность интерференционных импульсов, положение которых в шкале разности хода будет определяться толщинами слоев. Однако, в отличие от низкокогерентной интерферометрии, положение этих импульсов определяется не оптической толщиной слоев 2njhj (nj - показатель преломления слоя, hj - геометрическая толщина слоя), а отношением 2hjjnj .
Заключение
В рамках настоящей работы показано, что лазерный интерферометр с остросфокусированным зондирующим пучком и интегрирующей апертурой фотоприемника наравне с конфокальной микроскопией и низкокогерентной интерферометрией может служить надежным средством прецизионного контроля формы рельефа поверхности и параметров структуры слоистых объектов. Основные достоинства данного интерферометра заключаются: в возможности использования достаточно мощных лазерных источников и, соответственно, в получении сигнала достаточной величины при контроле сильно рассеивающих сред, когда имеют место большие энергетические потери; в существенно более низких требованиях к равенству оптических путей в плечах интерферометра; в более высокой поперечной и продольной разрешающих способностях. Немаловажным обстоятельством является более выгодное соотношение стоимостей источников излучения - суперлюминесцентные светодиоды, используемые в низкокогерентной интерферометрии, существенно дороже полупроводниковых лазерных модулей. Кроме этого, применение лазерных модулей коротковолнового излучения обеспечивает более высокую точность измерений.
Рассмотренная в настоящей работе интерференционная система достаточно надежна в эксплуатации вне лабораторных условий. Эта система может быть использована для размерного контроля формы и слоистой структуры объектов технического происхождения, а также в диагностических целях в биомедицине.
Работа выполнена при поддержке грантов РФФИ № 05-08-50318-а, CRDF № Y2-BEP-06-17 (REC-006, ANNEX BF4M06) и CRDF RUX0-006-SR-06 (ANNEX BP1M06).
Литература
1. Коломийцов Ю.В. Интерферометры. Основы инженерной теории. Применения // Л.: Машиностроение, 1976.
2. Wyant J.C. Interferometric optic metrology: basic principles and new systems // Laser Focus, 1982. V.18. P. 65-71.
3. Батраков А.С., Бутусов М.М., Гречка Г.П. и др. Лазерные измерительные системы // М.: Радио и связь, 1981.
4. Дубнищев Ю.Н., Ринкевичус Б.С. Методы лазерной доплеровской анемометрии // М.: Наука, 1982.
5. Коронкевич В.П., Соболев В.С., Дубинцев Ю.Н. Лазерная интерферометрия // Новосибирск: Наука, 1983.
6. Аблеков В.К., Колядин С.А., Фролов А.В. Высокоразрешающие оптические системы // М.: Машиностроение, 1985.
7. Чудов В. А. Измерения в машиностроении // Измерительная техника, 1990. №3. С.61-62.
8. Васильев В.Н., Гуров И.П. Оптическая интерферометрия и информационные технологии. Системы преци-
зионного технологического контроля // Известия вузов. Приборостроение, 1996. №5-6. С.13-20.
9. Васильев В.Н., Гуров И.П. Технология бесконтактного контроля объектов на основе когерентного и спектрального радаров в биомедицинских исследованиях и промышленности, в кн.: Оптические и лазерные технологии // СПб.: СПбГУ ИТМО, 2002.
10. Гуров И.П., Джабиев А.Н. Интерферометрические системы дистанционного контроля объектов // СПб.: СПбГУ ИТМО, 2000.
11. Wilson T., Sheppard C.J.R. Theory and Practice of Scanning Optical Microscopy // London: Academic Press, 1984.
12. Wilson T. ed. Confocal Microscopy // London: Academic Press, 1990.
13. Ченцов Ю.В. Конфокальная лазерная сканирующая микроскопия в биологии и медицине // Оптический журнал, 1994. №12. С.18-23.
14. Лужнев Э.И., Попова И.П., Кузьмин С.В., Сла-щев С. М. Конфокальная лазерная сканирующая микроскопия: принципы, устройство, применение (Часть 1) // Научное приборостроение, 2001. Т.11. №2. С.3-20.
15. Hamilton D.K., Sheppard C.J.R. A Confocal Interference Microscope // Optica Acta, 1982. V. 29. N. 12. P. 15731577.
16. Matthews H.J., Hamilton D.H., Sheppard C.J.R. Surface profiling by phase-locked interferometry // Applied Optics, 1986. V. 25. N. 14. P. 2372-2374.
17. Gu M. Time-resolved three-dimensional imaging based on confocal interferometry under ultrashot pulsed illumination // Optik, 1996. V. 104. N. 1. P. 32-34.
18. Fercher A.F. Optical coherence tomography // J. Biomed. Opt., 1996. V. 1. P. 157-173.
19. Schmitt J.M. Optical coherence tomography: a review // IEEE J. Select Topics Quant. Electron., 1999. V. 5. P. 1205-1215.
20. Masters B.R. ed. Optical low-coherence reflectometry and tomography // SPIE Milestone Series, 2001. V. MS165.
21. Fercher A. F., Drexler W., Hitzenberger C. K., Lasser T. Optical coherence tomography - principles and applications // Report Prog. Phys., 2003. N. 66. P. 239-303.
22. Bouma B.E. and Tearney G.J. eds. Handbook of optical coherence tomography // NY: Marcel Dekker Inc. 2002.
23. Гуров И.П. Оптическая когерентная томография: принципы, проблемы и перспективы // Проблемы когерентной и нелинейной оптики. Под ред. Гурова И. П. и Козлова С.А., СПб.: СПбГУ ИТМО. 2004. С. 6-30.
24. Рябухо В.П., Хомутов В.Л., Лякин Д.В., Константинов К.В. Лазерный интерферометр с остросфокусиро-ванными пучками для контроля пространственного положения объекта // Письма в ЖТФ, 1998. Т. 24. В. 4. С. 19-24.
25. Kempe M., Rudolph W. Scanning microscopy through thick layers based on linear correlation // Opt. Lett., 1994 V. 19. N. 23. P. 1919-1921.
26. Борн М., Вольф Э. Основы оптики // М.: Наука, 1973.
27. Лякин Д.В., Лобачев М.И., Резчиков А.Ф., Рябу-хо В.П., Миненкова И.Ф. Лазерная сканирующая интерференционная система для определения толщины слоистых микроструктур // Мехатроника, автоматизация, управление, 2003. №4. С.10-14.