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

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• バックアップ を表示
• 格子定数の計算 へ行く。
  • 1 (2015-06-19 (金) 21:41:21)
  • 2 (2015-06-19 (金) 22:52:51)
  • 3 (2015-06-20 (土) 05:09:29)
  • 4 (2017-11-30 (木) 20:22:39)

* 計算方法 [#y9329d23]
** 例1: bulk Cu [#b2afc8e1]
*** 入力ファイル [#v64bf6f8]
  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]
|CENTER:格子ベクトルの長さ (Bohr)|CENTER:全エネルギー (Hartree)|h
|CENTER:...|CENTER:...|
|CENTER:6.87|CENTER:-59.7655504296820|
|CENTER:&color(red){6.88};|CENTER:-59.76555&color(red){55602320};|
|CENTER:6.89|CENTER:-59.7655530960053|
|CENTER:...|CENTER:...|
全エネルギーが最小になる点が見つかったら、その周りの5~11点の計算結果を用いて最小二乗法により格子定数を求める。

* 構成定数の汎関数依存性 [#v9ee73f1]
||PBE|revPBE|vdW1|optPBE|optB88|optB86b|vdW2|rev-vdW2|備考|h
|Pt|3.946||4.007|3.966|3.957|3.926|4.090|3.929|ver.5.5.4, gm6gmp20|