Математические модели в сканирующей микроскопии ближнего поля и их реализация в виде комплекса программ

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Специальность 05.13.18 - Математическое моделирование, численные методы и комплексы программ

Автореферат

диссертации на соискание ученой степени кандидата физико-математических наук

Математические модели в сканирующей микроскопии ближнего поля и их реализация в виде комплекса программ

Беспалова Наталья Викторовна

Саратов 2010

Работа выполнена в ГОУ ВПО «Саратовский государственный технический университет»

Научный руководитель

доктор физико-математических наук, заслуженный деятель науки РФ,

профессор Байбурин Вил Бариевич

Официальные оппоненты:

доктор физико-математических наук, профессор Мельников Леонид Аркадьевич

доктор физико-математических наук, профессор Аникин Валерий Михайлович

Ведущая организация

Управление Российской академии наук Институт проблем точной механики и управления РАН

Защита состоится 19 февраля 2010 г. в 13 часов на заседании диссертационного совета Д212.242.08 при ГОУ ВПО «Саратовский государственный технический университет» по адресу: 410054, г. Саратов, ул. Политехническая, 77, Саратовский государственный технический университет, корп. 1, ауд.319.

С диссертацией можно ознакомиться в научно-технической библиотеке ГОУ ВПО «Саратовский государственный технический университет»или на сайте Саратовского государственного технического университета.

Автореферат разослан « 19 » января 2010 г.

Ученый секретарь диссертационного совета Терентьев А.А.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. В настоящее время одним из основных методов исследования и модификации наноразмерных структур является сканирующая зондовая микроскопия (СЗМ), в том числе, оптическая микроскопия ближнего поля (СОМБП), использующая световые поля рассеяния в ближней зоне. Указанным вопросам посвящены работы зарубежных и отечественных авторов: Binnig G., Rohrer H., Young R., Hansma P., Pohl D., Quate C.F., Чаплыгина Ю.А., Эдельмана В.С., Яминского И.В., Панова В.И., Неволина В.К., Логинова Б.А., Быкова В.А., Емельянова В.И., Булатова А.Н., Вернера В.Д., Байбурина В.Б., Волкова Ю.П. и др. Вместе с тем следует признать, что ряд задач, формулируемых далее, связанных с теорией и практикой применения СЗМ, требуют дополнительных исследований и пока далеки от полного решения.

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

Не менее важной является задача обработки полученных с помощью СЗМ изображений, с целью устранения дефектов и искажений, специфичных для данного типа СЗМ (вздутие, наклон поверхности и др.), вызванных, в частности, прилипанием зонда, тепловым дрейфом, вибрацией острия и др. Используемые для этих целей методы (вычитание плоскости или поверхностей второго порядка, преобразования Фурье с удалением низкочастотных составляющих и др.) зачастую вносят дополнительные искажения. Поэтому возникает необходимость выбора и применения математических методов, позволяющих устранять указанные дефекты.

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

Изложенное определило следующую цель работы.

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

Для достижения этой цели были поставлены и решены следующие задачи:

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

Разработка и реализация новых алгоритмов обработки СЗМ изображений поверхности (вейвлет-преобразования, фрактальный анализ и др.).

Создание комплекса программ для расчета распределения электрических и магнитных полей при рассеянии вблизи поверхности в СОМБП при различном расстоянии от сканирующего элемента (от ближнего до дальнего поля), различном направлении падения волны (из вакуума, из среды), произвольной диэлектрической проницаемости исследуемой поверхности и произвольной поляризации световой волны.

На основе численного моделирования выбор конструкций зондов с улучшенными фокусирующими свойствами.

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

В работе использовались методы решения задач волновой оптики, теория рассеяния света Ми, а также математические методы цифровой обработки сигналов, вейвлет-анализа, теории фракталов. В качестве оборудования использовались тестовые объекты и кантиливеры фирмы НТ - МДТ г. Зеленоград, электронный микроскоп HU-12A и сканирующий зондовый микроскоп (ООО «Пьезон» г. Саратов). Были использованы пакеты программ: Femlab, Matlab, а также «Программа управления универсальным комплексом сканирующей зондовой микроскопии» (авторы: Байбурин В.Б., Волков Ю.П., Якименко Р.А., Большаков А.А., свидетельство об официальной регистрации программы для ЭВМ №2003611643 от 10 июля 2003 г.).

Научная новизна работы:

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

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

Приведен сравнительный анализ различных методов расчета фрактальных размерностей изображений поверхностей, полученных с помощью СЗМ ближнего поля, и разработан математический алгоритм нахождения фрактальной размерности с помощью метода WTMM (метод модулей максимумов вейвлет-преобразования).

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

Проведен сравнительный анализ различных конструкций световодов (в случае идеального металла) и отмечены преимущества коаксиального световода (лучшие фокусирующие свойства).

Разработан комплекс программ, позволяющий:

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

Производить обработку СЗМ изображений поверхностей различными математическими методами (морфологическая фильтрация, Фурье-анализ, вейвлет-анализ и др.) с целью устранения дефектов, вызванных спецификой микроскопа.

Сравнивать различные методы расчета фрактальных размерностей и оценивать сложность изображений при помощи одномерного и двухмерного WTMM анализа.

Научные положения, выносимые на защиту:

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

По сравнению с другими математическими методами (дискретными преобразованиями Фурье, морфологическими фильтрами, вычитанием кривых 1-го и 2-го порядков) вейвлет-преобразования являются более эффективными, позволяя устранять наиболее специфичные для СЗМ искажения (наклон поверхности, вздутие и др.).

Нахождение фрактальной размерности с помощью метода WTMM дает возможность определить спектр фрактальных размерностей исследуемого объекта, позволяя оценить сложность рельефа поверхности объекта (гладкость поверхности, наличие дефектов и пор и др.).

Фокусирующие свойства коаксиальных острий СОМБП (в случае идеального металла) существенно превышают соответствующие параметры традиционных измерительных острий.

Практическая значимость.

Разработано программное обеспечение, позволяющее рассчитывать рассеянное поле любого порядка дифракции (проведены расчеты для 2-5 порядков, при этом увеличение порядка дифракции не приводит к существенным изменениям картины распределения светового поля) для произвольного рельефа поверхности с произвольной диэлектрической проницаемостью при различной поляризации падающего излучения.

Создан программный комплекс управления универсальным зондовым микроскопом, позволяющий производить сканирование, визуализацию изображений и обработку полученных объектов различными математическими методами. Указанное программное обеспечение используется при подготовке и проведении лабораторных работ по дисциплине «Компьютерное моделирование».

Создан программный продукт, исследующий фрактальные свойства заданной поверхности различными математическими методами.

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

Апробация работы. Материалы диссертационной работы доложены и обсуждены на Международной конференции SPIE «Fotonics West» (San Jose USA 2002), Международной конференции «Saratov fall meeting» (Саратов, 2002, 2009), Российском симпозиуме по растровой электронной микроскопии и аналитическим методам исследования твердых тел РЭМ (Черноголовка, 2002, 2007), ХХI Международной научной конференции «Математические методы в технике и технологиях» (Саратов, 2008), Международной научно-технической конференции «АПЭП-2008» (Саратов, 2008), Всероссийской научной школе-семинаре «Методы компьютерной диагностики в биологии и медицине - 2009» (Саратов, 2009), Международном форуме по нанотехнологиям «Rusnanotech - 09» (Москва, 2009).

Публикации. По теме диссертации опубликовано 15 работ (статьи в журналах, рекомендованных ВАК РФ - 2, статей в научных сборниках - 13), получено 2 свидетельства об официальной регистрации программы для ЭВМ.

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

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

Диссертация состоит из введения, четырех глав, заключения, списка использованной литературы.

Во введении обосновывается актуальность диссертационной работы, сформулированы цель и основные задачи исследований, указаны новизна, научная и практическая значимость результатов работы, приводятся основные положения, выносимые на защиту.

В первой главе описаны различные конструкции и принципы действия сканирующих зондовых микроскопов ближнего поля: туннельный микроскоп, атомно-силовой микроскоп и сканирующий оптический микроскоп ближнего поля (СОМБП).

Рассматриваются различные используемые на практике варианты зондов СОМБП и описываются известные способы моделирования формирования изображений, получаемых с помощью оптических микроскопов ближнего поля. Приводится описание различных программных продуктов, используемых для обработки и анализа изображений, получаемых с помощью зондовых микроскопов различного типа. Описаны некоторые перспективные методы обработки микроскопических данных (дискретное вейвлет-преобразование, фрактальный анализ, метод максимальных модулей вейвлет-преобразования и др.).

Во второй главе развиты математические модели расчета рассеяния световых полей в ближней зоне как на эталонных поверхностях, так и поверхностях различного типа.

Вначале, используя метод работы [Емельянов В.И., Семиногов В.Н., Соколов В.И. //Квантовая электроника. -1987. -Т.14. №1. -С.33-46], получим уравнения амплитуд проекций рассеянного оптического поля в ближней зоне для эталонной поверхности с синусоидальным рельефом в соответствии со схемой на (рис. 1,а), т.е. в случае падения света сверху на границу раздела вакуума и среды.

а б

Рис. 1. Схематическое изображение световой волны, падающей на синусоидальную поверхность: а - из вакуума; б - из среды

Пусть плоская электромагнитная волна распространяется в среде под углом и падает из вакуума на синусоидальную поверхность раздела изотропной и однородной среды с комплексным диэлектрическим показателем .

(1)

Синусоидальная по рельефу среда занимает нижнее полупространство

(2)

где - амплитуда Фурье-компоненты рельефа поверхности; q - вектор синусоидальной решетки (q=2, где T- период синусоиды); r={x,y}- радиус-вектор, лежащий в плоскости z=0; kt, kz - проекции волнового вектора падающей плоской волны k0 (k0=/c) на оси Z и Y соответственно; - угол между падающей волной и осью Z.

Необходимо найти решение для составляющих электрических и магнитных полей, прошедших через границу раздела и отраженных в вакуум. Для этого, предполагая гипотезу Релея справедливой, получаем для отраженного в среду и прошедшего в воздух полей выражения (3) и (4) соответственно:

(3)

(4)

(5)

где p - порядок дифракции, Re >0 и Re>0, если ; Im <0, и Im<0, если (p=0 описывает отраженную в воздух и прошедшую в среду Френелевские волны).

Амплитуды проекций рассеянных полей Epx, Epz и прошедших в среду полей , находятся из решения следующей линейной системы уравнений:

(6)

(7) (8)

(9)

где Jp(x)- функция Бесселя первого рода; ,. Для вычисления Y компонент поля и могут быть использованы следующие уравнения (полученные на основании уравнений Максвелла):

(10)

Магнитные компоненты электромагнитных полей, рассеянного в воздух и прошедшего в среду, могут быть вычислены из уравнений (3)-(5) и уравнений Максвелла: и . Уравнения (6)-(10) представляют математическую модель рассеяния оптического поля в ближней зоне для эталонной поверхности с синусоидальным рельефом.

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

Пусть плоская электромагнитная волна распространяется в среде под углом и падает изнутри на синусоидальную поверхность раздела с воздухом изотропной и однородной среды с комплексным диэлектрическим показателем .

(11)

Синусоидальная по рельефу среда занимает нижнее полупространство

(12)

Для отраженного в среду и прошедшего в воздух полей получаем выражения (13) и (14) соответственно:

(13)

(14)

где

, . (15)

Амплитуды проекций прошедших полей , и отраженных в среду полей Epx, Epz, входящих в выражения (13)-(15), находятся из решения следующей системы линейных уравнений:

(16)

(17)

(18)

(19)

Для получения Y - компоненты полей и используется следующая система уравнений, полученная из уравнений Максвелла:

(20)

Магнитные компоненты электромагнитного поля, рассеянного в воздухе и прошедшего в среду, вычисляются из уравнений (13), (14), (15) и уравнений Максвелла: и .

Для решения полученных систем линейных алгебраических уравнений (6)-(9) и (16)-(19), в отличие от метода [Емельянов В.И., Семиногов В.Н., Соколов В.И. //Квантовая электроника. -1987. -Т.14. №1. -С.33-46], используется численный метод исключения Гаусса, что позволяет получить решение при произвольной диэлектрической проницаемости поверхности и в случае больших амплитуд (порядка 10 нм). Полученные решения дают возможность рассчитать объемную плотность энергии светового ближнего поля в соответствии с соотношением:

. (21)

На основе данного метода создана программа расчета и визуализации изображения поверхности с нанометровым рельефом, выполненная на языке программирования Delphi (свидетельство о регистрации программы №2009615620 от 08 октября 2009 г.). При этом распределению плотности энергии от максимума до минимума (от высоких до низких точек рельефа) на изображении поверхности соответствует градация цвета от светлого до темного, соответственно.

На рис. 2 приведены изображение заданной эталонной поверхности (рис. 2,а) с синусоидальным рельефом и изображение, полученное в результате расчетов, проведенных по полученным выше соответствиям (рис. 2,б).

а б

Рис. 2. Синусоидальная поверхность: а - эталонная; б - полученная в результате расчета

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

а б в

Рис. 3. Картина плотности энергии рассеяния света на металлической синусоидальной поверхности с амплитудой 10 нм и периодом 50 нм. Поверхность освещается сверху под углом 30°, длина волны падающего излучении 630 нм, s-поляризация: а - расстояние 11 нм от поверхности; б - расстояние 30 нм от поверхности; в - расстояние 60 нм от поверхности

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

Следует отметить, что при получении изображений в ближнем поле с помощью СЗМ возникают дефекты (вздутие поверхности, наклон, шумы и др.), вызванные вибрацией острия, тепловым дрейфом и др. и искажающие вид рассматриваемой поверхности.

В третьей главе рассмотрены различные математические методы (вычитание наклона произвольной поверхности, преобразования Фурье, вейвлет-преобразования и т.д.), направленные на устранение указанных искажений.

Сравнение показало, что более удобным методом являются вейвлет-преобразования, способные адекватно обрабатывать как регулярные, так и сингулярные поверхности. В работе, в качестве рабочего вейвлета, был использован вейвлет Добеши 8-го порядка, с помощью которого обрабатывались различные изображения.

На рис. 4 представлены результаты обработки изображения полированной поликоровой пластины различными математическими методами. Аппроксимация полиномом второго порядка позволяет избавиться от наклона, однако не исключает вздутия; преобразования Фурье вообще не вносят значительной корректировки, а вейвлет-преобразования демонстрируют исчезновение дефектов вздутия и наклона.

А б

В г

Рис. 4. СЗМ изображение полированной поликоровой пластины (бугры на поверхности представляют собой увеличенные изображения следов полировки абразивным порошком): а - до обработки; б - после аппроксимации полиномом 2-го порядка; в - после обработки с помощью быстрого преобразования Фурье; г - после обработки вейвлетом Добеши 8-го порядка (очистка от дефектов)

В главе также представлены различные способы расчета фрактальных размерностей объектов, по величине которых можно судить о сложности рельефа поверхности (гладкость, наличие пор, дефектов).

В работе проведен расчет фрактальных размерностей существующими методами (подсчет кубов, треугольников, метод вариаций, спектральной мощности, WTMM). Расчеты показали, что наиболее удобным методом исследования структуры является метод модулей максимумов вейвлет-преобразования (wavelet transform modulus maxima - (WTMM)), поскольку позволяет не только рассчитать фрактальную размерность, но и определить, является ли изображение моно- или мультифракталом. Было создано соответствующее программное обеспечение, выполненное на языке программирования С++. При этом поверхность подвергается вейвлет-преобразованию (сглаживанию) по формуле

, (22)

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

, (23)

где L(a) - множество всех линий (l) локальных максимумов модулей вейвлет-коэффициентов, существующих на масштабе a, значение q выбирается в диапазоне (-n,n), где n - целое число.

При малых а, z можно представить как

~, (24)

где ф(q) - скейлинговые экспоненты, которые определяются для некоторого значения q путем вычисления наклона зависимости . Выбирая различные степени q, при построении частичных функций, можно получить линейную зависимость ф(q) с постоянным значением экспоненты Гёльдера Н=dф(q)/dq=const в случае монофрактальных объектов и нелинейную зависимость с большим числом гёльдеровских экспонент h(q)=dф(q)/dq?const в случае мультифракталов. На рис. 5,а в качестве примера задано монофрактальное изображение. Соответствующая ей, полученная из расчетов зависимость dф(q)/dq, как следует из рис. 5,б, является константой.

А б

Рис. 5: а - изображение «лестницы Кантора»; б - соответствующий график зависимости dф(q)/dq

На рис. 6,а представлено мультифрактальное изображение, поскольку из рис. 6,б видно, что соответствующая ему зависимость dф(q)/dq является не константой, а существенно нелинейной.

А б

Рис. 6: а - изображения поверхности полированного кремния, покрытого пленкой алюминия; б - соответствующий график зависимости dф(q)/dq

На языке C++ Builder было написано программное обеспечение, позволяющее производить расчеты фрактальных размерностей, как для одномерных, так и для двумерных данных различными методами (метод кубов, метод треугольников, метод Фурье, метод WTMM (одномерный и двумерный)).

Четвертая глава посвящена выбору конструкций оптических элементов (микролинз, зондов) в СОМБП. Для расчета фокусирующих характеристик микролинз в работе использовалась теория Ми. В соответствии с теорией предполагалось, что плоская линейно поляризованная монохроматическая волна дифрагирует на сфере радиуса a, помещенной в изотропную непроводящую и немагнитную среду. Сфера считается немагнитной, но может быть как проводящей, так и диэлектрической. Амплитуда электрического вектора падающей волны нормирована на 1.

В соответствии с теорией Ми, компоненты векторов поля рассеянной частицей волны в сферической системе координат записываются в виде бесконечных знакопеременных рядов, количество членов которых зависит от размера частиц и заданной точности.

Как видно из рис. 7, на котором приведены результаты расчета для диэлектрических сфер (стекло с показателем преломления 1,25) различного диаметра, наилучшие фокусирующие свойства имеют микролинзы с диаметром не менее 10 мкм.

Было создано соответствующее программное обеспечение, выполненное на языке программирования Delphi (имеется свидетельство о регистрации программы №2009615619 от 08 октября 2009 г.), позволяющее производить расчеты как для диэлектриков, так и проводников в ближнем и дальнем поле, при этом отдельно производится расчет для электрической и магнитной составляющих, а также для плотности энергии. Программа позволяет рассчитывать индикатрису рассеяния. Время счета по программе зависит от диаметра исследуемой микролинзы и составляет порядка 30 секунд для микролинзы диаметра 10 мкм.

а б в г

Рис. 7. Изображение распределения плотности энергии электромагнитного поля =0,63 мкм, рассеянного стеклянной сферической частицей: а - радиуса 0,01 мкм с показателем преломления 1,25; б- радиуса 3 мкм с показателем преломления 1,25; в - радиуса 6 мкм с показателем преломления 1,25; г - радиуса 30 мкм с показателем преломления 1,25

микроскопия математический световой поле нанометровый

На основе расчетов, с использованием волнового уравнения, по программе в среде Femlab было проведено сравнение конструкций как коаксиального, так и обычного световодов по коэффициенту пропускания и степени фокусировки. Расчеты были проведены для электрической, магнитной компонент поля и плотности энергии.

Как видно из рис. 8, коаксиальный световод имеет лучшие фокусирующие свойства (отношение диаметров фокальных пятен обычного световода ds и коаксиального световода dk на расстоянии 0,2 мкм , на расстоянии 2 мкм ), кроме этого, пучок света в коаксиальном световоде лучше локализован, что приводит к уменьшению потерь.

А б

Рис.8. Плотность электрической компоненты поля, рассеянного: а - коаксиальным световодом, покрытым металлом; б - стеклянным световодом, покрытым металлом

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

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Рис. 9. Комплекс программ для сканирующей зондовой микроскопии ближнего поля

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

ОСНОВНЫЕ ВЫВОДЫ ПО РАБОТЕ

Основные результаты работы заключаются в следующем:

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

Сравнительный анализ различных математических методов (морфологическая фильтрация, Фурье-анализ, вейвлет-анализ и др.) показал преимущество вейвлет-преобразования, позволяющего устранять типичные для СЗМ искажения изображений поверхностей.

Разработан рабочий алгоритм расчета фрактальных размерностей изображений поверхностей на основе метода WTMM (метод модулей максимумов вейвлет-преобразования), отличающийся от известных возможностью оценить мультифрактальность исследуемых изображений поверхностей.

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

На основе примененных математических моделей и рабочих алгоритмов разработан комплекс программ, позволяющий: рассчитывать плотности энергии рассеянной световой волны для получения изображения поверхностей; производить обработку СЗМ изображений поверхностей различными математическими методами с целью устранения дефектов изображений; производить расчет фрактальных размерностей различными методами и оценивать мультифрактальность изображений при помощи WTMM метода.

СПИСОК ОПУБЛИКОВАННЫХ РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Публикации в центральных изданиях, включенных в перечень периодических изданий ВАК РФ

Беспалова Н.В. Алгоритмы обработки изображений, используемые сканирующей зондовой микроскопией/ В.Б. Байбурин, Ю.П. Волков, Н.В. Беспалова// Вестник Саратовского государственного технического университета. 2007. №3(27). Вып. 2. С. 61-64.

Беспалова Н.В. Моделирование оптического поля в ближней зоне/ Н.В. Беспалова// Вестник Саратовского государственного технического университета. 2009. №4(43). Вып. 2. С. 20-24.

2. Публикации в других изданиях

Беспалова Н.В. Фокусирующие свойства сферических микролинз /Ю.П. Волков, В.Б. Байбурин, Ю.А. Аветисян, Н.В. Беспалова// Автоматизация и управление в машино- и приборостроении: межвуз. науч. сб. / Саратов: СГТУ, 2002. С. 31-34.

Bespalova N.V. Optical properties of spherical microlenses / U. P. Volkov, U. A. Avetisjan, V. B. Baibyrin, N.V. Bespalova // SPIE Proc. 2002. Vol. 4705. P. 145-155.

Bespalova N.V. Numerical simulation of image formation in near-filed optical microscopy / V. B. Baiburin, Y. P. Volkov, U.A. Avitisyan, I. V. Krasnikova, N.V. Bespalova // SPIE Proc. 2002. Vol. 4621. P. 117-120.

Беспалова Н.В. Фокусирующие свойства диэлектрических и металлических микролинз / Ю. П. Волков, В. Б. Байбурин, Ю. А. Аветисян, Н.В. Беспалова // Электронные приборы и устройства нового поколения: материалы науч.-техн. конф. / Саратов: Изд-во Сарат. ун-та, 2002. С. 149-152.

Беспалова Н.В. Обработка СЗМ изображений методами вейвлет анализа. / В.Б. Байбурин, Ю.П. Волков, Н.В. Беспалова//Тез. 15-го Рос. симп. по растровой электронной микроскопии и аналитическим методам исследования твердых тел (РЭМ-2007). Черноголовка, 2007. С. 21.

Беспалова Н.В. Измерение фрактальных размерностей в микроскопических изображениях / В.Б. Байбурин, Ю.П. Волков, Н.В. Беспалова// Математические методы в технике и технологиях - ММТТ-21 : сб. тр. XXI Междунар. науч. конф.: в 10 т. / Саратов: СГТУ, 2008. Т. 7. C. 67-70.

Беспалова Н.В. Гранулометрические методы анализа микроскопических данных / В.Б. Байбурин, Ю.П. Волков, Н.В. Беспалова, Е.А. Головешко// Математические методы в технике и технологиях - ММТТ-21 : сб. тр. XXI Междунар. науч. конф.: в 10 т. / Саратов: СГТУ, 2008. Т. 7. C. 65-67.

Беспалова Н.В. Использование ПИД - алгоритма регулирования в процессе исследования поверхности сканирующим зондовым микроскопом. / В.Б. Байбурин, Н.В. Беспалова, Ю.П. Волков// «АПЭП-2008»: сб. тр. Междунар. науч.-техн. конф. Саратов: СГТУ, 2008. С. 389-393.

Беспалова Н.В. Оборудование для преподавания нанотехнологий с использованием сети Интернет / В.Б. Байбурин, Н.В. Беспалова, Ю.П. Волков// Интернет и инновации: практические вопросы информационного обеспечения инновационной деятельности: сб. тр. Междунар. науч.-практ. конф. Саратов: СГТУ, 2008. С. 198-199.

Беспалова Н.В. Методы моделирования нанометровых оптических структур / В.Б. Байбурин, Н.В. Беспалова, В.С. Морев, Ю.П. Волков// Гетеромагнитная микроэлектроника. Вып. 6. Гетеромагнитная микро- и наноэлектроника. Системы информационной безопасности. Саратов: Изд-во Сарат. ун-та, 2009. С. 61-66.

Беспалова Н.В. Моделирование изображения поверхности в ближнем поле / Н.В.Беспалова, В.Б. Байбурин, Ю.П. Волков// Методы компьютерной диагностики в биологии и медицине - 2009: материалы ежегодной Всерос. науч. школы-семинара, под ред. проф. Д.А. Усанова. Саратов: Изд-во Сарат. ун-та, 2009. С. 31-35.

Беспалова Н.В. Моделирование изображения нанометровых структур в оптической микроскопии ближнего поля. / В.Б. Байбурин, Ю.П. Волков, Н.В. Беспалова, Н.П. Конов// Гетеромагнитная микроэлектроника. Выпуск 7. Гетеромагнитная микро- и наноэлектроника. Системы информационной безопасности. Саратов: Изд-во Сарат. ун-та, 2009. С. 71-83.

Bespalova N.V. Near-field optical image simulation / N.V. Bespalova, V.B. Bayburin, U.P. Volkov// Nanotechnology international forum Rusnanotech 09. Moscow, 2009. P. 232-235.

Беспалова Н.В. Программа моделирования изображения поверхности в оптике ближнего поля / В.Б. Байбурин, Н.В. Беспалова, Ю.П. Волков // Свидетельство об официальной регистрации программы для ЭВМ №2009615620 от 08 октября 2009 г.

Беспалова Н.В. Программа расчета геометрических характеристик фокального пятна для наноразмерных сферических линз / В.Б. Байбурин, Н.В. Беспалова, Ю.П. Волков// Свидетельство об официальной регистрации программы для ЭВМ №2009615619 от 08 октября 2009 г.

Размещено на Allbest.ru

...