Coverage for /builds/ase/ase/ase/data/g2_1.py: 26.09%
23 statements
« prev ^ index » next coverage.py v7.5.3, created at 2025-08-02 00:12 +0000
« prev ^ index » next coverage.py v7.5.3, created at 2025-08-02 00:12 +0000
1# fmt: off
3# flake8: noqa
4"""
5The following contains a database of small molecules
7Data for the G2/97 database are from
8Raghavachari, Redfern, and Pople, J. Chem. Phys. Vol. 106, 1063 (1997).
9See http://www.cse.anl.gov/Catalysis_and_Energy_Conversion/Computational_Thermochemistry.shtml for the original files.
11All numbers are experimental values, except for coordinates, which are
12MP2(full)/6-31G(d) optimized geometries (from http://www.cse.anl.gov/OldCHMwebsiteContent/compmat/G2-97.htm)
14Atomic species:
15ref: Curtiss et al. JCP 106, 1063 (1997).
16'Enthalpy' is the experimental enthalpies of formation at 0K
17'thermal correction' is the thermal corrections H(298)-H(0)
19Molecular species:
20ref: Staroverov et al. JCP 119, 12129 (2003)
21'Enthalpy' is the experimental enthalpies of formation at 298K
22'ZPE' is the zero-point energies
23'thermal correction' is the thermal enthalpy corrections H(298K) - H_exp(0K)
24ZPE and thermal corrections are estimated from B3LYP geometries and vibrations.
26For details about G2-1 and G2-2 sets see doi:10.1063/1.477422.
28Experimental ionization potentials are from http://srdata.nist.gov/cccbdb/
29Information presented on these pages is considered public information
30and may be distributed or copied http://www.nist.gov/public_affairs/disclaimer.cfm
32"""
34from ase.symbols import string2symbols
36atom_names = ['H', 'Li', 'Be', 'C', 'N', 'O', 'F', 'Na', 'Si', 'P', 'S', 'Cl']
38molecule_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',
39 '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']
41data = {
42 'Be': {'CAS No.': 7440417,
43 'charges': None,
44 'database': 'G2-1',
45 'description': 'Be atom',
46 'enthalpy': 76.48,
47 'ionization energy': 9.32,
48 'magmoms': None,
49 'name': 'Beryllium',
50 'positions': [[0.0, 0.0, 0.0]],
51 'symbols': 'Be',
52 'thermal correction': 0.46},
53 'BeH': {'CAS No.': 13597972,
54 'ZPE': 2.9073,
55 'charges': None,
56 'database': 'G2-1',
57 'description': 'Beryllium hydride (BeH), D*h symm.',
58 'enthalpy': 81.7,
59 'ionization energy': 8.21,
60 'magmoms': [0.8, 0.2],
61 'name': 'BeH (beryllium monohydride)',
62 'positions': [[0.0, 0.0, 0.269654], [0.0, 0.0, -1.078616]],
63 'symbols': 'BeH',
64 'thermal correction': 2.0739},
65 'C': {'CAS No.': 7440440,
66 'charges': None,
67 'database': 'G2-1',
68 'description': 'C atom',
69 'enthalpy': 169.98,
70 'ionization energy': 11.26,
71 'magmoms': [2.0],
72 'name': 'Carbon',
73 'positions': [[0.0, 0.0, 0.0]],
74 'symbols': 'C',
75 'thermal correction': 0.25},
76 'C2H2': {'CAS No.': 74862,
77 'ZPE': 16.6001,
78 'charges': None,
79 'database': 'G2-1',
80 'description': 'Acetylene (C2H2), D*h symm.',
81 'enthalpy': 54.2,
82 'ionization energy': 11.4,
83 'magmoms': None,
84 'name': 'C_2H_2',
85 'positions': [[0.0, 0.0, 0.60808],
86 [0.0, 0.0, -0.60808],
87 [0.0, 0.0, -1.67399],
88 [0.0, 0.0, 1.67399]],
89 'symbols': 'CCHH',
90 'thermal correction': 2.4228,
91 'vertical ionization energy': 11.49},
92 'C2H4': {'CAS No.': 74851,
93 'ZPE': 31.5267,
94 'charges': None,
95 'database': 'G2-1',
96 'description': 'Ethylene (H2C=CH2), D2h symm.',
97 'enthalpy': 12.5,
98 'ionization energy': 11.4,
99 'magmoms': None,
100 'name': 'C_2H_4',
101 'positions': [[0.0, 0.0, 0.66748],
102 [0.0, 0.0, -0.66748],
103 [0.0, 0.922832, 1.237695],
104 [0.0, -0.922832, 1.237695],
105 [0.0, 0.922832, -1.237695],
106 [0.0, -0.922832, -1.237695]],
107 'symbols': 'CCHHHH',
108 'thermal correction': 2.51,
109 'vertical ionization energy': 11.49},
110 'C2H6': {'CAS No.': 74840,
111 'ZPE': 46.095,
112 'charges': None,
113 'database': 'G2-1',
114 'description': 'Ethane (H3C-CH3), D3d symm.',
115 'enthalpy': -20.1,
116 'magmoms': None,
117 'name': 'C_2H_6',
118 'positions': [[0.0, 0.0, 0.762209],
119 [0.0, 0.0, -0.762209],
120 [0.0, 1.018957, 1.157229],
121 [-0.882443, -0.509479, 1.157229],
122 [0.882443, -0.509479, 1.157229],
123 [0.0, -1.018957, -1.157229],
124 [-0.882443, 0.509479, -1.157229],
125 [0.882443, 0.509479, -1.157229]],
126 'symbols': 'CCHHHHHH',
127 'thermal correction': 2.7912},
128 'CH': {'CAS No.': 3315375,
129 'ZPE': 3.9659,
130 'charges': None,
131 'database': 'G2-1',
132 'description': 'CH radical. Doublet, C*v symm.',
133 'enthalpy': 142.5,
134 'ionization energy': 10.64,
135 'magmoms': [1.0, 0.0],
136 'name': 'CH(Methylidyne)',
137 'positions': [[0.0, 0.0, 0.160074], [0.0, 0.0, -0.960446]],
138 'symbols': 'CH',
139 'thermal correction': 2.0739},
140 'CH2_s1A1d': {'CAS No.': 2465567,
141 'ZPE': 10.2422,
142 'charges': None,
143 'database': 'G2-1',
144 'description': 'Singlet methylene (CH2), C2v symm, 1-A1.',
145 'enthalpy': 102.8,
146 'magmoms': None,
147 'name': 'CH_2 (^1A_1)',
148 'positions': [[0.0, 0.0, 0.174343],
149 [0.0, 0.862232, -0.523029],
150 [0.0, -0.862232, -0.523029]],
151 'symbols': 'CHH',
152 'thermal correction': 2.3745},
153 'CH2_s3B1d': {'CAS No.': 2465567,
154 'ZPE': 10.6953,
155 'charges': None,
156 'database': 'G2-1',
157 'description': 'Triplet methylene (CH2), C2v symm, 3-B1.',
158 'enthalpy': 93.7,
159 'ionization energy': 10.4,
160 'magmoms': [2.0, 0.0, 0.0],
161 'name': 'CH_2 (^3B_1)',
162 'positions': [[0.0, 0.0, 0.110381],
163 [0.0, 0.982622, -0.331142],
164 [0.0, -0.982622, -0.331142]],
165 'symbols': 'CHH',
166 'thermal correction': 2.3877},
167 'CH3': {'CAS No.': 2229074,
168 'ZPE': 18.3383,
169 'charges': None,
170 'database': 'G2-1',
171 'description': 'Methyl radical (CH3), D3h symm.',
172 'enthalpy': 35.0,
173 'ionization energy': 9.84,
174 'magmoms': [1.0, 0.0, 0.0, 0.0],
175 'name': 'CH_3',
176 'positions': [[0.0, 0.0, 0.0],
177 [0.0, 1.07841, 0.0],
178 [0.93393, -0.539205, 0.0],
179 [-0.93393, -0.539205, 0.0]],
180 'symbols': 'CHHH',
181 'thermal correction': 2.5383},
182 'CH3Cl': {'CAS No.': 74873,
183 'ZPE': 23.3013,
184 'charges': None,
185 'database': 'G2-1',
186 'description': 'Methyl chloride (CH3Cl), C3v symm.',
187 'enthalpy': -19.6,
188 'ionization energy': 11.26,
189 'magmoms': None,
190 'name': 'CH_3Cl',
191 'positions': [[0.0, 0.0, -1.121389],
192 [0.0, 0.0, 0.655951],
193 [0.0, 1.029318, -1.47428],
194 [0.891415, -0.514659, -1.47428],
195 [-0.891415, -0.514659, -1.47428]],
196 'symbols': 'CClHHH',
197 'thermal correction': 2.4956,
198 'vertical ionization energy': 11.29},
199 'CH3OH': {'CAS No.': 67561,
200 'ZPE': 31.6635,
201 'charges': None,
202 'database': 'G2-1',
203 'description': 'Methanol (CH3-OH), Cs symm.',
204 'enthalpy': -48.0,
205 'ionization energy': 10.84,
206 'magmoms': None,
207 'name': 'H_3COH',
208 'positions': [[-0.047131, 0.664389, 0.0],
209 [-0.047131, -0.758551, 0.0],
210 [-1.092995, 0.969785, 0.0],
211 [0.878534, -1.048458, 0.0],
212 [0.437145, 1.080376, 0.891772],
213 [0.437145, 1.080376, -0.891772]],
214 'symbols': 'COHHHH',
215 'thermal correction': 2.6832,
216 'vertical ionization energy': 10.96},
217 'CH3SH': {'CAS No.': 74931,
218 'ZPE': 28.3973,
219 'charges': None,
220 'database': 'G2-1',
221 'description': 'Methanethiol (H3C-SH), Staggered, Cs symm.',
222 'enthalpy': -5.5,
223 'ionization energy': 9.44,
224 'magmoms': None,
225 'name': 'H_3CSH',
226 'positions': [[-0.047953, 1.149519, 0.0],
227 [-0.047953, -0.664856, 0.0],
228 [1.283076, -0.823249, 0.0],
229 [-1.092601, 1.461428, 0.0],
230 [0.432249, 1.551207, 0.892259],
231 [0.432249, 1.551207, -0.892259]],
232 'symbols': 'CSHHHH',
233 'thermal correction': 2.869,
234 'vertical ionization energy': 9.44},
235 'CH4': {'CAS No.': 74828,
236 'ZPE': 27.6744,
237 'charges': None,
238 'database': 'G2-1',
239 'description': 'Methane (CH4), Td symm.',
240 'enthalpy': -17.9,
241 'ionization energy': 12.64,
242 'magmoms': None,
243 'name': 'CH_4',
244 'positions': [[0.0, 0.0, 0.0],
245 [0.629118, 0.629118, 0.629118],
246 [-0.629118, -0.629118, 0.629118],
247 [0.629118, -0.629118, -0.629118],
248 [-0.629118, 0.629118, -0.629118]],
249 'symbols': 'CHHHH',
250 'thermal correction': 2.3939,
251 'vertical ionization energy': 13.6},
252 'CN': {'CAS No.': 2074875,
253 'ZPE': 3.0183,
254 'charges': None,
255 'database': 'G2-1',
256 'description': 'Cyano radical (CN), C*v symm, 2-Sigma+.',
257 'enthalpy': 104.9,
258 'ionization energy': 13.6,
259 'magmoms': [1.0, 0.0],
260 'name': 'CN (Cyano radical)',
261 'positions': [[0.0, 0.0, -0.611046], [0.0, 0.0, 0.523753]],
262 'symbols': 'CN',
263 'thermal correction': 2.0739},
264 'CO': {'CAS No.': 630080,
265 'ZPE': 3.1062,
266 'charges': None,
267 'database': 'G2-1',
268 'description': 'Carbon monoxide (CO), C*v symm.',
269 'enthalpy': -26.4,
270 'ionization energy': 14.01,
271 'magmoms': None,
272 'name': 'CO',
273 'positions': [[0.0, 0.0, 0.493003], [0.0, 0.0, -0.657337]],
274 'symbols': 'OC',
275 'thermal correction': 2.0739,
276 'vertical ionization energy': 14.01},
277 'CO2': {'CAS No.': 124389,
278 'ZPE': 7.313,
279 'charges': None,
280 'database': 'G2-1',
281 'description': 'Carbon dioxide (CO2), D*h symm.',
282 'enthalpy': -94.1,
283 'ionization energy': 13.78,
284 'magmoms': None,
285 'name': 'CO_2',
286 'positions': [[0.0, 0.0, 0.0],
287 [0.0, 0.0, 1.178658],
288 [0.0, 0.0, -1.178658]],
289 'symbols': 'COO',
290 'thermal correction': 2.2321,
291 'vertical ionization energy': 13.78},
292 'CS': {'CAS No.': 2944050,
293 'ZPE': 1.8242,
294 'charges': None,
295 'database': 'G2-1',
296 'description': 'Carbon monosulfide (CS), C*v symm.',
297 'enthalpy': 66.9,
298 'ionization energy': 11.33,
299 'magmoms': None,
300 'name': 'SC',
301 'positions': [[0.0, 0.0, -1.123382], [0.0, 0.0, 0.421268]],
302 'symbols': 'CS',
303 'thermal correction': 2.0814},
304 'Cl': {'CAS No.': 22537151,
305 'charges': None,
306 'database': 'G2-1',
307 'description': 'Cl atom',
308 'enthalpy': 28.59,
309 'ionization energy': 12.97,
310 'magmoms': [1.0],
311 'name': 'Chlorine',
312 'positions': [[0.0, 0.0, 0.0]],
313 'symbols': 'Cl',
314 'thermal correction': 1.1},
315 'Cl2': {'CAS No.': 7782505,
316 'ZPE': 0.7737,
317 'charges': None,
318 'database': 'G2-1',
319 'description': 'Cl2 molecule, D*h symm.',
320 'enthalpy': 0.0,
321 'ionization energy': 11.48,
322 'magmoms': None,
323 'name': 'Cl_2',
324 'positions': [[0.0, 0.0, 1.007541], [0.0, 0.0, -1.007541]],
325 'symbols': 'ClCl',
326 'thermal correction': 2.1963,
327 'vertical ionization energy': 11.49},
328 'ClF': {'CAS No.': 7790898,
329 'ZPE': 1.1113,
330 'charges': None,
331 'database': 'G2-1',
332 'description': 'ClF molecule, C*v symm, 1-SG.',
333 'enthalpy': -13.2,
334 'ionization energy': 12.66,
335 'magmoms': None,
336 'name': 'FCl',
337 'positions': [[0.0, 0.0, -1.084794], [0.0, 0.0, 0.574302]],
338 'symbols': 'FCl',
339 'thermal correction': 2.1273,
340 'vertical ionization energy': 12.77},
341 'ClO': {'CAS No.': 14989301,
342 'ZPE': 1.1923,
343 'charges': None,
344 'database': 'G2-1',
345 'description': 'ClO radical, C*v symm, 2-PI.',
346 'enthalpy': 24.2,
347 'ionization energy': 10.89,
348 'magmoms': [1.0, 0.0],
349 'name': 'ClO',
350 'positions': [[0.0, 0.0, 0.514172], [0.0, 0.0, -1.092615]],
351 'symbols': 'ClO',
352 'thermal correction': 2.1172,
353 'vertical ionization energy': 11.01},
354 'F': {'CAS No.': 14762948,
355 'charges': None,
356 'database': 'G2-1',
357 'description': 'F atom',
358 'enthalpy': 18.47,
359 'ionization energy': 17.42,
360 'magmoms': [1.0],
361 'name': 'Fluorine',
362 'positions': [[0.0, 0.0, 0.0]],
363 'symbols': 'F',
364 'thermal correction': 1.05},
365 'F2': {'CAS No.': 7782414,
366 'ZPE': 1.5179,
367 'charges': None,
368 'database': 'G2-1',
369 'description': 'F2 molecule, D*h symm.',
370 'enthalpy': 0.0,
371 'ionization energy': 15.7,
372 'magmoms': None,
373 'name': 'F_2',
374 'positions': [[0.0, 0.0, 0.710304], [0.0, 0.0, -0.710304]],
375 'symbols': 'FF',
376 'thermal correction': 2.0915,
377 'vertical ionization energy': 15.7},
378 'H': {'CAS No.': 12385136,
379 'charges': None,
380 'database': 'G2-1',
381 'description': 'H atom',
382 'enthalpy': 51.63,
383 'ionization energy': 13.6,
384 'magmoms': [1.0],
385 'name': 'Hydrogen',
386 'positions': [[0.0, 0.0, 0.0]],
387 'symbols': 'H',
388 'thermal correction': 1.01},
389 'H2CO': {'CAS No.': 50000,
390 'ZPE': 16.4502,
391 'charges': None,
392 'database': 'G2-1',
393 'description': 'Formaldehyde (H2C=O), C2v symm.',
394 'enthalpy': -26.0,
395 'ionization energy': 10.88,
396 'magmoms': None,
397 'name': 'H_2CO',
398 'positions': [[0.0, 0.0, 0.683501],
399 [0.0, 0.0, -0.536614],
400 [0.0, 0.93439, -1.124164],
401 [0.0, -0.93439, -1.124164]],
402 'symbols': 'OCHH',
403 'thermal correction': 2.3927,
404 'vertical ionization energy': 10.88},
405 'H2O': {'CAS No.': 7732185,
406 'ZPE': 13.2179,
407 'charges': None,
408 'database': 'G2-1',
409 'description': 'Water (H2O), C2v symm.',
410 'enthalpy': -57.8,
411 'ionization energy': 12.62,
412 'magmoms': None,
413 'name': 'H_2O',
414 'positions': [[0.0, 0.0, 0.119262],
415 [0.0, 0.763239, -0.477047],
416 [0.0, -0.763239, -0.477047]],
417 'symbols': 'OHH',
418 'thermal correction': 2.372},
419 'H2O2': {'CAS No.': 7722841,
420 'ZPE': 16.4081,
421 'charges': None,
422 'database': 'G2-1',
423 'description': 'Hydrogen peroxide (HO-OH), C2 symm.',
424 'enthalpy': -32.5,
425 'ionization energy': 10.58,
426 'magmoms': None,
427 'name': 'HOOH',
428 'positions': [[0.0, 0.734058, -0.05275],
429 [0.0, -0.734058, -0.05275],
430 [0.839547, 0.880752, 0.422001],
431 [-0.839547, -0.880752, 0.422001]],
432 'symbols': 'OOHH',
433 'thermal correction': 2.623,
434 'vertical ionization energy': 11.7},
435 'HCN': {'CAS No.': 74908,
436 'ZPE': 10.2654,
437 'charges': None,
438 'database': 'G2-1',
439 'description': 'Hydrogen cyanide (HCN), C*v symm.',
440 'enthalpy': 31.5,
441 'ionization energy': 13.6,
442 'magmoms': None,
443 'name': 'HCN',
444 'positions': [[0.0, 0.0, -0.511747],
445 [0.0, 0.0, 0.664461],
446 [0.0, 0.0, -1.580746]],
447 'symbols': 'CNH',
448 'thermal correction': 2.1768,
449 'vertical ionization energy': 13.61},
450 'HCO': {'CAS No.': 2597446,
451 'ZPE': 8.029,
452 'charges': None,
453 'database': 'G2-1',
454 'description': 'HCO radical, Bent Cs symm.',
455 'enthalpy': 10.0,
456 'ionization energy': 8.12,
457 'magmoms': [1.0, 0.0, 0.0],
458 'name': 'HCO',
459 'positions': [[0.06256, 0.593926, 0.0],
460 [0.06256, -0.596914, 0.0],
461 [-0.875835, 1.211755, 0.0]],
462 'symbols': 'COH',
463 'thermal correction': 2.3864,
464 'vertical ionization energy': 9.31},
465 'HCl': {'CAS No.': 7647010,
466 'ZPE': 4.1673,
467 'charges': None,
468 'database': 'G2-1',
469 'description': 'Hydrogen chloride (HCl), C*v symm.',
470 'enthalpy': -22.1,
471 'ionization energy': 12.74,
472 'magmoms': None,
473 'name': 'HCl',
474 'positions': [[0.0, 0.0, 0.07111], [0.0, 0.0, -1.208868]],
475 'symbols': 'ClH',
476 'thermal correction': 2.0739},
477 'HF': {'CAS No.': 7664393,
478 'ZPE': 5.7994,
479 'charges': None,
480 'database': 'G2-1',
481 'description': 'Hydrogen fluoride (HF), C*v symm.',
482 'enthalpy': -65.1,
483 'ionization energy': 16.03,
484 'magmoms': None,
485 'name': 'HF',
486 'positions': [[0.0, 0.0, 0.093389], [0.0, 0.0, -0.840502]],
487 'symbols': 'FH',
488 'thermal correction': 2.0733,
489 'vertical ionization energy': 16.12},
490 'HOCl': {'CAS No.': 7790923,
491 'ZPE': 8.1539,
492 'charges': None,
493 'database': 'G2-1',
494 'description': 'HOCl molecule, Cs symm.',
495 'enthalpy': -17.8,
496 'ionization energy': 11.12,
497 'magmoms': None,
498 'name': 'HOCl (hypochlorous acid)',
499 'positions': [[0.036702, 1.113517, 0.0],
500 [-0.917548, 1.328879, 0.0],
501 [0.036702, -0.602177, 0.0]],
502 'symbols': 'OHCl',
503 'thermal correction': 2.4416},
504 'Li': {'CAS No.': 7439932,
505 'charges': None,
506 'database': 'G2-1',
507 'description': 'Li atom',
508 'enthalpy': 37.69,
509 'ionization energy': 5.39,
510 'magmoms': [1.0],
511 'name': 'Lithium',
512 'positions': [[0.0, 0.0, 0.0]],
513 'symbols': 'Li',
514 'thermal correction': 1.1},
515 'Li2': {'CAS No.': 14452596,
516 'ZPE': 0.4838,
517 'charges': None,
518 'database': 'G2-1',
519 'description': 'Dilithium (Li2), D*h symm.',
520 'enthalpy': 51.6,
521 'ionization energy': 5.11,
522 'magmoms': None,
523 'name': 'Li_2',
524 'positions': [[0.0, 0.0, 1.38653], [0.0, 0.0, -1.38653]],
525 'symbols': 'LiLi',
526 'thermal correction': 2.3086},
527 'LiF': {'CAS No.': 7789244,
528 'ZPE': 1.4019,
529 'charges': None,
530 'database': 'G2-1',
531 'description': 'Lithium Fluoride (LiF), C*v symm.',
532 'enthalpy': -80.1,
533 'ionization energy': 11.3,
534 'magmoms': None,
535 'name': 'LiF',
536 'positions': [[0.0, 0.0, -1.174965], [0.0, 0.0, 0.391655]],
537 'symbols': 'LiF',
538 'thermal correction': 2.099},
539 'LiH': {'CAS No.': 7580678,
540 'ZPE': 2.0149,
541 'charges': None,
542 'database': 'G2-1',
543 'description': 'Lithium hydride (LiH), C*v symm.',
544 'enthalpy': 33.3,
545 'ionization energy': 7.9,
546 'magmoms': None,
547 'name': 'LiH',
548 'positions': [[0.0, 0.0, 0.41], [0.0, 0.0, -1.23]],
549 'symbols': 'LiH',
550 'thermal correction': 2.0783},
551 'N': {'CAS No.': 17778880,
552 'charges': None,
553 'database': 'G2-1',
554 'description': 'N atom',
555 'enthalpy': 112.53,
556 'ionization energy': 14.53,
557 'magmoms': [3.0],
558 'name': 'Nitrogen',
559 'positions': [[0.0, 0.0, 0.0]],
560 'symbols': 'N',
561 'thermal correction': 1.04},
562 'N2': {'CAS No.': 7727379,
563 'ZPE': 3.4243,
564 'charges': None,
565 'database': 'G2-1',
566 'description': 'N2 molecule, D*h symm.',
567 'enthalpy': 0.0,
568 'ionization energy': 15.58,
569 'magmoms': None,
570 'name': 'N_2',
571 'positions': [[0.0, 0.0, 0.56499], [0.0, 0.0, -0.56499]],
572 'symbols': 'NN',
573 'thermal correction': 2.0733,
574 'vertical ionization energy': 15.58},
575 'N2H4': {'CAS No.': 302012,
576 'ZPE': 32.9706,
577 'charges': None,
578 'database': 'G2-1',
579 'description': 'Hydrazine (H2N-NH2), C2 symm.',
580 'enthalpy': 22.8,
581 'ionization energy': 8.1,
582 'magmoms': None,
583 'name': 'H_2NNH_2',
584 'positions': [[0.0, 0.718959, -0.077687],
585 [0.0, -0.718959, -0.077687],
586 [0.211082, 1.092752, 0.847887],
587 [-0.948214, 1.005026, -0.304078],
588 [-0.211082, -1.092752, 0.847887],
589 [0.948214, -1.005026, -0.304078]],
590 'symbols': 'NNHHHH',
591 'thermal correction': 2.6531,
592 'vertical ionization energy': 8.98},
593 'NH': {'CAS No.': 13774920,
594 'ZPE': 4.5739,
595 'charges': None,
596 'database': 'G2-1',
597 'description': 'NH, triplet, C*v symm.',
598 'enthalpy': 85.2,
599 'ionization energy': 13.1,
600 'magmoms': [2.0, 0.0],
601 'name': 'NH',
602 'positions': [[0.0, 0.0, 0.129929], [0.0, 0.0, -0.909501]],
603 'symbols': 'NH',
604 'thermal correction': 2.0739,
605 'vertical ionization energy': 13.49},
606 'NH2': {'CAS No.': 13770406,
607 'ZPE': 11.742,
608 'charges': None,
609 'database': 'G2-1',
610 'description': 'NH2 radical, C2v symm, 2-B1.',
611 'enthalpy': 45.1,
612 'ionization energy': 10.78,
613 'magmoms': [1.0, 0.0, 0.0],
614 'name': 'NH_2',
615 'positions': [[0.0, 0.0, 0.14169],
616 [0.0, 0.806442, -0.495913],
617 [0.0, -0.806442, -0.495913]],
618 'symbols': 'NHH',
619 'thermal correction': 2.3726,
620 'vertical ionization energy': 12.0},
621 'NH3': {'CAS No.': 7664417,
622 'ZPE': 21.2462,
623 'charges': None,
624 'database': 'G2-1',
625 'description': 'Ammonia (NH3), C3v symm.',
626 'enthalpy': -11.0,
627 'ionization energy': 10.07,
628 'magmoms': None,
629 'name': 'NH_3',
630 'positions': [[0.0, 0.0, 0.116489],
631 [0.0, 0.939731, -0.271808],
632 [0.813831, -0.469865, -0.271808],
633 [-0.813831, -0.469865, -0.271808]],
634 'symbols': 'NHHH',
635 'thermal correction': 2.3896,
636 'vertical ionization energy': 10.82},
637 'NO': {'CAS No.': 10102439,
638 'ZPE': 2.7974,
639 'charges': None,
640 'database': 'G2-1',
641 'description': 'NO radical, C*v symm, 2-Pi.',
642 'enthalpy': 21.6,
643 'ionization energy': 9.26,
644 'magmoms': [0.6, 0.4],
645 'name': 'NO',
646 'positions': [[0.0, 0.0, -0.609442], [0.0, 0.0, 0.533261]],
647 'symbols': 'NO',
648 'thermal correction': 2.0745,
649 'vertical ionization energy': 9.26},
650 'Na': {'CAS No.': 7440235,
651 'charges': None,
652 'database': 'G2-1',
653 'description': 'Na atom',
654 'enthalpy': 25.69,
655 'ionization energy': 5.14,
656 'magmoms': [1.0],
657 'name': 'Sodium',
658 'positions': [[0.0, 0.0, 0.0]],
659 'symbols': 'Na',
660 'thermal correction': 1.54},
661 'Na2': {'CAS No.': 25681792,
662 'ZPE': 0.2246,
663 'charges': None,
664 'database': 'G2-1',
665 'description': 'Disodium (Na2), D*h symm.',
666 'enthalpy': 34.0,
667 'ionization energy': 4.89,
668 'magmoms': None,
669 'name': 'Na_2',
670 'positions': [[0.0, 0.0, 1.576262], [0.0, 0.0, -1.576262]],
671 'symbols': 'NaNa',
672 'thermal correction': 2.4699},
673 'NaCl': {'CAS No.': 7647145,
674 'ZPE': 0.5152,
675 'charges': None,
676 'database': 'G2-1',
677 'description': 'Sodium Chloride (NaCl), C*v symm.',
678 'enthalpy': -43.6,
679 'ionization energy': 9.2,
680 'magmoms': None,
681 'name': 'NaCl',
682 'positions': [[0.0, 0.0, -1.45166], [0.0, 0.0, 0.93931]],
683 'symbols': 'NaCl',
684 'thermal correction': 2.2935,
685 'vertical ionization energy': 9.8},
686 'O': {'CAS No.': 17778802,
687 'charges': None,
688 'database': 'G2-1',
689 'description': 'O atom',
690 'enthalpy': 58.99,
691 'ionization energy': 13.62,
692 'magmoms': [2.0],
693 'name': 'Oxygen',
694 'positions': [[0.0, 0.0, 0.0]],
695 'symbols': 'O',
696 'thermal correction': 1.04},
697 'O2': {'CAS No.': 7782447,
698 'ZPE': 2.3444,
699 'charges': None,
700 'database': 'G2-1',
701 'description': 'O2 molecule, D*h symm, Triplet.',
702 'enthalpy': 0.0,
703 'ionization energy': 12.07,
704 'magmoms': [1.0, 1.0],
705 'name': 'O_2',
706 'positions': [[0.0, 0.0, 0.622978], [0.0, 0.0, -0.622978]],
707 'symbols': 'OO',
708 'thermal correction': 2.0752,
709 'vertical ionization energy': 12.3},
710 'OH': {'CAS No.': 3352576,
711 'ZPE': 5.2039,
712 'charges': None,
713 'database': 'G2-1',
714 'description': 'OH radical, C*v symm.',
715 'enthalpy': 9.4,
716 'ionization energy': 13.02,
717 'magmoms': [0.5, 0.5],
718 'name': 'OH',
719 'positions': [[0.0, 0.0, 0.108786], [0.0, 0.0, -0.870284]],
720 'symbols': 'OH',
721 'thermal correction': 2.0739},
722 'P': {'CAS No.': 7723140,
723 'charges': None,
724 'database': 'G2-1',
725 'description': 'P atom',
726 'enthalpy': 75.42,
727 'ionization energy': 10.49,
728 'magmoms': [3.0],
729 'name': 'Phosphorus',
730 'positions': [[0.0, 0.0, 0.0]],
731 'symbols': 'P',
732 'thermal correction': 1.28},
733 'P2': {'CAS No.': 12185090,
734 'ZPE': 1.1358,
735 'charges': None,
736 'database': 'G2-1',
737 'description': 'P2 molecule, D*h symm.',
738 'enthalpy': 34.3,
739 'ionization energy': 10.53,
740 'magmoms': None,
741 'name': 'P_2',
742 'positions': [[0.0, 0.0, 0.966144], [0.0, 0.0, -0.966144]],
743 'symbols': 'PP',
744 'thermal correction': 2.1235,
745 'vertical ionization energy': 10.62},
746 'PH2': {'CAS No.': 13765430,
747 'ZPE': 8.2725,
748 'charges': None,
749 'database': 'G2-1',
750 'description': 'PH2 radical, C2v symm.',
751 'enthalpy': 33.1,
752 'ionization energy': 9.82,
753 'magmoms': [1.0, 0.0, 0.0],
754 'name': 'PH_2 (Phosphino radical)',
755 'positions': [[0.0, 0.0, 0.115396],
756 [0.0, 1.025642, -0.865468],
757 [0.0, -1.025642, -0.865468]],
758 'symbols': 'PHH',
759 'thermal correction': 2.3845},
760 'PH3': {'CAS No.': 7803512,
761 'ZPE': 14.7885,
762 'charges': None,
763 'database': 'G2-1',
764 'description': 'Phosphine (PH3), C3v symm.',
765 'enthalpy': 1.3,
766 'ionization energy': 9.87,
767 'magmoms': None,
768 'name': 'PH_3',
769 'positions': [[0.0, 0.0, 0.124619],
770 [0.0, 1.200647, -0.623095],
771 [1.039791, -0.600323, -0.623095],
772 [-1.039791, -0.600323, -0.623095]],
773 'symbols': 'PHHH',
774 'thermal correction': 2.4203,
775 'vertical ionization energy': 10.95},
776 'S': {'CAS No.': 7704349,
777 'charges': None,
778 'database': 'G2-1',
779 'description': 'S atom',
780 'enthalpy': 65.66,
781 'ionization energy': 10.36,
782 'magmoms': [2.0],
783 'name': 'Sulfur',
784 'positions': [[0.0, 0.0, 0.0]],
785 'symbols': 'S',
786 'thermal correction': 1.05},
787 'S2': {'CAS No.': 23550450,
788 'ZPE': 1.0078,
789 'charges': None,
790 'database': 'G2-1',
791 'description': 'S2 molecule, D*h symm, triplet.',
792 'enthalpy': 30.7,
793 'ionization energy': 9.36,
794 'magmoms': [1.0, 1.0],
795 'name': 'S_2',
796 'positions': [[0.0, 0.0, 0.960113], [0.0, 0.0, -0.960113]],
797 'symbols': 'SS',
798 'thermal correction': 2.1436,
799 'vertical ionization energy': 9.55},
800 'SH2': {'CAS No.': 7783064,
801 'ZPE': 9.3129,
802 'charges': None,
803 'database': 'G2-1',
804 'description': 'Hydrogen sulfide (H2S), C2v symm.',
805 'enthalpy': -4.9,
806 'ionization energy': 10.46,
807 'magmoms': None,
808 'name': 'SH_2',
809 'positions': [[0.0, 0.0, 0.102135],
810 [0.0, 0.974269, -0.817083],
811 [0.0, -0.974269, -0.817083]],
812 'symbols': 'SHH',
813 'thermal correction': 2.3808,
814 'vertical ionization energy': 10.5},
815 'SO': {'CAS No.': 13827322,
816 'ZPE': 1.6158,
817 'charges': None,
818 'database': 'G2-1',
819 'description': 'Sulfur monoxide (SO), C*v symm, triplet.',
820 'enthalpy': 1.2,
821 'ionization energy': 11.29,
822 'magmoms': [1.0, 1.0],
823 'name': 'SO',
824 'positions': [[0.0, 0.0, -1.015992], [0.0, 0.0, 0.507996]],
825 'symbols': 'OS',
826 'thermal correction': 2.0877},
827 'SO2': {'CAS No.': 7446095,
828 'ZPE': 4.3242,
829 'charges': None,
830 'database': 'G2-1',
831 'description': 'Sulfur dioxide (SO2), C2v symm.',
832 'enthalpy': -71.0,
833 'ionization energy': 12.35,
834 'magmoms': None,
835 'name': 'SO_2',
836 'positions': [[0.0, 0.0, 0.370268],
837 [0.0, 1.277617, -0.370268],
838 [0.0, -1.277617, -0.370268]],
839 'symbols': 'SOO',
840 'thermal correction': 2.5245,
841 'vertical ionization energy': 12.5},
842 'Si': {'CAS No.': 7440213,
843 'charges': None,
844 'database': 'G2-1',
845 'description': 'Si atom',
846 'enthalpy': 106.6,
847 'ionization energy': 8.15,
848 'magmoms': [2.0],
849 'name': 'Silicon',
850 'positions': [[0.0, 0.0, 0.0]],
851 'symbols': 'Si',
852 'thermal correction': 0.76},
853 'Si2': {'CAS No.': 12597352,
854 'ZPE': 0.7028,
855 'charges': None,
856 'database': 'G2-1',
857 'description': 'Si2 molecule, D*h symm, Triplet (3-Sigma-G-).',
858 'enthalpy': 139.9,
859 'ionization energy': 7.9,
860 'magmoms': [1.0, 1.0],
861 'name': 'Si_2 (Silicon diatomic)',
862 'positions': [[0.0, 0.0, 1.130054], [0.0, 0.0, -1.130054]],
863 'symbols': 'SiSi',
864 'thermal correction': 2.2182},
865 'Si2H6': {'CAS No.': 1590870,
866 'ZPE': 30.2265,
867 'charges': None,
868 'database': 'G2-1',
869 'description': 'Disilane (H3Si-SiH3), D3d symm.',
870 'enthalpy': 19.1,
871 'ionization energy': 9.74,
872 'magmoms': None,
873 'name': 'Si_2H_6',
874 'positions': [[0.0, 0.0, 1.167683],
875 [0.0, 0.0, -1.167683],
876 [0.0, 1.393286, 1.68602],
877 [-1.206621, -0.696643, 1.68602],
878 [1.206621, -0.696643, 1.68602],
879 [0.0, -1.393286, -1.68602],
880 [-1.206621, 0.696643, -1.68602],
881 [1.206621, 0.696643, -1.68602]],
882 'symbols': 'SiSiHHHHHH',
883 'thermal correction': 3.7927,
884 'vertical ionization energy': 10.53},
885 'SiH2_s1A1d': {'CAS No.': 13825906,
886 'ZPE': 7.1875,
887 'charges': None,
888 'database': 'G2-1',
889 'description': 'Singlet silylene (SiH2), C2v symm, 1-A1.',
890 'enthalpy': 65.2,
891 'ionization energy': 8.92,
892 'magmoms': None,
893 'name': 'SiH_2 (^1A_1)(silicon dihydride)',
894 'positions': [[0.0, 0.0, 0.131272],
895 [0.0, 1.096938, -0.918905],
896 [0.0, -1.096938, -0.918905]],
897 'symbols': 'SiHH',
898 'thermal correction': 2.3927},
899 'SiH2_s3B1d': {'CAS No.': 13825906,
900 'ZPE': 7.4203,
901 'charges': None,
902 'database': 'G2-1',
903 'description': 'Triplet silylene (SiH2), C2v symm, 3-B1.',
904 'enthalpy': 86.2,
905 'magmoms': [2.0, 0.0, 0.0],
906 'name': 'SiH_2 (^3B_1)(silicon dihydride)',
907 'positions': [[0.0, 0.0, 0.094869],
908 [0.0, 1.271862, -0.664083],
909 [0.0, -1.271862, -0.664083]],
910 'symbols': 'SiHH',
911 'thermal correction': 2.4078},
912 'SiH3': {'CAS No.': 13765441,
913 'ZPE': 13.0898,
914 'charges': None,
915 'database': 'G2-1',
916 'description': 'Silyl radical (SiH3), C3v symm.',
917 'enthalpy': 47.9,
918 'ionization energy': 8.14,
919 'magmoms': [1.0, 0.0, 0.0, 0.0],
920 'name': 'SiH_3',
921 'positions': [[0.0, 0.0, 0.079299],
922 [0.0, 1.41328, -0.370061],
923 [1.223937, -0.70664, -0.370061],
924 [-1.223937, -0.70664, -0.370061]],
925 'symbols': 'SiHHH',
926 'thermal correction': 2.4912,
927 'vertical ionization energy': 8.74},
928 'SiH4': {'CAS No.': 7803625,
929 'ZPE': 19.2664,
930 'charges': None,
931 'database': 'G2-1',
932 'description': 'Silane (SiH4), Td symm.',
933 'enthalpy': 8.2,
934 'ionization energy': 11.0,
935 'magmoms': None,
936 'name': 'SiH_4',
937 'positions': [[0.0, 0.0, 0.0],
938 [0.856135, 0.856135, 0.856135],
939 [-0.856135, -0.856135, 0.856135],
940 [-0.856135, 0.856135, -0.856135],
941 [0.856135, -0.856135, -0.856135]],
942 'symbols': 'SiHHHH',
943 'thermal correction': 2.5232,
944 'vertical ionization energy': 12.3},
945 'SiO': {'CAS No.': 10097286,
946 'ZPE': 1.7859,
947 'charges': None,
948 'database': 'G2-1',
949 'description': 'Silicon monoxide (SiO), C*v symm.',
950 'enthalpy': -24.6,
951 'ionization energy': 11.49,
952 'magmoms': None,
953 'name': 'SiO',
954 'positions': [[0.0, 0.0, 0.560846], [0.0, 0.0, -0.98148]],
955 'symbols': 'SiO',
956 'thermal correction': 2.0821}}
959def get_ionization_energy(name, vertical=True):
960 """Return the experimental ionization energy from the database.
962 If vertical is True, the vertical ionization energy is returned if
963 available.
964 """
965 if name not in data:
966 raise KeyError(f'System {name} not in database.')
967 elif 'ionization energy' not in data[name]:
968 raise KeyError(f'No data on ionization energy for system {name}.')
969 else:
970 if vertical and 'vertical ionization energy' in data[name]:
971 return data[name]['vertical ionization energy']
972 else:
973 return data[name]['ionization energy']
976def get_atomization_energy(name):
977 """Determine extrapolated experimental atomization energy from the database.
979 The atomization energy is extrapolated from experimental heats of
980 formation at room temperature, using calculated zero-point energies
981 and thermal corrections.
983 The atomization energy is returned in kcal/mol = 43.36 meV:
985 >>> from ase.units import *; print kcal / mol
986 0.0433641146392
988 """
989 d = data[name]
990 e = d['enthalpy']
991 z = d['ZPE']
992 dh = d['thermal correction']
993 ae = -e + z + dh
994 for a in string2symbols(d['symbols']):
995 h = data[a]['enthalpy']
996 dh = data[a]['thermal correction']
997 ae += h - dh
998 return ae