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

ВВЕДЕНИЕ.

1. КОМПЬЮТЕРНЫЙ ЭКСПЕРИМЕНТ В ФИЗИЧЕСКОЙ ХИМИИ.

1.1. Постановка вычислительного эксперимента.

1.2. Метод Кара-Паринелло (квантовая молекулярная динамика).

1.3. Кластерные методы квантовой химии.

1.3.1. Основные модельные приближения.

1.3.2. Неэмпирические методы.

1.3.3. Полуэмпирические методы.

1.4. Метод Монте-Карло и его модификации.

1.5. Метод молекулярной динамики.

1.6. Модели оксидных систем с разными потенциалами.

1.7. Статистико-геометрические методы.;.

1.8. Применение полимерной модели к исследованию оксидов.

1.9 Выводы.

2. МАТЕМАТИЧЕСКИЕ МОДЕЛИ СТРУКТУРЫ И СВОЙСТВ ОКИДНЫХ РАСПЛАВОВ.

2.1. Комплексная модель оксидного расплава.

2.2. Имитационная молекулярно-динамическая модель.

2.3.' Ионная и ионно-ковалентная модели.'.

2.4. Параметризация потенциальных функций.

2.5. Модель структуры.

2.5.1. Модель ближнего порядка.

2.5.2. Модель наноструктуры.

2.6. Модель процесса полимеризации.

2.7. Модель термодинамических свойств расплава.

2.8. Выводы. з ■ •

3. ПОДСИСТ ЕМА СТАТИСТИКО-ГЕОМЕТРИЧЕСКОГО • 64 МОДЕЛИРОВАНИЯ «SGR» ИИС "SLAGJMELT".:.

3.1. Постановка задачи создания ИИС.

3.2. Архитектура и состав системы.

3.3. Интеграция подсистемы «SGR» в ИИС.

3.3.1. Общая архитектура подсистемы SGR.

3.3.2. Использование XML-технологии для обмена данными.

3.3.3. Применение CORBA-технологии для управления работой SGR-приложения.

3.3.4. Структура и документация классов подсистемы SGR.

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

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

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

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

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

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

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

В связи с вышеизложенным, крайне актуальной является задача разработки математических моделей структуры и их компьютерной реализации, создания в рамках ИИС подсистемы статистико-геометрическое моделирования, позволяющей получать широкий набор структурных характеристик и структурочувствительных свойств для сильно взаимодействующих полимеризующихся систем, к которым относятся исследуемые объекты. С привлечением статистико-геометрических (СГ)" методов можно исследовать наиболее актуальную для оксидных расплавов проблему - фундаментальную проблему корреляций структура — свойство: определять параметры структурообразования, особенности ближней и протяженной структуры, выделять характерные структурные элементы и закономерности их взаимного расположения, .изучать процессы полимеризации/деполимеризации при смешении оксидов сеткообразователей и оксидов - модификаторов, а с привлечением полимерной теории можно рассчитать ряд практически важных термодинамических параметров. i

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

Поставленная цель включает следующие задачи:

• разработка двухуровневой математической модели для компьютерного эксперимента по исследованию структуры полимеризующихся оксидных расплавов;

• разработка математической модели процессов полимеризации в оксидных расплавах;

• разработка модели расчета термодинамических свойств оксидных расплавов, с привлечением подходов полимерной теории

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

• проектирование и реализация на базе новых информационных технологий подсистемы статистико-геометрического моделирования «SGR» и ее интеграция в ИИС «Slag Melt»

• проведение серии комплексных компьютерных экспериментов по статистико-геометрическому моделированию оксидных систем

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

На защиту выносятся следующие результаты:

1. модель структуры полимеризующихся оксидных расплавов;

2. модель процесса полимеризации в оксидных расплавах;

3. модель термодинамических свойств оксидных расплавов;

4. пакет программ для статистико-геометрического моделирования многокомпонентных оксидных систем большой размерности;

5. подсистема статистико-геометрического моделирования «СГМ» распределенной информационно-исследовательской системы «Шлаковые расплавы»;

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

Научную новизну работы определяют впервые полученные, перечисленные ниже результаты.

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

• Разработана модель термодинамических свойств полимеризующихся оксидных расплавов;

• разработана математическая модель процессов полимеризации;

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

• спроектирована и реализована на базе новых информационных технологий подсистема статистико-геометрического моделирования «SGR» и осуществлена ее интеграция в распределенную информационно-исследовательскую систему с удаленным доступом «Slag Melt»; •

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

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

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

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

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

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

• ' 10 л "

4.6. Выводы

1. С целью проверки • адекватности разработанных моделей * и тестирования подсистемы ЗвЯ для объектов исследования проведены серии комплексных компьютерных экспериментов по статистико-геометрическому моделированию в рамках ИИС.

2. В качестве объектов исследования выбраны четыре бинарных оксидных системы БЮг-НагО, БЮг-КгО, БЮг-СаО, БЮг-М^О, которые имеют большое практическое значение в металлургической и стекольной промышленности.

3. Для всех систем рассчитаны параметры потенциальных функций в рамках ионно-ковалентной модели.

4. В результате статистико-геометрического эксперимента для всех расплавов из заданных диапазонов каждой из систем:

• определены структурные характеристики ближнего порядка;

• проведено исследование особенностей наноструктуры с помощью выделения кремнекислородных комплексов различной сложности;

• исследованы процессы полимеризации с помощью функций распределения полианионных комплексов;

• получены термодинамические характеристики полимеризующихся оксидных расплавов

5. Получены и проанализированы концентрационные зависимости состав-структура-свойство для изучаемых бинарных оксидных систем;

6. Проведено сравнение модельных данных с экспериментом. Сравнение модельных результатов с экспериментальными данными показывает их удовлетворительное согласие;

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

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

ЗАКЛЮЧЕНИЕ

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

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

- модель процессов полимеризации в оксидных расплавах.

- модель термодинамических свойств оксидных расплавов.

• Разработан пакет программ статистико-геометрического моделирования многокомпонентных оксидных систем большой размерности;

• Спроектирована подсистема статистико-геометрического моделирования «СГМ» и осуществлена ее интеграция в информационно-исследовательскую систему (ИИС) «Шлаковые расплавы».

7. Для четырех бинарных оксидных систем 8Ю2-СаО, БЮг-МдО, 8Ю2-№20, 8Ю2-К20 в широком диапазоне составов проведены серии комплексных компьютерных экспериментов по статистико-геометрическому моделированию. Получены и проанализированы концентрационные зависимости состав-структура-свойство для бинарных оксидных систем. Расхождение между модельными и данными натурного эксперимента составляет не более 10%.

8. Определены:

- параметры ближнего порядка,

- функции распределения полианионных комплексов,

- конфигурационное время жизни комплексов,

- доли кислорода различного типа.

9. Получены термодинамические характеристики полимеризованных расплавов (степень и константа полимеризации, активности сеткообразующего и модифицирующего оксидов).

1. Аппен A.A. Химия стекла,- JI.: Химия, 1974.-351 с.

2. Атлас шлаков. Справ.изд. Пер.с нем./Под редакцией Куликова И.С.,

3. М:Металлургия, 1985, 208с.

4. Байдов В. В. Ультраакустические исследования и микроструктурасиликатных расплавов.- В кн.: Свойства и структура шлаковых расплавов,- М.: Наука, 1970, с. 23-38.

5. Белащенко Д.К. Компьютерное моделирование структуры и свойствнекристаллических оксидов. Успехи химии 66(9), 1997, с.811-844

6. Белащенко Д.К., Федько А.Д. Компьютерное моделирование структуры исвойств оксидов системы Na20-Si02. Расплавы. 1992. N 5. с.40-50

7. Бернал Дж., Кинг С. Экспериментальное моделирование простыхжидкостей. М.: Мир, 1971. с. 116-135.

8. Близнюк A.A., Войлок A.A. MNDO-параметры для атома Ca

9. Журн.структ.хим., 1988, т.29, №15, с.156-157.

10. Близнюк A.A., Войлок A.A. Комплекс программ MNDO-85 для расчетаэлектронной структуры, физико-химических свойств и реакционной способности молекулярных систем полуэмпирическими методами MNDO,MNDOChAM1. Журн. структ. химии, 1986, 27, N4, с. 190-191.

11. Близнюк A.A., Войтюк A.A. Программа поиска полуэмпирическихпараметров метода МПДП- Ж.структ.хим., 1987, т.28, №2, с. 172-173.

12. Бобылев И.Б., Анфилогов В.Н. Особенности кристаллизации силикатных расплавов и расчет кривых ликвидуса в бинарных системах.- В кн.: Исследование структуры магматических расплавов,- Свердловск: УНЦ АН СССР, 1981, с. 52-61.

13. Бухтояров О.И. Исследование оксидных расплавов методом Монте-Карло / Известия АН СССР. Металлы, 1991, № 4, с. 124-129.

14. Бухтояров О.И., Воронцов Б.С., Комогорова С.Г. Исследование структурных особенностей системы Mg0-Si02-P205 методом Монте

15. Карло: Тезисы докладов .XI конференции по физической химии и электрохимии расплавленных и твердых электролитов. Екатеринбург, 1998, - с.12-13.

16. Бухтояров О.И., Курлов С.П., Ленинских Б.М. Прогнозирование структуры и термодинамических свойств расплавов системы Ca0-Si02 методом Монте-Карло.- Известия ВУЗов. Черная металлургия, 1985, № 11, с. 1-4.

17. Бухтояров О.И., Школьник Я., Смирнов Д., Курлов С.П. Расчет теплоты смешения и структурных группировок в расплавах системы Са0-А1203-Si02 методом Монте-Карло / Расплавы, 1987, т. 1, № 6, с. 45-49.

18. Ватолин H.A., Пастухов Э.А. Дифракционные исследования строения высокотемпературных расплавов. М.: Наука, 1980. - 189 с.

19. Войлок A.A. Применение метода MNDO для исследования свойств и реакционной способности молекул. Журн. структ. хим., 1988, 29, №1, с.138-162.

20. Воронов В.И. Методы исследования структуры сильновзаимодейст вующих систем. В сб. науч. трудов «Математическое и программное обеспечение научных исследований и обучения», Курган: изд-во Курганского гос.ун-та, 2000г., с.58-65.

21. Воронов В.И. Программа многомерной оптимизации базовых данных, полученных MNDO-методом. Тезисы докладов 3-го Российского семинара «Компьютерное моделирование физико-химических свойств стекол и расплавов». - Курган: КГУ, 1996, с.37-38.

22. Воронов В.И., Судников A.B. Об СГМ объекте информационно-исследовательской системы «Шлаковый расплав».' Сборник тезисов докладов научной конференции студентов и аспирантов Курганского Государственного Университета, Курган, 2001.

23. Воронова Л.И., Бухтояров О.И. Прогнозирование физико- химических свойств борного ангидрида методом молекулярной динамики.- Физика и химия стекла, 1987, т.13, N 6, с.818-823.

24. Воронова Л.И., Глубокий Я.В., Воронов В.И., Гроховецкий Р.В. Расчет са мосогласованного набора потенциальных параметров для МЖ)0-М0 мо делирования бинарных оксидных расплавов.-Расплавы, 1999, №2, с.66-74.

25. Гроссман И. Н. Моделирование структуры стеклообразной двуокиси кремния методом Монте-Карло // Физика и химия стекла, 1985, т. 11, № 6, с. 641-644.

26. Губанов В.А., Кумраев Э.З., Ивановский А.Л. Квантовая химия твердого тела. М.: Наука, 1984. -304с.

27. Денисов В.М., Белоусова Н.В., Истомин С.А. и др. Строение и свойства расплавленных оксидов. Екатеринбург: УрО РАН, 1999. - 499с.

28. Денисов В.М., Пастухов Э.А., Ченцов В.П. и др. Расплавы металлургии полупроводников: строение и физико-химические свойства. М.: Наука, 1991,320 с.

29. Евсеев A.M. Уравнения движения в квантовой молекулярной динамике. -Журн.физ.хим., т. LXII, № 4, 1988, с.972-977.

30. Есин O.A. О полимерной модели расплавленных силикатов и других окислов. Сталь, 1979, № 7, с. 497-500.

31. Есин O.A. Уравнения полимерной модели расплавленных силикатов в при ближении регулярных растворов. -ЖФХ, 1974, т. 48, вып. 8, с. 2108-2110.

32. Есин O.A., Гельд П.В. Физическая химия пирометаллургических процессов. М.: Металлургия, 1966, 704с.

33. Жидомиров З.М., Багатурьянц A.A., Абронин H.A. Прикладная квантовая химия.-М.:Химия, 1979.-295с.

34. Жоголев Д.А., Волков В.Б. Методы, алгоритмы и программы для кван тово-химических расчетов молекул. -Киев; Наукова думка, 1976, 212с.

35. Кларк Т. Компьютерная химия: Пер. с англ. М.: Мир, 1990, 384 с.

36. Климонтович Ю.Л. Кинетическая теория неидеального газа и неидеальной плазмы. М.: Наука, 1975. -352 стр.

37. Кожеуров В. А. Термодинамика металлургических шлаков. Свердловск: Металлургиздат, 1955,- 164 с.

38. Купряков Ю.П. Шлаки медеплавильного производства и их переработка. М.: Металлургия, 1987, 200с.

39. Малиновский В.К., Гудаев O.A., Гусев В.А. и др. Фотоиндуцированные явления в силленитах. Новосибирск: Наука, 1990. 160с.

40. Медведев H.H., Лучников В.А., Наберухин Ю.И. Области «совершенной» структуры в аморфном аргоне,- ЖСХ. 1994. т. 35(1). С. 53-63.

41. Момчев В.П. Компьютерное моделирование структуры и свойств некоторых жидких и аморфных металлов и сплавов на основе Fe, Ag и Ni диссертация на соискание ученой степени кандидата физико-математических наук, М.:, МИСИС, 1994, 152с.

42. Наберухин Ю.И. Что такое структура жидкости,- ЖСХ. 1981. т. 22(6). С. 62-80.

43. Новые приоритеты науки и техники. Документы для руководства (нормативная основа). - Москва, 1996, 27с.

44. Роусон Г. Неорганические стеклообразующие системы. М.: Мир, 1970, 302 с.

45. Слетер Дж. Методы самосогласованного поля для молекул и твердых тел. -М.:Мир, 1978, 356 с.

46. Сокольский В.Э., Казимиров В.П., Баталин Г.И. и др. Некоторые закономерности строения расплавов бинарных силикатных систем, составляющих основу сварочных шлаков,- Изв. вузов, Черная металлургия, 1986, №3, с. 4-9.

47. Темкин М.И. Смеси расплавленных солей как ионные растворы. ЖФХ, 1946, т. 20, вып. 1, с. 105-110.

48. Хокни Р. Методы расчета потенциала и их предложения,- В кн.: Вычислительные методы в физике плазмы, под ред. Б.Олдера, С. Фернбаха и М. Ротенберга.- М.: Мир, 1974, с. 143.

49. Хокни Р., Иствуд Дж. Численное моделирование методом частиц -М.:Мир, 1987,638 с.

50. Шенк Г. Физико-химия металлургических процессов. Ч. 2. Производство стали. К.: ГНТИУ, 193, 306с. .

51. Эверестов Р.А. Квантовохимические методы в теории твердого тела. -Л.:ЛГУ, 1982, 280 с.

52. Alder В. J., Wainwright Т. Е. Studies in Molecular Dynamics. I. General Method // J. Chem. Phys., 1959, Vol.31, N 2, p.459-466.

53. Amini M., Mitra S.K., Hockney R.W. J.Phys.C, Solid State Phys., 14,3689 (1981)

54. Angell C.A., Clarke J.H.R., Woodcock L.V. Interaction potentials and glass formation: a survey of computer experiments Adv.Chem.Phys., 1981, V.48, p.397-453.

55. Barker J.A., Henderson D. What is "liquid" Understanding the states of matter. -Rev.Mod.Phys., 1976, v.48, N 4, p.587-671.

56. Beeman, D. Some multistep methods for use in molecular dynamics calculations, J. Comput. Phys., 1976, vol. 20, pp. 130-139.

57. Bernal J.G. A geometrical approach to the structure of liquids.- Nature. 1959. V. 183 (4655). P. 141-147.

58. Binkley J.S., Whiteside R.A., Hariharan P.C., Seeger R., Pople J.A., Hehre W.J., Newton M.D. QCPE Program No.368.

59. Binkley J.S., Whiteside R.A., Raghavachari K., Seeger R., DeFrees D.J., Schlegel H.B., Frisch M.J., Pople J.A., Kahn L.R. GAUSSUAN82 Release A., Carnegie-Mellon University, Pittsburg, 1982.

60. Borgianni C., Cranati P. Thermodynamic propertiies of silicates and almnino-silicates from Monte-Carlo calculations.-Met. Trans. B, 1977 v. 88, No. 1, p. 147-151.

61. Borgianni C., Granati P. Monte-Carlo calculations of ionic structure in silicate and alumino-silicate melts.- Met. Trans. B, 1979 v. 108, No. 1, p. 21-25.

62. Brawer S.A., WeberM.J. Molecular dynamics simulations of the structure of rare-earth doped beryllium-fluoride glasses J.Chem.Phys., 1981, V.75, N 7, p.3522-3541.

63. Buneman O. Time-reversible difference procedures, J. Comput. Phys., vol. 1, pp. 517-537.

64. Car R., Parrinello M. Unified approach for molecular dynamics and density -functional theory. -Phys.Rev.Lett., v.55, 1985, p.2471-2474.

65. Catlow C.A., Freeman C.M., Islam M.S. and others Philos.Mag. A,58, 123 (1988)

66. Catlow C.A., Thomas J.M, Parker S.C. and Jefferson D.A. Nature 295,1982, p.658

67. Dewar M.J.S. The Molecular Theory of Organic Chemistry // McGraw-Hill, New York, 1969.

68. Dewar M.J.S., Thiel W. Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters. //J.Am.Chem.Soc., 1977, v.99, N 15, p.4899-4907.

69. Eastwood J.W., Hockney R.W., Lawrence D.N. P3M3DP the three-dimentional periodic particle-particle / particle-mesh program. -Comp.Phys.Comm., 1980, v.19, N.2, p.215-261.

70. Fumi, F. G., and Tosi M. P. Ionic Sizes and born repulsion parameters in the NaCl type alkali halides -1.- J. Phys. Chem. Solids, 1964, vol. 25, pp. 31-43.

71. Gaskell D. R. The thermodynamic properties of the Masson polimerization models in liquid silicates.- Met. Trans., 1973, v. 4, No. 1, p. 185-192.

72. Gaskell D.R. Thermodynamic models of liquid silicates.- Canad. Met. Quart., 1981, 20, N1, p. 3-19.

73. Gaskell P.H., Tarrant I.D. Refinement of random network model for vitreous silicon dioxide. Phil.Mag. B, 1980, Vol.42, No.2, p.265-286.

74. Hehre W.J., Lathman W.A., Ditchfield R., Newton M.D., Pople J.A, QCPE Program No.326.

75. Hehre W.J.,Radom L., Schleyer P., Pople J.A Ab Initio Molecular Orbital Theory//Wiley-Intersience, New-York, 1985.

76. Hirao K., Soga N. Molecular dynamics studies of structure and atomic motion of BeF2 glass // J.Non-Cryst.Solids, 1983, V.57, N.l, p.109-117.

77. Inoue H., Aoki N., Yasui N. J.Am.Ceram.Soc., 70,622 (1987)4 143

78. Kawamura K. In Molecular Dynamics Simulations. Springer Series in 'SolidState Sciences. Vol.103. (Ed.F.Yonezawa). Springer-Verlag, Berlin,1990, p.88

79. Keating P.N. Phys.Rev., 145,637 (1966)

80. Kubicki J.D., Lasaga A.C. Am.J.Sci., 292, 159 (1992)

81. Lasaga A.C. and Gibbs G.V. Phys. Chem. Miner. 14, 1987, p.107

82. Lewis J. W. E., Singer K., Woodcock L. V. Thermodynamic and structural properties of liquid ionoc salt obtained by Monte-Carlo computation.- J. Chem. Soc., Faraday 11,1975, v. 71, pp. 301-312.

83. Masson C. R., Smith I. B., Whitemay S.G. Molecular size distribution in muttichain polymers: application of polymer theory of to silicate melts. Can. J. Chem., 1970, v. 48, p. 201-202.

84. Masson C.R. Anionic constitution of glassforming melts.-J. Non-Cryst. Solids, 1977, 25, N 1, p.1-41.

85. Mitra S.K. Molecular dynamics simulation of silicon dioxide glass. Phyl.Mag., B, 1982, v.45, N 5, p.529-548.

86. Mulliken R.S. Magic formula, structure of bond energies and isovalent hybridization. J.Phys.Chem, 1952, v.56, !3, p.295-317.

87. Murray R.A., SongL.W., Ching W.Y. J.Non-Cryst.Solids, 94, 133 (1987)

88. Murrel J.N., Harget A.J. Semiempirical Self-Consistent-Field Molecular Orbital Theory of Molecules Wiley-Interscience, London, 1972.

89. Mysen B.O., Fingert L.W., Seifert F.A., Virdo D. Curve fitting off Raman spectra of amorphous materials.-American Miner., 1982, 67, p. 686-696.

90. Newton M.D., Gibbs G.V. Ab initio calculated geometries and charge distributions for H4Si04 and H6Si207 compared with experimental values for silicates and silioxanes // Phys.Chem.Minerals. 1980,- V.6.- p.221-246.

91. Ohno K., Morokuma K. Quantum Chemistry in Literature Data Base // Elsevier, Amsterdam, 1982; THEOCHEM, 1982, v.8, p.1983.

92. Pople J.A., Beveridge D.L., Approximate Moleculsr Orbital Theory // McGraw-Hill, New York, 1970.

93. Rahman A., Fowler R. H., Narten A.H. Structure and motion in liquid BeF2, LiBeF3 and LiF from molecular dynamics calculations.- J. Chem. Phys., 1972, * vol.57, pp. 3010-3011.

94. Richards W.G., Walker T.E.H., Hinkley R.K. A bibliography of Ab initio Molecular Wave Functions // Oxford University Press, Oxford, 1971.

95. Sangster M.J.L., Dixon M. Interionic potentials in alkali halides and their use similation of molten salts. Adv. Phys., 1976, v.25, N3, p.247 - 342.

96. Sangster M.J.L., Stoneham A.M. Philos.Mag. B, 43, 597 (1981)

97. Selloni A., Carnevali P., Car R., Parrinello M. Localization, hopping and diffusion of electrons in molten salts. Phys.rev.lett., v.59, 1987, p.823-826.

98. Soules T.F. and Arun K.Varshneya Molecular Dynamic Calculations of A Sodium Borosilicate Glass Structure. J.Amer.Ceram.Soc., v.64, N.3, 1981, p.145-150.

99. Soules T.F. Molecular dynamics calculations of glass structure and diffusion in glass J. Non-Cryst. Solids, 1982, V.49, N 1-3, p.29-52.

100. Stewart J.J.P. Optimization of parameters for semiempirical methods J.Comp.Chem., 1989, v.10, p.209-220.

101. Stoneham A.M., Sangster M.J.L. Philos.Mag. B, 43, 609 (1981)

102. Thiel W. The MNDOC Method, a Correlated Version of the MNDO Model -J.Am.Chem.Soc., 1981, v.103, p.1413-1420.

103. Toop G. W., Samis C. S. Activities of ions in silicate melts.- Trans. Met. Soc. AIME, 1962, v. 229, p. 878-887.

104. Tosi M.P., Fumi F.G. J.Phys.Chem.Solids, 25,31 (1964)

105. Vashishta P., Kalia R.K. Electron transport in disordered systems: A nonequilibrium quantum-molecular-dynamics approach. Phys.Rev.B, v.43, no. 13, 1991, p.10928-10932.

106. Verlet, L. Computer experiments on classical fluids. I. Thermodynamic properties of Lennard-Jones molecules.- Phys. Rev., 1967, vol. 159, pp 98-103.

107. Voronov V.I. The analysis of complex-forming process in molecular-dynamic simulation of oxide systems. Eight Asian Chemical Congress (8ACC), Taipei, Taiwan, 1999.

108. Voronova L.I., Voronov V.I., Gluboky J.V. The Oxide Melts Nanostructure Research by "Covalent Bonds Network Covering" Method. Papers of International Symposium on Cluster and Nanostructure Interfaces, Richmond, Virginia, USA, 1999

109. Waseda Y. The structure of non crystalline materials. Liquid and amorphous solids. New-York: Mc. Graws-Hill International book Co, 1980.-350 p.

110. Whiteside R.A., Frish M.J., Binkley J.S, DeFreees, Schlegel H.B., Raghavachari K., Pople J.A. Carnegie-Melton Quantum Chemistry Archive, 2nd Ed. Carnegie-Melton University, Pittsburg, 1981.

111. Woodcock L.V., Angell K.A., Cheeseman P. Molecular dynamics studies of the vitreous state: simple ionic system and silica.- J. Chem. Phys., 1976, v.65, N.4, p.1565-1577.

112. Woodcock L.V., Singer K. Trans.Far.Soc., 67, 12(1971)

113. Wooten F., Weaire D. Solid State Phys., 40,1(1987)