Подготовка файлов с координатами и топологией

Индекс-файл 1.ndx с группой из одной молекулы этана.
et.gro файл с одной молекулой и заданной ячейкой.
Исправленный файл топологии et.top из прошлого задания (изменено количество молекул).


Визуальный анализ
Для каждой из пяти систем проводим конвертацию в pdb-файл. Получаем
et_an.pdb - молекула неподвижна, вращение по СС отсутствует, немного колеблются СН связи.
et_be.pdb - сначала вращение по связям СН, потом и по СС, молекула вращается.
et_nh.pdb - вращение только по связям С-Н.
et_sd.pdb - молекула движется, вращение по С-С связям, углы с водородами не изменяются.
et_vr.pdb - сначала изменяются углы между углеродом и водородом, потом немного вращение по С-С связи.

Сравнение потенциальной и кинетической энергии и распределения длины С-С связи во время моделирования

Метод Зависимость энергии Зависимость длины С-С-связи
Андерсена
Берендсена
Нуза-Хувера
Стохастической мол.динамики
Velocity rescale


Из полученных зависимостей более всего похожи на распределение Больцмана те, что получены методами Velocity rescale и стохастической молекулярной динамики. Однако в Pymol Velocity rescale выглядит адекватнее, без постоянного хаотичного движения молекулы, поэтому лучшим для контроля температуры методом я бы выбрала именно его.

 

 

© Дудина Дарья.