ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ЭЛЕМЕНТАРНОЙ

ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ЭЛЕМЕНТАРНОЙ
ЯЧЕЙКИ КРИСТАЛЛОВ
Учебно-методические указания к выполнению лабораторной работы
по курсу "Рентгеноструктурный анализ"
ОМСК
2004
Министерство образования и науки РФ
Омский государственный университет
ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ЭЛЕМЕНТАРНОЙ
ЯЧЕЙКИ КРИСТАЛЛОВ
Учебно-методические указания к выполнению лабораторной работы
по курсу "Рентгеноструктурный анализ"
___________________________________________________________________
Издание ОмГУ
Омск 2004
УДК: 539.26
Панова Т.В., Блинов В.И. Определение параметров элементарной ячейки
кристаллов: Учебно-методические указания к выполнению лабораторной работы
по курсу "Рентгеноструктурный анализ". Омск, 2004. 15с.
В работе даются основы методов прецизионных измерений параметров
элементарных ячеек поликристаллов.
Работа утверждена
в качестве учебно-методических указаний к
выполнению лабораторной работы на Ученом Совете физического факультета
23 апреля 2004 г.
 Омский госуниверситет, 2004.
Цель работы: ознакомиться с методами прецизионных измерений
параметров элементарных ячеек поликристаллов. Прецизионное определение
параметров элементарной ячейки с помощью экстраполяционного метода
наименьших квадратов.
Принадлежности: рентгеновский аппарат Дрон-3М; ионные кристаллы,
металлы ; ПЭВМ "PENTIUM".
ВВЕДЕНИЕ
Общие сведения
Для решения ряда задач, связанных с изучением твердого тела,
необходимо предельно точное определение периодов кристаллической
решетки /1/. Они зависят от температуры, от концентрации примеси,
напряжений, возникающих при упругой деформации. Измеряя с большой
точностью периоды решетки при постоянной температуре, можно определить
содержание растворенного элемента в твердом растворе, структурный тип
твердого раствора, измерить упругие напряжения в материале. Сопоставляя
периоды решетки одного и того же вещества, измеренные при разных
температурах,
находят
коэффициенты
термического
расширения.
По
периодам решетки кристаллов, закаленных с высоких температур, можно
оценить концентрацию вакансий при температуре нагрева под закалку.
Анализируя изменение периодов пересыщенного твердого раствора при его
распаде, можно установить закономерности кинетики этого процесса,
вызывающие существенные изменения свойств сплава. Это далеко не полный
перечень задач, которые можно решать путем точных измерений периодов
решетки.
Параметры ячейки определяются путем измерения межплоскостных
расстояний для ряда линий с известными индексами отражения hkl /2/. Число
линий должно быть, по крайней мере, равно числу неизвестных параметров.
Однако часто достаточная информация (например, о составе твердого
раствора) может быть получена путем измерений межплоскостного
расстояния по какой-либо одной линии. В любом случае важно измерить
межплоскостные расстояния с возможно большей точностью. Зная точное
значение межплоскостного расстояния и, воспользовавшись квадратичными
формами, связывающими d2 и периоды решеток (таблица 1 /5/), можно
вычислить периоды решеток для всех кристаллографических систем.
Таблица 1
Сингония
Межплоскостные расстояния dHKL
Кубическая
1
H 2  K 2  L2

d2
a2
Тетрагональная
1
H 2  K 2 L2

 2
2
2
c
d
a
Ромбическая
1
H 2 K 2 L2


 2
d2
a 2 b2 c
Ромбоэдрическая
1
d
Гексагональная
Моноклинная
Триклинная
Межплоскостные

(H
2
2
2
2
2
 K  L ) sin   2( HK  KL  HL )(cos   cos )
2
2
3
a (1  3 cos   2 cos  )
1
4 ( H 2  KH  K 2 ) L2
 
 2
c
d2 3
a2
1
H2
K2
L2
2 HL cos 




2
2
2
2
2
2
d
a sin  b
c sin 
ac sin 2 
1
1
 2 [ s11H 2  s22 K 2  s33L2  s12 HK  s23 KL  s13 HL ]
2
d
V
расстояния
определяют
экспериментально,
индицирование рентгенограммы (то есть приписывание индексов
а
hkl
дифракционным пикам) известных структурных типов проводят путем
сопоставления экспериментально полученной рентгенограммы с литературными
данными или при расчете теоретической дифрактограммы. Индицирование
неизвестных структур по порошковым данным является очень сложной и не
всегда однозначно решаемой задачей для специального исследования. Для
кристаллов всех сингоний, кроме кубической, межплоскостные расстояния в
общем случае зависят от всех линейных параметров решетки, и для определения
периодов необходимо использовать не менее стольких линий, сколько
различных линейных параметров в решетке данной сингонии. Для расчета
параметров ячейки выбирают соответствующее число проиндицированных,
наиболее
четких,
неперекрывающихся,
достаточно
интенсивных
линий.
Желательно выбирать линии в области больших углов, т.к. при одной и той же
точности
измерения

углов
точность
определения
межплоскостных
расстояний возрастает с увеличением угла отражения в соответствии с
уравнением:
d
 ctg  
d
d
d
Ниже приведены значения
при  = 0,001 радиан для измерений,
проведенных по линиям с разными углами:
0
20
40
50
60
70
80
85
d
d
0,275
0,12
0,084
0,058
0,036
0,018
0,009
Высокой точности определения периодов (погрешность 0,01- 0,001 %) можно
достигнуть, применяя особые методы съемки и обработки результатов
измерения рентгенограмм, так называемые прецизионные методы. Достижение
максимальной
точности
в
определении
периодов
решетки
возможно
следующими методами:
1) использование значений межплоскостных расстояний, определенных из
углов в прецизионной области;
2)
уменьшение
погрешности
в
результате
применения
точной
экспериментальной техники;
3) использование методов графической или аналитической экстраполяции.
Минимальная погрешность d/d получается при измерениях под углами  =
80830. К сожалению, далеко не все вещества дают на рентгенограмме линии
под такими большими углами. В этом случае для измерений следует
использовать линию под возможно большим углом . Увеличение точности
определения параметров ячейки связано также с уменьшением случайных
ошибок, которые можно учесть только усреднением, и с учетом
систематических погрешностей, если известны причины их возникновения.
Учет систематических погрешностей при определении параметров
решетки
сводится
к
нахождению
зависимости
систематических
погрешностей от брэгговского угла , что позволяет провести экстраполяцию
к углам  = 900, при которых погрешность определения межплоскостных
расстояний становится малой. К случайным погрешностям относятся
погрешности
измерения,
связанные
с
определением
положения
дифракционной линии на рентгенограмме. К систематическим относятся
погрешности, обусловленные геометрией съемки и физическими факторами.
Методика измерения углов отражения при дифрактометрической
регистрации
В дифрактометре щель счетчика движется по окружности. При
пересечении щелью счетчика конуса дифрагированных лучей возникает пик на
дифракционной диаграмме. На дифракционной диаграмме угол 2 является
углом поворота счетчика - отсчитывается как координата максимума пика по
оси абсцисс. Положение дифракционной линии может быть определено с
большой точностью, если ее запись проведена по точкам. Координата линии
определяется положением либо ее максимума, либо центра тяжести. В случае
симметричных линий максимум находится следующим образом. На записанном
профиле дифракционной линии проводят ряд прямых, соединяющих точки с
равной интенсивностью (измеряемой от линии фона) и лежащих по разные
стороны максимума. Полученные отрезки делят пополам и через середины
проводят прямую до пересечения с профилем линии, которое и определяет
положение максимума (2max) (рис. 1).
Рис. 1. К методике определения центра тяжести и максимума дифракционного
пика
Изменение фона вдоль линии приводит к смещению положения
максимума в сторону больших . Для асимметричных линий лучше определять
центр тяжести по формуле:

2
 I ( )    d

с = 1
,

2
 I ( )d

1
где I() - функция, описывающая профиль дифракционного пика; 1 и 2 границы сечения пика.
Практически с определяют следующим образом. Проводят линию фона.
Находят точки 21 и 22, в которых профиль линии сливается с линией фона.
Затем отрезок от 21 до 22 разбивают на n равных участков 2. Тогда:
 n
  k  Ik
2 c   k n1

  Ik
 k 1


 2  2 ,
1



где 2=(22-21)/(n-1), k – номер участка, Ik – интенсивность на данном
участке, промеренная от линии фона. Погрешность определения центра тяжести
связана как с точностью измерения интенсивности, так и с присутствием фона
на дифрактограмме. Более точно с учетом фона положение центра тяжести
можно вычислить по формуле:
n
 kI
 2 k  1 k rb0
2  2 


 1,
c
z 1 r
n
 2
 I
k
k 1
где b0 
n  1 2a  z

- центр тяжести трапециедального фона; r – отношение
3
az
интегрального фона к интегральной интенсивности:
r
(n / 2)( a  z )
,
n
 Ik
k 1
а- интенсивность при 2а; z- интенсивность при 2z (измерение
интенсивности ведется от нулевого уровня) (рис.2).
Рис. 2. Определение центра тяжести интерференционной линии с учетом
фона
Систематические погрешности можно учесть различными способами.
Рассмотрим некоторые из них:
1. Точное проведение эксперимента. Геометрические условия съемки
дифракционной линии нужно выбрать такие, чтобы величина систематической
погрешности не превышала случайной погрешности измерений. Для расчета
необходимо использовать дифракционные пики в прецизионной области. При
всей привлекательности этот метод имеет очень существенный недостаток, так
как
приводит
к
резкому
снижению
светосилы
дифрактометра
и
соответствующему увеличению экспозиций.
2. Аналитический метод. Наиболее просто систематические погрешности
могут быть учтены для центра тяжести дифракционного пика:
(с)ист= (с)изм+ сi,
где сi- положение центра тяжести i-й инструментальной функции. Чтобы
воспользоваться
этим
уравнением,
необходимо
точно
знать
значения
геометрических параметров съемки, что не всегда возможно. Неидеально
плоская поверхность образцов приводит к небольшому смещению центра
тяжести (до 1 при  30о), которое нельзя рассчитать заранее.
3. Экстраполяционные методы.
Эти
методы
применимы
главным
образом
к
высокосимметричным
веществам, относящимся к кубической, гексагональной или тетрагональной
сингониям /4/. Для кубических кристаллов с параметром элементарной ячейки а
можно записать:
 а/а0= -ctg  i =f(),
где i- сумма инструментальных смещений дифракционной линии. Тогда
величина параметра элементарной ячейки равна: а = а0[1+f()]. При  900,
ctg0 и, следовательно, f() 0. Кроме того, f() можно представить в виде
некоторой простой функции от угла . Например, для центров тяжести пиков
под углами  600 с достаточной точностью можно положить f()= cos2. Тогда
величина параметра а является линейной функцией сos2. Таким образом, нужно
точно измерить положение нескольких пиков под углами  600 (последний из
них должен иметь угол  = 78-820), для каждого из них определить значение
параметра элементарной ячейки а и отложить эти значения в зависимости от
cos2.
Пересечение
прямой,
соединяющей
экспериментальные
значения
параметра а с осью 900, определит значение параметра аэкстр, свободное от всех
систематических погрешностей, исключая погрешности из-за вертикальной
расходимости.
4. Использование эталонного
вещества. Этим методом получают
рентгенограмму смеси исследуемого вещества с эталоном (в качестве эталона
можно использовать NaCl, Si, CaF2), для которого известна величина параметров
элементарной ячейки. Для эталонного вещества рассчитывают величины
межплоскостных расстояний и затем величины углов отражения ’ист, которые
не будут совпадать с соответствующими ’изм, полученными из рентгенограммы.
Ясно, что истинные значения углов отражения для исследуемого вещества ист0
можно получить таким образом:
 ист0=0изм+ист - изм.
Важно лишь, чтобы углы  для линии исследуемого вещества и эталона были
близкими друг к другу. Это требование важно, потому что величина
систематических погрешностей зависит от угла , а для сравниваемых линий
исследуемого вещества и эталона она должна быть примерно равной.
Систематические смещения дифракционного пика вызываются и рядом
физических факторов, среди которых наиболее существенны изменения
дисперсии и множителя Лоренца, поляризации по ширине дифракционного
пика, изменение температуры образца во время съемки и преломление лучей.
Изменение температуры образца во время съемки приводит в результате
термического расширения к изменению параметров элементарной ячейки.
Учесть влияние колебаний температуры можно, используя специальную
термостатирующую приставку. Погрешности, связанные с вертикальной
расходимостью, дисперсией, множителем Лоренца и поляризацией, искажают
профиль дифракционного пика, однако значения их ничтожно малы и ими
можно пренебречь. Для получения максимальной точности необходимо вводить
поправку на преломление. Влияние преломления лучей таково. Как следует из
общей теории интерференции, в уравнение Вульфа-Брэгга должна входить
длина волны лучей, распространяющихся внутри исследуемого вещества, а не в
вакууме или воздухе, тогда как именно последняя фигурирует в таблицах и
используется в расчетах. Эта погрешность не устраняется экстраполяцией, так
как не зависит от угла . Погрешность, обусловленная преломлением, невелика
(0,003%), но если точность определения периода решетки превосходит 2 10 -4А,
то поправкой на преломление пренебрегать нельзя.
Современные
соответствующими
уточнения
рентгенографические
программами
параметров
решетки
для
методы
индицирования
методом
обеспечены
дифрактограмм
наименьших
квадратов,
и
что
существенно повышает достоверность результатов, т.к. при этом могут быть
использованы все измеренные рефлексы дифрактограммы.
Порядок выполнения работы
1.
числа
Подготовить образцы для съемки на дифрактометре ДРОН-3М (из
предложенных
преподавателем).
Ионные
кристаллы
тщательно
измельчить в ступке и завальцевать в кювету. Металлические образцы
отшлифовать и отполировать с одной стороны.
2.
Установить образец в центре держателя образцов, строго соблюдая
плоскопараллельность. С помощью преподавателя прописать дифракционные
спектры по точкам для образцов с кубической (ионные кристаллы) и
гексагональной (металлы) сингонией. Записать полученные данные на дискету,
создав два разных файла.
3.
С
помощью
программы
ORIGIN
или
FPEAK
построить
дифрактограммы для обоих образцов, сохранив их с соответствующим
расширением.
4.
Пользуясь таблицами приложения / 5 / или программой PowderCell
провести индицирование всех дифракционных пиков, присутствующих в
спектрах.
5.
Рассчитать межплоскостное расстояние для той длины волны
рентгеновского излучения, которая была использована при съемке.
6.
Рассчитать
параметры
элементарной
ячейки,
используя
квадратичные зависимости (таблица 1). Использовать пики, находящиеся в
прецизионной области или наиболее приближенные к ней. Для гексагональной
сингонии для расчета использовать пару дифракционных пиков с разными
индексами Миллера. В этом случае расчет производить по формулам,
получаемым из квадратичных зависимостей:
a2 
4

3
AB A B
1 2
2 1
;
B
B
2
1

2
2
d
d
H K L
H K L
1 11
2 2 2
c2 
AB A B
1 2
2 1
,
A
A
1
2

2
2
d
d
H 2K L
H K L
2 2
1 11
где А=H2+HK+K2 и B=L2.
7.
Провести
аппроксимацию
периода
решетки
а
для
ионных
кристаллов с кубической сингонией, путем подбора соответствующей функции с
помощью любого графического редактора (Origin или Excel). При этом расчет
параметра а произвести для всех пиков, относящихся к кубической сингонии.
Построение
проводить в последовательности
(а от cos 2) и т.д. Перебор
функций производить до тех пор пока не будет достигнуто максимальное
приближение к линейной
зависимости. После подбора функции определить
параметр а при угле 900.
8.
Рассчитать погрешности измерений.
9.
Произвести уточнение параметров элементарной ячейки методом
наименьших квадратов с помощью программы UnitCell. Запуск программы
осуществляется при загрузке файла UnitCell.exe. Для расчета предварительно
создать файл с расширением .dat
или .txt (пример созданием такого файла
приведен в описании программы readme.txt).
10.
Сравнить полученные данные с данными, рассчитанными без
уточнения.
Контрольные вопросы
1. Какие характеристики дифрактограммы необходимо брать для расчета
параметра решетки?
2. Как определяется угол при дифрактометрической регистрации спектра?
3. Какие существуют методики определения параметров элементарной ячейки?
4. Как учесть систематические и физические погрешности?
5. Каковы возможности программы UnitCell?
Рекомендуемый библиографический список
1. Русаков А.А. Рентгенография металлов: Учебник для вузов. М.: Атомиздат,
1977. 480 с.
2. Миркин
Л.И.
Рентгеноструктурный
контроль
машиностроительных
материалов: Справочник. М.: Изд-во МГУ, 1976. 140 с.
3. Миркин Л.И. Справочник по рентгеноструктурному анализу поликристаллов.
М.: Физматгиз, 1961. 863 с.
4. Уманский Я.С., Скаков Ю.А., Иванов А.Н., Расторгуев Л.Н. Кристаллография,
рентгенография и электронная микроскопия. М.: Металлургия, 1982. 632 с.
5. Горелик С.С., Расторгуев Л.Н., Скаков Ю.А. Рентгенографический и
электроннооптический анализ. Приложения. М: Металлургия, 1970. 107 с.
Панова Татьяна Викторовна
Блинов Василий Иванович
ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ЭЛЕМЕНТАРНОЙ
ЯЧЕЙКИ КРИСТАЛЛОВ
Учебно-методические указания к выполнению лабораторной работы
по курсу «Рентгеноструктурный анализ»
Редактор
Подписано в печать
Формат бумаги 6084 1/16
Печ.л. . Уч.-изд. л. .
Тираж
Издательско-полиграфический отдел ОмГУ
644077, Омск - 77, пр. Мира 55-а, госуниверситет
Заказ