- 追加された行はこの色です。
- 削除された行はこの色です。
* 例1: bulk Cu [#b2afc8e1]
** 入力ファイル [#a8b3f61e]
** 計算方法 [#a8b3f61e]
入力ファイル
0 0 0 0 0 0
6.00000 20.0000 1 1 1 : GMAX, GMAXP, NTYP, NATM, NATM2
225 2 : num_space_group, type
6.880000 6.880000 6.880000 90.000 90.000 90.000 : a,b,c,alpha,beta,gamma
16 16 16 2 2 2 : K_mesh
0 0 : NCORD, NINV, : IWEI, IMDTYP, ITYP
0.000000d0 0.000000d0 0.000000d0 1 0 1
29 0.5000 22104.64 6 1 0.d0 : TYPE 1IATOMN,ALFA,AMION,ILOC,IVAN
0 0 0 0 0 : ICOND 0-MD, 1-CONT.MD, 2-WAVE FN,, 3-CHARGE DEN.
0 1 : ipre, ipri
200 1000 0 84200.00 0 : nmd1,nmd2,last_iter,cpumax,ifstop
6 1 : SIMPLE=1,BROYD1=2,BROYD2=3,DFP=4,PULAY=5
0 20 0.60 : iter_start, KBXMIX, MIX_ALPHA
0.20 0.30 0.20 0.20 0.20 : dtim1,dtim2,dtim3,dtim4,dtim
300.00 4 1 0.50D-09 : dtio ,imdalg, iexpl, edelta
-10.0002 0.50D+03 0 : width,forccr,istress
ggapbe 1 : xctype,kspin
2.00 : destm
101 : nbztyp 0-sf, 1-bk, 2-sc, 3-bcc, 4-fcc, 5-dia, 6-hex
4 4 4
4 4 4
12 : keg
1 : nextst(mb)
0 : 0 random numbers 1: matrix
2 : imsd 1=MSD, 2=Preconditioned
0 : evaluate_eko_diff
0 : npdosao
0 0.000 : sm_n, dopping
格子ベクトルの長さを適当に増減させた場合も計算し、全エネルギーが最小になる点を求める。
** 結果 [#a861fedb]
|lattice const. (Bohr)|Etot (Hartree)|
|6.87|-59.7655504296820|
|6.88|-59.7655555602320|
|6.89|-59.7655530960053|
|CENTER:格子ベクトルの長さ (Bohr)|CENTER:全エネルギー (Hartree)|
|CENTER:6.87|CENTER:-59.7655504296820|
|CENTER:6.88|CENTER:-59.7655555602320|
|CENTER:6.89|CENTER:-59.7655530960053|
全エネルギーが最小になる点が見つかったら、その周りの5~11点の計算結果を用いて最小二乗法により格子定数を求める。