Постановка и решение коллимационной задачи при малоугловом рентгеновском рассеянии на анизотропных объектах Текст научной статьи по специальности «Математика»

УДК 517.968, 539.261

ПОСТАНОВКА И РЕШЕНИЕ КОЛЛИМАЦИОННОЙ ЗАДАЧИ ПРИ МАЛОУГЛОВОМ РЕНТГЕНОВСКОМ РАССЕЯНИИ НА АНИЗОТРОПНЫХ ОБЪЕКТАХ

Д. Д. Захаров1, А. В. Смирнов1, Б. А. Федоров1

1 Санкт-Петербургский государственный университет информационных технологий,

механики и оптики

smirnav@phd.ifmo.ru

PACS 61.05.C-, 02.30.Rz

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

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

Ключевые слова: малоугловое рентгеновское рассеяние, анизотропные объекты, метод базисных функций, итерационный метод Фридмана, высоко ориентированный полиэтилен.

1. Введение

Анализ и учет коллимационных искажений является одной из наиболее важных и сложных задач в методике малоуглового рентгеновского диффузного рассеяния [1,2]. Суть проблемы сводится к тому, что на практике для повышения светосилы рентгеновской установки используют, как правило, щелевые коллиматоры и щелевые приемники, в то время как почти вся развитая теория рентгеновского рассеяния основывается на точечном источнике и точечном приемнике. Очевидно, что использование щелей приводит к определенному искажению экспериментальных интенсивностей рассеяния, в результате которого затруднен их анализ с точки зрения общей теории дифракции. Возникает «коллимационная» задача — численный пересчет экспериментальных данных, полученных с использованием щелей («искаженная» интенсивность рассеяния), на данные, которые были бы получены при использовании точечного источника и точечного приемника («точечная» интенсивность рассеяния).

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

анизотропные свойства. В случае изотропного объекта его интенсивность рассеяния зависит только от угла рассеяния и является функцией одной переменной. Для изотропных объектов в различный лабораториях мира (в том числе, и в нашей) разработана целая серия алгоритмов введения коллимационных поправок [3-5], позволяющих с высокой точностью, доходящей до нескольких процентов [6, 7], осуществлять пересчет искаженной интенсивности на точечную.

В случае анизотропного объекта его интенсивность рассеяния зависит не только от угла рассеяния, но и от ориентации объекта в пространстве. В практическом отношении наиболее важным является случай рассеяния от анизотропных пленок, расположенных перпендикулярно первичному пучку. При этом интенсивность рассеянного излучения зависит от двух параметров: угла рассеяния и угла поворота пленки относительно коллимационной щели. Очевидно, задача коллимационного пересчета на точечную интенсивность в этом случае оказывается несравненно более сложной. К настоящему времени этот вопрос обсуждался только в работе [8], однако достаточно корректного алгоритма для введения коллимационных поправок описано не было.

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

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

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

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

2. Математическая постановка задачи

Пусть J = J(q, ф) — искаженная (экспериментальная) интенсивность рассеяния плоским анизотропным образцом с осью симметрии, лежащей в плоскости образца, I = I(q, ф) — искомая точечная интенсивность рассеяния того же образца. Здесь q — модуль вектора рассеяния, р — угол поворота образца относительно коллимационной щели. В дальнейшем в качестве первого аргумента обеих интенсивностей будем использовать угол рассеяния 9, связанный с модулем вектора рассеяния соотношением q = (4ж/Х) sin (0/2), где Л — длина волны рентгеновских лучей.

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

которого 2к1. Образец расположен перпендикулярно падающему излучению. Рассеянные лучи, прошедшие через образец, регистрируются щелью детектора длиной 2к2, расположенной под углом ф к первичному пучку (ось поворота совпадает с осью ув плоскости образца). Использование щелей повышает светосилу установки, но регистрируемая щелью детектора интенсивность ,1 будет отличаться от I, причем тем сильнее, чем больше длины 2к1 g и 2к2.

РИС. 1. Схема хода лучей при щелевой коллимации

/и

д(у) ¿у. (1)

Связь между интенсивностями 3 и I выражается в общем случае интегральным соотношением

У I(в, ф) и(у') йу'

В этом уравнении (см. рис. 1) Ь —расстояние между плоскостью образца и плоскостью регистрации; у — координата, отсчитываемая вдоль проекции первичного пучка в плоскости регистрации; у' — координата, отсчитываемая вдоль детекторной щели; 9 — фактический угол рассеяния; ф = ф + а — угол между плоскостью рассеяния и осью симметрии образца; а — угол между осью у и плоскостью рассеяния; д(х) — распределение интенсивности вдоль проекции первичного пучка в плоскости регистрации; и(у') —распределение чувствитель-

т? X тт( Л Г и(у'), при ¡у' \ ^

ности вдоль детекторной щели. Если ввести функции и (у ) = < ^ при ¡у' | > к и

\ \ я(у), при ¡у1 ^ к1 „ / , ч ,т

^(У) = | о при \у\ > к , а также перейти к новым переменным х = (У — У) /Ьи

( = у/Ь,то уравнение (1) можно привести к виду

—те

1{в,ф)= [ I ' \2- агс1у.;\ 0) • г:) II ( \ ) \ , (2)

где Ш (х)- «весовая функция», определяемая соотношением

н

^ (х) = У и ((х + С) Ь) С ((Ь) ¿С , при Н = (Иг + И2) /Ь, (3)

-н

и зависящая только от параметров рентгеновской установки.

Уравнение (2) представляет собой интегральное уравнение I рода типа Вольтерра, записанное в нестандартной форме, поскольку в нем нет в явном виде ядра, и искомая функция I зависит не от одного внутреннего аргумента х? а от комбинаций аргументов: л/в2 + х2 и аг<^(£/0) + <р.

Задача состоит в нахождении точечной интенсивности Iиз уравнения (2) по экспериментальной интенсивности ,1, измеренной в области значений аргументов [9тп,9тах] х [<£тт, <Апах]. Эта задача является некорректной, поскольку одной функции рассеяния 3 формально соответствует бесконечное множество решений I уравнения (2) [9]. Среди них необходимо найти решение, достаточно близкое к действительной точечной интенсивности. Поиску такого решения и посвящена настоящая работа.

В работе будет использована следующая терминология:

I(9,ф) —точечная интенсивность рассеяния, рассчитанная для выбранной анизотропной модели в условиях точечного источника и точечного приемника рентгеновского излучения;

,1 (9, ф) — искаженная интенсивность рассеяния, рассчитанная по формуле (2) с учетом параметров используемой коллимационной системы;

I'(9,ф) — восстановленная точечная интенсивность рассеяния, рассчитанная на основе развитых методов (метода базисных функций и модифицированного метода Фридмана);

,1'(9,ф) —восстановленная искаженная интенсивность рассеяния, рассчитанная по формуле (2) с использованием I'(в, ф) как новой точечной интенсивности.

В случае обработки экспериментальных данных для ,1 (в, ф) используется термин «экспериментальная» интенсивность; для ,1'(в, ф) — восстановленная «экспериментальная» интенсивность.

3. Метод базисных функций

Введем || — норму в пространстве квадратично интегрируемых функций двух пе-

/ ^тах Отах

ременных в и <р: Ц/1| = Ф / / /2(в, ф) ¿9 ¿<р. С использованием данной нормы относи-

У ^тт $тт

тельное отклонение функции Д от функции /2 можно определить как Пусть задано уравнение 3 = А(1), где интегральный оператор

те

А(1)= I I (л/92 + Х2,\агс1д(Х/9) + ф\)\УШх • (4)

—те

Выберем базисные функции (г = 1, 2,...,п) так, чтобы их линейная комбинация I' =

п

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

г=1

существовать такие коэффициенты сг, с2, ..., сп, для которых относительное отклонение

\\1-1'\\

——— является достаточно малой величинои.

Рассмотрим уравнение ,1' = А(Г). Исходя из линейности оператора А, ,1' = ^ с^,

г=1

где Рг = А (/г) (г = 1, 2,...,п). В поставленной задаче I является искомой функцией, и представить ее разложение по базисным функциям в явном виде нельзя. Однако можно найти для экспериментальной интенсивности ,1 ее приближение как разложение ,1' по функциям При корректном выборе функций можно найти коэффициенты с^, для

которых -————- будет достаточно малым. Эти коэффициенты будут использованы для

построения восстановленной точечной интенсивности I' = ^ сг1г.

г=1

\\1 - I' II

Очевидно, что если относительное отклонение ——— невелико, то невелико будет

\\J - ^ II

и относительное отклонение ——-—, поскольку оператор А в силу его линеиности непре-

рывен. Однако из-за некорректности задачи обратное утверждение несправедливо: одной экспериментальной функции ,1 соответствует множество решений I, для которых выполняется уравнение ,1 = А(1). Чтобы сузить класс решений, следует: а) выбрать базисные функции, наилучшим образом соответствующие данной задаче; б) ввести дополнительные ограничения на решение, используя его известные свойства.

В данной работе двухмерные базисные функции ¡г были построены как прямое произведение одномерных кубических В-сплайнов [10].

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

М / п \ 2

Д(С1, С2, . .., сп) = ^ I ,1 (9к, ) — ^ СгРг(вк, Рк) . (5)

к=1 \ г=1 )

Здесь вк—значения угла рассеяния и угла поворота образца, для которых получена экспериментальная интенсивность, М — общее число экспериментальных точек. В дальнейшем функцию Д будем называть «оценочной функцией».

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

М / п \ 2 п

Ai(ci, С2, сп) = ^ IJ(9к,фк) CiFi(®к)) + °2 (6)

к=1 \ г=1 J г=1

Здесь Х1 — параметр, регулирующий влияние введенного ограничения: слишком большие значения А1 приведут к почти нулевому решению (все q ~ 0), а слишком маленькие А1 не дадут желаемого эффекта.

Второе введенное ограничение связано с тем, что для нулевого угла рассеяния (в = 0) точечная интенсивность должна оставаться постоянной при изменении угла ц>, то есть I(0,р) = const. Если при проведении эксперимента минимальный угол рассеяния достаточно мал, то должно выполняться условие I(6min,p) ~ const. Это условие можно

учесть, если ввести в функцию оценки дополнительное слагаемое

^Е

кек

I п \ п

(®к, Рк М - ^ Су{9к, фк) \з = 1 / К 3=1

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

(^2 °3¡3(вк

1

^2^2 С31з (в к ).

к Л кек з=1

После алгебраических преобразований получаем окончательную функцию оценки в следующем виде:

2

Д2(С1, С2, ..., Сп)

М /

Е I

к=1 \

г=1

<1 (о к ,фк) — с^(вк +

+ А^ С2 +

г=1

кек

(вк )

.3=1

(7)

где введено обозначение

1

<Рк) = <Рк))К - ¡]{0к, <Рк) = У^ /){°к, <Рк) - ^{Ок, <Рк)-

ЦКик,фк) - д— / ,3з\ик,

К кек

Коэффициенты разложения Сг можно найти из условий системе линейных уравнений

дА2

дСг

0, которые преобразуются к

Е

3 = 1

(

м

мЕ

^=1

, >£к)Рэ (9к, фк) + ^2^2 ^ (®к, Vк )^з (®к, ^к)

кек

+ А

м

(8)

= > .1 (9к)Щ9к,Рк); г = 1, 2,

к=1

Главная задача при решении основного интегрального уравнения (2) состоит в подборе оптимальных значений регулировочных параметров Л1, Л2 и количества п базисных функций.

С этой целью:

а) Рассчитываются коэффициенты разложения сг путем минимизации оценочной функции Д2 для набора значений регуляризаторов А1 и Л2 в достаточно широком интервале при большом количестве базисных функций.

б) Выбирается оптимальное значение Л1. Очевидно, что при малых значениях А1 соответствующий член функции Д2 не оказывает влияния на решение уравнения (2), в то время как при больших значениях А1 минимизация функции Д2 приводит к занижению коэффициентов Сг и, следовательно, к искажению решения. При оптимальном выборе А1 подавляются указанные выше нежелательные осцилляции, и при этом не искажается форма искомой интенсивности рассеяния. Для выбора этого оптимального значения А1 в настоящей работе использован метод «точки перегиба» [11]. Суть метода состоит в том, что в форму-

М / п \ 2

ле (7) в окрестности оптимального значения А1 величины ^ I 7(9к, фк) — ^ ^Р^к, фк) )

&=1 V г=1 )

2

2

п п

и Е как функции А1 остаются приблизительно постоянными, причем величина ^ с2 име-

г=1 г=1

ет перегиб. Эти особенности поведения указанных величин в окрестности оптимального значения А1 проявляются только при выборе достаточно большого количества п базисных функций.

в) Если оптимальная величина А1 определена, то оптимальную величину А2 можно подобрать аналогичным способом. Как и в предыдущем случае, в окрестности оптималь-

М / п \ 2

ного значения А2 величины ^ 3(вк, ук) - ^ СгР{(9к, ) ) и ^

А:=1 V г=1 / кеК

Е С] $¡3 (°к

3 = 1

2

как

функции А2 остаются приблизительно постоянными. Отличие состоит в том, что величина

обнаруживает перегиб только в том случае, если при минимальных

Е

кек

Е сз(вк ) з=1

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

4. Модифицированный итерационный метод Фридмана

Второй подход к решению основного интегрального уравнения (2) использует итерационный метод Фридмана [8, 12], который, однако, необходимо существенно модифицировать.

Согласно классическому итерационному методу Фридмана выбирается начальное приближение 10 (например, 10 (9,ф) = 0), и каждое последующее приближение в каждой точке (д, ф) получается на основе рекуррентного соотношения

1п+ф9, ф = 1п(9, ф + и [3(в, ф — Зп(в, ф] , (9)

где Зп = А(1п), А(1) — интегральный оператор (4), а и — числовой параметр, регулирующий скорость сходимости решения. Очевидно, что если разность в квадратных скобках в уравнении (9) для любой точки (в, ф) равна нулю, тогда 1п является решением уравнения (2) . Последовательность приближений 1п с увеличение п сначала сходится, но затем начинает расходиться, не достигнув решения, достаточно близкого к точному решению [8].

Согласно (9), для приближения к точному решению нужно изменить значение 1п в каждой точке (в, ф) так, чтобы разность [3(в, ф) — Зп(в, ф] уменьшилась. Однако проблема заключается в том, что изменение значения 1п в точке (в, ф) может привести к изменению значения Зп не в точке (в, ф), а во множестве других точек. В этом случае изменение 1п в точке (в, ф) на основе разности [3(в, ф) — Зп(в, ф] теряет смысл.

Для исправления этого недостатка необходимо учесть влияние значения 1п(9, ф) на значения Зп в различных точках. С этой целью для оператора А вычисляется «функция влияния» ША : [0тт, ^тах] X [^тт, ^тах] X [9тт, ^тах] X [^тт, ^тах] ^ К, которая учитывает влияние аргумента оператора А (функции 1п) в точке (£, ц) на значение функции Зп в другой точке (в, ф).

Если «функция влияния» вычислена, то итерационную процедуру (9) можно переписать в виде:

^шах ^тах

У У и(фч,9,ф ■ (3(в,ф) — Зп(9,ф) АвАф . (10)

1п+1(^,'Ц) = 1п(Фч) + V

С помощью описанной модификации учитывается влияние значения функции 1п в точке

(фц) на все ,1п(д,ф).

Чтобы вычислить «функцию влияния» для оператора А, необходимо рассмотреть действие этого оператора на функции, отличные от нуля только в малой окрестности точки (£, ф). В качестве таких функций были выбраны «колокола» небольшой ширины с центром в точке (£, ц):

/ае, л') = (11)

где т — параметр, регулирующий ширину «колокола». Тогда в качестве «функции влияния» можно взять функцию

„({,„,*,„) =_ШШ_, (12)

/ / Шсф(е,ф ¿вйф

'Ушт ^шт

<?шах ^шах

для которой выполняется условие нормировки / / ш АвАф = 1.

<Ушт ^шт

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

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

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

Таким образом, модифицированный метод Фридмана предполагает: а) нахождение «функции влияния» (12) с использованием оператора (4) и функций (11) и б) применение итерационной процедуры (10) для поиска точного решения.

5. Оценка эффективности методов для модельных анизотропных систем

Для проверки эффективности методов был выбран ряд модельных анизотропных систем, для каждой из которых:

а) рассчитывалась исходная точечная интенсивность рассеяния I(в, ф);

б) с использованием параметров коллимационной системы, задаваемых при получении экспериментальных кривых (к1, Ь и т.д. —см. рис. 1.), по формуле (3) рассчитывалась весовая функция Ш (х), а затем по формуле (2) рассчитывалась искаженная интенсивность рассеяния ,1 (в, ф);

в) на основе развитых методов (метода базисных функций и модифицированного метода Фридмана) «восстанавливалась» точечная интенсивность I'(в, ф);

г) с использованием I'(в, ф) как новой точечной интенсивности, вновь рассчитывалась по формуле (2) искаженная интенсивность рассеяния ,1'(в, ф);

д) вычислялось относительное отклонение для точечных интенсивностей

\\i-n

1И1

\

м

Е(/ (ок ) - V (дк ))2 ^^--(13)

Е (! (Ок,4>к)? к=1

е) вычислялось относительное отклонение для искаженных интенсивностей

р - ? в

\

м 2

Е и (0к, Рк) - у (вк, Ук))2 к=1_

М

Е и (0к,^к))2 к=1

(14)

Для всех модельных структур интенсивность рассчитывалась в интервале углов 0.5 мрад ^ в ^ 10.5 мрад и 0 ^ ф ^ ж/2 при длине волны \сиКа = 0,1542 нм.

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

Вытянутые параллелепипеды. Пусть длинное ребро параллелепипеда ориентировано по оси ъ и два других ребра равны между собой. Ось симметрии образца выбрана параллельно длинному ребру. Интенсивность рассеяния от такого параллелепипеда [13] была усреднена по углу поворота относительно оси Z обратного пространства. С учетом указанного усреднения полученная окончательная формула для интенсивности рассеяния (в произвольных единицах) имеет вид

2

1(х,у,г) = -я (пег)

ж

тг/2

5

Ж

\.Т

(X2 + У2) 0С8 (у)

Б

IV

(X2 + У2) 8т(у) ¿у, (15)

где 5 (ж)

/ втаЛ V ж )

с — длина длинного ребра параллелепипеда, а — длина двух других

равных ребер.

Ламеллы. Ось симметрии образца выбрана перпендикулярно ламеллам. Как известно [13], системе бесконечных равноудаленных плоскостей в реальном пространстве, ориентированных перпендикулярно оси Z, соответствует в обратном пространстве интенсивность рассеяния, сконцентрированная в равноудаленных точках, расположенных по оси Z. Ограниченность плоскостей в реальном пространстве приводит к тому, что в каждой такой точке происходит уширение интенсивности рассеяния в плоскости, перпендикулярной оси Z. Считая это уширение радиально симметричным относительно оси Z, его можно учесть с помощью функции Гаусса. Учет толщины ламелл приводит к модуляции интенсив-

ности по оси Z в соответствии с формулой

Бт X

X

Таким образом, общая интенсивность

рассеяния (в произвольных единицах) системой ламелл может быть представлена как

I (X, У, г) = ехр

1 /X2 УЛ1 (

2 V

вт (^гbZ) \'

жьг

N

ехр

п=-N

24

(16)

где й — расстояние между ламеллами, Ь — толщина ламеллы. Выбранные значения параметров ах, оу, oz, N,а также и Ь будут указаны ниже. Величина выбиралась таким образом, чтобы интенсивность рассеяния имела помимо центрального максимума либо один, либо три пика в указанном выше диапазоне углов рассеяния.

При переходе в координаты в и ф в формулах (15) и (16) следует выбрать X = 0,

У

в Бт ф

и Z

в осе ф

Т

2

Ниже для оценки эффективности изложенных методов решения коллимационной задачи рассмотрены интенсивности рассеяния этими моделями.

5.1. Метод базисных функций

Согласно пунктам а)-г), для каждой модели были последовательно рассчитаны точечная и искаженная интенсивности рассеяния; методом «точки перегиба» из ряда значений Л1 и Л2 были выбраны оптимальные; определены восстановленная точечная и восстановленная искаженная интенсивности рассеяния.

Вытянутые параллелепипеды. В качестве параметров параллелепипеда были выбраны: а = 20 нм и с = 40 нм. На рисунках 2(а) и 3(а) приводятся, соответственно, точечная и искаженная интенсивности рассеяния. Рассчитанные оптимальные значения А1 = 10_8 и Л2= 3.162-10"4, количество сплайновых функций п = щ х п^, где щ = 20 (количество базисных функций по координате в) и п^ = 10 (количество базисных функций по координате ф). На рисунке 2(б) для данных величин А1 и А2 представлена восстановленная точечная интенсивность рассеяния.

а) б) в)

РИС. 2. Интенсивности рассеяния вытянутыми параллелепипедами: а) точечная; б) восстановленная точечная (метод базисных функций); в) восстановленная точечная (итерационный метод)

РИС. 3. Искаженные интенсивности рассеяния: а) вытянутыми параллелепипедами; б) ламеллами (в, = 25 нм); в) ламеллами (в, = 50 нм).

Ламеллы. Были выбраны: толщина ламеллы Ь =10 нм, параметры ах = 0, 02 пм,

ау = 0,1 пм, oz = 0, 05 пм, N = 10.

Расстояние между ламеллами й = 25 нм. При данном значении й в указанном выше диапазоне углов рассеяния точечная интенсивность рассеяния имеет помимо центрального максимума один пик. На рисунках 4(а) и 3(б) приводятся, соответственно, точечная и искаженная интенсивности рассеяния. Рассчитанные оптимальные значения А1 = 10_6 и Х2 = 10"5, количество сплайновых функций п = щ х п^, где щ = 20 и п^ = 10. На рисунке 4(б) для данных величин А1 и Л2 представлена восстановленная точечная интенсивность рассеяния.

Расстояние между ламеллами д, = 50 нм. При данном значении д, в указанном выше диапазоне углов рассеяния точечная интенсивность рассеяния имеет помимо центрального максимума три пика. На рисунках 5(а) и 3(в) приводятся, соответственно, точечная и искаженная интенсивности рассеяния. Рассчитанные оптимальные значения А1 =10"8 и Х2 = 1.778-10"5, количество сплайновых функций п =щ х п^, где щ = 20 и Пу = 10. На рисунке 5(б) для данных величин А1 и Х2 представлена, соответственно, восстановленная точечная интенсивность рассеяния.

а) б) в)

РИС. 4. Интенсивности рассеяния ламеллами (й = 25 нм): а) точечная; б) восстановленная точечная (метод базисных функций); в) восстановленная точечная (итерационный метод)

а) б) в)

РИС. 5. Интенсивности рассеяния ламеллами (й = 50 нм): а) точечная; б) восстановленная точечная (метод базисных функций); в) восстановленная точечная (итерационный метод)

5.2. Модифицированный метод Фридмана

Согласно пунктам а)—г), для каждой модели были последовательно рассчитаны точечная и искаженная интенсивности рассеяния; по формуле (12) с использованием оператора (4) и функций (11) была вычислена функция влияния; с помощью итерационной

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

Вытянутые параллелепипеды. Использовались те же геометрические параметры параллелепипедов: а = 20 нм и с = 40 нм. На рисунке 2(в) представлена восстановленная точечная интенсивность рассеяния, рассчитанная для оптимального количества итераций иопт= 20000.

Ламеллы. Были выбраны та же толщина ламеллы Ь = 10 нм и те же параметры ах = 0, 02 пм, ау = 0,1 пм, oz = 0,05 пм, N = 10.

Расстояние между ламеллами й = 25 нм. На рисунке 4(в) представлена восстановленная точечная интенсивность рассеяния, рассчитанная для оптимального количества итераций иопт = 80000.

Расстояние между ламеллами = 50 нм. На рисунке 5(в) представлена, восстановленная точечная интенсивность рассеяния, рассчитанная для оптимального количества итераций пош = 80000.

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

6. Применение развитых методов для коллимационного пересчета

экспериментальной интенсивности рассеяния анизотропным образцом

Были использованы экспериментальные индикатрисы для отожженных пленок из высокоориентированного полиэтилена [14], полученные в интервале углов рассеяния 9 от 0,8 до 23,5 миллирадиан и при семи различных углах ф поворота образца относительно коллимационной щели от 0° до 90° с шагом в 15°.

Следующий этап работы состоит в получении на основе указанных индикатрис рассеяния «экспериментальной» поверхности ,1 (в, ф). Трудность состоит в том, что имеющиеся исходные кривые рассеяния измерены с большим шагом по углу ф. Для получения достаточно объективной поверхности ,1 (9,ф) можно для промежуточных значений угла Ф рассчитать дополнительные кривые рассеяния. С этой целью каждая экспериментальная индикатриса нормировалась на ее значение в максимуме интенсивности, и в плоскости (9, ф) строились линии одинаковой интенсивности. Вдоль каждой такой линии проводилась интерполяция исходной (ненормированной) интенсивности к тем значениям углов 9 и ф, которые отвечают дополнительным кривым рассеяния. «Экспериментальная» поверхность 3 (9,ф), полученная с помощью описанной процедуры, представлена на рис. 6. Эта поверхность является достаточно гладкой, причем сечения, отвечающие экспериментальным значениям ф, полностью совпадают с соответствующими измеренными индикатрисами рассеяния.

Описанные методы были применены для получения восстановленной точечной интенсивности рассеяния I'(9,ф), и затем в соответствии с формулой (2) были получены восстановленные «экспериментальные» интенсивности ,1'(9,ф).

Метод базисных функций. Рассчитанные оптимальные значения А1 = 3,2-10-7 и \2 = 10-4, количество сплайновых функций п =щ х п^, где щ = 20 и п^ = 10. На рисунке 7(а) для данных величин А1 и А2 представлена восстановленная точечная интенсивность рассеяния.

Модифицированный метод Фридмана. По формуле (12) с использованием оператора (4) и функций (11) была вычислена функция влияния; с помощью итерационной процедуры (10) рассчитана восстановленная точечная интенсивность. Параметр, регулирующий скорость сходимости решения ^опт = 10. На рисунке 7(б) представлена восстановленная точечная интенсивность рассеяния, рассчитанная для оптимального количества итераций иопт = 3100.

Рис. 6. «Экспериментальная» интенсивность для отожженных пленок из высокоориентированного полиэтилена (см. текст)

РИС. 7. Восстановленные точечные интенсивности рассеяния для отожженных пленок из высокоориентированного полиэтилена: а) метод базисных функций; б) модифицированный итерационный метод Фридмана

7. Обсуждение результатов

Модельные системы. Как видно из рисунков 2, 4, 5, в рамках каждой модели наблюдается хорошее согласие между исходной точечной и восстановленной точечной ин-тенсивностями рассеянии, рассчитанными обоими методами. В то же время имеет место сильное различие между точечными и искаженными интенсивностями в рамках каждой модели (ср.: рис. 2 и рис. 3(а); рис. 4 и рис. 3(б); рис. 5 и рис. 3(в)).

Для количественного сравнения поверхностей соответствующих интенсивностей

было рассчитаны относительные отклонения -——-—- (13) и -—^—- (14). Результаты,

приводимые в таблице 1, свидетельствуют о том, что для рассмотренных модельных систем относительные отклонения -—цуц—- не превышают 3%, а относительные отклонения

^ - составляют 0,5% и менее. И хотя эти отклонения, как видно, сами по себе малы —

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

ТАБЛИЦА 1. Анизотропные модельные системы. Относительные отклонения для исходной точечной и восстановленной точечной интенсивностей рассеяния и относительные отклонения для искаженной и восстановленной искаженной интенсивностей рассеяния

Модель Метод базисных функций Модифицированный метод Фридмана

с 3 3' с 3 3' оу

I ' 3 ' /и I ' 3 '

Ориентированные параллелепипеды 2,6 0,55 0,8 0,08

Ламеллы, 3 =25 нм 1,8 0,24 1,27 0,07

Ламеллы, с1 = 50 нм 2,0 0,23 0,6 0,02

ТАБЛИЦА 2. Анизотропный образец. Относительные отклонения для «экспериментальной» искаженной и восстановленных искаженных интенсивностей рассеяния и относительное отклонение восстановленных точечных интенсив-ностей, рассчитанных двумя методами

Образец

Метод базисных функций

3

%

Модифицированный метод Фридмана

3

%

Сравнение методов

V _ т' 12

%

Отожженная пленка из высокоориентированного полиэтилена

3,0

4,8

19,9

РИС. 8. Сравнение экспериментальных индикатрис рассеяния (точки) при разных углах ориентации образца (ф =0, ^/4, ж/2) с сечениями «экспериментальной» интенсивности, восстановленной методом базисных функций (сплошная линия)

РИС. 9. Сравнение экспериментальных индикатрис рассеяния (точки) при разных углах ориентации образца (ф =0, ж/4, ж/2) с сечениями «экспериментальной» интенсивности, восстановленной модифицированным итерационным методом Фридмана (сплошная линия)

«Экспериментальная» интенсивность рассеяния на отожженной пленке из высокоориентированного полиэтилена. В случае «экспериментальной» интенсивности рассеяния, разумеется, отсутствует исходная точечная интенсивность, но возможно сравнение точечных интенсивностей, восстановленных с помощью двух развитых методов (рис. 7). С этой

\А\

новленная методом базисных функций, 1'2 — интенсивность, восстановленная итерационным методом. Кроме того, проводилось сравнение «экспериментальной» интенсивности (то есть поверхности, построенной на основе экспериментальных индикатрис рассеяния) с восстановленными искаженными интенсивностями, рассчитанными обоими методами. Для количественного сравнения поверхностей этих интенсивностей вновь рассчитывалось

\и - ^ I

относительное отклонение —-—-— для обоих методов. Результаты, представленные в таб-

целью вычислялось относительное отклонение —\ Т/ м 1 , где /( — интенсивность, восста-

лице 2, свидетельствуют о том, что относительное отклонение ^ ^ ^ составляет 3-5%,

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

Сравнение рис. 6 и рис. 7 обнаруживает разительное различие между «экспериментальной» интенсивностью и точечными интенсивностями, восстановленными с помощью обоих методов. Что касается сравнения восстановленных точечных интенсивностей

III' - I' II

рассеяния, то их относительное отклонение —* 2 оказывается ~20%. Это обстоятель-

НЛ11

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

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

Литература

[1] Guinier A., Fournet G. Small-angle Scattering of X-rays, New-York: Wiley, 1955. — 268 р.

[2] Porod G. General Theory in Small-angle X-ray scattering, London: Academic Press, 1983, P. 17-52.

[3] Kratky O. Instrumentation, Experimental Technique, Slit Collimation, in Small-angle X-ray scattering, London: Academic Press, 1983, P. 53-84.

[4] Федоров Б.А. Учет коллимационных искажений при малоугловом рассеянии рентгеновых лучей. Поправка на высоту щелей // Кристаллография, 1968, Т. 13, № 5, C. 763-769.

[5] Schelten J., Hossfeld F. Application of spline functions to the correction of resolution errors in small angle scattering // J. Appl. Cryst. 1971, V. 4(3), P. 210-223.

[6] Сизиков В.С., Смирнов А.В., Федоров Б.А. Решение одномерной коллимационной задачи оценки рентгеновского изотропного рассеяния излучения методом итераций // Изв. вузов. Приборостроение, 2005, Т. 48, № 10, C. 44-52.

[7] Смирнов А. В., Сизиков В. С., Федоров Б. А. Решение обратной коллимационной задачи для рентгеновского малоуглового изотропного рассеяния с помощью сплайновых функций // Изв. вузов. Приборостроение, 2006, Т. 49, № 1, C. 41-47.

[8] Захаров Д. Д., Сизиков В.С., Смирнов А. В., Федоров Б. А. Решение двумерной коллимационной задачи рассеяния рентгеновских лучей с использованием нестандартных интегральных уравнений // Научно-технический вестник СПбГУИТМО, 2006, Т. 32, C. 144-153.

[9] Тихонов А. Н., Арсенин В.Я. Методы решения некорректных задач. Наука, М., 1979, 284 с.

[10] Greville T. N. E. Introduction to spline functions // In: Theory and Applications of Spline Functions. Academic Press, New York, 1969, 1-35.

[11] Glatter O. Data Evaluation in Small Angle Scattering: Calculation Electron Density Distribution by Means of Indirect Fourier Transformation// Acta Physica Austriaca. 1977, V. 47, P. 83-102.

[12] Верлань А.Ф., Сизиков В.С. Интегральные уравнения: методы, алгоритмы, программы. Наук. думка, Киев, 1986. — 548 с.

[13] Вайнштейн Б. К. Дифракция рентгеновских лучей на цепных молекулах. Издательство АН СССР, М., 1963. — 372 с.

[14] Шаташвили М. В., Ресовский А. В., Смирнов А. В., Федоров Б. А., Курындин И. С., Ельяшевич Г.К. Исследование пористых полиэтиленовых пленок методом малоуглового рентгеновского рассеяния // Высокомолекулярные соединения, сер. А, 2005, Т. 47, № 6, C. 970-977.