Расчет энергии для C-C связи в этане с шагом 0.02 ангстрема
Для генерации входных файлов использовался скрипт, в результате получился выходной файл
Значения коэффициентов для аппроксимации полиномом f(x)=a + k*x*x - 2*k*x*b + k*b*b:
Final set of parameters Asymptotic Standard Error
======================= ==========================
a = -79.7649 +/- 0.000385 (0.0004826%)
k = 0.532259 +/- 0.02172 (4.08%)
b = 1.55305 +/- 0.002164 (0.1394%)
Расчет энергии для C-C связи в этане с шагом 0.1 ангстрема
Для генерации входных файлов использовался скрипт, в результате получился выходной файл
Логично попробовать применить наши знания с курса "Математическая обработка экспериментальных данных" и применить
полиномиальную аппроксимацию. Известно, что чем выше степень полинома, тем лучше результат но и выше вычислительная сложность. К сожалению, gnuplot на моей машине (а возможно и вообще) не осиливает полиномы степени выше 6. Результат фита приведен ниже.
Final set of parameters Asymptotic Standard Error
======================= ==========================
a = 1.77201 +/- 0.1437 (8.111%)
b = -18.7336 +/- 1.321 (7.053%)
c = 81.4527 +/- 4.867 (5.975%)
d = -186.945 +/- 9.147 (4.893%)
e = 239.989 +/- 9.202 (3.834%)
f = -164.388 +/- 4.67 (2.841%)
g = -32.4989 +/- 0.9295 (2.86%)
Расчет энергии для валентного угла H-C-H в этане с шагом 0.2 градуса
Для генерации входных файлов использовался скрипт, в результате получился выходной файл
Значения коэффициентов для аппроксимации полиномом f(x)=a + k*x*x - 2*k*x*b + k*b*b:
Final set of parameters Asymptotic Standard Error
======================= ==========================
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 в этане с шагом 12 градусов
Для генерации входных файлов использовался скрипт, в результате получился выходной файл
У данной функции 3 минимума, так как угол -Pi соответствуют углу Pi.
|