Научная статья на тему 'Рекуррентные соотношения в динамической теории рентгеновской дифракции в латерально ограниченном кристалле'

Рекуррентные соотношения в динамической теории рентгеновской дифракции в латерально ограниченном кристалле Текст научной статьи по специальности «Физика»

CC BY
234
44
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДИНАМИЧЕСКАЯ ТЕОРИЯ ДИФРАКЦИИ / ЛАТЕРАЛЬНО ОГРАНИЧЕННЫЙ КРИСТАЛЛ / РЕКУРРЕНТНЫЕ СООТНОШЕНИЯ / DYNAMICAL DIFFRACTION THEORY / LATERAL CRYSTAL / RECURRENCE RELATIONS

Аннотация научной статьи по физике, автор научной работы — Колосов С.И., Пунегов В.И.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

RECURRENCE RELATIONS IN THE DYNAMICAL THEORY OF X-RAY DIFFRACTION IN A LATERAL CRYSTAL

Darwin’s dynamical theory of X-ray diffraction is considered to the cases of plane-parallel and lateral (i.e., having a finite length in the lateral direction) crystals. Therefore this new approach can be widely used for non-destructive testing of lateral structures used in optoand microelectronic devices, nonvolatile memory devices and X-ray optics. We provided a short review of Darwin’s approach in the case of plane-parallel crystals and analyzed the exact solution of this problem. The main consideration is devoted to the dynamical X-ray diffraction in the lateral crystal. We restrict ourselves to the case of a symmetrical coplanar Bragg diffraction. Recurrence relations describing the dynamical X-ray diffraction in the lateral crystal is obtained. The angular distribution of X-ray scattering intensity in the reciprocal space is investigated. Phase changes of X-ray waves in the crystal are analyzed in detail. It is shown that this approach allows one to calculate rocking curves as well as Reciprocal Space Maps (RSMs) for lateral crystals having rectangular cross-section. An algorithm for computing of rocking curves and RSMs from the lateral crystal is developed. In accordance with the model, we used a rectangular lattice having a fixed distance between nodes to describe the dynamical diffraction process. The simulation procedure based on recurrence relations consists of the external and internal cycles. The external cycle starts with the first column and goes up to the last column. The internal cycle (from top to bottom) starts with the first reflecting plane and goes up to the last one. Then substituting the amplitudes calculated for the first column we can calculate the amplitudes for the second column and so on. At the end of the external cycle we obtain arrays of the reflected and transmitted amplitudes. We perform numerical modelling for crystals with different lateral sizes. The obtained simulation results are in good agreement with the solution obtained by integration of Takagi’s equations. It is demonstrated that the kinematical approximation is valid for thick crystals having small length in the lateral direction. It is shown that this approach, being simple and transparent, is faster than the one based on Takagi’s equations and allows simulations of RSMs. The obtained results can potentially be extended further to the three-dimensional case in both the Fourier space and real space.

Текст научной работы на тему «Рекуррентные соотношения в динамической теории рентгеновской дифракции в латерально ограниченном кристалле»

ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

УДК 548.732

РЕКУРРЕНТНЫЕ СООТНОШЕНИЯ В ДИНАМИЧЕСКОЙ ТЕОРИИ РЕНТГЕНОВСКОЙ ДИФРАКЦИИ В ЛАТЕРАЛЬНО ОГРАНИЧЕННОМ КРИСТАЛЛЕ

С.И. КОЛОСОВ, В.И. ПУНЕГОВ

Отдел математики, Коми НЦ УрО РАН, г. Сыктывкар kolos@dm.komisc.ru, vpunegov@dm.komisc.ru

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

Ключевые слова: динамическая теория дифракции, латерально ограниченный кристалл, рекуррентные соотношения

S.I. KOLOSOV, V.I. PUNEGOV. RECURRENCE RELATIONS IN THE DYNAMICAL THEORY OF X-RAY DIFFRACTION IN A LATERAL CRYSTAL

Darwin's dynamical theory of X-ray diffraction is considered to the cases of planeparallel and lateral (i.e., having a finite length in the lateral direction) crystals. Therefore this new approach can be widely used for non-destructive testing of lateral structures used in opto- and microelectronic devices, nonvolatile memory devices and X-ray optics. We provided a short review of Darwin's approach in the case of plane-parallel crystals and analyzed the exact solution of this problem. The main consideration is devoted to the dynamical X-ray diffraction in the lateral crystal. We restrict ourselves to the case of a symmetrical coplanar Bragg diffraction. Recurrence relations describing the dynamical X-ray diffraction in the lateral crystal is obtained. The angular distribution of X-ray scattering intensity in the reciprocal space is investigated. Phase changes of X-ray waves in the crystal are analyzed in detail. It is shown that this approach allows one to calculate rocking curves as well as Reciprocal Space Maps (RSMs) for lateral crystals having rectangular cross-section. An algorithm for computing of rocking curves and RSMs from the lateral crystal is developed. In accordance with the model, we used a rectangular lattice having a fixed distance between nodes to describe the dynamical diffraction process. The simulation procedure based on recurrence relations consists of the external and internal cycles. The external cycle starts with the first column and goes up to the last column. The internal cycle (from top to bottom) starts with the first reflecting plane and goes up to the last one. Then substituting the amplitudes calculated for the first column we can calculate the amplitudes for the second column and so on. At the end of the external cycle we obtain arrays of the reflected and transmitted amplitudes. We perform numerical modelling for crystals with different lateral sizes. The obtained simulation results are in good agreement with the solution obtained by integration of Takagi's equations. It is demonstrated that the kinematical approximation is valid for thick crystals having small length in the lateral direction. It is shown that this approach, being simple and transparent, is faster than the one based on Takagi's equations and allows simulations of RSMs. The obtained results can potentially be extended further to the three-dimensional case in both the Fourier space and real space.

Keywords: dynamical diffraction theory, lateral crystal, recurrence relations

Введение

В настоящее время существуют разные подходы для описания динамической дифракции рентгеновских лучей в кристаллах [1,2]. Динамическая тео-

рия дифракции Дарвина [3,4], основанная на методе рекуррентных соотношений, является наиболее простой и наглядной для решения конкретных задач рассеяния рентгеновских лучей и нейтронов от кристал-

лов. Продолжительное время теория Дарвина не использовалась, поскольку имела некоторые ограничения по сравнению с теорией Эвальда-Лауэ [1]. Бори [5] первым сделал попытку описать эффект Бормана с использованием системы разностных уравнений Дарвина. В рамках этого подхода Безирганяном и Навасардяном [6] показано, что эффект Бормана зависит от размера отражающих плоскостей. С помощью подхода Дарвина были исследованы важные особенности динамического рассеяния, включая дифракцию в асимметричной геометрии Лауэ[7,8], многоволновое рассеяние [9,10], а также свойства дисперсионной поверхности [7]. Рекуррентные соотношения, описывающие амплитуды проходящей и отраженной волны, применялись для исследования дифракции нейтронов в полиатомных кристаллах [11] и рассеянии света в жидких кристаллах [12]. Кроме того, процедура Дарвина была реализована в оптике слоистых сред [13].

Для исследования приповерхностных областей кристалла разработан формализм, получивший название crystal truncation rod (CTR) [14]. Схема динамической теории Дарвина была применена для исследования структуры поверхности кристалла методом CTR [15-19].

Яширо и Такахаши [20] исследовали коэффициенты отражения и прохождения рентгеновских лучей одной атомной плоскостью кристалла в общем случае для произвольной двумерной решетки Бра-вэ. Ими получены формулы для коэффициентов, отличающиеся от ранее опубликованных результатов Дурбина [21], но подобные формулам Бори [7]. Выражения коэффициентов отражения и прохождения, полученные Яширо и Такахаши [20], которые справедливы в случае произвольных направлений падающего и отраженного рентгеновского пучка, использованы в теории Дарвина для скользящей геометрии [22].

Известно, что теория Дарвина первоначально разработана для совершенного кристалла. Детальный анализ был выполнен для полубесконечного кристалла, включая эффекты преломления и поглощения [23]. Точное решение уравнений Дарвина для кристалла с произвольным числом отражающих плоскостей получено с использованием полиномов Чебышева [24]. Чен и Бхаттачаря [25] показали тождественность алгебраических уравнений Дарвина и дифференциальных уравнений Такаги [26] при переходе от модели массива дискретных атомных плоскостей к континуальной среде.

Отметим, что в природе идеальных кристаллов не существует. Поэтому особое внимание следует уделить работам, в которых формализм рекуррентных соотношений был обобщен для описания дифракции на искаженных структурах. Предприняты попытки для описания методом Дарвина дифракции рентгеновских лучей на кристаллах, решетка которых искажена под влиянием тепловых колебаний атомов [27] или из-за присутствия статистически распределенных дефектов [28]. Существенное влияние на дифракцию рентгеновских лучей оказывают деформации кристаллической решетки. Первая попытка численного моделирования рентгеновской ди-

фракции на кристалле с линейным изменением параметра решетки с использованием уравнений Дарвина была выполнена Фицгеральдом и Дарлингтоном [29]. Результаты численных расчетов кривых дифракционного отражения данной работы согласуются с позже полученными аналитическими решениями [30, 31]. Прудников аналитически обобщил подход Дарвина к описанию рентгеновской дифракции в деформированных кристаллах. В частности, с помощью рекуррентных соотношений исследована дифракция в кристаллах, промодулированных поверхностной акустической волной [32], и в несовершенных гетероструктурах [33].

Рекуррентные соотношения, аналогичные уравнениям Дарвина, получены для описания дифракции в многослойных структурах, при этом коэффициенты отражения и прохождения от атомных плоскостей заменены на соответствующие коэффициенты от отдельных слоев [34, 35]. Такой подход был обобщен на случай многоволновой дифракции [36, 37], резко асимметричной [38] и скользящей геометрии [39].

Все отмеченные выше теории, использующие формулы рекуррентных соотношений Дарвина, разработаны для планарных структур с отражающими атомными плоскостями бесконечных размеров в латеральном направлении. В последние годы заметно возрастает интерес к дифракции рентгеновских лучей на латеральных структурах [40]. Разные методы динамической дифракции в рамках уравнений Такаги [26] были использованы для вычисления кривых дифракционного отражения от совершенных кристаллов прямоугольного сечения [41-43]. Кинематическая теория дифракции применялась для описания дифракции в искаженных кристаллах трапецеидального [44, 45] и произвольного поперечного сечения [46].

Цель данной работы состоит в разработке подхода Дарвина для исследования дифракции рентгеновских лучей в латеральных кристаллах. На примере кристалла прямоугольного сечения будут рассчитаны кривые дифракционного отражения и карты распределения интенсивности рассеяния (Reciprocal Space Maps (RSM)) в зависимости от ширины и толщины кристалла.

В теории Дарвина кристалл рассматривается как совокупность бесконечно тонких атомных плоскостей, отстоящих друг от друга на расстоянии d. В отличие от теории Лауэ [1], в модели Дарвина вся электронная плотность вещества сосредоточена в этих плоскостях. Амплитудные коэффициенты отражения и прохождения отдельной атомной плоскостью вычисляются с использованием теории дифракции Френеля [7, 20]. С другой стороны, переходя от дискретной модели Дарвина к континуальной среде [25] амплитудные коэффициенты отражения и прохождения атомной плоскостью можно выразить через Фурье-коэффициенты рентгеновской поляризуемости xg = -r0\2Fg/(nVc), где Fg - структурный фактор (g = 0, h, h), X - длина волны падающего излучения Vc - объем элементарной ячейки, r0 = e2/(mc2) -классический радиус электрона, где e, m - заряд и масса электрона, соответственно.

В результате коэффициенты, описывающие симметричные отражение и прохождение одной атомной плоскостью в рекуррентных соотношениях Дарвина, запишутся как q0 = ndx0/(Asin0B), q = ndxh/(A sin 0b), q = ndx-h/(A sin 0b), где 0b - угол Брэгга. В этой статье мы рассматриваем случай а-поляризации.

Рассмотрим плоскопараллельный кристалл, имеющий конечную толщину Lz = Nd, где N - число отражающих атомных плоскостей. Пусть на этот кристалл падает плоская рентгеновская волна под углом 0 = 0B + А0, где А0 - малое отклонение от точного брэгговского положения. Согласно теории Дарвина, выражения для амплитуд, прошедшей Tn и отраженной Sn рентгеновских волн для атомной плоскости с номером n, где отсчёт идёт сверху вниз, запишутся в виде следующих рекуррентных соотношений [1,4]

Tn — (1 - iqo)егфTn-i - iqe2гфSn, Sn = (1 - iqo)егфSn+i - iqTn.

(1)

Здесь Tn-i, Sn+i - соответствующие амплитуды волн, относящиеся к плоскостям с номерами n - 1 и n +1, соответственно. Экспоненциальные множители в (1) учитывают набег фазы при распространении рентгеновского поля между соседними атомными плоскостями. Поскольку поле может распространяться под углом, не равным углу Брэгга 0B, то фаза Ф = (2nd/A) sin 0 при распространении рентгеновской волны от одной атомной плоскости до другой получает приращение, зависящее от угловой отстройки. Из-за малости отстройки А0 фаза волны может быть представлена в виде ф = (2nd/A)cos 0BА0. Решение системы (1) с граничными условиями T0 = 1 и SN = 0 описывается следующими аналитическими решениями [47]:

Sn — — В-

uNun - uNun

Tn —

(Aui - 1)uN - (Au2 - 1)uf (Aui - 1)uNNun - (Au2 - 1)uf un

(2)

(Aui - 1)uNN - (Au2 - 1)uf

N

где ui 2 — X ± л/X2 - 1, X —

1 + (1 - iqo)2е2гф + qqe2i*

2(1 - iqo)е1*

А = (1 - iqo)егф, В = —iq. Из решения (2) при значениях п = 0 и п = N получаем амплитудные коэффициенты отражения 50 и прохождения рентгеновских волн от кристаллической пластины

So — -В

TN —

uN

uN

(Aui - 1)uN - (Au2 - 1)uN ui - u2

N

(3)

(Aui - 1)uN - (Au2 - 1)uN

Выражения (3) согласуются сточными аналитическими решениями, полученными ранее с использованием матричного подхода [24]. Учитывая, что для всех п в (2), и™иП = 1, получаем

So —

TN —

iq

1 - (1 - iqo)егф sin([N - 1] arccos X)/ sin(N arccos X)'

(1 - iqo)егф2VX2-1

sin(N arccos X) - (1 - iqo)егф sin([N - 1] arccos X)

Формулы (4) совпадают с решениями, полученными ранее в работе [24]. Введя новый угловой пара-

1 — (1 — iqo )2 е2гф — С[(1е2гф

метр y —

и обозначив х —

2qqeгф

агееоБ ж, ф = агееоБ у, амплитудные коэффициенты отражения и прохождения можно также записать в виде

So —

TN —

q sin Nx q sin(Nx + ф) sin ф

exp(-ip),

(5)

sin(Nx + ф)

exp(-iNy>).

Решения в виде (5) являются компактными и по форме совпадают с решениями для идеальной одномерной сверхрешетки [48].

Рекуррентные соотношения для кристалла прямоугольного сечения

Используем подход Дарвина для описания дифракции от латерально ограниченного кристалла. Обозначим ширину кристалла Ьх, а ее толщину Ьг (рис.1). Выберем следующую систему координат: ось ж направлена вдоль входной поверхности кристалла, ось г - вглубь кристалла. Для простоты рассмотрим симметричную дифракцию Брэгга. Левая грань кристалла, на которую падает рентгеновская волна под углом в к оси ж, проходит через начало координат (рис.1).

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

Рассмотрим некоторый рентгеновский пучок, исходящий из начала координат. Этот пучок проходит расстояние d/ sin 0 и дифракционно отражается от ближайшей атомной плоскости. Проекция этого расстояния на латеральное направление равна Ax — dctg 0. Нанесем на атомные плоскости вдоль оси x разностную сетку с узлами xm — mAx, в которых этот рентгеновский пучок частично отражается и частично проходит к следующей плоскости.

Пусть тт - значение амплитуды проходящей волны и S¡m - соответствующее значение амплитуды отраженной волны непосредственно перед узлом (m; n). Здесь индексы m и n нумеруют координаты узлов в горизонтальном и вертикальном направлениях, соответственно (рис.2). Число узлов решетки вдоль осей X и z, соответственно равное ых и Nz, определяется размерами прямоугольного сечения кристалла Lx — Mx Ax и Lz — Nz d.

Рис. 2. Схематическое изображение отраженных и проходящих рентгеновских пучков в методе рекуррентных соотношений.

С учетом динамического рассеяния для отраженных 5 и проходящих Т волн можно записать следующие рекуррентные соотношения:

T

ш — 1

n+i = aTm + h S.

пШ am — i i i, п^ш — i

Sn = aSn+i + b2 Tn

(6)

Tn+1

где a = (1 — iqo)exp(i0i), Ъ\ = —igexp^i), 62 = -iqехр^ф1). Следует отметить, что мы используем здесь значения для q0, q и g, полученные для бесконечного кристалла, что, безусловно, является приближением. Параметр ф1 = i Л 2nd , входящий в выЛ sin вв

ражения для коэффициентов a и Ъ1?2, учитывает разность фаз, возникающую при распространении рентгеновского пучка в кристалле от одного узла до другого. В системе алгебраических уравнений (6) индексы m и n нумеруют координаты узлов в горизонтальном и вертикальном направлениях, соответственно (рис.2).

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

Рис. 3. Схематическое изображение рентгеновской дифракции в обратном пространстве.

Пусть на латеральный кристалл в направлении волнового вектора к падает рентгеновская волна под углом 01 = вВ + Д01 к поверхности кристалла (рис.3.) Отраженная волна регистрируется в обратном пространстве в направлении вектора к' под углом 02 = 9В + Д02 к поверхности кристалла, причем |к| = |к'| = 2п/Х. Углы Д01 и Д02 задают малые угловые отклонение от точного условия Брэгга падающего и отраженного пучка, соответственно. Отклонение вектора дифракции р = к' - к от конца вектора обратной решетки Ь определяется вектором ц. Проекции qx и qz вектора ц связаны с угловыми параметрами Д01?2 следующими соотношениями

qx = k sin 6b • (Л01 - Л02), qz = -k cos 0b • (Л01 + Л02).

(7)

Обозначив Л01 = и и Л02 = £ - и, соотношения (7) преобразуются в выражения, используемые в трех-кристальной рентгеновской дифрактометрии

qx = k sin 0B • (2и — £),

qz = —k cos 0B •

где и и £ - параметры, определяющие угловое положение образца и кристалла-анализатора [49].

Рассмотрим алгоритм вычислений интенсивности рассеяния в зависимости от значений (qx,qz), или, что одно и то же, от углов 01j2 (см. (7)). Угловые параметры Л01?2, согласно (7), связаны c проекциями вектора q следующими соотношениями: Л01?2 = (±2k cos 0B)-1 (qx ctg 0B T qz).

Граничные условия

При прохождении и отражении рентгеновских лучей в кристалле происходит изменение фаз рентгеновских волн. За начало отсчета фазовых изменений выберем начало системы координат (x = 0; 2 = 0), совмещенное с верхним левым углом прямоугольного сечения кристалла (рис.1). Для рентгеновских волн, падающих на левую боковую грань кристалла

(x = 0), разность хода растет с глубиной z по закону zsin<9i (рис. 4). Разность фаз для дискретных значений zn = nd вдоль вертикального направления равна = (2nd/\)nd sin 01. Граничное условие на левой боковой поверхности латерально ограниченного кристалла (x = о) запишется как тП = exp(^n,m) и #П = 0, где 7i = 0,1,2,...,iVz.

Рис. 4. К выводу фазовых изменений рентгеновских лучей в произвольной точке (x; z) за счет возникновения разности хода Ai = x cos в\ + z sin в\ и А2 =

— x cos в2 + Z sin в2 .

Для рентгеновских лучей, падающих на верхнюю поверхность латерального кристалла (z = 0), разность хода растет с ростом координаты x как mAx cos в!. Следовательно, имеет место изменение фазы падающей рентгеновской волны в горизонтальном направлении в зависимости от угла в!. Граничное условие на верхней поверхности кристалла задается выражением T^ = exp^^n), где ф™ги = (2n/A)mAx cos в!, m = 1, 2,..., Mx.

Для выходящей рентгеновской волны в направлении вектора k' под углом в2 = вв + 6в2 необходимо учитывать разность фаз =

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

—(2n/A)mAx cos в2 на верхней поверхности кристалла (рис.4). Амплитуда рентгеновской волны, выходящей из правой боковой поверхности кристалла, приобретает фазовый набег = —(2n/A)[ndsinв2 — Mx Ax cos в2]. Второе слагаемое в квадратных скобках этого выражения учитывает смещение правой грани кристалла на расстояние Lx = Mx Ax от начала координат.

Поскольку в рассматриваемой геометрии (рис.1) дифракционная волна выходит из верхней и правой боковой грани кристалла, граничные условия для отраженной волны S на нижней и левой боковой гранях кристалла запишутся, соответственно, как S%z =0, 5П = 0.

Амплитудный коэффициент отражения рентгеновской волны от кристалла прямоугольного сечения находится суммированием рентгеновских волн:

S (qx ,qz ) = S (в! ,в2 ) =

Mx

Nz

ЕаШ (• m \ , \ Л aMx (• n \ S0 ехр(1фх,вх) + ¿^ Sn X exP(iPz,ex),

(8)

m=0

n=!

где амплитуды рентгеновских волн Б™ = Б™№ ) и = Б^* (в1,02) находятся с использованием рекуррентных соотношений (6) и зависят от углов падения рентгеновского пучка в1 и угла регистрации в2 отраженной волны. Действительно, для падающей волны на границах кристалла амплитуды Т0 = Т0(в1) и

Tm = тт(в!) зависят только от угла падения. Однако в силу того, что фаза отраженной волны фПех = фПех (в2) зависит от угла регистрации пучка, то в процедуре вычислений в рамках рекуррентных соотношений (6) амплитуды волн в направлении прохождения тт = тт (в! ,в2) и отражения Sm = Sm (в! ,в2) будут зависеть уже от двух углов или вертикальной и латеральной проекций вектора q. Поэтому решение (8) позволяет рассчитывать карты распределения интенсивности рассеяния в обратном пространстве (RSM). Кривые дифракционного отражения (qx-сканы или CTR) могут непосредственно получены на основе решения (8) при условии в! = в2 (2и = г).

Численное моделирование

Численное моделирование кривых дифракционного (КДО, CTR) и RSM проводилось на основе решения (8) с использованием рекуррентных соотношений (б). Расчеты выполнены с использованием параметров, соответствующих (111) отражению CuKai излучения от монокристалла германия. Для всех вычислений толщина латерального кристалла Lz соответствовала 10000 межплоскостных расстояний (Lz = 10000dm = 3,27{ш).

Латеральный размер кристалла выбирался с учетом числа узлов Nx вдоль горизонтального направления (Nx = 1000, 4000, 10000, 40000; что соответствует латеральным размерам Lx =

1,35 {¿ш, 5,39 {¿ш, 13,5 {¿ш, 53,9 дш,). Максимумы кривых дифракционного отражения нормировались на величину максимума кривой Дарвина для плоскопараллельной пластинки той же толщины.

Кривые дифракционного отражения, изображенные на рис.5, вычислялись на основе решения (8) для случая и — 2в-сканирования (в! = в2). Поскольку в этом режиме сканирования полная разность хода A = A! + A2 = 2z sin в (в = в! = в2) (рис.4) не зависит от координаты x, то и полная разность фаз на верхней поверхности кристалла равна нулю: ^m = m + (pm,ex = 0. Отметим, что длина экстинк-ции на отражении (111) для полубесконечного кристалла кремния равна 0,67 дт, полуширина дарвиновской кривой - 15,4 угл. сек. (см. [50] и адрес сайта "X-Ray Server").

Для малой латеральной ширины кристалла (Lx = 1,35 дт) кривая отражения соответствует кинематической дифракции (рис.5а). С увеличением латеральной ширины наблюдается постепенный переход к динамической дифракции (рис.5Ь,с), при этом толщинные осцилляции накладываются на профиль КДО даже в области полного отражения.

На рис.6. представлены RSM от кристаллов разной ширины. В случае кристалла с малым латеральным размером RSM имеет вид, отвечающий кинематической дифракции (рис.ба), и распределение дифракционной интенсивности соответствует известному закону [1]

I (qx ,qz) = S (qx ,qz) rc

sin qx Lx/2 2 sin qz Lz /2

qx Lx/2 qz Lz/2

Рис. 5. Кривые дифракционного отражения для кристалла прямоугольного сечения толщиной Ь2 = 3,27 дш и разной ширины Ьх: а) 1,35 дт; б) 5,39 дт; в) 13,5 дт; г) 53,9 дт. Для сравнения тонкой красной линией 2 изображена кривая дифракционного отражения от плоскопараллельной пластинки (Ьх = то) той же толщины Ь2.

Рис. 6. RSM от кристалла прямоугольного сечения толщиной Ь2 = 3,27дш и разной ширины Ьх: а) 1,35 дт; б) 5,39 дт; в) 13,5 дт; г) 53,9 дт.

Рис. 7. Кривые качания ^х-сканирования) от кристалла прямоугольного сечения толщиной Ьг = 3,27 дш и разной ширины Ьх: а) 1,35 дт; б) 5,39 дт; в) 13,5 дт; г) 53,9 дт.

С увеличением латерального размера кристалла угловое распределение интенсивности рассеяния на RSM в qx-направлении становится все более узким, уменьшается также период латеральных осцилляций. Этот период обратно пропорционален

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

Результаты полностью согласуются с теми, что получены решением данной задачи путем интегрирования уравнений Такаги-Топена.

Заключение

Таким образом, мы показали, что теория Дарвина, в основе которой лежат алгебраические рекуррентные соотношения, в геометрии Брэгга может описывать не только дифракцию на плоскопараллельных пластинках, но и на латерально ограниченных кристаллах. Рассмотренный подход является не только наиболее простым и наглядным, но и требует существенно меньших временных затрат на вычисления кривых дифракционного отражения и RSMs. Поскольку в настоящее время латеральные структуры, создаваемые методами литографии и электрохимического травления, широко применяются в приборах микро- и оптоэлектроники, энергонезависимых запоминающих устройствах (nonvolatile memory devices), элементах рентгеновской оптики и т.д., разработанный подход может быть эффективно использован для неразрушающего рентгенодифрак-ционного анализа таких материалов.

Работа выполнена при частичной финансовой поддержке РФФИ (проект № 13-02-00272-а).

Литература

1. Authier A. Dynamical Theory of X-Ray Diffraction. New York: Oxford University Press, 2001. 661 p.

2. Pietsch U, Holy V., and Baumbach T. High resolution X-ray scattering: from thin films to lateral nanostructures. 2nd edition. New York: Springer-Verlag, 2004.

3. Darwin C.G. The Theory of X-Ray Reflexion // Philos. Mag. 1914. Vol. 27. P. 315-333.

4. Darwin C.G. The Theory of X-Ray Reflexion. Part II // Philos. Mag. 1914. Vol. 27. P. 675691.

5. Borie B. The Theory of the Borrmann Effect in Terms of Difference Equations // Acta Cryst. 1966. Vol. 21. P. 470-472.

6. Безирганян ПА., Навасардян MA. Теория эффекта Бормана для конечного кристалла // Изв. АН Арм. ССР. Физика. 1969. Т. 3. С. 269274.

7. Borie B. The Darwin dynamical theory of x-ray diffraction //Acta Cryst. 1967. Vol. 23. P. 210216.

8. Кузнецов A.B., Фофанов А.Д. Дарвиновская теория рассеяния идеальными кристаллами в асимметричном положении // Известия ВУЗов. Физика. 1970. № 10. C. 12-18.

9. Кузнецов A.B., Фофанов А.Д. О тождественности методов Дарвина и Эвальда-Лауэ (обобщение метода Дарвина на любое число сильных волн) // Известия ВУЗов. Физика. 1972. № 4. C. 108-113.

10. Игнатович В.К. Многоволновой алгебраический метод Дарвина в динамической теории дифракции // Кристаллография. 1992. Т. 37. № 5. С. 1100-1112.

11. Игнатович В.К. Алгебраическое описание динамической дифракции на многоатомных идеальных кристаллах // ЖЭТФ. 1990. Т. 97. № 5. C. 1616-1625.

12. Chandrasekhar S, Srinivasa Rao K. N. Optical rotatory power of Liquid crystals // Acta Cryst. 1968. Vol. A24. P. 445-451.

13. Dub P., Litzman O. The Darwin procedure in optics of layered media and the matrix theory // Acta Cryst. 1999. Vol. A55. P. 613-620.

14. Robinson I.K. Crystal truncation rods and surface roughness // Phys. Rev. B. 1986. Vol. 33. No. 6. P. 3830-3836.

15. Caticha A. Diffraction of x rays at the far tails of the Bragg peaks. II. Darwin dynamical theory // Phys. Rev. B. 1994. Vol. 49. No. 1. P. 3338.

16. Durbin S. M., Follis G.C. Darwin theory of heterostructure diffraction // Phys. Rev. B. 1995. Vol. 41. №15. P. 10127-10133.

17. Nakatani S., Takahashi T. Dynamical treatment of x-ray reflection from crystal surface // Sur. Sci. 1994. Vol. 311. P. 433-439.

18. Takahashi T., Nakatani S. Dynamical theory of X-ray diffraction for the study of crystal surfaces// Surface science, 1995. Vol. 326. P.347-360

19. Takahashi T., Yashiro W., Takahasi M., Kusano S. Effect of surface structure on crystal-truncation-rod scattering under the Bragg condition // Physical Review B. 2000. Vol. 62. P. 3630-3638.

20. Yashiro W., Takahashi T. X -ray diffraction from an atomic plan // Acta Cryst. 2000. Vol. A56. P. 163-167.

21. Durbin S. M. Darwin spherical-wave theory of kinematic surface diffraction // Acta Cryst. 1995. Vol. A51. P. 258-268.

22. Yashiro W., Ito Y., Takahasi M, Takahashi T. Darwin's theory for the grazing incidence geometry // Surface science. 2001. Vol. 490. P. 394-408.

23. Prins J. A. Die Reflexion von Rontgenstrahlen an absorbierenden idealen Kristallen // Zeitschrift fur Physik. 1930. Vol. 63. No. 7-8. P. 477-493.

24. Perkins R.T., Knight L.V. An exact analytic solution of Darwin's difference equations // Acta Cryst. 1984. Vol. A40. P. 617-619.

25. Chen, Y.C. Bhattacharya, P.K. Determination of critical layer thickness and strain tensor in InxGai-xAs/GaAs quantum-well structures by x-ray diffraction // Journal of Applied Physics. 1993. Vol. 73. No. 11. P. 7389-7394.

26. Takagi S. A Dynamical theory of diffraction applicable to crystals with any kind of small distortion // Acta Cryst. 1962. Vol. 15. № 12. P. 1311-1312.

27. Chung J.S., Durbin S.M. Temperature-dependent

X-ray dynamical diffraction: Darwin theory simulations // Acta Cryst. (1999). Vol. A55. P. 14-19.

28. Li M, Ress H, Gerhard T., Landwehr G., Cui S.F., Mai Z.H. A semi-kinematic approach to x-ray diffraction of real crystals with small defects // J. Appl. Phys. 1997. Vol. 81. P. 2143-2147.

29. Fitzgerald W.J., Darlington C.N.W. A numerical study of the dynamical theory of scattering from a distorted crystal // Acta Cryst. 1976. Vol. A32. P. 671-672.

30. Kolpakov A.V., Punegov V.I. X-ray diffraction in a crystal with a linear lattice parameter variation // Solid State Communications. 1985. Vol. 54. No. 7. P.573-578.

31. Punegov V.I., Vishnjakov Yu.V. Dynamical theory of X-ray diffraction from damaged epitaxial layers with constant strain gradient //J. Phys. D: Appl. Phys. 1995. Vol. 28. No. 4A. P. A184-A188.

32. Prudnikov I.R. Darwin's approach to X-ray dynamical diffraction in distorted crystals // Acta Cryst. 1998. Vol. A54. P. 1034-1036.

33. Prudnikov I.R. X-Ray Diffraction in Distorted Crystals and Imperfect Multilayered Structures. A Theoretical Approach Based on Resursion Equations // Phys. Stat. Sol. (b). 2000. Vol. 217. P. 725-735.

34. Bartels W.J., Hornstra J., Lobeek D.J.W. Acta Crystallogr., Sect. A: Found. Crystallogr. 42, 539 (1986).

35. Belyaev Yu.N., Kolpakov A.V. On the theory of x-ray diffraction in a perfect crystal with distorted surface layer // Phys.Stat.Sol.(a). 1983. Vol. 76. P. 641-645.

36. Ladanov A.V., Punegov V.I. A new approach to the problem of multiple X-ray diffraction in multilayer crystals // Twelfth Europen Crystallographic meeting. Moscow, USSR, August 20-29, 1989 / Collected Abstracts. Moscow. 1989. Vol. 3. P. 137.

37. Punegov V.I. X-ray diffraction from multilayer structures with statistically distributed microdefects. // Physica status solid (a). 1993. Vol. 136. No. 1. P. 9-19.

38. Пунегов В.И., Ладанов A.B. Теория дифракции рентгеновских лучей в многослойных кристаллах в условиях полного внешнего отражения // Поверхность. Физика. Химия. Механика. 1990. № 4. С. 45-50.

39. Пунегов В.И., Ладанов A.B. О дифракции рентгеновских лучей на многослойной пленочной структуре в скользящей геометрии // Журнал технической физики. 1989. Т. 59. № 11. С. 188-189.

40. Lee K, Yi H., Park W.-H, Kim Y.K. Baik S. Lateral size effects on domain structure in epitaxial PbTiO3 thin films // J. Appl. Phys. 2006. Vol. 100. P. 051615 (1-10).

41. Колосов С.И., Пунегов В.И. Методы численного интегрирования уравнений Та-каги-Топена для кристалла прямоугольного сечения // Кристаллография. 2005.

Т. 50. Вып. 3. С. 401-406.

42. Olekhnovich N.M., Olekhnovich A.I. Primary Extinction for Finite Crystals. Square-Section Parallelepiped // Acta Cryst. 1978. Vol. A34. № 2. P.321-326.

43. Thorkildsen G, Larsen H.B. X-ray diffraction in perfect crystals. Rocking curves // Acta Cryst. 1999. Vol. A55. P. 840-854.

44. Пунегов В.И., Колосов С.И., Павлов KM. К теории дифракции рентгеновских лучей на латеральном кристалле с упруго изогнутыми атомными плоскостями // Письма в ЖТФ. 2006. T. 32. Вып. 18. C. 65-72.

45. Пунегов В.И., Колосов С.И. Теория дифракции рентгеновских лучей на неидеальном кристалле трапецеидального сечения // Кристаллография. 2007. Т. 52. № 2. С. 215-222.

46. Пунегов В.И., Максимов А.И., Колосов С.И., Павлов К.М. К расчету дифракции рентгеновских лучей от многослойных латеральных кристаллических структур произвольных композиционного состава и формы // Письма в ЖТФ. 2007. Т. 33. Вып. 3. C. 64-71.

47. Пунегов В.И. Рентгеновская дифракция на полупроводниковой сверхрешетке с микродефектами // Письма в ЖТФ. 1992. Т. 18. № 4. С. 66-70.

48. Vardanyan D.M., Manoukyan H.M., Petrosyan H.M. The dynamic theory of x-ray diffraction by the one-dimensional ideal superlattice. I. Diffraction by the arbitrary superlattice // Acta Cryst. 1985. Vol. A41. P. 212-217.

49. Iida A., Kohra K. Separate measurements of dynamical and kinematical x-ray diffractions from silicon crystals with a triple crystal diffractometer // Phys. Stat. Sol. (a). 1979. Vol. 51. P. 533-542.

50. Stepanov S, Forrest R. Fitting dynamical X-ray diffraction data over the World Wide Web // Journal of Applied Crystallography. 2008. Vol. 41. P. 958962, http://sergey.gmca.aps.anl.gov/.

References

1. Authier A. Dynamical Theory of X-Ray Diffraction. New York: Oxford University Press, 2001. 661 p.

2. Pietsch U, Holy V., and Baumbach T. High resolution X-ray scattering: from thin films to lateral nanostructures, 2nd edition. New York: Springer-Verlag, 2004.

3. Darwin C.G. The Theory of X-Ray Reflexion // Philos. Mag. 1914. Vol. 27. P. 315-333.

4. Darwin C.G. The Theory of X-Ray Reflexion. Part II // Philos. Mag. 1914. Vol. 27. P. 675691.

5. Borie B. The Theory of the Borrmann Effect in Terms of Difference Equations // Acta Cryst. 1966. Vol. 21. P. 470-472.

6. Bezirganian P.H., Navasardian MA. Teoria ef-fekta Bormana dlja konechnogo kristalla [The theory of the Borrmann effect for the finite

crystal] // Izv. AN Arm. CCR, Fizika (1969) Vol. 3. P. 269-274.

7. Borie B. The Darwin dynamical theory of x-ray diffraction //Acta Cryst. 1967. Vol. 23. P. 210216.

8. Kuznetsov A.V., Fofanov A.D. Darvinovskaja teorija rassejanija ideal'nymi kristallami v asim-metrichnom polozhenii [The Darwin theory for scattering by an ideal crystal in the asymmetric case] // Izvestia VUZov, Fizika. 1970. No. 10. P. 12-18.

9. Kuznetsov A.V., Fofanov A.D. O tozhdestvennosti metodov Darvina i Jeval'da-Lauje (obobshhenie metoda Darvina na ljuboe chislo sil'nyh voln) [On the equivalence of methods of Darwin and Ewald-Laue] // Izvestia VUZov, Fizika. 1972. No. 4. P. 108-113.

10. Ignatovich V.K. Mnogovolnovoj algebraicheskij metod Darwina v dinamicheskoj teorii difrakcii [The multiwave algebraic method of Darwin in the dynamic theory of diffraction] // Kristallo-grafija. 1992. Vol. 37. No. 5. P. 1100-1112.

11. Ignatovich V.K. Algebraicheskoe opisanie di-namicheskoi difrakcii na mnogoatomnyh ide-alnyh kristallah [Algebraic description of dynamic diffraction by ideal polyatomic crystals] // JETF. 1990. Vol. 97. No. 5. P. 1616-1625.

12. Chandrasekhar S, Srinivasa Rao K.N. Optical rotatory power of Liquid crystals // Acta Cryst. Vol. A24. 1968. 445-451.

13. Dub P., Litzman O. The Darwin procedure in optics of layered media and the matrix theory // Acta Cryst. 1999. Vol. A55. P. 613-620.

14. Robinson I.K. Crystal truncation rods and surface roughness // Phys. Rev. B. 1986. Vol. 33. No. 6. P. 3830-3836.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

15. Caticha A. Diffraction of x rays at the far tails of the Bragg peaks. II. Darwin dynamical theory // Phys. Rev. B. 1994. Vol. 49. No. 1. P. 3338.

16. Durbin S.M., Follis G.C. Darwin theory of het-erostructure diffraction // Phys. Rev. B. 1995. Vol. 41, No. 15. P. 10127-10133.

17. Nakatani S., Takahashi T. Dynamical treatment of x-ray reflection from crystal surface // Sur. Sci. 1994. Vol. 311. P. 433-439.

18. Takahashi T., Nakatani S. Dynamical theory of X-ray diffraction for the study of crystal surfaces // Surface science. 1995. Vol. 326. P. 347360.

19. Takahashi T., Yashiro W., Takahasi M., Kusano S. Effect of surface structure on crystal-truncation-rod scattering under the Bragg condition // Phys. Rev. B. 2000. Vol. 62. P. 3630-3638.

20. Yashiro W., Takahashi T. X-ray diffraction from an atomic plan // Acta Cryst. 2000. Vol. A56. P. 163-167.

21. Durbin S.M. Darwin spherical-wave theory of kinematic surface diffraction // Acta Cryst. 1995. Vol. A51. 258-268.

22. Yashiro W., Ito Y., Takahasi M, Takahashi T. Darwin's theory for the grazing incidence geometry // Surface science. 2001. Vol. 490. P. 394-

408.

23. Prins JA. Die Reflexion von Röntgenstrahlen an absorbierenden idealen Kristallen // Zeitschrift fur Physik. 1930. Vol. 63. No. 7-8. P. 477-493.

24. Perkins R.T., Knight L.V. An exact analytic solution of Darwin's difference equations // Acta Cryst. 1984. Vol. A40. P. 617-619.

25. Chen Y.C., Bhattacharya P.K. Determination of critical layer thickness and strain tensor in InxGai-xAs/GaAs quantum-well structures by x-ray diffraction // J. Appl. Phys. 1993. Vol. 73. No. 11. P. 7389-7394.

26. Takagi S. A Dynamical theory of diffraction applicable to crystals with any kind of small distortion // Acta Cryst. 1962. Vol. 15. No. 12. P. 1311-1312.

27. Chung J.S., Durbin S.M. Temperature-dependent X-ray dynamical diffraction: Darwin theory simulations // Acta Cryst. 1999. Vol. A55, 14-19.

28. Li M., Ress H., Gerhard T., Landwehr G., Cui S.F., Mai Z.H. A semi-kinematic approach to x-ray diffraction of real crystals with small defects // J. Appl. Phys. 1997. Vol. 81. P. 21432147.

29. Fitzgerald W.J., Darlington C.N.W. A numerical study of the dynamical theory of scattering from a distorted crystal // Acta Cryst. 1976. Vol. A32. P. 671-672.

30. Kolpakov A.V., Punegov V.I. X-ray diffraction in a crystal with a linear lattice parameter variation // Solid State Communications. 1985. Vol. 54. No. 7. P. 573-578.

31. Punegov V.I., Vishnjakov Yu.V. Dynamical theory of X-ray diffraction from damaged epitaxial layers with constant strain gradient //J. Phys. D: Appl. Phys. 1995. Vol. 28. No. 4A. P. A184-A188.

32. Prudnikov I.R. Darwin's approach to X-ray dynamical diffraction in distorted crystals // Acta Cryst. 1998. Vol. A54. P. 1034-1036.

33. Prudnikov I.R. X-Ray Diffraction in Distorted Crystals and Imperfect Multilayered Structures. A Theoretical Approach Based on Resursion Equations // Phys. Stat. Sol. (b). 2000. Vol. 217. P. 725-735.

34. Bartels W.J., Hornstra J., Lobeek D.J.W. X-ray Diffraction of Multilayers and Superlattices // Acta Crystallogr., Sect. A: Found. Crystallogr. 1986. Vol. 42, P. 539.

35. Belyaev Yu.N, Kolpakov A.V. On the theory of x-ray diffraction in a perfect crystal with distorted surface layer // Phys.Stat.Sol.(a). 1983. Vol. 76. P. 641-645.

36. Ladanov A.V., Punegov V.I. A new approach to the problem of multiple X-ray diffraction in multilayer crystals. // Twelfth Europen Crystal-lographic meeting.- Moscow, USSR, August 2029,1989. / Collected Abstracts. Moscow. 1989. Vol. 3. P. 137.

37. Punegov V.I. X-ray diffraction from multilayer structures with statistically distributed microdefects // Physica status solidi (a). 1993.

Vol. 136. No. 1. P. 9-19.

38. Punegov V.I., Ladanov A.V. Teoria difrakcii rentgenovskih luchei v mnogosloinyh kristal-lah v usloviah polnogo vneshnego otrajenija [The theory of X-ray diffraction in multilayer crystals under total external reflection] // Poverhnost. Fizika. Himija. Mehanika. 1990. No. 4. P. 45-50.

39. Punegov V.I., Ladanov A.V. O difrakcii rentgen-ovskih luchey na mnogosloinoi plenochnoi struk-ture v skolzashei geometrii [X-ray diffraction under specular reflection conditions on multilayers] // J. techn. fiz. 1989. Vol. 59. No. 11. P.188-189.

40. Lee K, Yi H, Park W.-H, Kim Y.K. Baik S. Lateral size effects on domain structure in epitaxial PbTiO3 thin films // J. Appl. Phys. 2006. Vol. 100. P. 051615 (1-10)

41. Kolosov S.I., Punegov V.I. Metody chislennogo integrirovanija uravnenii Takagi-Topena dla kristalla pramougolnogo sechenia [Numerical integration methods for the Takagi-Taupin equations for crystals of rectangular cross section] // Kristallografija. 2005. Vol. 50. Vyp. 3. P. 401406.

42. Olekhnovich N.M., Olekhnovich A.I. Primary Extinction for Finite Crystals. Square-Section Parallelepiped // Acta Cryst. 1978. Vol. A34. No. 2. P. 321-326.

43. Thorkildsen G, Larsen H.B. X-ray diffraction in perfect crystals. Rocking curves // Acta Cryst. 1999. Vol. A55. P. 840-854.

44. Punegov V.I., Kolosov S.I., Pavlov K.M. K teorii difrakcii rentgenovskih luchej na lat-eral'nom kristalle s uprugo izognutymi atom-nymi ploskostjami [Theory of X-Ray Diffraction on a Lateral Crystal with Elastically Bent Lattice Planes] // Pis'ma v JTF. 2006. Vol. 32. Vyp. 18. P. 65-72.

45. Punegov V.I., Kolosov S.I. Theory of x-ray diffraction from a nonideal crystal with a trapezoidal cross section // Crystallography Reports. 2007. Vol. 52. No. 2. P. 191-198.

46. Punegov V.I, Maksimov A.I., Kolosov S.I., Pavlov K.M. K raschetu difrakcii rentgenovskih luchej ot mnogoslojnyh lateral'nyh kristallicheskih struktur proizvol'nyh kompozicionnogo sostava i formy [Calculating x-ray diffraction from multilayer lateral structures with arbitrary shapes and composition profiles] // Pis'ma v JTF. 2007. Vol. 33. Vyp. 3. No. 2. P. 125-127.

47. Punegov V.I. Rentgenovskaja difrakcija na poluprovodnikovoj sverhreshetke s mikrodefek-tami [X-ray diffraction by semiconductor superlattice with microdefects] // Pis'ma v JTF. 1992. Vol. 18. No. 4. P. 66-70.

48. Vardanyan D.M., Manoukyan H.M., Petrosyan H.M. The dynamic theory of x-ray diffraction by the one-dimensional ideal superlattice. I. Diffraction by the arbitrary superlattice // Acta Cryst. 1985. Vol. A41. P. 212-217.

49. Iida A., Kohra K. Separate measurements of dynamical and kinematical x-ray diffractions from

silicon crystals with a triple crystal diffractome-ter // Phys. Stat. Sol. (a). 1979. Vol. 51. P. 533542.

50. Stepanov S, Forrest R. Fitting dynamical X-ray

diffraction data over the World Wide Web //J. Appl. Crystallography. 2008. Vol. 41. P. 958-962.

Статья поступила в редакцию 06.04.2015.

i Надоели баннеры? Вы всегда можете отключить рекламу.