第一原理分子動力学プログラム STATE Senri Wiki

計算例: bulk Cu †

入力ファイル †

 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

格子ベクトルの長さを適当に増減させた場合も計算し、全エネルギーが最小になる点を求める。

結果 †

全エネルギーが最小になる点が見つかったら、その周りの5~11点の計算結果を用いて最小二乗法により格子定数を求める。

構成定数の汎関数依存性 †