Термохимия

[AlexKsen]

Друзья, подскажите, как в программе гауссиан посчитать Δ f H?

[Mugen37]

Термохимические данные программа выдает только после проведения анализа частот (frequency). Анализ частот следует проводить только после предварительной оптимизации структуры. Для новых молекул можно сразу задавать такую операцию # opt+freq. После проведения анализа частот в выходном файле появится раздел thermochemistry. В нем приведены необходимые данные для расчета Δ f H. Посмотрите также мануал по термохимии, выложенный на сайте разработчика программы: http://www.gaussian.com/g_whitepap/thermo.htm

[EvgeniX]

ΔfH

И как же в общем случае рассчитать теплоту образования вещества (не полуэмпирическими методами)?

[Mugen37]

ΔfH

И как же в общем случае рассчитать теплоту образования вещества (не полуэмпирическими методами)?

я сначала оптимизирую структуру методом B3LYP (работаю с порфиринами), затем выполняю анализ частот, строка в input файле такая: # freq B3LYP/(базисный набор) в результате расчета в .log файле появится раздел вида (привел не все, а только начало раздела):

-------------------
- Thermochemistry -
-------------------
Temperature 298.150 Kelvin. Pressure 1.00000 Atm.
Atom 1 has atomic number 6 and mass 12.00000
Atom 2 has atomic number 6 and mass 12.00000
Atom 3 has atomic number 1 and mass 1.00783
Atom 4 has atomic number 1 and mass 1.00783
Atom 5 has atomic number 1 and mass 1.00783
Atom 6 has atomic number 1 and mass 1.00783
Atom 7 has atomic number 1 and mass 1.00783
Atom 8 has atomic number 1 and mass 1.00783
Molecular mass: 30.04695 amu.

а дальше используя формулы из мануала (в самом конце пример разобран) http://www.gaussian.com/g_whitepap/thermo.htm, вычисляю теплоту образования вещества.

[EvgeniX]

Ну т.е. по изодесмическим реакциям...

[amge]

Ну т.е. по изодесмическим реакциям...

Можно и так сказать. Но вообще подход мне показался весьма разумным. Подумываю даже, не написать ли программку, в которую забить экспериментальные энтальпии и энтропии атомизации простых веществ (они есть по ссылке), а также расчетные энергии атомов для чаще всего используемых методов/базисов (не очень много счета), чтобы стандартные термодинамические параметры образования вычислялись автоматически.

[EvgeniX]

Хотя нет, не изодесмические реакции. Но думаю, что по изо(гомо)десмическим было бы точнее (ошибки учёта электронной корреляции компенсируются).

[VTur]

Если умеете подбирать гомодесмические реакции (одни из наиболее точных среди десмотических), то начните с них.
Если не получится, вот Вам небольшой текст, может и поможет.

fH.pdf

Но вопрос всегда один, с какой точностью это надо. От этого зависит комбинация метод/базис.

[Driver]

VTur'a с возвращением!

[Nbveh]

Вот статья

XXX.PDF

[VTur]

Спасибо.