#!/usr/bin/env python
# -*- coding: iso-8859-15 -*-
# この行（↑）を入れることにより日本語(UTF-8)でコメントを入れることができるようになる
 
# 単位格子ベクトル（長さの単位はÅ）
a1 = (5.0, 0.0, 0.0)
a2 = (0.0, 5.0, 0.0)
a3 = (0.0, 0.0, 5.0)
 
# 単位格子
box1 = [a1, a2, a3]
 
# 原子の位置（単位格子に対する相対座標）
r1 = (0.50, 0.50, 0.50)  # 炭素原子（単位格子の中央）
r2 = (0.72, 0.50, 0.50)  # 酸素原子（0.72 = 0.5 + 1.1Å/5.0Å）
 
from ASE import Atom
atom1 = Atom('C', r1)
atom2 = Atom('O', r2)
molecule1 = []
molecule1.append(atom1)  # molecule1 = [atom1, atom2]でもよいが、
molecule1.append(atom2)  # 原子数が増えるとこの方が便利
 
from ASE import ListOfAtoms
model1 = ListOfAtoms(molecule1)
model1.SetUnitCell(box1)
 
from Dacapo import Dacapo 
method1 = Dacapo()
method1.SetNumberOfBands(8) # 計算する電子状態の個数 = 総価電子数/2 + アルファ
 
model1.SetCalculator(method1)
method1.Calculate()

$ chmod +x lesson02.py 
$ ./lesson02.py 
Running /usr/bin//dacapo_2.7.7.run ./tmpXG27Yj out.nc -out out.txt -scratch . -stop stop
$