Особенности структуры и физических свойств некоторых суперионных кристаллов

УДК 548. 713: 548. 1. 025

ОСОБЕННОСТИ СТРУКТУРЫ И ФИЗИЧЕСКИХ СВОЙСТВ НЕКОТОРЫХ СУПЕРИОННЫХ КРИСТАЛЛОВ

Д. С. Нефёдова, Е. В. Николаева, А. С. Поплавной, Т. П. Федорова

STRUCTURE AND PHYSICAL PROPERTIES OF SOME SUPERIONIC CRYSTALS

D. S. Nefyodova, E. V. Nikolaeva, A. S. Poplavnoy, T. P. Fedorova

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

The paper discusses the mechanisms of superionic conductivity caused by features of the crystal structure, electronic and vibrational spectra of some series of crystals with anion and cation conductivity. Special attention is paid to the symmetry of the sublattices, which differs from the crystal structure in the implementation of superionic transfer.

Ключевые слова: суперионная проводимость, подрешетки, флюорит, перовскит, фононные спектры.

Keywords: superionic conductivity, sublattice, fluorite

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

Физический механизм, обеспечивающий супер-ионную проводимость, связан с высокой подвижностью какого-либо сорта ионов в решетке, образованной другим сортом (или сортами) ионов. Для реализации такого рода высокой подвижности необходимо выполнение одновременно нескольких условий. Во-первых, жесткая структура должна содержать заметно большее число вакантных позиций, чем ионов, которые могут их занимать. Во-вторых, позиции должны быть такими, чтобы энергия активации Еа между ними была не слишком велика, или точнее, чтобы отношение Еа/КвТ было порядка единицы при температурах Тс ниже температуры плавления (здесь КВ - постоянная Больцмана). В-третьих, должна существовать связная сетка траекторий движения ионов в жестком каркасе (каналов проводимости), пронизывающая всю структуру. В противном случае протекание постоянного тока оказывается невозможным. При определенном характере расположения этих каналов может реализоваться ионная проводимость пониженной размерности - одномерная, двухмерная.

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

, perovskite, phonon spectra.

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

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

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

1. Суперионные кристаллы с решеткой флюорита. Кристаллы MeF2 (Me = Ca, Cd, Sr, Pb, Ba) и MeO2 (Me = U, Pu, Th, Np) со структурой флюорита при высоких температурах обладают суперионной проводимостью [2, 3]. Они находят широкое практическое применение как в обычном, так и в суперион-

ном состояниях. Для точной оптики в области ультрафиолета используются кристаллы кубической структуры с коротковолновой границей пропускания, например, СаР2, ВаР2, 8гБ2. На основе С№2 и РЪР2 изготавливают химические источники тока, газовые датчики, твердотельные электрохимические датчики и сенсоры, различные оптоионные приборы. Кристаллы Ме02 применяются в качестве тепловыделяющих элементов (ТВЭЛ) ядерных реакторов.

Кристаллическая решетка флюорита - гранецен-трированная кубическая ГС, пространственная группа 01 (рис. 1а). Подрешетка металла по своему типу и пространственным периодам совпадает с кристал-

лической решеткой, а две Г^ подрешетки анионов можно рассматривать как одну Гс подрешетку с пространственными периодами вдоль всех осей, вдвое меньшими кристаллических. На рис. Іб представлены многогранники Дирихле-Вороного (МДВ) кристалла флюорита и подрешетки металлов (на рисунке обозначен цифрой І) и подрешетки анионов (2). МДВ кристаллической решетки Гf представляет собой ромбический додекаэдр (І), а МДВ Гс подрешетки аниона имеет вид куба (2) с объемом в 2 раза меньшим объема МДВ решетки.

Рис. 1. а) кристаллическая решетка флюорита; б) МДВ кристалла и подрешетки металлов

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

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

разования, миграции и взаимодействия дефектов относительно малы для суперионных проводников (десятые эВ). Вследствие тепловых возбуждений часть междоузлий разупорядоченным образом занята катионами (Френкелевские дефекты - ФД) или анионами (Антифренкелевские дефекты - АФД). Существование междоузельных ионов фтора в решетке этих соединений доказано опытами по дифракции нейтронов. При низких температурах концентрация дефектов мала, но имеются как экспериментальные, так и теоретические указания на то, что в суперион-ной фазе степень заполнения междоузлий может доходить до 0.3. Расчеты и эксперименты показывают, что энергия образования ФД в несколько раз выше энергий образования АФД и на этой основе сложились устойчивые представления об ионном транспорте в данном классе соединений. При этом различают два возможных механизма ионной проводимости: движение вакансий, образующихся при возникновении АФД, и движение междоузельных ионов. Как результаты обработки экспериментальных данных, так и молекулярно-динамическое моделирование показывают, что доминирует первый механизм ионной проводимости, для которого характерны меньшие значения энергии миграции [2; 3].

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

именно благодаря его наличию в теории объясняются быстрое возрастание числа дефектов при достижении температуры начала перехода в суперионное состояние и появление ряда аномалий в физико-химических свойствах фторидных суперпроводников. Существует общая тенденция плавного понижения энергии образования АФД от СаР2 кВаР2 и далее резкого уменьшения к РЪБ2, а также в зависимости от катиона в широком интервале температур наблюдается повышение проводимости от СаР2 к РЪБ2.

Влияние фононной подсистемы на переход кристаллов в суперионное состояние обсуждалось в ряде обзоров на качественном уровне, в том числе делались оценки параметров ангармонизма для рядов кристаллов с решеткой флюорита. Предлагаемый нами подход к описанию колебательных спектров МеБ2 и Ме02 с решеткой флюорита путем выделения из полного фононного спектра кристалла колебательных состояний подрешеток позволяет уже в гармоническом приближении исследовать различия в динамическом поведении подрешеток. В частности температурное поведение среднеквадратичных смещений атомов подрешеток показывает, почему именно анионная подрешетка будет разупорядочиваться с ростом температуры [3].

Особую роль в суперионной проводимости играют низкоэнергетические оптические фононы (НОФ). Вопрос о связи суперионной проводимости с особенностями НОФ исследовался в ряде работ [3]. Само существование таких фононов обусловлено некоторыми особенностями межионного взаимодействия. НОФ имеют большие значения коэффициентов затухания вблизи температуры суперионного перехода. В некоторых суперионных соединениях НОФ об-

ладают акустоподобным поведением. Особенности фононных спектров в значительной мере определяются структурой кристалла, которая также чрезвычайно важна и для суперионной проводимости. Проиллюстрируем особенности колебательных спектров кристаллов с решеткой флюорита на примере ТЮ2 [4].

Вычисления фононного спектраТЮ2 и его под-решеток, а также соответствующих плотностей частот выполнялось в модели жестких ионов с использованием потенциала Борна-Майера-Хаггинса в симметричных точках и направлениях ЗБ кристалла (1) (рис. 1е). На рис. 2 в центральной части приведены фононные спектры кристалла ТЮ2: светлыми кружками показаны экспериментальные данные, полученные методом нейтронного рассеяния, сплошными линиями результаты теоретических расчетов. В правой части рисунка представлены фононные спектры под-решетки кислорода (точечные линии) в «замороженной» подрешетке металла и подрешетки металла в «замороженной» подрешетке кислорода (сплошные линии). «Замораживание» означает обращение в ноль смещений соответствующей подрешетки. В фононном спектре подрешетки кислорода в точке Ь ЗБ имеет место вырождение частот, обусловленное «свертыванием» ветвей из ЗБ кислорода в ЗБ кристалла (рис. 1е). «Размораживание» подрешеток и взаимодействие их колебаний приводит к снятию вырождения, однако, как это видно из центральной части рисунка, соответствующие расщепления невелики, несмотря на большое различие масс тория и кислорода. В фонон-ных спектрах подрешеток отсутствуют акустические ветви, поскольку условия их реализации выполняются для полной динамической матрицы, а не ее отдельных блоков.

500

200

100

50 40 Цт), отн ед. 10 0

Рис. 2. Фононные спектры, плотности фононных частот кристалла ТН02 и его подрешеток

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

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

Обращает на себя внимание сильная дисперсия в оптической части фононного спектра кристалла и в спектре подрешетки кислорода в «замороженной» подрешетке металла. Отсюда широкая разрешенная полоса оптических колебаний, что не характерно для кристаллов с сильно различающимися массами анионов и катионов (например, в Тк02: шТ^ = 232 а.е., ш0 = 16 а.е.). В таких случаях обычно оптические и акустические фононы разделены широкой запрещенной полосой и спектры оптических фононов слабо дисперсны. Причина сильной дисперсии оптических ветвей может быть объяснена большим кулоновским взаимодействием между ионами кислорода в соответствующей подрешетке. Уширение полосы оптических фононов приводит к тому, что численные значения НОФ существенно отличаются от значений частот оптических фононов в верхней части спектра. Кроме того, значения НОФ по рядам соединений коррелируют со значениями температур перехода в суперионное состояние: Ме02 (Ме = и, Ри, ТИ, Мр) значения НОФ лежат в интервале частот 260 ^ 320 см-1 и Тс в интервале 2200 ^ 3000 К; МеБ2 (Ме = Са, 8г, С4 Ва) значения НОФ в интервале 160 + 190 см.'1 и Тс в интервале 1000 ^ 1400 К. Из плавных рядовых зависимостей выпадает РЪБ2 для которого резко уменьшается значение частоты НОФ до 96 см -1, так же как и Тс = 711 К.

Из изложенного выше ясно, что колебания подрешеток аниона (фтора и кислорода) в «замороженной» подрешетке металла достаточно хорошо описывают высокочастотную часть фононного спектра кристалла в целом. Это позволяет описывать статистические и термодинамические характеристики кристалл-лов МеБ2 и Ме02, относящиеся к анионам, рассматривая только эту подрешетку как простую решетку Браве, что существенно упрощает вычисления. Для подрешетки металла, хотя фононный спектр существенно отличается от кристаллического, однако пики плотности состояний коррелируют с кристаллическими, и по этой причине можно применять формулы для простой решетки Браве не только для вычисления термодинамических характеристик (например, факторов Дебая-Валлера) анионов, но и металлов, понимая, что здесь приближение будет более грубым, однако пригодным для качественных оценок [4].

О 400 800 1200 1600

Т, К

Рис. 3. Факторы Дебая-Валлера для анионов Б0 и катионов Вп в кристаллах ТН02

На рис. 3 представлены вычисленные и экспериментальные зависимости от температуры значений факторов Дебая-Валлерадля тория и кислорода при температурах ниже Тс . Результаты расчетов отображены сплошной и точечной линиями для кислорода и тория соответственно, экспериментальные данные представлены темными и светлыми кружками. Видно, что теория достаточно точно описывает эксперимент. Большие значения среднеквадратичных смещений аниона по сравнению с металлом и более быстрый рост с повышением температуры и приводит, в конце концов, к разупорядочению подрешетки анионов.

2. Суперионные кристаллы с решеткой перов-скита

К кристаллам с химической формулой АВХ3 относят соединения с разной кристаллической структурой, в том числе и переменного состава. Благодаря своим замечательным электрическим и магнитным свойствам особый интерес представляют составы со структурой перовскита. К этому структурному типу относятся многие сегнетоэлектрики, которые используются в самой быстро развивающейся отрасли -электронике. Некоторые кристаллы со структурой перовскита обладают достаточно высокой суперионной проводимостью [1].

Из всего многообразия соединений с перовскито-подобной структурой наиболее благоприятные условия для анионной подвижности имеются во фторидах. Впервые диффузионная подвижность ионов фтора была обнаружена среди фтористых перовскитов МеРЬБз (Ме = ЯЪ, Сб, Те). Наиболее полные исследования проведены на СбРЪР3. Прецизионные структурные исследования, проведенные на ряде кристаллов фтористых перовскитов, показали, что ионы фтора имеют аномально высокие амплитуды тепловых колебаний. В то же время эти исследования не позволили надежно разделить модель, связанную с появлением мягких мод за счет ангармонизма тепловых колебаний и модель порядок - беспорядок. Однако ряд косвенных данных (малое изменение энтропии при структурных превращениях, экспериментальное обнаружение мягких мод решеточных колебаний) свидетельствуют в пользу модели ангармонизма и полиморфных превращений типа смещения. Таким образом, модель разупорядочения, связанная с образованием АФД и характерная для структуры флюорита, в кристаллах перовскита, скорее всего, не работает.

Структурный тип перовскита АВХ3 (родоначальник СаТЮ3), где А, В - катионы, X- анионы хорошо известен. Элементарная ячейка может быть кубической, с пространственной группой 0\ (Рт3т), № 221, 7 = 1, В - (1а), А - (1Ъ) ‘/2 ‘/2 ‘/2, X - (3ф ‘/2 0 0, 0 / 0, 0 0 / (рис. 4а).

Рис. 4. Кубическая модификация перовскита: а) кристаллическая решетка; б) МДВ кристалла и подрешеток катионов

В кубическом перовските кристаллическая решетка относится к простому кубическому типу Браве (Гс). Структура перовскита составлена из 5 Гс подрешеток: 1А, 1В и 3Х. 3Х подрешетки расположены таким образом, что любые две из них можно представить одной базоцентрированной ромбической Г'о (/' = а, Ъ, с) подрешеткой, обладающей более высокой трансляционной симметрией, чем кубическая подрешетка кристалла. МДВ кубической модификации АВХ3 (1) и подрешетки Г X (2) представлены на рис. 46. В зависимости от того, какие две подрешетки X объединяются, центрирована будет грань, перпендикулярная к той или иной оси. Векторы элементарных трансляций для соответствующих типов Браве имеют вид:

Гс : а1 = (а, 0, 0), а2 = (0, а, 0), а3 = (0, 0, а);

Га: а = (а, 0, 0), а2 = (0, а/2, - а/2), а3 = (0, а/2, а/2);

Гьо : а = (а/2, 0, а/2), а^ = (0, а, 0), а3 = (-а / 2, 0, а/2);

Г : а = (а /2, - а/2, 0), а^ = (а /2, а/2, 0), а3 = (0, 0, а).

Подобное выделение подрешеток X приводит к кажущейся анизотропии, которая в реальности, конечно, отсутствует, если учесть полную структуру кристалла. Однако это выделение показывает, что в кубической модификации ABX3 имеются «потенциальная анизотропия» и «потенциальная неэквивалентность» атомов X.

В ряде работ высказывалась гипотеза, объясняющая возникновение ионной проводимости в СбРЪР3 резким увеличением числа дефектов. Принципиальная возможность обнаружения изменения концентрации дефектов при изменении объема образца обоснована в [1]. Появление подвижности ионов фтора при 220 К наблюдалось на основе метода ядерного магнитного резонанса (ЯМР). Измерения дифракции нейтронов в СбРЪР3 показали ромбоэдрические искажения ниже 187 К с изменением объема 0.113 %. По-

пытки экспериментальной проверки гипотезы о скачкообразном изменении концентрации точечных дефектов в кристаллах при возникновении ионной проводимости были предприняты путем измерения зависимости коэффициента теплового расширения СбРЪР3 от температуры в интервале 131 - 311 К. На этой зависимости обнаружена аномалия при 174.42 К, которая соответствует аномалии в теплоемкости, обусловленной полиморфным превращением при этой температуре. При 220 К аномалии теплового расширения не наблюдалось, хотя и был наибольший разброс экспериментальных точек. Таким образом, проводимость в СбРЪР3, по-видимому, не обусловлена механизмом умножения дефектов.

Для соединений ABX3 типа перовскита предполагается, что ионы X в высокотемпературной кубической фазе движутся внутри своего координационного полиэдра, а при охлаждении кристалла ниже некоторой критической температуры Тс ионы X «замерзают» в некоторой эффективной позиции. Такое изменение координат ионов X и соответственно симметрии решетки при охлаждении ниже Тс позволяет моделировать структурные превращения типа порядок-беспорядок с большим изменением энтропии, до Я 1п4. Для катионов в кубической модификации такая локальная подвижность не предполагается, так как для каждого из этих ионов имеется только по одной максимальной позиции. В соединениях КМпР3, TlCdF3, RЪCdР3, ЯЪСаР3 наблюдались превращения 0\ ^ , сопровождающиеся искажениями кристал-

лической решетки. Такая пространственная группа может реализовываться в случае разворота октаэдров ВР6 таким образом, что октаэдры повернуты в противоположные стороны в соседних слоях, ортогональных оси вращения, при этом оси вращения связаны с осями ромбических подрешеток Г'о (/' = а, Ъ, с). Измеренные углы поворотов элементарных октаэдров [ВХ6]4- являются структурными параметрами соединений. Причиной возникновения множественных позиций для ионов X является несоответствие межатомных расстояний А^, В^, соответствующих суммам

ионных радиусов. Параметры неизвестных кубических ячеек могут быть спрогнозированы на основании их зависимостей от размеров образующих атомов. Эти зависимости для фторидов приведены в [1]. Наличие таких зависимостей позволяет предложить для синтеза малоизученных веществ гипотетические структуры, полиморфные превращения в которых связаны с разворотами октаэдров. Действительно, такие соединения (ЯЪРЪР3, Сб8гР3, СбРЪР3) удалось обнаружить и исследовать. Если считать катионы А и В зафиксированными в позициях 1а и 1Ь, ионы Х в высокотемпературной фазе в результате тепловых колебаний перескакивают из одной позиции в другую; при понижении температуры ионы Х должны фиксироваться в одной из позиций. Изменение координат анионов Х приводит к искажению полиэдров, окружающих атомы А и В, что должно привести к изменению координат и этих атомов. При коррелированном распределении ионов для разных элементарных ячеек может оказаться, что среднее положение ионов Х изменилось. В этом случае можно говорить о происшедшем структурном полиморфном превращении типа порядок-беспорядок со смещениями Х по направлениям, соединяющими атомы В, с сильной анизотропией тепловых колебаний и с изменением энтропии.

Почти всегда размеры промежутков между атомами в кристаллической решетке меньше размеров подвижных ионов. Следовательно, для анализа возможности перемещения атомов или ионов по разным позициям необходимо учитывать способности ионов к деформации. Для более полного описания кристаллической решетки кроме ионного радиуса необходимо учитывать и физические свойства ионов. Каждый из ионов является сложной системой. Электронные оболочки могут смещаться относительно ядер. Сферическая форма ионов при этом искажается. Такое искажение в первом приближении можно учесть с помощью электронной поляризуемости. Поляризуемость количественно описывает деформацию, приводящую к образованию дипольного момента. Для нахождения электронной поляризуемости атомов используют различные феноменологические и перво-принципные методы [1].

Для количественного анализа способности разных ионов к подвижности необходимо учитывать размер промежутков в решетке и способность ионов к деформации, что выражается безразмерным параметром пропускания. Параметр пропускания численно равен отношению радиуса деформированной промежуточной позиции, образованной ионами одного заряда, к радиусу деформированного иона противоположного заряда. Для ряда соединений ABX3 параметры пропускания вычислены в [1] и сопоставлены с температурами фазовых переходов.

Вычисленные параметры пропускания для соединений CsPbX3, где X = Р, С1, Вг, I оказались соответственно равными 1.07, 0.75, 0.67 и 0.59. Рассчитанное тем же способом значение параметра пропускания для ионов серебра в известном ионном проводнике a-AgI составляет 0.57. Значения параметра пропускания для катионов оказываются значительно меньшими, чем для ионов фтора, но наблюдается тенденция роста их

параметров пропускания при увеличении размеров анионов Х и уменьшении радиусов катионов.

Впервые диффузионная подвижность ионов фтора была обнаружена среди фтористых перовскитов при исследовании методом ЯМР 19Р соединений АРЪР3 (А = ЯЪ, Сб, Те) [1]. Наиболее полные исследования проведены на СбРЪР3. Спектры ЯМР 19Р низкотемпературной фазы поликристаллического образца представляются одиночной широкой линией с заметной асимметрией. Возможная причина асимметрии спектра - структурная неэквивалентность фторов. Одновременное существование узкой и широкой компонент соответствует двум типам фторов, которые различаются динамическими свойствами. Различие химических сдвигов узкой и широкой компонент указывает на структурную неэквивалентность подвижных и малоподвижных фторов. Наблюдающееся поведение спектров можно интерпретировать следующим образом: в тетрагональной модификации выше 145 К начинает двигаться только часть ионов фтора. Кубическая модификация оказывается более благоприятной для движения, все фторы в ней диффундируют, и широкая компонента исчезает скачком. Это согласуется с выводами, полученными при исследовании электропроводности в СбРЪС13, где переход в кубическую модификацию сопровождается скачкообразным увеличением электропроводности и уменьшением энергии активации.

Энергия активации ионов фтора в СбРЪР3, найденная из релаксационных измерений, равна (0.23 ± 0.02) эВ для тетрагональной модификации и (0.11 ± 0.02) эВ для кубической. Имеется ряд факторов, указывающих на то, что в кубической фазе движение носит суперионный характер. Движение ионов фтора будет большей частью, по-видимому, по междоузлиям, также наличествует и вакансионный механизм. На примере СбРЪР3 и ЯЪРЪР3 сделано заключение, что кубические структуры оказываются более предпочтительными к диффузии фтора, чем тетрагональные.

Для выяснения влияния замещения ионов фтора на хлор проведено исследование систем СбРЪР2С1 и СбРЪРС12 [1]. Рентгеноструктурный анализ показал для этих соединений структуру типа перовскита со слабыми ромбическими искажениями. Эти искажения как раз имеют тип симметрии подрешеток фтора, представленный на рис. 46.

Первопринципные расчеты фононных спектров 10 кислородосодержащих соединений АТЮ3 (А = Са, 8г, Ва, Яа, Cd, 2и, Mg, ве, 8п, РЪ) выполнены в [5]. Качественно структура этих спектров подобна спектрам кристаллов с решеткой флюорита (рис. 2). Оптическая ветвь, отвечающая в основном колебаниям ионов кислорода, оказывается очень широкой, сливающейся с акустической. Это обусловлено, в частности, свертыванием фононных спектров из ЗБ подрешеток в ЗБ кристалла (рис. 4е). Отсюда можно сделать выводы о наличии НОФ в спектре этих кристаллов, аналогичные выводам для флюорита.

3. Серебросодержащие суперионные кристаллы

Серебросодержащие суперионные кристаллы (a-AgI, Ag2X (X = 8, 8е, Те), Ag3SI, Ag3SBr) обладают

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

Исследованию твердого электролита AgI посвящено много работ. Структура высокотемпературного йодистого серебра a-AgI от 420 К до температуры плавления 828 К описывается пространственной группой О[ - ]т3т. Ниже полиморфного превращения первого рода в-AgI принадлежит к структурному типу вюрцита с пространственной группой С64„ - Р63 тс. Кристаллическая структура высокотемпературной разупорядоченной а-фазы имеет объем-ноцентрированную кубическую Г анионную подре-шетку,а 2 катиона серебра в кубической ячейке могут занимать 3 типа позиций: 12ё - тетраэдрических, 24к - тригональных и 6Ь - октаэдрических. Молекулярно-динамические (МД) расчеты [2] показывают, что катионы серебра действительно распределены по этим позициям, при этом основной максимум плотности распределения катионов отвечает 12ё позициям, второй по интенсивности - 24к позициям, позициям 6Ь отвечает минимум распределения. По этой причине во многих работах при моделировании динамических и термодинамических свойств a-AgI ограничиваются рассмотрением только 12ё позиций для движения катионов Ag+.

На рис. 5 представлена модельная структура а-AgI. Тетраэдрические пустоты в Гс решетке можно

разделить на 6 ГС подрешеток, нумерация которых представлена на рис. 5. Анализ структуры этих 6 под-решеток показал, что из них можно образовать 3 простых тетрагональных Гподрешетки с более высокой

трансляционной симметрией, чем базовые Гс подрешетки, путем объединения подрешеток 1-4, 2-5, 3-6 (а, 6, в на рис. 5). Каждая Г подрешетки анизотропна, однако, в совокупности три Гподрешетки сохраняют кубическую симметрию составного кристалла.

В ряде работ, основанных на различных версиях МД, исследованы динамическая структура и термодинамические свойства суперионной подрешетки серебра. При высоких температурах среднеквадратичные смещения достаточно велики и изотропны. С понижением температуры среднеквадратичные смещения и диффузия оказываются анизотропными. В частности, при определенных моделях потенциала взаимодействия частицы Ag+ распределяются по некоторым подрешеткам, представленным на рис. 5а, 6, в, в то время как другие оказываются пустыми [6].

В [7] компьютерное моделирование беспорядка в AgI было выполнено с использованием более реалистичного потенциала, чем в [6], учтена также и перестройка подрешетки I. Показано, что при высоких температурах ионы Ag+ начинают конденсироваться по отдельным подрешеткам. В частности получено, что ионы Ag+ оказываются в локальном кристаллическом моноклинном окружении, которое не встречается в модели [6]. Получено также распределение ионов Ag+ по Г подрешеткам, представленным на рис. 5а,

6, в, которое отличается от модели [6]. Это говорит о том, что конкретный вид упорядочения Ag+ существенно зависит от потенциала взаимодействия МД моделирования. Однако следует отметить, что в любом случае при понижении температуры зарождается анизотропная фаза, наблюдаемая диффузия также оказывается анизотропной. Именно появление анизотропии в кубической фазе и создает предпосылки для а^в перехода.

Рис. 5. Анизотропные подрешетки в дробно занятых узлах соединений a-AgI, Ag2X

Кристаллическая структура некоторых кубических фаз Ag2X (X = 8, 8е, Те) оказывается аналогичной структуре высокотемпературной фазе AgI - а-AgI, только в отличие от a-AgI в междоузлиях Ag2X располагаются не 2, а 4 иона Ag+. Это приводит к увеличению кулоновского отталкивания в катионной подрешетке и распределение катионов по междоузли-

ям несколько отличается от того, которое имеет место в a-AgI. Усиливается неравномерность заселения серебра в ё и к позициях; так, в Ag2S серебро занимает главным образом тетраэдрические позиции, причем с ростом температуры эта тенденция усиливается. Вследствие сильной корреляции в катионной подсистеме кубических модификаций халькогенидов серебра

наблюдаются анизотропные флуктуации. Так, в Ag2S в Ag-подсистеме по данным рентгеновского и нейтронного диффузного рассеяния образуются динамические микродомены сигарообразной формы, диаметром ~ а, где а - параметры решетки, вытянуты вдоль направления <100> со временем жизни ~ 5х10-11 с [2].

Соединения Ag3SI и Ag3SBr называют обратными перовскитами (антиперовскитами), так как в них места анионов в перовскитах занимают катионы. Эти соединения являются одними из лучших ионных проводников. Схемы фазовых переходов для Ag3SI и Ag3SBr следующие [2]:

Ag3БІ (а - фаза)0^ (ІшЗт) 519К (в - фаза)0|1(Рт3т)157К (у - фаза)С^(ЛЗ),

AgSBr (фаза) разлагается выше '70:5К (в - фаза)О^(Рт3т)128К (у - фаза^1^ (Стст, С^16).

Структура а-Ag38I высокотемпературной полиморфной модификации образует Гс элементарную ячейку благодаря тому, что анионы 8 и I разупорядо-чены по позициям (0,0,0), (/г,1/,1/), ионы серебра распределены так же, как и в a-AgI по позициям 12ё (42т). Для этих позиций ближайшими соседями являются не только позиции на данной грани, но и на соседней. Вот почему проводимость в a-AgI и а-Ag3SI много больше, чем в в-Ag38I и в-Ag38Br.

Исследование в-Ag38I показало, что анионы 8 и I сосредотачиваются по позициям (0,0,0) и (/,/,/). Шесть ионов серебра распределены по 12к позициям с точечной симметрией (тт2): х / 0, 0 х /, / 0 х, х 0

/, / х 0, 0 / х, X / 0, 0 X /, / 0 X , X 0 /, / X 0,

0 / X, где х = 0.390, отклонение от центра граней около 0.5А. Подобную же структуру имеет и в^^гс, х = 0.385.

Размеры ионов серебра не позволяют расположить на грани более одного атома. Таким образом, в кубических фазах Ag3SI и Ag3SBr ионы серебра располагаются в позициях, близких к центрам граней, так что в итоге образуется структура, близкая к кубической перовскитной, что и оправдывает название «ан-типеровскит». В этом случае мы возвращаемся к приближенным подрешеткам, характерным для классического перовскита (рис. 4а, 6), только для СаТЮ3 эти подрешетки были точными.

4. Заключение

В суперионных соединениях МеБ2, Ме02 с решеткой флюорита результаты как обработки экспериментальных данных, так и молекулярно динамическое моделирование показывают, что доминирующий механизм суперионной проводимости обусловлен движением вакансий, образующихся при возникновении АФД. Вместе с тем существенный вклад в возникновение суперионной проводимости вносит и фононная подсистема, в том числе ангармонизм решеточных колебаний. Предлагаемый нами подход к описанию колебательных спектров МеБ2и Ме02 с решеткой флюорита путем выделения из полного фононного спектра кристалла колебательных состояний подре-шеток позволяет уже в гармоническом приближении исследовать различия в динамическом поведении подрешеток. В частности, температурное поведение среднеквадратичных смещений атомов подрешеток показывает, что смещения анионов имеют большие

значения, чем катионов и быстрее растут с ростом температуры, что и приводит в конце концов к разу-порядочению именно подрешетки анионов.

Особую роль в суперионной проводимости играют НОФ, появление которых обусловлено свертыванием оптических ветвей из ЗБ подрешетки аниона в ЗБ кристалла (рис. 1в). В результате свертывания и взаимодействия «свернутых» ветвей возникает сильная дисперсия и значительное уширение оптической части спектра. Уширение полосы оптических фононов приводит к тому, что численные значения НОФ существенно отличаются от значений частот оптических фононов в верхней части спектра. Кроме того, значения НОФ по рядам соединений коррелируют со значениями температур перехода в суперионное состояние.

В кубических модификациях кристаллов со структурой перовскита на основе анионных Гс под-решеток могут быть сформированы ромбические базоцентрированные Г'с ( = а, Ь, с) подрешетки, трансляционная симметрия которых выше трансляционной симметрии кристалла, однако, точечная симметрия ниже. При этом каждая подрешетка анизотропна, но симметрия анионной подрешетки в целом - кубическая. Такое формальное представление структуры перовскита выявляет внутреннюю (скрытую) симметрию анионной подрешетки, которая будет проявляться в ветвях спектров элементарных возбуждений, сформированных преимущественно из анионных состояний, если такие возбуждения реализуются. Со структурной точки зрения выделение анизотропных подрешеток связано с наличием потенциальной возможности в кубических кристаллах фазовых переходов в низкосимметричные структуры, которые для структуры перовскита представлены на рис. 46, в. Такие возможности проявляются также и в структурах замещения в анионной подрешетке, в частности, в системах СбРЬР2С1 и СбРЬРС12 рентгеноструктурный анализ показал структуру типа перовскита со слабыми ромбическими искажениями.

В кубических модификациях серебросодержащих суперионных кристаллов с катионной проводимостью анизотропные подрешетки могут быть выделены на основе междоузельных позиций, по которым происходит миграция ионов серебра. Эти анизотропные подрешетки для a-AgI, Ag2X (X = 8, 8е, Те) представлены на рис. 5а, 6, в. Для этих кристаллов в ряде работ, основанных на различных версиях МД, исследованы динамическая структура и термодинамические

свойства суперионной подрешетки серебра. При высоких температурах среднеквадратичные смещения достаточно высоки и изотропны. С понижением температуры среднеквадратичные смещения и диффузия оказываются анизотропными и при определенных моделях потенциала взаимодействия частицы Ag+ распределяются по некоторым подрешеткам, представленным на рис. 5а, б, в, в то время как другие оказы-

ваются пустыми. Появление анизотропии в кубической фазе и создает предпосылки для а^в перехода.

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

Литература

1. Воронов, В. Н. Ионная подвижность и свойства соединений ABX3 типа перовскита I В. Н. Воронов II Препринит № ОООФ. - Красноярск: Институт физики СО РАН. - 2006. - 64 с.

2. Иванов-Шиц, А. Ионика твердого тела: в 2 т. I А. Иванов-Шиц, И. В. Мурин. - Т. 2. - СПб.: Изд-во С. - Петерб. ун-та, 2010. - 1000 с.

3. Поплавной, А. С. Механизмы суперионного переноса в кристаллах I А. С. Поплавной. - Кемерово: Изд. КемГУ. - 2009. - 182 с.

4. Теория физических и физико-химических свойств сложных кристаллических соединений с различным типом химической связи; под общ. ред. А. С. Поплавного. - Кемерово: Изд. КемГУ. - 2012. - 400 с.

5. Лебедев, А. И. Abinitio расчеты фононных спектров в кристаллах перовскитов ATiO3 (А = Ca, Sr, Ba, Cd, Ra, Zn, Mg, Ge, Sn, Pb) I А. И. Лебедев II Физика твердого тела.- 2009. - Т. 51. - № 2.- С. 341 - 350.

6. Szabo, G. The lattice gas model on tetrahedral sites of a bcc lattice: anisotropic diffusion in the imtermediate phase I G. Szabo, J. Kertesz II J. Phys. C: Solid State Phys. - 1986. - V. 19. - P. L273 - L277.

7. Madden, P. A. Ordering of the Silver ions in a-AgI: A mechanism for the a^P phase transition I P. A. Madden, K. F. O'Sullivan, G Chiarotti II Phys. Rev. B. - 1992. - V. 45. - № 18. - P. 10206 - 10212.

Информация об авторах:

Нефёдова Дарья Сергеевна - магистрант физического факультета КемГУ, 8(3842)58-31-95,

darinuca@mail.ru.

Daria S. Nefyodova - Master’s Degree student at the Faculty of Physics, Kemerovo State University.

Николаева Елена Владимировна - кандидат физико-математических наук, руководитель бизнес-инкубатора КемГУ, 8(3842)58-31-95, elena eruk@mail.ru.

Elena V. Nikolaeva - Candidate of Physics and Mathematics, Head of Business Incubator, Kemerovo State University.

Поплавной Анатолий Степанович - доктор физико-математических наук, профессор, заведующий кафедрой теоретической физики КемГУ, 8(3842)58-31-95, popl@kemsu.ru.

Anatoly S. Poplavnoy - Doctor of Physics and Mathematics, Professor, Head of the Department of Theoretical Physics, Kemerovo State University.

Федорова Татьяна Петровна - кандидат физико-математический наук, доцент кафедры общей физики КемГУ, 8(3842)58-31-95, kirienko@kemsu.ru.

Tatiana P. Fedorova - Candidate of Physics and Mathematics, Assistant Professor at the Department of General Physics, Kemerovo State University.