Вычисление валентных параметров для молекулярной механики
Файл-заготовкаScript
В результате полулся 21 inp файл и в каждом разное значение для переменной сс.C помощью GAMESS получены значения энергий. Выходной файл содержит эти значения и соотвутсвующие длины связей
Зависимость строим в Gnuplot:
Значения коэффициентов:
a = -79.7652 +/- 0.0004522 (0.000567%) k = 0.563608 +/- 0.02335 (4.142%) b = 1.55432 +/- 0.002455 (0.1579%)Aналогичные операции для валентного угла HCH:
Заготовка
Script
Результат
Значения коэффициентов:
a = -79.7647 +/- 1.21e-08 (1.517e-08%) k = 3.56076e-05 +/- 6.229e-09 (0.01749%) b = 111.38 +/- 9.954e-05 (8.937e-05%)Торсионный угол d3:
Заготовка
Script
Результат
Значения коэффициентов:
a = 0.00234519 +/- 8.891e-07 (0.03791%) k = 3.00014 +/- 0.0002247 (0.007491%) b = -79.7623 +/- 6.577e-07 (8.245e-07%)