CC BY
79
УДК 538.9:548.1
СКРЫТАЯ АНИЗОТРОПИЯ, МНОГОМЕРНЫЕ РЕШЕТКИ, ЭЛЕКТРОННЫЕ И КОЛЕБАТЕЛЬНЫЕ СПЕКТРЫ СЛОЖНЫХ КРИСТАЛЛИЧЕСКИХ СОЕДИНЕНИЙ
А. С. Поплавной
LATENT ANISOTROPY, MULTIDIMENSIONAL LATTICES, ELECTRONIC AND VIBRATIONAL SPECTRA OF COMPLEX CRYSTALLINE COMPOUNDS
A. S. Poplavnoi
Работа выполнена при поддержке государственного задания № 2014/64 (2014 - 2016 гг.)
Дан обзор описания сложных кристаллических структур на основе их представления совокупностью под-решеток Браве. При наличии в структуре подрешеток с более высокой симметрией, чем кристаллическая, сложный кристалл обладает скрытой симметрией, или псевдосимметрией, которая проявляется в его физических и физико-химических свойствах. Представлены методы исследования этой скрытой симметрии.
The paper reviews the methods of complex crystal structure symmetry description based on representing by a set of Bravaissublattices. If sublattices in the structure are of higher symmetry than that of crystalline one, the composite crystal is of latent symmetry or pseudo symmetry that is revealed in its physical and physico-chemical properties. The investigation methods of suchsymmetry arepresented.
Ключевые слова: типы Браве, сорта Делоне, подрешетки, многомерные кристаллические пространства, фо-нонные спектры, электронные спектры.
Keywords: Bravais types, Delaunay sorten, sublattices, multidimensional crystal spaces, phonon spectra, electron spectra.
Введение. В случае если все подрешетки кристалла обладают одинаковой симметрией, то его описание дается в трехмерном пространстве с использованием решетки с базисом. В случае же если несколько подрешеток имеют различающую симметрию, то полная симметрия таких кристаллов уже не описывается в трехмерном пространстве, необходимо перейти к пространству с размерностью, равной сумме размерностей пространств подрешеток с различающейся симметрией. Многомерные операторы симметрии строятся на основе трехмерных операторов симметрии под-решеток. В случае если имеются подрешетки из одинаковых химических элементов, но различающиеся симметрией, то появляется не унитарный точечный элемент симметрии, связанный с матрицей трансляционной совместимости этих подрешеток, при действии которого происходит изменение масштаба. Пригодными для описания симметрии кристаллов с под-решетками являются семейства групп в 3п (п - число подрешеток) пространствах, представляющие собой прямое произведение подгрупп, действующих в Я3 взаимно ортогональных инвариантных пространствах [1; 2]. Для кристаллов с двумя подрешетками в Я6 пространствах установлено число пригодных для описания семейств, для которых определены число свободных параметров в метрическом тензоре, порядки максимальных голоэдрических групп в семействах, число классов Браве, характер приводимости, записаны метрические тензоры [3].
Развитый кристаллографический метод описания сложных кристаллических соединений применялся при моделировании их электронных и колебательных спектров. При моделировании электронных спектров кристаллический потенциал представлялся в форме суммы потенциалов, приписываемых отдельным под-решеткам и потенциала, отвечающего за гибридиза-
цию подрешеток. Из-за трансляционной совместимости всех подрешеток с решеткой кристалла зоны Бриллюэна (ЗБ) подрешеток могут быть перестроены в ЗБ кристалла. Размещение энергетических зон подрешеток в ЗБ кристалла приводит к появлению вырождений, обусловленных двумя факторами: а) трансляционные вырождения, возникающие при свертывании спектров подрешеток; б) вырождения, обусловленные более высокой точечной симметрией подре-шеток, если таковая имеется. Эти вырождения снимаются при учете гибридизации состояний [4].
При наличии в кристалле подрешеток с симметрией выше кристаллографической группы соединения в фононных спектрах появляются такие же особенно -сти, как и в электронных, и они могут быть исследованы путем свертывания фононных спектров подре-шеток в ЗБ кристалла. Вместе с тем уравнения динамики кристаллической решетки принципиально отличаются от зонного уравнения. Если электронная подсистема может рассматриваться как единая из-за равенства масс и эквивалентности электронов, то уравнения динамики решетки с самого начала записываются с учетом индивидуальных свойств атомов. Таким образом, возникает многомерное пространство степеней свободы, симметрию которого и нужно исследовать. В силу этого особенности фононных спектров существенно зависят от соотношения масс атомов, входящих в различные подрешетки [5].
Результаты применения метода подрешеток к исследованию особенностей зонных и колебательных спектров многих сложных диэлектрических и полупроводниковых кристаллов, к анализу их химической связи, изучению генезиса спектров из состояний под-решеток представлены в [1]. В настоящей статье обсуждаются методические аспекты кристаллографического подхода с позиции подрешеток, особое внима-
©А. С. Поплавной, 52014- 177
ние уделяется преимуществам метода подрешеток, возможностям получения новых результатов при изучении физических и физико-химических свойств сложных кристаллических соединений, обусловленных естественным учетом «скрытой» симметрии сложного соединения.
1. Анизотропные подрешетки в кубических кристаллах. В сложных кристаллических соединениях можно выделить подрешетки, трансляционная симметрия которых оказывается выше кристаллической, однако точечные группы симметрии каждой подрешетки принадлежат к низшей сингонии, по сравнению с кристаллической точечной группой. При этом точечная группа симметрии полной структуры будет составлена из подгрупп, отвечающих точечной симметрии каждой подрешетки с добавлением элементов симметрии совмещающих низкосимметричные (в смысле точечных групп) подрешетки между собой. Такая структура приводит к специфическим особенностям спектров элементарных возбуждений и содержит потенциальную возможность понижения симметрии кристалла под влиянием колебаний решетки или внешних воздействий. Проиллюстрируем эту общую идею на некоторых примерах.
К кристаллам с химической формулой АВХ3 относят соединения с разной кристаллической структурой, в том числе и переменного состава. Благодаря своим замечательным электрическим и магнитным свойствам особый интерес представляют составы со структурой перовскита. К этому структурному типу относят многие сегнетоэлектрики, которые используются в самой быстро развивающейся отрасли - электронике. Некоторые кристаллы со структурой перов-скита обладают достаточно высокой суперионной проводимостью, прежде всего, это соединения, содержащие в качестве аниона фтор.
Структурный тип перовскита (родоначальник -минерал CaTiO3) хорошо известен (рис. 1). Элементарная ячейка может быть кубической, с ПГ № 221 (Оь). Кубической структурой CaTiO3 обладает при температурах T >1580K, при понижении температуры претерпевает последовательность фазовых переходов и при T <1380 K стабилизируется и описывается ромбической структурой с симметрией D2h16. На рис. 2 представлены многогранники Дирихле-Вороного (МДВ) различных кристаллических модификаций CaTiO3, а также их точных и приближенных подреше-ток.
Са
Щ
а)
у г * и
• •
в)
в)
Рис. 1. Элементарные ячейки: а) кубической, б) моноклинной, в) ромбической модификаций СаТЮ3
Кристаллическая решетка кубической модификации CaTiO3 простая кубическая (рис. 1a) с параметрами a = b = с = 3.79 А составлена из 5 Гс подрешеток: 1 кальция, 1 титана и 3 кислорода. Три Гс подрешетки кислорода расположены таким образом, что любые две из них можно рассматривать как простую тетрагональную (ПТ) Tq подрешетку, обладающую более высокой трансляционной симметрией, чем кубическая решетка кристалла. МДВ кубической модификации CaTiO3 (1) и подрешетки Tq кислорода (2) представлены на рис. 2a. В зависимости от того, какие две подрешетки кислорода объединяются, ось симметрии четвертого порядка будет направлена по той или иной координатной оси. Подобное выделение подрешеток кислорода приводит к кажущейся анизотропии, которая в реальности, конечно, отсутствует, если учесть полную структуру кристалла. Однако это выделение показывает, что в кубической модификации соединений ABX3 имеются «потенциальная анизотропия» и «потенциальная неэквивалентность» атомов кислорода. Реализация этих потенциальных возможностей наблюдалась при исследовании физических и физико-химических свойств систем CsPbF2Cl, CsPbFCl2. Рентгеноструктурный анализ показал для этих соединений структуру типа перовски-та со слабыми ромбическими искажениями.
Моноклинная модификация CaTiO3 (рис. 1б) имеет кристаллическую решетку Гт, ПГС^2, параметры решетки a = b = с = 3.79 А, а = в = Y = 90°. Видно, что кристалл обладает трансляционной симметрией простой кубической решетки с удвоенными периодами во всех направлениях, по сравнению с кубической модификацией. МДВ кристалла (1) на рис. 2б в 8 раз больше МДВ кубической модификации (1) на рис. 2a и (2) на рис. 2б. Фазовый переход в моноклинную модификацию произошел только за счет малых смещений атомов кислорода из симметричных позиций. Подрешетки кальция и титана не претерпели изменений и относятся к МДВ (2) на рис. 2б. Если атомы кислорода сместить в симметричные позиции на величины не более 6 % от параметров решетки, то структура «вернется» в кубическую модификацию с МДВ (2) и (3) на рис. 2б, которые эквивалентны (1) и (2) на рис. 2a.
Ромбическая модификация CaTiO3 (рис. 1в) отвечает кристаллической решетке Г0 с параметрами a = 5.3881 А, b = 5.4471 А, с = 7.6541 А, а = в = Y = 90°, ПГ D2h16. Видно, что параметр с - удвоенный по сравнению с кубической модификацией. Эта структура обладает как точной, так и приближенной внутренней симметрией, которая иллюстрируется на рис. 2в. На этом рисунке (1) - МДВ кристаллической решетки. Четыре подрешетки титана объединяются в одну Г0с подрешетку с пространственным периодом вдоль оси z, вдвое меньшим соответствующего периода кристаллической решетки и совпадающим с периодом кубической решетки, а периоды вдоль осей x и y совпадают с кристаллическими. МДВ для Г0с подрешетки - прямая призма с шестиугольным основанием (2), представляющая собой прямоугольник, у которого два противолежащих угла слегка «срезаны», при этом объем данного МДВ в 4 раза меньше МДВ решетки. Атомы кальция можно сместить на величины ~ 3 % периода
решетки в симметричные позиции так, что они также образуют подрешетку Г0с с МДВ (2). Атомы кислорода разделены на два неэквивалентных множества О1 и О2.
При их смещении на величины ~ 6 % от периода ре-
Г с
0с
МДВ (2), а другая в подрешетку Г0с пространственными периодами вдоль всех осей вдвое меньшей кристаллических. Таким образом, в ромбической модификации CaTiO3 проявляются как неэквивалентность, так и анизотропия ромбического типа атомов кислорода, искусственно введенные для кубической модификации.
Другую интересную группу кубических кристаллов с анизотропными подрешетками представляют серебросодержащие суперионные кристаллы (a-AgI, Ag2X (X = S, Se, Te), Ag3SI, Ag3SBr), которые обладают рекордными для твердых тел катионной проводимостью и коэффициентом диффузии. Некоторые из них (халькогениды серебра), кроме того, имеют интересные полупроводниковые свойства, что резко выделяет их среди классических суперионных проводников, проявляющих в большинстве своем чисто ионную проводимость.
Исследованию твердого электролита AgI посвящено много работ. Структура высокотемпературного йодистого серебра a-AgI от 420К до температуры плавления 828 К описывается пространственной группой
Oh9. Ниже полиморфного превращения первого рода в-AgI принадлежит к структурному типу вюртцита с пространственной группой C6v4. Кристаллическая структура высокотемпературной разупорядоченной а-фазы имеет ОЦК анионную подрешетку, а два катиона серебра в кубической ячейке могут занимать 3 типа позиций: 12d - тетраэдрических, 24h - тригональных и 6b - октаэдрических. Молекулярно-динамические (МД) расчеты показывают, что катионы серебра действительно распределены по этим позициям, при этом основной максимум плотности распределения катионов отвечает 12d позициям, второй по интенсивности -24h позициям, позициям 6b отвечает минимум распределения. По этой причине во многих работах при моделировании динамических и термодинамических свойств a-AgI ограничиваются рассмотрением только
12d позиций для движения катионов Ag+.
На рис. 3 представлена модельная структура a-AgI. Тетраэдрические пустоты в ОЦК решетке можно разделить на 6 ОЦК подрешеток, нумерация которых представлена на рис. 3. Анализ структуры этих 6 под-решеток показал, что из них можно образовать 3 ПТ подрешетки с более высокой трансляционной симметрией, чем базовые ОЦК подрешетки, путем объединения подрешеток 1-4, 2-5, 3-6 (а, б, в на рис. 3). Каждая ПТ подрешетка анизотропна, однако, в совокупности три ПТ подрешетки сохраняют кубическую симметрию составного кристалла. Результаты МД моделирования показали, что при определенных температурах ионы Ag+ могут конденсироваться по отдельным ПТ подре-шеткам, представленным на рис. 3. Это может быть причиной зарождения анизотропной фазы и создания предпосылок для а^в перехода.
,_Р 0^
а о в
Рис. 3. Анизотропные подрешетки в дробно занятых узлах соединений a—AgI, Ag2X (X = Б, Бе, Те)
Кристаллическая структура некоторых кубических фаз Ag2X (X = 8, Бе, Те) оказывается аналогичной структуре высокотемпературной фазы AgI -a-AgI, только в отличие от a-AgI в междоузлиях Ag2X располагаются не 2, а 4 иона Ag+. Это приводит к увеличению кулоновского отталкивания в катионной под-решетке и распределение катионов по междоузлиям несколько отличается от того, которое имеет место в а— AgI. Усиливается неравномерность заселения серебра в ё и И позициях; так, в Ag2S серебро занимает главным образом тетраэдрические позиции, причем с ростом температуры эта тенденция усиливается. Вследствие сильной корреляции в катионной подсистеме кубических модификаций халькогенидов серебра наблюдаются анизотропные флуктуации. Так, в Ag2S в Ag-подсистеме по данным рентгеновского и нейтронного диффузного рассеяния образуются динамические микродомены сигарообразной формы, диаметром ~ а, где а - параметр решетки, вытянутые вдоль направления <100>, со временем жизни ~ 5*10 11с.
Соединения Ag3SI и Ag3SBr называют обратными перовскитами (антиперовскитами), так как в них места анионов в перовскитах занимают катионы. Эти соединения являются одними из лучших ионных проводников. Структура a-Ag3SI высокотемпературной полиморфной модификации образует объемноцентрирован-ную элементарную ячейку благодаря тому, что анионы Б и I разупорядочены по позициям (0,0,0), (/г,1/,1/), ионы серебра распределены так же, как и в a-AgI по позициям 12ё(42ш). Для этих позиций ближайшими соседями являются не только позиции на данной грани, но и на соседней. Вот почему проводимость в a-AgI и а-Ag3SI много больше, чем в в-Ag3SI и в-Ag3SBr. Исследование в-Ag3SI показало, что анионы 8 и I сосредотачиваются по позициям (0,0,0) и (/,/,/). Шесть ионов серебра распределены по 12И позициям с точечной симметрией (шш2): х / 0, 0 х /, / 0 х, х 0 /, / х
0, 0 / х, х / 0, 0 х /, / 0 х, х 0 /, / х 0, 0 / х, где х = 0.390, отклонение от центра граней около 0.5 А. Подобную же структуру имеет и в-Ag3SBr с х = 0.385. Ag3SBr разлагается при температурах выше 703 К и фаза a-Ag3SBr не обнаружена. Размеры ионов серебра не позволяют расположить на грани более одного атома. Таким образом, в кубических фазах Ag3SI и Ag3SBr ионы серебра располагаются в позициях, близких к центрам граней, так что в итоге образуется структура, близкая к кубической перовскитной, что и оправдывает
название «антиперовскит». В этом случае мы возвращаемся к приближенным подрешеткам, характерным для классического перовскита (рис. 1а), только для СаТЮ3 эти подрешетки были точными.
Приведенные примеры показывают, что формальное введение анизотропных подрешеток в кубических кристаллах отражает поведение таких кристаллов в процессах упорядочения (разупорядочения), фазовых переходах. Так, в кубических модификациях кристаллов со структурой перовскита на основе анионных Гс подрешеток могут быть сформированы ПГ подрешет-ки, трансляционная симметрия которых выше трансляционной симметрии кристалла, однако, точечная симметрия ниже. При этом каждая подрешетка анизотропна, но симметрия анионной подрешетки в целом — кубическая. Такое представление структуры перовскита выявляет внутреннею (скрытую) симметрию анионной подрешетки, которая будет проявляться в ветвях спектров элементарных возбуждений, сформированных преимущественно из анионных состояний, если такие возбуждения реализуются. Со структурной точки зрения выделение анизотропных подрешеток связано с наличием потенциальной возможности в кубических кристаллах фазовых переходов в низкосимметричные структуры, которые для структуры перовскита представлены на рис. 16, в. Такие возможности также проявляются и в кристаллах замещения в анионной под-решетке, в частности, в системах СбРЬР2С1 и СбРЬРС12 рентгеноструктурный анализ показал структуру перов-скита со слабыми ромбическими искажениями. В кубических фазах соединений Ag3SI и Ag3SBr ионы серебра располагаются в позициях близких к центрам граней, так что в итоге образуется структура, близкая к кубической перовскитной. В этих соединениях можно ввести приближенные тетраэдрические подрешетки, аналогичные точным, представленным на рис. 1а.
В кубических модификациях серебросодержащих суперионных кристаллов с катионной проводимостью анизотропные подрешетки могут быть выделены на основе междоузельных позиций, по которым происходит миграция ионов серебра. Эти анизотропные подрешетки для a-AgI, Ag2X (X = S, Se, Те) представлены на рис. 3а, 6, в. Для этих кристаллов в ряде работ, основанных на различных версиях МД исследованы динамическая структура и термодинамические свойства суперионной подрешетки серебра. При высоких температурах среднеквадратичные смещения достаточно высоки и изотропны. С понижением температуры
среднеквадратичные смещения и диффузия оказываются анизотропными и при определенных моделях потенциала взаимодействия частицы Ag+ распределяются по некоторым подрешеткам, представленным на рис. 3а, б, в, в то время как другие оказываются пустыми. Появление анизотропии в кубической фазе и создает предпосылки для а^в перехода.
2. Симметрия кристаллов с подрешетками в многомерных кристаллических пространствах. Симметрия реальных физических объектов не всегда может быть адекватно описана на языке трехмерного кристаллического пространства. Достаточно простым физически прозрачным примером являются кристаллы с подрешетками. Кристаллографическая группа симметрии таких кристаллов содержит элементы симметрии, отвечающие подрешеткам с наименьшей симметрией. Полная симметрия таких кристаллов описывается только в пространстве, размерность которого равна суммарной размерности подрешеток, относящихся к разным системам Браве или имеющих различные периоды. Что касается размерности самих подрешеточ-ных кристаллических пространств, то при изучении симметрии кристаллов некоторых материалов (например, ферромагнитных) имеет значение не только геометрическое положение атомов, но и их состояния (например, электрический или магнитный заряд). Тогда необходимо рассматривать преобразования, которые могут менять не только геометрическое положение объекта, но и его состояние, то есть увеличивать размерность пространства, в котором они действуют. Это
( а11(Г(« * )
а(Г( „),) =
делается при нахождении кристаллических классов шубниковских групп. Метод, представленный ниже, наиболее важен для кристаллов, имеющих сложное строение и химический состав. Прежде всего, именно в сложных системах наиболее вероятно существование дополнительной симметрии.
Свойства многомерных кристаллических пространств с трансляционно совместимыми подпространствами исследовались в [1 - 3]. Мы будем рассматривать кристаллические пространства, размерность которых N кратна некоторому числу к (Ы = пк), и кристаллографические группы которых распадаются на подгруппы, элементы симметрии которых действуют только в подпространствах размерности п. Иными словами, мы рассматриваем многомерное кристаллическое пространство Я(пк), составленное из прямой суммы к подпространств размерности п:
Япк) = Ът ® Яп)2 ® ••• ® Яп)к• (1)
В каждом Я(ф (я = 1, 2,..., к) может реализоваться лишь конечное число кристаллографических групп С(„)., (^ = 1, 2,., к), каждая из которых содержит в себе подгруппу параллельных переносов Г^ (я = 1, 2,., к). В случае п = 3 неэквивалентные группы исчерпываются 230 федоровскими кристаллографическими группами; подгруппы параллельных переносов определяются 14 решетками Браве. Введем в подпространствах Я(ф базисы решетки из п линейно независимых векторов и,(Г(ф) (/ = 1, 2,., п), которым отвечают генерирующие матрицы:
I (Г(п)я )
а
а
(Г(п)я ) а22(Г(п)я )
(Г(п)я ) ап2 (Г(п)я )
а
а
Л
а
г (Г(п)5 )
г (Г(п)5 )
(2)
Примитивная ячейка решетки в Я^ представляет собой п-мерный параллелепипед, вершинами которого являются точки решетки, и имеет объем, равный абсолютной величине определителя генерирующей
матрицы |ёе1 а(Г(п)я)| . Минимальные ячейки полностью заполняют все пространство Я^, не пересекаясь друг с другом. Сопряженное или обратное пространство определяется системой векторов Ь/(Г(п)!!) (/ = 1,2,...,п), которые называются биортогональными векторам а,(Г(п)я) и удовлетворяют соотношению (^(Г^Ь/Г^)) = 8/. Генерирующая матрица сопряженного пространства и метрические тензоры прямого и сопряженного пространств даются соответствующими соотношениями:
*(Г(п)я) = (аТ(Г(п)я г1;
ёа = а(Г(п)я )аТ (Г(п)я ); (3)
ёЬ = Ь(Г(п)я Ь (Г(п)я X
где Т означает транспонирование, ёа и ё - матрицы Грамма прямого и сопряженного пространств.
Введем нумерацию пространств Я^ в соответствии с убыванием объемов их примитивных ячеек:
|det а(Г(п)1 ^ ^ |det а(Г(п)2 ^ ^ |det а(Г(п)к ^ , по-
добным же образом будем нумеровать группы О(п)я, подгруппы Г(п), и векторы а,(Г(п)х), А/Г^). Определим понятие трансляционной совместимости подпространств Я(п)я под которым мы будем понимать наложение на подпространства Я^ (я = 2, 3,., к) таких условий, которые позволяют вложить эти подпространства в подпространство Я(п)1, имеющее максимальный объем примитивной ячейки, с сохранением в таком составном пространстве трансляционной симметрии Г(п)1. Соотношения трансляционной совместимости между базисными векторами а^Г^) и «/(Г(п>) (я = 2, 3,., к) запишем:
п
а (Г(п)1) = Е (Г(п)1 I Г(п)я )/ а/ (Г(п)я X (4)
/=1
где (Г(п)^ Г(п)я )г/ - матрицы связи, обеспечивающие
трансляционную совместимость пространства Я(п)1 с пространствами Я(п)х (я = 2, 3,., к). Чтобы пространства Я(п)я были инвариантны относительно операций группы параллельных переносов пространства Я(п)1,
необходимо и достаточно, чтобы все элементы матрицы связи были целыми числами. Реально соотношение (4) налагает ограничения на компоненты векторов аДГ^) пространств R(n)s(s = 2, 3,..., к). Записывая формулу (4) в терминах генерирующих матриц (2)
а(Г(n)l) = (Г(n)11 Г(n)s )а(Г(n)s ), (5)
получим выражение для матриц связи:
(Г(n)11 Г(n)s ) = а(Г(n)1 )a(Г(n)s ). (6)
Решим далее задачу конструирования кристаллического пространства R(„l, отвечающего гипотетическому сложному кристаллическому соединению путем вложения в пространство R(n)1 пространств R(n)s(s = 2, 3,., к). Для этого необходимо задать расположение всех подрешеток s = 1, 2,., к друг относительно друга. Если за общее начало отчета выбрать узел в пространстве R(n)1, то начала отчета оставшихся подрешеток смещаются относительно общего начала на некоторые векторы с^Г^) (s = 2, 3,., к). Построенное таким образом кристаллическое пространство R(n)L оказывается составленным из подрешеток, каждой из которых отвечает свое кристаллическое пространство. Подгруппа трансляций R(n)L совпадает с Г(„)Ь как это следует из (4) - в сущности, из самого определения трансляционной совместимости. Что касается точечной симметрии пространства R(n)L, то она будет зависеть как от точечной симметрии всех пространств R(n)s (s = 1,2,.,n), так и от способа их размещения, т. е. от выбора векторов с,(Г(п)х). В итоге определяется и пространственная группа симметрии пространства R(n)L -
G(n)L.
В пространстве R(n)L имеется связь между элементами симметрии и векторами с,(Г(ф) - под действием элементов симметрии из группы G(n)L эти векторы должны переходить в эквивалентные положения:
Zg (у) s-(Г („) s)+т?\- = j=1
n (7)
= ^Г(,) s-) + £ n Г а, (Г(„ )l),
j=1
/V (у)
где g (у = 1, 2, ..., M) - точечный элемент группы
(у)
симметрии G(nL, состоящей из M элементов, - не-
целый параллельный перенос для несимморфных
(у) w
групп, n j - набор целых чисел, соответствующий каждому элементу симметрии, С^(Г( n)s,) совпадает
либо с с,(Г(п)х), либо с началом отсчета симметрично эквивалентной подрешетки.
Симметрия сложного трехмерного кристалла, составленного из двух подрешеток Браве, обладающих различной, но совместимой с кристаллографической трансляционной симметрией соединения, может быть определена в 6-мерном пространстве R6 [3]. В этом пространстве имеется 91 семейство, 841 класс Браве, 7104 геометрических, 85311 арифметических классов и 28 927 922 ПГ [4]. В [5] представлена информация о 91 семействе в R6, для которых предложена нумерация,
обозначения в системе Германна-Могена, приведены некоторые метрические тензоры.
Из 91 семейства в Я6 пригодными для описания симметрии кристаллов с двумя подрешетками являются те, которые представляют собой прямые произведения подгрупп, действующих в Я3 взаимно ортогональных инвариантных подпространствах. Таких семейств оказывается 21, необходимая информация о них представлена в таблице 1.
В этой таблице первый столбец - предлагаемый в [4; 5] номер семейства, №р - число свободных параметров в метрическом тензоре, |Н| - порядок максимальной голоэдрической группы в семействе, |В| - число классов Браве в семействе. Более детально обозначения разъясняются в [4; 5]; информация необходимая для наших целей, будет сообщаться по мере необходимости. В таблице 2 приведены метрические тензоры семейств, представленных в таблице 1. Номера в нижних левых углах тензоров в таблице 3 обозначают номера семейств. Часть таблицы 1 представляет собой фрагмент таблицы 8 из [5]. Нами добавлен столбец с числом классов Браве в семействах, кроме того, для семейств 13, 14, 15, 18, 29, 30, 40, 42, 53, 54, 57 введена другая нумерация векторов базиса, чем принятая в [5], для того, чтобы, не выходя за рамки семейства, перейти к прямым произведениям подгрупп, действующих в Я3. Например, для семейства 14, вместо символа 4т ^
13 ± т, принятого в [5], нужно ввести 4т ± т ± 1 3. Подобным же образом можно перенумеровать векторы базиса и переобозначить символы для других перечисленных семейств.
Как видно из таблицы 1, семейство 4 отвечает двум ортогональным, относящимся к триклинной сингонии подпространствам в Я3, образующим один класс Браве. Четыре элемента голоэдрической симметрии включают единичный элемент, два пространственных отражения в каждом отдельном Я3 и одновременное отражение в обоих Я3. Семейство 13 отвечает ортогональным подпространствам Я3, относящимся к моноклинной синго-нии. В этой сингонии имеется две решетки Браве -простая и базоцентрированная. Следовательно, в Я6 будет 3 класса Браве, два из которых относятся к одинаковым решеткам Браве в Я3 и один - к различающимся. Семейство 8 относится к двум ортогональным Я3, одно из которых отвечает триклинной, другое моноклинной решеткам, что реализует в Я6 два класса Браве. Подобным же образом реализуются классы Браве всех семейств, представленных в таблице 1.
Рассмотрим подробнее инвариантные подпространства Я3 кубической сингонии. В кубической сингонии имеется три типа Браве - кубическая простая Гс, кубическая гранецентрированная Г/ и кубическая объ-емноцентрированная Г/. Группа преобразований для
всех трех типов (0Ь - Шенфлис, т 3 т- Германн-Моген) содержит 48 элементов симметрии. В Я3 голоэдрическая группа кубической сингонии 3-непро-водима. Семейства, содержащие такую группу: 21, 44, 55, 71, 74, 80. Те семейства, в которые 3-неприводимая группа входит однократно - 21, 44, 55, 71, 74 и действует в Я3, имеют ортогональные подпространства, отвечающие триклинной, моноклинной, ромбической, тетрагональной и гексагональной сингониям соответ-
ственно. Семейство, в которое 3-неприводимая группа входит дважды - 80. В этом семействе группа т3 т реализуется в двух ортогональных подпространствах Я3, голоэдрическая точечная группа является прямым
произведением групп т 3 т, число элементов симметрии 2 304 = 482, число классов Браве равно шести.
Далее необходимо учесть то обстоятельство, что Я3 подпространства, формирующие пространство Я6, должны быть трансляционно совместимы. При установлении трансляционной совместимости, естественно, необходимо учитывать схему подчинения сингоний. Обратимся к таблице 1 с этих позиций. Семейства 54 и 74 отвечают подпространству Я3 гексагональной син-гонии и ортогональным подпространствам тетрагональной и кубической сингоний соответственно. Согласно схеме подчинения сингоний, Я3 подпространства этих сингоний трансляционно не совместимы, что
не позволяет рассматривать эти семейства как пригодные для описания сложных кристаллов с двумя подре-шетками.
Помещение подрешеток высших сингоний в триклинную или моноклинную решетки возможно только при сохранении трансляционной симметрии высокосимметричной подрешетки. Это относится к семействам 8, 12, 14, 15, 18, 21, 29, 30, 44. В этом случае понижение симметрии до триклинной или моноклинной происходит за счет размещения подрешеток в низко -симметричных точках примитивной ячейки, то есть за счет понижения точечной симметрии. В качестве примера можно привести FeS2 и NiS2, в одной из своих фаз кристаллизующихся с ПГ С1, полностью теряющей точечную симметрию, но сохраняющей трансляционную симметрию простой кубической решетки.
Таблица 1
Семейства в Я6, содержащие прямые произведения подгрупп, действующих в Я3 ортогональных инвариантных подпространствах
№ Nfp \H\ \B\ R-приводимость Символ
4 12 4 1 (1+1+1)+(1+1+1) 13 113
8 10 8 2 (1+1+1)+(1+1)+1 13 1 2 1 m
12 9 16 4 (1+1+1)+1+1+1 131mmm
13 8 16 3 (1+1)+1+(1+1)+1 2 1 m 1 2 1 m
14 8 32 2 2+1+(1+1+1) 4m 1 m 113
15 8 48 1 2+1+(1+1+1) 6m 1 m 1 13
18 7 32 8 (1+1)+1+1+1+1 2 1 m 1 mmm
21 7 96 3 3+(1+1+1) m3m 1 13
28 6 64 10 1+1+1+1+1+1 mmm 1 mmm
29 6 64 4 2+1+(1+1)+1 4m 1 m 1 2 1 m
30 6 96 2 2+1+(1+1)+1 6m 1 m 121 m
40 5 128 8 2+1+1+1+1 4m 1 m 1 mmm
42 5 192 4 2+1+1+1+1 6m 1 m 1 mmm
44 5 192 6 3+(1+1)+1 m3m 1 2 1 m
53 4 256 3 2+1+2+1 4m 1 m 1 4m 1 m
54 4 384 2 2+1+2+1 6m 1 m 1 4m 1 m
55 4 384 12 3+1+1+1 m3m 1 mmm
57 4 576 1 2+1+2+1 6m 1 m 1 6m 1 m
71 3 768 6 3+2+1 m3m 1 4m 1 m
74 3 1152 3 3+2+1 m3m 1 6m 1 m
80 2 2304 6 3+3 m3m 1 m3m
Таблица 2
___________Метрические тензоры семейств, приведенных в таблице 1_____________
а й Ь 0 0 0 а е Ь 0 0 0 а g Ь 0 0 0 а g 0 0 0 0 а 0 0 0 0 0
Ъ 1 0 0 0 Ъ 1 0 0 0 Ъ 1 0 0 0 Ъ 0 0 0 0 а 0 0 0 0
с 0 0 0 с 0 0 0 с 0 0 0 с 0 0 0 Ь 0 0 0
<1 j к (1 ) 0 <1 0 0 <1 ь 0 с г ь
4 е 1 8 е 0 12 е 0 13 е 0 14 <1 1
Г Т т Г 15 е
а & 0 0 0 0 а 0 0 0 0 0 а 0 0 0 0 0 а 0 0 0 0 0 а 0 0 0 0 0
Ъ 0 0 0 0 а 0 0 0 0 Ъ 0 0 0 0 а 0 0 0 0 а 0 0 0 0
с 0 0 0 а 0 0 0 с 0 0 0 Ъ 0 0 0 Ъ 0 0 0
<1 0 0 ь г ь (1 0 0 с ь 0 с 0 0
18 е 0 21 с 1 28 е 0 29 (1 0 40 (1 0
г (1 г 30 е 42 е
а 0 0 0 0 0 а 0 0 0 0 0 а 0 0 0 0 0 а 0 0 0 0 0 а 0 0 0 0 0
а 0 0 0 0 а 0 0 0 0 а 0 0 0 0 а 0 0 0 0 а 0 0 0 0
а 0 0 0 Ъ 0 0 0 а 0 0 0 а 0 0 0 а 0 0 0
ь 11 0 53 с 0 0 ь 0 0 ь 0 0 ъ 0 0
44 с 0 54 с 0 55 с 0 71 ь 0 80 ъ 0
(1 57 (1 (1 74 с ъ
Рассмотрим в качестве примера кристаллы с решеткой флюорита СаР2, имеющие Г/ решетку, ПГ симметрии № 225 О5. К этой же симметрии относится и структура антифлюорита, являющаяся обратной структуре флюорита в смысле взаимной замены анионов и катионов. Структура флюорита может быть сформирована из двух кубических подрешеток анионов и катионов, при этом катионы занимают позиции в узлах Г/ подрешетки с периодом а, атомы аниона -
а
в узлах Гс подрешетки с периодом а/2, сдвинутой относительно катионной подрешетки на четверть пространственной диагонали кубической элементарной ячейки.
На рис. 4 представлены МДВ (а) и ЗБ (б) кристалла и подрешеток структуры флюорита. В таблице 3 приведены генерирующая матрица и метрический тензор в Я6 для этой структуры.
* кг
б
Рис. 4. МДВ (а) и ЗБ (б) подрешетки кристалла со структурой флюорита
ФИЗИКА | Таблица 3
Генерирующая матрица и метрический тензор в Я6 для кристаллов с решеткой флюорита
Генерирующая матрица Метрический тензор
0 а/2 а/2 0 0 0 а2/2 а2/4 а2/4 0 0 0
а/2 0 а/2 0 0 0 а2/4 а2/2 а2/4 0 0 0
а/2 а/2 0 0 0 0 а2/4 а2/4 а2/2 0 0 0
0 0 0 а/2 0 0 0 0 0 а2/4 0 0
0 0 0 0 а/2 0 0 0 0 0 а2/4 0
0 0 0 0 0 а/2 0 0 0 0 0 а2/4
В Я6 структура флюорита относится к семейству 80, имеет (48)2 точечных элементов симметрии и
предложенный символ т 3 т ^ т 3 т (таблица 1).
Точечные элементы симметрии группы т 3 т действуют независимо в ортогональных Я3 подпространствах - подпространствах анионной и катионной под-решеток. Однако трансляционная симметрия этих подрешеток оказывается разной как в прямом, так и обратном пространствах, при этом трансляционная симметрия кристалла и катионной подрешетки одинакова. Объем ЗБ кристалла и катионной подрешетки в два раза меньше объема ЗБ анионной подрешетки, поэтому ЗБ последней необходимо перестроить в ЗБ кристалла - то есть провести разложение неприводимых звезд волновых векторов анионной подрешетки по неприводимым звездам волновых векторов кристалла. На этой основе можно выполнять качественный анализ спектров элементарных возбуждений кристаллов с решеткой флюорита и антифлюорита.
3. Исследование электронных и колебательных спектров на основе метода подрешеток. Метод исследования генезиса зонных спектров кристаллов, основанный на анализе симметрии подрешеток, изложен в [6]. Проанализируем строение энергетических зон кристалла с позиции зонных состояний подреше-ток. Кристаллический потенциал Уь (г) представим в форме:
»»=Хад+4»^(г). (8)
где У5 (г) - потенциалы подрешеток, АЖЬЗ (г)-
потенциал, отвечающий за гибридизацию подреше-точных состояний. Введем далее векторы обратной решетки кристалла и подрешеток, которые обозначим Ь,(ГА) и Ь/Г) (,, у = 1, 2, 3) соответственно. Разложения потенциалов в ряды Фурье запишутся:
V(г) = X V(Ь,.(Гь))ехр[,Ь,(Гь)• г],
ь, (ГI)
'V(г) = X V(Ь,.(Г))ехр[,Ь,(Гх)• г], (9)
Ь„(Гя)
А»„ (г) = X АГ(Ь,.(^))ехр[,'Ь,.(^) • г].
ь, (Г1)
Суммирование ведется по всем соответствующим векторам обратной решетки, У(Ь^ь)), К(Ь,(Г5)),
ДW(Ьi(ГL)) - форм-факторы потенциалов. Величины форм-факторов определяют расщепления энергетических зон, рассеяние волн, взаимодействующих с электронами кристаллической решетки. Действие формфакторов в обратном пространстве определяется границами ЗБ, которые будут разными для кристалла и подрешеток при наличии подрешеток с отличными от кристаллической решетками Браве.
Поскольку примитивная ячейка кристалла вмещает в себя примитивные ячейки всех подрешеток, объем ее ЗБ меньше либо равен объему ЗБ любой подре-шетки, при этом наибольшей оказывается ЗБ подре-шетки с наименьшим объемом элементарной ячейки. Из-за трансляционной совместимости всех подреше-ток с кристаллической решеткой их ЗБ могут быть перестроены в ЗБ кристалла. Перестройка векторов ЗБ подрешеток к в векторы ЗБ кристалла кь производится согласно соотношению:
к і=к х+Е п(гь і г ж(г ь ), (10)
г=1
где и,(Гь|Гх) - целые числа, Ь,(Гь) - векторы обратной решетки кристалла. При перестройке ЗБ подрешеток в ЗБ кристалла необходимо принять во внимание точечную симметрию, которая может оказаться выше для подрешетки по отношению к кристаллической. В таком случае необходимо провести разложение неприводимых представлений групп симметрии подре-шетки по кристаллическим. Это разложение начинается с установления соответствия между неприводимыми звездами представлений подрешеток и кристаллической решетки с использованием соотношения (19), а далее ведется разложение неприводимых представлений групп волновых векторов. В итоге устанавливается генезис зонного спектра кристалла из под-решеточных состояний.
Размещение энергетических зон подрешеток в ЗБ кристалла приводит к появлению вырождений, обусловленных двумя факторами:
а) трансляционные вырождения, возникающие при свертывании спектров подрешеток;
б) вырождения, обусловленные более высокой точечной симметрией подрешеток, если таковая имеется.
Эти вырождения снимаются при учете гибридизации подрешеточных состояний, величины расщеплений определяются форм-факторами ДW(Ьi(ГL)). В
случае слабой гибридизации подрешеточных состояний хорошим нулевым приближением к зонному спектру кристалла будет спектр «свернутых» подрешеточных состояний, расщепления в котором могут быть вычислены, например, по теории возмущений.
Сильная гибридизация каких-либо подрешеток означает наличие между ними связи типа ковалентной, что проявляется в образовании молекулярных комплексов либо ковалентно связанных каркасов. Здесь, в зависимости от ситуации, подрешетки могут выбираться не в форме простых решеток Браве, а решеток с базисом, либо могут быть введены некоторые «усредненные», виртуальные подрешетки, в которых уже учтена наиболее сильная, ковалентная составляющая химической связи.
В качестве примера приведем кристаллы с решеткой флюорита. Подрешеточная структура этого кристалла описана в предыдущем разделе, МДВ и ЗБ кристалла и подрешеток представлены на рис. 4. Разложение неприводимых звезд волновых векторов анионной подрешетки на неприводимые звезды представлений группы симметрии кристалла представлено в таблице 4. Из таблицы видно, что в симметричных точках и на линиях ЗБ кристалла возникают трансляционные вырождения, обусловленные свертыванием зонного спектра анионной подрешетки.
Элементарная ячейка флюорита СаР2 содержит два семивалентных аниона (Б) и один двухвалентный катион (Са), число валентных электронов в ЗБ равно 16 и энергетический спектр валентной зоны СаБ2 состоит из 8 ветвей. Валентный зонный спектр 2-кратно заряженной анионной подрешетки аниона в ЗБ ПК решетки содержит 4 ветви, что естественно, так как ЗБ ПК подрешетки в два раза больше ЗБ кристалла.
Зонный спектр флюорита, рассчитанный в рамках теории функционала плотности методом псевдопотенциала, представлен на рис. 5. На этом же рисунке изображены зонные спектры заряженной подрешетки фтора - исходный спектр в ЗБ анионной подрешетки и свернутый спектр в результате помещения спектра подрешетки фтора в ЗБ кристалла.
Как видно из рисунка, свернутый зонный спектр заряженной подрешетки фтора повторяет все особенности кристаллического зонного спектра в области валентной зоны. Двукратное вырождение в точках Ь и X, возникшее при свертывании, может быть снято с применением теории возмущений. Сложные топологические особенности верхней валентной зоны СаР2 можно воспринять как чисто случайные, возникшие за счет пересечения различных ветвей спектра. Однако проведенный анализ показывает, что эти особенности не случайные, а являются результатом свертывания спектра анионной подрешетки.
Уравнения динамики кристаллической решетки принципиально отличаются от зонного уравнения. Если электронная подсистема может рассматриваться как единая из-за принципа эквивалентности частиц (все электроны имеют одинаковые массы, заряд и спин) и отнесение электронов к какой-либо подре-шетке заведомо условно, то уравнения динамики решетки с самого начала записываются с учетом конкретных характеристик структурных элементов - масс и параметров сил связи. Фактически здесь возникает многомерное пространство степеней свободы, симметрию которого и нужно исследовать. Такое исследование детально проведено в [7].
Таблица 4
Разложение неприводимых звезд представлений группы симметрии анионной подрешетки по неприводимым звездам кристалла со структурой флюорита
Решетка кристалла Подрешетка
тип звезды векторы звезды (ед. 2л/а) тип звезды векторы звезды (ед. 2л/а)
Г (0, о, 0) ГЯ (0, 0, 0) (1, 1, 1)
X (0, 1, 0) X (0, 1, 0)
W (1/2, 1, 0) Ъ0 (1/2, 1, 0)
Ь (1/2, 1/2, 1/2) Л0 (1/2, 1/2, 1/2)
£о (1, 1, 0) м (1, 1, 0)
Д (0, 2ц, 0), -1/2<ц<1/2 ДТ’ (0, 2ц, 0), -1/2<ц<1/2 (1, 2ц, 1), -1/2<ц<1/2
ъ (2ц, 1, 0), -1/4<ц<1/4 ъ (1+2ц, 0, 1), -1/2<ц<-1/4 (2ц, 1, 0), -1/4<ц<1/4 (-1+2ц, 0, 1), 1/4<ц<1/2
Е (2ц, 2ц, 0), -3/8<ц<3/8 Е 8’ (1+2ц, 1+2 ц, 1), -1/2<ц<-3/8 (2ц, 2ц, 0), -3/8<ц<3/8 (-1+2ц, -1+2ц, 1), 3/8<ц<1/2 (1+2ц, 1+2 ц, 1), -1/2<ц<-1/8 (2ц, 2ц, 1), -1/8<ц<1/8 (-1+2ц, -1+2ц, 1), 1/8<ц<1/2
Л (ц, ц, ц), -1/2<ц<1/2 Л (1+2ц, 1+2ц, 1+2ц), -1/2<ц<-1/4 (2ц, 2ц, 2ц), -1/4<ц<1/4 (-1+2ц, -1+2ц, -1+2ц), 1/4<ц<1/2
Рис. 5. Зонные спектры кристалла СаЖ2 и заряженной подрешетки V1 (значения энергии даны в Эв).
На крайне правом рисунке изображен зонный спектр подрешетки V1 в ЗБ ПК-подрешетки фтора, на среднем рисунке изображен зонный спектр подрешетки V1 в ЗБ кристалла СаЖ2 в результате свертки
Рассмотрим сложную кристаллическую структуру с г частицами в элементарной ячейке, координаты которых обозначим Я(у) (V = 1, 2, ..., г). Предположим далее, что рассматриваемую кристаллическую структуру можно представить как совокупность к вставленных друг в друга подрешеток Браве. Если в элементарной ячейке кристалла все атомы различны, то число подрешеток к будет равно числу частиц в элементарной ячейке г, трансляционная симметрия под-решеток и кристалла будет одинаковой. Что касается точечной симметрии, то она может быть выше у под-решеток, чем у кристалла, если их начала отсчета помещены в низкосимметричные точки элементарной ячейки. Если же в элементарной ячейке имеются одинаковые частицы, расстояния между которыми находятся в рациональном отношении к параметрам кристаллической решетки, то появляется возможность введения подрешеток с объемами элементарных ячеек
меньшими, чем объем элементарной ячейки кристалла.
Начала отсчета для каждой подрешетки с® (£ = 1,2, ..., к) целесообразно совместить с векторами с(5)еК(у), при этом часть одинаковых атомов в элементарной ячейке окажутся связанными друг с другом векторами трансляций подрешеток. Здесь целесообразно ввести двойную индексацию частиц - наряду с номером частицы в элементарной ячейке V ввести номер подрешетки 5”:
Я(у5) = о(у5 ) + и(у5 ), (11)
где векторы п(у5) равны нулю для случая, когда число подрешеток к равно числу атомов в примитивной ячейке г, и отличны от нуля при к<г и тогда они связывают начала отсчета 5-й подрешетки с v-м атомом элементарной ячейки.
В соответствии с нумерацией частиц (20) динамическая матрица размерности 3гх Згможет быть записана в форме блоков:
В(ь)(к) =
_М 1" 2
г - 1 г
12
В,
В к
В1
В к
г - 1 г \ к
В1
В к
(12)
диагональные блоки В11, . , В55, . ,Вкк отвечают взаимодействию атомов в отдельных подрешетках, недиагональные В5к . - взаимодействию подрешеток друг с другом.
Динамическая матрица подчиняется известным соотношениям, вытекающим из однородности и изотропности пространства, симметрии кристаллической решетки, а также условием совместимости микроскопической теории и теории упругости. При наличии в
кристаллическом соединении высокосимметричных подрешеток динамическую матрицу (12) необходимо еще подчинить дополнительным требованиям, вытекающим из этой подрешеточной симметрии. Учет дополнительной симметрии динамической матрицы позволяет провести качественный анализ колебательных спектров без проведения вычислений и интерпретировать некоторые их особенности, если спектры вычислены. К таким особенностям относятся прежде всего квазивырождения, возникающие за счет более высокой симметрии подрешеток. Трансляционные квазивырождения возникают при условии, если объем примитивной ячейки подрешетки меньше объема примитивной ячейки кристалла, и такие вырождения находятся путем свертывания ветвей спектров, относящимся к подрешеткам в ЗБ кристалла. Квазивырождения могут быть также обусловлены более высокой, чем кристаллическая, точечной симметрией подреше-ток. Такие квазивырождения устанавливаются путем разложения неприводимых представлений кристаллографических групп подрешеток по неприводимым представлениям группы симметрии кристалла. Что касается величины расщепления в области квазивырождений, то она будет тем меньше, чем больше разница атомных масс подрешеток и чем меньше силовые постоянные взаимодействия подрешеток.
Проиллюстрируем изложенное на примере кристалла с решеткой флюорита. Динамическая матрица (12) в случае структуры флюорита (г = 3) имеет размерность 9x9 и может быть записана в виде блоков трехмерных матриц, каждая из которых отвечает определенному типу атомов или взаимодействию между ними:
Аі А2 | Аз
Рл. _р22 _[ р23_
= Р31 Р32 ! Р33
(13)
где блок Р11 отвечает атомам металла, двумерный блок - атомам фтора; блоки Р12, Р13, Р21, Р31 - взаимодействию металла и фтора. Двумерный блок пред-
ставляет собой динамическую матрицу простой кубической решетки, составленную из атомов аниона, с удвоенной элементарной ячейкой. При этом Б22 = Б33, а в точке Ь ЗБ кристалла Б23 =032 = 0. Таким образом, в этой точке для подрешетки фтора имеет место вырождение фононных ветвей, обусловленное «свертыванием» ветви из ЗБ кислорода в ЗБ кристалла. Далее в расчете использовалась модель «замороженных» подрешеток. В этой модели Б12 = Б13 = Б21 = = Б31 = 0, но это не означает, что исключается взаимодействие металла и аниона, взаимодействие между ними входит через силовые матрицы в диагональных элементах динамической матрицы. Фактически при «замораживании» подрешеток взаимодействие металл-анион учитывается в первом порядке теории возмущений, второй порядок учитывается при «размораживании» - учете недиагональных элементов Б12, Б13, Б21, Б31. Дальнейшие вычисления фононных спектров кристалла и подрешеток выполнялись в изложенной модели.
На рис. 6 представлены результаты вычислений колебательных спектров кристалла РЪБ2 и его «замороженных» подрешеток в модели жестких ионов [8]. В центральной части рисунка приведен фононный спектр кристалла, светлыми и темными кружками обозначены экспериментальные данные по нейтронному рассеянию. В правой части рисунка представлены фононные спектры подрешетки фтора (точечные линии) в «замороженной» подрешетке металла и под-решеток металла в «замороженной» подрешетке фтора (сплошные линии). В фононном спектре подрешетки фтора в точке Ь ЗБ имеет место вырождение частот, обусловленное «свертыванием» ветвей из ЗБ фтора в ЗБ кристалла. «Размораживание» подрешеток и взаимодействие их колебаний приводит к снятию вырождения, однако, как это видно из центральной части рисунка, соответствующие расщепления невелики, что обусловлено существенным различием масс анионов и катионов. В фононных спектрах подреше-ток отсутствуют акустические ветви, поскольку условия их реализации выполняются для полной динамической матрицы (21), а не ее отдельных блоков.
Рис. 6. Фононные спектры, плотности фононных частот кристалла РЪЖ2 и его подрешеток
188 I Вестник Кемеровского государственного университета. 2014 № 2 (58) Т. 1
Сравнение фононных спектров кристаллов (центральная часть рисунка) со спектрами подрешеток (правая часть рисунка) показывает, что оптические ветви кристаллов чрезвычайно близки к спектрам под-решетки фтора в «замороженной» подрешетке металла.
В левой части рисунка приведены плотности фонон-ных частот кристалла и подрешеток металла и фтора.
Кристаллическая плотность представлена жирной линией, подрешетки металла - тонкой, а фтора - точечной. Из рисунка видно, что в высокочастотной области плотность фононных частот кристалла оказывается достаточно близкой к плотности частот подрешетки фтора, а в низкочастотной основные особенности кристаллической плотности частот коррелируют с особенностями плотности частот подрешетки металла.
Полученные результаты показывают, что метод подрешеток позволяет установить некоторые интересные особенности частотных спектров сложных кристаллических соединений. В случае соединений с существенно различающимися массами частиц в подре-шетках (на примере PbF2) показано, что полный фо-нонный спектр кристалла достаточно хорошо количественно (оптические ветви) или качественно (акустические ветви) описывается спектрами в поле «замороженных» подрешеток. Это позволяет оценивать различные статистические и термодинамические характеристики кристалла в названной простой модели.
4. Заключение. Развиваемый нами в последние годы метод структурных подрешеток показал себя эффективным при решении многих задач кристаллографии, кристаллохимии и физики конденсированного состояния.
1. В области кристаллографии представление сложных кристаллических структур совокупностью подрешеток Браве позволяет исследовать внутреннюю «скрытую» симметрию кристаллического соединения, обусловленную возможной более высокой симметрией подрешеток по сравнению с федоровской группой
Литература
1. Поплавной, А. С. Многомерные кристаллические пространства с трансляционно совместимыми подпространствами I А. С. Поплавной II Известия вузов. Физика. - 200б. - Т. 49. - № б. - С. 10 - 14.
2. Поплавной, А. С. Кристаллы с подрешетками и их симметрия в многомерных кристаллических пространствах
I А. С. Поплавной II Кристаллография. - 200б. - Т. 49. - № 10. - С. 277 - 2В0.
3. Поплавной, А. С. Симметрия кристаллов, составленных из двух подрешеток Браве в шестимерных кристаллических пространствах I А. С. Поплавной II Кристаллография. - 2010. - Т. ЗЗ. - № 2. - С. 1В1 - 1Вб.
4. Report of a Subcommittee on the Nomenclature of n-Dimensional Crystallography. II. Symbols for geometric crystal classes, Bravais classes and space groups I T. Janssen [et al.] II ActaCryst. - 2002. - V. AЗВ. - P. б0З - б21.
З. Report of a Subcommittee on the Nomenclature of n-Dimensional Crystallography. I. Symbols for point-group transformations, families, systems and geometric crystal classes I T. Janssen [et al.] II ActaCryst. - 1999. - V. A55. - P. 7б1 - 7В2.
6. Поплавной, А. С. Симметрия подрешеток и генезис спектров элементарных возбуждений в кристаллах I А. С. Поплавной II Материаловедение. - 200З. - № 9. - С. 2 - 7.
7. Поплавной, А. С. Симметрия фононных спектров в кристаллах с подрешетками I А. С. Поплавной II Известия вузов. Физика. - 200В. - Т. З1. - № 7. - С. 31 - ЗВ.
В. Теория физических и физико-химических свойств кристаллических соединений с различным типом химической связи I под общ. ред. А. С. Поплавного. - Кемерово: КемГУ, 2012. - 400 с.
Информация об авторе:
Поплавной Анатолий Степанович - доктор физико-математических наук, профессор, заведующий кафедрой теоретической физики КемГУ, В(3В42)ЗВ-31-9З, popl@kemsu.ru.
Anatoly S. Poplavnoy - Doctor of Physics and Mathematics, Professor, Head of the Department of Theoretical Physics, Kemerovo State University.
Статья поступила в редколлегию 26.02.2014 г.
сложного кристаллического соединения. Несмотря на то, что эта симметрия за счет взаимодействия подре-шеток не является точной (псевдосимметрия), она проявляется в физических и физико-химических свойствах кристалла. Развитием кристаллографического подхода к исследованию свойств твердых тел является применение многомерных кристаллических пространств для описания структуры соединения, построенного из под-решеток Браве, использование классификации Делоне кристаллических решеток, на основе которой сформулирован метод построения трансляционно совместимых МДВ, найден граф подчинения сортов Делоне по принципу повышения симметрии.
2. В области кристаллохимии разработан новый подход к анализу химической связи на основе исследования вкладов электронной плотности подрешеток в полную электронную плотность кристалла и введения понятия разностной электронной плотности, учитывающий эффект гибридизации подрешеток. С применением этого подхода установлены новые особенности химической связи сложных кристаллических соединений с преимущественно ионной, ионно-молекулярной и ионно-ковалентной типами химической связи.
3. В области физики конденсированного состояния развит метод исследования генезиса электронных и колебательных спектров сложных кристаллических соединений из состояний их подрешеток. Установлены особенности электронных и колебательных спектров, обусловленные более высокой симметрией подреше-ток: вырождения и квазивырождения спектров, возникающие за счет свертывания ветвей спектров из ЗБ подрешеток в ЗБ кристалла, а также вызванные более высокой точечной симметрией подрешеток. Новые возможности для анализа фазовых переходов в сложных кристаллических соединениях представляет найденный граф подчинения сортов Делоне по принципу повышения симметрии.
