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

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

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

Введение

Глава 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. Систематические - расчеты энергий и позиций собственных точечных дефектов в минералах и неорганических соединениях разных классов.

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

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

Работа выполнена на кафедре кристаллографии и кристаллохимии Геологического факультета МГУ им. М.В.Ломоносова. Отдельные этапы работы были поддержаны грантами 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 проекта А0109 "Науки о Земле". Научные заслуги автора также отмечались грантами Дж. Сороса (1996 г.) и фондом содействия отечественной науки (2004 г.).

Личный вклад автора: Все результаты, представленные в диссертации, получены самим автором, либо при его непосредственном участии, либо* под его руководством. Создание программного продукта 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 наименований.

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

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

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

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

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

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

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

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

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

заключение

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

Список литературы диссертационного исследования доктор химических наук Еремин, Николай Николаевич, 2009 год

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