第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* 銅(111)表面 [#w9763304]
ここではPBEで最適化したバルクの格子定数を用いて作成した6...
スラブの構造(原子位置)は例えば[[slabgen>計算コードの使...
** SCF計算 [#k70d7070]
- 入力ファイル (nfinp_scf)
0 0 0 0 0 0
6.00000 20.00000 1 6 6 : GMAX GMAXP NTYP...
1 0 : num_space_group...
Cartesian
4.865948314378 0.000000000000 0.00000000...
-2.432974157189 4.214034853753 0.00000000...
0.000000000000 0.000000000000 79.46060323...
12 12 1 1 1 1 : K_mesh
1 0 : NCORD NINV : IW...
0.000000000000 0.000000000000 9.93257540...
0.000000000000 2.809356569169 5.95954524...
0.000000000000 -2.809356569169 1.98651508...
0.000000000000 0.000000000000 -1.98651508...
0.000000000000 2.809356569169 -5.95954524...
0.000000000000 -2.809356569169 -9.93257540...
29 0.5000 22104.64 6 1 0.2 : IATOMN ALFA AMI...
0 0 0 0 0 : ICOND INIPOS IN...
0 1 : ipre ipri
200 200 0 3200.00 0 : nmd1 nmd2 last_...
6 1 : WAY_MIX MIX_WHAT
0 20 0.30 : iter_start KBX...
0.20 0.30 0.20 0.20 0.20 : dtim1 dtim2 dti...
200.00 4 1 5.00D-12 : dtio imdalg i...
-0.0010 1.00D+03 0 : width forccr is...
ggapbe 1 : xctype kspin
2.00 : destm
101 : nbztyp
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : keg
1 : nextst
0 : 0 random number...
2 : imsd
0 : evaluate_eko_diff
0 : npdosao
0 0.000 : sm_n dopping
STATEの実行
mpirun -np 8 ./STATE < nfinp_scf > nfout_scf
** バンド計算 [#hcb8fa4b]
- 入力ファイル (nfinp_band)
0 0 0 0 0 0
6.00000 20.00000 1 6 6 : GMAX, GMAXP, NT...
1 0 : num_space_group...
Cartesian
4.865948314378 0.000000000000 0.00000000...
-2.432974157189 4.214034853753 0.00000000...
0.000000000000 0.000000000000 79.46060323...
12 12 1 1 1 1 : K_mesh
1 0 : NCORD NINV : IW...
0.000000000000 0.000000000000 9.93257540...
0.000000000000 2.809356569169 5.95954524...
0.000000000000 -2.809356569169 1.98651508...
0.000000000000 0.000000000000 -1.98651508...
0.000000000000 2.809356569169 -5.95954524...
0.000000000000 -2.809356569169 -9.93257540...
29 0.5000 22104.64 6 1 0.2 : IATOMN,ALFA,AMI...
22 0 0 0 0 : ICOND 0-MD, 1-C...
0 1 : ipre, ipri
200 200 0 3200.00 0 : nmd1,nmd2,last_...
6 1 : SIMPLE=1,BROYD1...
0 20 0.60 : iter_start, KBX...
0.20 0.30 0.20 0.20 0.20 : dtim1,dtim2,dti...
200.00 4 1 5.0D-12 : dtio ,imdalg, i...
-0.0010 1.00D+03 0 : width,forccr,is...
ggapbe 1 : xctype,kspin
2.00 : destm
101 : nbztyp 0-sf, 1-...
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : keg
1 : nextst
0 : 0 random number...
2 : imsd
0 : evaluate_eko_diff
0 : npdosao
0 0.000 : sm_n, dopping
&KPOINTS_BAND
NKSEG 3
KMESH 40 20 40
KPOINTS
0.0 0.0 0.0
0.5 0.0 0.0
0.33333333 0.33333333 0.0
0.0 0.0 0.0
&END
ここでバンドはGamma (0,0,0) -> M (1/2,0) -> K (1/3,1/3) -...
STATEの実行
mpirun -np 8 ./STATE < nfinp_band > nfout_band
energy.dataが得られるので
energy2band
を実行してバンドをXYデータ形式にしてプロットすると以下の...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
ここでBrillouin zoneは以下のようになっている。
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
得られたバンドを文献のものと比較してみる。例えばEuceda et...
終了行:
* 銅(111)表面 [#w9763304]
ここではPBEで最適化したバルクの格子定数を用いて作成した6...
スラブの構造(原子位置)は例えば[[slabgen>計算コードの使...
** SCF計算 [#k70d7070]
- 入力ファイル (nfinp_scf)
0 0 0 0 0 0
6.00000 20.00000 1 6 6 : GMAX GMAXP NTYP...
1 0 : num_space_group...
Cartesian
4.865948314378 0.000000000000 0.00000000...
-2.432974157189 4.214034853753 0.00000000...
0.000000000000 0.000000000000 79.46060323...
12 12 1 1 1 1 : K_mesh
1 0 : NCORD NINV : IW...
0.000000000000 0.000000000000 9.93257540...
0.000000000000 2.809356569169 5.95954524...
0.000000000000 -2.809356569169 1.98651508...
0.000000000000 0.000000000000 -1.98651508...
0.000000000000 2.809356569169 -5.95954524...
0.000000000000 -2.809356569169 -9.93257540...
29 0.5000 22104.64 6 1 0.2 : IATOMN ALFA AMI...
0 0 0 0 0 : ICOND INIPOS IN...
0 1 : ipre ipri
200 200 0 3200.00 0 : nmd1 nmd2 last_...
6 1 : WAY_MIX MIX_WHAT
0 20 0.30 : iter_start KBX...
0.20 0.30 0.20 0.20 0.20 : dtim1 dtim2 dti...
200.00 4 1 5.00D-12 : dtio imdalg i...
-0.0010 1.00D+03 0 : width forccr is...
ggapbe 1 : xctype kspin
2.00 : destm
101 : nbztyp
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : keg
1 : nextst
0 : 0 random number...
2 : imsd
0 : evaluate_eko_diff
0 : npdosao
0 0.000 : sm_n dopping
STATEの実行
mpirun -np 8 ./STATE < nfinp_scf > nfout_scf
** バンド計算 [#hcb8fa4b]
- 入力ファイル (nfinp_band)
0 0 0 0 0 0
6.00000 20.00000 1 6 6 : GMAX, GMAXP, NT...
1 0 : num_space_group...
Cartesian
4.865948314378 0.000000000000 0.00000000...
-2.432974157189 4.214034853753 0.00000000...
0.000000000000 0.000000000000 79.46060323...
12 12 1 1 1 1 : K_mesh
1 0 : NCORD NINV : IW...
0.000000000000 0.000000000000 9.93257540...
0.000000000000 2.809356569169 5.95954524...
0.000000000000 -2.809356569169 1.98651508...
0.000000000000 0.000000000000 -1.98651508...
0.000000000000 2.809356569169 -5.95954524...
0.000000000000 -2.809356569169 -9.93257540...
29 0.5000 22104.64 6 1 0.2 : IATOMN,ALFA,AMI...
22 0 0 0 0 : ICOND 0-MD, 1-C...
0 1 : ipre, ipri
200 200 0 3200.00 0 : nmd1,nmd2,last_...
6 1 : SIMPLE=1,BROYD1...
0 20 0.60 : iter_start, KBX...
0.20 0.30 0.20 0.20 0.20 : dtim1,dtim2,dti...
200.00 4 1 5.0D-12 : dtio ,imdalg, i...
-0.0010 1.00D+03 0 : width,forccr,is...
ggapbe 1 : xctype,kspin
2.00 : destm
101 : nbztyp 0-sf, 1-...
4 4 4 : DUMMY
4 4 4 : DUMMY
48 : keg
1 : nextst
0 : 0 random number...
2 : imsd
0 : evaluate_eko_diff
0 : npdosao
0 0.000 : sm_n, dopping
&KPOINTS_BAND
NKSEG 3
KMESH 40 20 40
KPOINTS
0.0 0.0 0.0
0.5 0.0 0.0
0.33333333 0.33333333 0.0
0.0 0.0 0.0
&END
ここでバンドはGamma (0,0,0) -> M (1/2,0) -> K (1/3,1/3) -...
STATEの実行
mpirun -np 8 ./STATE < nfinp_band > nfout_band
energy.dataが得られるので
energy2band
を実行してバンドをXYデータ形式にしてプロットすると以下の...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
ここでBrillouin zoneは以下のようになっている。
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
得られたバンドを文献のものと比較してみる。例えばEuceda et...
ページ名: