CC BY
18
УДК 548.73
В.Г. Кон1, М.А. Орлов2 1 Национальный исследовательский центр «Курчатовский институт»
2 Московский государственный университет им. М.В. Ломоносова
Корректировка изображения в методе фазового контраста Цернике для жесткого рентгеновского излучения
Проведено численное моделирование эксперимента по визуализации прозрачных (фазовых) микрообъектов с использованием метода фазового контраста Цернике в жестком рентгеновском излучении. В расчетах использованы параметры пучков, характерные для источников синхротронного излучения третьего поколения. В качестве фокусирующего элемента рассмотрена преломляющая параболическая линза. Фазосдвигающая четвертьволновая пластинка расположена в области изображения точечного источника. Показано, что изображение объектов может быть скорректировано путем деления регистрируемого детектором профиля интенсивности на профиль интенсивности, формируемый сразу за линзой из-за поглощения.
Ключевые слова: численное моделирование, жесткое рентгеновское излучение, фазовый контраст Цернике, преломляющая параболическая линза, фазосдвигающая четвертьволновая пластинка, корректировка изображения.
I. Введение
В последние 15 лет активно развивается метод рентгеновского фазового контраста, позволяющий неразрушающим образом исследовать внутреннюю структуру слабопоглощающих некристаллических объектов. Различают две основных разновидности этого метода: метод фазового контраста с кристаллом-анализатором [1], применяющийся для исследования объектов субмил-лиметровых размеров и имеющий принципиальное ограничение по разрешению, составляющее несколько микрон, и метод фазового контраста в in-line схеме [2], применяющийся для исследования объектов с меньшими размерами и требующий высокой степени когерентности падающего пучка рентгеновского излучения. Принципиально другим подходом является использование метода фазового контраста Цернике для жесткого рентгеновского излучения. На данный момент не сделано ни одного эксперимента, в котором этот метод применялся бы с преломляющей линзой в качестве объектива. В опубликованных работах объективом служила более распространенная в рентгенооптических исследованиях зонная пластинка [3-5]. В работе [6] был проведен теоретический анализ соответствующих экспериментов, показавший, что зонная пластинка обладает худшим разрешением (около 3 мкм), чем преломляющая линза, позволяющая получить разрешение меньше 1 мкм. Это объясняется наличием у зонной пластинки большого количества порядков фокусировки. Тем не менее из результатов работы следовало, что рентгеновская преломляющая линза обладает существенным недостатком, а именно плохо изображает объекты, удаленные от центра апертуры. Это связано с тем, что толщина линзы на краях существенно больше, чем в центре, и значительная часть излучения, проходящего через края апертуры, поглощается. Однако, как оказалось, изображение легко скорректировать путем деления регистрируемого детектором профиля интенсивности на профиль интенсивности, формируемый сразу за линзой из-за поглощения. Удивительный факт заключается в том, что при этом не только выравнивается фон, но и контраст выравнивается сам собой. Такая корректировка позволит, в частности, наблюдать с помощью метода Цернике движущиеся объекты субмиллиметрового размера. Таким образом, линза в качестве объектива превосходит зонную пластинку, поскольку обладает лучшим разрешением и после корректировки позволяет получить неискаженное изображение вдоль практически всей апертуры.
II. Метод расчета
При решении задач распространения жесткого рентгеновского излучения с высокой точностью выполняется параксиальное приближение. При этом перенос волновой функции излучения в воздухе описывается с помощью интегральной формы Кирхгофа для решения уравнения Максвелла. Пусть ось г декартовой системы координат направлена вдоль оптической оси, то есть вдоль распространения пучка излучения. Задача состоит в вычислении зависимости амплитуды волнового поля от поперечных координат х и у в каждой точке на оси г. При этом характерный интервал существенного изменения волнового поля вдоль оси г равен десяткам сантиметров, в то время как в поперечных направлениях волновое поле меняется на расстояниях в несколько микрон и меньше. Так как поляризация в рассматриваемых нами процессах не изменяется, то достаточно ограничиться скалярной волновой функцией поля. Более того, для упрощения задачи мы будем рассматривать только одномерные объекты, однородные вдоль оси у.
Пусть Е1(х) — волновая функция в точке гь Тогда волновая функция Е(х) в точке г при условии, что между г1 и г нет объектов, определяется следующим образом:
Е(х) =
1 Л х
^Р(х - хьг - г1)Е1(х1), Р(х,г) = (^дг)1/2 ехР ^ ) • (1)
„2
Здесь Р(х, г) — пропагатор Кирхгофа в параксиальном приближении, Л — длина волны излучения. Что касается объектов, то если их продольные размеры не слишком сильно отличаются от поперечных, мы можем пренебречь изменением траектории лучей при прохождении объекта, поскольку углы рассеяния весьма малы (обычно — десятки микрорадиан). При этом можно пренебречь длиной объекта вдоль оптической оси и считать его плоским в плоскости его середины. Тогда взаимодействие излучения с таким объектом описывается путем умножения волновой функции поля излучения на так называемую трансмиссионную функцию:
Т(х) = ехр ^—г[5 — гв]“Л^(х)^ , (2)
где £(х) — локальное изменение толщины объекта вдоль оси г, при условии, что он состоит из одного материала, а 5 и в — действительная и мнимая части комплексного показателя преломления п = 1 — 5 + гв материала, из которого состоит объект.
Вычисления начинаются с точечного источника. Волновая функция в плоскости источника формально может быть взята в виде дельта-функции Ео(х) = 5(х). Подставляя это в (1), сразу получаем волновую функцию перед исследуемым объектом в виде пропагатора Кирхгофа Е1(х) = Р(х, п), где Г1 — расстояние от источника до объекта. Далее необходимо умножить волновую функцию на трансмиссионную функцию объекта, в которой учесть конкретный вид зависимости £(х), и применить формулу (1) еще раз. Однако в формуле (1) пределы интегрирования бесконечные, в то время как в конкретных численных расчетах область интегрирования конечна. Поэтому для проведения расчетов мы будем предполагать, что трансмиссионная функция не меняется за пределами конкретной области интегрирования и равна константе, то есть Т(х) = То, где Т0 описывается формулой (2) при £(х) = £0. В частности, может быть £0 = 0. Тогда интеграл (1) можно преобразовать к виду
Е2(х) = ТоР (х,Г2) +
^Р (х — х1,Г2 — п)[Т (х)Е1(х1) — ТоР (хьп)]. (3)
В этой формуле Г2 — расстояние от точки наблюдения до источника. Кроме того, в общем случае предполагается, что Е1(х) = Р(х, г1) за пределами области интегрирования. В первом расчете это равенство выполняется и внутри области интегрирования, но при следующих расчетах ситуация более сложная.
В первом расчете Г2 равно расстоянию от источника до фокусирующей линзы (рис. 1). Затем расчет по формуле (3) нужно повторить с новым объектом в виде фокусирующей линзы. При
этом надо сделать замену т2 ^ т3, п ^ т2, где т3 — расстояние от источника до плоскости его фокусировки, в которой стоит фазосдвигающий элемент. В случае двояковогнутой преломляющей параболической линзы £(х) = х2/К, где К — радиус кривизны у вершины параболы. В данном случае бесконечные пределы интегрирования эффективно обрезаются из-за поглощения в линзе, если линза на краях апертуры достаточно толстая. В третьем расчете объектом является фазосдвигающий элемент с размерами порядка размера пучка в фокусе линзы. В этом случае в формуле (3) надо сделать замену г2 ^ т4, т1 ^ г3, где т4 — расстояние от источника до координатного детектора.
Интеграл (3) представляет собой свертку двух комплексных функций, и его прямой расчет является достаточно сложной задачей. Поэтому расчет проводился с использованием преобразования Фурье. Сначала вычислялся фурье-образ функции в квадратных скобках, затем он умножался на фурье-образ пропагатора Кирхгофа, который имеет аналитический вид Р(д, т) = ехр(—гЛт д2/4п), и вычислялось обратное преобразование Фурье. Преобразование Фурье рассчитывалось с использованием известной процедуры быстрого преобразования Фурье [7].
Рис. 1. Схема фазового контраста Цернике при использовании преломляющей линзы. Почти параллельный пучок СИ падает слева. 1 — объект, 2 — преломляющая линза, 3 — фазосдвигающая пластинка, 4 — координатный детектор. При изображении без увеличения объект и детектор ставятся на двойном фокусном расстоянии от линзы, а фазосдвигающая пластинка — в точке фокусировки источника
III. Схема численного эксперимента, результаты и выводы
Численный эксперимент был проведен в соответствии со схемой, изображенной на рис. 1. Источник жесткого рентгеновского излучения с размером Бо удален от объекта наблюдения 1 на расстояние Т1. Объект 1 находится на двойном фокусном расстоянии от рентгеновской преломляющей линзы 2, то есть Т2 — Т1 = 2¥. Фазосдвигающая пластинка 3 помещается в точке изображения источника Тз = Т2 + ¥/(1 — ¥/т2), определяемой по формуле линзы. Детектор 4 располагается на расстоянии Т4 — Т2 = 2¥ от линзы, вследствие чего изображение получается перевернутым и без увеличения.
Для расчетов были выбраны стандартные параметры синхротронных источников 3-го поколения, а именно расстояние от источника до объекта п = 50 ми размер источника £о = 50 мкм. В качестве объектива была рассмотрена бериллиевая линза со следующими параметрами: фокусное расстояние ¥ = 20 см, длина линзы р = 5 см, радиус кривизны линзы Ко = 0.88 мкм, геометрическая апертура Ао = 400 мкм. Расчеты проводились для энергии падающего излучения Е = 12.4 кэВ, что соответствует длине волны Л = 0.1 нм. Источник является некогерентным в том смысле, что каждая его точка светит независимо по фазе. Так как источник имеет конечный размер Бо, то его изображение в фокусе также будет иметь определенный размер Б = Б0(т3 — Т2) / (т1 + Т2), и в нашем случае этот размер равен 0,2 мкм, что превышает размер пучка при фокусировке точечного источника (дифракционный предел).
Учет размеров источника проводился вычислением свертки распределения конечной интенсивности от точечного источника с гауссовой кривой, полуширина которой равна Б^ = Бо(Т4 — Т2) / (т1 + Т2). Радиус фазосдвигающей пластинки также подбирался с учетом значения Б. Чтобы выделить изображение объекта в чистом виде, в экспериментах часто применяют специальный прием, а именно записывают изображение в отсутствие объекта и с объектом и затем вычита-
ют первое изображение из второго. При построении графиков распределения интенсивности мы использовали такой же прием, то есть для получения изображения объектов из полного профиля интенсивности проводилось вычитание профиля интенсивности, формируемого линзой в отсутствие объекта.
В работе [6] утверждалось, что качество изображения объектов заметно ухудшается при удалении объектов от центра апертуры. Однако оказалось, что изображение легко восстановить путем деления профиля интенсивности, формируемого линзой, на функцию поглощения в линзе exp
(-A[4nx2^Ro]e), где в = 5,19 • 10
-10
индекс поглощения для бериллия. Из рис. 2 видпо, что
при этом изображение восстанавливается на значительной части апертуры. Поясним значение коэффициента А. Как оказалось, внесение фазовой пластинки само по себе приводит к некоторому искажению профиля интенсивности в плоскости изображения (рис. 3). Он опускается в центре и поднимается на краях. Такое поведение не имеет простой интерпретации, однако ясно, что в силу закона сохранения энергии интегральная интенсивность должна сохраняться. Поскольку нашей целью является получение правильного изображения объектов, в показатель экспоненты следует добавить множитель А, немного меньший 1 (в нашем случае — 0,95). Понятно, что реальные линзы будут неидеальными. И если отклонение от идеальности будет велико, то наилучший вид корректирующей функции можно будет определять непосредственно из эксперимента, беря тестовый образец с известными параметрами и определяя, как искажается его изображение в разных точках х.
•100 0 100 Координата х, мкм
Рис. 2. Изображение параболической преломляю щей линзой с апертурой 400 мкм серии кремни- Координата х, мкм
евых объектов прямоугольной и эллиптической Рис. 3. Изображение, даваемое: (а) «чистой» лин-формы, с продольными размерами 0,5 и 0,1 мкм: зой; (б) линзой при наличии фазовой пластинки.
(а) до корректировки; (б) после корректировки. На краях видны классические брэгговские осцилляции, возникающие из-за наличия резкой границы линзы
Во втором случае все объекты, кроме крайних, изображаются без существенных искажений
Таким образом, линза, будучи более точным оптическим инструментом, чем зонная пластинка, обладает несомненными преимуществами перед зонной пластинкой и в рентгеновском диапазоне, однако здесь необходимо учитывать ее специфику, а именно сильное поглощение на краях, и делить результирующий профиль интенсивности на функцию поглощения в линзе.
Литература
1. Ингал В.В., Беляевская Е.А. Метод фазодисперсионной интроскопии // ЖТФ. — 1997. — Т. 61, вып. 1. — С. 68-77.
2. Snigirev A., Snigireva I., Kohn V., Kuznetsov S., Schelokov I. On the possibilities of x-ray phase contrast microimaging by coherent high-energy synchrotron radiation // Rev. Sci. Instrum. — 1995. — V. 66, N. 12. — P. 5486-5492.
3. Awaji M., Suzuki Y., Takeuchi A., Takano H., Kamijo N., Tamuraa S., Yasumoto M.. Zernike-type X-ray imaging microscopy at 25 keV with Fresnel zone plate optics // J. Synchrotron Rad. — 2002. — N. 9. — P. 125--127.
4. Hiroki Yokosuka, Norio Watanabe, Takuji Ohigashi, Yasutoshi Yoshida, Shunichi Maeda, Sadao Aoki, Yoshio Suzuki, Akihisa Takeuchi, Hidekazu Takano. Zernike-type phase-contrast hard X-ray microscope with a zone plate at the Photon Factory // J. Synchrotron Rad. — 2002. — N. 9. — P. 179-181.
5. Chu Y.S., Yi J.M., De Carlo F., Shen Q., Lee Wah-Keat, Wu H.J., Wang C.L., Wang J.Y., Liu C.J., Wang C.H., Wu S.R., Chien C.C., Hwu Y., Tkachuk A., Yun W., Feser M., Liang K.S., Yang C.S., Je J.H., Margaritondo G. Hard-x-ray microscopy with Fresnel zone plates reaches 40 nm Rayleigh resolution // Appl. Phys. Lett. — 2008. — N. 92, 103119-103121.
6. Кон В.Г., Орлов М.А. Численное моделирование фазового контраста Цернике в жестком рентгеновском излучении с использованием преломляющих линз и зонных пластинок // Поверхность. Рентгеновские, синхротронные и нейтронные исследования. — 2010. — Вып. 11. — С. 76-81.
7. Интернет-ресурс : http://alglib.sources.ru/fft
Поступила в редакцию 27.04.2011.
