Attachment 'CO_in_a_box.py'
Download 1 """CO dimer in a box."""
2
3 import os
4 from Dacapo import Dacapo
5 from ASE import Atom, ListOfAtoms
6
7 atoms = ListOfAtoms([Atom('C', (2, 2, 2),magmom=1.),
8 Atom('O', (3.6, 2, 2),magmom=1,)],
9 cell=(5,3,3), periodic=1)
10
11 calc = Dacapo(planewavecutoff=350, # in eV
12 nbands=6, # 5 extra empty bands
13 spinpol=False,
14 out='CO_in_a_box.nc',
15 txtout='CO_in_a_box.txt')
16
17 atoms.SetCalculator(calc)
18 calc.SetNetCDFEntry('ConvergenceControl',att='MaxNumberOfSteps',value=3)
19 energy = atoms.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.