Hello everyone,
I am calculating reaction energies for various small molecules and thermodynamic corrections for high temperatures (650 - 950 K). On the MP2 level I am using numforce to calculate the vibr. force constants with the option -central. I used the default displacement of 0.02. All geometries have only real modes and the SCF and density convergence criteria are tight, respectively.
Here comes my problem: For 950 K the thermodynamic corrections (freeh) are pretty off. For some reactions, the Gibbs free energy rises at high temperature, as opposed to expectation (and GGA-level aoforce calculations). In most cases, however, the Gibbs free energies are way too low, so there is no similar behaviour or trend through all the reactions. This puzzles me. Some of the dG values are completely unrealistic (although corresponding MP2 energies and dG for lower temperatures from the same calculation are fine).
Is this odd behaviour of numforce thermodynamic properties at high temperatures known? Is there a solution or a way around? Or do I have to stick to lower temperatures when applying numforce?
Best regards and many thanks for suggestions!
Stegi