- 追加された行はこの色です。
- 削除された行はこの色です。
*Molecular dynamics [#c3820ec9]
分子動力学計算の実行方法
分子動力学シミュレーションの実行方法
**HOWTO [#refd1c79]
- 分子動力学の実行を指定
--IMDALG=-1 or IMDALG=-2を選択。
~通常はIMDALG=-1を指定する。
--IMDALG=-1を指定するとTRAJECTORY,ENERGIES等の有用なデータが出力される。
--IMDALG=2はquenching
--IMDALG=3はVibration analysis
--IMDALG=4はGDIIS
--IMDALG=5はGDIISTS 制約条件の入れ方を要追記
--IMDALG=6はNEB
--IMDALG=7はNEBwithClimbingImage
*** 分子動力学計算 [#f5253999]
- 分子動力学シミュレーションの実行を指定
--IMDALG=-1あるいはIMDALG=-2を選択
~通常はIMDALG=-1を指定する。TRAJECTORY、ENERGIES等が生成される。
- [[計算結果のアニメーション化>計算コードの使い方/STATE/原子構造のアニメーション化]]
-IMDALG=-1(有限温度分子動力学)の場合、固定フォーマット形式の最終行(...SM_dopping)の後に、有限温度分子動力学の条件を付け加える。
2 : mvelsc
300.0D0 30.0D0 10 1.0D0 : tempw,tolp,nroll,anneal
300.0D0 8 15 1 : wnosep,nhc,nosy,ndrt
5.0D-8 : frict
0 : mcnstr
- MVELSC: 温度(速度)の制御法を指定
| mvelsc | velocity scaling |
| 0 | micro cano |
| 1 | Simulated annealing |
| 2 | simple velocity rescaling |
| 10 | Nose-Hoover |
- TEMPW: 温度 (K)
- TOLP: 温度の閾値 (K) Velocity rescalingの場合に目的の温度よりTOLPずれた場合に温度を目標の温度にスケールする
- NROLL: Rolling averageのインターバル 通常10が推奨される
- ANNEAL: アニーリングの因子 Simulated annealingの際に用いられる
- WNOSEP: 特徴的な系の振動数 (cm^-1)
- NHC: Nose-Hoover chainの長さ (NHCの場合4, GGMTの場合2が推奨される)
- NOSY: Suzuki-Yoshida積分の次数
- NDRT: 熱浴変数についての積分の数 通常NDRT=1で十分である
- FRICT: Langevin分子動力学シミュレーションの際の摩擦係数
- MCNSTR: 拘束条件の数
MCNSTRが0よりも大きい場合、拘束条件の数だけ行を追加する
DIST 1 2 1.83401~
BEND 1 2 3 104.23195~
DIHED 1 2 3 4 60.0~
ここで拘束条件は
-- DIST: 距離
-- BEND: 角度
-- DIHED: 二面角
が実装されており、2-5 columnの数字(column数は拘束条件に依存する)で拘束に関係する原子の番号を指定し、最後の実数で距離あるいは角度を指定する。
- 制約条件付きMD
*** 拘束条件付き分子動力学シミュレーション [#vdad45f2]
以下のオプションを追加することによっても拘束条件付き分子動力学シミュレーションを実行することが可能である。
&COSTRAINT
NMCNSTR 1 (3の誤りか?)
DIST 6 7 4.102
BEND 1 2 3 104.23195
DIHED 1 2 3 4 60.0
NMCNSTR 3
DIST 6 7 4.102
BEND 1 2 3 104.23195
DIHED 1 2 3 4 60.0
&END
DIST は原子間距離~
BENDは結合間の角度~
DIHEDは2面角(要確認)~
&CONTRAINT~&ENDと書けばあらゆるMDの方法(quenching等)でも機能する。
NMCNSTR: 拘束条件の数
DIST: 原子間距離~
BEND: 結合間の角度~
DIHED: 2面角~
このオプションは拘束条件付き構造最適化を念頭に置いており、有限温度分子動力学シミュレーションを実行する際にCONSTRAINTファイルが生成されない(要修正)。
**分子動力学シミュレーション結果の可視化 [#ye0162b9]
- [[計算結果のアニメーション化>計算コードの使い方/STATE/原子構造のアニメーション化]]
- 温度制御の方法の指定
--MVELSCの設定。頻繁に使われるのは以下のものである。~
0: ミクロカノニカル アンサンブル(温度制御なし)~
2: velocity rescaling~
10: Nose-Hoover
~
- 有限温度のMDの場合、&CONTRAINT~&ENDがなくても機能する。
nfinpの最後に~
DIST 1 2 1.83401~
BEND 1 2 3 104.23195~
DIHED 1 2 3 4 60.0~
等と記述する。~
DIST は原子間距離~
BENDは結合間の角度~
DIHEDは2面角(要確認)~