Вычисление параметров для валентных углов
Шестой семестр
Построение зависимости энергии от длины связи С-С в молекуле этана
1). Расчет энергии молекулы этана с помощью GAMESS при изменении длины связи С-С с шагом 0.2 Å.
+ - результат расчетов GAMESS,
зеленая линия - подобранная (средствами gnuplot) функция f(x) = a + k(x-b)^2 , описывающая зависимость.
исходные значения коэффициентов: a=-80, k=1, b=1.5.
оптимизированные значения коэффициентов: a=-79.198, b=1.5514, k=0.6170.
коэффициент k подгоняется хуже всего.
На следующем графике видно, что точки удовлетворительно ложатся на параболу только в узком диапазоне значений.
2). Расчет энергии молекулы этана при изменении длины связи С-С с шагом 0.1 Å.
+ - результат расчетов GAMESS,
красная линия - подобранная функция f(x)=t*exp(s*x)/x+r , описывающая зависимость.
исходные значения коэффициентов: t=-0.1, s=-1, r=-80.
оптимизированные значения коэффициентов: t=22.5302, s=22.5302, r=-79.1715.
коэффициент k подгоняется хуже всего.
Построение зависимости энергии от величины валентного угла H-C-H в молекуле этана
3). Расчет энергии молекулы этана при изменении величины валентного угла H-C-H с шагом 0.1°.
+ - результат расчетов GAMESS,
зеленая линия - подобранная функция f(x) = a + k(x-b)^2 , описывающая зависимость.
исходные значения коэффициентов: a=-80, k=1, b=1.5.
оптимизированные значения коэффициентов: a=-79.1976, k=4.06217e-05, b=111.195.
Построение зависимости энергии от величины торсионного угла H-C-H в молекуле этана
4). Расчет энергии молекулы этана при изменении величины торсионного угла d3 от 180° до -180° с шагом 12°.
+ - результат расчетов GAMESS,
Функция имеет минимумы в: 60°, 180° и -60° на отрезке (-180;180].
Возможные причины расхождений с результатами J. Wang et al
1. Оптимизированные длины связей вычислялись другими методами, поэтому они несколько отличаются от полученных нами (длина С-С - 1.526).2. Подгонка коэффициентов делалась не с помощью gnuplot, а в Parmscan.