Оптимизируя результат прошлой работы была получена модель структуры et.pdb
Следуя инструкциям, был с нуля создан файл, описывающий молекулы этана - et.top
Используяскрипт и используя файлы жидкой и газовой фазы было получено следующие данные о среднем значении энергии для каждого значения эпсилон водорода
Для жидкой фазы:
epsilon | LJ | Coulomb |
1.00000 | -253.746 | 0.00479425 |
0.12500 | -31.1671 | 0.00120514 |
0.03703 | -5.85892 | 0.0065481 |
0.01562 | -3.37202 | 0.000458128 |
0.00800 | -3.14446 | 0.00056079 |
0.00462 | -1.78802 | 0.00042166 |
0.00291 | -1.96088 | 0.000483273 |
Видно, что влияние кулоновских сил на энергию крайне мало - это вызвано нейтральностью заряда этана. В то же самое время влияние Ван-дер-Ваальсовых сил значительно выше - и потому для оценки энтальпии испарения можно использовать эначения именно ВдВ взаимодействий.
Для газовой фазы:
epsilon | LJ | Coulomb |
1.00000 | -0.000226648 | 0.000157056 |
0.12500 | -0.000157936 | 0.00015705 |
0.03703 | -0.000141262 | 0.00015704 |
0.01562 | -0.000134335 | 0.000157045 |
0.00800 | -0.000130698 | 0.000157043 |
0.00462 | -0.000128505 | 0.000157045 |
0.00291 | -0.000127074 | 0.000157045 |
И вклад Ван-дер-Ваальсовых взаимодействий, и вклад кулоновских сил мал для оценки энтальпии. Лучше использовать все-таки значения VdW
При Т=25 значение энтальпии испарения этана равно 5.4 кДж/моль. Это значение лежит меж epsilon 0.03703 (Δ значения составляет 5.8588 кДж/моль) и epsilon 0.01562 (Δ значения составляет 3.3720 кДж/моль)