A toolkit for ab initio thermochemistry
1. G4(MP2) enthalpy of formation for the PPE1694 dataset
2. G4(MP2)-XP enthalpy of formation for the PPE1694 dataset
3. G4(MP2) enthalpy of formation for 270 entries in the G3/05 dataset
4. G4(MP2)-XP enthalpy of formation for 270 entries in the G3/05 dataset
5. G4(MP2)-XP enthalpy of formation for C60
The enclosed file PPE1694_G4MP2.csv
contain geometries, frequencies and thermochemical energies of all 1694 molecules in the PPE1694 dataset calculated at the G4MP2 level. The file ATOMS_G4MP2.csv
contains energies for the atoms using which atomization energy and standard formation enthalpy are calculated. Here is a screenshot of the output printed when running the python code enclosed in the zip.
$ python3.9 calc_hof.py | tee outfile
#Mol. Atm. E Form. H Exp. form. H [all in kcal/mol]
1 22.1762 53.3112 51.6000
2 136.5204 -80.4132 -80.7800
3 257.4797 -27.7258 -26.4100
4 225.0425 0.0114 0.0000
5 35.6105 1.3222 0.0000
...
...
...
1690 2449.2282 12.4709 12.3000
1691 2449.0343 12.8907 14.8000
1692 2111.7479 11.3013 14.2000
1693 1519.7924 195.8740 198.6000
1694 1208.4291 -54.2314 -59.6000
To print geometry and frequencies, please edit the following lines in calc_hof.py
print_geom='false' # Set to 'true' for printing geometries in g4mp2_geom.xyz
print_freq='false' # Set to 'true' for printing frequencies in g4mp2_freq.dat
The enclosed file PPE1694_G4MP2-XP.csv
contain geometries, frequencies and thermochemical energies of all 1694 molecules in the PPE1694 dataset calculated at the G4MP2-XP level. The file ATOMS_G4MP2-XP.csv
contains energies for the atoms using which atomization energy and standard formation enthalpy are calculated. Here is a screenshot of the output printed when running the python code enclosed in the zip
$ python3.9 calc_hof.py | tee outfile
#Mol. Atm. E Form. H Exp. form. H [all in kcal/mol]
1 22.2452 53.2421 51.6000
2 136.7751 -80.6680 -80.7800
3 258.0407 -28.2868 -26.4100
4 225.3437 -0.2898 0.0000
5 36.0511 0.8817 0.0000
...
...
...
1690 2450.5240 11.1751 12.3000
1691 2450.3301 11.5949 14.8000
1692 2112.8919 10.1574 14.2000
1693 1520.6502 195.0162 198.6000
1694 1209.9533 -55.7556 -59.6000
To print geometry and frequencies, please edit the following lines in calc_hof.py
print_geom='false' # Set to 'true' for printing geometries in g4mp2_geom.xyz
print_freq='false' # Set to 'true' for printing frequencies in g4mp2_freq.dat
The python code contains G4(MP2)-XP level geometry, harmonic frequencies, energies from CCSD(T), HF, and MP2 methods, and demonstrates how the atomization energy and standard formation enthalpy of buckminsterfullerene, C60 can be calculated.
After installing the pople code, running the above C60_enthalpy.py
with python returns the following output
Atomization energy is 15.2943 hartree
Atomization energy is 9597.3248 kcal/mol
ZPVE is 230.0659 kcal/mol
HOF_atoms(C) HOF molecule [all in kcal/mol]
169.9800 600.4915
170.0270 603.3115
170.1100 608.2915