Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

Молекулярное моделирование играет важную роль в развитии современной химии. Для каждого круга задач предлагались свои модели, как, например, шаро-стержневые, Стюарта-Бриглеба и т.д. С развитием новых компьютерных технологий возникли потребности в новых представлениях молекул в трехмерном пространстве. В частности, интересные результаты получены в результате анализа поверхности постоянной электронной плотности Р.Бейдером [1].

В настоящей работе предлагается новый, более простой и доступный химикам подход, базирующийся на модели DENSON [2]. Новизна его состоит в том что, в конечном итоге, атомы в молекуле не имеют сферической симметрии, в сравнении с более ранним подходом [3].

Для получения пространственной формы молекулы требуется знание декартовых координат атомов ее составляющих. Такую информацию можно получить, например, в рамках вычислительного комплекса "Hyper Chemistry" [4] посредством оптимизации геометрии молекулы одним из доступных методов, вычисляя полную энергию системы методом РМ3. Далее каждому типу атомов Периодической системы Д.И.Менделееева приписывается собственный атомный радиус (r0), являющийся параметром модели. В частности, нами предлагается в качестве r0 для атома водорода принять величину в 0,3325 ангстрем. Соответственно для атома углерода - r0 (С) =0,6500 . Предполагается, что в результате атом-атомного взаимодействия в молекуле сферические оболочки атомов должны деформироваться. Величина деформации (∆r0) для атомов, находящихся на валентных расстояниях, вычисляется по формуле (1), где Ri µ евклидово расстояние между атомами i и μ ; n - главное квантовое число атома μ. Для невалентных взаимодействий такая поправка имеет вид (2). Очевидно, что если рассматриваемый атом имеет несколько соседних атомов, с которыми есть взаимодействие, то исходная сферическая оболочка будет деформироваться в разных направлениях и на разную величину. Предлагается простой механизм трансформации исходной атомной оболочки.

Δri = rµ0 exp (Ri µ)/ Ri µ n-1                                                (1)

Δri = - r0μ exp (Rμ j)/ Rμ j 6 - r0μ exp (Rμ j)/ Rμ j 12                      (2)

В начале поверхность каждого из атомов молекулы разбивается на N- точек, координаты которых запоминаются. Далее для каждой из точек поверхности атома рассматриваются все возможные взаимодействия и строятся вектора различного направления и значения согласно вычисленным ∆r. В итоге имеем для каждой точки несколько векторов для которых потом находится результирующий. Координаты конца результирующего вектора запоминаются и определяют местоположение новой точки поверхности атома с учетом всех взаимодействий.

На рисунке видны результаты трансформации исходной атомной оболочки каждого из атомов молекулы на примере пропилена. На сколько новый молекулярный объем получил удовлетворительную оценку можно определить путем сравнения экспериментальных и вычисленных величин относительных плотностей веществ при 20оС (d420), как это и делалось ранее [2,3]. Отдельные результаты расчетов приведены в таблице.

p

p

    без учета атомных взаимодействий                                    с учетом атомных взаимодействий

Рисунок 1. Вид молекулы пропилена без учета атомных взаимодействийи с их учетом

Таблица 1. Экспериментальные d420(эксп.) [5] и вычисленные d420(расч.) относительные плотности некоторых углеводородов

Соединения

d420(эксп.)

d420(расч.)

Соединения

d420(эксп.)

d420(расч.)

Гексан

0,6548

0,6966

Циклобутан

0,703

0,712

Гептан

0,6838

0,6985

Циклопентан

0,7554

0,7105

Октан

0,7025

0,6999

Пентин-2

0,7127

0,7503

Нонан

0,7176

0,701

Гексадиен-1,5

0,688

0,742

Декан

0,7301

0,7019

2-метилбутен-2

0,668

0,718

Додекан

0,7437

0,7033

2-метилпентан

0,6599

0,6969

Тетрадекан

0,764

0,704

2,2-диметилбутан

0,6485

0,6973

Пентадекан

0,769

0,705

2,3-диметилбутадиен-1,3

0,745

0,741

Циклопропан

0,6886

0,7061

2,5-диметилгексадиен-2,4

0,716

0,731

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

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

СПИСОК ЛИТЕРАТУРЫ

  1. Бейдер Р. Атомы в молекулях: Квантовая теория.-М.:Мир, 2001.-532с.
  2. Белик А.В. Компьютерное моделирование атомных радиусов для прогноза плотности органических веществ //Изв. вузов. Сер. Химия и хим. технология.1992. Т. 35, №4. С.51-55.
  3. Белик А.В., Потемкин В.А. Новый подход к компьютерному моделированию атомных радиусов //ЖОХ. 1993. Т.63. Вып.6. С.1201-1203.
  4. http: //www.hyper.com.
  5. Свойства органических соединений. Справочник / Под ред. А. А. Потехина. Л.: Химия, 1984. 520 С.