Coverage for /builds/kinetik161/ase/ase/data/g2_1.py: 29.17%

24 statements  

« prev     ^ index     » next       coverage.py v7.2.7, created at 2023-12-10 11:04 +0000

1# flake8: noqa 

2""" 

3The following contains a database of small molecules 

4 

5Data for the G2/97 database are from 

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

7See http://www.cse.anl.gov/Catalysis_and_Energy_Conversion/Computational_Thermochemistry.shtml for the original files. 

8 

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

10MP2(full)/6-31G(d) optimized geometries (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 

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

25 

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

27Information presented on these pages is considered public information 

28and may be distributed or copied http://www.nist.gov/public_affairs/disclaimer.cfm 

29 

30""" 

31 

32from ase.symbols import string2symbols 

33 

34atom_names = ['H', 'Li', 'Be', 'C', 'N', 'O', 'F', 'Na', 'Si', 'P', 'S', 'Cl'] 

35 

36molecule_names = ['LiH', 'BeH', 'CH', 'CH2_s3B1d', 'CH2_s1A1d', 'CH3', 'CH4', 'NH', 'NH2', 'NH3', 'OH', 'H2O', 'HF', 'SiH2_s1A1d', 'SiH2_s3B1d', 'SiH3', 'SiH4', 'PH2', 'PH3', 'SH2', 'HCl', 'Li2', 'LiF', 'C2H2', 'C2H4', 

37 'C2H6', 'CN', 'HCN', 'CO', 'HCO', 'H2CO', 'CH3OH', 'N2', 'N2H4', 'NO', 'O2', 'H2O2', 'F2', 'CO2', 'Na2', 'Si2', 'P2', 'S2', 'Cl2', 'NaCl', 'SiO', 'CS', 'SO', 'ClO', 'ClF', 'Si2H6', 'CH3Cl', 'CH3SH', 'HOCl', 'SO2'] 

38 

39data = { 

40 'Be': {'CAS No.': 7440417, 

41 'charges': None, 

42 'database': 'G2-1', 

43 'description': 'Be atom', 

44 'enthalpy': 76.48, 

45 'ionization energy': 9.32, 

46 'magmoms': None, 

47 'name': 'Beryllium', 

48 'positions': [[0.0, 0.0, 0.0]], 

49 'symbols': 'Be', 

50 'thermal correction': 0.46}, 

51 'BeH': {'CAS No.': 13597972, 

52 'ZPE': 2.9073, 

53 'charges': None, 

54 'database': 'G2-1', 

55 'description': 'Beryllium hydride (BeH), D*h symm.', 

56 'enthalpy': 81.7, 

57 'ionization energy': 8.21, 

58 'magmoms': [0.8, 0.2], 

59 'name': 'BeH (beryllium monohydride)', 

60 'positions': [[0.0, 0.0, 0.269654], [0.0, 0.0, -1.078616]], 

61 'symbols': 'BeH', 

62 'thermal correction': 2.0739}, 

63 'C': {'CAS No.': 7440440, 

64 'charges': None, 

65 'database': 'G2-1', 

66 'description': 'C atom', 

67 'enthalpy': 169.98, 

68 'ionization energy': 11.26, 

69 'magmoms': [2.0], 

70 'name': 'Carbon', 

71 'positions': [[0.0, 0.0, 0.0]], 

72 'symbols': 'C', 

73 'thermal correction': 0.25}, 

74 'C2H2': {'CAS No.': 74862, 

75 'ZPE': 16.6001, 

76 'charges': None, 

77 'database': 'G2-1', 

78 'description': 'Acetylene (C2H2), D*h symm.', 

79 'enthalpy': 54.2, 

80 'ionization energy': 11.4, 

81 'magmoms': None, 

82 'name': 'C_2H_2', 

83 'positions': [[0.0, 0.0, 0.60808], 

84 [0.0, 0.0, -0.60808], 

85 [0.0, 0.0, -1.67399], 

86 [0.0, 0.0, 1.67399]], 

87 'symbols': 'CCHH', 

88 'thermal correction': 2.4228, 

89 'vertical ionization energy': 11.49}, 

90 'C2H4': {'CAS No.': 74851, 

91 'ZPE': 31.5267, 

92 'charges': None, 

93 'database': 'G2-1', 

94 'description': 'Ethylene (H2C=CH2), D2h symm.', 

95 'enthalpy': 12.5, 

96 'ionization energy': 11.4, 

97 'magmoms': None, 

98 'name': 'C_2H_4', 

99 'positions': [[0.0, 0.0, 0.66748], 

100 [0.0, 0.0, -0.66748], 

101 [0.0, 0.922832, 1.237695], 

102 [0.0, -0.922832, 1.237695], 

103 [0.0, 0.922832, -1.237695], 

104 [0.0, -0.922832, -1.237695]], 

105 'symbols': 'CCHHHH', 

106 'thermal correction': 2.51, 

107 'vertical ionization energy': 11.49}, 

108 'C2H6': {'CAS No.': 74840, 

109 'ZPE': 46.095, 

110 'charges': None, 

111 'database': 'G2-1', 

112 'description': 'Ethane (H3C-CH3), D3d symm.', 

113 'enthalpy': -20.1, 

114 'magmoms': None, 

115 'name': 'C_2H_6', 

116 'positions': [[0.0, 0.0, 0.762209], 

117 [0.0, 0.0, -0.762209], 

118 [0.0, 1.018957, 1.157229], 

119 [-0.882443, -0.509479, 1.157229], 

120 [0.882443, -0.509479, 1.157229], 

121 [0.0, -1.018957, -1.157229], 

122 [-0.882443, 0.509479, -1.157229], 

123 [0.882443, 0.509479, -1.157229]], 

124 'symbols': 'CCHHHHHH', 

125 'thermal correction': 2.7912}, 

126 'CH': {'CAS No.': 3315375, 

127 'ZPE': 3.9659, 

128 'charges': None, 

129 'database': 'G2-1', 

130 'description': 'CH radical. Doublet, C*v symm.', 

131 'enthalpy': 142.5, 

132 'ionization energy': 10.64, 

133 'magmoms': [1.0, 0.0], 

134 'name': 'CH(Methylidyne)', 

135 'positions': [[0.0, 0.0, 0.160074], [0.0, 0.0, -0.960446]], 

136 'symbols': 'CH', 

137 'thermal correction': 2.0739}, 

138 'CH2_s1A1d': {'CAS No.': 2465567, 

139 'ZPE': 10.2422, 

140 'charges': None, 

141 'database': 'G2-1', 

142 'description': 'Singlet methylene (CH2), C2v symm, 1-A1.', 

143 'enthalpy': 102.8, 

144 'magmoms': None, 

145 'name': 'CH_2 (^1A_1)', 

146 'positions': [[0.0, 0.0, 0.174343], 

147 [0.0, 0.862232, -0.523029], 

148 [0.0, -0.862232, -0.523029]], 

149 'symbols': 'CHH', 

150 'thermal correction': 2.3745}, 

151 'CH2_s3B1d': {'CAS No.': 2465567, 

152 'ZPE': 10.6953, 

153 'charges': None, 

154 'database': 'G2-1', 

155 'description': 'Triplet methylene (CH2), C2v symm, 3-B1.', 

156 'enthalpy': 93.7, 

157 'ionization energy': 10.4, 

158 'magmoms': [2.0, 0.0, 0.0], 

159 'name': 'CH_2 (^3B_1)', 

160 'positions': [[0.0, 0.0, 0.110381], 

161 [0.0, 0.982622, -0.331142], 

162 [0.0, -0.982622, -0.331142]], 

163 'symbols': 'CHH', 

164 'thermal correction': 2.3877}, 

165 'CH3': {'CAS No.': 2229074, 

166 'ZPE': 18.3383, 

167 'charges': None, 

168 'database': 'G2-1', 

169 'description': 'Methyl radical (CH3), D3h symm.', 

170 'enthalpy': 35.0, 

171 'ionization energy': 9.84, 

172 'magmoms': [1.0, 0.0, 0.0, 0.0], 

173 'name': 'CH_3', 

174 'positions': [[0.0, 0.0, 0.0], 

175 [0.0, 1.07841, 0.0], 

176 [0.93393, -0.539205, 0.0], 

177 [-0.93393, -0.539205, 0.0]], 

178 'symbols': 'CHHH', 

179 'thermal correction': 2.5383}, 

180 'CH3Cl': {'CAS No.': 74873, 

181 'ZPE': 23.3013, 

182 'charges': None, 

183 'database': 'G2-1', 

184 'description': 'Methyl chloride (CH3Cl), C3v symm.', 

185 'enthalpy': -19.6, 

186 'ionization energy': 11.26, 

187 'magmoms': None, 

188 'name': 'CH_3Cl', 

189 'positions': [[0.0, 0.0, -1.121389], 

190 [0.0, 0.0, 0.655951], 

191 [0.0, 1.029318, -1.47428], 

192 [0.891415, -0.514659, -1.47428], 

193 [-0.891415, -0.514659, -1.47428]], 

194 'symbols': 'CClHHH', 

195 'thermal correction': 2.4956, 

196 'vertical ionization energy': 11.29}, 

197 'CH3OH': {'CAS No.': 67561, 

198 'ZPE': 31.6635, 

199 'charges': None, 

200 'database': 'G2-1', 

201 'description': 'Methanol (CH3-OH), Cs symm.', 

202 'enthalpy': -48.0, 

203 'ionization energy': 10.84, 

204 'magmoms': None, 

205 'name': 'H_3COH', 

206 'positions': [[-0.047131, 0.664389, 0.0], 

207 [-0.047131, -0.758551, 0.0], 

208 [-1.092995, 0.969785, 0.0], 

209 [0.878534, -1.048458, 0.0], 

210 [0.437145, 1.080376, 0.891772], 

211 [0.437145, 1.080376, -0.891772]], 

212 'symbols': 'COHHHH', 

213 'thermal correction': 2.6832, 

214 'vertical ionization energy': 10.96}, 

215 'CH3SH': {'CAS No.': 74931, 

216 'ZPE': 28.3973, 

217 'charges': None, 

218 'database': 'G2-1', 

219 'description': 'Methanethiol (H3C-SH), Staggered, Cs symm.', 

220 'enthalpy': -5.5, 

221 'ionization energy': 9.44, 

222 'magmoms': None, 

223 'name': 'H_3CSH', 

224 'positions': [[-0.047953, 1.149519, 0.0], 

225 [-0.047953, -0.664856, 0.0], 

226 [1.283076, -0.823249, 0.0], 

227 [-1.092601, 1.461428, 0.0], 

228 [0.432249, 1.551207, 0.892259], 

229 [0.432249, 1.551207, -0.892259]], 

230 'symbols': 'CSHHHH', 

231 'thermal correction': 2.869, 

232 'vertical ionization energy': 9.44}, 

233 'CH4': {'CAS No.': 74828, 

234 'ZPE': 27.6744, 

235 'charges': None, 

236 'database': 'G2-1', 

237 'description': 'Methane (CH4), Td symm.', 

238 'enthalpy': -17.9, 

239 'ionization energy': 12.64, 

240 'magmoms': None, 

241 'name': 'CH_4', 

242 'positions': [[0.0, 0.0, 0.0], 

243 [0.629118, 0.629118, 0.629118], 

244 [-0.629118, -0.629118, 0.629118], 

245 [0.629118, -0.629118, -0.629118], 

246 [-0.629118, 0.629118, -0.629118]], 

247 'symbols': 'CHHHH', 

248 'thermal correction': 2.3939, 

249 'vertical ionization energy': 13.6}, 

250 'CN': {'CAS No.': 2074875, 

251 'ZPE': 3.0183, 

252 'charges': None, 

253 'database': 'G2-1', 

254 'description': 'Cyano radical (CN), C*v symm, 2-Sigma+.', 

255 'enthalpy': 104.9, 

256 'ionization energy': 13.6, 

257 'magmoms': [1.0, 0.0], 

258 'name': 'CN (Cyano radical)', 

259 'positions': [[0.0, 0.0, -0.611046], [0.0, 0.0, 0.523753]], 

260 'symbols': 'CN', 

261 'thermal correction': 2.0739}, 

262 'CO': {'CAS No.': 630080, 

263 'ZPE': 3.1062, 

264 'charges': None, 

265 'database': 'G2-1', 

266 'description': 'Carbon monoxide (CO), C*v symm.', 

267 'enthalpy': -26.4, 

268 'ionization energy': 14.01, 

269 'magmoms': None, 

270 'name': 'CO', 

271 'positions': [[0.0, 0.0, 0.493003], [0.0, 0.0, -0.657337]], 

272 'symbols': 'OC', 

273 'thermal correction': 2.0739, 

274 'vertical ionization energy': 14.01}, 

275 'CO2': {'CAS No.': 124389, 

276 'ZPE': 7.313, 

277 'charges': None, 

278 'database': 'G2-1', 

279 'description': 'Carbon dioxide (CO2), D*h symm.', 

280 'enthalpy': -94.1, 

281 'ionization energy': 13.78, 

282 'magmoms': None, 

283 'name': 'CO_2', 

284 'positions': [[0.0, 0.0, 0.0], 

285 [0.0, 0.0, 1.178658], 

286 [0.0, 0.0, -1.178658]], 

287 'symbols': 'COO', 

288 'thermal correction': 2.2321, 

289 'vertical ionization energy': 13.78}, 

290 'CS': {'CAS No.': 2944050, 

291 'ZPE': 1.8242, 

292 'charges': None, 

293 'database': 'G2-1', 

294 'description': 'Carbon monosulfide (CS), C*v symm.', 

295 'enthalpy': 66.9, 

296 'ionization energy': 11.33, 

297 'magmoms': None, 

298 'name': 'SC', 

299 'positions': [[0.0, 0.0, -1.123382], [0.0, 0.0, 0.421268]], 

300 'symbols': 'CS', 

301 'thermal correction': 2.0814}, 

302 'Cl': {'CAS No.': 22537151, 

303 'charges': None, 

304 'database': 'G2-1', 

305 'description': 'Cl atom', 

306 'enthalpy': 28.59, 

307 'ionization energy': 12.97, 

308 'magmoms': [1.0], 

309 'name': 'Chlorine', 

310 'positions': [[0.0, 0.0, 0.0]], 

311 'symbols': 'Cl', 

312 'thermal correction': 1.1}, 

313 'Cl2': {'CAS No.': 7782505, 

314 'ZPE': 0.7737, 

315 'charges': None, 

316 'database': 'G2-1', 

317 'description': 'Cl2 molecule, D*h symm.', 

318 'enthalpy': 0.0, 

319 'ionization energy': 11.48, 

320 'magmoms': None, 

321 'name': 'Cl_2', 

322 'positions': [[0.0, 0.0, 1.007541], [0.0, 0.0, -1.007541]], 

323 'symbols': 'ClCl', 

324 'thermal correction': 2.1963, 

325 'vertical ionization energy': 11.49}, 

326 'ClF': {'CAS No.': 7790898, 

327 'ZPE': 1.1113, 

328 'charges': None, 

329 'database': 'G2-1', 

330 'description': 'ClF molecule, C*v symm, 1-SG.', 

331 'enthalpy': -13.2, 

332 'ionization energy': 12.66, 

333 'magmoms': None, 

334 'name': 'FCl', 

335 'positions': [[0.0, 0.0, -1.084794], [0.0, 0.0, 0.574302]], 

336 'symbols': 'FCl', 

337 'thermal correction': 2.1273, 

338 'vertical ionization energy': 12.77}, 

339 'ClO': {'CAS No.': 14989301, 

340 'ZPE': 1.1923, 

341 'charges': None, 

342 'database': 'G2-1', 

343 'description': 'ClO radical, C*v symm, 2-PI.', 

344 'enthalpy': 24.2, 

345 'ionization energy': 10.89, 

346 'magmoms': [1.0, 0.0], 

347 'name': 'ClO', 

348 'positions': [[0.0, 0.0, 0.514172], [0.0, 0.0, -1.092615]], 

349 'symbols': 'ClO', 

350 'thermal correction': 2.1172, 

351 'vertical ionization energy': 11.01}, 

352 'F': {'CAS No.': 14762948, 

353 'charges': None, 

354 'database': 'G2-1', 

355 'description': 'F atom', 

356 'enthalpy': 18.47, 

357 'ionization energy': 17.42, 

358 'magmoms': [1.0], 

359 'name': 'Fluorine', 

360 'positions': [[0.0, 0.0, 0.0]], 

361 'symbols': 'F', 

362 'thermal correction': 1.05}, 

363 'F2': {'CAS No.': 7782414, 

364 'ZPE': 1.5179, 

365 'charges': None, 

366 'database': 'G2-1', 

367 'description': 'F2 molecule, D*h symm.', 

368 'enthalpy': 0.0, 

369 'ionization energy': 15.7, 

370 'magmoms': None, 

371 'name': 'F_2', 

372 'positions': [[0.0, 0.0, 0.710304], [0.0, 0.0, -0.710304]], 

373 'symbols': 'FF', 

374 'thermal correction': 2.0915, 

375 'vertical ionization energy': 15.7}, 

376 'H': {'CAS No.': 12385136, 

377 'charges': None, 

378 'database': 'G2-1', 

379 'description': 'H atom', 

380 'enthalpy': 51.63, 

381 'ionization energy': 13.6, 

382 'magmoms': [1.0], 

383 'name': 'Hydrogen', 

384 'positions': [[0.0, 0.0, 0.0]], 

385 'symbols': 'H', 

386 'thermal correction': 1.01}, 

387 'H2CO': {'CAS No.': 50000, 

388 'ZPE': 16.4502, 

389 'charges': None, 

390 'database': 'G2-1', 

391 'description': 'Formaldehyde (H2C=O), C2v symm.', 

392 'enthalpy': -26.0, 

393 'ionization energy': 10.88, 

394 'magmoms': None, 

395 'name': 'H_2CO', 

396 'positions': [[0.0, 0.0, 0.683501], 

397 [0.0, 0.0, -0.536614], 

398 [0.0, 0.93439, -1.124164], 

399 [0.0, -0.93439, -1.124164]], 

400 'symbols': 'OCHH', 

401 'thermal correction': 2.3927, 

402 'vertical ionization energy': 10.88}, 

403 'H2O': {'CAS No.': 7732185, 

404 'ZPE': 13.2179, 

405 'charges': None, 

406 'database': 'G2-1', 

407 'description': 'Water (H2O), C2v symm.', 

408 'enthalpy': -57.8, 

409 'ionization energy': 12.62, 

410 'magmoms': None, 

411 'name': 'H_2O', 

412 'positions': [[0.0, 0.0, 0.119262], 

413 [0.0, 0.763239, -0.477047], 

414 [0.0, -0.763239, -0.477047]], 

415 'symbols': 'OHH', 

416 'thermal correction': 2.372}, 

417 'H2O2': {'CAS No.': 7722841, 

418 'ZPE': 16.4081, 

419 'charges': None, 

420 'database': 'G2-1', 

421 'description': 'Hydrogen peroxide (HO-OH), C2 symm.', 

422 'enthalpy': -32.5, 

423 'ionization energy': 10.58, 

424 'magmoms': None, 

425 'name': 'HOOH', 

426 'positions': [[0.0, 0.734058, -0.05275], 

427 [0.0, -0.734058, -0.05275], 

428 [0.839547, 0.880752, 0.422001], 

429 [-0.839547, -0.880752, 0.422001]], 

430 'symbols': 'OOHH', 

431 'thermal correction': 2.623, 

432 'vertical ionization energy': 11.7}, 

433 'HCN': {'CAS No.': 74908, 

434 'ZPE': 10.2654, 

435 'charges': None, 

436 'database': 'G2-1', 

437 'description': 'Hydrogen cyanide (HCN), C*v symm.', 

438 'enthalpy': 31.5, 

439 'ionization energy': 13.6, 

440 'magmoms': None, 

441 'name': 'HCN', 

442 'positions': [[0.0, 0.0, -0.511747], 

443 [0.0, 0.0, 0.664461], 

444 [0.0, 0.0, -1.580746]], 

445 'symbols': 'CNH', 

446 'thermal correction': 2.1768, 

447 'vertical ionization energy': 13.61}, 

448 'HCO': {'CAS No.': 2597446, 

449 'ZPE': 8.029, 

450 'charges': None, 

451 'database': 'G2-1', 

452 'description': 'HCO radical, Bent Cs symm.', 

453 'enthalpy': 10.0, 

454 'ionization energy': 8.12, 

455 'magmoms': [1.0, 0.0, 0.0], 

456 'name': 'HCO', 

457 'positions': [[0.06256, 0.593926, 0.0], 

458 [0.06256, -0.596914, 0.0], 

459 [-0.875835, 1.211755, 0.0]], 

460 'symbols': 'COH', 

461 'thermal correction': 2.3864, 

462 'vertical ionization energy': 9.31}, 

463 'HCl': {'CAS No.': 7647010, 

464 'ZPE': 4.1673, 

465 'charges': None, 

466 'database': 'G2-1', 

467 'description': 'Hydrogen chloride (HCl), C*v symm.', 

468 'enthalpy': -22.1, 

469 'ionization energy': 12.74, 

470 'magmoms': None, 

471 'name': 'HCl', 

472 'positions': [[0.0, 0.0, 0.07111], [0.0, 0.0, -1.208868]], 

473 'symbols': 'ClH', 

474 'thermal correction': 2.0739}, 

475 'HF': {'CAS No.': 7664393, 

476 'ZPE': 5.7994, 

477 'charges': None, 

478 'database': 'G2-1', 

479 'description': 'Hydrogen fluoride (HF), C*v symm.', 

480 'enthalpy': -65.1, 

481 'ionization energy': 16.03, 

482 'magmoms': None, 

483 'name': 'HF', 

484 'positions': [[0.0, 0.0, 0.093389], [0.0, 0.0, -0.840502]], 

485 'symbols': 'FH', 

486 'thermal correction': 2.0733, 

487 'vertical ionization energy': 16.12}, 

488 'HOCl': {'CAS No.': 7790923, 

489 'ZPE': 8.1539, 

490 'charges': None, 

491 'database': 'G2-1', 

492 'description': 'HOCl molecule, Cs symm.', 

493 'enthalpy': -17.8, 

494 'ionization energy': 11.12, 

495 'magmoms': None, 

496 'name': 'HOCl (hypochlorous acid)', 

497 'positions': [[0.036702, 1.113517, 0.0], 

498 [-0.917548, 1.328879, 0.0], 

499 [0.036702, -0.602177, 0.0]], 

500 'symbols': 'OHCl', 

501 'thermal correction': 2.4416}, 

502 'Li': {'CAS No.': 7439932, 

503 'charges': None, 

504 'database': 'G2-1', 

505 'description': 'Li atom', 

506 'enthalpy': 37.69, 

507 'ionization energy': 5.39, 

508 'magmoms': [1.0], 

509 'name': 'Lithium', 

510 'positions': [[0.0, 0.0, 0.0]], 

511 'symbols': 'Li', 

512 'thermal correction': 1.1}, 

513 'Li2': {'CAS No.': 14452596, 

514 'ZPE': 0.4838, 

515 'charges': None, 

516 'database': 'G2-1', 

517 'description': 'Dilithium (Li2), D*h symm.', 

518 'enthalpy': 51.6, 

519 'ionization energy': 5.11, 

520 'magmoms': None, 

521 'name': 'Li_2', 

522 'positions': [[0.0, 0.0, 1.38653], [0.0, 0.0, -1.38653]], 

523 'symbols': 'LiLi', 

524 'thermal correction': 2.3086}, 

525 'LiF': {'CAS No.': 7789244, 

526 'ZPE': 1.4019, 

527 'charges': None, 

528 'database': 'G2-1', 

529 'description': 'Lithium Fluoride (LiF), C*v symm.', 

530 'enthalpy': -80.1, 

531 'ionization energy': 11.3, 

532 'magmoms': None, 

533 'name': 'LiF', 

534 'positions': [[0.0, 0.0, -1.174965], [0.0, 0.0, 0.391655]], 

535 'symbols': 'LiF', 

536 'thermal correction': 2.099}, 

537 'LiH': {'CAS No.': 7580678, 

538 'ZPE': 2.0149, 

539 'charges': None, 

540 'database': 'G2-1', 

541 'description': 'Lithium hydride (LiH), C*v symm.', 

542 'enthalpy': 33.3, 

543 'ionization energy': 7.9, 

544 'magmoms': None, 

545 'name': 'LiH', 

546 'positions': [[0.0, 0.0, 0.41], [0.0, 0.0, -1.23]], 

547 'symbols': 'LiH', 

548 'thermal correction': 2.0783}, 

549 'N': {'CAS No.': 17778880, 

550 'charges': None, 

551 'database': 'G2-1', 

552 'description': 'N atom', 

553 'enthalpy': 112.53, 

554 'ionization energy': 14.53, 

555 'magmoms': [3.0], 

556 'name': 'Nitrogen', 

557 'positions': [[0.0, 0.0, 0.0]], 

558 'symbols': 'N', 

559 'thermal correction': 1.04}, 

560 'N2': {'CAS No.': 7727379, 

561 'ZPE': 3.4243, 

562 'charges': None, 

563 'database': 'G2-1', 

564 'description': 'N2 molecule, D*h symm.', 

565 'enthalpy': 0.0, 

566 'ionization energy': 15.58, 

567 'magmoms': None, 

568 'name': 'N_2', 

569 'positions': [[0.0, 0.0, 0.56499], [0.0, 0.0, -0.56499]], 

570 'symbols': 'NN', 

571 'thermal correction': 2.0733, 

572 'vertical ionization energy': 15.58}, 

573 'N2H4': {'CAS No.': 302012, 

574 'ZPE': 32.9706, 

575 'charges': None, 

576 'database': 'G2-1', 

577 'description': 'Hydrazine (H2N-NH2), C2 symm.', 

578 'enthalpy': 22.8, 

579 'ionization energy': 8.1, 

580 'magmoms': None, 

581 'name': 'H_2NNH_2', 

582 'positions': [[0.0, 0.718959, -0.077687], 

583 [0.0, -0.718959, -0.077687], 

584 [0.211082, 1.092752, 0.847887], 

585 [-0.948214, 1.005026, -0.304078], 

586 [-0.211082, -1.092752, 0.847887], 

587 [0.948214, -1.005026, -0.304078]], 

588 'symbols': 'NNHHHH', 

589 'thermal correction': 2.6531, 

590 'vertical ionization energy': 8.98}, 

591 'NH': {'CAS No.': 13774920, 

592 'ZPE': 4.5739, 

593 'charges': None, 

594 'database': 'G2-1', 

595 'description': 'NH, triplet, C*v symm.', 

596 'enthalpy': 85.2, 

597 'ionization energy': 13.1, 

598 'magmoms': [2.0, 0.0], 

599 'name': 'NH', 

600 'positions': [[0.0, 0.0, 0.129929], [0.0, 0.0, -0.909501]], 

601 'symbols': 'NH', 

602 'thermal correction': 2.0739, 

603 'vertical ionization energy': 13.49}, 

604 'NH2': {'CAS No.': 13770406, 

605 'ZPE': 11.742, 

606 'charges': None, 

607 'database': 'G2-1', 

608 'description': 'NH2 radical, C2v symm, 2-B1.', 

609 'enthalpy': 45.1, 

610 'ionization energy': 10.78, 

611 'magmoms': [1.0, 0.0, 0.0], 

612 'name': 'NH_2', 

613 'positions': [[0.0, 0.0, 0.14169], 

614 [0.0, 0.806442, -0.495913], 

615 [0.0, -0.806442, -0.495913]], 

616 'symbols': 'NHH', 

617 'thermal correction': 2.3726, 

618 'vertical ionization energy': 12.0}, 

619 'NH3': {'CAS No.': 7664417, 

620 'ZPE': 21.2462, 

621 'charges': None, 

622 'database': 'G2-1', 

623 'description': 'Ammonia (NH3), C3v symm.', 

624 'enthalpy': -11.0, 

625 'ionization energy': 10.07, 

626 'magmoms': None, 

627 'name': 'NH_3', 

628 'positions': [[0.0, 0.0, 0.116489], 

629 [0.0, 0.939731, -0.271808], 

630 [0.813831, -0.469865, -0.271808], 

631 [-0.813831, -0.469865, -0.271808]], 

632 'symbols': 'NHHH', 

633 'thermal correction': 2.3896, 

634 'vertical ionization energy': 10.82}, 

635 'NO': {'CAS No.': 10102439, 

636 'ZPE': 2.7974, 

637 'charges': None, 

638 'database': 'G2-1', 

639 'description': 'NO radical, C*v symm, 2-Pi.', 

640 'enthalpy': 21.6, 

641 'ionization energy': 9.26, 

642 'magmoms': [0.6, 0.4], 

643 'name': 'NO', 

644 'positions': [[0.0, 0.0, -0.609442], [0.0, 0.0, 0.533261]], 

645 'symbols': 'NO', 

646 'thermal correction': 2.0745, 

647 'vertical ionization energy': 9.26}, 

648 'Na': {'CAS No.': 7440235, 

649 'charges': None, 

650 'database': 'G2-1', 

651 'description': 'Na atom', 

652 'enthalpy': 25.69, 

653 'ionization energy': 5.14, 

654 'magmoms': [1.0], 

655 'name': 'Sodium', 

656 'positions': [[0.0, 0.0, 0.0]], 

657 'symbols': 'Na', 

658 'thermal correction': 1.54}, 

659 'Na2': {'CAS No.': 25681792, 

660 'ZPE': 0.2246, 

661 'charges': None, 

662 'database': 'G2-1', 

663 'description': 'Disodium (Na2), D*h symm.', 

664 'enthalpy': 34.0, 

665 'ionization energy': 4.89, 

666 'magmoms': None, 

667 'name': 'Na_2', 

668 'positions': [[0.0, 0.0, 1.576262], [0.0, 0.0, -1.576262]], 

669 'symbols': 'NaNa', 

670 'thermal correction': 2.4699}, 

671 'NaCl': {'CAS No.': 7647145, 

672 'ZPE': 0.5152, 

673 'charges': None, 

674 'database': 'G2-1', 

675 'description': 'Sodium Chloride (NaCl), C*v symm.', 

676 'enthalpy': -43.6, 

677 'ionization energy': 9.2, 

678 'magmoms': None, 

679 'name': 'NaCl', 

680 'positions': [[0.0, 0.0, -1.45166], [0.0, 0.0, 0.93931]], 

681 'symbols': 'NaCl', 

682 'thermal correction': 2.2935, 

683 'vertical ionization energy': 9.8}, 

684 'O': {'CAS No.': 17778802, 

685 'charges': None, 

686 'database': 'G2-1', 

687 'description': 'O atom', 

688 'enthalpy': 58.99, 

689 'ionization energy': 13.62, 

690 'magmoms': [2.0], 

691 'name': 'Oxygen', 

692 'positions': [[0.0, 0.0, 0.0]], 

693 'symbols': 'O', 

694 'thermal correction': 1.04}, 

695 'O2': {'CAS No.': 7782447, 

696 'ZPE': 2.3444, 

697 'charges': None, 

698 'database': 'G2-1', 

699 'description': 'O2 molecule, D*h symm, Triplet.', 

700 'enthalpy': 0.0, 

701 'ionization energy': 12.07, 

702 'magmoms': [1.0, 1.0], 

703 'name': 'O_2', 

704 'positions': [[0.0, 0.0, 0.622978], [0.0, 0.0, -0.622978]], 

705 'symbols': 'OO', 

706 'thermal correction': 2.0752, 

707 'vertical ionization energy': 12.3}, 

708 'OH': {'CAS No.': 3352576, 

709 'ZPE': 5.2039, 

710 'charges': None, 

711 'database': 'G2-1', 

712 'description': 'OH radical, C*v symm.', 

713 'enthalpy': 9.4, 

714 'ionization energy': 13.02, 

715 'magmoms': [0.5, 0.5], 

716 'name': 'OH', 

717 'positions': [[0.0, 0.0, 0.108786], [0.0, 0.0, -0.870284]], 

718 'symbols': 'OH', 

719 'thermal correction': 2.0739}, 

720 'P': {'CAS No.': 7723140, 

721 'charges': None, 

722 'database': 'G2-1', 

723 'description': 'P atom', 

724 'enthalpy': 75.42, 

725 'ionization energy': 10.49, 

726 'magmoms': [3.0], 

727 'name': 'Phosphorus', 

728 'positions': [[0.0, 0.0, 0.0]], 

729 'symbols': 'P', 

730 'thermal correction': 1.28}, 

731 'P2': {'CAS No.': 12185090, 

732 'ZPE': 1.1358, 

733 'charges': None, 

734 'database': 'G2-1', 

735 'description': 'P2 molecule, D*h symm.', 

736 'enthalpy': 34.3, 

737 'ionization energy': 10.53, 

738 'magmoms': None, 

739 'name': 'P_2', 

740 'positions': [[0.0, 0.0, 0.966144], [0.0, 0.0, -0.966144]], 

741 'symbols': 'PP', 

742 'thermal correction': 2.1235, 

743 'vertical ionization energy': 10.62}, 

744 'PH2': {'CAS No.': 13765430, 

745 'ZPE': 8.2725, 

746 'charges': None, 

747 'database': 'G2-1', 

748 'description': 'PH2 radical, C2v symm.', 

749 'enthalpy': 33.1, 

750 'ionization energy': 9.82, 

751 'magmoms': [1.0, 0.0, 0.0], 

752 'name': 'PH_2 (Phosphino radical)', 

753 'positions': [[0.0, 0.0, 0.115396], 

754 [0.0, 1.025642, -0.865468], 

755 [0.0, -1.025642, -0.865468]], 

756 'symbols': 'PHH', 

757 'thermal correction': 2.3845}, 

758 'PH3': {'CAS No.': 7803512, 

759 'ZPE': 14.7885, 

760 'charges': None, 

761 'database': 'G2-1', 

762 'description': 'Phosphine (PH3), C3v symm.', 

763 'enthalpy': 1.3, 

764 'ionization energy': 9.87, 

765 'magmoms': None, 

766 'name': 'PH_3', 

767 'positions': [[0.0, 0.0, 0.124619], 

768 [0.0, 1.200647, -0.623095], 

769 [1.039791, -0.600323, -0.623095], 

770 [-1.039791, -0.600323, -0.623095]], 

771 'symbols': 'PHHH', 

772 'thermal correction': 2.4203, 

773 'vertical ionization energy': 10.95}, 

774 'S': {'CAS No.': 7704349, 

775 'charges': None, 

776 'database': 'G2-1', 

777 'description': 'S atom', 

778 'enthalpy': 65.66, 

779 'ionization energy': 10.36, 

780 'magmoms': [2.0], 

781 'name': 'Sulfur', 

782 'positions': [[0.0, 0.0, 0.0]], 

783 'symbols': 'S', 

784 'thermal correction': 1.05}, 

785 'S2': {'CAS No.': 23550450, 

786 'ZPE': 1.0078, 

787 'charges': None, 

788 'database': 'G2-1', 

789 'description': 'S2 molecule, D*h symm, triplet.', 

790 'enthalpy': 30.7, 

791 'ionization energy': 9.36, 

792 'magmoms': [1.0, 1.0], 

793 'name': 'S_2', 

794 'positions': [[0.0, 0.0, 0.960113], [0.0, 0.0, -0.960113]], 

795 'symbols': 'SS', 

796 'thermal correction': 2.1436, 

797 'vertical ionization energy': 9.55}, 

798 'SH2': {'CAS No.': 7783064, 

799 'ZPE': 9.3129, 

800 'charges': None, 

801 'database': 'G2-1', 

802 'description': 'Hydrogen sulfide (H2S), C2v symm.', 

803 'enthalpy': -4.9, 

804 'ionization energy': 10.46, 

805 'magmoms': None, 

806 'name': 'SH_2', 

807 'positions': [[0.0, 0.0, 0.102135], 

808 [0.0, 0.974269, -0.817083], 

809 [0.0, -0.974269, -0.817083]], 

810 'symbols': 'SHH', 

811 'thermal correction': 2.3808, 

812 'vertical ionization energy': 10.5}, 

813 'SO': {'CAS No.': 13827322, 

814 'ZPE': 1.6158, 

815 'charges': None, 

816 'database': 'G2-1', 

817 'description': 'Sulfur monoxide (SO), C*v symm, triplet.', 

818 'enthalpy': 1.2, 

819 'ionization energy': 11.29, 

820 'magmoms': [1.0, 1.0], 

821 'name': 'SO', 

822 'positions': [[0.0, 0.0, -1.015992], [0.0, 0.0, 0.507996]], 

823 'symbols': 'OS', 

824 'thermal correction': 2.0877}, 

825 'SO2': {'CAS No.': 7446095, 

826 'ZPE': 4.3242, 

827 'charges': None, 

828 'database': 'G2-1', 

829 'description': 'Sulfur dioxide (SO2), C2v symm.', 

830 'enthalpy': -71.0, 

831 'ionization energy': 12.35, 

832 'magmoms': None, 

833 'name': 'SO_2', 

834 'positions': [[0.0, 0.0, 0.370268], 

835 [0.0, 1.277617, -0.370268], 

836 [0.0, -1.277617, -0.370268]], 

837 'symbols': 'SOO', 

838 'thermal correction': 2.5245, 

839 'vertical ionization energy': 12.5}, 

840 'Si': {'CAS No.': 7440213, 

841 'charges': None, 

842 'database': 'G2-1', 

843 'description': 'Si atom', 

844 'enthalpy': 106.6, 

845 'ionization energy': 8.15, 

846 'magmoms': [2.0], 

847 'name': 'Silicon', 

848 'positions': [[0.0, 0.0, 0.0]], 

849 'symbols': 'Si', 

850 'thermal correction': 0.76}, 

851 'Si2': {'CAS No.': 12597352, 

852 'ZPE': 0.7028, 

853 'charges': None, 

854 'database': 'G2-1', 

855 'description': 'Si2 molecule, D*h symm, Triplet (3-Sigma-G-).', 

856 'enthalpy': 139.9, 

857 'ionization energy': 7.9, 

858 'magmoms': [1.0, 1.0], 

859 'name': 'Si_2 (Silicon diatomic)', 

860 'positions': [[0.0, 0.0, 1.130054], [0.0, 0.0, -1.130054]], 

861 'symbols': 'SiSi', 

862 'thermal correction': 2.2182}, 

863 'Si2H6': {'CAS No.': 1590870, 

864 'ZPE': 30.2265, 

865 'charges': None, 

866 'database': 'G2-1', 

867 'description': 'Disilane (H3Si-SiH3), D3d symm.', 

868 'enthalpy': 19.1, 

869 'ionization energy': 9.74, 

870 'magmoms': None, 

871 'name': 'Si_2H_6', 

872 'positions': [[0.0, 0.0, 1.167683], 

873 [0.0, 0.0, -1.167683], 

874 [0.0, 1.393286, 1.68602], 

875 [-1.206621, -0.696643, 1.68602], 

876 [1.206621, -0.696643, 1.68602], 

877 [0.0, -1.393286, -1.68602], 

878 [-1.206621, 0.696643, -1.68602], 

879 [1.206621, 0.696643, -1.68602]], 

880 'symbols': 'SiSiHHHHHH', 

881 'thermal correction': 3.7927, 

882 'vertical ionization energy': 10.53}, 

883 'SiH2_s1A1d': {'CAS No.': 13825906, 

884 'ZPE': 7.1875, 

885 'charges': None, 

886 'database': 'G2-1', 

887 'description': 'Singlet silylene (SiH2), C2v symm, 1-A1.', 

888 'enthalpy': 65.2, 

889 'ionization energy': 8.92, 

890 'magmoms': None, 

891 'name': 'SiH_2 (^1A_1)(silicon dihydride)', 

892 'positions': [[0.0, 0.0, 0.131272], 

893 [0.0, 1.096938, -0.918905], 

894 [0.0, -1.096938, -0.918905]], 

895 'symbols': 'SiHH', 

896 'thermal correction': 2.3927}, 

897 'SiH2_s3B1d': {'CAS No.': 13825906, 

898 'ZPE': 7.4203, 

899 'charges': None, 

900 'database': 'G2-1', 

901 'description': 'Triplet silylene (SiH2), C2v symm, 3-B1.', 

902 'enthalpy': 86.2, 

903 'magmoms': [2.0, 0.0, 0.0], 

904 'name': 'SiH_2 (^3B_1)(silicon dihydride)', 

905 'positions': [[0.0, 0.0, 0.094869], 

906 [0.0, 1.271862, -0.664083], 

907 [0.0, -1.271862, -0.664083]], 

908 'symbols': 'SiHH', 

909 'thermal correction': 2.4078}, 

910 'SiH3': {'CAS No.': 13765441, 

911 'ZPE': 13.0898, 

912 'charges': None, 

913 'database': 'G2-1', 

914 'description': 'Silyl radical (SiH3), C3v symm.', 

915 'enthalpy': 47.9, 

916 'ionization energy': 8.14, 

917 'magmoms': [1.0, 0.0, 0.0, 0.0], 

918 'name': 'SiH_3', 

919 'positions': [[0.0, 0.0, 0.079299], 

920 [0.0, 1.41328, -0.370061], 

921 [1.223937, -0.70664, -0.370061], 

922 [-1.223937, -0.70664, -0.370061]], 

923 'symbols': 'SiHHH', 

924 'thermal correction': 2.4912, 

925 'vertical ionization energy': 8.74}, 

926 'SiH4': {'CAS No.': 7803625, 

927 'ZPE': 19.2664, 

928 'charges': None, 

929 'database': 'G2-1', 

930 'description': 'Silane (SiH4), Td symm.', 

931 'enthalpy': 8.2, 

932 'ionization energy': 11.0, 

933 'magmoms': None, 

934 'name': 'SiH_4', 

935 'positions': [[0.0, 0.0, 0.0], 

936 [0.856135, 0.856135, 0.856135], 

937 [-0.856135, -0.856135, 0.856135], 

938 [-0.856135, 0.856135, -0.856135], 

939 [0.856135, -0.856135, -0.856135]], 

940 'symbols': 'SiHHHH', 

941 'thermal correction': 2.5232, 

942 'vertical ionization energy': 12.3}, 

943 'SiO': {'CAS No.': 10097286, 

944 'ZPE': 1.7859, 

945 'charges': None, 

946 'database': 'G2-1', 

947 'description': 'Silicon monoxide (SiO), C*v symm.', 

948 'enthalpy': -24.6, 

949 'ionization energy': 11.49, 

950 'magmoms': None, 

951 'name': 'SiO', 

952 'positions': [[0.0, 0.0, 0.560846], [0.0, 0.0, -0.98148]], 

953 'symbols': 'SiO', 

954 'thermal correction': 2.0821}} 

955 

956 

957def get_ionization_energy(name, vertical=True): 

958 """Return the experimental ionization energy from the database. 

959 

960 If vertical is True, the vertical ionization energy is returned if 

961 available. 

962 """ 

963 if name not in data: 

964 raise KeyError(f'System {name} not in database.') 

965 elif 'ionization energy' not in data[name]: 

966 raise KeyError(f'No data on ionization energy for system {name}.') 

967 else: 

968 if vertical and 'vertical ionization energy' in data[name]: 

969 return data[name]['vertical ionization energy'] 

970 else: 

971 return data[name]['ionization energy'] 

972 

973 

974def get_atomization_energy(name): 

975 """Determine extrapolated experimental atomization energy from the database. 

976 

977 The atomization energy is extrapolated from experimental heats of 

978 formation at room temperature, using calculated zero-point energies 

979 and thermal corrections. 

980 

981 The atomization energy is returned in kcal/mol = 43.36 meV: 

982 

983 >>> from ase.units import *; print kcal / mol 

984 0.0433641146392 

985 

986 """ 

987 d = data[name] 

988 e = d['enthalpy'] 

989 z = d['ZPE'] 

990 dh = d['thermal correction'] 

991 ae = -e + z + dh 

992 for a in string2symbols(d['symbols']): 

993 h = data[a]['enthalpy'] 

994 dh = data[a]['thermal correction'] 

995 ae += h - dh 

996 return ae