第一原理分子動力学プログラム STATE Senri Wiki
開始行:
* 水素分子 [#s026a0c2]
このチュートリアルでは分子系について平面波カットオフとユ...
周期境界条件の下で孤立分子のシミュレーションを行う場合、...
ここでは水素分子を取り上げ、局所密度近似を用いて計算を行...
使用するポテンシャルはpot.H_lda1です。
** SCF 計算 [#l471888c]
- 入力ファイル (nfinp_1)
0 0 0 0 0 0 : H2 molecul...
5.0 10.0 1 2 2 : GMAX GMAXP...
1 0 : NUM_SPACE_...
10.0 10.0 10.0 90.0 90.0 90.0 : A B C ALPH...
1 1 1 1 1 1 : N1 N2 N3 M...
1 0 : NCORD NINV
-0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
1 0.15 1.00794 3 1 0.d0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 57200.00 0 : NMD1 NMD2 ...
3 1 : MIX_WAY MI...
0 8 0.7 : STARTING_M...
0.60 0.50 0.60 0.70 1.00 : DTIM1 DTIM...
100.0 4 1 0.10D-08 : DTIO IMDAL...
0.0010 0.50D+03 0 : WIDTH FORC...
ldapw91 1 : XCTYPE NSPIN
1.00 : DESTM
101 : NBZTYP
0 0 0 : DUMMY ...
0 0 0 : DUMMY ...
2 : NEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.0 : SM_N dopping
- 実行~
$ mpirun -np 4 ./STATE < nfinp_1 > nfout_1
- 結果~
STATEを実行すると以下のメッセージが表示されます。
*******************************************************...
* ...
* ...
* ...
* ****** ******** ** ******** *****...
* ******** ******** **** ******** *****...
* ** ** ** ** ** ** ...
* *** ** ******** ** *****...
* *** ** ********** ** *****...
* ** ** ** ** ** ** ...
* ******** ** ** ** ** *****...
* ****** ** VERSION 5.6.5 ** *****...
* RICS-AIST ...
* OSAKA UNIVERSITY ...
* ...
*******************************************************...
SCF計算が収束すると以下の全エネルギーおよび運動エネルギー...
TOTAL ENERGY AND ITS COMPONENTS
TOTAL ENERGY = -1.1375878...
FREE ENERGY = -1.1375878...
KINETIC ENERGY = 0.9682457...
HARTREE ENERGY = 0.7357007...
XC ENERGY = -0.6455573...
LOCAL ENERGY = -2.3137234...
NONLOCAL ENERGY = -0.0233694...
EWALD ENERGY = 0.1411159...
PC ENERGY = 0.0000000...
ENTROPIC ENERGY = 0.0000000...
全エネルギー以降にフェルミエネルギー、Kohn-Sham固有値と占...
FERMI ENERGY = ...
NKP= 1 NGP= 2103 K=( 0.00000 0.00000 0.000...
EIGEN VALUE
-0.37074 -0.00963
OCCUPATION
1.00000 0.00000
さらに原子にかかる力が以下のように表示されます。
ATOM COORDINATES ...
MD: 1
MD: 1 H -0.710000 0.000000 0.000000 -0.012...
MD: 2 H 0.710000 0.000000 0.000000 0.012...
** 収束性の検証 [#f30e2cb1]
平面波基底を用いた計算の場合、波動関数の平面波展開のため...
ノルム保存擬ポテンシャルを用い、なおかつ非線形内殻補正を...
しかしながらウルトラソフト擬ポテンシャルの場合は局在した...
STATEでは波動関数を展開する平面波のカットオフ波数ベクトル...
電荷密度のカットオフ波数ベクトルはGMAXPで指定され、GMAXP...
GMAXとGMAXPを変えながら全エネルギー計算を実行することで以...
[GMAXP=15 (E_cut^den=225 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12781490
4.5 -1.13418820
5.0 -1.13697688
5.5 -1.13816998
6.0 -1.13847300
6.5 -1.13850127
7.0 -1.13851867
7.5 -1.13857102
[GMAXP=20 (E_cut^den=400 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12780183
4.5 -1.13417130
5.0 -1.13695809
5.5 -1.13815036
6.0 -1.13845313
6.5 -1.13848136
7.0 -1.13849876
7.5 -1.13855109
[GMAXP=25 (E_cut^den=625 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12780130
4.5 -1.13417074
5.0 -1.13695751
5.5 -1.13814979
6.0 -1.13845255
6.5 -1.13848078
7.0 -1.13849818
7.5 -1.13855052
[GMAXP=30 (E_cut^den=900 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12780182
4.5 -1.13417124
5.0 -1.13695801
5.5 -1.13815027
6.0 -1.13845303
6.5 -1.13848127
7.0 -1.13849867
7.5 -1.13855100
GMAX=6とGMAXP=20で全エネルギーは0.1 mRyの精度で収束してい...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
~
孤立分子のシミュレーションを実行する場合、カットオフエネ...
GMAX=6とGMAXP=20を用いて全エネルギーをユニットセル(立方...
#A(Bohr) ETOT(Ha)
5.00 -1.31378028
10.00 -1.13845313
15.00 -1.13662806
20.00 -1.13661187
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
上の結果から立方体のユニットセルの辺の長さが15 Bohrとなっ...
大きい分子の場合はより大きいユニットセルを用いる必要があ...
** 構造最適化 (1) [#d2e5e2e9]
二原子分子の場合、原子間距離(ボンド長)を変えながら全エ...
ここではGMAX=6、GMAXP=20、辺の長さが10 Bohrの立方体のユニ...
# Bond length (Bohr) Total energy (Hartree)
1.28000000 -1.13291015
1.30000000 -1.13434083
1.32000000 -1.13555229
1.34000000 -1.13655937
1.36000000 -1.13737569
1.38000000 -1.13801370
1.40000000 -1.13848494
1.42000000 -1.13880026
1.44000000 -1.13896985
1.46000000 -1.13900337
1.48000000 -1.13890994
1.50000000 -1.13869810
1.52000000 -1.13837576
1.54000000 -1.13795024
1.56000000 -1.13742831
1.58000000 -1.13681622
1.60000000 -1.13611994
1.62000000 -1.13534513
1.64000000 -1.13449712
1.66000000 -1.13358088
得られた全エネルギーを6次の多項式にフィットすることで平衡...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
** 構造最適化 (2) [#v236650a]
より複雑な系では手動で上記のような構造最適化を行うのは大...
そのような場合にはHellmann-Feynman力を求めて構造最適化を...
Hellmann-Feynmann力を用いて構造最適化を行う場合、力の閾値...
この例では最適化のアルゴリムとしてgeneralized direct inve...
0 0 0 0 0 0 : H2 molecul...
6.0 20.0 1 2 2 : GMAX GMAXP...
1 0 : NUM_SPACE_...
10.0 10.0 10.0 90.0 90.0 90.0 : A B C ALPH...
1 1 1 1 1 1 : N1 N2 N3 M...
1 0 : NCORD NINV
-0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
1 0.15 1.00794 3 1 0.d0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 57200.00 0 : NMD1 NMD2 ...
3 1 : MIX_WAY MI...
0 8 0.7 : STARTING_M...
0.60 0.50 0.60 0.70 1.00 : DTIM1 DTIM...
100.0 4 1 0.10D-08 : DTIO IMDAL...
0.0010 0.50D-03 0 : WIDTH FORC...
ldapw91 1 : XCTYPE NSPIN
1.00 : DESTM
101 : NBZTYP
0 0 0 : DUMMY
0 0 0 : DUMMY
2 : NEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.0 : SM_N dopping
得られた平衡ボンド長は1.455 Bohr (0.77 Angstrom)で手動最...
** 演習問題 [#dc85ed87]
- ボンド長のユニットセルサイズおよびカットオフ依存生を調...
終了行:
* 水素分子 [#s026a0c2]
このチュートリアルでは分子系について平面波カットオフとユ...
周期境界条件の下で孤立分子のシミュレーションを行う場合、...
ここでは水素分子を取り上げ、局所密度近似を用いて計算を行...
使用するポテンシャルはpot.H_lda1です。
** SCF 計算 [#l471888c]
- 入力ファイル (nfinp_1)
0 0 0 0 0 0 : H2 molecul...
5.0 10.0 1 2 2 : GMAX GMAXP...
1 0 : NUM_SPACE_...
10.0 10.0 10.0 90.0 90.0 90.0 : A B C ALPH...
1 1 1 1 1 1 : N1 N2 N3 M...
1 0 : NCORD NINV
-0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
1 0.15 1.00794 3 1 0.d0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 57200.00 0 : NMD1 NMD2 ...
3 1 : MIX_WAY MI...
0 8 0.7 : STARTING_M...
0.60 0.50 0.60 0.70 1.00 : DTIM1 DTIM...
100.0 4 1 0.10D-08 : DTIO IMDAL...
0.0010 0.50D+03 0 : WIDTH FORC...
ldapw91 1 : XCTYPE NSPIN
1.00 : DESTM
101 : NBZTYP
0 0 0 : DUMMY ...
0 0 0 : DUMMY ...
2 : NEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.0 : SM_N dopping
- 実行~
$ mpirun -np 4 ./STATE < nfinp_1 > nfout_1
- 結果~
STATEを実行すると以下のメッセージが表示されます。
*******************************************************...
* ...
* ...
* ...
* ****** ******** ** ******** *****...
* ******** ******** **** ******** *****...
* ** ** ** ** ** ** ...
* *** ** ******** ** *****...
* *** ** ********** ** *****...
* ** ** ** ** ** ** ...
* ******** ** ** ** ** *****...
* ****** ** VERSION 5.6.5 ** *****...
* RICS-AIST ...
* OSAKA UNIVERSITY ...
* ...
*******************************************************...
SCF計算が収束すると以下の全エネルギーおよび運動エネルギー...
TOTAL ENERGY AND ITS COMPONENTS
TOTAL ENERGY = -1.1375878...
FREE ENERGY = -1.1375878...
KINETIC ENERGY = 0.9682457...
HARTREE ENERGY = 0.7357007...
XC ENERGY = -0.6455573...
LOCAL ENERGY = -2.3137234...
NONLOCAL ENERGY = -0.0233694...
EWALD ENERGY = 0.1411159...
PC ENERGY = 0.0000000...
ENTROPIC ENERGY = 0.0000000...
全エネルギー以降にフェルミエネルギー、Kohn-Sham固有値と占...
FERMI ENERGY = ...
NKP= 1 NGP= 2103 K=( 0.00000 0.00000 0.000...
EIGEN VALUE
-0.37074 -0.00963
OCCUPATION
1.00000 0.00000
さらに原子にかかる力が以下のように表示されます。
ATOM COORDINATES ...
MD: 1
MD: 1 H -0.710000 0.000000 0.000000 -0.012...
MD: 2 H 0.710000 0.000000 0.000000 0.012...
** 収束性の検証 [#f30e2cb1]
平面波基底を用いた計算の場合、波動関数の平面波展開のため...
ノルム保存擬ポテンシャルを用い、なおかつ非線形内殻補正を...
しかしながらウルトラソフト擬ポテンシャルの場合は局在した...
STATEでは波動関数を展開する平面波のカットオフ波数ベクトル...
電荷密度のカットオフ波数ベクトルはGMAXPで指定され、GMAXP...
GMAXとGMAXPを変えながら全エネルギー計算を実行することで以...
[GMAXP=15 (E_cut^den=225 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12781490
4.5 -1.13418820
5.0 -1.13697688
5.5 -1.13816998
6.0 -1.13847300
6.5 -1.13850127
7.0 -1.13851867
7.5 -1.13857102
[GMAXP=20 (E_cut^den=400 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12780183
4.5 -1.13417130
5.0 -1.13695809
5.5 -1.13815036
6.0 -1.13845313
6.5 -1.13848136
7.0 -1.13849876
7.5 -1.13855109
[GMAXP=25 (E_cut^den=625 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12780130
4.5 -1.13417074
5.0 -1.13695751
5.5 -1.13814979
6.0 -1.13845255
6.5 -1.13848078
7.0 -1.13849818
7.5 -1.13855052
[GMAXP=30 (E_cut^den=900 Ry)]
#GMAX ETOT (Ha)
4.0 -1.12780182
4.5 -1.13417124
5.0 -1.13695801
5.5 -1.13815027
6.0 -1.13845303
6.5 -1.13848127
7.0 -1.13849867
7.5 -1.13855100
GMAX=6とGMAXP=20で全エネルギーは0.1 mRyの精度で収束してい...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
~
孤立分子のシミュレーションを実行する場合、カットオフエネ...
GMAX=6とGMAXP=20を用いて全エネルギーをユニットセル(立方...
#A(Bohr) ETOT(Ha)
5.00 -1.31378028
10.00 -1.13845313
15.00 -1.13662806
20.00 -1.13661187
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
上の結果から立方体のユニットセルの辺の長さが15 Bohrとなっ...
大きい分子の場合はより大きいユニットセルを用いる必要があ...
** 構造最適化 (1) [#d2e5e2e9]
二原子分子の場合、原子間距離(ボンド長)を変えながら全エ...
ここではGMAX=6、GMAXP=20、辺の長さが10 Bohrの立方体のユニ...
# Bond length (Bohr) Total energy (Hartree)
1.28000000 -1.13291015
1.30000000 -1.13434083
1.32000000 -1.13555229
1.34000000 -1.13655937
1.36000000 -1.13737569
1.38000000 -1.13801370
1.40000000 -1.13848494
1.42000000 -1.13880026
1.44000000 -1.13896985
1.46000000 -1.13900337
1.48000000 -1.13890994
1.50000000 -1.13869810
1.52000000 -1.13837576
1.54000000 -1.13795024
1.56000000 -1.13742831
1.58000000 -1.13681622
1.60000000 -1.13611994
1.62000000 -1.13534513
1.64000000 -1.13449712
1.66000000 -1.13358088
得られた全エネルギーを6次の多項式にフィットすることで平衡...
#ref(http://www-cp.prec.eng.osaka-u.ac.jp/puki_state/grap...
** 構造最適化 (2) [#v236650a]
より複雑な系では手動で上記のような構造最適化を行うのは大...
そのような場合にはHellmann-Feynman力を求めて構造最適化を...
Hellmann-Feynmann力を用いて構造最適化を行う場合、力の閾値...
この例では最適化のアルゴリムとしてgeneralized direct inve...
0 0 0 0 0 0 : H2 molecul...
6.0 20.0 1 2 2 : GMAX GMAXP...
1 0 : NUM_SPACE_...
10.0 10.0 10.0 90.0 90.0 90.0 : A B C ALPH...
1 1 1 1 1 1 : N1 N2 N3 M...
1 0 : NCORD NINV
-0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
0.710 0.000 0.000 1 1 1 : CPS(1:3) I...
1 0.15 1.00794 3 1 0.d0 : IATOMN ALF...
0 0 0 0 0 : ICOND INIP...
0 1 : IPRE IPRI
200 200 0 57200.00 0 : NMD1 NMD2 ...
3 1 : MIX_WAY MI...
0 8 0.7 : STARTING_M...
0.60 0.50 0.60 0.70 1.00 : DTIM1 DTIM...
100.0 4 1 0.10D-08 : DTIO IMDAL...
0.0010 0.50D-03 0 : WIDTH FORC...
ldapw91 1 : XCTYPE NSPIN
1.00 : DESTM
101 : NBZTYP
0 0 0 : DUMMY
0 0 0 : DUMMY
2 : NEG
1 : NEXTST
0 : DUMMY
2 : IMSD
0 : EVALUATE_E...
0 : NPDOSAO
0 0.0 : SM_N dopping
得られた平衡ボンド長は1.455 Bohr (0.77 Angstrom)で手動最...
** 演習問題 [#dc85ed87]
- ボンド長のユニットセルサイズおよびカットオフ依存生を調...
ページ名: