Атомистическое компьютерное моделирование кристаллической структуры и свойств минералов, их дефектов и твердых растворов

Еремин, Николай Николаевич

Бесплатный автореферат и диссертация по наукам о земле на тему
Атомистическое компьютерное моделирование кристаллической структуры и свойств минералов, их дефектов и твердых растворов
ВАК РФ 25.00.05, Минералогия, кристаллография

Автореферат диссертации по теме "Атомистическое компьютерное моделирование кристаллической структуры и свойств минералов, их дефектов и твердых растворов"

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

ЕРЕМИН Николай Николаевич

АТОМИСТИЧЕСКОЕ КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ КРИСТАЛЛИЧЕСКОЙ СТРУКТУРЫ И СВОЙСТВ МИНЕРАЛОВ, ИХ ДЕФЕКТОВ И ТВЕРДЫХ РАСТВОРОВ

25.00.05 - минералогия, кристаллография

Автореферат диссертации на соискание ученой степени доктора химических наук

□□3478225

Москва - 2009

003478225

Работа выполнена на кафедре кристаллографии и кристаллохимии Геологического факультета Московского государственного университета

имени М.В. Ломоносова.

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

Урусов Вадим Сергеевич

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

наук, профессор Гусаров Виктор Владимирович

член-корреспондент РАН, доктор химических наук, профессор Кусков Олег Львович

доктор физико-математических наук, профессор Лебедев Александр Иванович

Ведущая организация: Институт геологии и геохимии УрО РАН

Защита состоится « _ часов в аудитории 415

на заседании диссертационного совета Д.501.002.06 при Московском государственном университете имени М.В. Ломоносова по адресу: 119991, ГСП-1, Москва, Ленинские горы, МГУ, Геологический факультет.

С диссертацией можно ознакомиться в библиотеке Геологического факультета МГУ имени М.В. Ломоносова (зона А, 6 этаж).

Автореферат разослан «.

:£УС» 2009 г.

Ученый секретарь , Киселева И. А.

диссертационного совета '"Ус/¿¿¿¿м- у*

доктор геолого-минералогических наук

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

Актуальность проблемы. Два последних десятилетия прошлого века и первое - нынешнего ознаменовались гигантским прогрессом в области непосредственного теоретического предвидения структуры и свойств кристаллических веществ методами, основанными на поиске минимума структурной энергии кристалла. Такая задача решается, в принципе, строго с помощью квантовохимических расчетов и приближенно - с помощью атомистических полуэмпирических методов моделирования. На сегодняшний день, использование полуэмпирических подходов способно принести успех там, где квантовохимические методы пока не в состоянии решить задачу поиска устойчивой атомной конфигурации. Метод атомистического моделирования, ориентируясь на понятийный аппарат кристаллохимии, не требует рассмотрения внутриатомных взаимодействий, что сразу существенно облегчает вычислительную задачу поиска минимума структурной энергии системы. Наиболее заманчивые возможности метода заключаются в предсказании структуры и свойств еще неизвестных, но прогнозируемых веществ с необычными, особенно важными для технологии свойствами. В области наук о Земле предсказательный потенциал методов компьютерного моделирования играет большую роль при воссоздании реалистичной картины минерального строения и физических свойств глубинных недр планет, которые не могут наблюдаться прямыми методами. Однако, в рамках атомистического подхода существуют определенные трудности: недостаточная транеферабельность (переносимость) наборов межатомных потенциалов от одного класса соединений к другим, отсутствие строгих энергетических критериев моделирования с неиелочислегшыми зарядами атомов, недостаточная разработка методик успешного моделирования реальных кристаллов, содержащих различные структурные несовершенства, вплоть до протяженных изоморфных смесей.

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

Иа защиту выносятся:

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

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

3. Систематические расчеты энергий и позиций собственных точечных дефектов в минералах и неорганических соединениях разных классов.

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

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

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

Практическая значимость работы;

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

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

3) Разработанный алгоритм нахождения максимально неупорядоченной атомной конфигурации произвольного состава для ячейки любого размера может быть рекомендован как для ab-initio расчетов, так и для моделирования полуэмпирическими методами. Кроме того, этот подход может быть использован при интерпретации сложных экспериментальных спектров изоморфных смесей.

4) Разработанные программные продукты адаптированы для пользователей, снабжены руководством пользователя с подробным описанием, инструкцией по эксплуатации и примерами файлов информации, которые находятся в открытом доступе по адресам: http://ciyst.geol.msu.ru/odss/ и http://ciyst.geol.msu.ru/values/. Программы востребованы для расчета научными группами, как в нашей стране, так и за рубежом.

5) Материалы диссертации явились основой для создания нового спецкурса «Современные методы моделирования твердых тел» читаемого автором на геологическом факультете МГУ начиная с 2007 года, некоторые материалы настоящей работы используются в курсе «Компьютерное программное обеспечение в кристаллохимии», читаемого для студентов 3-го курса геохимического потока геологического факультета.

Личный вклад автора: Все результаты, представленные в диссертации, получены самим автором, либо при его непосредственном участии, либо под его руководством. Создание программного продукта VALUES осуществлялось автором совместно со студентом Оболенским О.В. Создание пакета компьютерных программ ODSS, разработка алгоритмов для расчета осуществлялось автором совместно с Урусовым B.C. (научные консультации) и Деяновым Р.З. (программирование). Компьютерное моделирование структур оксидов щелочноземельных металлов и их твердых растворов осуществлялось совместно с Петровой Т.Г., длительное время работавшей сначала под руководством автора, а затем и в качестве соавтора. Часть результатов, представленных в работе получена совместно со студентами и аспирантами кафедры кристаллографии Леоненко Е.В. (системы NaCl-KCl и Ti02-Sn02) и Талисом P.A. (системы на основе корунда), защитившими под руководством автора свои магистерские и бакалаврские работы. Изоморфные системы на основе хризоберилла исследовались автором совместно с аспиранткой кафедры кристаллографии Громаловой H.A.

Работа выполнена на кафедре кристаллографии и кристаллохимии Геологического факультета МГУ им. М.В.Ломоносова. Отдельные этапы работы были поддержаны грантами INTAS № 94-2826, 97-32174, РФФИ №№ 96-05-64567, 99-05-65139-а, 00-15-98582-л, 01-05-06164, 02-05-64845-а, 05-05-64721-а, 09-05-00403-а, Ведущих научных школ Российской федерации №№ НШ-1955.2003.5, НШ-8091.2006.5, НШ-1880.2008.5, ФЦП «Интеграция» подпроект №683 проекта АО 109 «Науки о Земле».

Апробация работы: Материалы, включенные в диссертацию, докладывались на следующих всероссийских и международных конференциях: 16 конгресс Международного союза кристаллографов, Пекин, 1993; 5 Международное совещание по ядерно-спектроскопическим исследованиям сверхтонких взаимодействий. Дубна, 1993; 11 Мегвдународное кристаллографическое совещание «Sagamore», 1994; 2 Европейское совещание по резонансной абсорбции и спектроскопии в минералогии, Берлин, 1995; Международное совещание АСА, Монреаль, 1995; 3 Европейское совещание по спектроскопическим методам в минералогии, Киев, 1996; Международная конференция «Спектроскопия, рентгенография и кристаллохимия минералов», Казань, 1997; 1 Национальная кристаллохимическая конференция, Черноголовка, 1998; 10 Международное совещание кристаллографов, Страсбург, 1999; 18 Конгресс международного союза кристаллографов, Глазго, 1999; 19 Европейское совещание кристаллографов, Нанси, 2000; 19 Чтения

имени академика Белова, Нижний Новгород, 2000; 2 Национальная конференция по кристаллохимии, Черноголовка, 2000; Международная конференция «Ломоносовские чтения», Москва, 2000, 2007, 2009; Конференция французского общества минералогов и кристаллографов, 2001; 15 Международное совещание «Рентгенография и кристаллохимия минералов», Санкт-Петербург, 2003; 16 Международная конференция по химии фосфора, Бирмингем, 2004; Всероссийская конференция «Химия твердого тела и функциональные материалы-2004», Екатеринбург, 2004; Международная конференция «Спектроскопия, рентгенография и кристаллохимия минералов», Казань, 2005; 20 Конгресс международного союза кристаллографов, Флоренция, 2005; 4 Национальная кристаллохимическая конференция, Черноголовка, 2006; Международная конференция «Спектроскопия и кристаллохимия минералов», Екатеринбург, 2007; 16 Международное совещание «Кристаллохимия и рентгенография минералов», Миасс, 2007; 27 Научные чтения имени академика Н.В.Белова, Нижний Новгород, 2008; 4 Всероссийская конференция «Математические исследования в кристаллографии, минералогии и петрографии», Апатиты, 2008; Федоровская сессия РМО, Санкт-Петербург, 2008; Всероссийская конференция «Минералы: строение, свойства, методы исследования», Миасс, 2009.

Публикации. Результаты исследований, изложенные в диссертации, отражены в 82 публикациях: 31 статье в отечественных и международных журналах, 3 обзорных статьях в сборниках, 3 препринтах и 45 тезисах докладов, а также в отчетах по проектам РФФИ, ИНТАС, грантам Президента РФ «Ведущие научные школы». Основных работ: 35.

Структура и объем работы. Диссертация состоит из введения, шести глав, заключения, выводов, списка цитированной литературы и приложений. Общий объем работы - 352 страницы, включая 74 таблицы (плюс 5 таблиц в приложениях), 114 рисунков (плюс 1 рисунок в приложениях). Список литературы состоит из 444 наименований.

Благодарности. При выполнении работы автору оказывалась постоянная помощь и поддержка со стороны своего научного консультанта академика РАН Урусова B.C. с которым автору посчастливилось работать вместе уже около 20-ти лет. При проведении теоретических расчетов автору помогали его коллеги по научной группе В.С.Урусова: проф. Оганов А.Р., н.с. Петрова Т.Г. и Кучериненко Я.В., асп. Леоненко Е.В. и Громалова H.A., магистрант Талис P.A., студенты Курыкина М.М., Оболенский О.В., Горяева A.M. Автор также благодарит всех своих коллег-соавторов, внесших свой вклад в его

диссертацию: Фурманову Н.Г., Кузнецова В.А., Охрименко Т.М., свою жену и соавтора ряда статей Еремину Т.А., Якубович О.В. Дубровинского Л.С., Вассермана Е.А., Храмова Д.А., Глазкову М.А., Ефремову Е.П., Игнатьева В.Д., [Леонюк Л.И.), Жарикова Е.В., Дудникову В.Б., Суханова М.В., Веремейчик Т.Ф., Петькова В.И., Куражковскую B.C. Отдельно хотелось бы поблагодарить Деянова Р.З. за воплощение авторских идей и алгоритмов в законченный программный продукт, Вассермана Е.А. и Товбиса А.Б. за адаптацию программных продуктов METAPOCS и GULP к персональным компьютерам, а Агошкова В.М. за предоставление пакета программ SPECTRA для термодинамических расчетов. Огромную признательность хотелось бы высказать зарубежным коллегам - профессорам С. Паркеру и Дж. Гэйлу, предоставившим в распоряжение автора программные продукты METAPOCS, GULP и METADISE, а также профессору Дж. Прайсу предоставившему интернет-доступ к расчетному кластеру в Великобритании. Существенную помощь автору оказали высококвалифицированные консультации профессора Русакова B.C. (физический факультет МГУ), доцентов кафедры

кристаллографии геологического факультета МГУ рЕгорова-Тисменко Ю.К.

Дороховой Г.И. и многих других сотрудников кафедры кристаллографии и кристаллохимии геологического факультета МГУ, лаборатории геохимии твердого тела ГЕОХИ РАН и кафедры общей физики физического факультета МГУ. Всем этим лицам автор выражает свою глубокую признательность и сердечную благодарность.

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

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

Глава 1. Общая методология атомистического моделирования структур и свойств неорганических кристаллов.

Представлен обзор современного состояния вопроса теоретического моделирования кристаллических структур и свойств минералов с использованием метода атомистического моделирования, основанного на поиске минимума энергии межатомного взаимодействия в кристалле путем варьирования координат атомов и параметрами элементарной ячейки. В найденной точке энергетического минимума силы притяжения и отталкивания на оптимальном для каждой пары атомов расстоянии Яо выровнены (рис. 1а).

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

более короткодействующих взаимодействий, входит еще и АЕ - энергия переноса заряда - внутриатомный эффект, описывающий работу, необходимую для отрыва внешних электронов от катионов и для присоединения их к анионам:

Еат = Естр + ДЕ = Е"* + YUp +... + АЕ, (1)

где Uij- энергия парных, Ui]k - трехчастичных межатомных взаимодействий и так далее.

расстояния R: 1 - притяжение, 2 - потенциалы ионизации катионов к, F -отталкивание, 3 - результирующая сродство к электрону анионов /). кривая.

Из-за нестабильности в свободном виде (в газовой фазе) таких многозарядных ионов, как О2", S2", N3" и т.д. термохимическая оценка энергии решетки U, выделяющейся при образовании кристалла из газа свободных ионов, невозможна для многих минеральных классов соединений [Урусов B.C. (1975)]. Поэтому единственной экспериментальной оценкой энергии сцепления для таких веществ становится энергия атомизации Еат, выделяющаяся при образовании кристалла из газа нейтральных атомов (рис. 16). Однако,

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

В работе впервые показано, что ионизация в кристалле многовалентного атома в гибридном валентном состоянии существенно отличается от последовательной ионизации свободного изолированного атома. Для описания такого процесса ионизации атома в кристалле функция переноса заряда должна быть существенно модернизирована. Попытки такой модификации неоднократно предпринимались рядом авторов [Watson R.E. et а/ (1983), Дубровинский Л.С. (1986), Урусов B.C. и др., (1994)], однако их синтезированные кривые учета энергии переноса заряда разумно описывают процесс катионной ионизации только в некотором ограниченном интервале степеней ионности связи /.' В связи с этим данная проблема требует более строгого решения.

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

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

атомов, параметры потенциалов отталкивания, дисперсионного взаимодействия и др.) с привлечением известных из эксперимента свойств кристалла.

Проведен сравнительный анализ статического моделирования в ионном и частично ионном приближении, отмечены достоинства и недостатки каждого. Огромное число работ, относящихся к предсказанию структур и свойств минералов, основывается на ионной модели, в рамках которой кристалл состоит из дискретных ионов (катионов и анионов) с интегральными зарядами Z, , Zj и.т.д. [Burnham C.W. (1985), Post J.E. and Burnham С. IV. (1986)]. Неоспоримым преимуществом и удобством ионной модели является ее универсальность - в любом соединении (ионном, ковалентном и т.д.) каждому атому можно приписать заряд, соответствующий его формальной валентности. Никаких трудностей не возникает и при проведении более сложных вычислений, таких как расчеты дефектов, включая гетеровалентные замещения - недостаток или избыток целочисленного заряда компенсируется по одной из схем гетеровалентного изоморфизма для соблюдения принципа полного электронного баланса. Однако рассчитанные значения ряда свойств моделируемых кристаллов не очень хорошо согласуются с экспериментальными величинами. Так, значения упругих констант обычно существенно завышены и теоретические «ионные кристаллы» гораздо более «жесткие», чем следует из эксперимента.

Как хорошо известно [Урусов B.C. (1975), Catlow C.R.A. and Price G.D. (1990), Урусов B.C., Дубровинский Jl.С. (1989)], подавляющее большинство минералов, например, оксиды, оксисоли, халькогениды и др., не могут быть корректно описаны ни как чисто ионные, ни как чисто ковалентные соединения. В общем случае связи в таких кристаллах имеют промежуточный характер между ионными и ковалентными. Опыт многих исследователей [Price G.D. and Parker S.C. (1984), CattiM. (1986),MatsuiM. (1986) и др.] показывает, что с учетом эффектов ковалентности в наборах межатомных потенциалов, удается добиться очень хорошего предсказания кристаллических структур и физических свойств таких минералов, как оксиды и силикаты.

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

программного обеспечения атомистического моделирования, при этом обращено особое внимание на возможности ведущег о программного продукта в этой области - программы GULP [Gale J.D., RohlA.L. (2003)].

Глава 2. Разработка и усовершенствование методов атомистического компьютерного моделирования.

2.1. Метод минимизации энергии атомизации, основанный на поиске оптимальных зарядов на атомах.

Проблема учета энергии переноса заряда при ионизации катионов, содержащих в валентной электронной оболочке электроны разных типов, потребовала разработки определенной физически обоснованной процедуры смешения энергий отдельных валентных уровней. Такая процедура учета энергии переноса нецелочисленного заряда катиона, находящегося в гибридном валентном состоянии, при любом значении / разрабатывались при непосредственном участии автора на протяжении ряда лет [Urusov V.S. and Eremin N.N. (1995), (1996), (1997), Еремин H.H (1996), Урусов B.C., Еремин Н.Н., Огонов А.Р. (1997), Урусов B.C., Оганов А.Р., Еремин ДД (1999)]. Основным принципом процедуры является метод выравнивания орбитальных электроотрицательностей (ЕЕМ) отдельных электронных уровней одного ионизируемого атома. Под орбитальной электроотрицательностью понимается

первая производная энергии орбита™ по ее заселенности: х, = /¿q > гДе Ч, -

число электронов. Принцип выравнивания сводится к требованию: X, = ХР = Xd - Xt, что является основанием для определения заселенности s, р, d,

/-орбиталей при заданной величине Q,, т.е. для восстановления хода ионизации. На рис. 2. в качестве примера приведена синтезированная кривая ионизации атома Сг в .«/^-состоянии.

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

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

Рис. 2. Энергия переноса заряда для атома Сг в .чрс!- валентном

ЯП

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

Расчет структурных, упругих ' ■ «термодинамическихсвойств соединения в оптимальной точке

Рис. 3. Схема структурного моделирования путем минимизации энергии атомизации соединения.

2.2. Разработка методик структурного моделирования твердых растворов минералов.

Теоретические подходы к моделированию структур и свойств концентрированных твердых растворов основываются на феноменологических моделях, либо на атомистических расчетах или расчетах «ab-initio». Для учета всех локальных конфигураций, которые возникают в изоморфной смеси требуется применение моделей, принципиально отличных от методики вложенных сфер, используемой для расчетов изолированных точечных дефектов [Mott N. F., Littleton M.J. (1938)]. На рис. 4 приведена общая схема оригинальной методики моделирования концентрированных твердых растворов замещения во всем диапазоне составов твердого раствора.

1 этап 5 этап

Рис. 4. Общая схема структурного моделирования твердого раствора замещения

по предлагаемой методике.

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

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

На втором этапе выбирается сверхструктура, в пределах которой задаются конфигурации различного состава, максимально приближенные к неупорядоченному твердому раствору. Описаны проблемы, связанные с выбором неупорядоченной конфигурации, и предлагается новый способ их решения, подробно описанный в недавней работе автора [Еремин H.H., Деянов Р.З., Урусов B.C. (2008)]. В качестве критерия степени неупорядоченности конфигурации использовалась величина квадратов отклонений числа разнородных пар атомов во второй координационной сфере для случайной конфигурации от статистической теоретической гистограммы (критерий согласия Пирсона х2)- Множество случайных конфигураций анализировалось по величине отклонения от идеальной статистической гистограммы частоты встречаемости разнородных вторых соседей для каждого состава (рис. 5). Результирующая гистограмма распределения (серые столбцы на рис. 5) получается как суперпозиция двух и более вкладов от атомов различных сортов, изоморфно замещающих друг друга в одной атомной позиции. К достоинствам предложенного способа нахождения оптимальной конфигурации атомов, реализованного в программе BINAR [Еремин H.H., Деянов Р.З, Урусов B.C. (2008)] можно отнести: наличие количественной оценки качества конфигурации, возможность в рамках сравнительно небольшой сверхячейки максимально приблизиться к неупорядоченному статистическому распределению в кристалле, возможность в дальнейших расчетах обойтись одной конфигу рацией, что немаловажно для минимизации расчетного времени.

Относительное количество

пас в свеохячейке

Рис. 5. Статистические

(черные столбцы) и полученные по программе BINAR (серые столбцы) гистограммы распределения разнородных пар М'-М2 в сверхячейке 4x4x8 состава 1:2 системы tío2-s11o2.

23456789 10 Число разнородных пар М'-М*

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

а) вторые соседи должны окружать центральный атом как можно более равномерно со всех сторон (образовывать координационный полиэдр таким образом, чтобы центральный атом находился внутри него и вблизи от его «центра тяжести» с минимальной несферичностью);

б) число вторых соседей должно находиться, как правило, в интервале 10-14 (12 - для ГЦК и ГПУ и 14 - для ОЦК-упаковок). Отметим, что число 10 является минимально рекомендованной величиной для критерия Пирсона %2.

Так, для структурного типа рутила ТЮ2 при построении замкнутого полиэдра вокруг центрального атома необходимо учитывать, помимо ближайших, 8 следующих вторых соседей. Следовательно, общее число вторых соседей катиона равно 10 (2+8) (рис. 6а). В структурном типе корунда ситуация более сложная, и, как показал проведенный геометрический анализ, наиболее обоснованное число вторых соседей катиона равно 13 (рис. 66).

Рис. 6. Вторая координационная сфера катиона в структурном типе рутила -содержит 10 атомов, корунда - 13 атомов.

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

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

Для практической реализации предложенного метода моделирования твердых растворов замещения автором совместно с Деяновым Р.З. и Урусовым B.C. был написан пакет компьютерных программ ODSS, состоящий, помимо программы BINAR, еще из трех программ: RELAX, GISTOGRAMMA и VOLUME, используемых на пятом этапе для анализа локатьной структуры твердого раствора. Так как сверхъячейка в рамках вышеописанного подхода описывается симметрией Р\, расположение атомов не закреплено элементами симметрии и при оптимизации они смещаются со своих первоначальных позиций. Смещения атомов характеризуют релаксацию структуры твердого раствора. Автором предложена величина «сдвигаемость атома» <оь определяемая как усредненная по всем парам i-j среднеквадратичная разность их начального и конечного расстояний:

со, = £(Я„-Ryf/n. (2)

./ -i."

Здесь и - число частиц в сфере, внутри которой анализируются атомные смещения, Щ и Rначальное и конечное расстояние от атома i до атома /', соответственно. Помимо расчета сдвигаемостей атомов, в программе RELAX можно оценить податливость каждой кагаонной позиции C¡, определяемая по аналогии с определением Долласа [Dolíase W.A. (1980)] как реальная доля увеличения (или уменьшения) конкретных длин связей в полиэдре относительно разности длин связи чистых компонентов:

С, ^ " R . (3)

Здесь Л среднее значения межатомного расстояния дая конкретного катиона со своим ближайшим окружением в твердом растворе, -

аддитивное значение межатомного расстояния но правилу Вегарда (приближение виртуального кристалла), а Я - расстояние в чистом беспримесном кристалле. Отметим, что понятия «сдвигаемость атома» и

«податливость позиции» описывают несколько различные типы искажении в кристаллической структуре. Проиллюстрируем это на примере искажения октаэдра М06 (рис. 7).

9 ^Рс ЛР Рис. 7. Демонстрация

I £ = 0 ■ А-'отличий «податливости

'Я

л, ' I я, Ш О

позиции» С.«,- и

г" «сдвигаемости атома» со.

% ^

Ш Лк

Так, при смещении центрального катиона из геометрического центра октаэдра может возникнуть вариант, когда три связи катион-анион удлинятся, а три укоротятся практически на одну и ту же величину АЛ. Такая картина может наблюдаться, например, в корундовых системах при смещении катиона А1 (Сг, Ре и т.д.) вдоль оси с элементарной ячейки. В этом случае величина податливости позиции С,, будет либо равна нулю, либо крайне мала, что может быть ошибочно интерпретировано как отсутствие атомных смещений вообще. Напротив, значения со в этом случае будут отличны от нуля, со центрального катиона будет существенно превышать величины сдвигаемостей анионов, что позволит сделать вывод о том, что из стартовой позиции сдвигается именно катион.

Программа вКТОСКАММА позволяет провести анализ конечного

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

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

V -V

СГ = т"'р (4)

V -V

аддит

Здесь Утв р_р - усредненное по всем катионам одного сорта в сверхячеике значение объема координационного полиэдра, - аддитивное значение

объема полиэдра в соответствии с правилом Ретгерса, а V - объем полиэдра в

чистом беспримесном кристалле. Для оценки объемов тетраэдров и октаэдров на основами алгоритма Робинсона [Robinson К. et al (1971)] автором была написана программа VOLUME, проводящая геометрический анализ (включая расчет объемов) полиэдров. Подробное описание всех программ доступно по адресу http://cryst.geol.msu.ro/odss/.

Разработанный пакет программ был использован при компьютерном моделировании ряда изоморфных систем минералов [Урусов B.C. и др. (2007), Eremin N.N., Deyanov R.Z., Urusov V.S. (2008), Урусов B.C., Петрова Т.Г., Еремин H.ff. (2008), Еремин Н.Н., Талис Р.А., Урусов B.C. (2008)]. Конкретные результаты расчетов свойств смешения и локальной структуры твердых растворов представлены в главах 5 и 6.

Глава 3. Уточнение параметров межатомных потенциалов, структурное моделирование, расчет и предсказание свойств кристаллов различных минералогических классов.

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

3-1-1. Моделирование купратов щелочноземельных металлов.

Благодаря исследованиям высокотемпературной сверхпроводимости, открытой в 1986 году, для класса купратов в последние два десятилетия был накоплен значительный объем кристаллохимических данных [Leonyuk L.I. et al (1998)]. Однако работы по структурному моделированию этих соединений были крайне малочисленны, что связано с тем, что С11О4 квадратные группировки практически не поддаются корректному описанию с использованием традиционных приемов структурного моделирования. В работе в рамках модели формальных зарядов был разработан набор потенциалов, способный корректно описать кристаллические структуры купратов Са, Sr и Ва. Проведенное структурное моделирование Са2СиОз, CaCu02, SrCu03, (Sri 19Сао.7з)Си204 и ВаСиСь позволило предсказать ряд их неизвестных энергетических и упругих свойств. Разработанный набор потенциалов оказался хорошо переносимым в схожие структуры, несмотря на различия в мотивах соединения группировок С11-О4. Впоследствии он был рекомендован авторами

программы GULP для структурного моделирования купратов, содержащих С11О4 квадратные группировки [GaleJ.D., RohlA.L. (2003)].

3-1-2. Система кианит-андалузит-силлимапит.

На примере полиморфных модификаций AI2S1O5 (кианита, андалузита, силлиманита и двух возможных фаз высокого давления) было показано, как с помощью ионной модели можно относительно успешно предсказывать структуру, упругие свойства минералов, уравнение состояния, фононные спектры и термодинамических свойства кристаллических веществ. Однако, показан целый ряд систематических ошибок ионной модели, связанный с завышением вклада электростатического взаимодействия. В частности, показаны заметные искажения структуры силлиманита и гипотетической фазы высокого давления со структурой V3O5. Для упругих свойств кристаллов получены хорошие (в целом несколько завышенные) полуколичественные оценки, а для диэлектрических свойств - качественные. С удовлетворительной точностью были рассчитаны фононные спектры, в частности, фононный спектр кианита. Вычисленные термодинамические функции неплохо согласуются с экспериментом. Вместе с тем, попытка воспроизведения фазовой диаграммы AbSiOs в рамках ионной модели оказалась неудачной. Для исправления этой и других ошибок следовало учесть частичную ковалентность химических связей, что было сделано в §3.3.2 для кианита.

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

3-2-1. Структурное моделирование фосфатов структурного типа насикона.

Структурный тип NaZr2(PO,t)3 (NZP, насикон) широко распространен среди неорганических солей с тетраэдрическими анионами РО43". Благодаря своим свойствам (термической, химической и радиационной устойчивости, малому тепловому расширению, низкой теплопроводности, хорошей ионной проводимости и т.д.) они являются предметом интенсивных теоретических и экспериментальных исследований. Однако исследования по их структурному моделированию до последнего времени не проводились. В работе была разработана модель межатомных потенциалов и проведено компьютерное моделирование кристаллической структуры щелочно-циркониевых фосфатов MZr2(P04)3, где М = Li, Na, К, Rb, Cs и Fr, с использованием частично-ковалентного приближения. Был предсказан ряд их неизвестных энергетических и упругих свойств.

3-2-2. Разработка набора межатомных потенциалов для структурного моделирования дигидрофосфата калия КОР.

Дигидрофосфат калия КН2РО4 (КОР) является одним из важнейших сегнетоэлектриков. Простота получения из низкотемпературных водных растворов кристаллов большого размера, высокое оптическое качество и стойкость к лазерному излучению ставит КГ)Р на одно из первых мест при использовании в квантовой электронике и лазерной технике. Для моделирования структуры КОР было выбрано частично ковалентное приближение, как соответствующее реальному характеру химической связи в соединении. Параметры парных потенциалов межатомного взаимодействия в форме Букингема и Морзе оптимизировались с использованием экспериментальной информации о кристаллической структуре и упругих свойствах соединения. Разработанная модель позволила воспроизвести параметры элементарной ячейки КОР с ошибкой не более 1.8%, а атомные координаты - не более 2.5%. Добавочным параметром оценки качества проведенного моделирования и предсказательной возможности оптимизированного набора потенциалов явилось согласие рассчитанной и экспериментальной величины энтальпии образования соединения. Это позволило применить разработанный набор потенциалов в расчетах дефектных областей структуры КОР (глава 4).

3.2.3. Особенности структурного моделирования соединений, содержащих органические комплексы на примере бифталата калия С8Н5О4К.

Один из известнейших нелинейно-оптических кристаллов бифталат калия С8Н504К (БФК) (рис. 8) содержит в своей кристаллической структуре как органический комплекс С8Н504, гак и неорганический ион К. В связи с этим, роль кулоновского взаимодействия достаточно велика по сравнению с чисто органическими соединениями. Следовательно, использование потенциалов, разработанных для органических кристаллов, не приводит (и не может привести) к удовлетворительному воспроизведению кристаллической структуры. С другой стороны, очевидно, что атом калия в структуре БФК не несет на себе полный положительный заряд +1, а находится в некотором промежуточном состоянии, и использование известных из моделирования структур неорганических соединений потенциалов парного взаимодействия КО также не может привести к положительному результату. В связи с этим структурное моделирование БФК осуществлялось в рамках частично ковалентного приближения. Помимо кулоновского взаимодействия, использовались потенциалы Морзе и Букингема, а также угловые

трехчастичные гармонические потенциалы и потенциалы кручения. Использование разработанного оптимизированного набора потенциалов позволило превосходно воспроизвести структуру БФК: максимальные расхождения в параметрах элементарной ячейки составили всего лишь 0.23%, а атомные координаты отличались от экспериментальных значений не более чем на 0.012 А (в среднем - 0.001 А).

В §3.3 представлены результаты зарядовой оптимизации методом минимизации энергии атомизации, описанном во второй главе, продемонстрирована трансферабельность разработанных моделей потенциалов и показаны предсказательные возможности приближения на примере кварца, корунда, стишовита, кианита, малаяита Са5п08Ю4 и ряда оксидов переходных металлов.

3.3.1. Минимизация энергии атомизации кварца, корунда и стишовита

В качестве объектов исследования для отработки методики минимизации энергии атомизации были выбраны следующие соединения кремния и алюминия: а-кварц и стишовит 8Ю2 и корунд А120з. Для их структурного моделирования была выбрана модель парных потенциалов с введением квадратичных весовых функций для потенциалов Букингема и Морзе. Оказалось, что структуру кварца можно удовлетворительно воспроизвести лишь в очень узком интервале степеней данности /=50(6)%, в то время как структура стишовита воспроизводима при степенях ионности приблизительно от 45 до 70%. Структура корунда может быть с успехом рассчитана в пределах у=60(±15)%. Оптимальные величины атомных зарядов, соответствующие

Рис. 8. Мотив кристаллической структуры БФК. Полиэдрами показано окружение иона К.

минимумам на кривых зависимости энергии атомизации от степени данности связи равны: Z(Sirv) = +1,92е, Z(О) = -0.96е для кварца; Z(SiVi) = +2.04е, Z(О) = -1.02е для стишовита и Z(Al)= + 1.86е, Z(0) = -1.24е для корунда. Особенно убедительно преимущества модели с варьируемыми зарядами проявляются в предсказании энергии атомизации кристаллов, которая при оптимальных зарядах атомов оказывается весьма близкой к эксперименту. Неплохо предсказывается и энтальпия перехода кварц-стишовиг (расчет - 0.34 эВ, эксперимент - 0.51 эВ).

3.3.2. Минимизация энергии атомизации кианита А!2SiOs

Одним из критериев разумности и результативности предложенной модели с варьируемыми атомными зарядами является ее предсказательная способность. Для проверки было осуществлено структурное моделирование кианита, в котором варьирование зарядов осуществлялось одновременно на атомах Si и AI. Структура и свойства кианита были рассчитаны с использованием потенциалов, разработанных для кварца и корунда. Полученные для параметров ячейки, координат атомов и энергии сцепления результаты показали превосходное согласие с экспериментом. Более того, на момент проведения расчетов экспериментально упругие свойства кианита не изучались. Однако в 1997 году в работе [Yang el al. (1997)] сжимаемость кианита была экспериментально изучена и экспериментальная величина модуля всестороннего сжатия Кэкси, равная 193(1) гПа, оказалась в великолепном согласии с предсказанной по оптимальной модели Кщсч= 191 гПа.

/(AI-O)

Рис. 9. Поверхность энергии атомизации кианита Е как функция параметров ионности / связей 81-0 и А1-0.

На рис. 9 показан ландшафт поверхности энергии атомизации кианита А125Ю5 в координатах степеней ионности связей А1-0 и 81-0. Обратим внимание на то, что оптимальные значения ДА1-0) иДБьО) близко совпадают с соответствующими значениями для корунда и кварца, а значение эффективного

заряда кислорода (-1.15е) является промежуточным между его зарядами в корунде и кварце. Величина эффективного заряда А1 (+1.83е) крайне незначительно отличается от заряда А1 в корунде (+1.86е). Кроме того, нужно отметить чрезвычайно близкое совпадение рассчитанного (-50.96 эВ) и экспериментального (-51.27 эВ) значений энергии атомизации кианита.

3.3.3. Минимизация энергии атомизации малаяита CaStiOSiO*

В отличие от изоструктурного титанита CaTïOSiO,», в структуре малаяита атомы олова уже при комнатной температуре находятся в центре SnCV октаэдров, что обуславливает кристаллизацию минерала в пространственной группе Alla. Моделирование малаяита осуществлялось в частично ковалентном приближении методом минимизации энергии атомизации. Такое моделирование представляло самостоятельный интерес, так как упругие и диэлектрические свойства минерала экспериментально не исследовались. Проведенные расчеты позволили успешно воспроизвести структурные особенности малаяита, предсказать значения упругих, диэлектрических и энергетических характеристик соединения, которые согласуются с имеющейся ограниченной экспериментальной информацией о его свойствах. Поверхность £(ZS„,ZS1) имеет четко выраженный минимум при зарядах Si=+2.40(4)e и Sn=+3.20(4)e, что согласуется с мессбауэровскими данными, согласно которым заряд на ядре олова в малаяите выше, чем в S11O2 (для последнего расчет дает ZSn=+2.96e).

3.3.4. Минимизация энергии атомизации оксидов переходных металлов

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

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

3.4.1. Оксиды щелочноземельных металлов со структурой галита.

В структурном типе NaCl кристаллизуются несколько десятков природных и синтетических оксидов с формулой МО. Для них феноменологическая теория твердых растворов замещения разработана наиболее полно (Урусов, 1977). Естественно, что эти соединения являются одними из наиболее изученных объектов методом межатомных потенциалов, тем более что для них наиболее простая ионная модель может быть применена совершенно обоснованно. Однако, парные потенциалы с учетом неполной ионности химической связи гораздо лучше описывают упругие свойства кристаллов, поэтому в нашей работе предпочтение было отдано модели потенциалов со степенью ионности /=0.85.

3.4.2. Модель меясионных потенциалов для галита и сильвина.

Так как связь в галогенидах носит еще более ионный характер, чем в изоструктурных им оксидах, было принято решение не вводить в межатомные потенциалы для этой системы ковалентную составляющую. Пробные расчеты показали, что отталкивание во второй координационной сфере Ме-Ме и С1-С1 может быть описано исключительно электростатическими силами без использования короткодействующих потенциалов, поэтому в модели короткодействующих потенциалов были использованы только потенциалы Букингема Ме-С1. Разработанная модель обеспечивает точное описание структуры и упругих свойств NaCl и KCl. Температурные зависимости теплоемкости хорошо согласуются с экспериментом, однако зависимости энтропии от температуры испытывают систематические отклонения от экспериментальных величин, что связано, вероятно, с чрезмерной простотой набора парных потенциалов. Тем не менее, разность ^ксг^каа почти в точности соответствует экспериментально наблюдаемой, что позволяет использовать предлагаемую модель ионных потенциалов для дальнейшего моделирования твердых растворов в системе галит-сильвин.

3.4.3. Разработка согласованной модели потенциалов для полуторных оксидов М2О3 со структурой корунда.

Несмотря на то, что корунд неоднократно подвергался структурному моделированию ранее, критический анализ предлагаемых моделей показал необходимость разработать оригинальную согласованную модель для последующего моделирования твердых растворов состава VbOj. Заряд металла в разработанной модели был принят равным величине +1.89е, а совместное использование потенциала Букингема и Морзе позволило существенно

улучшить описание взаимодействия катион-анион. Некоторые рассчитанные структурные, упругие и термодинамические свойства корунда, гематита и эсколаита приведены в таблице 1 в сравнении с экспериментальными данными и расчетными данными из работы \Catlow С.НА. е( а/ (1993)].

Таблица 1. Рассчитанные свойства корунда, гематита и эсколаита в сравнении с _экспериментальными величинами._

Соединение а, к с. А К А' К, ГПа

А120з (эксперимент) \GaleJ.D. е1 а\ (1992)1 4.760 12.993 254.977 253.60 [Goto Т. et al (1989)1

АЬОз (расчет) 4.753 13.016 254.613 296.58

А1203 (расчет) \Ca1lowC.RA. е? а/(1993)1 4.7688 12.9869 255.77 333.22

СггОз (эксперимент) \IshimataN. etal(ШЩ 4.954 13.566 288.356 238 [Finger L, HazenRM. (1980)]

Сг20з (расчет) 4.965 13.538 288.993 230.26

БегОз (эксперимент) [АпИрт М.1и. е! ей (1985)]) 5.050 13.720 303.017 225 [Finger L„ HazenRM (1980)]

РегОз (расчет) 5.053 13.730 303.537 213.82

Как видно из таблицы, наша модель существенно успешнее описывает упругие свойства корунда, чем модель работы Кэтлоу, которая, более того, оказалась совершенно непригодной для моделирования термодинамических свойств соединения. С нашей точки зрения, это вызвано тем обстоятельством, что авторы вышеназванной работы выбрали несколько завышенные значения эффективных атомных зарядов +2.08 (А1) и -1.39 (О) (взятых из квантовомеханического расчета небольшого кластера [Gale J.D. et al (1992)]).

Рис. 10. Рассчитанные зависимости энтропии от температуры (сплошные линии):

(а) корунда,

(б) эсколаита,

(в) гематита

в сравнении с экспериментом (квадраты - корунд, треугольники эсколаит,

кружки - гематит)

Как видно из рис. 10 модель с зарядами +1.89 (Me) и -1.26 (О) успешно описывает энтропию корунда, во всем интересующем нас диапазоне

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

3.4.4. Разработка согласованного набора межатомных потенциалов для структурного моделирования оксидов и силикатов бериллия.

Методом атомистических потенциалов соединения бериллия практически не изучались, возможно, из-за отсутствия до сих пор рекомендованного потенциала межатомного взаимодействия Ве-0 в известных банках межатомных потенциалов. Оптимизация параметров короткодействующих потенциалов Ве-0 и О-Ве-О в рамках частично-ковалентного приближения проводилась по экспериментальным значениям параметров ячейки, координат атомов и упругих констант бромеллита ВеО, что позволило точнее воспроизвести кристаллическую структуру ВеО, чем расчет ab-initio [Шеин И.Р. и др. (2007)]. Рассчитанный модуль всестороннего сжатия отличается от экспериментальной величины всего на 5% (расчет ab-initio - 21%). С целью проверки трансферабельности разработанного набора параметров потенциалов Ве-0 он был использован при моделировании структур и свойств хризоберилла ВеА1204 и его изоструктурных аналогов: ВеСг204 и неизученного до сих пор рентгеновскими методами BeFe204. Рассчитанные параметры элементарной ячейки хризобериллов отличаются от экспериментальных величин не более чем на 1%, а координаты атомов не более чем на 1.75%. В целом, хорошо воспроизвелся и тензор упругих констант BeAb04. Как видно из рис. 11 энтропия соединений согласуется с экспериментальными значениями в пределах ±3% в интервале температур от 298 до 1800 К.

450 «1С 3S0 300 250 20С 150 »ОС 50

Энтропия, -0

Дж1(Моль'К)

........ ....... , ^

_......: ....... .................._........

УХ

/ ................•*©

т,к

Рис. 11. Рассчитанные зависимости энтропии от температуры

(сплошные линии):

(а) ВеА1204,

(б) ВеСг204,

(в) ВеРе204,

(г) ВеО в сравнении с экспериментальными данными для ВеА1204 и ВеО (кружки).

400

600

1£>00 1200 1400 1600 1800

госо

Таблица 2. Рассчитанные структурные и механические свойства BeFe2C>4

Предсказанная кристаллическая структура

Объем, А3 263.883

а, Ь, с, А 4.580, 9.921, 5.808

Fei х,у, z 0, 0,0

Fe2 х, у, z 0.9926, 0.2774, 0.25

В ex,y,z 0.4319, 0.0945,0.25

Ol x,y,z 0.7759, 0.0910,0.25

02 х, у, z 0.2400, 0.4388, 0.25

02 х, у, z 0.2668,0.1599, -0.9751

Расстояния Ве-О, А 1x1.5759,2x1.6439, 1x1.7339

Расстояния Fe 1-0, А 2x1.9735, 2x1.9942, 2x2.0079

Расстояния Fe2-0, А 1x1.9612, 2x2.0011, 1x2.0987, 2x2.1552

Предсказанные упругие свойства, гПа (усреднение Хилла)

Сц, С22, Сзз, С44, С55, Сбб 409.61, 323.45, 364.08, 97.49, 113.62, 112.42

Ст2, Стз, С23, К, G 119.50, 125.47, 119.37, 201.90, 112.91

Скорость прохождения акустических колебаний, км/сек

S-продольные, /'-поперечные 15.584, 27.533

Рассчитанные структурные и упругие свойства минералов также хорошо согласуются с имеющейся экспериментальной информацией, что позволяет считать ряд рассчитанных свойств этих соединений и модельную кристаллическую структуру Fe-хризоберилла (табл. 2) достоверным предсказанием. Это дает основание, во-первых, использовать разработанный набор потенциалов для моделирования свойств смешения и локальной структуры твердых растворов в хризобериллах (главы 5 и 6), а во-вторых, использовать его для моделирования других соединений, содержащих ВеС>4 тетраэдрические группировки, например фенакита ВегБЮ.» и берилла ВезAbiSiftOje] • Кристаллические структуры этих силикатов хорошо известны и могут быть описаны как смешанные каркасы из тетраэдров ВеО.» и SiC>4. Для их моделирования с использованием разработанного набора потенциалов потребовалось оптимизировать потенциал взаимодействия Si-O для заряда на атоме Si +2.52е (для соблюдения электронейтральности элементарных ячеек). Потенциалы Si-О в форме Букингема и Si-O-O в трехчастичной гармонической форме оптимизировались на кристаллической структуре «-кварца.

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

воспроизвел кристаллические структуру берилла и фенакита (отклонения параметров элементарных ячеек менее 1%) при разумном описании его упругих и термодинамических свойств. Для фенакита модуль всестороннего сжатия расходится с экспериментом не более 16%, а энтропия - не более 3% (Т= 3001800 К). Для берилла модуль всестороннего сжатия отличается от экспериментальной величины на 15%, а отклонения энтропии от экспериментальных значений не превышают 5% при температурах выше 900 К, и лишь при более низких температурах значения энтропии занижены по сравнению с экспериментальными значениями.

3.4.5. Разработка модели межатомных взаимодействий для оксидов группы рутила

Рутил и касситерит неоднократно и с переменным успехом подвергались структурному моделированию методом атомных потенциалов. Проведенный в работе анализ предыдущих исследований продемонстрировал необходимость разработки самостоятельной согласованной модели потенциалов. Наилучшее согласие экспериментальных и расчетных величин структурных, упругих и термодинамических характеристик этих оксидов удалось достичь при зарядах И и Бп равным +2.8, а кислорода - -1.4. Как видно из таблицы 3, разработанная для рутила и касситерита согласованная модель межатомных потенциалов неплохо описывает их структурные, упругие (а также термодинамические) свойства и может быть использована для проведения дальнейших расчетов в изоморфной системе (8п,Т1)02.

Таблица 3. Сравнение вычисленных и экспериментальных [Park el al. (1985), _Ahrens (1995)1 структурных и упругих свойств ТЮ2 и Sn02._

Параметр тю2 Sn02

эксперимент расчет эксперимент расчет

а, с (А) 4.5943, 2.9586 4.5993, 2.9857 4.737, 3.187 4.7623,3.1591

*(0) 0.305 0.3084 0.306 0.3059

v,ks 62.449 63.158 71.559 71.645

Си, С,2 (ГПа) 26.9, 17.7 32.0, 10.8 26.3, 17.8 31.7, 16.6

С13, Сзз (ГПа) 14.8, 48.2 18.6, 56.4 15.6, 45.0 14.4, 49.2

С 44, Сб6 (ГПа) 12.4, 19.3 13.7, 15.9 10.3, 20.8 13.5, 17.6

а:, ГПа 218.6 226.6 217.7 223.2

Глава 4. Предсказание возможных форм вхождения и локализация структурных позиций примесей в моделируемых соединениях.

4.1. Локализация позиций примесей железа в структуре касситерита SnC>2

Для точного задания координат атомов в

области кристаллического дефекта с целью интерпретации мессбауэровской информации о величинах ГЭП на ядрах Fe в касситерите было проведено структурное моделирование S11O2 с точечным дефектом 2Sn4<" = 2Fe3+ + V(O). Структурная оптимизация области дефекта показала (рис. 12), что ионы Fe3+ смещаются в противоположных друг от друга направлениях, а находящийся между ними анион стремится к линии, соединяющей катионы Fe, но не достигает ее. Оптимизированные координаты дефектного кластера были использованы для расчета градиентов электрического поля на ядрах 57Fe в структуре S11O2 по программе LATTICE [Русаков B.C. (2000)]. Было отмечено, что использование такого оптимизированного кластера существенно улучшает согласие между рассчитанной и экспериментальной величиной ГЭП на ядре Fe в SnOi.

4.2. Моделирование дефектных центров в структуре ильменита.

Задача провести моделирование кристаллической структуры и локальных структурных искажений частично окисленного ильменита встала после того, как в работе [Глазкова МЛ. (1996)] было обнаружено, что по результатам обработки мессбауэровских спектров ионы Fe2 по величине изоморфных сдвигов можно разделить на две группы. К первой группе Fe2+(1) относятся ионы со сдвигами, характерными для октаэдрического кислородного окружения (5= 1.04 - 1.12 мм/с). Для второй группы Fe2+(2) значения 5 ~ 0.9 мм/с близки к значениям сдвигов, характерным для ионов Fe2+ в кислородном тетраэдрическом окружении. Проведенное структурное моделирование показало, что такие ионы Fe2+(2) находятся не в соседних октаэдрических слоях, а в одном с дефектом Fe3+ слое и связаны с ним общими атомами кислорода. Отметим также, что соотношение содержаний Fe2+(2)/Fe3+ . отвечает экспериментально обнаруженному отношению 3:2. При изоморфном замещении Fe2" + Ti4+ на 2Fe + в структуре ильменита происходит образование

Ф f

Рис. 12. Смещение атомов в дефектном центре РегОю

гематитового кластера, вытянутого вдоль оси с элементарной ячейки. Появление такого кластера способствует уменьшению этого параметра. Внутри гематитового кластера два Fe3+ - октаэдра геометрически неэквивалентны. В связи с этим, можно рекомендовать описывать Fe3+ при обработке мессбауэровского спектра не одним, а двумя дублетами. В области, смежной с центром дефекта, возникают три октаэдра Fe2+(2)-0 с двумя значительно удлиненными расстояниями металл-кислород, что отражается на мессбауэровском спектре в виде появления аномальной позиции Fe2t(2).

4.3. Моделирование дефектных центров, образованных двух- и трехвалентными металлами в структуре KDP.

Существенное влияние на физические свойства технически важных кристаллов дигидрофосфата калия (KDP) оказывают примеси двух- и трехвалентных металлов, в частности двух- и трехвалентного железа. Для определения количественных характеристик изменений в дефектной области в результате внедрения примесных ионов в кристаллы KDP было проведено структурное моделирование методом атомных потенциалов, которое установило, что примеси двух- и трехвалентных металлов формируют в структуре KDP три типа дефектных центров. Первый тип представлен изолированными дефектными центрами, формирующими «каналы» вдоль оси с кристалла (катионы М + и второй тип образован примесными цепочками «кластерного» типа, ориентированные под углом 60 - 75° к оси с кристалла (катионы Fe2 , Мп2+, частично Со2 ). Третий тип представляет собой дефектные центры в позициях калия, сформированные как с образованием вакансий, так и по механизму сопряженного гетеровалентного замещения. Был установлен факт меньшей деформации структуры KDP ионами М + по сравнению с что согласуется с экспериментально наблюдаемым характером изменения рельефа граней KDP в присутствии примесей.

4.4. Расчет энергии различных точечных дефектов в рутиле и касситерите Разработанный в § 3.4.5 согласованный потенциальный набор был

использован для расчета энергии образования вакансий и интерстиций атомов, дефектов Шоттки, Френкеля и вхождения одиночных примесей в ТЮ2 и Sn02. Полученные в настоящей работе оценки в целом не противоречат более ранним расчетам и экспериментальной величине энергии образования кислородной вакансии в рутиле (7.74 эВ - [Picard С., Gerdanian Р. (1975)]; настоящая работа 8.76 эВ), хотя разброс результатов расчета различными методами и моделями в ряде случаев очень существенен. Таким образом, можно сделать вывод, что

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

4.5, Расчет энергии различных точечных дефектов в полуторных оксидах М2О3 со структурой корунда

Методом атомистических потенциалов энергетика образования точечных дефектов в структуре корунда неоднократно изучалась различными исследователями на протяжении многих лет, но лишь в рамках формальной ионной модели. Проведенные расчеты по модели потенциалов из § 3-4-3 позволили оценить энергию образования таких структурных несовершенств, как пятерки Шоттки, катионных и анионных пар Френкеля, а также значения параметров взаимодействия О, используемых далее для расчета свойств смешения. Наиболее предпочтительными типами дефектов для всех структур являются пятерки Шоттки, а рассчитанная энергия образования дефектов согласуется существенно лучше предыдущих теоретических расчетов с экспериментальными оценками [Mohapatra S.K and Kroger F.A. (1978)]. Так, для А12Оз энергия образования дефекта Шоттки (на дефект) равна 3.83 эВ (эксперимент) и 4.50 эВ (расчет); энергия образования катионной пары Френкеля (на дефект) равна 4.45 эВ (эксперимент) и 5.08 эВ (расчет). Рассчитанные значения энергий бесконечного разбавления для системы корунд-эсколаит равны 18.78 и 17.53 кДж/моль, что хорошо согласуется с данными [Jacob К.Т. (1987)], полученными методом ЭДС (()[ 18.7 кДж/моль, Q2 15.8 кДж/моль) и теоретическим термодинамическим расчетом [Гусаров В.В. (1996)] {О, 18.3 кДж/моль, (Л 14.1 кДж/моль). Для системы корунд-гематит рассчитанные значения 0\ 32.97 кДж/моль и Q?= 27.81 кДж/моль также отлично согласуются с данными Гусарова, полученными теоретическим термодинамическим расчетом (Qr 32.2 кДж/моль, Qrr 21.0 кДж/моль) и термодинамическим анализом фазовых диаграмм (Qr 40.0±2.0 кДж/моль, Q2 23.6±1.8 кДж/моль). Такое согласие позволяет использовать разработанную согласованную модель межатомных потетщиалов для расчета свойств смешения в протяженных твердых растворах системы корунд-гематит-эсколаит.

4.6. Расчет энергии образования точечных дефектов в бромеллите и хризобериллах

Для расчетов собственных и примесных дефектов с помощью разработанной (§ 3.4.4) модели потенциалов была использована модель Мотта-Литтлтона с радиусами зоны дефекта в 6.5 А (около 160 атомов) и экранирующего слоя до 16.5 А (около 2700 атомов). Для хризоберилла и его изоструктурных аналогов наиболее предпочтительными дефектами оказались бериллиевые пары Френкеля, тогда как для бромеллита механизм образования дефектов Шоттки является на 1.3 эВ энергетически более выгодным, чем образование катионных пар Френкеля. Отметим, что результаты настоящей работы по энергии Ве-пар Френкеля в бромеллите в целом согласуются с ab-initio расчетом [Софронов А.А. и др. (2006)]. В этой работе в зависимости от выбора расчетного кластера, энергия образования дефекта составила от 2.04 до 4.42 эВ, при этом октаэдрическая позиция оказалась более выгодной для вхождения междоузельного Be, чем тетраэдрическая. Это согласуется с проведенным методом атомистических потенциалов расчетом: в результате оптимизации области дефекта атом Be мигрирует из стартовой тетраэдрической позиции в незаполненные октаэдрические пустоты, при этом энергия образования дефекта составляет 4.17 эВ. Во всех соединениях анионные пары Френкеля являются наиболее энергетически затратными и их образование маловероятно.

Глава 5. Компьютерное моделирование свойств смешения в бинарных и тройных системах твердых растворов.

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

5.1. Компьютерное моделирование свойств смешения и стабильности твердых растворов оксидов щелочноземельных металлов.

Расчеты велись в сверхячейке 4><4х4 без нетрансляционной симметрии. Катионы разного сорта распределялись с максимальной степенью неупорядоченности (по критерию Брэгга-Вильямса) по 256 катионным позициям в требуемом стехиометрическом соотношении. Структурные и термодинамические свойства неупорядоченных твердых растворов Mg^Ca^O, CavSr(i_X)0, SrtBa< 1-^)0 были рассчитаны для пяти составов х = 0.125, 0.25, 0.50,

0.75, 0.875 в интервале температуры 298 - 1800 К. Графическое построение зависимостей энергий Гиббса от состава в диапазоне температур от 298 до 1800 К (через 50 К) и поиск минимумов на этих кривых позволили приближенно оценить границы устойчивости твердых растворов (рис. 13). Они находятся в хорошем согласии как с экспериментальными данными, так и с результатами других теоретических расчетов. Особо следует подчеркнуть, что такие рассчитанные зависимости как Д У(х), АК(х), а также Д^ол еще недоступны в большинстве случаев для прямого экспериментального определения.

2200 р \ 2000 ' j 1 1800 • I

1S0Ö

i400 -г-—г——---1-

О 0.2 0.4* 0.6 0.8 1 MtO

0 0-2 0.4 0.6 0.6 1

ХСяО

Рис. 13. Кривые сольвуса в системе Т^Са^О \а), Са^Гп^О (б), Бг^Ва^О (в). Сплошная линия - атомистический расчет, пунктир - эксперимент \Jacob К.Т. е/ о/ (1995), (1998)], точечная линия - расчет Монте-Карло \LavrentievМ.У. е1 а\ (2001)].

5.2. Моделирование свойств смешения и стабильности твердого раствора галит - сильвин

Моделирование неупорядоченных твердых растворов галит-сильвин также осуществлялось в сверхячейке с учетверенными параметрами структуры типа NaCl, содержащей 512 атомов. Однако, в отличие от моделирования твердых растворов оксидов, описанных в предыдущем параграфе, для составов с x(Na) = 0.125, 0.25, 0.5, 0.75 и 0.875, неупорядоченные конфигурации ионов Na+ и К+ в сверхячейке определялись по программе BINAR. Сравнение с экспериментальными данными (рис. 14-15) показывает неплохое согласие различных методов оценки энтальпии смешения и границ устойчивости твердого раствора, хотя экспериментальная критическая температура распада несколько выше. По результатам моделирования значение критической температуры Гкр= 720±25 К при Хмаа = 0.57±0.05 (эксперимент: Ткр = 780 К).

Рис. 14. Энтальпия смешения твердого Рис. 15. Границы устойчивости раствора NaCl-KCl. Белые кружки - твердого раствора. Штриховая атомистический расчет, черные - полоса - атомистический расчет, эксперимент [Barett W.T., Wallace W.E. черные кружки - эксперимент 1954]; крестики - расчет Дика и Даса [Dick [Barett W.T., Wallace W.E. В. G., Das Т. P. (1962)], линия - феноменоло- (1954)]. гическая теория [Урусов B.C. (1977)].

5.3. Расчет свойств смешения и стабильности твердых растворов замещения в системе корунд-эсколаит-гематит

Для расчетов свойств смешения твердых растворов в системе корунд-эсколаит-гематит была выбрана сверхячейка 4><4><1 структурного типа корунда пр. гр. PI, содержащая 672 частицы (288 атомов кислорода и 192 катиона, состоящих из остова и оболочки). Расчеты свойств смешения проводились для различных бинарных составов (М'хМ^.^Оз с катионным соотношением М':М2 = 1:5, 1:2, 1:1, 2:1, 5:1. Для расчетов свойств смешения в тройной системе корунд-гематит-эсколаит дополнительно рассчитывались еще 10 промежуточных составов. Из рис. 16 видно, что для всех трех бинарных систем энтальпия смешения описывается параболами с достаточно малым параметром асимметрии. Для системы корунд-эсколаит рассчитанные значения А#см отлично согласуются с экспериментом [Chatterjee N.D. et al (1985)]. Как и следовало ожидать, минимальные величины энтальпии (менее 1 кДж во всем диапазоне составов) зафиксированы для системы эсколаит-гематит, что связано с близкими величинами ионных радиусов и электроотрицательностей атомов и предполагает отсутствие области распада твердого раствора при всех температурах, вплоть до комнатных.

Рис. 16. Л#см бинарных твердых растворов. Сплошная линия - расчет, пунктир и квадраты - экспериментальные данные для А120з-Ре20з треугольники - для А12Оз-Сг203.

Рис. 17. Кривые сольвуса в системах А12Оэ-Сг2Оз - 1 и А1203-Ре203 - 2. Расчет -сплошные линии, пунктирные линии - эксперимент.

Содержание СггО, {Fe А)

Несколько хуже согласуются с экспериментальной информацией результаты расчета энтальпии смешения для системы корунд - гематит. Расчетные данные систематически занижены по сравнению с данными [Majzlan J. et al (2002), Feenstra A. eí al (2005)]. Вероятно, это может быть связано с тем, что в энтальпию смешения этой системы вносят свой вклад более тонкие эффекты, например некоторое различие эффективных зарядов атомов Al и Fe, а также различные магнитные свойства компонентов, которые не учтены в примененной модели. Любопытен факт отрицательных, хотя и крайне малых, отклонений объемов от правила Ретгерса для всех бинарных и тройных составов твердого раствора. Отклонения модулей всестороннего сжатия от аддитивности также отрицательны для всех бинарных и тройных составов. Однако, если для системы гематит-эсколаит отклонения очень малы, то для двух других систем они составляют до 7 (А1-Сг) и даже 12 ГПа (Al-Fe), что может быть подвергнуто экспериментальной проверке. Теоретические диаграммы распада твердых растворов корунд - эсколаит и корунд - гематит приведены на рис. 17.

Для системы корунд-эсколаит критическая температура оценена как 1215±20 К. что хорошо согласуется с другими оценками этой величины, в том числе экспериментальными. Отметим, что кривая сольвуса относительно симметрична и её максимум близок к среднему составу л*(Сг)-0.47, в согласии с экспериментальными данными. В системе корунд-гематит по данным [Feensira A. el al (2005)] при 7-1400 К и выше образуется ромбическая фаза А1РеОз, поэтому критическую температуру распада можно только приближенно оценить как 1550 (расчет) и 1800 К (эксперимент).

5.4. Моделирование свойств смешения в системе ВеА120гВеСг204 с учетом катионногораспределения катионов по октаэдрическим позициям

Для расчетов кристаллических структур и свойств смешения твердых растворов хризобериллов была выбрана сверхячейка 4*2x2 структурного типа оливина, содержащая 128 атомов Al (Сг), 64 атома Ве и 256 атомов О. Многочисленные экспериментальные факты указывают на существенное предпочтение ионами Сг более крупных и менее симметричных позиций М2. Анализ экспериментальных данных вынудил автора усложнить процедуру выбора атомной конфигурации для различных составов твердого раствора BeAl204-BeCr204. Для составов, обогащенных алюминием, катионное распределение ионов хрома по 64 позициям М1 и 64 позициям М2 в сверхячейке выдерживалось в соответствии с коэффициентом 1:2, определенным в работе [Rabadanov М.К. and Dudka А.Р. (1998)]. Соответственно, для составов, обогащенных хромом, распределение атомов алюминия выдерживалось с коэффициентом А1^(М2)/А13+(М1), равным 1:2. Программа BINAR применялась для размещения определенного количества катионов одного copra не по 128, а по 64 позициям М1 и но 64 позициям М2 раздельно. Рассчитанная зависимость энтальпии смешения от состава с учетом катионного распределения представлена на рис. 18.

Как и в случае корундовых систем наблюдается крайне малое отрицательное отклонение объема от аддитивности: наибольшее отклонение AV—0.15 Á\ Такое отклонение, по-видимому, не может быть сопоставлено с экспериментальными оценками. Отклонения модулей всестороннего сжатия от аддитивности также отрицательны для всех составов твердого раствора. Величины отклонений малы и не превышают 3.5 ITIa.

На основе графических построений были найдены составы сосуществующих твердых растворов, отвечающих, минимумам на кривых AG(x) при Т> 1000 К. Критическая температура оценена как 2150±20 К при а- Сг=0.51 (рис. 19), т.е. практически соответствует температуре плавления хризоберилла

Гпл=2143 К [Гапахов Ф.Я. (1957)], следовательно купол сольвуса в действительности не замыкается. Более того, согласно данным [Гусаров В.В., Семин Е.Г. (1992)] при температуре 2016 К хризоберилл претерпевает фазовый переход в фазу с тетраэдрической координацией А1.

Рис. 18. Энтальпия смешения твердого Рис. 19. Кривая сольвуса в системе

раствора ВеА1204-ВеСг204 с учетом ВеА1204-ВеСг204. Пунктиром

катионного распределения ионов А1 и показано положение критической

Сг по октаэдрическим позициям. точки.

Глава 6. Анализ локальной структуры твердых растворов замещения по данным компьютерного моделирования.

Одним из самых непростых для экспериментального изучения любого твердого раствора является вопрос о его реальной локальной структуре. Рентгеноструктурный анализ дает сведения только о так называемой «средней» структуре твердого раствора. Экспериментальные сведения о реальной локальной структуре могут быть получены только с помощью спектроскопических методов, если они доступны для рассматриваемой системы. Одним го них является метод ЕХАББ, который требует применения синхротро мпо го излучения. Поэтому теоретический анализ локальной структуры интересен как сам по себе, так и как инструмент для расшифровки ЕХАРБ спектров.

6.1. Локальная структура твердых растворов Mgo.5Cao.5O и Cao.5Sro.5O.

Методом атомистических потенциалов локальная структура твердых растворов оксидов со структурой №С1 изучалась для эквимолярных составов систем Са0-М§0 и СаО-БтО. Как видно на рис. 20, межатомные расстояния М§—О в твердом растворе в целом увеличиваются относительно расстояний в структуре чистого периклаза, а расстояния Са-О соответственно уменьшаются относительно расстояний в структуре извести. Аналогичная картина наблюдается и для максимально неупорядоченного твердого раствора

Cao.5Sro.5O. По данным моделирования возвращение к исходной структуре чистых компонент (релаксация структуры твердого раствора) больше, чем это предписывается геометрической ДЛС моделью [Dollace W.A. (1980)]. Отметим, что дисперсия межатомных расстояний 0-0 намного выше дисперсии расстояний Mg-O, Са-О и Sr-O. Это значит, что именно анионы сильнее смещаются из своих регулярных позиций, что находится в согласии с

Рис. 20. Частотная диаграмма межатомных расстояний в твердом растворе эквимолярного состава Mgo.5Cao.5O. В ертикальными линиями показаны расстояния между

ионами первой и второй координационной сферы в чистых компонентах MgO и СаО.

6.2. Локальная структура твердого раствора галит-сильвин

Расчет локальной структуры твердого раствора галит-сильвин проводился в сверхячейке 4x4x4, содержащей 128 ионов К и 128 ионов Na. Расчет показал, что межатомные расстояния Na-Cl увеличиваются относительно расстояний в галите примерно на 1/4 AR, а расстояния К-С1 уменьшаются относительно расстояний в сильвине на ту же величину. Следовательно, релаксация структуры твердого раствора несколько больше, чем это предписывается феноменологической моделью [Урусов B.C. и др. (1997)], согласно которой с наибольшей вероятностью должны встречаться изменения расстояний ± AR/2. Рассчитанный параметр релаксации X = 0.62 находится в хорошем согласии с данными исследования локальной структуры твердых растворов (K,Rb)Br и Rb(Br,I) методом EXAFS: X для первого составляет 0.55, а для второго 0.65 [Воусе J.B., Mikkelsen J.C. (1985)]. Отметим, что дисперсия межатомных расстояний С1-С1 намного выше дисперсии расстояний Na-Cl и К-С1. Это значит, что именно анионы СГ сильнее смещаются из своих регулярных позиций, что согласуется с представлениями феноменологической теории, которая не учитывает смещений катионов, допуская, что искажения происходят только в анионной подрешетке.

представлениями феноменологической теории.

Доля связей

0.8

0.6

0.4

0.2

о 1.8

1-зякоординационная сфера

2.102 2.406

/Л

2-ая координационная сфера

2.973 3,168 3.403

)-ОДта-р-р;

R, А

2.3

3.3

3.8

6.3. Анализ локальной структуры эквимолярных твердых растворов (А110,Сг10)О3, (Те1МСг1Л)03, (А1,л,ГсГп)Оз

Анализ локальной структуры твердых растворов замещения на основе корунда проводился для максимально неупорядоченных конфигураций эквимолярных составов в рамках сверхъячейки 4x4^1. Анализ атомных смещений показал, что во всех случаях вклад анионной подрешетки в искажение идеальной структуры существенно выше катионной. На рис. 21 видно, что три более длинных расстояний металл-кислород (М-ЗОЛЛШ!) испытывают существенно большие изменения по сравнению с изменениями трех более коротких расстояний в октаэдре (М-30КоР). Этот факт говорит о том, что в твердом растворе октаэдры в основном искажаются за счет длинных связей М-О, а группировка М-ЗОкор подвержена меньшим искажениям.

Рис. 21. Частотное распределение расстояний М-0 в твердом растворе (А10.5Сг0.5)2О3. Средние положения пиков А1-0 и Сг-О в твердом растворе -вертикальные пунктирные линии; соответствующие расстояния в исходных оксидах - сплошные линии.

Отметим, что перекрывание пиков межатомных расстояний даже для первой координационной сферы металл-кислород для разных катионов усложняет интерпретацию наблюдаемой картины релаксации координационных полиэдров. В этом случае логично перейти от рассмотрения конкретных межатомных расстояний в полиэдре к анализу изменений полиэдрических объемов, что может давать более наглядные гистограммы распределений. Так, гистограмма, представленная на рис. 22, отчетливо показывает, что истинная релаксация объемов октаэдров занимает промежуточное положение между моделями виртуального кристалла и чередования связей. Расчет по программе Relax эквимолярного состава (А1о.5Сго.5)2Оз показал, что ионы Сг и A1 практически одинаковы по сдвигаемости: со(Сг)=0.00196(63) А2, а со(А1)=0.00191(74) А2. Вдвое большими значениями сдвигаемости характеризуется кислородная подрешетка: со(0)=

:Чисяо расстояний |» сюрхачсйхс

J Сг-О.,,

i i / ckktt'kwt)

Г I / ...

■ i 7 <И,Р>,1Д'

.....А!-0 .

-а-о

W> кор, -.......

■ ■ \ п

1.82 1.94 «.86 1.88 j..9c

.1.93 1.9в .2.00 2.112 2.с* 2,0?. 208 2.1 с! • 2.12

0.00380(60) А2, что подтверждает вывод о большем вкладе анионов в деформацию кристаллической структуры.

,,, Число октаэдров.........................................

всверхячейке

I сю, \

; V \>ра'лфа. Á

Рис. 22. Частотное распределение объемов октаэдров МОб в твердом растворе (Alo5Сго.5)203. Сплошными вертикальными линиями показаны

соответствующие объемы в чистых оксидах, аддитивный объем отмечен вертикальной пунктирной линией.

6.4. Анализ локальной структуры эквимолярных твердых растворов на основе хризоберилла.

Для твердых растворов на основе хризоберилла анализ локальной структуры проводился в рамках сверхячейки 4x2x2 структурного типа оливина для эквимолярного состава бинарной системы ВеА1204-ВеСг204 с учетом катионного распределения катионов А1 и Сг по позициям по позициям М1 и М2, описанного в предыдущей главе. Были построены гистограммы распределения расстояний М-О, М-М и 0-0 в твердом растворе, а также частотное распределение объемов октаэдров М06. Как видно из таблицы 4 средняя степень релаксации Ст06 октаэдра в хризобериллах >. = 67%. Эта величина несколько меньше, чем для гранатов, но существенно больше, чем в системе корунд-эсколаит. Данный факт является неплохой иллюстрацией известного правила содействия (ассистирования) изоморфизма: общая структурная единица изоморфной смеси не является нейтральной, увеличение ее относительного размера расширяет область смесимости.

Таблица 4. Оценки степени релаксации СгОб октаэдров в изоморфных системах

Система Метод

Ca3(Al1_xCrx)2SÍ3012 гроссуляр-уваровит 0.82 Анализ электронно-абсорбционных спектров [Taran M.N. et al (2004)1

Mg3(Al1.xCrx)2Si3012 пироп-кноррингит 0.77 Анализ электронно-абсорбционных спектров [Langer К. et al (2004)]

Be(Al0.67Cr033)M1(Al0.33Cr067)M2O4 0.67 Атомистическое моделирование

Ali оСг] 0О3 корунд-эсколаит 0.45 Атомистическое моделирование

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

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

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

ЗАКЛЮЧЕНИЕ

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

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

ОСНОВНЫЕ ВЫВОДЫ

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

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

3. Использование частично ковалентной модели при моделировании точечных дефектов по процедуре Мотта-Литтлтона позволило сделать более обоснованные выводы о позиционном предпочтении и механизме вхождения примесей в кристаллические структуры ряда неорганических соединений.

4. Разработанные автором оптимизированные наборы потенциалов были успешно применены для структурного моделирования и предсказания физических свойств купратов, содержащих Си04 квадратные группировки, фосфатов структурного типа насикона, оксидов и силикатов, содержащих катионы А1, Бе, и Сг в октаэдрическом окружении а также минералов, содержащих Ве04 - тетраэдры.

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

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

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ.

1 Урусов B.C., Дубровииский Л.С., Вассерман Е.А., Еремин H.H. (1994) Моделирование структуры и свойств оксидов типа рутила методом минимизации энергии межатомного взаимодействия // Кристаллография, Т. 39, стр. 446-456.

2 Urusov V.S., Eremin N.N. (1995) Energy Minimum Criteria in Modeling Structures and Properties of Minerals // Phys. Chem. Minerals, T. 22, стр. 151-158.

3. Храмов Д.А., Русаков B.C., Еремин H.H., Урусов B.C. (1994) Градиенты электрических полей на ядрах Sn и Ti в оксидах со структурой рутила // Известия РАН, сер. физическая. Т. 58, №4, стр. 30-34.

4. Урусов B.C., Якубович О.В., Еремин H.H. (1994) Деформационная электронная плотность в касситерите SnCb // Доклады РАН, Т. 336, №3, стр. 330-334.

5. Урусов B.C., Еремин H.H., Якубович О.В. (1995) Распределение электронной плотности в титаните CafTiOfSiO.»]} // Кристаллография, Т. 40, стр. 485-492.

6. Еремин H.H. (1996) Теоретическое моделирование и прецизионные экспериментальные исследования оксидов и силикатов олова и титана // Диссертация на соискание ученой степени кандидата химических наук, Институт геохимии и аналитической химии РАН, Москва, 146 стр.

7. Urusov V.S., Eremin N.N. (1996) Structure modeling of inorganic crystals using charge transfer energy term // Phys. Chem. Minerals V.23, p. 216.

8. Urusov V.S., Eremin N.N. (1997) Charge-transfer energy in computer modeling of structure and properties of minerals // Phys. Chem. Minerals V.24, №5, P. 374-383.

9. Урусов B.C., Оганов A.P., Еремин H.H. (1997) Компьютерное моделирование структуры, свойств и устойчивости модификаций AbSiOs. I. Ионное приближение // Геохимия, №5, стр. 456-474.

10. Еремина 'Г.А., Кузнецов В.А., Охримеико Т.М., Фурманова Н.Г., Еремин H.H., Урусов B.C. (1998) Моделирование дефектной области в кристаллах КДР с примесью трехвалентного железа // Кристаллография т.43, №5, стр. 906-911.

11. Урусов B.C., Еремин H.H. (1999) Компьютерное моделирование структуры и свойств кристаллов - современные достижения и возможности. // Сборник «Проблемы кристаллологии». Вып. 5. Москва, Изд-во ГЕОС, стр. 228-266.

12. Урусов B.C., Еремин H.H., Оганов А.Р. (1999) Моделирование структуры и свойств оксидных кристаллов методом минимизации энергии атомизации // Кристаллография, Т. 44, стр. 396-405.

13. Игнатьев В.Д., Еремин H.H. (1999) Теоретическое моделирование дефектов в структуре ильменита // Вестник Института Геологии (РАН, Уральское отд. Коми научный центр) Лг»10, стр. 3-7.

14. Кузнецов В.А., Еремина Т.А., Еремин H.H., Охрименко Т.М., Фурманова Н.Г., Ефремова E.H. (2000) Структурный механизм влияния примесей на кинетику роста кристаллов. Кристаллохимический аспект // Кристаллография, Т. 45, Мб, стр. 1.120-1125.

15. Еремина Т.А. , Еремин H.H., Фурманова Н.Г.. Кузнецов В.А., Охрименко Т.М., Урусов B.C. (2001) Моделирование дефектной области в кристаллах KDP с примесью ионов двухвалентного железа. Сравнение дефектов, создаваемых двух- и трехвалентными металлами // Кристаллография, Т. 46, №1, 2001, стр. 82-87.

16. Дудникова В.Б. Жариков Е.В. Еремин H.H., Жаркова Е.В., Лебедев В.Ф., Урусов B.C. (2001) Распределение ванадия между кристаллом и расплавом форстерита, структурное и валентное состояние ванадия /7 Геохимия №3, стр. 734-743.

17. Еремин H.H., Игнатьев БД., Урусов B.C. (2001) Компьютерное моделирование дефектов в структуре измененного ильменита // Вестник МГУ, сер. 4 - Геология, №6, сгр. 25-28.

18. Eremin N. N., Leonyuk L. /., Uriisov V. S. (2001) Interatomic Potentials for Structure Simulation of Alkaline-Earth Cuprates // J. Solid State Chem. V. 158, pp. 162-168.

19. Урусов B.C., Еремин H.H. (2002) Развитие кристаллохимических методов предсказания и моделирования структуры, дефектов и свойств минералов и их твердых растворов. В сборнике: «Геология, геохимия и геофизика на рубеже XX и XXI веков». Т.2. С.338-387.

20. Еремин H.H., Урусов B.C.; Русаков В.С, Якубович О.В. (2002) Прецизионные рентгеновские, мессбаузровские исследования и компьютерное моделирование малаяита CaSnOSiO.<, /У Кристаллография, Т.47 №5, стр. 825-833.

21. Veremeichik Т. F., Zharikov Е. V., Dudnikova V. В., Gaisier А. V., Eremin N. N., Smirnov V. A., Soubbotin K. A., Urusov V. S., Shcherbakov J. A. (2002) Vanadium in forsterite: oxidation states and stmctural localization in the crystals grown by the Czochralski technique // Opt, Mater. V.19, № 3, p. 319-328.

22. Еремина Т.А., Еремин H.H., Кузнецов В.А., Фурманова Н.Г., Урусов B.C. (2002) Моделирование дефектов, образованных катионами двух- и трехвалентных металлов в структуре дигидрофосфата калия. Методика расчетов // Кристаллография, Т. 47, №5 стр. 819-824.

23. Урусов B.C., Петрова Т.Е., Еремин H.H. (2002) Компьютерное моделирование свойств твердых растворов MgO-CaO с учетом ближнего порядка // Доклады Академии Наук, 'Г. 387, №2, стр. 191-195.

24. Урусов B.C., Петрова Т.Г., Еремин H.H. (2003) Моделирование локальной структуры и свойств твердых растворов CaO-SrO и SrO-BaO // Доклады Академии Наук, Т.392, №1, стр. 1-6.

25. Еремин H.H., Суханов М. В., Петьков В. И., Урусов В. С, (2004) Межатомные потенциалы для структурного моделирования двойных

ортофосфатов циркония и щелочных элементов // Доклады Академии Наук. Т. 396. № 4. стр. 496-499.

26. Урусов B.C., Еремин ПЛ. (2004) Компьютерное моделирование структур и свойств минералов с использованием атомистических потенциалов межатомного взаимодействия // Вестник МГУ сер. 4 - Геология, №5, стр. 37-44.

27. Rak М., Eremin N.N., Eremina Т.А., Kuznetsav V.A., Okhrimenko Т.At, Furmanova N.G., Efremova E.P. (2005) On the mechanism of impurity influence on growth kinetics and surface morphology of KDP crystals—I: defect centres formed

by bivalent and trivalent impurity ions incorporated in KDP structure.......theoretical

study // Journal of Crystal Growth V.273, pp. 577-585.

28. Eremina T.A., Kuznelsov V.A., Eremin N.N., Okhrimenko T.M., Furmanova N.G., Efremova E.P., Rak M. (2005) On the mechanism of impurity influence on growth kinetics and surface morphology of KDP crystals — II: experimental study of influence of bivalent and trivalent impurity ions on growth kinetics and surface morphology of KDP crystals /'/' Journal of Crystal Growth V.273, pp. 586-593.

29. Суханов M. В., Петьков В. И., Куражковская B.C., Еремин Н.Н., Урусов В. С. (2006) Компьютерное моделирование структуры, синтез и изучение фазообразования молибдатофосфатов Ai.yZr2(P04)3.x(Mo04)x (А-щелочной металл) // Журнал неорг анической химии, Т. 51, №5, стр. 773-779.

30. Урусов B.C., Петрова Т.Е., Леоненко Е. В. Еремин Н.Н. (2007) Компьютерное моделирование локальной структуры, свойств и стабильности твердого раствора галит NaCl - сильвин КО. Вестник Московского университета сер. 4 - Геология №2, стр. 58-63.

31. Еремин П.Н., Талис Р.А., Урусов B.C. (2008) Компьютерное моделирование локальной структуры, свойств смешения и стабильности бинарных оксидных твердых растворов замещения системы корунд-эсколаит-гематит // Кристаллмрафия, Т.53, №5, стр. 802-810.

32. Еремин ПЛ., Деянов Р.З, Урусов B.C. (2008) Выбор сверхячейки с оптимальной атомной конфигурацией при моделировании неупорядоченных твердых растворов // Физика и химия стекла, Т 34, №1, стр. 9-18.

33. Урусов B.C., Петрова Т.Е., Еремин Н.Н. (2008) Компьютерное моделирование локальной структуры, свойств смешения и стабильности твердых растворов оксидов щелочноземельных металлов. Кристаллография, Т.53 №>6, стр. 1086- 1095.

34. Еремин Н.Н., Урусов B.C. (2009) Разработка и усовершенствование методов атомистического компьютерного моделирования твердых растворов замещения // В сборнике «Проблемы кристаллологии». Вып. 6. ГЕОС, Москва, стр. 29-83.

35. Еремин Н.Н., Громалова НА., Урусов B.C. (2009) Атомистическое моделирование и предсказание структуры, энергетики точечных дефектов, термодинамических и упругих свойств простых и сложных оксидов бериллия // Физика и химия стекла (в печати).

Подписано в печать 05 августа 2009 г.

Формат 60x90/16

Объём 2,75 п.л.

Тираж 200 экз.

Заказ Л® 070809227

Оттиражировано на ризографе в ООО «УниверПринт» ИНН/КПП 77285 72912\772801001

Адрес: 119333, г. Москва, Университетский проспект, д. 6, кор. 3.

Тел. 740-76-47, 125-22-73.

Ьйр://\\л\л\'.ип1\'егрпп1.ги

Содержание диссертации, доктора химических наук, Еремин, Николай Николаевич

Введение

Глава 1. ОБЩАЯ МЕТОДОЛОГИЯ АТОМИСТИЧЕСКОГО МОДЕЛИРОВАНИЯ 16 СТРУКТУР И СВОЙСТВ НЕОРГАНИЧЕСКИХ КРИСТАЛЛОВ

1.1. Энергия межатомного взаимодействия в кристалле

1.2. Потенциалы межионного и межатомного взаимодействия

1.2.1. Кулоновское взаимодействие

1.2.2. Парные короткодействующие потенциалы

1.2.3. Учет дисперсионного взаимодействия и поляризуемости атомов

1.2.4. Трех- и четырех- частичные взаимодействия

1.3. Процедуры поиска и уточнения параметров потенциала

1.4. Статическое моделирование в ионном приближении

1.5. Статическое моделирование в частично-ионном приближении

1.6. Возможности метода для расчета и предсказания ряда физических свойств 51 кристаллов

1.6.1. Расчет упругих свойств кристалла

1.6.2. Скорость прохождения акустических колебаний 53 1.6.3 Статические и высокочастотные диэлектрические константы. 54 Показатели преломления

1.6.4. Пьезоэлектрические константы

1.6.5. Электростатический потенциал и градиент электрического поля

1.6.6. Термодинамические свойства кристалла и фононные спектры 59 1.7 Методика расчета точечных, одно- и двумерных дефектов в кристаллах, включая поверхности

1.7.1. Классификация дефектов структуры

1.7.2. Формализм Мотта-Литтлтона

1.7.3. Особенности расчета кристаллической поверхности, как 66 разновидности двумерного дефекта кристалла

1.8. Методики расчета твердых растворов

1.8.1 .Энергетические аспекты образования твердых растворов замещения

1.8.2. Локальная структура твердых растворов замещения

Г.8.3. Современные теоретические подходы к моделированию изоморфных 78 смесей

1.9. Современное программное обеспечение структурного моделирования

Глава 2. РАЗРАБОТКА И УСОВЕРШЕНСТВОВАНИЕ МЕТОДОВ 86 АТОМИСТИЧЕСКОГО КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ

2.1. Метод минимизации энергии атомизации, основанный на поиске 86 оптимальных зарядов на атомах

2.1.1. Разработка обоснованных методик учета энергии переноса заряда

2.1.2. Программный продукт VALUES

2.1.3. Общая схема проведения моделирования путем оптимизации энергии 95 атомизации соединения»

2.2. Разработка методик моделирования твердых растворов минералов»

2.2.1. Особенности оригинальной методики структурного моделирования- 97 твердых растворов замещения

2.2.2. Выбор сверхячейки и возможные способы учета порядка-беспорядка 100 внутри нее

2.2.3. Критерии выбора второй координационной сферы

2.2.4. Общая методология расчета свойств смешения!

2.2.5. Пакет компьютерных программ ODSS' для выбора и анализа 112 оптимальной атомной конфигурации в рамках заданной сверхячейки

Программа Binar

Анионная'координация. Возможности задания различных зарядов, на 117 атомах

Программа Relax

Программа Gistogramma

Программа Volume

Глава 3. УТОЧНЕНИЕ ПАРАМЕТРОВ МЕЖАТОМНЫХ ПОТЕНЦИАЛОВ, 126 СТРУКТУРНОЕ МОДЕЛИРОВАНИЕ, РАСЧЕТ И* ПРЕДСКАЗАНИЕ СВОЙСТВ МИНЕРАЛОВ И НЕОРГАНИЧЕСКИХ КРИСТАЛЛОВ РАЗЛИЧНЫХ КЛАССОВ

3.1 Структурное моделирование в ионном приближении. Трансферабельность 126 параметров потенциалов, предсказательные возможности

3.1.1. Структурное моделирование купратов щелочноземельных металлов

3.1.2. Компьютерное моделирование кристаллической структуры, свойств 139 и устойчивости полиморфных модификаций'AbSiOj

3.2. Структурное моделирование в частично-ионном приближении. 157 Трансферабельность разработанных моделей потенциалов, предсказательные возможности приближения

3.2.1. Моделирование фосфатов структурного типа насикона

3.2.2. Межатомные потенциалы и структурное моделирование 163 дигидрофосфата калия KDP

3.2.3. Особенности разработки модели межатомных потенциалов и 167 структурного моделирования соединений, содержащих органические комплексы на примере бифталата калия КС8Н5О

3.3. Зарядовая оптимизация в структурном моделировании методом минимизации энергии'атомизации

3.3.1. Минимизация энергии атомизации кварца, корунда и стишовита

3.3.2. Минимизация энергии атомизации кианита 181 3.3.3: Минимизация энергии атомизации малаяита CaSn0Si04 184 3.3.4. Минимизация энергии атомизации оксидов переходных металлов

3.4. Особенности разработки согласованных потенциальных наборовv для 196 изострукгурных соединений! с целью последующего изучения их твердых растворов

3.4.1. Разработка согласованной модели потенциалов для полуторных 197 оксидов М2О3 со структурой*корунда

3.4.2. Разработка согласованного набора межатомных потенциалов для 200 структурного моделирования оксидов и силикатов бериллия

3.4.3. Разработка модели- потенциалов и моделирование структуры и 210 свойств оксидов группы рутила

3.4.4. Разработка согласованной модели для оксидов щелочноземельных 217 металлов со структурой галита

3.4.5. Моделирование структуры и свойств галита и сильвина

Глава 4. ПРЕДСКАЗАНИЕ ФОРМ ВХОЖДЕНИЯ И ЛОКАЛИЗАЦИЯ 222 СТРУКТУРНЫХ ПОЗИЦИЙ ПРИМЕСНЫХ И СОБСТВЕННЫХ ДЕФЕКТОВ

4.1. Локализация позиций примесей железа в структуре касситерита S11O

4.2. Моделирование дефектных центров в структуре ильменита

4.3. Моделирование дефектных центров, образованных двух- и 231 трехвалентными металлами в структуре KDP

4.4. Расчет энергии различных точечных дефектов в оксидах МО2 (M=Ti, Sn) 238 со структурой рутила.

4.5. Расчет энергии различных точечных дефектов в полуторных оксидах М2О3 242 со структурой корунда

4.6. Расчет энергии различных точечных дефектов в хризобериллах и бромеллите

Глава 5. КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ СВОЙСТВ СМЕШЕНИЯ В 247 БИНАРНЫХ И ТРОЙНЫХ СИСТЕМАХ ТВЕРДЫХ РАСТВОРОВ ЗАМЕЩЕНИЯ

5.1. Современные подходы к моделированию свойств смешения твердых 247 растворов замещения

5.2. Компьютерное моделирование свойств смешения и стабильности твердых 249 растворов оксидов щелочноземельных металлов

5.3. Моделирование свойств смешения и стабильности твердого раствора галит 256 NaCl - сильвин КС

5.4. Расчет свойств смешения и стабильности твердых растворов замещения в 261 системе корунд-эсколаит-гематит

5.5. Моделирование свойств смешения в системе ВеА1г04-ВеСг204 с учетом 269 распределения катионов А1(Сг) по различным октаэдрическим позициям

Глава1 6. АНАЛИЗ ЛОКАЛЬНОЙ СТРУКТУРЫ' ТВЕРДЫХ РАСТВОРОВ 277 ЗАМЕЩЕНИЯ ПО ДАННЫМ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ

6.1. Современное состояние вопроса изучения локальной^ структуры твердых 277 растворов*

6.2. Локальная структура твердых растворов Mgo jCao.jO и Cao.5Sro.5O.

6.3. Анализ локальной структуры твердого раствора галит-сильвин

6.4 Анализ локальной структуры эквимолярных твердых растворов 284 (All о,Сп о)03, (Fei о,Сп.о)Оз, (Ali.0,Fei.0)O

6.5 Анализ локальной структуры эквимолярных твердых растворов на основе 291 хризоберилла

Введение Диссертация по наукам о земле, на тему "Атомистическое компьютерное моделирование кристаллической структуры и свойств минералов, их дефектов и твердых растворов"

Актуальность проблемы: Со времени оформления себя в качестве науки в двадцатые годы прошлого века кристаллохимия стала вырабатывать собственные методы предвидения и предсказания кристаллических структур. Хорошим примером являются известные пять правил строения существенно ионных кристаллов знаменитого американского кристаллохимика Лайнуса Полинга [Полинг JI. (1947)]. Предложенные им еще в 1929 году правила до сих пор сохранили свое значение в качестве первичных рецептов построения кристаллических структур, к которым в настоящее время можно отнести методы начального предсказания топологии (геометрической схемы) структуры, а также способы уточнения этой схемы на основании некоторых критериев К первой группе относится серия приемов, которые позволяют конструировать основные геометрические черты строения' кристаллического вещества. Прежде всего, надо напомнить, что главным законом такого строения* является подчинение геометрического узора распределения атомов в кристаллическом пространстве одной из 230 пространственных групп симметрии. Другой важный принцип, который ограничивает возможное многообразие кристаллических структур, заключается в стремлении атомов или-молекул заполнить пространство не только наиболее симметрично, но и наиболее плотно. Поэтому значительную роль в описании^ и топологическом конструировании кристаллических структур, особенно интерметаллических и ионных соединений, играет теория к плотнейших упаковок. Второе правило Полинга позволяет наиболее рационально распределить в структуре атомы разной валентности. Оно требует, чтобы в стабильной структуре соблюдался локальный баланс валентностей, то есть сумма валентных усилий (отношение валентности атома к его координационному числу) всех ближайших соседей данного атома стремилась быть равной его собственной валентности Это правило делает маловероятным заполнение позиций с большим отклонением от локального баланса (больше 10-20%). Достаточно подробно эта группа методов описана в [Урусов и др, (1990)].

После того как пробная топологическая схема строения кристалла заданного состава сконструирована, возникает задача уточнения этой общей схемы путем определения размеров, отдельных деталей и элементарной ячейки кристалла Это может быть достигнуто либо с помощью процедуры DLS> (Distance-Least-Squares) - оптимизации пробной, геометрии структуры с помощью заданного оптимального набора межатомных расстояний [Meier W.M. and Villiger Н. (1969), Baur W.H. (1977)], либо с помощью метода валентностей связи - современной математической интерпретации модифицированного второго правила Полинга [Brown I. D. (1977), Урусов B.C. Орлов И.П. (1999)]. Отметим, что даже эти, самые усовершенствованные из геометрических методов конструирования кристаллических структур, используют весьма приближенные критерии устойчивости атомных конфигураций.

Два последних десятилетия прошлого века и первое — нынешнего ознаменовались гигантским прогрессом в области непосредственного теоретического предвидения структуры и свойств кристаллических веществ на основании единственно верного критерия нахождения устойчивого варианта кристаллической структуры: условия минимума структурной энергии кристаллического ансамбля, что связано с быстрым ростом вычислительных мощностей. Эта задача решается в принципе точно квантовохимически и приближенно - с помощью кристаллохимических атомистических полуэмпирических методов моделирования. Квантовая химия твердого тела пытается решить задачу нахождения стабильной атомной конфигурации заданного состава поиском минимума энергии взаимодействия всех электронов* и ядер в рассматриваемой системе. Основная принципиальная трудность при квантовохимическом моделировании давно известна - невозможность точного решения уравнения Шредингера с 1023 аргументами. Это обуславливает существование различных приближенных моделей, прогностические возможности которых на сегодняшний день еще ограничены [.ХаррисотУ'., (1983)]. Это связано с тем, что межатомные взаимодействия, ответственные за1 пространственное размещение атомов в элементарной ячейке, на несколько порядков слабее внутриатомных. Следовательно, для правильного предсказания энергии межатомного взаимодействия кристаллической структуры и функционально связанных с ней всех важнейших производных свойств кристалла (упругих, механических, диэлектрических, термодинамических и др.) квантовомеханические расчеты должны обладать огромной, и не всегда достижимой на сегодняшний день точностью. Пока удовлетворительная точность расчета достигается только для достаточно простых кристаллических систем (небольшое число атомов в элементарной ячейке, отсутствие элементов с большим атомным весом) и на настоящее время этот метод еще крайне «компьютероемкий». Однако, если вспомнить, с какой скоростью повышается быстродействие компьютеров, можно утверждать, что прямые методы расчета кристаллических структур окажутся на лидирующих позициях уже в самом ближайшем будущем.

На сегодняшний день использование кристаллохимических подходов способно, безусловно, принести успех там, где квантовохимические методы пока не в состоянии решить задачу поиска устойчивой атомной конфигурации. К ним относится структурное атомистическое моделирование методом межатомных потенциалов. Этот метод, ориентируясь на понятийный аппарат кристаллохимии, не требует рассмотрения внутриатомных взаимодействий, а ограничивается лишь межатомными, что сразу существенно облегчает вычислительную задачу поиска минимума структурной энергии системы (хотя также требует применения достаточно мощной вычислительной техники). В связи с этим, он стал быстро развиваться только в последние 30 - 35 лет, хотя кристаллохимия уже достаточно давно создала методику парных потенциалов, которые действуют между молекулами [Китайгородский А.И. (1971), Pertsin A.J., Kitaigorodsky A.I. (1986)] или ионами [Борн М, Хуан Кунъ (1958)]. В последние годы в кристаллохимической и минералогической литературе резко возросло количество публикаций, относящихся к методологии и результатам теоретического моделирования структуры и свойств минералов и неорганических веществ этим методом. Появилось много обзорных и монографических работ [Catlow C.R.A., Mackrodt W.C. Ed. (1982), Catti М. (1986), Burnham C.W. (1985), (1990), Урусов B.C., Дубровинский Л.С. (1989), Catlow C.RA., Price G.D. (1990), Catlow C.RA., Thomas J.M. et al (1993), Catlow C.R.A., Bell R.G., Gale J.D. (1994)], в которых с той или иной полнотой излагаются его принципы и детали. Разработка атомистических потенциалов' для основных типов, химического (металлического, ковалентного, ионного) и нехимического (ван-дер-ваальсова, индукционного, ориентационного) взаимодействий, позволило достаточно' корректно описать промежуточные по характеру типььхимической,связи!в кристаллах. Это, в свою очередь, позволило удачно воспроизвести структуры- и свойства, многих минералов и неорганических соединений, включая и относительно сложные, такие, как силикаты [Урусов B.C., Дубровинский JI.C. (1989)], в том числе цеолиты [Catlow C.RA., Thomas J.M., etal. (1993)].

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

Однако, нельзя утверждать, что в рамках атомистического подхода преодолены все трудности, и фраза, сказанная около двадцати лет назад редактором известного научного журнала «Nature» Дж.Мэддоксом [Maddox J. (1988)]: «Одним из продолжающихся скандалов в физических науках является то, что остается в общем невозможно предсказать структуру даже простейших кристаллических веществ из знания их химического состава» может быть окончательно опровергнута. К трудностям относятся недостаточная трансферабельность (переносимость) наборов межатомных потенциалов от одного класса соединений к другим, отсутствие строгих энергетических критериев моделирования с нецелочисленными зарядами атомов, недостаточная разработка результативных методик, для успешного моделирования^ реальных кристаллов, содержащих различные структурные несовершенства, вплоть до протяженных изоморфных- смесей. Все это определило основную цель работы: дальнейшее совершенствование метода атомистического моделирования структур и свойств неорганических кристаллов применительно к реальным, минеральным объектам, разработку новых алгоритмов и методик повышения** предсказательной, способности метода, внедрение новых, более физически обоснованных и результативных наборов межатомных потенциалов.

Научная новизна работы заключается, во-первых, в разработке оригинального и наиболее последовательного способа учета энергии переноса заряда от катиона к аниону для атомов в различных смешанных валентных конфигурациях. В работе продемонстрировано, что таким путем достигается заметное улучшение результатов компьютерного моделирования структуры и свойств неорганических кристаллов. Во-вторых, предложена оригинальная методика моделирования твердого раствора замещения, в рамках которой' используется предложенный автором новый критерий степени неупорядоченности расположения атомов в сверхячейке, что позволило оценить свойства смешения ряда бинарных и тройных твердых растворов и сравнить полученные результаты с имеющейся экспериментальной информацией. В-третьих, предложены и реализованы в оригинальных программах новые способы анализа локальной' структуры твердого раствора замещения, где помимо традиционных характеристик, таких как «релаксация» и «податливость катионной позиции» предложено использовать дополнительно еще предлагаемые автором характеристики: «объемную податливость» и «сдвигаемость атома». Это позволило провести более глубокий анализ локальной структуры в ряде изоморфных систем, где экспериментальная информация о ближнем атомном окружении ограничена или вообще отсутствует. В четвертых, разработанные в работе новые наборы межатомных потенциалов позволили в целом, обеспечить лучшее согласие рассчитанных и экспериментально наблюдаемых структурных, упругих, диэлектрических и термодинамических характеристик изучаемых природных и синтетических неорганических соединений, а также неизученных экспериментально гипотетических фаз, по сравнению с более ранними расчетами, где использовались известные наборы потенциалов, такие как библиотека Буша [Bush T.S. et al (1994)], Кэтлоу [iSchroder К.Р. et al (1992)] и Льюиса [Lewis G.V., Catlow C.RA. (1985)]. Это дает основание считать предсказание ряда неизвестных свойств этих соединений и модельные кристаллические структуры, неизученных фаз достаточно достоверными, а также рекомендовать разработанные наборы потенциалов межатомного взаимодействия для, моделирования других родственных соединений:

Необходимо'отметить еще тот факт, что за времяфаботы автора по теме настоящей' работы* произошли почти революционные изменения в быстродействии вычислительной техники: зашериод с. 1995 по 2008 г. быстродействие доступных автору вычислительных средств^ выросло как минимум наг два, порядка. Это; безусловно; наложило некоторый, «эволюционный» отпечаток на1, представленные, материалы диссертации: в конкретный, период времени-,, автор оперировал доступным на текущий момент программным обеспечением- и скоростью расчета; следовательно, вычисления, которые могли быть легко осуществлены непосредственно в период написания настоящей! работы, были невозможны по соображениям'быстродействия за несколько лет до этого. Следовательно,-многие, из полученных в ходе работы потенциалы, размерности сверхячеек, радиусы областей дефектов, и т.д. могут быть модернизированы с помощью компьютерной техники, доступной.на текущий момент.

Практическая значимость работы. Некоторые прикладные аспекты проведенного исследования состоят в следующем:

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

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

3) Созданный^ в ходе работы алгоритм нахождения максимально неупорядоченной атомной конфигурации произвольного состава для ячейки любого размера может быть рекомендован, как для ab-initio расчетов, так и для моделирования полуэмпирическими методами. Кроме того, этот подход может быть использован при разложении сложных экспериментальных спектров изоморфных смесей, которые с трудом поддаются интерпретации* без привлечения модельных представлений, о локальной структуре многокомпонентных систем.

4) Разработанные с участием автора программные продукты адаптированы для конечных пользователей, снабжены руководством пользователя* с подробным описанием, инструкцией по эксплуатации и примерами, файлов: информации, которые находятся в открытом доступе по* адресам: http://cryst.geol.msu.ru/odss/ и http://cryst.geol.rnsu.ru/values/. Программы,востребованы для расчета научными группами: как в нашей стране, так и за. рубежом.

5) Материалы диссертации* явились основой для. создания» нового спецкурса «Современные методы-моделирования твердых тел», читаемого автором на геологическом факультете МГУ, начиная: с 2007 года для специализаций^011304 (Кристаллография и кристаллохимия) по специальности 511010 (Геохимия).

6) Некоторые материалы, настоящей работы, используются автором в курсе «Компьютерное программное обеспечение в кристаллохимии» для студентов-3-его курса геохимического потока геологического факультета МГУ специализации1 011304 (Кристаллография и кристаллохимия).

Назащиту выносятся:

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

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

3. Систематические - расчеты энергий и позиций собственных точечных дефектов в минералах и неорганических соединениях разных классов.

Личный вклад автора: Все результаты, представленные в диссертации, получены самим автором, либо при его непосредственном участии, либо* под его руководством. Создание программного продукта VALUES- осуществлялось автором совместно со студентом Оболенским О.В. Создание пакета компьютерных программ ODSS, разработка, алгоритмов для расчета осуществлялось автором совместно с Урусовым B.C. (научные консультации) и< Деяновым- Р.З:- (программирование). Компьютерное моделирование структур-оксидов щелочноземельных металлов; и их твердых растворов осуществлялось совместно с Петровой): Т.Г., длительное время- работавшей! сначала под руководством автора, а затем и в качестве соавтора. Часть результатов, представленных в работе получена совместно со студентами и аспирантами-кафедры, кристаллографии Леоненко Е.В. (системы, NaCl-KCl и Ti02-Sn02) и Талисом Р.А. (системы, на основе корунда), защитившими, под руководством* автора свои- магистерские и бакалаврские работы. Изоморфные системы на основе хризоберилла исследовались автором совместно с аспиранткой кафедры кристаллографии Громаловой Н.А.

Апробация- работы и публикации: Материалы, включенные в диссертацию, докладывались на следующих всероссийских и международных конференциях: 16-ый конгресс Международного союза кристаллографов, Пекин, 1993; 5-ое Международное совещание по ядерно-спектроскопическим исследованиям сверхтонких взаимодействий. Дубна, 1993; 11-ое Международное кристаллографическое совещание «Sagamore», 1994; 2-ое Европейское совещание по резонансной абсорбции* и спектроскопии в минералогии, Берлин, 1995; Международное совещание АСА, Монреаль, 1995; 3-ее Европейское совещание по спектроскопическим, методам в минералогии, Киев, 1996; Международная конференция.^ «Спектроскопия, рентгенография и кристаллохимия* минералов», Казань, 1997; 1-ая Национальная кристаллохимическая конференция, Черноголовка, 1998; 10-ое Международное совещание кристаллографов, Страсбург, 1999; 18-ый Конгресс международного союза кристаллографов, Глазго, 1999; 19-ое Европейское совещание кристаллографов, Нанси, 2000; 19-ые Чтения имени академика Белова, Нижний Новгород,

2000; 2-ая Национальная конференция по кристаллохимии, Черноголовка, 2000; Международная конференция «Ломоносовские чтения», Москва, 2000, 2007, 2009; Конференция французского общества минералогов и кристаллографов, 2001; 15-ое Международное совещание «Рентгенография? и кристаллохимия минералов», Санкт-Петербург, 2003; 16-ая Международная конференция по химии фосфора, Бирмингем, 2004; Всероссийская конференция «Химия твердого тела и функциональные материалы-2004», Екатеринбург, 2004; Международная, конференция «Спектроскопия, рентгенография и кристаллохимия, минералов», Казань, 2005; 20-ый Конгресс международного союза кристаллографов, Флоренция, 2005; 4-ая Национальная кристаллохимическая конференция, Черноголовка, 2006; Международная конференция «Спектроскопия и кристаллохимия минералов», Екатеринбург, 2007; 16-ое Международное совещание «Кристаллохимия и рентгенография минералов», Миасс, 2007; 27-ые Научные чтения- имени академика Н.В.Белова, Нижний' Новгород, 2008; 4-ая, Всероссийская конференция «Математические исследования, в кристаллографии, минералогии и петрографии», Апатиты, 2008; Федоровская- сессия РМО; Санкт-Петербург, 2008; Всероссийская конференция^ «Минералы: строение,* свойства, методы исследования», Миасс, 2009.

Результаты, исследований, изложенные в диссертации,» отражены в 82 публикациях: 31 статье в отечественных и международных журналах, 3 обзорных* статьях в сборниках, 3 препринтах и 45* тезисах докладов, а также в отчетах по f проектам РФФИ, ИНТАС, грантам Президента РФ «Ведущие научные школы». Основных работ: 35.

Структура и объем работы. Диссертация состоит из введения, шести глав, заключения, выводов, списка цитированной литературы и приложений. Общий объем работы — 352 страницы, включая 74 таблицы (плюс 5 таблиц в приложении), 114 рисунков (плюс 1 рисунок в приложении) и список литературы из 444 наименований.

Заключение Диссертация по теме "Минералогия, кристаллография", Еремин, Николай Николаевич

основные вывода

299 оптимальных эффективных зарядов атомов в кристаллах и находить оптимальные значения ряда физических свойств кристалла

3. Использование частично ковалентной» модели при моделировании точечных дефектов по процедуре Мотга-Литглтона позволило сделать более обоснованные выводы о позиционном предпочтении, и механизме вхождения примесей в кристаллические структуры ряда неорганических соединений.

4. Разработанные автором оптимизированные наборы потенциалов были успешно применены для структурного моделирования и предсказания^ физических свойств купратов, содержащих С1Ю4 квадратные группировки, фосфатов структурного типа насикона, оксидов и силикатов, содержащих катионы Al, Fe, и Сг в октаэдрическом окружении а также минералов, содержащих Ве04 — тетраэдры.

заключение

В диссертационной работе проведен систематический обзор метода атомистического моделирования кристаллической структуры и свойств минералов и неорганических соединений, показаны его достоинства в сравнении с другими способами предсказания атомной структуры и свойств кристаллов и обозначены границы применимости изучаемого подхода. Отметим, что метод базируется на крисггаллохимической основе. Такие понятия кристаллохимии, относящиеся к эффективным характеристикам атомов и химических связей, как заряды атомов, эффективные радиусы атомов, степень ионности связи, несмотря на отсутствие физической строгости, оказываются чрезвычайно полезными не только для.описания-эмпирического материала, но и для.количественного прогнозирования. В- работе показано, что в рамках атомистического моделирования, возможно с использованием сравнительно простых полуэмпирических потенциалов межатомного взаимодействия предсказать основные упругие, структурные и термодинамические характеристики даже относительно сложных кристаллов. Структурное моделирование не ограничено в настоящее время описанием структуры и свойств совершенных кристаллов*- оно все больше используется для расчета структуры и свойств реальных кристаллов, содержащих дефекты в объеме и на поверхности, для моделирования локальной структуры твердых растворов' и оценки их стабильности. Необходимо отметить и то, что в рамках этого метода еще есть существенные резервы. Во-первых, это постоянно' растущие компьютерные возможности. Во-вторых -продолжение теоретических исследований в области разработки более реалистичных способов описания межатомных взаимодействий, в частности, путем включения в модели результатов ab-initio расчетов. Продолжение работ в этом направлении, без сомнения, повысит предсказательную способность метода и позволит еще с большим успехом применять его в качестве возможной альтернативы экспериментальных методов исследования кристаллов, путем параллельного развития и сближения возможностей атомистического подхода и методов «из первых принципов».

Библиография Диссертация по наукам о земле, доктора химических наук, Еремин, Николай Николаевич, Москва

1. Агошков В.М. (1985) Расчет термодинамических свойств минералов по модели колебательного спектра кристаллической решетки силикатов и окислов // Геохимия Т. 10, стр. 1441-1454.

2. Альберт А. (1977) Регрессия, псевдоинверсия и рекуррентное оценивание // М.: Наука, 223 стр.

3. Аракчеева А.В., Шамрай В.Ф., Леонюк Л.И., Лубман Г.У. (1994) Упорядочение атомов Са и Sr катионов в структуре (Cao5Sr05)(Cuo.75Bio.25)02 //Кристаллография, Т.39, стр. 928-931.

4. Асланов ЛА. (1989) Структуры веществ //М.: Изд. МГУ, 160 стр.

5. Белоконева ЕЛ. Щербакова Ю.К. (2003) Электронная плотность в синтетическом эсколаите СГ2О3 со структурой корунда и ее связь с антиферромагнитными свойствами // Журнал неорг. химии, Т.48, №6, стр. 960-968.

6. Беляев Р.А. (1980) Окись бериллия. М.: Атомиздат, 224 стр.

7. Берн Фр. (1969) Физические константы горных пород // М.: Мир, Гл.7.

8. Бетгер X. (1986) Принципы динамической теории решеток // М. Мир, 392 стр.

9. Бетехтин А.Г. (1951) Курс минералогии // М.: Государственное Издательство геологической литературы, 542 стр.

10. Борисанова Л.М., Засурская Л.А. (1990) Расчет пустот в кристаллической структуре ZSM-5 пентасила//Вестник МГУ. Сер. 2. Химия, Т.2, стр. 454-457.

11. Борисов С.В., Клевцова Р.Ф., Магарилл С.А., Первухина Н.В., Подберезская Н.В. (2002) Экспериментальная кристаллография от атомной до супрамолекулярной // Журнал структурной химии, Т.43, №4, стр. 664-670.

12. Борн М., Хуан К. (1958) Динамическая теория кристаллических решеток // М.: ИЛ, 488 стр.

13. Брэгг У.Л., Кларингбулл Г. Кристаллическая структура минералов // М.: Мир, 1967, 390 стр.

14. Вайполин А.А. (1990) Атомные смещения и сверхструктура точечных дефектов в тетраэдрических структурах полупроводников // Физика твердого тела Т.32, №12, стр. 3637-3641.

15. Веснин Ю.И., Заковряшин С.И (1979) О распаде твердых растворов NaCl КС1 // Solid State Communications. Vol.31, pp. 635-639.

16. Вотяков С.Л., Щапова Ю.В. (2008) Компьютерное моделирование атомной и электронной структуры дефектов в ряде оксидных и силикатных минералов. В17

Информация о работе

Еремин, Николай Николаевич
доктора химических наук
Москва, 2009
ВАК 25.00.05

Диссертация

Атомистическое компьютерное моделирование кристаллической структуры и свойств минералов, их дефектов и твердых растворов - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Атомистическое компьютерное моделирование кристаллической структуры и свойств минералов, их дефектов и твердых растворов - тема автореферата по наукам о земле, скачайте бесплатно автореферат диссертации

Похожие работы