CC BY
25
НАНОСИСТЕМЫ: ФИЗИКА, ХИМИЯ, МАТЕМАТИКА, 2012, 3 (4), С. 20-27 УДК 535.361.1
ОПРЕДЕЛЕНИЕ ВЕЛИЧИНЫ НАНОРАЗМЕРНЫХ НЕОДНОРОДНОСТЕЙ В РАССЕИВАЮЩИХ И ПОГЛОЩАЮЩИХ СРЕДАХ
Д.Н. Вавулин, А. В. Альфимов, О. В. Андреева, Е. М. Арысланова, А. В. Пантелеев, С. А. Чивилихин
Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, Санкт-Петербург, Россия
dima-vavulin@mail.ru, alfimov.anton@gmail.com, olga_andreeva@mail.ru, elizabeth.aryslanova@gmail.com, pantandrej@yandex.ru, sergey.chivilikhin@gmail.com
PACS 45.25.Dd
В работе предложен оригинальный метод определения размеров наноразмерных неоднородностей в рассеивающей и поглощающей среде путём сравнения теоретической и экспериментальной зависимостей коэффициента пропускания среды от длины волны света. Метод был опробован на образце нанопористого стекла и дал удовлетворительные результаты. Также были проведены численные оценки явления с помощью программного пакета CST Microwave Studio.
Ключевые слова: Рэлеевское рассеяние, наноразмерные неоднородности, нанопористые стёкла, диффузия излучения.
1. Введение
Определение характерного размера неоднородностей среды является важной научно-технической задачей, которая привлекает внимание многих исследователей [1,3,4,11]. Известные на данный момент способы решения этой задачи зачастую являются непригодными ввиду своей сложности или высокой стоимости (например, малоугловое рассеяние рентгеновских лучей). Нами был предложен простой и дешёвый метод определения размеров наноразмерных неоднородностей в рассеивающих и поглощающих средах. Метод был опробован на примере пластины из поглощающего нанопористого стекла, в котором неоднородности (поры со средним диаметром 17 нм, заполненные водой) служили рассеи-вателями.
Для реализации данного метода была получена теоретическая зависимость коэффициента пропускания от коэффициентов рассеяния и поглощения и экспериментальная зависимость коэффициента пропускания от длины волны света.
Основная идея метода заключается в следующем: при больших длинах волн поглощение среды становится пренебрежимо малым и ослабление света происходит за счёт рассеяния на неоднородностях среды. Таким образом, возможно определить экспериментально коэффициент рассеяния, исследуя пропускание пластины в диапазоне больших длин волн.
Так как характерный размер неоднородностей среды много меньше длины световой волны, рассеяние имеет рэлеевский характер. Следовательно, коэффициент рассеяния
обратно пропорционален четвёртой степени длины волны. Далее, коэффициент пропорциональности в выражении для коэффициента рассеяния и спектральная зависимость коэффициента поглощения подбираются таким образом, чтобы обеспечить наилучшее согласование между теоретической и экспериментальной зависимостями коэффициента пропускания от длины волны. Так как рассеяние имеет рэлеевский характер, полученный коэффициент пропорциональности зависит от шестой степени диаметра пор. Из полученных результатов определяется средний размер пор.
2. Вывод теоретической зависимости коэффициента пропускания от параметров системы
2.1. Вывод уравнения диффузии излучения
Для получения интересующей нас зависимости коэффициента пропускания от коэффициентов рассеяния и поглощения, необходимо описать процесс распространения света в рассеивающей и поглощающей среде. Выведем сначала для этого уравнение диффузии излучения.
В предположении, что размер рассеивающих частиц много меньше длины волны света [4], уравнение переноса излучения имеет вид:
ОУ/ = - (та + ъ) I + ^Ш, (1)
где 1 — интенсивность излучения, О — единичный вектор в направлении луча, 7а и 7« — коэффициенты поглощения и рассеяния среды (которые считаются постоянными), ¿ш — элемент телесного угла. Величина (^а + 7«)-1 представляет собой характерный пространственный масштаб изменения интенсивности излучения. Будем считать, что эта величина мала по сравнению с характерными размерами рассматриваемой системы. Введем суммарную интенсивность излучения, падающего на данный элемент среды:
( = ! Ыш (2)
и применим к (1) преобразование Фурье. Тогда имеем:
¿Ш/к = - (7а + 7з) 1к + у1 (к,
4п
где 1к = / 1е гкх^3ж. Откуда
4 = — (к . (3)
7а + 7« + ¿кО'
Интегрируя (3) по полному телесному углу и используя (2), получаем:
7« ^ [ ¿ш
(к = -¡— (к
или
4п У та + 7« + ¿кО '
7« , ( к
1 — у1 аг^ (
+ , (к = 0, (4)
к \ 7а + 7«
где к = |к|. Выражение (4) представляет собой Фурье-представление псевдодифферен-
к
циального уравнения [7] для величины (. В длинноволновом приближении - ^ 1
7а + 7«
это уравнение приобретает вид:
к2Ъ
я = о.
(5)
37а (1а + 7«)2
Ограничимся рассмотрением сред, в которых коэффициент рассеяния 7з существенно превышает коэффициент поглощения 7а. Тогда, переводя (5) в координатное представление, т.е. производя обратное преобразование Фурье, получаем искомое уравнение диффузии излучения [4]:
[А - 37а7*] Я = 0. (6)
В тех же предположениях может быть получено граничное условие к уравнению (6):
о - ^ = 2Я
27з дп
(7)
д
где —--производная вдоль внешней нормали к границе области, Яе — суммарная интен-
дп
сивность излучения, падающая снаружи на границу области.
2.2. Вывод коэффициента пропускания нанопористого стекла
Теперь получим коэффициент пропускания среды для случая плоскопараллельной пластины нанопористого стекла толщины к. Тогда мы можем рассматривать изменение интенсивности света только в одном пространственном направлении (вдоль оси х).
Согласно (6), (7), уравнение диффузии излучения и граничные условия к нему имеют вид:
где
й2Я
йх2
- 3^а!зЯ = 0,
Я - — йЯ
йх
2Яе
х=0
Я + 27^ йх
0.
х=Н
Интегрируя (8) с использованием (9), получаем:
А =„ 8Ь(а) + всЬ (а) , , 2Яе, В
Я = А сЬ (ах) + В вЬ (ах),
сЬ (а) + в вЬ (а)
(1 + в2)вЬ(а) + 2в сЬ (а) (1 + в2) вЬ (а) + 2в сЬ (а)
а = к^/3тОтЗ, в = 0, 5а/7з.
2Яе
(8) (9)
(10)
Зная величину Я, мы можем определить интенсивность излучения, распространяющегося под углом
1 = / - с08 (I) йО 4п \ 7з йх
(11)
и суммарную интенсивность излучения, распространяющегося в положительном направлении вдоль оси х:
Г п/2
Я+ = 2И I яп(^) ио
Определение величины наноразмерных неоднородностей 23
(+ (Л)
Используя (10), (11), (12), рассчитаем искомый коэффициент пропускания Т пластины нанопористого стекла:
(+ (0)
Т =----^- . (13)
еЬ (Л^37а7«) + У (ЛУЗтОт!)
Из полученного соотношения в совокупности с экспериментальными результатами предложенный нами метод позволяет определить характерный размер пор.
3. Проведение эксперимента
Нанопористые матрицы на основе силикатного стекла представляют собой уникальный инструмент исследования физико-химических процессов в ограниченном объеме, соизмеримом с масштабом протекания процессов и размерами изучаемых объектов. Ограниченность пространства и эффективный контакт со стенками пор обуславливают существенные особенности состояния и свойств заполняющего вещества по сравнению со случаем его нахождения в свободном объеме.
Применение нанопористых стекол в настоящее время, главным образом, связано с прозрачностью таких стекол в видимой области спектра и возможностью получения образцов оптического качества [8,9]. Оптическая плотность образцов в ближней УФ области спектра обусловлена рассеянием на структуре и поглощением материала. При рассмотрении оптических свойств пористых стекол в видимой области спектра и ближнем ИК диапазоне поглощением в образцах можно пренебречь ввиду малости этой величины, а эффективные оптические постоянные образцов определяются только ослаблением излучения за счет рассеяния на пористой структуре и ее неоднородностях. Это свойство нанопористых стёкол и легло в основу данного метода.
Для экспериментального измерения зависимости коэффициента пропускания образца рассеивающей и поглощающей среды, в качестве модельного объекта был взят полированный диск нанопористого стекла НПС-17 диаметром 15 мм и толщиной 1 мм, со средним размером пор 17 нм и относительным объемом пор 58%.
Перед проведением измерений поры были заполнены водой. Для этого образец был погружен на сутки в дистиллированную воду. Измерения проводились на стандартном спектрофотометре Еуо1и1юп-300. Далее измерялся коэффициент пропускания образца на различных длинах волн в диапазоне 350-1100 нм. Результаты измерений приведены на Рис. 1.
4. Результаты и обсуждение
Полученные экспериментальные и расчётные зависимости коэффициентов пропускания и поглощения приведены на Рис. 1.
Коэффициент рассеяния обратно пропорционален четвертой степени длины волны
С
падающего света, 7« = —. Коэффициент С и спектральная зависимость коэффициента
поглощения 7а подобраны так, чтобы обеспечить наилучшее согласование между теоретической и экспериментальной зависимостями коэффициента пропускания от длины волны света. Полученное при этом значение коэффициента С позволяет приближенно оценить характерный размер неоднородности среды В ~ С1/6.
Расчетное значение среднего диаметра пор В = 13,5 нм близко к среднему размеру пор (17 нм) нанопористого стекла. Это говорит о том, что предложенный метод обладает сравнительно высокой точностью при своей простоте и наглядности.
Рис. 1. Экспериментальная и расчетная зависимости коэффициента пропускания (а), а также коэффициента поглощения нанопористого стекла (б), от длины волны света
5. Вывод коэффициента пропускания нанопористого стекла (второй порядок точности)
Во втором порядке точности уравнение диффузии излучения будет иметь, соответственно, вид:
5
А2 + ^78 А - 57а7,
Я = 0.
(14)
Граничные условия к нему:
Я + ^ дЯ + ^
275 дп 37^ дп2
= 2Я
+ АдО + п д2я
4 дп 167? дп2
= 2и
(15)
(16)
где суммарная интенсивность и некоторая мера расходимости пучка равны соответственно:
п/2
Qe = 2п J I sin(tf) dtf, (17)
0
п/2
Ue = 2п j I sin2 (tf) dtf. (18)
0
Выражение для коэффициента пропускания T = Q+ ( ) в данном случае в неявном
Q+ (0)
виде выглядит следующим образом:
T = Q (h) = (aisi + s2) c3 (19)
Qe Qe
где:
5
s1 = sh(Aih), S2 = sh(A2h), t1 = ch(A1h), ¿2 = ch(A2h),
a2,2 = 1 (-«± va^+ib)
a
c =
d
3
5ТаТ3, 1
3Y2,
а1
Ai
A3;
«2 =
A
Ал
аз =
Аз
A4
и имеют место соотношения:
3
п
167s2
п
m = —,
n =
х
37s'
(20)
(21)
3f
2
3
2
3
sJ 1 - -m ) + s1A1 d —— c2 + sW 1 - -m + s2A2 d —— c4 = 2Qe — 3Ue
2
3/
2
s^ 1 — 2mj + S1A2 ^—d — 31 j j C1 + ^2 ^ 1 — 2m j + S2A2 ^—d — 31 ) ) C3 —nA1c2 + fc1fc2c3 — nA2c4 = 0,
c3 c3
—nA1«1C3 + кза2--ПА2С3 + k1— = 2Ue
аз «3
(22)
где:
k2 = -
k1 = ms2 + nA2¿2 + /А^2, ms1 + nA111 + /A2s1
s1 (1 — §m) + s1A2 (—d — §/)
k3 = ms1 + nA1t1 + /А^ь
/
0
6. Численное моделирование процесса
Численное моделирование с использованием программного пакета CST Microwave Studio качественно подтвердили эти результаты в приближении экспоненциальной зависимости коэффициента поглощения среды от длины волны. В этом пакете моделировался
Рис. 2. Небольшой элемент среды (стекла), заполненный порами с водой на 58% в объёмном отношении
небольшой элемент среды (стекла), заполненный порами с водой на 58% в объёмном отношении. На этот элемент падала электромагнитная волна (порт 1 на Рис. 2). Часть электромагнитной волны, которая прошла через этот элемент, регистрировалась с другой стороны элемента. Далее, полученные результаты обрабатывались и сравнивались с теоретическими и экспериментальными зависимостями.
На Рис. 3 изображена численная зависимость коэффициента пропускания от длины волны, с исключёнными длинами волн 432 нм (Т=0,80), 1017 нм (Т=0,86) и 1064 нм (Т=0,67), так как на этих длинах волн наблюдаются провалы в спектре пропускания, связанные с выполнением условий Брэгга для определённых длин волн в данном фотонном
Рис. 3. Качественное подтверждение результатов в приближении экспоненциальной зависимости коэффициента поглощения среды от длины волны
кристалле. Эти провалы наблюдались при численном моделировании в связи с упорядоченностью структуры пор, чего в реальном эксперименте наблюдаться не может, а потому не представляет интереса для нашего случая.
7. Заключение
Предложен метод определения среднего размера характерных неоднородностей среды. Метод был опробован на образце нанопористого стекла и дал удовлетворительные результаты. Было проведено численное моделирование процесса, которое также оказалось в хорошем согласии с опытом. Таким образом, метод обладает сравнительно хорошей точностью, при низкой стоимости и трудоёмкости.
Работа поддержана ФЦП «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы» (контракты NK-526P, 14.740.11.0879, 16.740.11.0030) и грантом 11-08-00267 РФФИ.
Литература
[1] Ван де Хюлст. Рассеяние света малыми частицами — М: ИЛ, 1961.
[2] Чандрасекар С. Перенос лучистой энергии — М: ИЛ, 1953
[3] Шифрин К.С. Рассеяние света в мутной среде — М-Л: Гос. Изд. Техн.-теор. лит., 1951.
[4] Исимару А. Распространение и рассеяние волн в случайно-неоднородных средах. Т.1 — М: Мир, 1981.
[5] Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений — М: Наука, 1966.
[6] Schartl W. Light scattering frompolymer solutions and nanoparticle dispersions — Springer Laboratory. V. XIV, 2007.
[7] Маслов В.П. Операторные методы — M.: Наука, 1973.
[8] Андреева О.В., Обыкновенная Е.И. Нанопористые матрицы — основа композитов различного назначения для оптических экспериментов // Труды научно-исследовательского центра фотоники и Оптоинформа-тики СПбГУ ИТМО: сборник статей под редакцией И.П. Гурова и С.А. Козлова. — 2009. — C. 303-322.
[9] Андреева О.В., Обыкновенная И.Е. Нанопористые матрицы НПС-7 и НПС-17 — возможность использования в оптическом эксперименте // Наносистемы: физика, химия, математика. — 2010. — Т. 1, № 1. — С. 37-53.
[10] J M Pitarke, V M Silkin, E V Chulkov, P M Echenique. Theory of surface plasmons and surface-plasmon polaritons // Institute of physics publishing reports on progress in physics. — 2007. — V. 70. — P. 1-87.
[11] Альфимов А.В., Арысланова Е.М., Вавулин Д.Н., Андреева О.В., Темнова Д.Д., Лесничий В.В., Альмя-шев В.И., Кириллова С.А., Чивилихин С.А., Попов И.Ю., Гусаров В.В. Диффузия света в среде с нано-размерными неоднородностями // Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики. — 2011. —Т. 72. — С. 105-109.
[12] Rubem M. F. Vargas, Cynthia F. Segatto and Marco Tullio Vilhena. On the Analytical Solution of the SN Radiative Transport Equation in a Slab for a Space-dependent Albedo Coefficient // Eurotherm Conference No. 95: Computational Thermal Radiation in Participating Media IV IOP Publishing. Journal of Physics: Conference Series. — 2012. — V. 369. — P. 012018.
[13] Michael F. Modest and Shenghui Lei. The Simplified Spherical Harmonics Method For Radiative Heat Transfer // Eurotherm Conference No. 95: Computational Thermal Radiation in Participating Media IV IOP Publishing. Journal of Physics: Conference Series. — 2012. — V. 369. — P. 012019.
INVESTIGATION OF NANOINHOMOGENEITIES IN SCATTERING
AND ABSORBING MEDIA
D.N. Vavulin, A.V. Alfimov, O.V. Andreeva, E.M. Aryslanova, A.V. Panteleev, S.A. Chivilikhin
In this paper, we propose an original method for determining the size of nanoscale inhomogeneities in scattering and absorbing medium by comparing the theoretical and experimental dependences of the transmission medium on the wavelength of light. The method was tested on a sample of nanoporous glasses and gave satisfactory results. The numerical simulations of light propagation process were also carried out using software package CST Microwave Studio.
Keywords: Rayleigh scattering, nanoscale inhomogeneities, nanoporous glass, the diffusion of radiation.
Dmitrii Vavulin - Saint Petersburg National Research University of Information Technologies Mechanics and Optics, Saint Petersburg, Russia, student, bachelor, dima-vavulin@mail.ru Anton Alfimov - Saint Petersburg National Research University of Information Technologies Mechanics and Optics, Saint Petersburg, Russia, student, bachelor, Alfimov. anton@gmail. com
Olga Andreeva - Saint Petersburg National Research University of Information Technologies Mechanics and Optics, Saint Petersburg, Russia, associate professor, PhD, senior scientist, olga_andreeva@mail. ru
Elizabeth Aryslanova - Saint Petersburg National Research University of Information Technologies Mechanics and Optics, Saint Petersburg, Russia, student, bachelor, Elizabeth. aryslanova@gmail. com
Andrey Panteleev - Saint Petersburg National Research University of Information Technologies Mechanics and Optics, Saint Petersburg, Russia, student, bachelor, pantandrej@yandex.ru
Sergey Chivilikhin - Saint Petersburg National Research University of Information Technologies Mechanics and Optics, Saint Petersburg, Russia, associate professor, PhD, senior scientist, sergey.chivilikhin@gmail.com
