Destio небольшой блог

Термохимия в Gaussian

Для того чтобы Gaussian рассчитал всю термохимию достаточно указать параметр freq во входном файле.

При этом, дополнительно к электронной энергии, будут рассчитаны такие параметры как энергия нулевых колебаний ZPE, энтропия и термические поправки к энергии. Из этих чисел несложно получить энергию Гиббса и энтальпию, которые тоже появляются в результатах.

По умолчанию расчет идет для температуры 298.15 Кельвинов, со стандартными массами атомов как смесей каноничных изотопов.

В результате, среди прочих параметров мы получим вывод нужных нам параметров:

 Zero-point correction=                           0.139798 (Hartree/Particle)
 Thermal correction to Energy=                    0.149467
 Thermal correction to Enthalpy=                  0.150412
 Thermal correction to Gibbs Free Energy=         0.105081
 Sum of electronic and zero-point Energies=           -493.527050
 Sum of electronic and thermal Energies=              -493.517380
 Sum of electronic and thermal Enthalpies=            -493.516436
 Sum of electronic and thermal Free Energies=         -493.561767

Важно не забывать что ВСЕ 8 значений энергии включают в себя ZPE (даже сам ZPE!). Бывало ээ что порой я учитывал ZPE дважды. Термическая поправка (Et) к энергии складывается из составных частей, указанных в следующей таблице. Энтальпия из энергии получается нехитро: H = E + kT. Энергия Гиббса складывается из энтальпии и энтропии: G = H - TS.

Следующие четыре строки, отвечающие за суммы, получаются из предыдущих путем прибавления электронной энергии ко всем найденным значениям и поправкам. В данном примере электронная энергия равна -493.666847632 хартри.

Чуть ниже можно увидеть подробные данные по теплоемкости и энтропии. Почему-то все в дурацких калориях:

                    E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
 Total                   93.792             35.443             95.406
 Electronic               0.000              0.000              0.000
 Translational            0.889              2.981             40.722
 Rotational               0.889              2.981             29.441
 Vibrational             92.015             29.482             25.243

Численное значение энтропии можно взять только отсюда.

Чтобы не ворошить каждый выходной файл в поисках хартри и калорий, можно использовать Quantum Wizard

Ярлыки: ,


Ближайшие сообщения:

Похожие сообщения:

Комментарии к этой записи отключены.