Здравствуйте,
Мне очень нужна помощь вот по какому вопросу:
Можно ли с помощью МД воспроизвести химпотенциалы для молекул в газовой фазе? По логике вещей, изменение потенциала от температуры будет разницей от энергии релаксированной молекулы и средней энергией термостатированной. Я попробовала сделать такое для углеводородов, но соотношение тангенсов угла наклона для нескольких молекул не соответствует эксперименту (сверялась с JANAF), по экспериментальным данным химпотенциал последовательно растет в ряду метан-ацетилен-этилен, у меня выходит другой порядок. Вопроса два:
1. Реально ли в принципе такое сделать
2. Какой ансамбль будет правильным? NVT или NPT?
Или, вероятно, я не учитываю еще какие-то важные параметры, которые нужны?
Заранее огромное спасибо
Химпотенциалы для молекул
-
- Сообщения: 2
- Зарегистрирован: Пт мар 10, 2023 12:30 pm
Re: Химпотенциалы для молекул
В полный химпотенциал основной вклад вносит внутренняя электронно-ядерная энергия молекулы, вычислить которую методами силового поля, по моему, не реально (и в любом случае, она будет зависеть не от термостата, а от параметров силового поля). Однако, эта внутренняя энергия, насколько я понимаю, не зависит от температуры. Поэтому с температурной зависимостью химпотенциала, наверно, работать можно. Мне кажется, что нужно сравнивать не просто темп. зависимости химпотенциалов отдельных соединений, а компонентов изодесмических реакций типа ацетилен + Н2 = этилен.Любовь Варламова писал(а): ↑Пн мар 13, 2023 4:37 pmхимпотенциал последовательно растет в ряду метан-ацетилен-этилен, у меня выходит другой порядок.
-
- Сообщения: 2
- Зарегистрирован: Пт мар 10, 2023 12:30 pm
Re: Химпотенциалы для молекул
Да, но, по сути, эта энергия - глобальный минимум, который можно найти из релаксации и принять за ноль, смотреть только на изменение.
Ради этих реакций и пишется. В конечном итоге это должны быть температурные зависимости изменения энергии Гиббса для присоединения водорода к довольно большим кластерам графена - в конечном итоге мне нужно найти как этот самый химпотенциал будет меняться для кластеров графена разного размера и хотя бы как-то оценить это изменение количественно
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 11 гостей