Занятие 6
Подготовка
Проводим DFT оптимизацию, получаем идеально симметричную молекулу в заторможенной конформации.
После обработки Ante_RED - заряд 0, мультиплетность - 1, как и должно быть.
После обработки RED - почему-то ассиметричный этан получился, разные заряды на атомах.
Расчет epsilon для водорода
Файл топологии для системы
Скрипт
Файлы с результатами расчета:
Для жидкой фазы:
v_1.log,
v_2.log,
v_3.log,
v_4.log,
v_5.log,
v_6.log,
v_7.log,
Для газа:
vb_1.log,
vb_2.log,
vb_3.log,
vb_4.log,
vb_5.log,
vb_6.log,
vb_7.log.
Результат
Вклад В-д-В взаимодействий оказался намного больше, чем вклад кулоновских для жидкой фазы. Для газа вклад В-д-В также больше, чем кулоновский, но сам он значительно уменьшился и стал соизмерим с ним (отличие "всего" на пару порядков).
По результатам, epsilon должно лежать в интервале от 1/27 до 1/64.