1# flake8: noqa 

2# atomization energies in kcal / mol (= 43.364 meV) 

3# All values evaluated with PBE xc-orbitals and densities at 

4# experimental geometries. Zero-point vibration has been removed 

5# from experimental energies (from [1]). 

6atomization = { 

7# Molec expt LSD PBE RPBE BLYP 

8'H2' : ( 109.5, 113.2, 104.6, 105.5, 109.4), 

9'LiH' : ( 57.8, 61.0, 53.5, 53.4, 58.1), 

10'CH4' : ( 419.3, 462.3, 419.8, 410.6, 416.6), 

11'NH3' : ( 297.4, 337.3, 301.7, 293.2, 301.4), 

12'OH' : ( 106.4, 124.1, 109.8, 106.3, 109.6), 

13'H2O' : ( 232.2, 266.5, 234.2, 226.6, 232.5), 

14'HF' : ( 140.8, 162.2, 142.0, 137.5, 141.0), 

15'Li2' : ( 24.4, 23.9, 19.9, 20.2, 20.5), 

16'LiF' : ( 138.9, 156.1, 138.6, 132.9, 140.1), 

17'Be2' : ( 3.0, 12.8, 9.8, 7.9, 6.1), 

18'C2H2': ( 405.4, 460.3, 414.9, 400.4, 405.3), 

19'C2H4': ( 562.6, 632.6, 571.5, 554.5, 560.7), 

20'HCN' : ( 311.9, 361.0, 326.1, 313.6, 320.3), 

21'CO' : ( 259.3, 299.1, 268.8, 257.9, 261.8), 

22'N2' : ( 228.5, 267.4, 243.2, 232.7, 239.8), 

23'NO' : ( 152.9, 198.7, 171.9, 161.6, 166.0), 

24'O2' : ( 120.5, 175.0, 143.7, 133.3, 135.3), 

25'F2' : ( 38.5, 78.2, 53.4, 45.6, 49.4), 

26'P2' : ( 117.3, 143.8, 121.1, 114.1, 121.0), 

27'Cl2' : ( 58.0, 83.0, 65.1, 58.9, 57.2) 

28} 

29 

30# exchange-only atomization energies in kcal / mol (= 43.364 meV) 

31# All values evaluated with PBE xc-orbitals and densities at 

32# experimental geometries. (from [1]). 

33ex_atomization = { 

34# Molec exact LSD PBE RPBE BLYP 

35'H2' : ( 84.0, 81.5, 84.8, 85.8, 85.4), 

36'LiH' : ( 33.9, 33.6, 36.9, 36.8, 36.2), 

37'CH4' : ( 327.2, 369.9, 336.0, 326.9, 331.2), 

38'NH3' : ( 199.5, 255.0, 227.4, 218.9, 222.6), 

39'OH' : ( 67.3, 96.2, 84.5, 80.9, 82.7), 

40'H2O' : ( 154.6, 212.9, 183.9, 176.4, 180.5), 

41'HF' : ( 96.1, 136.1, 117.1, 112.6, 115.4), 

42'Li2' : ( 3.5, 6.5, 6.4, 6.7, 3.9), 

43'LiF' : ( 86.8, 129.7, 116.5, 110.8, 113.6), 

44'Be2' : ( -11.0, 9.4, 3.1, 1.2, 1.6), 

45'C2H2': ( 290.6, 382.7, 333.0, 318.5, 325.5), 

46'C2H4': ( 423.9, 517.7, 456.5, 439.6, 447.4), 

47'HCN' : ( 194.5, 294.0, 256.1, 243.5, 249.1), 

48'CO' : ( 169.2, 261.9, 224.0, 213.1, 218.7), 

49'N2' : ( 110.2, 211.4, 184.1, 173.6, 177.6), 

50'NO' : ( 45.6, 156.9, 122.8, 112.5, 117.0), 

51'O2' : ( 24.9, 147.5, 104.4, 94.1, 99.3), 

52'F2' : ( -43.3, 64.0, 32.5, 24.7, 28.8), 

53'P2' : ( 31.8, 98.4, 73.1, 66.1, 70.1), 

54'Cl2' : ( 15.5, 68.2, 39.8, 33.7, 37.0) 

55} 

56 

57# Exchange energy of some spherical atoms in Hartrees (= 27.211 eV) 

58# All functionals were evaluated with self-consistent exchange-only 

59# OEP orbitals and densities (from [1]). 

60ex_energy = { 

61# Atom exact LSD PBE RPBE BLYP 

62'H' : ( 0.3125, 0.2680, 0.3059, 0.3112, 0.3098), 

63'He' : ( 1.0258, 0.8840, 1.0136, 1.0313, 1.0255), 

64'Li' : ( 1.7807, 1.5379, 1.7572, 1.7876, 1.7753), 

65'Be' : ( 2.6658, 2.3124, 2.6358, 2.6801, 2.6578), 

66'N' : ( 6.6044, 5.9008, 6.5521, 6.6252, 6.5961), 

67'Ne' : ( 12.1050, 11.0335, 12.0667, 12.1593, 12.1378), 

68'Na' : ( 14.0131, 12.7859, 13.9506, 14.0528, 14.0304), 

69'Mg' : ( 15.9884, 14.6117, 15.9147, 16.0260, 16.0005), 

70'P' : ( 22.6341, 20.7931, 22.5028, 22.6369, 22.6221), 

71'Ar' : ( 30.1747, 27.8632, 29.9961, 30.1494, 30.1535), 

72'Kr' : ( 93.8330, 88.6245, 93.4257, 93.6645, 93.8721), 

73'Xe' : ( 179.0635, 170.5660, 178.2450, 178.5649, 179.0427) 

74} 

75 

76# Correlation energy of some spherical atoms in Hartrees (= 27.211 eV) 

77# All functionals were evaluated with self-consistent exchange-only 

78# OEP orbitals and densities (from [1]). 

79# 'Exact' values are from reference [2] 

80ec_energy = { 

81# Atom exact LSD PBE BLYP 

82'H' : ( 0.0000, 0.0222, 0.0060, 0.0000), 

83'He' : ( 0.0420, 0.1125, 0.0420, 0.0438), 

84'Li' : ( 0.0455, 0.1508, 0.0514, 0.0534), 

85'Be' : ( 0.0950, 0.2240, 0.0856, 0.0945), 

86'N' : ( 0.1858, 0.4268, 0.1799, 0.1919), 

87'Ne' : ( 0.3929, 0.7428, 0.3513, 0.3835), 

88'Na' : ( 0.3988, 0.8010, 0.3715, 0.4083), 

89'Mg' : ( 0.4424, 0.8874, 0.4110, 0.4594), 

90'P' : ( 0.5446, 1.1127, 0.5265, 0.5664), 

91'Ar' : ( 0.7314, 1.4242, 0.7067, 0.7508), 

92'Kr' : ( 3.2693, 1.7672, 1.7486, 2.0788), 

93'Xe' : ( 5.1773, 2.9184, 2.7440, 3.1789) 

94} 

95 

96# atomization energies in kcal / mol (= 43.364 meV). 

97# All values evaluated with self-consistent orbitals and densities. 

98# Geometry optimization to within 10meV/Ang. 

99# Data from reference [3]. 

100atomization_vasp = { 

101# Molecule Expt. PBE_VASP PBE_G03 PBE0_VASP PBE0_G03 

102'LiH' : ( 58, 53.5, 53.5, 52.6, 52.9,), 

103'BeH' : ( 48, 55.5, 55.6, 55.4, 56.0,), 

104'CH' : ( 84, 84.7, 84.8, 83.3, 83.0,), 

105'CH2_s3B1d' : ( 189, 194.4, 194.6, 193.9, 193.8,), 

106'CH2_s1A1d' : ( 182, 178.8, 179.1, 176.3, 176.5,), 

107'CH3' : ( 306, 309.7, 310.1, 308.3, 308.6,), 

108'CH4' : ( 420, 419.6, 420.2, 417.2, 417.9,), 

109'NH' : ( 82, 88.6, 88.6, 85.2, 85.3,), 

110'NH2' : ( 182, 188.7, 188.9, 183.4, 183.3,), 

111'NH3' : ( 297, 301.7, 302.3, 294.7, 295.3,), 

112'OH' : ( 107, 109.7, 110.1, 105.4, 105.8,), 

113'H2O' : ( 233, 233.7, 234.5, 226.4, 227.3,), 

114'HF' : ( 142, 141.5, 142.2, 136.2, 137.0,), 

115'SiH2_s1A1d': ( 154, 147.9, 148.0, 147.2, 147.4,), 

116'SiH2_s3B1d': ( 131, 131.3, 131.8, 132.2, 132.5,), 

117'SiH3' : ( 226, 222.2, 222.6, 223.5, 223.8,), 

118'SiH4' : ( 324, 313.3, 313.7, 315.0, 315.7,), 

119'PH2' : ( 153, 154.5, 154.6, 153.0, 153.2,), 

120'PH3' : ( 241, 239.0, 239.3, 237.3, 237.5,), 

121'SH2' : ( 182, 182.0, 182.2, 179.8, 180.0,), 

122'HCl' : ( 107, 106.3, 106.5, 104.4, 105.0,), 

123'Li2' : ( 26, 19.9, 20.1, 19.3, 19.3,), 

124'LiF' : ( 139, 138.4, 139.0, 131.0, 131.9,), 

125'C2H2' : ( 404, 414.5, 415.1, 404.5, 404.7,), 

126'C2H4' : ( 562, 571.0, 571.9, 563.8, 564.2,), 

127'C2H6' : ( 711, 716.0, 717.1, 711.6, 712.4,), 

128'CN' : ( 179, 197.5, 197.7, 179.1, 179.1,), 

129'HCN' : ( 313, 326.3, 326.5, 311.1, 311.5,), 

130'CO' : ( 261, 268.6, 269.1, 255.3, 255.8,), 

131'HCO' : ( 279, 294.9, 295.5, 280.5, 280.9,), 

132'H2CO' : ( 376, 385.5, 386.3, 371.9, 372.8,), 

133'CH3OH' : ( 513, 519.3, 520.4, 509.0, 510.3,), 

134'N2' : ( 227, 243.7, 243.9, 225.3, 225.9,), 

135'N2H4' : ( 437, 452.7, 453.7, 437.9, 438.8,), 

136'NO' : ( 153, 172.0, 172.5, 153.3, 153.8,), 

137'O2' : ( 118, 143.3, 144.0, 124.1, 124.9,), 

138'H2O2' : ( 268, 281.6, 282.6, 262.7, 263.8,), 

139'F2' : ( 38, 52.6, 53.0, 35.2, 35.3,), 

140'CO2' : ( 392, 415.4, 416.5, 390.8, 392.0,), 

141'Na2' : ( 19, 17.7, 18.1, 15.6, 15.9,), 

142'Si2' : ( 74, 81.3, 81.4, 76.5, 77.3,), 

143'P2' : ( 116, 121.5, 121.7, 111.8, 111.7,), 

144'S2' : ( 98, 115.4, 115.2, 107.3, 107.0,), 

145'Cl2' : ( 57, 65.8, 65.8, 60.1, 59.9,), 

146'NaCl' : ( 99, 93.6, 94.5, 92.1, 93.6,), 

147'SiO' : ( 191, 195.6, 196.6, 182.2, 183.3,), 

148'CS' : ( 172, 179.5, 179.6, 168.0, 168.2,), 

149'SO' : ( 122, 141.5, 141.3, 127.9, 127.3,), 

150'ClO' : ( 62, 81.6, 81.5, 67.4, 67.6,), 

151'ClF' : ( 62, 72.3, 72.5, 61.3, 61.3,), 

152'Si2H6' : ( 533, 519.5, 520.4, 522.2, 523.3,), 

153'CH3Cl' : ( 395, 399.4, 400.2, 395.0, 395.7,), 

154'CH3SH' : ( 473, 477.8, 478.6, 472.7, 473.5,), 

155'HOCl' : ( 165, 175.2, 175.7, 162.9, 163.3,), 

156'SO2' : ( 253, 281.1, 280.7, 254.1, 253.5,), 

157} 

158 

159# Experimental, and calculated bindinglengths of 16 diatomic molecules 

160# of the G2-1 test set. In Angstroms. 

161# Data from reference [3]. 

162diatomic = { 

163#System Expt. PBEVASP PBEG03 PBE0VASP PBE0G03 

164'BeH': (1.343, 1.354, 1.353, 1.350, 1.348), 

165'CH' : (1.120, 1.136, 1.136, 1.124, 1.124), 

166'Cl2': (1.988, 1.999, 2.004, 1.973, 1.978), 

167'ClF': (1.628, 1.648, 1.650, 1.614, 1.617), 

168'ClO': (1.570, 1.576, 1.577, 1.554, 1.555), 

169'CN' : (1.172, 1.173, 1.174, 1.159, 1.159), 

170'CO' : (1.128, 1.136, 1.135, 1.122, 1.122), 

171'F2' : (1.412, 1.414, 1.413, 1.377, 1.376), 

172'FH' : (0.917, 0.932, 0.930, 0.919, 0.918), 

173'HCl': (1.275, 1.287, 1.288, 1.276, 1.278), 

174'Li2': (2.673, 2.728, 2.728, 2.727, 2.727), 

175'LiF': (1.564, 1.583, 1.575, 1.571, 1.561), 

176'LiH': (1.595, 1.604, 1.604, 1.602, 1.597), 

177'N2' : (1.098, 1.103, 1.102, 1.089, 1.089), 

178'O2' : (1.208, 1.218, 1.218, 1.193, 1.192), 

179'Na2': (3.079, 3.087, 3.076, 3.086, 3.086), 

180} 

181 

182## Note the difference between the experimental data from Blaha [1] and 

183## from Kresse [3]: 

184## Mol BLAHA KRESSE 

185## LiH : 57.8 58.0 

186## CH4 : 419.3 420.0 

187## NH3 : 297.4 297.0 

188## OH : 106.4 107.0 

189## H2O : 232.2 233.0 

190## HF : 140.8 142.0 

191## Li2 : 24.4 26.0 

192## LiF : 138.9 139.0 

193## C2H2: 405.4 404.0 

194## C2H4: 562.6 562.0 

195## HCN : 311.9 313.0 

196## CO : 259.3 261.0 

197## N2 : 228.5 227.0 

198## NO : 152.9 153.0 

199## O2 : 120.5 118.0 

200## F2 : 38.5 38.0 

201## P2 : 117.3 116.0 

202## Cl2 : 58.0 57.0 

203## ----------------- 

204## MAE : 0.0 1.0 

205 

206## and the difference between the PBE results of Blaha [1] and those from 

207## VASP and GAUSSIAN-03 [3] 

208## Mol BLAHA VASP G03 

209## LiH : 53.5 53.5 53.5 

210## CH4 : 419.8 419.6 420.2 

211## NH3 : 301.7 301.7 302.3 

212## OH : 109.8 109.7 110.1 

213## H2O : 234.2 233.7 234.5 

214## HF : 142.0 141.5 142.2 

215## Li2 : 19.9 19.9 20.1 

216## LiF : 138.6 138.4 139.0 

217## C2H2: 414.9 414.5 415.1 

218## C2H4: 571.5 571.0 571.9 

219## HCN : 326.1 326.3 326.5 

220## CO : 268.8 268.6 269.1 

221## N2 : 243.2 243.7 243.9 

222## NO : 171.9 172.0 172.5 

223## O2 : 143.7 143.3 144.0 

224## F2 : 53.4 52.6 53.0 

225## P2 : 121.1 121.5 121.7 

226## Cl2 : 65.1 65.8 65.8 

227## ----------------------- 

228## MAE : 0.0 0.3 0.4 

229## 

230## Where in the last two, geometry optimization has been performed, but not in 

231## the first. 

232 

233# References: 

234# [1]: 

235# Kurth, Perdew, and Blaha 

236# Molecular and Solid-State Tests of Density Functional Approximations 

237# International Journal of Quantum Chemistry, Vol. 85, 889-909 (1999) 

238# [2]: 

239# Krieger, Chen, Iafrate, and Savin 

240# In Electron Correlations and Materials Properties 

241# Gonis and Kioussis; eds. 

242# Plenum: New York, 1999. 

243# [3]: 

244# Paier, Hirschl, Marsman, and Kresse 

245# The Perdew-Burke-Ernzerhof exchange-correlation functional applied to the 

246# G2-1 test set using a plane wave basis set 

247# The Journal of Chemical Physics, Vol 122, 234102 (2005) 

248 

249def convert(input, column): 

250 keys = sorted(input.keys()) 

251 data = {} 

252 for k in keys: 

253 data[k] = input[k][column] 

254 return data