search

patonlab / GoodVibes

Calculate quasi-harmonic free energies from Gaussian output files with temperature and other corrections
http://www.patonlab.colostate.edu
MIT License
137 stars 51 forks source link

SPC energy with double hybrid functional #51

Closed npe1011 closed 2 years ago

npe1011 commented 2 years ago

Dear developers,

I am using a double hybrid functional method (DSDPBEP86/def2TZVPP) in Gaussian 16 rev.C for single point calculation after opt and freq calculation using a standard hybrid functional (M062X/def2SVP). When I run goodvibes 3.0.1 to combine the results (using --spc option) to calculate the final free energies, the goodvibes program seemed to read the SCF energy of the DSDPBEP86 results as printed: SCF Done: E(RDSDPBEP86) = @@@@@@@@ A.U. after 8 cycles But the correct final EE of double hybrid methods would be the sum of SCF energy and E2 correction after MP2 calculation, as printed: E2(DSDPBEP86) = @@@@@@@@@@@@@@@@@@ E(DSDPBEP86) = @@@@@@@@@@@@@@@

I would really appreciate if you could check this point.

Accutally, the official GaussView ver 6 seems to read the SCF energy without E2 corrections.

jvalegre commented 2 years ago

Hi npe1011 - I just added some lines to account for calculations with double hybrid functionals, could you make sure that GoodVibes is working correctly for you too? The changes are in the latest commit (download + installation through setup.py).

Thanks for reporting this issue!

npe1011 commented 2 years ago

Thank you very much for a quick reply and change. The new version works fine with my outputs.