Определение зависимости энергии молекулы этана от длины связи C-C
файл-заготовка, на основе которого будут происходить дальнейшие вычисления
Первый и второй скрипты. Первый предназначен для формирования inp и log файлов, второй же - для создания bond-файла, используемого в дальнейшем в gnuplot.
Не аппроксимированная зависимость энергии от длины связи. Говоря об аппроксимации f(x)=a+k*(x-b)^2; стартовые значения коэффициентов:
a = -80
k = 1
b = 1.5
Параметры аппроксимирующей функции:
a = -79,7652 +/- 0,0004522 (0,0006% ошибка)
k = 0,563608 +/- 0,02335 (4,142% ошибка)
b = 1,55432 +/- 0,002455 (0,16% ошибка)
График
Сильное отклонение в k - возможно, функция сложнее, чем кажется
Для валентного угла HCH: Значения углов; Скрипт
И снова парабола. Попробуем ту же функцию с теми же коэффициентами, кроме b, для него поставим стартовый параметр 111
Итоговые коэффициенты:
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
значение изменялись от -180 до 180. Тут лежат данные; Скрипт
Крайние точки, что логично, совпадают/
Назад