Hide keyboard shortcuts

Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

1"""The following contains a database of small molecules 

2 

3Data for the G2/97 database are from 

4Raghavachari, Redfern, and Pople, J. Chem. Phys. Vol. 106, 1063 (1997). 

5See http://www.cse.anl.gov/Catalysis_and_Energy_Conversion/ 

6Computational_Thermochemistry.shtml for the original files. 

7 

8All numbers are experimental values, except for coordinates, which are 

9MP2(full)/6-31G(d) optimized geometries 

10(from http://www.cse.anl.gov/OldCHMwebsiteContent/compmat/G2-97.htm) 

11 

12Atomic species: 

13ref: Curtiss et al. JCP 106, 1063 (1997). 

14'Enthalpy' is the experimental enthalpies of formation at 0K 

15'thermal correction' is the thermal corrections H(298)-H(0) 

16 

17Molecular species: 

18ref: Staroverov et al. JCP 119, 12129 (2003) 

19'Enthalpy' is the experimental enthalpies of formation at 298K 

20'ZPE' is the zero-point energies 

21'thermal correction' is the thermal enthalpy corrections H(298K) - H_exp(0K) 

22ZPE and thermal corrections are estimated from B3LYP geometries and vibrations. 

23 

24Experimental ionization potentials are from http://srdata.nist.gov/cccbdb/. 

25 

26For details about G2-1 and G2-2 sets see doi:10.1063/1.477422. 

27""" 

28 

29from ase.data.g2_1 import data as data_g2_1 

30from ase.data.g2_2 import data as data_g2_2 

31from ase.data.g2_1 import atom_names as atom_names_g2_1 

32from ase.data.g2_1 import molecule_names as molecule_names_g2_1 

33from ase.data.g2_2 import atom_names as atom_names_g2_2 

34from ase.data.g2_2 import molecule_names as molecule_names_g2_2 

35from ase.data.g2_1 import get_ionization_energy 

36from ase.data.g2_1 import get_atomization_energy 

37 

38__all__ = ['data', 'molecule_names', 'atom_names', 'get_ionization_energy', 

39 'get_atomization_energy'] 

40 

41data = data_g2_1.copy() 

42data.update(data_g2_2) 

43 

44atom_names = [] 

45for a in atom_names_g2_1 + atom_names_g2_2: 

46 if a not in atom_names: 

47 atom_names.append(a) 

48molecule_names = molecule_names_g2_1 + molecule_names_g2_2