Вычислением параметров молекулярной механики.
1. Зависимость энергии молекулы от длины одной сс связи в этане
График зависимости энергии от длины связи cc
файл с результатами
коэффициенты аппроксимирующей функции f(x)=a + k*x*x - 2*k*x*b + k*b*b
Final set of parameters Asymptotic Standard Error
======================= ==========================
a = -79.7652 +/- 0.0004522 (0.000567%)
k = 0.563608 +/- 0.02335 (4.142%)
b = 1.55432 +/- 0.002455 (0.1579%)
на графике видна типичная ситуация для полиномов, когда на краях аппроксимация не очень точная, чтобы этого избежать можно брать точки на краях почаще, полином более высокой степени или использовать сплайн.
2. Зависимость энергии молекулы от длины одной нсн связи в этане
График зависимости энергии от длины связи нсн
файл с результатами
коэффициенты аппроксимирующей функции f(x)=a*x*x + b*x+c
Final set of parameters Asymptotic Standard Error
======================= ==========================
a = 3.56076e-005 +/- 6.227e-009 (0.01749%)
c = -79.3229 +/- 7.7e-005 (9.707e-005%)
b = -0.00793195 +/- 1.385e-006 (0.01746%)
Расхождения со статьей.
В оптимальном значении связи сс наши результаты практически сошлись у них 1.526, а у нас около 1.52986, расхождение может быть связано с различной шириной шага или обновлением программы, подсчитывающей энергию.
3. Зависимость энергии молекулы от значения торсионного угла в этане
График зависимости энергии от значения торсионного угла
файл с результатами
На графике видно 3 минимума, которые видимо соответствуют заторможенным конформациям этана.
© Garanina Irina