Авторефераты по всем темам  >>  Авторефераты по физике  

На правах рукописи

Ивлиев Михаил Петрович

ТЕОРИЯ СТРУКТУРНЫХ ФАЗОВЫХ ПЕРЕХОДОВ

С НЕСКОЛЬКИМИ ПАРАМЕТРАМИ ПОРЯДКА

В КРИСТАЛЛАХ С ОКТАЭДРИЧЕСКОЙ

АНИОННОЙ ПОДРЕШЁТКОЙ

Специальность:

01.04.07 - физика конденсированного состояния

Автореферат диссертации

на соискание ученой степени

доктора физико-математических наук

Ростов-на-Дону

2011 

        Работа выполнена в отделе кристаллофизики Научно-исследовательского института физики Южного федерального университета

Научный консультант:  доктор физико-математических наук, 

                              профессор САХНЕНКО Владимир Павлович

Официальные оппоненты:  доктор физико-математических наук,

              профессор ЗИНЕНКО Виктор Иванович 

                              Институт физики имени Л. В. Киренского РАН

         доктор физико-математических наук,

  профессор ДАРИНСКИЙ Борис Михайлович

  Воронежский государственный университет 

  доктор физико-математических наук, старший

  научный сотрудник ТОРГАШЕВ Виктор Иванович

  Южный федеральный университет

Ведущая организация:       Воронежский государственный

  технический университет 

       Защита диссертации состоится 23 декабря 2011 года в 1400 часов на заседании диссертационного совета Да212.208.05 по физико-математическим наукам (специальность 01.04.07) при Южном федеральном университете в здании НИИ физики ЮФУ по адресу: 344090, г. Ростов-на-Дону, пр. Стачки, 194, ауд. 411

       С диссертацией можно ознакомиться в Зональной научной библиотеке ЮФУ по адресу: г. Ростов-на-Дону, ул. Пушкинская, 148

       Автореферат разослан л__ ноября 2011 года

               Отзывы на автореферат, заверенные подписью рецензента и печатью учреждения,  просим  направлять  ученому  секретарю диссертационного совета Да212.208.05 при ЮФУ по адресу: 344090, г. Ростов-на-Дону, пр. Стачки, 194, НИИ физики ЮФУ

Ученый секретарь диссертационного

совета Д 212.208.05  при ЮФУ  Гегузина Г.А.

                       ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Одним из наиболее распространённых типов структурных превращений в кристаллах являются ФП порядок - беспорядок, обусловленные упорядочением частиц по нескольким кристаллографически эквивалентным позициям (КЭП) [1]. В таких соединениях, как правило, наблюдается несколько упорядоченных фаз, возникаюших вследствие перераспределения частиц по набору КЭП. В этом случае описание фазовых превращений сводится к поиску оптимального распределения частиц по упомянутым позициям (модель Френкеля [2]). Такая модель использовалась для интерпретации структурных ФП в простых ионно-ковалентных кристаллах, например, в перовскитах [3, 4].

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

В рамках метода ПСП упорядочения, возникающие в системах с большим числом КЭП, характеризуются несколькими многокомпонентными  параметрами порядка (ПП), равновесные значения которых определяются системой трансцендентных уравнений самосогласования [7]. Сложность такой системы уравнений заставляет прибегать к дополнительным упрощениям уже внутри самого приближения. В большинстве работ либо рассматриваются конкретные вещества и исследуются лишь отдельные упорядочения, наблюдаемые именно в этих кристаллах [7], либо применяются дополнительные, трудно контролируемые приближения, позволяющие упростить систему уравнений самосогласования. Кроме того области в пространстве термодинамических и модельных параметров, расположенные вблизи линий ФП второго рода из симметричной фазы, удаётся исследовать методами феноменологической теории [8].

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

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

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

Объекты исследования: кристаллы с октаэдрической анионной подрешёткой. Основу этой обширной группы кристаллов составляют перовскиты, упорядоченные перовскиты (включая эльпасолиты, криолиты), слоистые перовскиты и ряд других соединений.

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

Научная новизна. В работе впервые:

1. Получена диаграмма фазовых состояний 8-мин. кубической модели: установлено наличие трёх сегнетоэлектрических фаз, определены факторы, влияющие на расположение фаз, и показано, что последовательность ФП, наблюдаемых в BaTiO3 и KNbO3 может быть описана на основе этой модели. 

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

3. Получена диаграмма фазовых состояний 3-мин. двухподрешёточной модели: на границе с симметричной фазой помимо трикритической имеются две четырехфазные точки, а в области Т 0 межфазные границы трёх упорядоченных фаз сходятся асимптотически к двум разным линиям.

4. Для параметра порядка,  относящегося к представлению М2+, на различных сечениях получены диаграммы фазовых состояний 3-мин. восьмиподрешёточной модели: на границе с симметричной фазой помимо трикритической имеются несколько трехфазных точек, число которых изменяется от нуля до двух.

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

6. Установлен и описан в рамках статистической теории один из возможных механизмов изоструктурных (изосимметричных) фазовых переходов для кристаллов с несколькими наборами кристаллографически эквивалентных позиций. 

7. Установлены условия, при которых возникают наиболее распространённые последовательности ротационных упорядочений в кристаллах со структурой перовскита (пр. гр. Oh1), упорядоченных перовскитов (включая эльпасолиты, криолиты), антифлюоритах (пр. гр.Oh5) и слоистых перовскитах (пр. гр. D4h1).

8. Выявлено влияние фактора толерантности и заряда катиона А на формирование фазовых состояний при ротационных фазовых переходов в кубических перовскитах АВО3. 

9. Показано, что упорядоченное фазовое состояние в MnF3 образуется вследствие совместной конденсации ротационного и орбитального параметров порядка при фазовом переходе триггерного типа.

Основные положения, выносимые на защиту

1. Диаграммы фазовых состояний 4-минимумной квадратной, 6- и 8- минимумных кубических моделей, включающие области вблизи Т = 0, содержат сегнетоэластическую и три сегнетоэлектрические фазы. На большинстве межфазных границ имеются трикритические точки, разделяющие области фазовых переходов первого и второго рода.

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

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

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

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

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

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

8. На базе феноменологической теории разработана термодинамическая модель, способная описать всю совокупность ротационных и сегнетоэлектрических фазовых переходов, подобным наблюдаемым в твёрдом растворе Na1-xKxNbO3 при Т > 300 С, 0.03 < x < 0.2. Из анализа данных по структуре и температурных зависимостей диэлектрической проницаемости для твёрдых растворов Na1-xKxNbO3 различного состава определены основные параметры модели, что позволило на качественном уровне воспроизвести x - T диаграмму фазовых состояний, наблюдаемую в Na1-xKxNbO3 в области Т > 300 С, 0.03 < x < 0.2.

Научная и практическая значимость

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

Апробация работы происходила на Междунар. сем. по физике сенетоэластиков, 1994, Воронеж; Междунар. научно-практической конф. Фундаментальные проблемы пьезоэлектроники, 1995, Азов; XV Всерос. конф. по физике сегнетоэлектриков, 1999, Азов; Междунар. симп. Фазовые превращения в твёрдых растворах и сплавах, 2003, Сочи; 7-ом  Междунар. симп. Порядок, беспорядок и свойства оксидов, 2004, Сочи; XVII Всерос. конф. по физике сегнетоэлектриков, 2005, Пенза; 8-ом  Междунар. симп. Порядок, беспорядок и свойства оксидов, 2005, Сочи; 5 Междунар. сем. по физике сенетоэластиков, 2006, Воронеж; XVIII Всерос. конф. по физике сегнетоэлектриков, 2008, Санкт-Петербург; 6 Междунар. сем. по физике сенетоэластиков, 2009, Воронеж; XIX Всерос. конф. по физике сегнетоэлектриков, 2011, Москва.

ичный вклад автора. Основные темы и направления работы обсуждались и конкретизировались с сотрудниками Департамента физики Южного федерального университета: д-ром физ-мат. наук, проф. Сахненко В.П. и канд. физ-мат. наук Тимониным П.Н.

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

Благодарности. Автор глубоко признателен профессору Сахненко В.П. за помощь и ценные советы при работе над диссертацией. Искренне благодарит коллег Тимонина П.Н., Раевского И.П. и Резниченко Л.А. за участие в обсуждениях тем, затронутых в диссертации. Их помощь способствовала плодотворной работе, в процессе которой и были получены  все результаты, выносимые на защиту.

       Структура и объем работы.  Диссертация состоит из введения, 5 разделов, заключения и списка цитируемой литературы из 145 наименований, изложенных на 263 страницах, включая 87 рисунков и 3 таблицы.  В конце автореферата приведен список основных публикаций автора по теме диссертации с литером А при каждом наименовании. 

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первом разделе приводятся примеры объектов, в которых реализуются многоминимумные модели.  Отмечено, что за последние годы число таких соединений значительно возросло. В частности, экспериментально установлено, что в кубической ячейке параэлектрической фазы целого ряда перовскитов с общей структурной формулой ABO3 катионы А и (или) В занимают нецентросимметричное положение. Ряд работ свидетельствует о том, что в кубической фазе перовскитов KNbO3, NaNbO3 и BaTiO3 катионы Nb и Ti распределены по восьми [10, 11], а в PbTiO3 катион Ti - по шести позициям [12].

Ранее, за редким исключением, предполагалось, что  в кристалле имеется только один набор КЭП, однако исследования, выполненные за последние годы, показали, что в ряде кристаллов для упорядочивающихся при фазовом переходе атомов имеется несколько наборов минимумов потенциальной энергии, т.е. несколько наборов КЭП, характеризуемых разной энергией. Существование нескольких наборов КЭП обнаружено (или предполагается по косвенным признакам) в перовскитах, в суперионных проводниках CuI, AgI и др. [1], в LiTaO3, в Hg2Cl2 [5], в кристаллах A3MOxF6-x (x = 1,2,3; M = Ti, Nb, W, Mo и др.) и A2MO2F4 (M = W, Mo), где А - одновалентный катион [6].

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

Обсуждаются специфические особенности физических свойств, наблюдаемые при сегнетоэлектрических фазовых переходах в перовскитах, часть из которых может свидетельствовать пользу того, что сегнетоэлектрические ФП в перовскитах - ФП типа порядок - беспорядок [3, 4, 11, 12] однако другая часть свидетельствует  пользу того, что сегнетоэлектрические ФП в перовскитах - ФП типа смещения. Отмечено, что в последние годы сформировалась точка зрения, согласно которой в перовскитах имеет место своеобразная гибридизация процессов размягчения кристалла по одной из мод решёточных колебаний и перераспределения частиц по КЭП [10]. При этом формирование упорядоченного состояния будут определять оба процесса.

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

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

Для описания упорядочений, возникающих в таких системах, требуется, как правило, несколько многокомпонентных параметров порядка. Их равновесные значения можно найти, решив систему уравнений самосогласования (СУС) и выбрать тот набор компонент ПП, которому соответствует абсолютный минимум термодинамического потенциала (ТП) при данных значениях модельных и термодинамических параметров [7, A1]. Сложность такой системы уравнений, а также отсутствие общих методов решения, побуждают прибегать к дополнительным упрощениям. Наиболее распространённые способы упрощённого описания следующие: при исследовании упорядочений учитывают только тот ПП, по которому система теряет устойчивость (критический), отбрасывая те параметры из конфигурационного пространства модели, которые могли бы появиться вследствие нелинейного взаимодействия с критическим, то есть вторичные:

- при рассмотрении конкретных веществ исследуют лишь отдельные упорядочения, описываемые небольшим числом компонент критического и вторичного ПП (обычно одной компонентой), наблюдаемые именно в этих веществах [7];

-  пренебрегают реальной кристаллической анизотропией. 

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

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

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

Для решения этой задачи в диссертации используется следующая схема:

- для критического ПП, с помощью теоретико-групповых методов выделяется полный конденсат, то есть совокупность ПП из конфигурационного пространства модели, порождаемых нелинейным взаимодействием с критическим (вторичных, несобственных ПП) [13];

- из анализа асимптотических областей, находящихся вблизи линии ФП второго рода из симметричной фазы [9], а также в окрестности Т 0, определяется набор фаз, описываемых критическим ПП, существующих в широкой области значений термодинамических и модельных параметров и составляющих основу диаграммы ФС (в большинстве случаев эти фазы характеризуются небольшим числом независимых параметров);

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

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

Такой подход даёт возможность построить диаграмму фазовых состояний (ФС) модели в широкой области значений термондинамических и модельных параметров. Сопоставляя условия появления тех или иных ФС, можно определить механизмы, ответственные за их возникновение, а также выявить закономерности формирования картины ФС в целом [A2 - A4].

Наличие нескольких наборов КЭП с различной глубиной потенциальных ям существенно усложняет и без того непростую задачу по исследованию статистических свойств модели, поскольку необходимо учитывать возможность перераспределения частиц между КЭП разного типа и формирования ФС набором КЭП с менее глубокими потенциальными ямами. Теоретические  работы, посвящённые исследованию моделей такого типа, практически отсутствуют.

Применяя вышеупомянутый подход, в широкой области значений термодинамических и модельных параметров исследованы статистические свойства простых многоминимумных моделей: в главе 2 - четырёхминимумной  (4-мин.) квадратной,  шести-, восьми- и двенадцатиминимуных (6-, 8- и 12-мин.) кубических с одной подрешёткой; в главе 3 - трёхминимумной (3-мин.) кубической с двумя, четырьмя и восьмью подрешётками. В главе 4 проанализированы статистические свойства сложных составных  многоминимумных моделей: 1 + 2-минимуных и 1+4-минимумных, а также исследованы специфические особенности термодинамических свойств кристаллов, описываемых этими моделями. 

В качестве объектов теоретического исследования структурных фазовых переходов в диссертации рассмотрены соединения с октаэдрической анионной подрешёткой. Основу обширной группы таких кристаллов составляют перовскиты, упорядоченные перовскиты (включая эльпасолиты и криолиты), слоистые перовскиты и ряд других соединений [9]. В кристаллах этой группы одним из наиболее распространённых типов структурных превращений являются ФП, обусловленные смещениями анионов, которые можно истолковать как повороты октаэдров (ротационные ФП, "смятие", tilt). В настоящее время накоплен большой объём экспериментального материала по ротационным ФП в этих кристаллах, проведена его систематизация и предварительный теоретико-групповой анализ, позволивший установить ПП, описывающие ротационно упорядоченные фазы. Также был выделен набор наиболее распространённых (базовых) последовательностей ротационных упорядочений [9].

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

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

Для кристаллов со структурой перовскита сопоставление условий реализации различных базовых последовательностей позволило оценить влияние кристаллохимических параметров компонент соединения таких, как заряд катиона А и соотношение ионных радиусов, характеризуемое фактором толерантности, на формирование ФС [A5].

Также исследовано взаимное влияние поляризации и ротационных упорядочений различного типа на формирование упорядоченных ФС. В качестве модельного объекта выбраны твёрдые растворы Na1-xKxNbO3 (NKN), поскольку в них наблюдаются в различных сочетаниях как ротационные, так и сегнетоэлектрические ФП [14]. В диссертации предложена термодинамическая модель [A6], способная описать всю совокупность ротационных и сегнетоэлектрических ФП, подобных тем, что наблюдаются в твёрдом растворе NKN при Т > 300 C, 0.03 < x < 0.2. Её основные параметры получены из анализа данных по структуре и температурных зависимостей диэлектрической проницаемости для твёрдых растворов NKN различного состава. На базе предложенной модели проанализированы особенности диэлектрических свойств и исследовано влияние внешних воздействий - давления и электрического поля на ФС.

Во втором разделе исследуются статистические свойства простых многоминимумных моделей с одной подрешёткой (без мультипликации ячейки): 4-мин. квадратной, 6-, 8- и 12-мин. кубических [A1, A2, A4, A7 - A9]. Такой выбор обусловлен тем, что 6- и 8-мин. кубические модели являются наиболее известными и часто используемыми многоминимумными моделями [3, 4, 7].  Четырехминимумная квадратная модель - это, фактически, более простой (плоский) вариант 6-мин. кубической модели. Она является хорошим модельным объектом для демонстрации в деталях всего процесса исследования статистических свойств таких моделей. Статистические свойства 12-мин. кубической модели [7] в широкой области значений термодинамических и модельных параметров ранее никогда не исследовались. Эта модель интересна тем, что её конфигурационное пространство содержит ПП, которые по трансформационным свойствам идентичны поляризации и тензорам одноосных и сдвиговых напряжений.

       Исследование в широкой области значений термодинамических и модельных параметров статистических свойств упомянутых моделей проведено в приближении Горского - Брэгга - Вильямса, при этом выявлен ряд общих закономерностей формирования  картины фазовых состояний (ФС).

В подразделе 2.1 проведён анализ 6-мин. кубической модели [A1, A7, A8]. В качестве объекта исследования рассматривается кристалл с  кубической решёткой  симметрии Оh1, содержащий подсистему частиц, каждая из которых с одинаковой вероятностью может находиться в одной из шести КЭП, отвечающих минимумам их потенциальной энергии в кристаллической решётке. Минимумы немного смещены от узла вдоль направлений типа [001].

Вероятности заполнения минимумов, характеризуются следующим набором

n1,2 = (1 + 1 1); n3,4 = (1 Ц  1 + 2 2)                                 n5,6 = (1 Ц  1 - 2 3).                                                        

Функции , выполняют роль параметров порядка (ПП),  преобразуется по неприводимому представлению Eg, а - по Т1u. Равновесные значения параметров порядка , можно найти решив систему уравнений состояния (СУС):

  (i = 1,2) ,    (j =1,2,3)

и выбрав то решение, которому соответствует абсолютный минимум F при  заданных значениях  A, B и Т.

Неравновесный термодинамический потенциал (ТП) 6-мин. кубической модели в приближении среднего поля в варианте Горского - Брэгга - Вильямса (ГБВ) в расчёте на одну частицу имеет вид

= [() + ()]  - ln 6 + [(1 + 1 1) ln (1 + 1 1)+

  + (1 Ц  1 + 2 2) ln (1 - 1 + 2 2) +                        

+ (1 Ц  1 - 2 3) ln (1 Ц  1 - 2 3)] .                        (1)

Знак л предполагает суммирование выражений сначала со знаком л+, а затем со знаком л-.  Константы А и В - феноменологические параметры теории. Функции , выполняют роль параметров порядка (ПП),  преобразуется по неприводимому представлению Eg, а - по Т1u.                

Общий вид диаграммы ФС, описываемой термодинамическим потенциалом (1), приведён на рис. 1, где I - кубическая фаза с симметрией Oh1; II - тетрагональная сегнетоэластическая фаза с симметрией D4h1; III, IV, V - сегнетоэлектрические (СЭ) фазы, соответственно, тетрагональная - C4v1, ромбоэдрическая - C3v5 и моноклинная - Сi1.        Диаграмма ФС (рис. 1) содержит сегнетоэластическую и три СЭ фазы. В области, где критическим является ПП , при В < 0 реализуется фаза III Ц(1,0, 1,0,0), при В > 0 реализуется фаза IV Ц(0,0,,,). Таким образом, с ростом В фаза III либо через промежуточную фазу V - (1,0,1,2,2), либо непосредственно заменяется фазой IV.

                       

Рисунок 1 - Диаграмма ФС 6-мин. модели, описываемой ТП (1). Здесь и далее сплошными обозначены линии ФП первого рода (ФП 1), штриховыми - второго рода. Линия В = - А - асимптота, к которой сходятся межфазные границы при Т 0

Диаграмма ФС (рис. 1), это диаграмма для недеформируемой, зажатой решётки. Учёт деформируемости решётки относительно сдвиговых деформаций даже при слабом стрикционном взаимодействии γ2/c = 0.04.B ( - константа стрикции, с Цсоответствующий модуль упругости), приводит к тому, что вместо моноклинной СЭ фазы на диаграмме фазовых состояний появляется ромбическая СЭ фаза (рис. 2).

Также в подразделе 2.1 проведён анализ четырёминимумной (4-мин.) квадратной модели [A1]. В качестве объекта исследования рассматривается слоистый кристалл с тетрагональной решёткой  симметрии D4h, содержащий подсистему частиц, каждая из которых с одинаковой вероятностью может находиться в одной из четырёх КЭП, отвечающих минимумам их потенциальной энергии в кристаллической решётке. Минимумы находятся в слоях ортогональных оси 4-го порядка, они немного смещены от узла в направлениях [ 1,0,0]  и [0, 1,0].

               

Рисунок 2 - Диаграмма ФС 6-мин. модели при учёте деформируемости решётки относительно сдвиговых деформаций. Тонкая сплошная линия - термодинамический путь, воспроизводящий последовательность  ФП, наблюдаемых в BaTiO3 и КNbО3

       Неравновесный ТП 4-мин. модели имеет вид:   

где А, В - параметры модели. Функции , выполняют роль параметров порядка,  преобразуется по неприводимому представлению В1g, а 1,2 - по Eu. Система уравнений состояния имеет пять типов решений, отличающихся по симметрии: 1.( =1= 2 = 0) - D4h; 2. ( , 1= 2 = 0) - D2h;  3. (, 1,0) - С2v;  4.( 0, , ) - С2v; 5. (, 1, 2) - Сs.

Диаграмма ФС 4-мин. квадратной модели имеет вид подобный тому, что изображён на рис. 1. Она содержит все пять возможных по симметрии фазовых состояний, т.е. фазы I - V. В области, где критическим является ПП , при В < 0 реализуется фаза III - (,1,0), при В > 0 реализуется фаза IV - (0,,). Таким образом, с ростом В фаза III либо через промежуточную фазу V Ц  (,1,2), либо непосредственно заменяется фазой IV.

В подразделе 2.2 проведён анализ 8-мин. кубической модели [A2, A9]. Рассматривается кристалл с простой кубической решёткой (пр. гр. Oh1 (Z=1)), содержащий подсистему частиц, каждая из которых с одинаковой вероятностью может находиться в одной из восьми КЭП, отвечающих минимумам их потенциальной энергии в кристаллической решётке. Минимумы немного смещены в направлениях типа [111] от узла в центросимметричной позиции.

Неравновесный ТП 8-мин. модели в приближении ГБВ имеет вид: 

       F =  A() + D() + C 2 +                                 

+[(1+ e1 + e2 +e3 ( + 1 + 2 + 3)) ln(1+ e1 + e2 +e3 ( + 1 + 2 + 3)+

+ (1+ e1 - e2 - e3 ( + 1 - 2 - 3)) ln(1+ e1 - e2 - e3 ( + 1 - 2 - 3))+

+ (1 - e1 + e2 - e3 ( - 1 + 2 - 3)) ln(1 - e1 + e2 - e3 ( - 1 + 2 - 3)+

+(1 - e1 - e2 +e3 ( - 1 - 2 +3)) ln(1 - e1 - e2 +e3 ( - 1 - 2 + 3))],  (2)

где  A, B, C - модельные константы. Функции , , е выступают в качестве параметров порядка, ОT1u, eОT2g, ОA2u.

Диаграмма ФС 8-мин. кубической модели, описываемая ТП (5) при С = 0 имеет вид, изображённый на рис. 3, где I - кубическая фаза с симметрией Oh1; II - (0,0,0,e,e,e,0) ромбоэдрическая сегнетоэластическая фаза с симметрией D3d5; III, IV, VI - СЭ фазы, соответственно: ромбоэдрическая - (,,,e,e,e,), C3v5; тетрагональная - (1,0,0,0,0,0,0), C4v1; ромбическая - (0,2,2,e1,0,0,0), (C2v14) Двигаясь вдоль термодинамического пути 1, обозначенного тонкой сплошной линией (рис. 3), можно воспроизвести последовательность ФП, наблюдаемых в KNbO3: IIV VI III. Отличие от реальной ситуации лишь в том, что на диаграмме ФП I Ц  IV - переход второго рода, а в KNbO3 и в BaTiO3 это переход первого рода.

       

  Рисунок 3 Ц  Диаграмма ФС 8-мин. модели, описываемая ТП (2) при С = 0

Учёт деформируемости решётки относительно одноосных деформаций приводит к тому, что на диаграмме фазовых состояний, на границе с симметричной фазой появляется вторая трикритическая точка и линия ФП первого рода между фазами I - IV. Показано, что при стрикционном взаимодействии 2/c0 = 2D  ( Цконстанта стрикции, с0=(1/2)(c11 - c12) - модуль упругости) имеется  термодинамический путь, двигаясь вдоль которого удаётся  воспроизвести всю последовательность ФП, которая наблюдается в KNbO3, причём все ФП - переходы первого рода.

       В подразделе 2.3 проведён анализ статистических свойств 12-мин. кубической модели [7, A4].        Конфигурационное пространство этой модели содержит ПП, которые по своим трансформационным свойствам идентичны тензорам одноосных и сдвиговых деформаций, а также вектору поляризации. Можно предположить, что такое сочетание ПП позволит в рамках единого подхода описывать ФП первого и второго рода из кубической в наиболее распространённые сегнетоэлектрические фазы - С2v, С3v и  С4v. Фазовые переходы такого типа наблюдаются в твёрдых растворах PbZr1-xTixO3 (ЦТС), (1-х)PbMg1/3Nb2/3O3 - xPbTiO3 (x>0.25), (1-х)PbZn1/3Nb2/3O3 - xPbTiO3 и других. Для оценки возможностей 12-мин. модели, исследована её применимость для описания х-Т диаграммы фазовых состояний, аналогичной наблюдаемой в ЦТС.  Выбор ЦТС обусловлен тем, что его х-Т диаграмма хорошо известна, а кроме того она содержит ряд особенностей, описать которые  на основе других многоминимумных моделей, не прибегая к введению дополнительных взаимодействий, невозможно.

Рассматривается кристалл с простой кубической решёткой (пр. гр. Oh1 (Z=1)), содержащий подсистему частиц, каждая из которых с одинаковой вероятностью может находиться в одной из двенадцати КЭП, отвечающих минимумам их потенциальной энергии в кристаллической решётке. Минимумы немного смещены в направлениях типа [011] от узла в центросимметричной позиции, в котором пересекаются три оси С4, их расположение и нумерация приведены на рис. 4.        

Неравновесный ТП 12-мин. модели имеет вид:

       12F =  A + B(6f12 + 2f22) + D + G +                 

+ T[(1 + 2f1 + e3 (1 + 2 + 1 - 2)) ln(1 + 2f1 + e3 (1 + 2 + 1 - 2)) +

+ (1 + 2f1 - e3 (1 - 2 + 1 + 2)) ln(1 + 2f1 - e3 (1 - 2 + 1 + 2)) +

+ (1 - f1 + f2 + e2 (1 + 3 - 1 + 3)) ln(1 - f1 + f2 + e2 (1 + 3 - 1 + 3)) +

+ (1 - f1 + f2 - e2 (1 - 3 - 1 - 3)) ln(1 - f1 + f2 - e2 (1 - 3 - 1 - 3)) +

+ (1 - f1 - f2 + e1 (2 + 3 + 2 - 3))ln(1 - f1 - f2 + e1 (2 + 3 + 2 - 3)) +

+ (1 - f1 - f2 - e1 (2 - 3 + 2 + 3))ln(1 - f1 - f2 - e1 (2 - 3 + 2 + 3))],         (3)        

означает суммирование по верхнему и по нижнему знакам. A, B, D, G - модельные константы. Функции f, , e, выполняют роль ПП и преобразуются по неприводимым представлениям соответственно  Eg, T1u, T2g и T2u. Параметры порядка f и e можно считать сегнетоэластическими, а ПП - сегнетоэлектрическим. Далее полагаем, что критическим является параметр .        Из сопоставления результатов исследования ТП (3) вблизи границы с симметричной фазой с экспериментальными данными, получен набор условий, который позволил подобрать модельные константы A(x), B(x), D(x) таким образом, чтобы вид диаграммы фазовых состояний, полученной теоретически, вблизи границы с симметричной фазой соответствовал по основным параметрам экспериментальной х-Т диаграмме твёрдых растворов цирконата-титаната свинца.

Рисунок 4  - Расположение и нумерация 12-ти  КЭП

Полагая, что 

Тс(х) = 450 + 450х - 200х2;  А(х) = - 4 Тс(х);                                 В(х) = (0.4 - 2.4х) Тс(х);         D(x) = (- 1.5 + 3.5x) Тс(х)                 (4)

при G(x) = 0.8(4 - 3x)Tc(x) и G(x) = (4 - 3x)Tc(x) получаем диаграммы ФС (рис. 5).  На диаграммах ФС вблизи  линии ФП второго рода из симметричной фазы картина ФС подобна той, что наблюдается у всех рассмотренных ранее моделей. Однако в области x < ~ 1 появляется необычная сегнетоэлектрическая ромбическая фаза Orth II. В области T ~ 0 картина совершенно меняется. Тетрагональная сегнетоэлектрическая фаза либо присутствует в ограниченной области значений модельных параметров, либо отсутствует вообще. Вместо неё присутствуют либо одна Orth II, либо две ромбические сегнетоэлектрические фазы Orth I и Orth II.

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

         

а)                                                        б)

Рисунок 5 - Диаграммы фазовых состояний, описываемые ТП(3) с учётом (4) при: G(x) = 0.8(4 - 3x)Tc(x) - (а); G(x) = (4 - 3x)Tc(x) -(б)

Появление такой необычной фазы как Orth II обусловлено тем, что она возникает не из кубической, а из тетрагональной сегнетоэластической (D4h1). В этой фазе имеет место преимущественное заполнение четырёх КЭП в средней, экваториальной плоскости, ортогональной оси 4-го порядка. Поэтому, когда критическим становится ПП , то перераспределение частиц по КЭП происходит по этим четырём позициям. В результате, как следует из анализа 4-мин. квадратной модели, могут появиться две сегнетоэлектрические фазы. При G(x) = 0.8(4-3x)Tc(x) появляются обе. Одна из них фаза Orth II, другая,  Orth I по структуре подобна ромбической сегнетоэлектрической фазе, возникающей из кубической (в данном случае она неустойчива), однако по происхождению это совсем другая фаза.

В третьем разделе на примере трёхминимумной (3-мин.) кубической модели рассмотрены ФП, при которых возникают фазы с несколькими по-разному упорядоченными подрешётками, т.е. фазовые переходы, сопровождаемые мультипликацией ячейки [A3, A10]. Эта модель является одной из наиболее распространённых многоминимумных многоподрешёточных моделей. Известна большая группа ионно-молекулярных кристаллов (CsSH, CsO2), кристаллов, содержащих октаэдры, искажённые вследствие эффекта Яна-Теллера (KCuF3, KCrF3), упорядочения в которых описываются именно этой моделью. Деформирующие ЯТ ионы переходного элемента можно представить в виде частиц в форме вытянутого цилиндра, которые в симметричной (кубической) фазе могут быть ориентированы вдоль трёх направлений типа [001]. На таком уровне исследований ориентационные ФП в ионно-молекулярных кристаллах и ФП, обусловленные статическим эффектом Яна - Теллера, описываются на основе одной и той же статистической модели. В качестве модельной системы рассматривается кристалл, имеющий простую кубическую решётку с симметрией O(Z=1), в её узлах расположены частицы в форме вытянутого цилиндра, которые с одинаковой вероятностью могут быть ориентированы вдоль одного из трёх направлений типа [001], то есть имеют три КЭП.

В подразделе 3.1 анализируется 3-мин. двухподрешёточная  модель. Исследуются возможные упорядочения, относящиеся к вектору  k = ()(b1 + b2 + b3), (bi - базисные векторы обратной решётки), то есть к точке R пространства векторов обратной решётки [A10]. В этом случае кристалл можно представить в виде двух одинаковых, взаимопроникающих, по-разному упорядоченных подрешёток, узлы которых чередуются по трём направлениям типа [001].

Функции, характеризующие заполнение КЭП для каждой из подрешёток, имеют вид:                 

n1, (I,II) = , n2, (I,II) =

                               n3, (I,II) =  ,                        

I, II - номер подрешётки, i = (1,2,3) -  номер КЭП, относящийся к той, или иной подрешётке.  В столбце л знак л+ соответствует первой подрешётке, а знак л - - второй.

Неравновесный ТП двухподрешёточной трёхминимумной (3-мин.) модели в приближении Горского - Брэгга - Вильямса можно представить:

       F = + + T.        (5)

Функции и выполняют роль параметров порядка, где Eg(k = 0), R2+, k = (1/2,1/2,1/2). Диаграмма ФС, описываемая ТП (5), изображена на рис. 6. Фаза I- (1,0,0,0), фазы II и II1 - (1,0,0,2), фаза III- (1,0,1,0), фаза IV- (1,2,1,2).

               

Рисунок 6 - Диаграмма ФС 3-мин. двухподрешёточной модели, описываемая ТП (5). Линии В = - А и B = -  3A - асимптоты, к которым сходятся межфазные границы при Т 0

Система уравнений состояния может иметь помимо симметричного || = || = 0 (пр. гр. Oh1(Z = 1)) два набора решений. В первый набор входят два решения типа 1, 2 = 0 (пр. гр. D4h1(Z = 1)) и 1, 2 (пр. гр. D2h1(Z = 1)), описывающие сегнетоэластические фазы с упорядочением ферродисторсионного типа [A7]. Во второй набор - три решения типа:  1, 2 = 0, 1 = 0, 2 (пр. гр. D4h18(Z = 4));  1, 2 = 0, 1, 2 = 0 (пр. гр. D4h17(Z = 4)) и 1, 2, 1, 2 (пр. гр. D2h23(Z = 4)), описывающие фазы с упорядочением антиферродисторсионного типа с k = (1/2,1/2,1/2).        

Характерная особенность диаграммы (рис. 6) в том, что фазы  II и II1 имеют одинаковую симметрию. На границе с симметричной фазой С помимо трикритической точки K имеются две четырёх фазные мультикритические точки L и M, а в области низких температур, границы между фазами II - IV,  IV - III и III - IV,  IV Ц  II1 асимптотически сходятся к двум разным линиям.

В подразделе 3.2 анализируются 3-мин. четырёх- и восьмиподрешёточные  модели [A3]. Для исследования упорядочений, описываемых 3-мин. восьмиподрешёточной моделью предложен следующий метод:

- исходная кристаллическая решётка (пр. гр. Oh1) разбивается на две эквивалентные, взаимопроникающие подрешётки с симметрией Oh5;

- исследуются симметрийные и статистические свойства 3-мин. четырёхподрешёточной модели для кристалла с симметрией Oh5.

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

Для параметра порядка,  относящегося к представлению М2+, на различных сечениях получены диаграммы фазовых состояний 3-мин. восьмиподрешёточной модели. Показано, что на границе между симметричной (Oh1) и фазами (D4h5 (Z = 2)) и (Th5 (Z = 8)) соответственно могут  существовать  трикритические точки, также  на границе с симметричной фазой имеются мультикритические 3-х фазные точки, в которых граничат симметричная, тетрагональная и ромбоэдрическая фазы. Количество 3-х фазных точек может меняться от нуля до двух.

       Используя приближение взаимодействия ближайщих соседей, определены условия, при которых наиболее распространённые упорядочения антиферродисторсионного типа (пр. гр. D4h18 и D4h5) могут реализоваться в одном и том же кристалле.

В четвёртом разделе исследуются статистические свойства многоминимумных моделей с несколькими наборами КЭП, отличающимися по энергии (составными  многоминимумными моделями) [A11, A12]. Одной из основных целей этого раздела диссертации является выявление на примере наиболее простых Ц  1 + 2-мин. и 1 + 4-мин. моделей, сохраняющих основные черты, присущие кристаллам с несколькими наборами КЭП, особенностей термодинамического поведения, обусловленных перераспределением атомов между различными наборами КЭП.

В подразделе 4.1 исследованы свойства 1 + 2-мин. модели [A11]. Рассмотрена трёхминимумная модель, состоящая из двух подсистем минимумов - центрального и двух боковых (1+2-мин.). На рис. 7 приведены два варианта рельефа потенциальной энергии частицы для  такой 1+2-мин. модели.

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

Функции, характеризующие вероятности заполнения минимумов:  центрального - p = (1+ ) /2,  боковых - n1 =(1 - )(1+ )/4, n2 =(1 - )(1 - )/4. Функции , выполняют роль параметров порядка, преобразуется по неприводимому представлению A1g , а  - по A1u.

               

Рисунок 7 - Варианты рельефа потенциальной энергии частицы

для  1+2- мин.  модели

Неравновесный термодинамический потенциал трёхминимумной (1+2-мин.) модели в приближении Горского - Брэгга - Вильямса имеет вид: 

                      (6)      

Параметр в основном зависит от разницы глубин потенциальных ям -  центральной и боковых, если > 0 то центральная глубже боковых. Параметр В характеризует соотношение между взаимодействиями внутри подсистем минимумов и между этими подсистемами (одна подсистема состоит из  центральных ям, другая - из боковых). Упрощая ситуацию, можно сказать, что при В < 0 частицы стремятся расположиться в одной подсистеме минимумов, а при В  > 0 - по разным. Параметр А  характеризует взаимодействие между частицами, находящимися в боковых ямах. Так, если А < 0, то частицы стремятся расположиться в боковых ямах одного типа (т.е. находящихся с одной стороны от центральной), а если А > 0, то - в разных.        

Система уравнений состояния

,                                                (7)

имеет два типа решений: 1. (, = 0) - пр. гр. D4h; 2. (, ) - пр. гр.C4v.

Показано, что в случае одного ПП в области > 0, B < 0 в кристаллах, описываемых этими моделями, наблюдаются следующие особенности термодинамического поведения: 

1. В отсутствие фазовых переходов могут наблюдаться максимумы на температурных зависимостях ряда физических характеристик, в частности упругой податливости и избыточной теплоёмкости Ср (рис. 8).

2. Появляются состояния, в которых при Т ~ 0 оказываются частично заполнены КЭП с более высокой энергией, то есть возбуждённые состояния. 

3. В симметричной фазе имеют место изоструктурные фазовые переходы двух типов: когда фазовый переход происходит между фазами, параметры порядка которых принадлежат одной ветви решений системы уравнений состояния, и фазовый переход между фазами, параметры порядка которых принадлежат разным ветвям системы уравнений состояния. Фазовый переход первого типа связан с появлением локальной неустойчивости у каждой из фаз. Фазовый переход второго типа связан с потерей устойчивости лишь одной из фаз. Другая остаётся устойчивой во всём температурном интервале. 

       В случае двух ПП в области > 0, B > 0, А < 0 для k = 0.4 (k = B/) при

0 < s < 1.207 (s = -A/), несмотря на то, что А < 0 ФП не происходит и ПП не появляется. Хотя в области s < ~ 1.207 наблюдаются хорошо заметные аномалии восприимчивости по ПП (рис. 9). В  области  1.207 < s < 1.8 СУС (7) имеет только одно решение (,)1, ответвляющееся от решения (,0) в двух точках. Фаза (,)1 возникает вследствие ФП второго рода (ФП2) из симметричной и существует в ограниченном интервале температур. Затем снова посредством ФП2 происходит возврат в симметричную фазу, которая остаётся устойчивой при Т ~ 0.

       Показано, что при s = 1+ 2k = 1.8 (k = 0.4) в СУС появляется вторая ветвь решений - (,)2, возникающая в точке = -1,  = 1 при Т= 0. В области 1.8 < s 2 с ростом s ветви решений (,) сближаются, область сосуществования фаз (,)1 и (,)2 расширяется, а разница между ТП этих фаз уменьшается. Однако фаза (,)2 так и остаётся метастабильной, а последовательность упорядочений остаётся той же, что была раньше.

Рисунок 8 - Зависимости от температуры в отсутствие ФП податливости   относительно деформации по оси z и избыточной теплоёмкости Ср

                       

Рисунок 9 Ц  Температурная зависимость восприимчивости 

при k = 0.4, s = 1.15

       При s > ~ 2 обе ветви решений смыкаются и уже при s = 2.002 расположение ветвей коренным образом меняется. Теперь каждая из них ответвляется от ветви, описывающей симметричную фазу, и заканчивается в области  Т ~ 0 при ~ 1, одна в точке = -1 - первая ветвь, другая в точке > 0 - вторая ветвь. При s > 2 реализуется следующая последовательность упорядочений: при высоких Т устойчива фаза (,0), затем происходит ФП 2 в фазу (,). Эта упорядоченная фаза остаётся устойчивой и при Т ~ 0, где ~ -1,  ~ 1. Именно такой сценарий обычно реализуется в простых многоминимумных моделях.

       В подразделе 4.2 исследованы свойства 1 + 4-мин. модели [A12]. Рассматривается модель, состоящая из двух подсистем минимумов - центрального и четырёх одинаковых боковых (1+4-мин.). Расположение минимумов приведено на рис. 10.

                                       

Рисунок 10 - Расположение и нумерация минимумов в 1 + 4-мин. модели.

Варианты рельефа потенциальной энергии частицы на сечениях по линиям у = 0 и х = 0 для  такой 1+4-мин. модели будут иметь тот же вид, изображён на рис. 7. В качестве объекта исследования рассматриваем слоистый кристалл с простой тетрагональной решёткой  симметрии D4h, в узлах которого расположены частицы, каждая из которых имеет пять минимумов в плоскостях ортогональных оси 4-го порядка (рис.11). Центральный минимум находится в самом узле, а четыре боковых - немного смещены.

Неравновесный ТП такой пятиминимумной (1+4-мин.) модели в приближении ГБВ имеет вид: (8)        

где , В - константы, характеризующие те же энергии и взаимодействия, что и в модели 1 + 2-мин. Параметры А и D характеризуют взаимодействие между частицами, находящимися в боковых ямах. D - между подсистемами минимумов, одна из которых состоит из минимумов 1 и 2, а другая - из 3 и 4. Если А < 0, то частицы, стремятся расположиться в боковых ямах одного типа (с одним номером), а если А > 0, то - в разных. Функции ,, выполняют роль параметров порядка (ПП),  преобразуется по неприводимому представлению A1g , - по В1g, 1,2 - по E1u. Система уравнений состояния имеет пять типов решений: 1.(, =1= 2 = 0) - D4h; 2. ( ,, 1= 2 = 0) - D2h;  3. (,, 1,0) - С2v;  4.( ,0, , ) - С2v; 5. (,, 1, 2) - Сs. 

       В случае одного и двух ПП особенности термодинамического поведения кристаллов, описываемых моделью 1 + 4-мин. аналогичны особенностям термодинамического поведения кристаллов, описываемых 1 + 2-мин. моделью.

Показано, что в случае трёх ПП в области >0, B >0, D >0, A< 0, где k = B/, s = -A/, с = - В/D тип упорядочения в пространстве ПП , определяется исключительно соотношением между константами A, D и температурой. Однако появление этих ПП зависит от всех модельных и термодинамических параметров. Необходимым условием возникновения ПП является А < 0, но это условие не достаточно, в случае k = 0.4 этот ПП появляется при s s* = 7.475. При с = 0.95 решение (,,1,0), возникшее из симметричной фазы, при s*   s s1 (s1 10.5) существует в ограниченном сверху и снизу интервале температур. И сверху, и снизу соседом является та же самая фаза. При s > s0 (s0 = 9.443) в области Т ~ 0, с параметрами ~ - 1, ~ 1, 1 ~ 2 появляется вторая ветвь решения (,,1,0)2. В области Т ~ 0 при s s1 (s1 = 10.5) термодинамически устойчива фаза (,0,0,0), а при s > s1 термодинамически устойчивой становится вторая ветвь фазы (,,1,0)2. 

Для решений СУС типа (,0,,) вторая ветвь появляется при s > 8(1+2k). При s > 16 фаза, отвечающая второй ветви решения (,0,,)2, становится термодинамически устойчивой при Т ~ 0.

При понижении температуры реализуются следующие последовательности упорядочений:

для k = 0.4, s = 8, с = 0.95 -  (,0,0,0) - (,,1,0) - (,0,0,0);

для k = 0.4, s = 10.4, с = 0.95 -  (,0,0,0) - (,,1,0)1 - (,,1,0)2 - (,0,0,0);

для k = 0.4, s = 12, с = 0.95 -  (,0,0,0) - (,,1,0)1;

для k = 0.4, s = 9.6, с = 1.1 -  (,0,0,0) - (,0,,) - (,,1,2) - (,,1,0) - (,,1,2)  - (,0,,) - (,0,0,0).

В первой, второй и  четвёртой последовательностях имеет место возврат в фазу, причём в четвёртой последовательный  возврат в фазу, при котором фактически возникает целый каскад ФП.

В пятом разделе исследуются ротационные фазовые переходы, один из типов структурных превращений в кристаллах с октаэдрической анионной подрешёткой, в частности в кристаллах семейства перовскита, упорядоченных перовскитах (включая эльпасолиты и  криолиты), слоистых перовскитах и в некоторых других соединениях [A5, A10]. Отмечено, что при таких переходах меняются координация ионов, межионные расстояния и т.п. Эти изменения определяются характером ротационного упорядочения и могут существенно повлиять на физические характеристики соединения.

В подразделе 5.1 исследуются ротационные фазовые переходы в кубических перовскитах и родственных им кристаллах. Вся совокупность ротационных искажений в перовскитах описывается двумя трехкомпонентными параметрами порядка, один из которых - относится к трёхлучевой звезде вектора k = (1/2)(b1 + b2), другой - принадлежит однолучевой звезде вектора k = (1/2)(b1 + b2 + b3), где bi - базисные векторы обратной решетки. Параметры и преобразуются соответственно по неприводимым представлениям М3 и R9.

Одна из основных целей этого подраздела выяснение на основе феноменологической теории условий, при которых могли бы быть реализованы указанные ниже наиболее распространённые последовательности ротационных упорядочений [9], а также исследование факторов, влияющих на формирование ФС.

1. 00000012 ;  1a. 000000;

2. 00000;                 3. 000000 1,2,3 + 3;        

4. 0000;                 5. 000000;

6. 000;                 7. 000.                                 (9)

Модельный термодинамический потенциал для описания всех возможных ротационных ФП, обусловленных конденсацией ПП и , имеет вид:

                       (,) = М +  R + MR;                               (10)

М = 1MG1M + 2MG21M +  1MG2M + MG1MG2M +  MG3M + 2MG22M + MG4M Е

R  = 1RG1R + 2RG21R +  1RG2R + RG1RG2R +  R G3R + 2RG22R Е;

MR = 1G1MG1R + (2 - 1) + 1 + 2G1MG2R +

               + 3 + 4G1RG2M + Е;

       G1M = , G2M = , G3M = ,

G4M =  123 [14 (22 - 32) + 24 (32 - 12) + 34 (12 - 22)], M3,

       G1R = ,  G2R = , G3R = , R25

везде берется либо i<j, i =1,2, j =2,3, либо i j, i k, j<k, i =1, 2, 3 j =1,2, k=2, 3. 

Он содержит члены, характеризующие вклад каждого ПП - М (), R(), а также смешанные инварианты, описывающие взаимодействие между этими параметрами порядка - МR (,). Диаграммы ФС при различных соотношениях между коэффициентами ТП (10) приведены на рис. 11,12. Диаграммы ФС получены при условии: 1M >> 0, 1<0, 2>0, 2R > 0, 2M > 0, R< 0, = 1 + 2 < 0, R < 0, M < 0. На них показаны термодинамические пути, на которых реализуются последовательности фаз из набора (9) под тем же номером. На пути 1 реализуется последовательность фаз, наблюдаемых в CsSrCl3 и CsPbCl3, последовательность фаз 1а, наблюдается в RbCdCl3 и NaTaO3 (рис. 12а); на пути 2 (рис. 11, а и 12, а и б ) - в KCdF3, CsPbBr3; на пути 3  - в NaNbO3; на пути 4 (рис. 12, б) - в KСaF3; на пути 5 (рис. 11б, 12а,б) - в SrZrO3*; на пути 6 (рис. 11, б) - в SmAlO3; на пути 7 (рис. 12, б) - в CaTiO3. Условия реализации последовательностей упорядочений 1, 4, 5 из (9) были получены ранее в работе [14]. Эти результаты согласуются с результатами настоящей работы.

 

Рисунок 11 - Диаграммы фазовых состояний ТП(10) на 1,1 - плоскости:

1М 1R, Г > 0, R< 0 (а);   1R1M, Г > 0, R< 0 (б), где R = 42R 2R - R2,

Г = 42M2R - .  Линии 1, 2, 5, 6 - термодинамические пути, на которых реализуются последовательности фаз из набора (9) под тем же номером

Из проведённого анализа было получено, что ключевым для формирования последовательности фаз является: во-первых, последовательность конденсации ПП, определяемых соотношением между 1М и 1R; во-вторых, расположение термодинамического пути относительно мультикритической N-фазной точки 1= 0, 1 = 0, определяемое взаимодействием различных компонент каждого ПП, в-третьих, характером взаимодействия параметров порядка между собой.

Аналогичным образом, с учётом упомянутых ключевых моментов, получены диаграммы фазовых состояний для кристаллов со структурой упорядоченных перовскитов (включая эльпасолиты и криолиты) (пр.гр. Oh5) и слоистых перовскитов (пр. гр. D4h1) [9].

Модельный ТП, описывающий ротационные упорядочения в кристаллах со структурой упорядоченных перовскитов (включая эльпасолиты и  криолиты) (пр. гр Oh5), имеет тот же вид, что и ТП (10), в который помимо инварианта G4 добавлен инвариант такого же вида, но относящийся к другому ПП. Как и в случае кубических перовскитов, наличие таких инвариантов 9-ой степени практически не влияет на вид диаграмм фазовых состояний, но может изменить набор ПП, характеризующих ФС. Поэтому диаграммы ротационно-упорядоченных ФС в кристаллах со структурой упорядоченных перовскитов (пр. гр Oh5) будут иметь тот же вид, что и диаграммы рис. 11, 12, однако соответствующий этим кристаллам набор ПП, характеризующих ФС, может отличаться. 

Рисунок 12 - Диаграммы фазовых состояний ТП (21) на 1М -  1R плоскости: 1R> ~ 0, R < 0, Г > 0 (а);  1R> ~ 0, R < 0, Г< 0 (б)

Ротационно-упорядоченные ФС в слоистых перовскитах описываются большим числом ПП, но меньшей размерности. Поэтому исследование ротационных ФП в этих кристаллах имеет свою специфику. Несмотря на это, использование подхода, который применялся при анализе ротационных ФП в кубических перовскитах, позволило получить диаграммы ФС и определить условия, при которых возникают наиболее распространённые последовательности ротационных упорядочений в кристаллах со структурой слоистых перовскитов [А13].

Из сопоставления характера ротационных упорядочений в кубических перовскитах с величиной фактора толерантности t и зарядом катиона А - ZA получено, что  уменьшение t при ZA= const 0 стимулирует мягкость системы по ПП в большей степени, чем по ПП , аналогично, уменьшение ZA (ZA 0) при t=const  также стимулирует мягкость преимущественно по ПП [A5].

В диссетртации предложена термодинамическая модель [A6, A14], способная на основе феноменологической теории описать всю совокупность ротационных и сегнетоэлектрических ФП, подобных тем, что наблюдаются в твёрдом растворе Na1-xKxNbO3 (NKN) при Т > 300 C, 0.03 < x < 0.2 [15]. Её основные параметры получены из анализа данных по структуре и температурных зависимостей диэлектрической проницаемости для твёрдых растворов NKN различного состава. Это дало возможность определить расположение фазовых состояний на x-T плоскости и на качественном уровне воспроизвести x-T диаграмму фазовых состояний, наблюдаемую в NKN при Т>300 C, 0.03< x< 0.2.

На базе этой модели проанализированы особенности диэлектрических свойств и исследовано влияние внешних воздействий - давления и электрического поля на ФС [A6, A15]. Показано, что вблизи перехода в полярную фазу аномалии диэлектрических свойств, обусловленные ФП между параэлектрическими фазами (ротационные ФП), должны проявляться значительно сильнее чем те, что наблюдаются вдали от него. Это хорошо заметно из сравнения аномалий (изломов) на температурных зависимостях диэлектрической проницаемости при х ~ 0,04÷0,06 и при x ~ 0.1. 

Отмечено, что особенности структуры перовскита позволяют посредством изовалентных и гетеровалентных замещений компонент соединений существенным образом влиять на степень их ротационной и поляризационной неустойчивости [А16 - А18].

Показано, что формирование упорядоченного ФС в MnF3 обусловлено совместной конденсацией ротационного ПП и орбитального ПП, характеризующего ориентационное упорядочение октаэдров MnF6, искажённых вследствие эффекта Яна-Теллера. Вклад от орбитального ПП описывается на основе 3-мин. двухподрешёточной модели, вклад от ротационного - на основе феноменологической теории и при этом учитывается нелинейное взаимодействие между обоими параметрами порядка [A10, A19]. Это дало возможность определить условия, при которых фазовый переход из симметричной фазы MnF3 становится фазовым переходом триггерного типа, при котором появление одного параметра порядка стимулирует появление другого.

                               ЗАКЛЮЧЕНИЕ

       Основные результаты и выводы:

1. Диаграммы фазовых состояний 4-мин. квадратной, 6- и 8-мин. кубических моделей содержат сегнетоэластическую и три сегнетоэлектрические (полярные) фазы. На границе между сегнетоэластической и одной из полярных фаз, а также на границе между симметричной и этой же полярной фазами имеются трикритические точки, в которых фазовые переходы I рода становятся фазовыми переходами II рода. На линии фазовых переходов II рода из симметричной фазы расположена мультикритическая трёхфазная точка, в которой с симметричной граничат две полярные фазы. В области низких температур на диаграммах фазовых состояний 4-, 6- и 8-мин. моделей имеются по три полярные фазы. В случае 6- и 8-мин. моделей, между  ромбоэдрической и тетрагональной фазами расположена либо моноклинная (6-мин.), либо ромбическая (8-мин.) фаза. В случае 4-мин. модели между двумя ромбическими фазами разного типа расположена моноклинная фаза.

2. Учёт одноосных деформаций в 8-мин. кубической модели, приводит к появлению на границе с симметричной фазой ещё одной трикритической точки, в которой фазовые переходы II рода в тетрагональную полярную фазу  становятся фазовыми переходами I рода. Установлено, что последовательность ФП, наблюдаемых в BaTiO3 и KNbO3 может быть описана на основе этой модели. 

3. Учёт сдвиговых деформаций в 6-мин. кубической модели, даже при относительно слабом стрикционном взаимодействии, приводит к тому, что вместо моноклинной СЭ фазы на диаграмме фазовых состояний появляется ромбическая СЭ фаза.

4. Совпадение температур Кюри - Вейсса кристаллов ANbO3 (A = Na, K, Ag), экстраполированных из данных по температурной зависимости диэлектрической проницаемости в кубической фазе, объяснено на основе восьмиминимуной модели для катионов Nb.

5. Получены диаграммы фазовых состояний 12-мин. кубической модели. Они содержат либо три, либо четыре полярные фазы. На границе с симметричной имеются: две трикритические точки, в которых фазовый переход I рода превращается в фазовый переход II рода; мультикритическая трёхфазная точка, в которой сходятся две линии фазовых переходов II рода и линия фазовых переходов I рода между полярными фазами; трёхфазная точка, в которой сходятся три линии фазовых переходов I рода. В области низких температур существуют три полярные фазы: либо ромбоэдрическая и две, отличающиеся по типу, ромбические, либо ромбоэдрическая, тетрагональная и ромбическая.

6. Предложен и применён при исследования статистических свойств многоминимумных, многоподрешёточных моделей метод, включающий следующие этапы:

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

б. Исследование симметрийных и статистических свойств  рассматриваемой многоминимумной модели одной из подрешёток. 

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

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

7. В приближении Горского - Брэгга - Вильямса получены базовые диаграммы фазовых состояний для 3-мин. восьмиподрешёточной кубической (пр. гр. Oh1) модели, описывающие ориентационные упорядочения антиферродисторсионного типа, характеризуемые параметрами порядка, относящимися к представлениям R3+ и M2+.

8. Диаграмма фазовых состояний 3-мин. двухподрешёточной кубической модели содержит одну сегнетоэластическую  фазу (фаза I, D4h1 (Z= 1)) и четыре фазы с антиферродисторсионным типом упорядочения: фазы II и II1 с симметрией  D4h18 (Z = 4), фазу III (D4h17 (Z = 4)) и фазу IV (D2h23 (Z = 4)). Фаза I граничит с симметричной (С) и фазой II по линии фазовых переходов I рода. На границе с симметричной (Oh) и фазой  II имеется трикритическая точка, в которой фазовый переход I рода становятся фазовым переходом II рода. На линии фазовых переходов II рода из симметричной фазы имеются  две мультикритические четырёхфазные точки, в которых с фазой С граничат фазы II и III, либо III и II1, а также ромбическая фаза IV. В области низких температур, границы между фазами II - IV,  IV - III и III - IV,  IV Ц  II1 асимптотически сходятся к двум разным линиям, подобно тому, как они сходились к критическим точкам на границе с симметричной фазой. Поскольку фазы II и II1 имеют одинаковую симметрию и существуют при разных  термодинамических условиях, то может иметь место возврат в фазу, сопровождающийся целым каскадом ФП: С - II - IV - III - IV - II1.

9. Показано, что в кристаллах с несколькими наборами кристаллографически эквивалентных позиций появляются состояния, в которых при Т ~ 0 оказываются частично заполнены позиции с более высокой энергией, то есть  возбуждённые состояния. 

10. Установлен и описан в рамках статистической теории один из возможных механизмов изоструктурных (изосимметричных) фазовых переходов для кристаллов с несколькими наборами кристаллографически эквивалентных позиций. 

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

12. Влияние кристаллохимических параметров компонент соединений со структурой перовскита (пр. гр. Oh1) таких, как заряд катиона А (ZA) и соотношение ионных радиусов, характеризуемое фактором толерантности t, на формирование фазовых состояний при ротационных фазовых переходах проявляется следующим образом: при t < ~ 1 температура конденсации моды R25 выше чем моды М3. При фиксированной величине ZA уменьшение t сопровождается более интенсивным ростом  температуры конденсации моды М3 по сравнению с температурой конденсации моды R25. При t (t~0,9, ZА=1) моды М3 и R25 конденсируются одновременно. При t ~ 0,87 сначала конденсируется М3, затем - R25. С ростом ZA при фиксированном t  в большей степени стимулируется конденсация моды R25. Поэтому для кристаллов с ZA  = 2 при t = 0.88 сначала конденсируется мода R25, затем М3. Совместная конденсация мод М3 и R25 имеет место только при t ~ 0.85. Вариант, при котором сначала конденсируется мода М3, представляется маловероятным, так как при t < 0.85 проблематично само существование структуры перовскита.

13. Нелинейные взаимодействия ротационных степеней свободы различного типа с деформационными в потенциале Ландау приводят к появлению промежуточных псевдособственных деформационных фазовых переходов в кубических перовскитах: RbCdCl3, NaTaO3, KСaF3, SrZrO3; в упорядоченных перовскитах: Pb2LuSbO6, Pb2CoTeO6, K2SeCl6, Rb2KMF6 (M = Sc, In, Lu), Cs2RbDyF6; а также в слоистых перовскитах: K1-xRbxAlF4 и TlAlF4.

14. Предложен модельный термодинамический потенциал, учитывающий вклады от параметров порядка, характеризующих ротационное и  ориентационное упорядочение октаэдров MnF6, искажённых вследствие эффекта Яна-Теллера, а также  нелинейное взаимодействие между обоими параметрами порядка. На его основе определены  условия, при которых фазовый переход из симметричной фазы MnF3 становится фазовым переходом триггерного типа, при котором появление одного параметра порядка стимулирует появление другого.

Список цитируемой литературы

1. Парсонидж Н., Стейвли Л. Беспорядок в кристаллах. Т.1.-М.: Мир, 1982.-434с.

  2. Френкель Я.И. Кинетическая теория жидкости.- Л.: Наука, 1975. - 592с.

  3. Mason W.P., Mattias В.Т. Theoretical model explaining the ferroelectric effect in Barium titanate // Phys.Rev.-1948.- V.74.- P.1622 - 1636.

  4. Comes R., Lambert M., Guinier A. The chain structure of BaTiO3and KNbO3 // Sol. St. Commun.-1968.-V.6.- P.715 - 719.

  5. Вихнин В.С., Зайцев О.А. Фазовые переходы и динамические эффекты в кристаллах, обладающих одноячеечными потенциалами с многоямным возбужденным состоянием // ФТТ.- 1997.-Т.39.- С.547- 552.

6. Udovenko A.A., Laptash N. M. Disorder in crystals of dioxofluorotungstates, (NH4)2 and Rb2WO2F4 //Acta Cryst. B.- 2008.-V. 64.-P. 645 - 651. 

7. Блат Д. Х., Зиненко В. И. Теория структурных фазовых переходов в цианидах щелочных металлов // ЖЭТФ.- 1980.- Т.79.- С. 974 - 986.

8. Ландау Л.Д., Лифшиц Е.М. Статистическая физика. ч.1.- М.: Наука, 1976.-584с.

9. Александров К.С., Безносиков Б.В. Перовскиты. Настоящее и будущее.- Новосибирск: Изд-во СО РАН,  2004.- 231с.

10. Pirc R., Blinc R. Off-center Ti model of barium titanate // Phys. Rev. B.-  2004.- V.70.- 134107-1 - 134107-8. 

11. Lemeshko M.P., Nazarenko E.S., Gonchar A. A., Reznichenko L. A., Nedoseykina N. I., Novakovich A.A., Mathon O., Joly Y., Vedrinskii R.V. EXAFS studies of  the local atomic structure of lead- free piezoelectric ceramics KxNa1-x NbO3 over  temperature range 10 - 1023 K // Phys. Rev. B.- 2007.- V.76.- 134106-1 - 134106-11.

12. Sicron N., Ravel B., Yacoby Y., Stern E. A., Dogan F., Rehr J. J. Nature of ferroelectric phase transition in PbTiO3 // Phys. Rev.B.- 1994.-V.50.- P. 13168 - 13180. 

  13. Сахненко В.П., Таланов В.М., Чечин Г.М. Теоретико-групповой анализ полного конденсата, возникающего при структурных фанзовых переходах // Физика металлов и металловедение. -1986.- Т.62.-  С. 847 - 856.

  14. Ларин Е.С. О фазовой диаграмме CsPbCl3 // ФТТ.- 1984.-Т.26.- С.3019 - 3023.

  15. Ahtee M., Glazer A.M. Lattice parametrs and tilted octahedra in sodium - potassium niobate solid solutions //Acta. Cryst. A.- 1976.-V. 32.-P. 434 - 446.

Список основных публикаций автора по теме диссертации

А1. Ивлиев, М.П. Фазовые переходы типа порядок - беспорядок, описываемые двумя многоминимумными параметрами порядка (приближение среднего поля)./ М.П. Ивлиев, В.П. Сахненко // Деп. ВИНИТИ № 3418 - 89 от 19.05.89 - Ростовск. гос. ун-т, 1989. - 35с.

  А2. Ивлиев, М.П. Формирование сегнетоэлектрических фазовых состояний в KNbO3 и других ниобатах со структурой перовскитов/ М.П. Ивлиев, И.П. Раевский, С.И. Раевская, В.А. Шуваева, И.Н. Пирог // ФТТ.- 2007.- Т.49.- С. 731-740.

А3. Ивлиев, М.П. Сегнетоэластические фазовые переходы в перовскитах, обусловленные орбитальным и ротационным упорядочением / М.П. Ивлиев // ФТТ.- 2009.- С. 1472-1476.

А4. Ивлиев, М.П. Влияние сегнетоэластических параметров порядка на формирование фазовых состояний титаната - цирконата свинца / М.П. Ивлиев // Изв. РАН. сер. физ.- 2010.-  Т.74.- С. 1257 - 1259.

А5. Ивлиев, М.П. Фазовые состояния в ротационно-искажённых перовскитах / М.П. Ивлиев // Кристаллография.- 2002.- Т.47.- С.1065-1071.

А6. Ивлиев, М.П. Фазовые состояния и особенности диэлектрических свойств твёрдых растворов ниобата натрия-калия / М.П. Ивлиев, И.П. Раевский, Л.А. Резниченко, С.И. Раевская, В.П. Сахненко // ФТТ.- 2003.- Т.45.- С. 1886-1891.

А7. Ивлиев, М.П. Возможность существования квадрупольной сегнетоэластической фазы в KTaO3 : Li / М.П. Ивлиев, В.П. Сахненко // ФТТ.- 1986.- Т.28.- С. 632-634.

А8. Ивлиев, М.П. Особенности сегнетоэлектрических упорядочений в перовскитах АNbO3 (А - Na, Ag, K) / М.П. Ивлиев, С.И. Раевская, И.П. Раевский // Труды VIII Междунар. Симпоз. Порядок, беспорядок и свойства оксидов (ODPO - 2005). Сочи, 2005. Ч.2.- С.150 - 153.

А9. Ивлиев, М.П. Теория сегнетоэластических фазовых переходов в KCN / М.П. Ивлиев, В.П. Сахненко // Изв. АНСССР. сер. физ.- 1979. -Т.43. - С.1606-1610.

  А10. Ивлиев, М.П. Фазовые состояния трёхминимумной двухподрешёточной модели (приближение среднего поля) / М.П. Ивлиев // Труды Междунар. симпоз. Фазовые превращения в твёрдых растворах и сплавах (OMA - 2003). Сочи, 2003. - С.133-136.

А11. Ивлиев, М.П. Особенности термодинамических свойств кристаллов, описываемых составными многминимумными моделями / М.П. Ивлиев, В.П. Сахненко // Известия ВУЗов. Сев.Кавказский регион.- 2010.- Т.23.- С.41-45.

А12. Ивлиев, М.П. Complex Aggregative Multiminima Models for Ferroelastic Transitions in Crystals / M.P. Ivliev, V.P. Sakhnenko // Ferroelectrics.-2010.-Р.16-21.

А13. Ивлиев, М.П. Soft modes and martensitic transformations in K1-xRbxAlF4 / M.P. Ivliev, V.P. Sakhnenko // Ferroelectrics.- 1996.- V.175.- P.65-71.

А14. Ивлиев, М.П. Получение и исследование монокристаллических твёрдых растворов (Na,K)NbO3 / И.П. Раевский, Л.А. Резниченко, М.П. Ивлиев, В.Г. Смотраков, В.В. Ерёмкин, М.А. Малицкая, Л.А. Шилкина, С.И. Шевцова, А.В. Бородин // Кристаллография.- 2003.- Т.48.- С.531-535.

А15. Ивлиев, М.П. Кристаллохимические аспекты влияния термодинамической предыстории на вид фазовых диаграмм температура - состав твёрдых растворов ниобата натрия - лития и ниобата натрия - калия / М.П. Ивлиев, И.П. Раевский, Л.А. Резниченко, М.Н. Палатников, Л.Е. Балюнис, М.А. Малицкая // ЖТФ.- 2002.- Т. 72.- С. 120 - 124.

А16. Ивлиев, М.П. Влияние модифицирования литием на устойчивость сегнетоэлектрического состояния в пьезоэлектрических, керамических материалах / М.П. Ивлиев, С.И. Раевская,  О.Ю. Кравченко, Л.А. Резниченко, И.П. Раевский//Конструкции из композитных материалов.-2008.-Вып.4.- С.80-87.

  А17. Ивлиев, М.П. Влияние несоответствия кристаллохимических параметров катионов Na и Li на диэлектрические свойства твёрдых растворов NaNbO3 - LiNbO3 / М.П. Ивлиев, С.И. Раевская,  О.Ю. Кравченко, Л.А. Резниченко, И.П. Раевский // ФТТ.- 2009.- С. 1988-1993.

  А18. Ивлиев, М.П. Влияние модифицирования изовалентными и гетеровалентными ионами на диэлектрические свойства пьезоэлектрических керамических материалов на основе ниобата натрия / М.П. Ивлиев, С.И. Раевская,  О.Ю. Кравченко, Л.А. Резниченко, И.П. Раевский // Конструкции из композитных материалов.- 2009.- Вып.1.- С.61-68. 

А19. Ивлиев, М.П. Взаимодействие ротационного и орбитального упорядочений при формировании фазового состояния в MnF3 / М.П. Ивлиев // Труды VII Междунар. Симпоз. Порядок, беспорядок и свойства оксидов (ODPO - 2004). Сочи, 2004.- С.268 - 269.

  Авторефераты по всем темам  >>  Авторефераты по физике