Часть полного текста документа:Построение 3D-моделей нециклических молекул в естественных переменных Е.Г. Атавин, В.О. Тихоненко, Омский государственный университет, кафедра органической химии 1. Введение По мере накопления химической информации роль данных о пространственном геометрическом строении молекул возрастает. Устанавливать его можно как экспериментальными, так и теоретическими методами, а описывать принято либо в декартовой системе координат, либо в естественных (внутренних, молекулярных) переменных. Первый способ предполагает знание 3N декартовых координат N атомов, позволяет легко строить графическое изображение молекулы, вычислять значения всех естественных переменных и используется в большинстве современных программ квантовой механики, молекулярной механики и колебательной спектроскопии. Однако произвол в выборе положения начала координат и ориентации координатных осей затрудняет сравнение результатов, полученных разными авторами. Кроме того, наличие у молекулы трех поступательных и трех вращательных степеней свободы приводит к появлению шести нулевых собственных значений у матрицы вторых производных энергии по координатам и к дополнительным осложнениям вычислительного характера [1]. Наконец, само задание декартовых координат атомов - нетривиальная задача, поскольку они не являются справочными данными. Описание (и анализ) геометрического строения в естественных переменных (ниже - межъядерные расстояния R, валентные углы ?, ?, ? и углы внутреннего вращения ?, F) проще, поскольку задание их не представляет проблемы и менее зависит от произвола исследователя, благодаря имеющимся эмпирическим закономерностям [2]. При оптимизации геометрии молекулы можно упрощать задачу, фиксируя значения хорошо известных параметров. Легко организовать поиск глобального минимума энергии путем перебора допустимых значений всех или некоторых параметров. При работе же с декартовыми координатами реализация этих возможностей сопряжена со значительными трудностями. Однако непосредственно по значениям естественных переменных невозможно в общем случае построить графическое изображение молекулы. Также затруднительно выполнять любые вычислительные операции с моделью молекулы, например, определять вандерваальсовые расстояния между атомами. Таким образом, оба способа описания молекулярной геометрии обладают рядом практически важных достоинств и весьма существенных недостатков. Совмещение достоинств достигается вычислением декартовых координат атомов по заданным естественным переменным, что представляет собой в общем случае весьма громоздкую стереометрическую задачу. Цель настояшей работы и состоит в рассмотрении алгоритмов вычисления декартовых координат атомов по заданным естественным переменным то есть построения 3D-моделей молекул. 2. Метод Эйринга [3] Обычно систему координат связывают с положением первых трех атомов (рис. а), координаты которых, таким образом, определяются по формулам: x1 = R12 cos?, x2 = 0, x3 = R23 y1 = R12 sin?, y2 = 0, y3 = 0 (1) z1 = 0, z2 = 0, z3 = 0 Легко также вычислить координаты четвертого атома: x4=x3-R34cos? y4=R34sin?cos? (2) z4=R34sin?sin? Рис. 1. Ориентация молекулы в системе координат Далее, построенный фрагмент с помощью переноса и двух поворотов переводится в положение, показанное на рис. ............ |