Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Калин Физическое материаловедение Том 1 Физика твердого тела 2007.pdf
Скачиваний:
2377
Добавлен:
16.08.2013
Размер:
7.64 Mб
Скачать

+∞ +∞ +∞

 

Ш(r) = ∑ ∑ ∑δ(r Rmnp ) ,

(1.54)

m=−∞ n=−∞ p=−∞

 

+∞ +∞ +∞

 

Ш(h) = ∑ ∑ ∑δ(h Hhkl ) ,

(1.55)

h=−∞ k =−∞ l =−∞

где r и h – текущие векторы в прямом и обратном пространствах соответственно.

Функции Ш(r) и Ш(h) связаны взаимным преобразованием Фурье, поэтому обратное пространство называют фурье-пространст- вом. Так, при рассмотрении дифракции рентгеновских лучей на кристаллах в прямом пространстве задается электронная плотность ρ(r), а в обратном – амплитуда рассеянного рентгеновского излучения A(h). По этой причине обратное пространство называют ди-

фракционным пространством.

1.2.5. Матрица ортогонального преобразования

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

 

 

3

 

 

 

aik ail = δkl (k, l = 1, 2, 3),

(1.56)

 

 

i =1

 

где aij − элементы матрицы A, а δkl = 1 при k = l и δkl

= 0 при k ≠ l.

Если A − матрица ортогонального преобразования, то

detA = ±1 и

A

–1

~

 

=

A . Матрица A с detA = +1 соответствует собственному пово-

роту, а с detA = −1 − несобственному повороту.

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

Таким образом, собственные значения матрицы A должны иметь вид: 1, exp(−iα), exp(+iα). При этом фаза α комплексного собственного значения называется углом поворота, а собственный вектор u с собственным значением λ = 1 называется осью поворота.

51

Угол поворота α находят из линейного инварианта матрицы поворота A, т.е. из ее следа (обозначают trA от английского слова

«trace» или SpA − от немецкого «Spur»)

trA = a11 + a22 + a33 = 1 + exp(−iα) + exp(+iα) = 1 + 2cosα. (1.57)

Компоненты оси поворота определяются из элементов матрицы поворота A как

u1 : u2 : u3 = (a32 a23) : (a13 a31) : (a21 a12).

(1.58)

Направляющие косинусы оси вращения u определяются как

 

l = u (u2

+u2

+u2 )

1

 

 

2

 

 

 

1

1

1

2

3

 

1

 

 

 

 

 

 

(u2

+u2

+ u2 )

 

 

l

2

= u

2

2 .

(1.59)

 

 

 

1

2

3

 

1

 

 

 

 

 

 

(u2

+u2

+u2 )

 

 

l

3

= u

3

2

 

 

 

 

 

 

 

1

2

3

 

 

 

 

В общем случае матрица поворота на угол φ вокруг оси u с направляющими косинусами l1, l2, l3 имеет вид

 

 

 

 

 

 

R(u, ϕ) =

 

 

l2

(1cosϕ) +cosϕ

l l

(1cosϕ) l cosϕ

 

 

1

 

 

1 2

3

=

 

 

(1cosϕ) +l3sinϕ

2

(1cosϕ) +cosϕ

l1l2

l2

 

l l

(1cosϕ) l

sinϕ

l l

(1cosϕ) +l sinϕ

 

 

1 3

2

 

2 3

1

l1l3(1cosϕ) +l2sinϕ

l2l3(1cosϕ) l1sinϕ .(1.60) l32 (1cosϕ) +cosϕ

Если поворот задан эйлеровскими углами φ, θ, ψ, то каждая из матриц поворота имеет простой вид в соответствующей повернутой системе координат (рис. 1.24). Так матрицы первого поворота на угол φ вокруг оси z Rz(φ), второго поворота на угол θ вокруг оси x′ − Rx′(θ) и третьего поворота на угол ψ вокруг оси z′ − Rz′(ψ) имеют вид:

cosϕ

sinϕ

0

 

1

0

0

Rz(φ) = sinϕ

cosϕ

0

,

Rx(θ) =

0

cosθ

sinθ ,

 

0

1

 

 

0

sinθ

 

0

 

 

cosθ

 

 

cosψ

sinψ

0

 

 

 

 

Rz(ψ) = sinψ

cosψ

0

.

 

 

 

 

 

0

0

1

 

 

 

 

 

 

 

 

 

Матрица поворота R(φ, θ, ψ) обычно записывается в неподвижной системе координат xyz, причем если ее выражать через матри-

52

цы Rz(φ), RX′(θ) и RZ′(ψ), то порядок их перемножения обратен по-

рядку выполняемых поворотов, т.е. R(φ, θ, ψ) = RZ(φ) RX′(θ) RZ′(ψ), так что

R(ϕ,θ,ψ) =

 

cosϕcosψ −sinϕsinψ

cosϕsinψ −cosθsinϕcosψ

sinθsinϕ

 

 

 

sinϕsinψ + cosθcosϕcosψ

 

 

.(1.61)

= sinϕcosψ + cosθcosϕsinψ

sinθcosϕ

 

sinθsinψ

sinθcosψ

cosθ

 

 

 

 

 

Углы φ и ψ изменяются от 0 до 2π, а угол θ − от 0 до π. Эйлеровские углы поворота часто используются в физике твер-

дого тела, в кристаллографии интересуются кристаллографическими индексами оси поворота и значением угла поворота.

В результате поворота, определяемого матрицей A, произвольный вектор r переходит в r′, что в матричном виде записывается как

r= Ar,

т.е.

 

 

 

r

a

a

1

 

11

12

r2

 

= a21 a22

 

 

 

a32

r3

 

a31

a13 a23 a33

(1.62)

r1

r2 , (1.63) r3

где aij − элементы матрицы A, (r1, r2, r3) − компоненты вектора r,

(r1, r2, r3) − компоненты вектора r′.

Вообще говоря, матрица A, как и любая матрица ортогонального преобразования, допускает двойную интерпретацию: «активную» и «пассивную». В первом случае матрица A поворачивает вектор r на угол α вокруг оси u, переводя его в r′, причем оба вектора задаются в одной системе координат. Во втором случае − при пассивной интерпретации − можно считать, что матрица A поворачивает систему координат вокруг оси u на угол α в противоположном направлении, причем компоненты вектора r не изменяются (рис. 1.36).

Компоненты матрицы A легко определить, исходя из соотношения (1.63). Действительно, направление [100] с координатами (r1, r2, r3) после поворота имеет координаты (r1, r2, r3). Таким образом, матрицу A можно представить состоящей из матрицстолбцов a*1, a*2, a*3, компонентами которых являются координаты векторов [100], [010], [001] после поворота.

53

Рис. 1.36. Активная (а) и пассивная (б) интерпретация поворота на угол α вокруг оси x

Если матрица A задает операцию в системе координат xyz

r2 = Ar1,

(1.64)

то в системе координат x'y'z' эта операция задается матрицей A'

r'2 = A'r'1.

Переход из системы координат xyz в систему координат x'y'z' осу-

ществляется матрицей C, т.е.

r'1 = Cr1 и r'2 = Cr2, r1 = C–1r'1 и r2 = C–1r'2.

Подставляя r1 и r2 в (1.64), получим C–1r'2 = AC–1r'1 или, умножая слева и справа на C, имеем r'2 = CAC–1 r'1. Матрица CAC–1 = A' осуществляет в системе координат x'y'z' такую же операцию, как и

A в системе xyz. Матрицы A и A', связанные соотношением

 

A' = CAC–1

(1.65)

или, если D = C–1, A' = D–1AD, называют матрицами подобия. У матриц подобия равны все инварианты, в том числе и линейный инвариант (след матрицы), т.е.

trA = trA' = 1 + 2 cosα.

(1.66)

1.2.6.Преобразование индексов направлений

иплоскостей при изменении системы координат

Произвольное направление [mnp] задается вектором Rmnp = ma1 + + na2 + pa3 или вектором Rт1т2 т3 = m1a1 + m2a2 + m3a3. В матричном виде Rт1т2 т3 Rm:

54

( ) m1 = ~
Rm = a1 a2 a3 m2 am ,
m3

(1.67)

где ~ − матрица-строка из базисных векторов, а − матрица- a m

столбец из компонентов вектора.

При переходе из системы координат xyz к системе x'y'z' базисные

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

 

 

 

 

a

 

 

A a A a A a

 

 

1

 

 

11 1

12 1

13 1

 

 

a2

 

=

A21a2

A22a2

A23a2

 

, т.е. a' = Aa

 

 

 

A31a3

A32a3

A33a3

 

 

a3

 

 

 

 

или

~~

 

 

~

 

(1.68)

a

′ = aA .

 

 

~

= Rm

Вектор Rm в новой системе координат переходит в Rm' = a m

 

и с учетом (1.68)

 

 

 

~~

~

 

(1.69)

aAm′ = am ,

 

~

 

 

 

откуда Am′ = m и

~

 

 

 

 

(1.70)

m′ = A1 m.

 

Таким образом, при изменении системы координат базисные векторы преобразуются как ковариантные векторы с помощью

матрицы

~

 

 

 

 

 

 

 

 

A (1.68), а компоненты вектора R преобразуются как

контравариантные векторы с помощью матрицы

~

1

. Иногда

A

 

матрицу

A~1 называют матрицей контраградиентного преобразо-

вания.

 

 

 

 

 

 

 

 

 

 

Скалярное произведение между векторами Rт т

т

и Hhkl равно

 

 

 

 

 

 

 

1

2

 

3

 

 

Rт т

т

Hhkl = m1h + m2k + m3 l или в матричной записи

 

 

 

1

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

~

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

(h k l) m2

 

= hm .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m3

 

 

 

 

 

 

 

При изменении системы координат скалярное произведение сохра-

~

=

~′ ′

и с учетом (1.70)

~

~

~

1

m и

~ ~

~

1

,

няетсяhm

h m

hm =

h

A

 

h = h

A

 

откуда

55

~

~

~

(1.71)

h

= h

A .

Следовательно, при изменении системы координат компоненты вектора обратной решетки и базисные векторы преобразуются, как

~

ковариантные векторы с помощью матрицы A .

1.2.7. Основные формулы структурной кристаллографии

Период идентичности. Период идентичности I[mnp] вдоль данного направления [mnp] – расстояние между ближайшими узлами вдоль данной прямой – равен абсолютной величине вектора Rmnp.

Длину вектора проще вычислить через скалярное

произведение

I[mnp] = (Rmnp Rmnp)1/2, записанное в матричном виде.

 

 

 

Переходя от Rmnp к Rт т

т

и учитывая соотношения (1.67) и

 

 

 

1

 

2

 

3

 

 

 

 

(1.71), запишем Rт т

т

в прямом и обратном базисах как

 

1

2

 

3

 

 

 

 

 

 

 

 

 

 

 

m

 

~

h

 

~

 

Rт1т2 т3 = (a1 a2

 

1

 

1

 

(1.72)

a3 ) m2

 

= am = (b1 b2

b3 ) h2

= bh

 

 

 

m

3

 

 

h

 

 

 

 

 

 

 

 

 

 

3

 

 

 

~ ~

иam = bh . Умножая слева левую и правую часть последнего соот-

~

~

~

~

ношения на a, получим aam = abh и, поскольку aa

=G, ab =I, то

 

Gm = h.

~

(1.73)

Аналогично, умножая на b

~

 

соотношение am = bh , получим

 

m = G 1h .

 

(1.74)

Соотношение (1.73) показывает, что использование метрической

матрицы G позволяет при

известных компонентах m вектора

 

Rm m

m

3

в прямом пространстве найти его компоненты h в обратном

1

2

 

 

 

 

 

 

 

 

пространстве. Матрица G–1 аналогично переводит компоненты из

обратного пространства в прямое пространство. (1.74).

 

 

Теперь выражение для периода идентичности имеет вид

 

 

 

 

 

 

 

Im m

m

=

~

(1.75)

 

 

 

 

 

 

mGm .

 

 

 

 

 

 

1

2

 

3

 

 

 

Квадратичная форма.

Длина

вектора обратной

решетки

 

Hhkl

 

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

 

 

 

 

56

 

 

 

 

 

 

 

Hhkl

 

=

~

1

h .

(1.76)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hG

 

Выражение

1

=

 

Hhkl

 

2 = h~G1h

, которое является функцией па-

 

 

dhkl2

 

 

 

 

 

 

 

 

 

 

 

 

 

раметров решетки a, b, c, α, β, γ, называют квадратичной формой соответствующей сингонии.

Квадратичные формы для некоторых сингоний имеют вид: для кубической сингонии

 

 

 

 

1

 

 

 

=

h

2 + k 2 + l 2

,

 

 

(1.77)

 

 

 

dhkl2

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для ромбической сингонии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=

h2

+

k 2

+

l2

 

,

 

 

(1.78)

 

 

dhkl2

a2

b2

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

для гексагональной сингонии

 

 

 

 

+ k 2 )

 

 

 

 

1

 

 

=

 

 

4(h2

+ hk

+

l2

.

(1.79)

 

dhkl2

 

 

 

 

 

 

3a2

 

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Угол между направлениями. Угол φ между направлениями

[m1(1) m2(1)

m3(1)]и [m1(2) m2(2) m3(2)]определяют из соотношения

 

 

 

 

 

 

~

(2)

 

(1)

 

 

 

 

 

 

 

 

 

 

cosφ =

 

 

m

 

 

 

Gm

 

 

 

 

 

 

,

 

(1.80)

 

 

 

 

~ (1)

 

 

(1)

~

(2)

 

(2)

 

 

 

 

 

 

m Gm

 

 

 

m

 

Gm

 

 

 

 

 

~ (1)

~

(2)

− матрицы-строки; m

(1)

, m

(2)

− матрицы-столбцы.

где m

, m

 

 

 

Угол

между плоскостями.

 

Угол

 

 

φ

между

плоскостями

(h1(1) h2(1) h3(1)) и (h1(2) h2(2) h3(2))

равен углу между нормалями к этим

плоскостям:

 

 

 

h~(2)G1h(1)

 

 

 

 

 

 

 

 

 

cosφ =

 

 

 

 

 

 

.

(1.81)

 

 

 

 

~(1)

1

h

(1)

~(2)

1

h

(2)

 

 

 

 

 

h G

 

 

 

h

 

G

 

 

 

 

 

Угол между плоскостью и направлением. Угол между плос-

костью

(h1 h2 h3 )

и направлением

[m1 m2 m3

]

определяют

через

угол φ между этим направлением и нормалью к плоскости:

 

 

 

 

 

~

 

 

 

 

 

 

 

cosφ =

 

hm

 

.

 

 

(1.82)

 

 

~

1

h

~

 

 

 

 

 

 

 

 

 

 

 

 

hG

 

mGm

 

 

 

 

57

Объем элементарной ячейки. Объем элементарной ячейки, построенной на базисных векторах прямой решетки a1, a2, a3, выражается смешанным произведением

V = a1[a2×a3], (1.83)

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

V = detG .

(1.84)

Для обратной решетки, построенной на базисных векторах b1, b2, b3, объем элементарной ячейки определяется аналогичными выражениями

V* = b1[b2×b3] = det G1 .

(1.85)

1.2.8. Области Вороного, ячейки ВигнераЗейтца, зоны Бриллюэна

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

Эта ячейка примитивная, так как содержит всего один узел. Область, определенная таким образом, называется в кристаллографии

областью Вороного, а в физике твердого тела ячейкой Вигнера– Зейтца.

ВГЦК решетке ближайшие узлы расположены вдоль направлений <110>, поэтому ячейка Вигнера–Зейтца огранена плоскостями {110}. Поскольку множитель повторяемости для {110} равен 12, то ячейка Вигнера–Зейтца для ГЦК решетки образована двенадцатью гранями {110}, имеющими форму ромба. Таким образом, ячейка Вигнера–Зейтца для ГЦК решетки является правильным ромбическим додекаэдром (рис. 1.37,а).

Вслучае ОЦК решетки ячейка ВигнераЗейтца огранена плоскостями {111} (правильные шестиугольники) и {100} (квадраты). Множители повторяемости для {111} и {100} равны соответственно 8 и 6, и ячейка ВигнераЗейтца изображается многогранником с

14 гранями (рис. 1.37,б).

58