Attachment 'Fe-ferro_compl.py'
Download 1 from Dacapo import Dacapo
2 from ASE import Atom, ListOfAtoms
3
4
5 a = 2.87
6 bulk = ListOfAtoms([Atom('Fe', (0, 0, 0), magmom=2.2),
7 Atom('Fe', (a/2, a/2, a/2), magmom=2.2)],
8 cell=(a, a, a))
9
10 calc = Dacapo(spinpol=True,
11 kpts=(8, 8, 8),
12 planewavecutoff=300,
13 densitycutoff=500,
14 nbands=17,
15 usesymm=True,
16 out='ferro.nc',
17 txtout='ferro.txt')
18
19 bulk.SetCalculator(calc)
20
21 print bulk.GetPotentialEnergy()
Attached Files
To refer to attachments on a page, use attachment:filename, as shown below in the list of files. Do NOT use the URL of the [get] link, since this is subject to change and can break easily.You are not allowed to attach a file to this page.