Леоненко Е.В., Еремин Н.Н., Урусов В.С.
КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ СТРУКТУРЫ, СВОЙСТВ И ТОЧЕЧНЫХ ДЕФЕКТОВ РУТИЛА И КАССИТЕРИТА



РЕФЕРАТ

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

Методом атомистических парных потенциалов проведено компьютерное моделирование структур, упругих и термодинамических свойств рутила TiO2 и касситерита SnO2 при помощи программного комплекса GULP (GaleRohl, 2003), в основе которого лежит поиск минимума структурной энергии. Моделирование проводилось в рамках ионной и ионно-ковалентной моделей. Параметры потенциалов оптимизированы путем расчетов параметров и и объема элементарной ячейки V, упругих констант Cij и модуля упругости K, и энтропии ST чистых компонентов TiO2 и SnO2.

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

Выбор степени ионности рутила TiO2 и касситерита SnO2 в рамках ионно-ковалентной модели основывался на литературных данных по расчету зарядов ионов Ti и Sn, было выбрано значение = 0.7. Для учета ковалентного связывания был использован потенциал Морзе. Ионно-ковалентная модель хорошо описывает как структурные и термодинамические, так и упругие свойства рутила и касситерита.

В таблице сравниваются структурные и  упругие свойства TiO2 и SnO2, полученные экспериментально и вычисленные в данной работе в рамках ионно-ковалентной модели.

 

 

TiO2

SnO2

эксперимент

расчет

эксперимент

расчет

a, Å

4.5943

4.5993

4.737

4.7623

c, Å

2.9586

2.9857

3.187

3.1591

uостова

0.305

0.3048

0.306

0.3036

uоболочки

0.3084

0.3059

V, Å3

62.449

63.158

71.559

71.645

C11, ГПа

26.9

32.0

26.3

31.7

C 12, ГПа

17.7

10.8

17.8

16.6

C 13, ГПа

14.8

18.6

15.6

14.4

C 33, ГПа

48.2

56.4

45.0

49.2

C 44, ГПа

12.4

13.7

10.3

13.5

C 66, ГПа

19.3

15.9

20.8

17.6

K, ГПа

218.6

226.6

217.7

223.2

Eстр, кДж/моль

-60.770

-58.058

 

Моделирование точечных дефектов в рутиле и касситерите было проведено в рамках модели Мотта-Литтлтона как по ионной, так и по ионно-ковалентной моделям. Вычислены энергии образования вакансий и интерстиций анионов и катионов, примесного дефекта Ti в SnO2 и Sn в TiO2, тройных дефектов Шоттки и катионных и анионных дефектов Френкеля. Расчет по ионно-ковалентной модели лучше согласуется с экспериментальными данными, чем расчеты по ионной модели и «из первых принципов».

Литература

 Gale J.D., Rohl A.L.The General Utility Lattice Program // Mol. Simul., 2003. V. 29. №5. P. 291-341.