This is an old revision of the document!
いっとき流行ったホウ素ドープダイアモンドの仮想結晶近似を例にVCAを説明します。
プロジェクト名は Diamond とします
仮想結晶近似を行う原子の電子配置を調べておきます。
structファイルは次のものを使用しました。
blebleble F LATTICE,NONEQUIV.ATOMS: 1 227 Fd-3m MODE OF CALC=RELA unit=bohr 6.740880 6.740880 6.740880 90.000000 90.000000 90.000000 ATOM 1: X=0.87500000 Y=0.37500000 Z=0.37500000 MULT= 2 ISPLIT=15 1: X=0.12500000 Y=0.12500000 Z=0.12500000 C 1 NPT= 781 R0=.000100000 RMT= 1.45 Z: 6.00000 LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0 NUMBER OF SYMMETRY OPERATIONS
$ init_lapw -b -numk 1000 -rkmax 5.5 $ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -p
blebleble F LATTICE,NONEQUIV.ATOMS: 1 227 Fd-3m MODE OF CALC=RELA unit=bohr 6.740880 6.740880 6.740880 90.000000 90.000000 90.000000 ATOM 1: X=0.87500000 Y=0.37500000 Z=0.37500000 MULT= 2 ISPLIT= 2 1: X=0.12500000 Y=0.12500000 Z=0.12500000 C 1 NPT= 781 R0=.000100000 RMT= 1.45 Z: 5.95 LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 0.0000000 1.0000000 48 NUMBER OF SYMMETRY OPERATIONS ...
TOT (TOT,FOR,QTL,EFG,FERMI) -12.0 7.9 0.50 0.05 1 EMIN, NE, ESEPERMIN, ESEPER0, iqtlsave TETRA 0.000 (GAUSS,ROOT,TEMP,TETRA,ALL eval) 0 0 4 0 4 4 6 0 6 4 -3 2 12.00 GMAX NOFILE FILE/NOFILE write recprlist
C He 3 2,-1,1.0 N 2,-1,1.0 N 2, 1,1.0 N 2, 1,0.0 N 2,-2,0.95 N 2,-2,0.0 N **** **** END of input (instgen_lapw)
$ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -p -NI
$ grep ":GAP" Diamond.scf