Вычисление валентных параметров для молекулярной механики

Длина С-С связи

Пусть длина С-С связи варьируется от 1.32986 до 1.72986 ангстрем с шагом 0,02 ангстрема.
Скрипт для вычисления 21 значения энергии.
Итоговые значения энергии.

Зависимость энергии от длины связи:
Подобранная кривая описывается уравнением f(x)=a + k*x*x - 2*k*x*b + k*b*b
Где

a = -79.7652

k = 0.563608

b = 1.55432
Кривая не везде совпадает со значениями энерии, так как квадратный трехчлен - симметричен относительно минимума, а энергия меняется по-разному при уменьшении/увеличении длины связи.

Валентный угол HCH

Пусть угол HCH варьируется от 109.2 до 113.2 градусов с шагом 0,2 градуса.
Скрипт для вычисления 21 значения энергии.
Итоговые значения энергии.

Зависимость энергии от длины связи:
Подобранная кривая описывается уравнением f(x)=a*x*x + b*x +c
Где


a = 3.56076e-05

b = -0.00793196

c = -79.3229


Торсионный угол d3

Пусть угол варьируется от -180 до 180 с шагом 12 градусов.
Скрипт для вычисления 31 значения энергии.
Итоговые значения энергии.

Зависимость энергии от длины связи:

Кривая похожа на синусоиду. Имеютя три точки минимума (60, -60 и 180 градусов).