1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
# fmt: off
"""Module to read atoms in chemical json file format.
https://wiki.openchemistry.org/Chemical_JSON
"""
import json
import numpy as np
from ase import Atoms
from ase.cell import Cell
# contract and lower case string
def contract(dictionary):
dcopy = {key.replace(' ', '').lower(): dictionary[key] for key in
dictionary}
return dcopy
def read_cjson(fileobj):
"""Read a Chemical Json file as written by avogadro2 (>=1.93.0)
See https://wiki.openchemistry.org/Chemical_JSON
"""
data = contract(json.load(fileobj))
atoms = Atoms()
datoms = data['atoms']
atoms = Atoms(datoms['elements']['number'])
if 'unitcell' in data:
cell = data['unitcell']
a = cell['a']
b = cell['b']
c = cell['c']
alpha = cell['alpha']
beta = cell['beta']
gamma = cell['gamma']
atoms.cell = Cell.fromcellpar([a, b, c, alpha, beta, gamma])
atoms.pbc = True
coords = contract(datoms['coords'])
if '3d' in coords:
positions = np.array(coords['3d']).reshape(len(atoms), 3)
atoms.set_positions(positions)
else:
positions = np.array(coords['3dfractional']).reshape(len(atoms), 3)
atoms.set_scaled_positions(positions)
yield atoms
|
