Coverage for gpaw/utilities/acwf.py: 25%

1"""ACWF-benchmark. 

2 

3Do equation-of-state calculation for the 10 reference systems from 

4the AiiDA common workflows (ACWF) benchmark: 

5DIAMOND, FCC, SC, BCC, XO3, XO, X4O6, XO2, X4O10, X2O. 

6""" 

7from __future__ import annotations 

8import json 

9from pathlib import Path 

10from typing import Any 

11 

12import numpy as np 

13from ase import Atoms 

14from ase.data import atomic_numbers 

15from gpaw.atom.check import all_names 

16from gpaw.mpi import world 

17from gpaw.new.ase_interface import GPAW 

18from gpaw.dft import Parameters 

19 

20 

21def eos(atoms: Atoms, 

22 npoints: int = 7, 

23 maxstrain: float = 0.02) -> dict: 

24 """Do EOS calculation and fit result.""" 

25 cell_cv = atoms.get_cell() 

26 

27 energies = [] 

28 strains = np.linspace(-maxstrain, maxstrain, npoints) 

29 for i, strain in enumerate(strains): 

30 atoms.set_cell((1.0 + strain) * cell_cv, scale_atoms=True) 

31 e = atoms.get_potential_energy() 

32 energies.append(e) 

33 

34 strain0, e0, d2eds2 = fit(strains, energies) 

35 v0 = abs(np.linalg.det(cell_cv)) 

36 volume = v0 * (1 + strain0)**3 

37 

38 return dict( 

39 volume=volume, 

40 strains=strains.tolist(), 

41 energies=energies, 

42 energy=e0, 

43 strain=strain0, 

44 d2eds2=d2eds2) 

45 

46 

47def fit(strains: np.ndarray, 

48 energies: np.ndarray | list[float]) -> tuple[float, float, float]: 

49 """Fit energies. 

50 

51 :: 

52 

53 3 

54 --- n 1 

55 f(s) = p(x) = > c x, x = -----, 

56 --- n 1 + s 

57 n=0 

58 

59 df/ds = dp/dx dx/ds, 

60 

61 2 2 2 2 2 

62 d f/ds = d p/dx (dx/ds). 

63 

64 """ 

65 

66 x = 1 / (1 + strains) 

67 p0 = np.poly1d(np.polyfit(x, energies, 3)) 

68 p1 = np.polyder(p0, 1) 

69 p2 = np.polyder(p1, 1) 

70 for x0 in np.roots(p1): 

71 if isinstance(x0, float) and x0 > 0 and p2(x0) > 0: 

72 break 

73 else: # no break 

74 raise ValueError 

75 strain = 1 / x0 - 1 

76 e0 = p0(x0) 

77 d2eds2 = p2(x0) * x0**2 

78 return strain, e0, d2eds2 

79 

80 

81def run_eos_calculation(structure: str, symbol: str, **params) -> dict: 

82 """Create GPAW calculator and do EOS calculation.""" 

83 atoms = reference_structure(symbol, structure) 

84 builder = Parameters(**params).dft_component_builder(atoms, log=None) 

85 params['kpts'] = builder.ibz.bz.size_c # type: ignore 

86 params['gpts'] = builder.grid.size 

87 params.pop('h', None) 

88 atoms.calc = GPAW(txt=f'{builder.mode}-{structure}.txt', **params) 

89 data = eos(atoms) 

90 return data 

91 

92 

93def workflow() -> None: 

94 """Submit PW and LCAO calculations. 

95 

96 101 PAW-potentials, 10 crystals, PW and LCAO: 2020 calculations. 

97 """ 

98 from myqueue.workflow import run 

99 symbol, _, setup_name = Path().absolute().name.partition('.') 

100 for x in acwf_structures: 

101 run(function=work, args=[x, symbol, setup_name, 'pw'], 

102 cores=24, tmax='1h', name=f'pw-{x}') 

103 run(function=work, args=[x, symbol, setup_name, 'lcao'], 

104 cores=24, tmax='5h', name=f'lcao-{x}', restart=2) 

105 

106 

107def work(structure: str, 

108 symbol: str, 

109 setup_name: str = '', 

110 mode: str = 'pw'): 

111 """Do single EOS calculations with PBE.""" 

112 params: dict[str, Any] = dict( 

113 xc='PBE', 

114 occupations=dict(name='fermi-dirac', width=0.0612), 

115 kpts=dict(density=6.0), 

116 convergence=dict(density=1e-6)) 

117 if mode == 'pw': 

118 params['mode'] = dict(name='pw', ecut=1000) 

119 else: 

120 params['mode'] = 'lcao' 

121 params['basis'] = 'dzp' 

122 params['h'] = 0.12 

123 if setup_name: 

124 params['setups'] = {symbol: setup_name} 

125 data = run_eos_calculation( 

126 structure, 

127 symbol, 

128 **params) 

129 if world.rank == 0: 

130 Path(f'{mode}-{structure}.json').write_text( 

131 json.dumps(data, indent=1)) 

132 

133 

134def create_folders() -> None: 

135 """Create folder for all 101 PAW-potentials.""" 

136 for name in all_names: 

137 Path(name).mkdir(exist_ok=True) 

138 

139 

140def collect_data() -> None: 

141 """Collect everything need for a table on the web-page.""" 

142 results = {} 

143 for name in all_names: 

144 result: dict[str, list[tuple[str, float]]] = {'pw': [], 'lcao': []} 

145 for mode in ['pw', 'lcao']: 

146 for path in Path(name).glob(f'{mode}-*.json'): 

147 structure = path.stem.split('-')[1] 

148 strain = json.loads(path.read_text())['strain'] 

149 result[mode].append((structure, strain)) 

150 results[name] = result 

151 Path('acwf-results.json').write_text(json.dumps(results)) 

152 

153 

154def reference_structure(symbol: str, 

155 name: str) -> Atoms: 

156 """Create one of the ACWF structures with the WIEN2K volume.""" 

157 atoms = acwf_structures[name].copy() 

158 x = (volumes[symbol][name] / atoms.get_volume())**(1 / 3) 

159 if symbol == 'O' and name in ['X4O10', 'XO']: 

160 x *= 2**(1 / 3) 

161 atoms.set_cell(atoms.cell * x, scale_atoms=True) 

162 atoms.numbers[atoms.numbers == 0] = atomic_numbers[symbol] 

163 return atoms 

164 

165 

166acwf_structures = { 

167 'XO3': Atoms( 

168 'XO3', 

169 cell=[[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]], 

170 scaled_positions=[[0.0, 0.0, 0.0], [0.5, 0.0, 0.0], 

171 [0.0, 0.5, 0.0], [0.0, 0.0, 0.5]], 

172 pbc=True), 

173 'XO': Atoms( 

174 'XO', 

175 cell=[[0.0, 1.0, 1.0], [1.0, 0.0, 1.0], [1.0, 1.0, 0.0]], 

176 scaled_positions=[[0.0, 0.0, -0.0], [0.5, -0.5, 0.5]], 

177 pbc=True), 

178 'X4O6': Atoms( 

179 'X4O6', 

180 cell=[[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]], 

181 scaled_positions=[[0.25, 0.25, 0.25], [0.75, 0.75, 0.25], 

182 [0.75, 0.25, 0.75], [0.25, 0.75, 0.75], 

183 [0.5, 0.0, 0.0], [0.0, 0.5, 0.0], 

184 [0.0, 0.0, 0.5], [0.5, 0.5, 0.0], 

185 [0.5, 0.0, 0.5], [0.0, 0.5, 0.5]], 

186 pbc=True), 

187 'XO2': Atoms( 

188 'XO2', 

189 cell=[[0.0, 1.0, 1.0], [1.0, 0.0, 1.0], [1.0, 1.0, 0.0]], 

190 scaled_positions=[[0.0, 0.0, -0.0], 

191 [0.25, 1.25, 0.25], 

192 [-0.25, 0.75, 0.75]], 

193 pbc=True), 

194 'X4O10': Atoms( 

195 'X4O10', 

196 cell=[[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]], 

197 scaled_positions=[[0.75, 0.75, 0.75], [0.25, 0.25, 0.75], 

198 [0.25, 0.75, 0.25], [0.75, 0.25, 0.25], 

199 [0.25, 0.75, 0.75], [0.75, 0.25, 0.75], 

200 [0.75, 0.75, 0.25], [0.25, 0.25, 0.25], 

201 [0.5, 0.5, 0.0], [0.5, 0.0, 0.5], 

202 [0.0, 0.5, 0.5], [0.5, 0.0, 0.0], 

203 [0.0, 0.5, 0.0], [0.0, 0.0, 0.5]], 

204 pbc=True), 

205 'X2O': Atoms( 

206 'X2O', 

207 cell=[[0.0, 1.0, 1.0], [1.0, 0.0, 1.0], [1.0, 1.0, 0.0]], 

208 scaled_positions=[[0.25, 1.25, 0.25], 

209 [-0.25, 0.75, 0.75], 

210 [0.0, 0.0, -0.0]], 

211 pbc=True), 

212 'Diamond': Atoms( 

213 'X2', 

214 cell=[[0.0, 1.0, 1.0], [1.0, 0.0, 1.0], [1.0, 1.0, 0.0]], 

215 scaled_positions=[[0.0, 0.0, -0.0], [0.25, 0.25, 0.25]], 

216 pbc=True), 

217 'FCC': Atoms( 

218 'X', 

219 cell=[[0.0, 1.0, 1.0], [1.0, 0.0, 1.0], [1.0, 1.0, 0.0]], 

220 scaled_positions=[[0.0, 0.0, -0.0]], 

221 pbc=True), 

222 'SC': Atoms( 

223 'X', 

224 cell=[[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]], 

225 scaled_positions=[[0.0, 0.0, 0.0]], 

226 pbc=True), 

227 'BCC': Atoms( 

228 'X', 

229 cell=[[-1.0, 1.0, 1.0], [1.0, -1.0, 1.0], [1.0, 1.0, -1.0]], 

230 scaled_positions=[[-0.0, 0.0, 0.0]], 

231 pbc=True)} 

232 

233 

234volumes: dict[str, dict[str, float]] = { 

235 'Ac': { 

236 'BCC': 45.943689065180855, 

237 'Diamond': 129.76193763607355, 

238 'FCC': 45.55065392302246, 

239 'SC': 49.82906214069011, 

240 'X2O': 80.3766777041316, 

241 'X4O6': 203.326117907875, 

242 'X4O10': 207.76797985303918, 

243 'XO': 38.992350681601565, 

244 'XO2': 52.96179138555256, 

245 'XO3': 91.44195279647013 

246 }, 

247 'Ag': { 

248 'BCC': 17.98159943412582, 

249 'Diamond': 60.138447753538145, 

250 'FCC': 17.83856038211966, 

251 'SC': 20.807816653747608, 

252 'X2O': 48.38808819777706, 

253 'X4O6': 131.2721383029882, 

254 'X4O10': 149.92799548230082, 

255 'XO': 25.50923480620743, 

256 'XO2': 34.38335555693807, 

257 'XO3': 65.23416320682753 

258 }, 

259 'Al': { 

260 'BCC': 16.925656872192757, 

261 'Diamond': 55.25812791788162, 

262 'FCC': 16.49535905981626, 

263 'SC': 20.155139416208016, 

264 'X2O': 46.21663447693852, 

265 'X4O6': 98.4792078105453, 

266 'X4O10': 126.14348555770839, 

267 'XO': 22.456524937543783, 

268 'XO2': 26.300916999476847, 

269 'XO3': 49.74390043966775 

270 }, 

271 'Am': { 

272 'BCC': 16.191064699425894, 

273 'Diamond': 38.92606088376118, 

274 'FCC': 17.363741775723295, 

275 'SC': 16.11661629141215, 

276 'X2O': 45.419636119060456, 

277 'X4O6': 144.59436691493505, 

278 'X4O10': 148.57639001361287, 

279 'XO': 26.66553580766982, 

280 'XO2': 36.96323694007846, 

281 'XO3': 67.85355676723577 

282 }, 

283 'Ar': { 

284 'BCC': 53.35453138568846, 

285 'Diamond': 197.20372578528713, 

286 'FCC': 52.276383334536746, 

287 'SC': 65.24068819085674, 

288 'X2O': 113.25843565850016, 

289 'X4O6': 206.36298900234982, 

290 'X4O10': 185.38977453731096, 

291 'XO': 39.43160913365185, 

292 'XO2': 51.33857264225451, 

293 'XO3': 68.32594446494795 

294 }, 

295 'As': { 

296 'BCC': 19.052151090336626, 

297 'Diamond': 57.03345638956617, 

298 'FCC': 19.317577518973586, 

299 'SC': 20.367749732901732, 

300 'X2O': 47.81293031010519, 

301 'X4O6': 130.57910836008148, 

302 'X4O10': 144.03087762987815, 

303 'XO': 27.191332237135178, 

304 'XO2': 31.49881160785104, 

305 'XO3': 47.319197509869724 

306 }, 

307 'At': { 

308 'BCC': 40.007256561713156, 

309 'Diamond': 133.82441316731246, 

310 'FCC': 39.030666540798215, 

311 'SC': 46.14406680206717, 

312 'X2O': 77.79112798187565, 

313 'X4O6': 177.486626152432, 

314 'X4O10': 169.0902164530073, 

315 'XO': 40.774272420246774, 

316 'XO2': 43.83458777704999, 

317 'XO3': 71.71609632447499 

318 }, 

319 'Au': { 

320 'BCC': 18.042083697768916, 

321 'Diamond': 58.541509746305394, 

322 'FCC': 17.97894940070536, 

323 'SC': 20.769602790294257, 

324 'X2O': 49.45298956153372, 

325 'X4O6': 134.26817861354488, 

326 'X4O10': 150.36198987822257, 

327 'XO': 26.881873150878082, 

328 'XO2': 35.00356102786805, 

329 'XO3': 63.961952949423605 

330 }, 

331 'B': { 

332 'BCC': 6.139355211328122, 

333 'Diamond': 16.626316908817323, 

334 'FCC': 5.891550054407846, 

335 'SC': 6.7004880836453475, 

336 'X2O': 30.02852147513494, 

337 'X4O6': 71.97017487676789, 

338 'X4O10': 108.7466438012355, 

339 'XO': 14.710476771224293, 

340 'XO2': 20.373538745231247, 

341 'XO3': 33.628829190051064 

342 }, 

343 'Ba': { 

344 'BCC': 63.30541055553059, 

345 'Diamond': 113.16998237960502, 

346 'FCC': 64.11404100611098, 

347 'SC': 61.60561802807037, 

348 'X2O': 79.89102772639173, 

349 'X4O6': 233.84319295075028, 

350 'X4O10': 227.20804659070285, 

351 'XO': 43.30590522313118, 

352 'XO2': 60.342062817349294, 

353 'XO3': 91.21675807848992 

354 }, 

355 'Be': { 

356 'BCC': 7.815756001821475, 

357 'Diamond': 29.37398729097643, 

358 'FCC': 7.8716206738707015, 

359 'SC': 10.26723560161619, 

360 'X2O': 26.833686313645984, 

361 'X4O6': 77.59929916270295, 

362 'X4O10': 114.64908157675103, 

363 'XO': 12.124123532413858, 

364 'XO2': 22.18537049044565, 

365 'XO3': 38.99206686558694 

366 }, 

367 'Bi': { 

368 'BCC': 31.634667577446024, 

369 'Diamond': 96.90223944114533, 

370 'FCC': 31.810468127145906, 

371 'SC': 35.20505822838321, 

372 'X2O': 69.24164021274835, 

373 'X4O6': 174.87946592197846, 

374 'X4O10': 177.94746787410855, 

375 'XO': 36.02504757390335, 

376 'XO2': 42.07312443213904, 

377 'XO3': 74.04342790643858 

378 }, 

379 'Br': { 

380 'BCC': 26.784212478063857, 

381 'Diamond': 86.17112439244427, 

382 'FCC': 26.417536335741566, 

383 'SC': 29.8278511781182, 

384 'X2O': 61.628473750026465, 

385 'X4O6': 147.33738792781116, 

386 'X4O10': 155.97272615705282, 

387 'XO': 31.106626065259217, 

388 'XO2': 37.01371722315619, 

389 'XO3': 58.256274759446185 

390 }, 

391 'C': { 

392 'BCC': 6.685722025660118, 

393 'Diamond': 11.391522946768983, 

394 'FCC': 7.321634771212951, 

395 'SC': 5.582211299597029, 

396 'X2O': 28.554339698563062, 

397 'X4O6': 84.53807260317583, 

398 'X4O10': 117.2300823463996, 

399 'XO': 15.70175180655379, 

400 'XO2': 22.756394307135132, 

401 'XO3': 32.919218490588555 

402 }, 

403 'Ca': { 

404 'BCC': 42.15047346676994, 

405 'Diamond': 159.96939544107403, 

406 'FCC': 42.19429685750488, 

407 'SC': 43.584913194130294, 

408 'X2O': 56.093056556486744, 

409 'X4O6': 159.87070574443365, 

410 'X4O10': 168.84121978259935, 

411 'XO': 28.190280638134226, 

412 'XO2': 42.30717797818566, 

413 'XO3': 94.28216713499612 

414 }, 

415 'Cd': { 

416 'BCC': 23.419572636611186, 

417 'Diamond': 74.9080478398366, 

418 'FCC': 22.84127817506259, 

419 'SC': 26.923600245740374, 

420 'X2O': 53.08274544890885, 

421 'X4O6': 146.50420325772748, 

422 'X4O10': 161.23783826032457, 

423 'XO': 27.093952692307997, 

424 'XO2': 38.53660786052451, 

425 'XO3': 73.15212767492528 

426 }, 

427 'Ce': { 

428 'BCC': 27.323964391559745, 

429 'Diamond': 60.36962802029463, 

430 'FCC': 26.52235577369828, 

431 'SC': 24.921013988065226, 

432 'X2O': 56.10337207266545, 

433 'X4O6': 164.25418463410853, 

434 'X4O10': 168.09408301376885, 

435 'XO': 31.03931067267397, 

436 'XO2': 40.79349906605082, 

437 'XO3': 81.65266123145898 

438 }, 

439 'Cl': { 

440 'BCC': 21.45474447156215, 

441 'Diamond': 67.51485775600094, 

442 'FCC': 21.28818952390059, 

443 'SC': 23.464029289304705, 

444 'X2O': 54.559995010249025, 

445 'X4O6': 135.2505227617596, 

446 'X4O10': 150.35467864680692, 

447 'XO': 26.958660419342714, 

448 'XO2': 34.69762129543773, 

449 'XO3': 52.00560808001245 

450 }, 

451 'Cm': { 

452 'BCC': 16.520782253696993, 

453 'Diamond': 38.65236260941347, 

454 'FCC': 17.49246649598698, 

455 'SC': 16.399002476823348, 

456 'X2O': 45.8647310715384, 

457 'X4O6': 144.5823454157204, 

458 'X4O10': 149.2067378771397, 

459 'XO': 26.794240945658323, 

460 'XO2': 36.65760496498912, 

461 'XO3': 67.29565446076828 

462 }, 

463 'Co': { 

464 'BCC': 10.544792682061145, 

465 'Diamond': 29.76943354558406, 

466 'FCC': 10.30839278272375, 

467 'SC': 11.89326679933492, 

468 'X2O': 29.835635132747164, 

469 'X4O6': 90.62563815518891, 

470 'X4O10': 117.91158291691508, 

471 'XO': 17.2707951704295, 

472 'XO2': 25.033682445053497, 

473 'XO3': 47.13776514479808 

474 }, 

475 'Cr': { 

476 'BCC': 11.548199339585498, 

477 'Diamond': 33.07091175174273, 

478 'FCC': 11.88591646485507, 

479 'SC': 12.806958012708076, 

480 'X2O': 30.441620462552265, 

481 'X4O6': 96.55853606378332, 

482 'X4O10': 116.26560686420035, 

483 'XO': 17.664208258348324, 

484 'XO2': 25.425018327785768, 

485 'XO3': 46.87447267496256 

486 }, 

487 'Cs': { 

488 'BCC': 116.84172285179223, 

489 'Diamond': 377.51195120825054, 

490 'FCC': 117.36057894318543, 

491 'SC': 128.35905522569587, 

492 'X2O': 96.49631431852791, 

493 'X4O6': 283.36831128941446, 

494 'X4O10': 244.2227704566822, 

495 'XO': 53.296471975757754, 

496 'XO2': 63.55584224078437, 

497 'XO3': 76.20321967007675 

498 }, 

499 'Cu': { 

500 'BCC': 12.004516294003697, 

501 'Diamond': 38.35431876920354, 

502 'FCC': 11.952224652322, 

503 'SC': 13.93497312424489, 

504 'X2O': 34.538824246939555, 

505 'X4O6': 103.80050324183955, 

506 'X4O10': 128.49045023552725, 

507 'XO': 19.107179475952414, 

508 'XO2': 28.233135239309007, 

509 'XO3': 52.01329300394603 

510 }, 

511 'Dy': { 

512 'BCC': 32.288569204367704, 

513 'Diamond': 46.02037644692374, 

514 'FCC': 32.47663596342035, 

515 'SC': 31.857193083077554, 

516 'X2O': 51.791906354845, 

517 'X4O6': 148.57156262807263, 

518 'X4O10': 155.3517188345445, 

519 'XO': 28.17588955662974, 

520 'XO2': 36.94493400809891, 

521 'XO3': 71.77986544273953 

522 }, 

523 'Er': { 

524 'BCC': 33.930706759113306, 

525 'Diamond': 160.71821137423967, 

526 'FCC': 34.82288485206389, 

527 'SC': 35.94303590731403, 

528 'X2O': 53.81459081043589, 

529 'X4O6': 147.37472200530317, 

530 'X4O10': 156.02502208600697, 

531 'XO': 28.024032432884447, 

532 'XO2': 36.90354261761303, 

533 'XO3': 73.1595986199416 

534 }, 

535 'Eu': { 

536 'BCC': 26.134083220040758, 

537 'Diamond': 41.36507024860562, 

538 'FCC': 24.991976993289317, 

539 'SC': 17.795148296296187, 

540 'X2O': 50.001667417550664, 

541 'X4O6': 152.32836518004666, 

542 'X4O10': 156.76478360639607, 

543 'XO': 28.504144398060703, 

544 'XO2': 37.69178097711213, 

545 'XO3': 71.09394505182284 

546 }, 

547 'F': { 

548 'BCC': 10.084127330253507, 

549 'Diamond': 29.00480135135658, 

550 'FCC': 10.146987257565563, 

551 'SC': 10.520941986117641, 

552 'X2O': 30.965700436027305, 

553 'X4O6': 112.64541976601427, 

554 'X4O10': 121.58537901423888, 

555 'XO': 18.85666827649016, 

556 'XO2': 30.8381241671117, 

557 'XO3': 53.32461728600295 

558 }, 

559 'Fe': { 

560 'BCC': 10.500484196384424, 

561 'Diamond': 28.92444480806717, 

562 'FCC': 10.260211776246214, 

563 'SC': 11.651347771934487, 

564 'X2O': 29.393596433158027, 

565 'X4O6': 90.8708646557079, 

566 'X4O10': 116.01663708535033, 

567 'XO': 17.13602302251458, 

568 'XO2': 24.193568321607785, 

569 'XO3': 45.93568036284212 

570 }, 

571 'Fr': { 

572 'BCC': 116.49231248676648, 

573 'Diamond': 384.032159640575, 

574 'FCC': 117.16295360986017, 

575 'SC': 132.17408585059354, 

576 'X2O': 106.7081347574964, 

577 'X4O6': 262.19120504163266, 

578 'X4O10': 229.91803402594206, 

579 'XO': 55.39235301269504, 

580 'XO2': 58.91457051825229, 

581 'XO3': 79.51922007961696 

582 }, 

583 'Ga': { 

584 'BCC': 19.205587610065592, 

585 'Diamond': 50.83804519831408, 

586 'FCC': 18.94650008463183, 

587 'SC': 20.117312153714767, 

588 'X2O': 53.04283169078497, 

589 'X4O6': 110.9910906198663, 

590 'X4O10': 135.6777983405775, 

591 'XO': 24.3519260047855, 

592 'XO2': 29.261674771756013, 

593 'XO3': 57.47020127815971 

594 }, 

595 'Gd': { 

596 'BCC': 28.947403790341237, 

597 'Diamond': 41.89427015507523, 

598 'FCC': 27.993870497573496, 

599 'SC': 20.809598933525876, 

600 'X2O': 50.44857514811846, 

601 'X4O6': 151.01960637372918, 

602 'X4O10': 156.1541161775155, 

603 'XO': 28.361967230056884, 

604 'XO2': 37.36132396877068, 

605 'XO3': 71.11221246370377 

606 }, 

607 'Ge': { 

608 'BCC': 19.269495309541668, 

609 'Diamond': 47.82658983637545, 

610 'FCC': 19.582450801516337, 

611 'SC': 19.941748311438904, 

612 'X2O': 49.682664655338414, 

613 'X4O6': 121.20682589314904, 

614 'X4O10': 139.51685630367012, 

615 'XO': 27.122474603810435, 

616 'XO2': 28.190777443226906, 

617 'XO3': 51.07396911183831 

618 }, 

619 'H': { 

620 'BCC': 2.9667742929491983, 

621 'Diamond': 6.831111656170396, 

622 'FCC': 2.964809112631019, 

623 'SC': 3.0866686938954198, 

624 'X2O': 11.96347891754409, 

625 'X4O6': 61.84619062082082, 

626 'X4O10': 103.34760585605244, 

627 'XO': 10.027848555766598, 

628 'XO2': 19.10459642792833, 

629 'XO3': 31.08928579616201 

630 }, 

631 'He': { 

632 'BCC': 18.03038805667906, 

633 'Diamond': 64.21487893655794, 

634 'FCC': 17.77258072097029, 

635 'SC': 21.485271132346007, 

636 'X2O': 92.1784561523737, 

637 'X4O6': 94.28925599148977, 

638 'X4O10': 112.79696925975088, 

639 'XO': 31.528419933938714, 

640 'XO2': 24.776988443791783, 

641 'XO3': 43.05953909874314 

642 }, 

643 'Hf': { 

644 'BCC': 22.304719470255996, 

645 'Diamond': 70.15774064881346, 

646 'FCC': 22.56760341919224, 

647 'SC': 24.781139308615224, 

648 'X2O': 48.646400121650814, 

649 'X4O6': 128.35931068254416, 

650 'X4O10': 142.6630326850143, 

651 'XO': 24.23287799836855, 

652 'XO2': 33.11711229980498, 

653 'XO3': 70.82641284222711 

654 }, 

655 'Hg': { 

656 'BCC': 29.237160212480383, 

657 'Diamond': 112.58234465071779, 

658 'FCC': 32.347797693222105, 

659 'SC': 29.854323349822568, 

660 'X2O': 56.18140344248944, 

661 'X4O6': 152.67994259762438, 

662 'X4O10': 166.5429926971196, 

663 'XO': 29.815337246318833, 

664 'XO2': 39.30320403679202, 

665 'XO3': 72.63086986279107 

666 }, 

667 'Ho': { 

668 'BCC': 33.267285286514536, 

669 'Diamond': 50.850112182697885, 

670 'FCC': 33.89173938695238, 

671 'SC': 34.29406972261268, 

672 'X2O': 52.720584538970684, 

673 'X4O6': 147.77237454797466, 

674 'X4O10': 155.498638197736, 

675 'XO': 28.09751593211625, 

676 'XO2': 36.88059139506703, 

677 'XO3': 72.36955554798217 

678 }, 

679 'I': { 

680 'BCC': 35.98676703422804, 

681 'Diamond': 121.14527750115774, 

682 'FCC': 35.10488450303778, 

683 'SC': 41.563154436885426, 

684 'X2O': 72.18837519600658, 

685 'X4O6': 167.70157634361635, 

686 'X4O10': 163.00329788775545, 

687 'XO': 38.012187732710515, 

688 'XO2': 41.44308732251686, 

689 'XO3': 60.8593853544075 

690 }, 

691 'In': { 

692 'BCC': 27.78058028977925, 

693 'Diamond': 76.43205868239434, 

694 'FCC': 27.510098861849173, 

695 'SC': 29.550210982224606, 

696 'X2O': 65.93128247399382, 

697 'X4O6': 140.37972957567774, 

698 'X4O10': 156.84754366426756, 

699 'XO': 30.493020046452415, 

700 'XO2': 36.4728530560971, 

701 'XO3': 74.76771621358446 

702 }, 

703 'Ir': { 

704 'BCC': 15.055641509624866, 

705 'Diamond': 43.19355479899218, 

706 'FCC': 14.504994129130639, 

707 'SC': 16.994504125715814, 

708 'X2O': 40.335773558548205, 

709 'X4O6': 111.03802506574334, 

710 'X4O10': 130.89009942137127, 

711 'XO': 23.407091787120702, 

712 'XO2': 30.412061386006737, 

713 'XO3': 53.86205787750703 

714 }, 

715 'K': { 

716 'BCC': 73.77951798925906, 

717 'Diamond': 223.8044852857056, 

718 'FCC': 74.0044361285296, 

719 'SC': 79.35441292282184, 

720 'X2O': 68.00231952386676, 

721 'X4O6': 227.58585562272344, 

722 'X4O10': 198.7688451835494, 

723 'XO': 42.74542841231468, 

724 'XO2': 58.364291316549775, 

725 'XO3': 136.53458069146527 

726 }, 

727 'Kr': { 

728 'BCC': 67.46343111231373, 

729 'Diamond': 249.70798358307468, 

730 'FCC': 66.04176534194272, 

731 'SC': 82.66454249647283, 

732 'X2O': 118.9585508090093, 

733 'X4O6': 194.32198413188252, 

734 'X4O10': 192.3183698028015, 

735 'XO': 40.074708805542606, 

736 'XO2': 47.54241675063361, 

737 'XO3': 68.05698604893871 

738 }, 

739 'La': { 

740 'BCC': 37.81757025987239, 

741 'Diamond': 74.6380534647987, 

742 'FCC': 36.94689152187624, 

743 'SC': 36.74506294698075, 

744 'X2O': 65.1355386334241, 

745 'X4O6': 182.51664661894716, 

746 'X4O10': 189.9565090452286, 

747 'XO': 34.431148430436664, 

748 'XO2': 47.879356819541435, 

749 'XO3': 90.87226882578761 

750 }, 

751 'Li': { 

752 'BCC': 20.26746521326914, 

753 'Diamond': 51.386718035200765, 

754 'FCC': 20.224478592944862, 

755 'SC': 20.408801974058356, 

756 'X2O': 24.721960544273422, 

757 'X4O6': 89.47029219795057, 

758 'X4O10': 123.69216300739882, 

759 'XO': 16.825820190091974, 

760 'XO2': 24.922284658946587, 

761 'XO3': 46.08915629110866 

762 }, 

763 'Lu': { 

764 'BCC': 29.625706403168756, 

765 'Diamond': 101.21824904812954, 

766 'FCC': 28.971404244986218, 

767 'SC': 32.94002697403406, 

768 'X2O': 54.07084024181938, 

769 'X4O6': 144.42641168708099, 

770 'X4O10': 157.87453325697595, 

771 'XO': 26.56021013564612, 

772 'XO2': 37.2635986858799, 

773 'XO3': 78.08376717694293 

774 }, 

775 'Mg': { 

776 'BCC': 22.917256964511758, 

777 'Diamond': 80.85167299497583, 

778 'FCC': 23.125246144035756, 

779 'SC': 27.580585068743265, 

780 'X2O': 44.86183000035521, 

781 'X4O6': 113.54760895185626, 

782 'X4O10': 138.89600330910122, 

783 'XO': 19.249038606577802, 

784 'XO2': 30.59038256739605, 

785 'XO3': 62.138771825273366 

786 }, 

787 'Mn': { 

788 'BCC': 10.780826839746656, 

789 'Diamond': 30.343273943677755, 

790 'FCC': 10.74714431323143, 

791 'SC': 11.898840334365449, 

792 'X2O': 29.56938020768425, 

793 'X4O6': 92.93093174774003, 

794 'X4O10': 115.32569138735641, 

795 'XO': 17.29322893846185, 

796 'XO2': 24.65355351858858, 

797 'XO3': 45.60581266709946 

798 }, 

799 'Mo': { 

800 'BCC': 15.79258104697389, 

801 'Diamond': 46.0075591929068, 

802 'FCC': 16.03513025909276, 

803 'SC': 17.596261003631483, 

804 'X2O': 39.71697359069131, 

805 'X4O6': 113.5385735057236, 

806 'X4O10': 125.80601908145752, 

807 'XO': 21.517172018665377, 

808 'XO2': 29.728405877601176, 

809 'XO3': 55.619143262267016 

810 }, 

811 'N': { 

812 'BCC': 7.234713384391083, 

813 'Diamond': 18.353307328655553, 

814 'FCC': 7.601290987311529, 

815 'SC': 6.4799301718134945, 

816 'X2O': 26.668195825122606, 

817 'X4O6': 91.34940280257845, 

818 'X4O10': 114.46774472474029, 

819 'XO': 15.338064037074709,