Вычисление точечных зарядов и VdW параметров для молекулярной механики
Файл et.top - описание молекулы в формате пакета программ GROMACS.
script
|
Жидкость |
|
Газ |
|
| Epsilon |
LJ(kJ/mol) |
Coulomb(kJ/mol) |
LJ(kJ/mol) |
Coulomb(kJ/mol) |
| 1.00000 |
-241.252 |
0.104627 |
-0.000226489 |
0.00563674 |
| 0.12500 |
-33.4987 |
0.03734 |
-0.000157829 |
0.00563675 |
| 0.03703 |
-2.17679 |
0.0167636 |
-0.000141172 |
0.00563674 |
| 0.01562 |
-2.28293 |
0.0148439 |
-0.000134252 |
0.00563673 |
| 0.00800 |
-2.52502 |
0.0145191 |
-0.000130619 |
0.00563673 |
| 0.00462 |
-2.73929 |
0.0144559 |
-0.000128428 |
0.00563673 |
| 0.00291 |
-1.07137 |
0.0141016 |
-0.000126999 |
0.00563673 |
Как видно из полученных данных, вклад кулоновских сил мал по сравнению с Ван-дер-Ваальсовыми.
Для газа значения энергий VdW взаимодействий все порядка 10-4,
тогда как у жидкости минимальный порядок 1.
Т.к. для обеих фаз кулоновские взаимодействия намного меньше, чем VdW,
то для оценки энтальпии испарения можно использовать значения только для VdW взаимодействий в жидкой фазе.
Чтобы воспроизводилась энтальпия испарения этана, epsilon водорода должна лежать в диапозоне от 0.01562 до 0.03703.
© Шерстюк Александра, MSU 2012