АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ ПРИ МАЛОУГЛОВОМ РЕНТГЕНОВСКОМ РАССЕЯНИИ НА ОРИЕНТИРОВАННЫХ ОДНОРОДНЫХ ЦИЛИНДРАХ. ОЦЕНКА ВЛИЯНИЯ ЭКСПЕРИМЕНТАЛЬНОЙ ПОГРЕШНОСТИ НА КАЧЕСТВО КОЛЛИМАЦИОННОГО ПЕРЕСЧЕТА
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
6 МАТЕРИАЛОВЕДЕНИЕ И НАНОТЕХНОЛОГИИ
УДК 517.968, 539.261
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ ПРИ МАЛОУГЛОВОМ РЕНТГЕНОВСКОМ РАССЕЯНИИ НА ОРИЕНТИРОВАННЫХ ОДНОРОДНЫХ ЦИЛИНДРАХ. ОЦЕНКА ВЛИЯНИЯ ЭКСПЕРИМЕНТАЛЬНОЙ ПОГРЕШНОСТИ
НА КАЧЕСТВО КОЛЛИМАЦИОННОГО ПЕРЕСЧЕТА
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
Для ориентированных однородных цилиндров получены интенсивности малоуглового рентгеновского рассеяния, в них внесены коллимационные искажения. На основе развитых ранее методов коллимационного пересчета рассчитаны восстановленные интенсивности рассеяния, которые сравнивались с неискаженными интенсивностями. Показано, что возникающая при этом относительная погрешность составляет несколько процентов, что свидетельствует о высоком качестве коллимационного пересчета. Показано также, что внесение в искаженные интенсивности рассеяния статистического шума в 1% и 4% – что соответствует реальным условиям рентгеновского эксперимента – лишь незначительно ухудшает совпадение исходных и восстановленных интенсивностей рассеяния. Ключевые слова: малоугловое рентгеновское рассеяние, ориентированные однородные цилиндры, метод базисных функций, итерационный метод Фридмана.
Введение
Использование коллимационных щелей при рентгеновском малоугловом рассеянии приводит к искажению экспериментальных индикатрис, что существенно затрудняет их интерпретацию с точки зрения структуры рассеивающего объекта [1]. Для устранения этого недостатка необходимо ввести коллимационные поправки, позволяющие перейти от «искаженной» интенсивности рассеяния J при коллимационной системе «щель–щель» к «точечной» интенсивности рассеяния I при коллимационной системе «точка–точка». В случае изотропных объектов, для которых интенсивность зависит только от угла рассеяния , такие поправки вносятся сравнительно легко [2–5]. Для анизотропных объектов интенсивности J (,) и I (,) являются двухпараметрическими функциями, зависящими как от угла рассеяния , так
и от угла поворота рассеивающего образца относительно коллимационной щели. При этом задача
внесения коллимационных поправок становится несравненно более сложной [6]. В недавней работе [7] для решения этой задачи были развиты два метода. Первый из них – метод ба-
зисных функций (МБФ) [8] – состоит в том, что экспериментальная интенсивность рассеяния представляется в виде линейной комбинацией двухмерных базисных функций и с помощью модифицированного метода наименьших квадратов достигается наилучшее согласие между исходной экспериментальной интенсивностью рассеяния и интенсивностью, полученной с помощью разложения по базисным функциям. Второй метод – модифицированный метод Фридмана (ММФ), основанный на итерационном методе Фридмана [9], существенно дополненном в соответствии с особенностями настоящей задачи. Показано, что оба метода дают достаточно высокую точность коллимационного пересчета как в случае рассеяния на ряде анизотропных модельных объектов, так и в случае экспериментальных индикатрис рассеяния при анизотропном рассеянии на отожженных пленках из высокоориентированного полиэтилена.
В настоящей работе с помощью указанных методов рассмотрен еще один модельный объект – ориентированные цилиндры разной вытянутости, и на его основе проведена оценка влияния экспериментальных погрешностей на относительную ошибку, возникающую при внесении коллимационных поправок.
Основы метода базисных функций и модифицированного метода Фридмана [7]
J (,) I 2 2 , | arctg |W d ,
(1)
где W – «весовая функция», зависящая только от параметров рентгеновской установки; I (,) – то-
чечная интенсивность рассеяния, рассчитанная для выбранной анизотропной модели в условиях точеч-
ного источника и точечного приемника рентгеновского излучения; J (,) – искаженная интенсивность
рассеяния, рассчитанная по формуле (1) с учетом параметров используемой коллимационной системы; J~(,) – зашумленная («псевдоэкспериментальная») интенсивность рассеяния; I~(,) – восстановлен-
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 4 (74)
77
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ…
ная точечная интенсивность рассеяния, рассчитанная на основе развитых методов (МБФ и ММФ); J~(,) – восстановленная искаженная интенсивность рассеяния, рассчитанная по формуле (1) с использованием I~(,) как новой точечной интенсивности.
МБФ. Введем f – норму в пространстве квадратично интегрируемых функций двух перемен-
max max
ных и : f
f 2(,) d d . С использованием данной нормы относительное отклонение
min min
функции
f1 от функции
f2 можно определить как f1, f2
f1 f2 f2
.
n
Выбираются базисные функции fi (i 1,2,, n) так, чтобы их линейная комбинация I ci fi
i 1
могла хорошо описывать I . Для этих функций должны существовать такие коэффициенты c1,c2,,cn ,
для которых относительное отклонение I, I является достаточно малой величиной.
n
Рассмотрим уравнение J A(I) . Исходя из линейности оператора A , J ciFi , где
i 1
Fi A fi (i 1,2,...,n) . В поставленной задаче I является искомой функцией, и представить ее разложе-
ние по базисным функциям в явном виде нельзя. Однако можно найти для экспериментальной интенсив-
ности J ее приближение как разложение J по функциям Fi . При корректном выборе функций fi
можно определить коэффициенты, для которых J , J будет достаточно малым. Эти коэффициенты
будут использованы для построения восстановленной точечной интенсивности
n
I ci fi .
(2)
i 1
Так как исходная задача некорректна, то, очевидно, одной экспериментальной функции J соот-
ветствует множество решений I , для которых выполняется уравнение J A(I ) .Для того чтобы сузить
класс решений, следует: а) выбрать базисные функции, наилучшим образом соответствующие данной
задаче; б) ввести дополнительные ограничения на решение, используя его известные свойства. Двухмер-
ные базисные функции fi были построены как прямое произведение одномерных кубических
B-сплайнов [10]. Для получения коэффициентов разложения по базисным функциям использовался метод
наименьших квадратов (МНК), который был модифицирован исходя из двух положений. Первое состоит
в том, что при применении МНК к быстро спадающим функциям возникают коэффициенты ci , имеющие большие значения и противоположные знаки. При этом в решении появляются нежелательные осцилля-
ции, которых не должно быть в искомой функции. Второе положение связано с тем, что для нулевого
угла рассеяния ( 0 ) точечная интенсивность должна оставаться постоянной при изменении угла ,
т.е. I (0,) const .
Математический формализм, развитый с учетом указанных положений, приводит к системе линейных уравнений относительно коэффициентов ci :
n c j 1
j
M k 1
Fi
(k
,
k
)Fj
(k
,
k
)
2
fi
kK
(k
,
k
)f
j
(k
,
k
)
1ci
,
i
1,2,...,n ,
M
J (k ,k )Fi (k ,k )
k 1
(3)
где f j (k ,k )
f j (k ,k )
K
f j (k ,k )
1 NK
f j (k ,k )
kK
f j (k ,k ) .
Здесь k , k – значения угла рассеяния и угла поворота образца, для которых получена экспери-
ментальная интенсивность; M – общее число экспериментальных точек; K – множество точек, соответ-
ствующих min ; NK – количество точек в множестве K ; а 1 и 2 – параметры, выбор которых
детально рассмотрен в [4].
78 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 4 (74)
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
ММФ. Согласно классическому итерационному методу Фридмана выбирается начальное прибли-
жение I0 (например, I0, 0 ), и каждое последующее приближение в каждой точке (q,) получает-
ся на основе рекуррентного соотношения
In1(,) In(,) J (,) Jn(,) ,
(4)
где Jn A(In) , A(I ) – интегральный оператор, а – числовой параметр, регулирующий скорость сходимости решения.
Согласно (4), для приближения к точному решению нужно изменить значение In в каждой точке
(,) так, чтобы разность J (,) Jn(,) уменьшилась. Однако проблема заключается в том, что
изменение значения In в точке (,) может привести к изменению значения Jn не в точке (,) , а во
множестве других точек. В этом случае изменение In в точке (,) на основе разности
J (,) Jn(,) теряет смысл. Для исправления этого недостатка необходимо учесть влияние значения
In(,) на значения Jn в различных точках. Вводится и обосновывается «функция влияния»
(,,,) , которая учитывает влияние аргумента оператора A (функции In ) в точке (,) на значение
функции Jn в другой точке (,) . Итерационная процедура (4) с учетом «функции влияния» может быть представлена в виде
max max
In1(,)
In
(, )
(,,,) J (,) Jn (,) dd .
min min
(5)
Выбор оптимального количества итераций nопт , а также регулировочного параметра описан в
работе [7].
Оценка эффективности методов коллимационного пересчета на примере однородных ориентированных цилиндров
Из общей теории дифракции [11] следует, что если в реальном пространстве (x,y,z) цилиндр ориентирован по оси z, то в обратном пространстве (X,Y,Z), его интенсивность рассеяния имеет вид
I
X
,Y
,
Z
d
sinbZ J1 d
2Z X 2
X2 Y2
Y
2
2
,
(6)
где d – диаметр цилиндра, b – его высота. При переходе в координаты и в формуле (6) следует вы-
брать X 0 , Y sin и Z cos .
В работе рассмотрены два цилиндра с одинаковым радиусом, но разной высотой. Интенсивность рассчитывалась в интервале углов 0,5 мрад 10,5 мрад и 0 / 2 при длине волны
CuKα = 0,1542 нм. Для оценки влияния экспериментальной погрешности на качество коллимационного пересчета в
искаженную (экспериментальную) интенсивность J рассеяния в каждой точке k , k вносился статистический «шум», распределенный по нормальному закону со средним квадратичным отклонением
k ,k J k ,k 104 J 0,0,
причем в качестве общей относительной погрешности интенсивности (в %) выбиралось J , J~ , где J~ –
зашумленная интенсивность. Уровень шума можно менять, умножая интенсивность J во всех точках на один и тот же коэффициент. Для зашумленных («псевдоэкспериментальных») интенсивностей, получен-
ных с разным уровнем шума J , J~ , рассчитывались обоими методами восстановленные точечные
I~(,) интенсивности, и на их основе с использованием формул (2) или (5) рассчитывались соответствующие восстановленные искаженные интенсивности рассеяния J~(,) .
Параметры цилиндра: d = 20 нм, b = 40 нм. Исходя из выражения (6), была рассчитана точечная интенсивность рассеяния цилиндра I (,) (рис. 1, а), а по формуле (1) – искаженная интенсивность
J (,) , в которую был внесен «шум» в 1% и 4%, и получены зашумленные («псевдоэкспериментальные») интенсивности J~(,) (рис. 2, а). На приводимых ниже рисунках представлены интенсивности,
относящиеся только к шуму в 4%.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 4 (74)
79
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ… МБФ. По формулам (2) и (3) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 1, b). Далее на основе полученных функций I~(,) по формуле (1) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) как для шума в 1%, так и для шума в 4% (рис. 2, b). ММФ. По формуле (5) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) (рис. 2, а) по формуле (5) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 1, с). Далее на основе полученных функций I~(,) по формуле (5) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) для шума в 1% и 4% (рис. 2, с).
abc Рис. 1. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 40 нм: (а) – точечная интенсивность I (,) ; (b) – восстановленная точечная интенсивность I~(,) ,
МБФ, «шум» в 4%; (с) – восстановленная точечная интенсивность I~(,) , ММФ, «шум» в 4%
abc Рис. 2. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 40 нм: (а) – зашумленная («псевдоэкспериментальная») интенсивность J~(,) ; (b) – восстановленная искаженная интенсивность J~(,) , МБФ, «шум» в 4%; (с) – восстановленная искаженная интенсивность J~(,) ,ММФ, «шум» в 4%
Параметры цилиндра: d = 20 нм, b = 80 нм. Исходя из выражения (6), была рассчитана точечная интенсивность рассеяния цилиндра I (,) (рис. 3, а), а по формуле (1) – искаженная интенсивность J (,) , в которую был внесен «шум» в 1% и 4%, и получены зашумленные («псевдоэкспериментальные») интенсивности J~(,) (рис. 4, а). На приводимых ниже рисунках представлены интенсивности, относящиеся только к шуму в 4%.
МБФ. По формулам (2) и (3) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 3, b). Далее на основе полученных функций I~(,) по формуле (1) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) как для шума в 1%, так и для шума в 4% (рис. 4, b).
80 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 4 (74)
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
ММФ. По формуле (5) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) по формуле (5) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 3, с). Далее на основе полученных функций I~(,) по формуле (5) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) для шума в 1% и 4% (рис. 4, с).
abc Рис. 3. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 80 нм: (а) – точечная интенсивность I (,) ; (b) – восстановленная точечная интенсивность I~(,) ,
МБФ, «шум» в 4%; (с) – восстановленная точечная интенсивность I~(,) , ММФ, «шум» в 4%
abc Рис. 4. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 80 нм: (а) – зашумленная («псевдоэкспериментальная») интенсивность J~(,) ; (b) – восстановленная искаженная интенсивность J~(,) , МБФ, «шум» в 4%; (с) – восстановленная искаженная интенсивность J~(,) , ММФ, «шум» в 4%
Заключение
Результаты, отражающие качество коллимационного пересчета, представлены в таблице.
Уровень шума 1
J, J~, %
0 1,0 4,0
0 1,0 4,0
Метод базисных функций
Модифицированный метод Фридмана
Расхождение результатов
23
45 67
8
I, I~, % J~, J~ ,% J, J~,% I, I~,%
J~, J~ ,%
J, J~ ,%
I~М БФ, I~М МФ ,%
Однородный ориентированный цилиндр с параметрами d = 20 нм, b = 40 нм
1,0 0,6 0,6
3,3 0,8 0,08
2,6
1,4 1,1 0,6
3,1 0,9 0,3
2,9
2,5 4,0 1,0
3,4 3,8 0,9
2,8
Однородный ориентированный цилиндр с параметрами d = 20 нм, b = 80 нм
2,7 0,8 0,8
3,8 0,3 0,3
2,6
3,4 1,4 1,1
3,9 0,9 0,4
5,3
5,6 4,1 1,5
4,8 3,8 1,1
5,2
Таблица. Относительное расхождение между интенсивностями: точечной и точечной восстановленной; искаженной и искаженной восстановленной для рассеяния на однородном ориентированном цилиндре
при различном уровне вносимого «шума»
Из анализа этой таблицы можно сделать следующие заключения.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 4 (74)
81
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ…
1. Погрешность восстановления как точечных, так и искаженных интенсивностей рассеяния нигде не превышает нескольких процентов даже с учетом «псевдоэкспериментальной» зашумленности.
2. Оба метода дают приблизительно одинаковый результат внесения коллимационных поправок, хотя, как можно видеть из сравнения столбцов 2 и 5, МБФ «работает» несколько лучше, если «шум» не слишком большой и цилиндры не слишком вытянутые. Из сравнения точечных интенсивностей рассеяния, восстановленных методами МБФ и ММФ (столбец 8), следует, что результаты применения обоих методов согласуются с точностью в несколько процентов.
3. Зашумление искаженной интенсивности, т.е. переход от J к J~ , естественно, приводит к появлению расхождения между этими функциями (столбец 1). Интересно, однако, отметить, что если эту зашумленную искаженную интенсивность J~ «пропустить» через двойной коллимационный пересчет J~ I~ J~ , то полученная восстановленная искаженная интенсивность J~ оказывается заметно ближе к незашумленной искаженной интенсивности J , чем просто зашумленная искаженная интенсивность J~ : ср. столбцы 1 и 4 (МБФ) и 1 и 7 (ММФ). Эту закономерность можно также проследить при сравнении рис. 2, а–с, (цилиндр 2040) и рис. 4, а–с, (цилиндр 2080). Как видно, поверхности J~ на рис. 2, а, и рис. 4, а, оказываются существенно более шероховатыми по сравнению с поверхностями J~ на рис. 2, b, с, и рис. 4, b, с, соответственно.
Литература
1. Симаков А.П., Федоров Б.А., Смирнов А.В. Исследование поверхности тонкодисперсных фракций ряда минералов с помощью рентгеновского малоуглового рассеяния // Научно-технический вестник СПбГУ ИТМО. – 2009. – № 2(60). – С. 59–66.
2. Kratky O. Instrumentation, Experimental Technique, Slit Collimation, in Small-angle X-ray scattering. – London: Academic Press, 1983. – Р. 53–84.
3. Федоров Б.А. Учет коллимационных искажений при малоугловом рассеянии рентгеновых лучей. Поправка на высоту щелей // Кристаллография. – 1968. – Т. 13. – № 5. – С. 763–769.
4. Schelten J., Hossfeld F. Application of spline functions to the correction of resolution errors in small angle scattering // J. Appl. Cryst. – 1971. – V. 4. – Р. 210–223.
5. Фомичева Е.Е., Темнов Д.Э., Смирнов А.В., Федоров Б.А. Влияние дисперсионного наполнителя на основе алюминия на структуру и свойства полипропилена // Научно-технический вестник СПбГУ ИТМО. – 2009. – № 6(64). – С. 59–65.
6. Захаров Д.Д., Сизиков В.С., Смирнов А.В., Федоров Б.А. Решение двухмерной коллимационной задачи рассеяния рентгеновских лучей с использованием нестандартных интегральных уравнений // Научно-технический вестник СПбГУ ИТМО. – 2006. – № 32. – С. 144–153.
7. Захаров Д.Д., Смирнов А.В., Федоров Б.А. Решение коллимационной задачи при малоугловом рентгеновском рассеянии на анизотропных объектах // Наносистемы: физика, химия, математика. – 2011. – Т. 2. – № 3. – С. 26–44.
8. Glatter O. Data Evaluation in Small Angle Scattering: Calculation Electron Density Distribution by Means of Indirect Fourier Transformation // Acta Physica Austriaca. – 1977. – V. 47. – Р. 83–102.
9. Верлань А.Ф., Сизиков В.С. Интегральные уравнения: методы, алгоритмы, программы. – Киев: Наукова думка, 1986. – 543 с.
10. Greville T.N.E. Introduction to spline functions. In Theory and Applications of Spline Functions. – New
York: Academic Press, 1969. – Р. 1–35.
11. Вайнштейн Б.К. Дифракция рентгеновских лучей на цепных молекулах. – М.: Изд-во АН СССР, 1963. – 372 с.
Смирнов Александр Витальевич
Захаров Денис Дмитриевич Федоров Борис Александрович
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат физ.-мат. наук, доцент,
smirnav@phd.ifmo.ru – Санкт-Петербургский государственный университет информационных тех-
нологий, механики и оптики, аспирант, denzakharov@gmail.com – Санкт-Петербургский государственный университет информационных тех-
нологий, механики и оптики, доктор физ.-мат. наук, профессор,
borfedorov@rambler.ru
82 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 4 (74)
6 МАТЕРИАЛОВЕДЕНИЕ И НАНОТЕХНОЛОГИИ
УДК 517.968, 539.261
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ ПРИ МАЛОУГЛОВОМ РЕНТГЕНОВСКОМ РАССЕЯНИИ НА ОРИЕНТИРОВАННЫХ ОДНОРОДНЫХ ЦИЛИНДРАХ. ОЦЕНКА ВЛИЯНИЯ ЭКСПЕРИМЕНТАЛЬНОЙ ПОГРЕШНОСТИ
НА КАЧЕСТВО КОЛЛИМАЦИОННОГО ПЕРЕСЧЕТА
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
Для ориентированных однородных цилиндров получены интенсивности малоуглового рентгеновского рассеяния, в них внесены коллимационные искажения. На основе развитых ранее методов коллимационного пересчета рассчитаны восстановленные интенсивности рассеяния, которые сравнивались с неискаженными интенсивностями. Показано, что возникающая при этом относительная погрешность составляет несколько процентов, что свидетельствует о высоком качестве коллимационного пересчета. Показано также, что внесение в искаженные интенсивности рассеяния статистического шума в 1% и 4% – что соответствует реальным условиям рентгеновского эксперимента – лишь незначительно ухудшает совпадение исходных и восстановленных интенсивностей рассеяния. Ключевые слова: малоугловое рентгеновское рассеяние, ориентированные однородные цилиндры, метод базисных функций, итерационный метод Фридмана.
Введение
Использование коллимационных щелей при рентгеновском малоугловом рассеянии приводит к искажению экспериментальных индикатрис, что существенно затрудняет их интерпретацию с точки зрения структуры рассеивающего объекта [1]. Для устранения этого недостатка необходимо ввести коллимационные поправки, позволяющие перейти от «искаженной» интенсивности рассеяния J при коллимационной системе «щель–щель» к «точечной» интенсивности рассеяния I при коллимационной системе «точка–точка». В случае изотропных объектов, для которых интенсивность зависит только от угла рассеяния , такие поправки вносятся сравнительно легко [2–5]. Для анизотропных объектов интенсивности J (,) и I (,) являются двухпараметрическими функциями, зависящими как от угла рассеяния , так
и от угла поворота рассеивающего образца относительно коллимационной щели. При этом задача
внесения коллимационных поправок становится несравненно более сложной [6]. В недавней работе [7] для решения этой задачи были развиты два метода. Первый из них – метод ба-
зисных функций (МБФ) [8] – состоит в том, что экспериментальная интенсивность рассеяния представляется в виде линейной комбинацией двухмерных базисных функций и с помощью модифицированного метода наименьших квадратов достигается наилучшее согласие между исходной экспериментальной интенсивностью рассеяния и интенсивностью, полученной с помощью разложения по базисным функциям. Второй метод – модифицированный метод Фридмана (ММФ), основанный на итерационном методе Фридмана [9], существенно дополненном в соответствии с особенностями настоящей задачи. Показано, что оба метода дают достаточно высокую точность коллимационного пересчета как в случае рассеяния на ряде анизотропных модельных объектов, так и в случае экспериментальных индикатрис рассеяния при анизотропном рассеянии на отожженных пленках из высокоориентированного полиэтилена.
В настоящей работе с помощью указанных методов рассмотрен еще один модельный объект – ориентированные цилиндры разной вытянутости, и на его основе проведена оценка влияния экспериментальных погрешностей на относительную ошибку, возникающую при внесении коллимационных поправок.
Основы метода базисных функций и модифицированного метода Фридмана [7]
J (,) I 2 2 , | arctg |W d ,
(1)
где W – «весовая функция», зависящая только от параметров рентгеновской установки; I (,) – то-
чечная интенсивность рассеяния, рассчитанная для выбранной анизотропной модели в условиях точеч-
ного источника и точечного приемника рентгеновского излучения; J (,) – искаженная интенсивность
рассеяния, рассчитанная по формуле (1) с учетом параметров используемой коллимационной системы; J~(,) – зашумленная («псевдоэкспериментальная») интенсивность рассеяния; I~(,) – восстановлен-
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 4 (74)
77
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ…
ная точечная интенсивность рассеяния, рассчитанная на основе развитых методов (МБФ и ММФ); J~(,) – восстановленная искаженная интенсивность рассеяния, рассчитанная по формуле (1) с использованием I~(,) как новой точечной интенсивности.
МБФ. Введем f – норму в пространстве квадратично интегрируемых функций двух перемен-
max max
ных и : f
f 2(,) d d . С использованием данной нормы относительное отклонение
min min
функции
f1 от функции
f2 можно определить как f1, f2
f1 f2 f2
.
n
Выбираются базисные функции fi (i 1,2,, n) так, чтобы их линейная комбинация I ci fi
i 1
могла хорошо описывать I . Для этих функций должны существовать такие коэффициенты c1,c2,,cn ,
для которых относительное отклонение I, I является достаточно малой величиной.
n
Рассмотрим уравнение J A(I) . Исходя из линейности оператора A , J ciFi , где
i 1
Fi A fi (i 1,2,...,n) . В поставленной задаче I является искомой функцией, и представить ее разложе-
ние по базисным функциям в явном виде нельзя. Однако можно найти для экспериментальной интенсив-
ности J ее приближение как разложение J по функциям Fi . При корректном выборе функций fi
можно определить коэффициенты, для которых J , J будет достаточно малым. Эти коэффициенты
будут использованы для построения восстановленной точечной интенсивности
n
I ci fi .
(2)
i 1
Так как исходная задача некорректна, то, очевидно, одной экспериментальной функции J соот-
ветствует множество решений I , для которых выполняется уравнение J A(I ) .Для того чтобы сузить
класс решений, следует: а) выбрать базисные функции, наилучшим образом соответствующие данной
задаче; б) ввести дополнительные ограничения на решение, используя его известные свойства. Двухмер-
ные базисные функции fi были построены как прямое произведение одномерных кубических
B-сплайнов [10]. Для получения коэффициентов разложения по базисным функциям использовался метод
наименьших квадратов (МНК), который был модифицирован исходя из двух положений. Первое состоит
в том, что при применении МНК к быстро спадающим функциям возникают коэффициенты ci , имеющие большие значения и противоположные знаки. При этом в решении появляются нежелательные осцилля-
ции, которых не должно быть в искомой функции. Второе положение связано с тем, что для нулевого
угла рассеяния ( 0 ) точечная интенсивность должна оставаться постоянной при изменении угла ,
т.е. I (0,) const .
Математический формализм, развитый с учетом указанных положений, приводит к системе линейных уравнений относительно коэффициентов ci :
n c j 1
j
M k 1
Fi
(k
,
k
)Fj
(k
,
k
)
2
fi
kK
(k
,
k
)f
j
(k
,
k
)
1ci
,
i
1,2,...,n ,
M
J (k ,k )Fi (k ,k )
k 1
(3)
где f j (k ,k )
f j (k ,k )
K
f j (k ,k )
1 NK
f j (k ,k )
kK
f j (k ,k ) .
Здесь k , k – значения угла рассеяния и угла поворота образца, для которых получена экспери-
ментальная интенсивность; M – общее число экспериментальных точек; K – множество точек, соответ-
ствующих min ; NK – количество точек в множестве K ; а 1 и 2 – параметры, выбор которых
детально рассмотрен в [4].
78 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 4 (74)
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
ММФ. Согласно классическому итерационному методу Фридмана выбирается начальное прибли-
жение I0 (например, I0, 0 ), и каждое последующее приближение в каждой точке (q,) получает-
ся на основе рекуррентного соотношения
In1(,) In(,) J (,) Jn(,) ,
(4)
где Jn A(In) , A(I ) – интегральный оператор, а – числовой параметр, регулирующий скорость сходимости решения.
Согласно (4), для приближения к точному решению нужно изменить значение In в каждой точке
(,) так, чтобы разность J (,) Jn(,) уменьшилась. Однако проблема заключается в том, что
изменение значения In в точке (,) может привести к изменению значения Jn не в точке (,) , а во
множестве других точек. В этом случае изменение In в точке (,) на основе разности
J (,) Jn(,) теряет смысл. Для исправления этого недостатка необходимо учесть влияние значения
In(,) на значения Jn в различных точках. Вводится и обосновывается «функция влияния»
(,,,) , которая учитывает влияние аргумента оператора A (функции In ) в точке (,) на значение
функции Jn в другой точке (,) . Итерационная процедура (4) с учетом «функции влияния» может быть представлена в виде
max max
In1(,)
In
(, )
(,,,) J (,) Jn (,) dd .
min min
(5)
Выбор оптимального количества итераций nопт , а также регулировочного параметра описан в
работе [7].
Оценка эффективности методов коллимационного пересчета на примере однородных ориентированных цилиндров
Из общей теории дифракции [11] следует, что если в реальном пространстве (x,y,z) цилиндр ориентирован по оси z, то в обратном пространстве (X,Y,Z), его интенсивность рассеяния имеет вид
I
X
,Y
,
Z
d
sinbZ J1 d
2Z X 2
X2 Y2
Y
2
2
,
(6)
где d – диаметр цилиндра, b – его высота. При переходе в координаты и в формуле (6) следует вы-
брать X 0 , Y sin и Z cos .
В работе рассмотрены два цилиндра с одинаковым радиусом, но разной высотой. Интенсивность рассчитывалась в интервале углов 0,5 мрад 10,5 мрад и 0 / 2 при длине волны
CuKα = 0,1542 нм. Для оценки влияния экспериментальной погрешности на качество коллимационного пересчета в
искаженную (экспериментальную) интенсивность J рассеяния в каждой точке k , k вносился статистический «шум», распределенный по нормальному закону со средним квадратичным отклонением
k ,k J k ,k 104 J 0,0,
причем в качестве общей относительной погрешности интенсивности (в %) выбиралось J , J~ , где J~ –
зашумленная интенсивность. Уровень шума можно менять, умножая интенсивность J во всех точках на один и тот же коэффициент. Для зашумленных («псевдоэкспериментальных») интенсивностей, получен-
ных с разным уровнем шума J , J~ , рассчитывались обоими методами восстановленные точечные
I~(,) интенсивности, и на их основе с использованием формул (2) или (5) рассчитывались соответствующие восстановленные искаженные интенсивности рассеяния J~(,) .
Параметры цилиндра: d = 20 нм, b = 40 нм. Исходя из выражения (6), была рассчитана точечная интенсивность рассеяния цилиндра I (,) (рис. 1, а), а по формуле (1) – искаженная интенсивность
J (,) , в которую был внесен «шум» в 1% и 4%, и получены зашумленные («псевдоэкспериментальные») интенсивности J~(,) (рис. 2, а). На приводимых ниже рисунках представлены интенсивности,
относящиеся только к шуму в 4%.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 4 (74)
79
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ… МБФ. По формулам (2) и (3) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 1, b). Далее на основе полученных функций I~(,) по формуле (1) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) как для шума в 1%, так и для шума в 4% (рис. 2, b). ММФ. По формуле (5) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) (рис. 2, а) по формуле (5) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 1, с). Далее на основе полученных функций I~(,) по формуле (5) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) для шума в 1% и 4% (рис. 2, с).
abc Рис. 1. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 40 нм: (а) – точечная интенсивность I (,) ; (b) – восстановленная точечная интенсивность I~(,) ,
МБФ, «шум» в 4%; (с) – восстановленная точечная интенсивность I~(,) , ММФ, «шум» в 4%
abc Рис. 2. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 40 нм: (а) – зашумленная («псевдоэкспериментальная») интенсивность J~(,) ; (b) – восстановленная искаженная интенсивность J~(,) , МБФ, «шум» в 4%; (с) – восстановленная искаженная интенсивность J~(,) ,ММФ, «шум» в 4%
Параметры цилиндра: d = 20 нм, b = 80 нм. Исходя из выражения (6), была рассчитана точечная интенсивность рассеяния цилиндра I (,) (рис. 3, а), а по формуле (1) – искаженная интенсивность J (,) , в которую был внесен «шум» в 1% и 4%, и получены зашумленные («псевдоэкспериментальные») интенсивности J~(,) (рис. 4, а). На приводимых ниже рисунках представлены интенсивности, относящиеся только к шуму в 4%.
МБФ. По формулам (2) и (3) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 3, b). Далее на основе полученных функций I~(,) по формуле (1) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) как для шума в 1%, так и для шума в 4% (рис. 4, b).
80 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 4 (74)
А.В. Смирнов, Д.Д. Захаров, Б.А. Федоров
ММФ. По формуле (5) для J (,) были рассчитаны восстановленная точечная I(,) и на ее основе – восстановленная искаженная J(,) интенсивности рассеяния. Для зашумленных интенсивностей J~(,) по формуле (5) рассчитывались соответствующие восстановленные точечные интенсивности I~(,) (рис. 3, с). Далее на основе полученных функций I~(,) по формуле (5) вновь рассчитывались восстановленные искаженные интенсивности рассеяния J~(,) для шума в 1% и 4% (рис. 4, с).
abc Рис. 3. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 80 нм: (а) – точечная интенсивность I (,) ; (b) – восстановленная точечная интенсивность I~(,) ,
МБФ, «шум» в 4%; (с) – восстановленная точечная интенсивность I~(,) , ММФ, «шум» в 4%
abc Рис. 4. Интенсивности рассеяния однородным ориентированным цилиндром с параметрами d = 20 нм,
b = 80 нм: (а) – зашумленная («псевдоэкспериментальная») интенсивность J~(,) ; (b) – восстановленная искаженная интенсивность J~(,) , МБФ, «шум» в 4%; (с) – восстановленная искаженная интенсивность J~(,) , ММФ, «шум» в 4%
Заключение
Результаты, отражающие качество коллимационного пересчета, представлены в таблице.
Уровень шума 1
J, J~, %
0 1,0 4,0
0 1,0 4,0
Метод базисных функций
Модифицированный метод Фридмана
Расхождение результатов
23
45 67
8
I, I~, % J~, J~ ,% J, J~,% I, I~,%
J~, J~ ,%
J, J~ ,%
I~М БФ, I~М МФ ,%
Однородный ориентированный цилиндр с параметрами d = 20 нм, b = 40 нм
1,0 0,6 0,6
3,3 0,8 0,08
2,6
1,4 1,1 0,6
3,1 0,9 0,3
2,9
2,5 4,0 1,0
3,4 3,8 0,9
2,8
Однородный ориентированный цилиндр с параметрами d = 20 нм, b = 80 нм
2,7 0,8 0,8
3,8 0,3 0,3
2,6
3,4 1,4 1,1
3,9 0,9 0,4
5,3
5,6 4,1 1,5
4,8 3,8 1,1
5,2
Таблица. Относительное расхождение между интенсивностями: точечной и точечной восстановленной; искаженной и искаженной восстановленной для рассеяния на однородном ориентированном цилиндре
при различном уровне вносимого «шума»
Из анализа этой таблицы можно сделать следующие заключения.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 4 (74)
81
АНАЛИЗ И УЧЕТ КОЛЛИМАЦИОННЫХ ИСКАЖЕНИЙ…
1. Погрешность восстановления как точечных, так и искаженных интенсивностей рассеяния нигде не превышает нескольких процентов даже с учетом «псевдоэкспериментальной» зашумленности.
2. Оба метода дают приблизительно одинаковый результат внесения коллимационных поправок, хотя, как можно видеть из сравнения столбцов 2 и 5, МБФ «работает» несколько лучше, если «шум» не слишком большой и цилиндры не слишком вытянутые. Из сравнения точечных интенсивностей рассеяния, восстановленных методами МБФ и ММФ (столбец 8), следует, что результаты применения обоих методов согласуются с точностью в несколько процентов.
3. Зашумление искаженной интенсивности, т.е. переход от J к J~ , естественно, приводит к появлению расхождения между этими функциями (столбец 1). Интересно, однако, отметить, что если эту зашумленную искаженную интенсивность J~ «пропустить» через двойной коллимационный пересчет J~ I~ J~ , то полученная восстановленная искаженная интенсивность J~ оказывается заметно ближе к незашумленной искаженной интенсивности J , чем просто зашумленная искаженная интенсивность J~ : ср. столбцы 1 и 4 (МБФ) и 1 и 7 (ММФ). Эту закономерность можно также проследить при сравнении рис. 2, а–с, (цилиндр 2040) и рис. 4, а–с, (цилиндр 2080). Как видно, поверхности J~ на рис. 2, а, и рис. 4, а, оказываются существенно более шероховатыми по сравнению с поверхностями J~ на рис. 2, b, с, и рис. 4, b, с, соответственно.
Литература
1. Симаков А.П., Федоров Б.А., Смирнов А.В. Исследование поверхности тонкодисперсных фракций ряда минералов с помощью рентгеновского малоуглового рассеяния // Научно-технический вестник СПбГУ ИТМО. – 2009. – № 2(60). – С. 59–66.
2. Kratky O. Instrumentation, Experimental Technique, Slit Collimation, in Small-angle X-ray scattering. – London: Academic Press, 1983. – Р. 53–84.
3. Федоров Б.А. Учет коллимационных искажений при малоугловом рассеянии рентгеновых лучей. Поправка на высоту щелей // Кристаллография. – 1968. – Т. 13. – № 5. – С. 763–769.
4. Schelten J., Hossfeld F. Application of spline functions to the correction of resolution errors in small angle scattering // J. Appl. Cryst. – 1971. – V. 4. – Р. 210–223.
5. Фомичева Е.Е., Темнов Д.Э., Смирнов А.В., Федоров Б.А. Влияние дисперсионного наполнителя на основе алюминия на структуру и свойства полипропилена // Научно-технический вестник СПбГУ ИТМО. – 2009. – № 6(64). – С. 59–65.
6. Захаров Д.Д., Сизиков В.С., Смирнов А.В., Федоров Б.А. Решение двухмерной коллимационной задачи рассеяния рентгеновских лучей с использованием нестандартных интегральных уравнений // Научно-технический вестник СПбГУ ИТМО. – 2006. – № 32. – С. 144–153.
7. Захаров Д.Д., Смирнов А.В., Федоров Б.А. Решение коллимационной задачи при малоугловом рентгеновском рассеянии на анизотропных объектах // Наносистемы: физика, химия, математика. – 2011. – Т. 2. – № 3. – С. 26–44.
8. Glatter O. Data Evaluation in Small Angle Scattering: Calculation Electron Density Distribution by Means of Indirect Fourier Transformation // Acta Physica Austriaca. – 1977. – V. 47. – Р. 83–102.
9. Верлань А.Ф., Сизиков В.С. Интегральные уравнения: методы, алгоритмы, программы. – Киев: Наукова думка, 1986. – 543 с.
10. Greville T.N.E. Introduction to spline functions. In Theory and Applications of Spline Functions. – New
York: Academic Press, 1969. – Р. 1–35.
11. Вайнштейн Б.К. Дифракция рентгеновских лучей на цепных молекулах. – М.: Изд-во АН СССР, 1963. – 372 с.
Смирнов Александр Витальевич
Захаров Денис Дмитриевич Федоров Борис Александрович
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат физ.-мат. наук, доцент,
smirnav@phd.ifmo.ru – Санкт-Петербургский государственный университет информационных тех-
нологий, механики и оптики, аспирант, denzakharov@gmail.com – Санкт-Петербургский государственный университет информационных тех-
нологий, механики и оптики, доктор физ.-мат. наук, профессор,
borfedorov@rambler.ru
82 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 4 (74)