* 計算方法 [#y9329d23] ** 例1: bulk Cu [#b2afc8e1] *** 入力ファイル [#v64bf6f8] * 計算例: 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] ** 構成定数の汎関数依存性 [#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|