Stefan Grimme, Jens Antony, Stephan Ehrlich, and Helge Krieg, J. Chem. Phys. 132, 154104 (2010)
Jan 18, 2016
Stefan Grimme, Jens Antony, Stephan Ehrlich, and Helge Krieg, J. Chem. Phys. 132, 154104 (2010)
Forces de dispersió: forces intermoleculars febles entre sistemes no polars (multipols temporals).
Problemes amb l’ús de la teoria DFT: errors importants.
Existència de diferents mètodes per corregir els errors: entre ells DFT-D3.
Presenta avantatges respecte altres mètodes: menor grau d’empiricisme, major acuracitat, etc.
Es basa en:
Dades necessàries: Geometria òptima (coord. xyz) Funcional DFT
Hpiqc, script: “GrimmeD3”
S’obté la correcció de l’Energia de Dispersió per al cas concret (molècula+funcional)
Opt Freq Total E & optim geometry
Hpiqc/users/josepm/FEIV/DFTD3 estructura1_bp
Angstroms
GrimmeD3 estructura1_bp b3-lypD3 Grime E:
Edisp /kcal,au: -50.2160 -0.08002435Total E + D3 Grime E = DFT-D3 Energy
*Exemples: Hpcc/users/mgarcia/DFT-D3_exemples
Hpiqc/users/josepm/FEIV/DFTD3 estructura1_bp.gd3
Structure Total E (au)Total E
Kcal/mol)D3 Grime (Kcal/mol)
DFTD3 (Kcal/mol)
DFTD3 (Kcal/mol)
Structure1 -2258.727585 -1416626.906 -50.216 -1416677.122 2.99TS12 -2258.709860 -1416615.789 -50.9954 -1416666.785 13.33
Structure2 -2258.726607 -1416626.293 -49.4517 -1416675.745 4.37TS23 -2258.727180 -1416626.652 -49.7535 -1416676.405 3.71
Structure3 -2258.733515 -1416630.625 -49.4872 -1416680.113 0.00TS34 -2258.712781 -1416617.621 -51.3037 -1416668.925 11.19
Structure4 -2258.730673 -1416628.843 -50.7911 -1416679.634 0.48
Structure Total E (au) E (au)E
Kcal/mol)D3 Grime (Kcal/mol)
E+D3 Grime (Kcal/mol)
DFTD3 (Kcal/mol)
Structure1 -2258.727585 0.0059305400 3.72 -50.216 -46.49 2.99TS12 -2258.709860 0.0236555000 14.84 -50.9954 -36.15 13.34
Structure2 -2258.726607 0.0069078700 4.33 -49.4517 -45.12 4.37TS23 -2258.727180 0.0063354800 3.98 -49.7535 -45.78 3.71
Structure3 -2258.733515 0.0000000000 0.00 -49.4872 -49.49 0.00TS34 -2258.712781 0.0207342100 13.01 -51.3037 -38.29 11.19
Structure4 -2258.730673 0.0028422200 1.78 -50.7911 -49.01 0.48