Например, Бобцов

Оценка основных параметров волноводного микрорезонаторного преобразователя ИК изображений.

ФИЗИЧЕСКАЯ ОПТИКА
УДК 681.783.323: (621.397.61.029.71/73 : 681.772.7)
ОЦЕНКА ОСНОВНЫХ ПАРАМЕТРОВ ВОЛНОВОДНОГО МИКРОРЕЗОНАТОРНОГО ПРЕОБРАЗОВАТЕЛЯ ИК ИЗОБРАЖЕНИЙ
© 2009 г. В. А. Пилипович, доктор физ.-мат. наук; А. К. Есман, доктор физ.-мат. наук; В. К. Кулешов, канд. техн. наук; Г. Л. Зыков, канд. физ.-мат. наук Институт физики имени Б.И. Степанова Национальной академии наук Беларуси, г. Минск, Белоруссия
E-mail: lomoi@inel.bas-net.by
Показано, что использование нового метода преобразования изображений инфракрасного (ИК) диапазона в видимое изображение на основе волноводного кремниевого микрорезонатора при построении матричных ИК преобразователей тепловизионной техники, воспринимающих контраст температур порядка 10–3 K, позволяет получить пространственное разрешение 25 мкм при радиусе замкнутого оптического кремниевого волновода 10 мкм и энергии падающего на него ИК излучения 5,6×10–10 Дж с частотой полнокадрового режима 50 Гц.
Коды OCIS: 110.0110, 110.3080, 3000.33080
Поступила в редакцию 13.10.2008.

В работе [1] рассмотрены структурная схема и принцип работы матричного преобразователя изображений инфракрасного диапазона в видимый или ближний ИК диапазоны спектра электромагнитных волн с использованием в качестве чувствительного элемента микрорезонатора в виде замкнутого оптического волновода на диэлектрической подложке. Зависимость резонансной длины волны микрорезонатора от его температуры положена в основу работы преобразователя.
Данная работа посвящена оценке основных параметров микрорезонатора преобразователя инфракрасного электромагнитного излучения, физической основой работы которого являются не фотонные эффекты, а интегральные изменения оптических свойств материала чувствительного элемента, происходящие в результате изменений его температуры под действием падающего излучения.
Анализ работы волноводного микрорезонаторного преобразователя ИК изображений может быть выполнен на основе его структурной схемы, приведенной на рис. 1. Замкнутый оптический волновод 3 расположен на кольцевом выступе 5 диэлектрической подложки 6 и выполнен из кремния. Для получения условий полного внутреннего отражения сверху и снизу волновода 3 расположены буферные слои 2 и 4 из материала

PИвхК 0 y 1 hэ
2 xz
3 4
5

rэ 2 3 4
5

d2 d3 d4
d5

d

6

Pвых

8 rэ 3

Pвсхч 7

Рис. 1. Структурная схема волноводного микрорезонаторного преобразователя ИК
изображений. 1 – поглощающая пленка, 2 и 4 – буферные слои, 3 – замкнутый оптический волновод, 5 – кольцевой выступ подложки, 6 – подложка.

“Оптический журнал”, 76, 5, 2009

3

с меньшим показателем преломления – окиси кремния. На верхнем буферном слое 2 нанесена пленка 1 из материала, эффективно поглощающего тепловое излучение, например, платиновая или золотая чернь, обычно используемые в болометрах [2]. Для считывания информации на входной волновод 7 подается излучение с мощностью Pвсхч видимого или ближнего ИК диапазонов. В рассматриваемом случае – это излучение на длине волны 1,5 мкм, которое через спадающие поля поступает через замкнутый оптический волновод 3 в выходной волновод 8.
Количество энергии электромагнитного излучения, поглощенной чувствительным элементом для изменения его температуры, зависит от следующих основных параметров:
плотности потока мощности ИК излучения PИвхК, падающего на поглощающую пленку 1,
параметров поглощающей пленки 1, эффективности теплоотдачи от поглощающей пленки 1 к замкнутому оптическому волноводу 3, теплофизических параметров материала замкнутого оптического волновода 3. Эффективность преобразования падающего ИК излучения PИвхК и другие параметры преобразователя получены численным моделированием рассматриваемой структуры чувствительного элемента с учетом ее геометрических размеров и всех вышеуказанных факторов. В общем виде температурное поле в многослойной тепловоспринимающей структуре определяется из нестационарного уравнения теплопроводности [3]

ρc

¶T ¶t

=

¶ ¶x

çççèæk

¶T ¶x

÷÷÷øö

+

¶ ¶y

æççèçk

¶T ¶y

÷ø÷ö÷÷

+

+

¶ ¶z

ççèæçk

¶T ¶z

÷÷÷øö

+

Q(x,y,z,t),

(1)

где ρ – плотность материала, с – его теплоемкость, T – температура, t – время, k – теплопроводность, x, y и z – пространственные переменные, Q(x, y, z, t) – функция теплового источника.
В данном случае инфракрасное излучение направлено по нормали к поверхности поглощающей пленки 1, выполненной в виде диска с радиусом rэ, равным радиусам элементов 2–5 (рис. 1). Температура нижней поверхности подложки 6 принималась постоянной и равной начальной температуре T0.
Анализ характеристик отдельного элемента матричного преобразователя ИК изображения проводился методом численного моделирования на основе одномерного уравнения теплопроводности

ρici

¶T ¶t

=

¶ ¶x

çæèççki

¶T ¶x

ö÷÷ø÷ +

Q(x,t),

i

=

1, 2…N

при граничных и начальном условиях

(2)

¶T ¶x

= 0,
x=0

T(x = d,t) = T0,

T(x,t = 0) = T0,

(3)

где i – номер текущего слоя чувствительного

элемента матричного преобразователя изображе-

ний, d – толщина рассматриваемого чувствитель-

ного элемента, T0 – начальная температура. Так как температура в плоскости уoz (рис. 1)

предполагалась одинаковой, то при моделиро-

вании температурных полей рассматривался

одномерный случай (2). Кроме этого, темпера-

турное поле в преобразователе ИК изображений

удовлетворяет условиям сопряжения на грани-

цах раздела ξ его элементов j, выполненных из

различных материалов

T(j) = T(j + 1) при x = ξ, j = 1, 2…5,

(4)

k(j)

¶T(j) ¶x

- k(j+1)
x=ξ

¶T(j+1) ¶x

= 0,
x=ξ

(5)

где j = 1 – поглощающая пленка, j = 2 – верхний буферный слой, j = 3 – замкнутый оптический волновод, j = 4 – нижний буферный слой, j = 5 – кольцевой выступ подложки, j = 6 – подложка, T(j) и T(j+1) – температуры контактирующих поверхностей j-го и j + 1-го элементов, выполненных из различных материалов, k(j) и k(j+1) – теплопроводности рассматриваемых элементов.
Количество тепла Q(x, t), выделяемое в i-ом слое j-го элемента при поглощении ИК излучения, задавалось выражением

òQ(x,

t)

=

(1-

R)

PâÈõÊ Si

αi (x)expçççæççèçç-

x 0

αi (x¢)dx¢øö÷÷÷÷÷÷÷,

(6)

где R – коэффициент оптического отражения поглощающей пленки, αi(x) – коэффициент поглощения i-го слоя j-го элемента, PИвхК – мощность падающего ИК излучения, Si – площадь i-го слоя рассматриваемого элемента.
Уравнение (2) решалось методом конечных
разностей. Суть метода заключается в покрытии расчетной области (x, t) сеткой из N×M точек. Тем самым определяются узлы, в которых будет
осуществляться поиск решения. Шаги по вре-
мени и координате выбирались из требования получения точности счета, равной 1,5%.
Численное моделирование проводилось при следующих параметрах элементов j: толщина поглощающей пленки 1 равна одному шагу по координате x, а теплоемкость единицы ее площади –

4 “Оптический журнал”, 76, 5, 2009

8,368×10-3 Дж см–2 K–1 [2], ширина верхнего 2 и нижнего 4 буферных слоев из SiO2 и замкнутого оптического волновода 3 из Si одинаковы и равны hэ = 0,5 мкм, а их толщина составляет d2 = d3 = = d4 = 0,5мкм. Толщина кольцевого выступа подложки 5, выполненного из кремния, составляет
d5 = 20 мкм. Теплофизические и оптические параметры материалов чувствительного элемента,
использованные при расчетах, представлены в
таблице.
При расчетах выбирались размеры чувстви-
тельного элемента 10÷40 мкм (rэ = 5÷20 мкм) с учетом теоретических пределов дифракционного
рассеяния для дальнего ИК диапазона [4].
Считывание информации с чувствительного
элемента осуществляется при достижении мак-
симальных изменений температуры в замкнутом
оптическом волноводе 3. Зависимости значений этих изменений температуры замкнутого оптического волновода 3 от входной мощности PâÈõÊ для различных длительностей воздействия ИК
излучения приведены на рис. 2. Из характери-
стик, полученных для радиуса замкнутого опти-
ческого волновода 3 rэ = 10 мкм, следует, что при значении его температурного контраста порядка 5×10-4 K энергия на один чувствительный элемент составляет 3,1×10–11, 1,6×10–10 и 1,5×10–11 Дж при временах воздействия ИК излучения соот-
ветственно 1, 20 и 30 мс или при входных мощностях 3,1×10–8, 8×10–9 и 5×10–10 Вт/элемент. В разрабатываемых современных преобразова-
телях ИК изображений в видимое изображение
значение температурного контраста составляет порядка 10–3 K [5]. Таким образом, для достижения указанного значения необходимо на вход рас-
сматриваемого преобразователя подать мощности порядка 8,4×10–8, 2,8×10–8 и 1,8×10–8 Вт/элемент при длительностях воздействия соответственно
1, 20 и 30 мс.
Время реагирования чувствительного эле-
мента на изменение ИК сигнала, т.е. временные
зависимости изменения температуры замкнутого
оптического волновода 3 с различными радиусами rэ при минимальном времени воздействия – 1 мс и плотности входной мощности ИК излучения ~3,2×10–2 Вт/см2, падающей на него, при-

Параметры ρ, г/см3
c, Дж/(г K) k, Вт/(см K)
α, см–1

Si 2,33 1,648 0,4 418,9

SiO2 2,65 0,89 0,0155 0,16755

T, 10–3 °C

ведены на рис. 3. Указанная плотность мощно-
сти соответствует уровню теплового воздействия 10–7 Вт/элемент при rэ = 10мкм. Как видно из этих зависимостей, сумма времен срабатывания и вос-
становления после снятия теплового воздействия
для радиусов замкнутого оптического волновода
3, равных 5, 10 и 20 мкм, соответственно составляют примерно 14, 20 и 35 мс.
Таким образом, при радиусе замкнутого
оптического волновода 3 rэ = 10 мкм возможно построение рассматриваемого матричного ИК

Tmax, 10–3 °C

3
3 2
2
11

0 10–9

10–8

10–7

PИвхК, Вт/элемент

Рис. 2. Зависимости максимальных изменений температуры замкнутого оптического
волновода от мощности PИвхК для разных длительностей воздействия. 1 – 1 мс, 2 – 20 мс, 3 – 30 мс.

1,2
1,0
0,9 3
0,6 2
0,3 1
0 0 30 60 t, мс
Рис. 3. Зависимости изменения температуры замкнутого оптического волновода ΔT от времени t. 1 – rэ = 5 мкм, 2 – rэ = 10 мкм, 3 – rэ = = 20 мкм.

“Оптический журнал”, 76, 5, 2009

5

преобразователя с частотой полнокадровой развертки 50 Гц.
Эффективность преобразования падающего ИК излучения PИвхК с учетом того, что золотая чернь во всем диапазоне длин волн поглощает 90% энергии [2], оценивалась по отношению максимальных температур замкнутого оптического волновода 3 и поглощающей пленки 1. Для кривых 1–3, приведенных на рис. 3, такая эффективность преобразования соответственно составила 0,32; 0,35; 0,404, т. е. с увеличением радиуса эффективность преобразования увеличивается.
Зависимость максимального изменения температуры ΔT замкнутого оптического волновода 3 от радиуса rэ при плотности мощности входного ИК излучения порядка 3,2×10–2 Вт/см2 приведена на рис. 4. С увеличением времени теплового воздействия и радиуса замкнутого оптического волновода изменение температуры ΔT возрастает.
На рис. 5 представлены распределения температурных изменений в слоях элементов j = 2, 3, 4, 5 преобразователя в различные моменты времени, из которых следует, что на стадии нагревания замкнутого оптического волновода 3, имеющего радиус 10 мкм, градиент температуры по его глубине незначителен, и, следовательно, влиянием пьезоэффекта на оптические свойства чувствительного элемента можно пренебречь.
Для получения номинального значения оптического сигнала на выходе чувствительного элемента на основе кремния с радиусом замкнутого оптического волновода 10 мкм необходимо время 100 пс [6]. Поэтому при времени считывания информации порядка 200 пс и мощности входного считывающего сигнала Pвсхч = 1 мВт в замкнутом оптическом волноводе 3 выделится примерно 10–13 Дж тепловой энергии. Это значение получено с учетом того, что экспериментально измеренное затухание оптического сигнала в кремниевых резонансных волноводных структурах составляет от 2,5 до 5,5 дБ [7]. При использовании мощности считывания Pвсхч = 10–7 Вт, что вполне реально для современных фотоприемников, тепловыделение в замкнутом оптическом волноводе снижается на 4 порядка. Следовательно, тепловыделением при считывании можно пренебречь.
Математическое моделирование основных параметров чувствительного элемента преобразователя ИК изображений на основе замкнутых волноводных кремниевых микрорезонаторов показывает, что использование нового метода

T, 10–3 °C

Tmax, 10–3 °C

3,3 3 2
3,0
2,7 1,4
1 1,2

1,0 5

10 15 rэ, мкм

Рис. 4. Зависимости максимального изменения температуры замкнутого оптического волновода от его радиуса для разных длительностей воздействия ИК излучения. 1 – 1 мс, 2 – 20 мс, 3 – 30 мс.

1,20 1,05

2 мс

0,90 1 мс

3 мс 4 мс

0,75
2 34 0,60
0 0,3 0,6 0,9

5 1,2 x, мкм

Рис. 5. Профиль изменения температуры в отдельном чувствительном элементе преоб-
разователя ИК изображения на глубине до
1,5 мкм с радиусом чувствительного элемента
rэ = 10мкм в различные моменты времени на стадии нагревания. Области 2 – верхний буферный слой, 3 – замкнутый оптический волновод, 4 – нижний буферный слой, 5 – кольцевой выступ подложки.

преобразования изображений инфракрасного диапазона в видимое изображение в современных тепловизионных приборах, воспринимающих контраст температур порядка 10–3 K, позволяет достичь следующих основных параметров:
пространственное разрешение – 25 мкм, частота полнокадрового режима – 50 Гц, время считывания одной строки – 0,2 нс, энергия ИК излучения, падающего на один элемент – 5,6×10–10 Дж,

6 “Оптический журнал”, 76, 5, 2009

радиус замкнутого оптического волновода – 10 мкм.
Таким образом, на основе предложенного чувствительного элемента возможно построение нового типа матричного преобразователя ИК изображений, удовлетворяющего современным требованиям тепловизионной техники и неподверженного электромагнитным наводкам.
ЛИТЕРАТУРА
11. Пилипович В.А., Есман А.К., Гончаренко И.А., Кулешов В.К. Матричный преобразователь инфракрасных изображений в видимые на основе волноводных микрорезонаторов // Оптический журнал. 2008. Т. 75. № 5. С. 59–64.
12. Синцов В.Н. Исследование свойств золотой черни // ЖПС. 1966. Т. 4, вып. 6. С. 503–508.
13. Мельников А.А. Расчет температурных полей в многослойных фотоприемных структурах //

Нано- и микросистемная техника. 2000. № 2. С. 21–26.
14. Певцов Е., Чернокожин В. Матричные ИК приемники для малогабаритных тепловизионных камер // Электронные компоненты. 2001. № 3. С. 12–20.
15. Таубкин И.И., Тришенков М.А. Предельные значения чувствительности тепловизоров и информативности систем технического зрения // Прикладная физика. 2001. № 6. С. 48–60.
16. Пилипович В.А., Есман А.К., Гончаренко И.А., Кулешов В.К., Зыков Г.Л. Управляемый интегрально-оптический фильтр высокого разрешения на волноводных микрорезонаторах // Доклады НАН Беларуси. 2008. Т. 52. № 3. C. 48–51.
17. Barwicz T., Popovic M.A., Watts M.R., Rakich P.T., Ippen E.P., Smith H.I. Fabrication of add-drop filters based on frequency-matched microring resonators // Journal of Lightwave Technology. 2006. V. 24. № 5. P. 2207–2218.

“Оптический журнал”, 76, 5, 2009

7