グラフェンを例に、より詳しい計算の演習を行おう。

a0 = 2.461 # Lattice constant of Graphene
 
from math import sqrt
c1 = ( a0*sqrt(3.0)/2.0, a0/2.0, 0.0)
c2 = (-a0*sqrt(3.0)/2.0, a0/2.0, 0.0)
c3 = (0.0, 0.0, a0*8.0)              # Height of 3D unit cell
 
from ase import Atom
atom1 = Atom('C', (1.0/3.0, 2.0/3.0, 0.5))
atom2 = Atom('C', (2.0/3.0, 1.0/3.0, 0.5))
 
from ase import Atoms
box1 = Atoms([atom1, atom2], pbc = True)
box1.set_cell([c1, c2, c3], scale_atoms = True)
 
from ase.calculators.jacapo import Jacapo
solver1 = Jacapo(nbands = 8, kpts = (2, 2, 1), pw = 150, dw = 150)
box1.set_calculator(solver1)
solver1.calculate()

a0 = 2.461 # Lattice constant of Graphene
 
from math import sqrt
c1 = ( a0*sqrt(3.0)/2.0, a0/2.0, 0.0)
c2 = (-a0*sqrt(3.0)/2.0, a0/2.0, 0.0)
c3 = (0.0, 0.0, a0*8.0)              # Height of 3D unit cell
 
unit_cell = [c1, c2, c3] # Unit cell

from ase import Atom
atom1 = Atom('C', (1.0/3.0, 2.0/3.0, 0.5))
atom2 = Atom('C', (2.0/3.0, 1.0/3.0, 0.5))
 
list_of_atoms = [atom1, atom2]

# -------------------------
# Calculational conditions
# -------------------------
kwargs = {}
kwargs['nbands'] = 8        # Number of Valence electrons
kwargs['kpts'] = (2, 2, 1)  # k-point sampling in the 1st BZ
kwargs['pw'] = 150          # Energy cutoff for wave functions
kwargs['dw'] = 150          # Energy cutoff for charge density

# ---------------------------
# Common part of calculations
# ---------------------------
from ase import Atoms
box1 = Atoms(list_of_atoms, pbc = True)
box1.set_cell(unit_cell, scale_atoms = True)
from ase.calculators.jacapo import Jacapo
solver1 = Jacapo(**kwargs)
box1.set_calculator(solver1)
solver1.calculate()

# ------------------
# Output on screen
# ------------------ 
print box1.get_total_energy()