Lesson 3 : スクリプトとして実行する

#!/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
$
ab_initio/dacapo_tutorial/lesson3.txt · 最終更新: 2009/08/26 14:23 by kimi
www.chimeric.de Creative Commons License Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0