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
|
from __future__ import print_function
import os
from ase.build import molecule
from ase.io.bader import attach_charges
fname = 'ACF.dat'
f = open(fname, 'w')
print("""
# X Y Z CHARGE MIN DIST
----------------------------------------------------------------
1 7.0865 8.5038 9.0672 9.0852 1.3250
2 7.0865 9.9461 7.9403 0.4574 0.3159
3 7.0865 7.0615 7.9403 0.4574 0.3159
----------------------------------------------------------------
NUMBER OF ELECTRONS: 9.99999
""", file=f)
f.close()
atoms = molecule('H2O')
atoms.set_cell([7.5, 9, 9])
atoms.center()
attach_charges(atoms)
attach_charges(atoms, fname)
os.remove(fname)
for atom in atoms:
print('Atom', atom.symbol, 'Bader charge', atom.charge)
|
