Gibbs free Energy Calculations

[Habib]

Hello,

I'm trying to figure out the Gibbs free energy for a reaction on the LaFeO3(100) surface. I read some forum posts and found out that the "IdealGasThermochemistry" object can help with this, but my QuantumATK version (Q-2019.12) doesn't have it. We're planning to upgrade, but it'll take some time.
To calculate Δ𝐺, we need to find Δ𝛦, Δ𝐸𝑍𝑃𝐸, and Δ𝑆. I already have Δ𝛦, which is -0.06 eV, and I can get 𝐸𝑍𝑃𝐸 and S for the isolated molecules from a database. But I'm not sure how to calculate the Zero-Point Energy (ZPE) corrections and the entropy change for the LaFeO3(100)@H system in QuantumATK Q-2019.12.
Is there an easy way to do this or a script that can help me extract ZPE and S from the optimized file of LaFeO3@H? Your guidance would be much appreciated!

Kind Regards
Habib

[Habib]

If someone can help?

[Anders Blom]

The ZPE can be obtained as a query method on a PhononDensityOfStates (https://docs.quantumatk.com/manual/Types/PhononDensityOfStates/PhononDensityOfStates.html) , but computing the entropy is too complex to provide a simple answer to here, as it includes translational, rotational, vibrational, pressure and electronic contributions.