差分

このページの2つのバージョン間の差分を表示します。

--- 分子シミュレーション関連:amber:解析方法 [2018/02/05 10:27] – [水素結合解析] mumeiyamibito
+++ 分子シミュレーション関連:amber:解析方法 [2022/08/18 13:48] (現在) – mumeiyamibito
@@ 行 1: / 行 1: @@
 ====== 解析方法 ======
 ===== RMSD 解析 =====
-  * 指定された対象の RMSD の経時変化を解析する
+  * 指定された対象の RMSD の経時変化を解析する。
-  * ''cpptraj'' に同梱のパッケージで解析ができる
+  * $RMSD = \sqrt{\frac {\sum_{i = 1}^{N} (a_i - b_i)^2}{N}}$
+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要
-  * 対象 ''TARGET_MASK'' の RMSD を測定する\\ <code bash>
+  * 対象 ''TARGET_MASK'' の RMSD を測定する。\\ <code bash>
 $ cpptraj
 > parm test.prmtop
@@ 行 12: / 行 13: @@
 </code>
     * ''TARGET_MASK'': 対象のマスク
-    * ''first'': 最初のフレームをリファレンスにする
+    * ''first'': 最初のフレームをリファレンスにする。
     * ''OUTPUT'': 出力ファイル
     * ''time'': 各時間の刻み幅
@@ 行 18: / 行 19: @@
 ===== RMSF解析 =====
-  * 指定された対象の RMSF を解析する
+  * 指定された対象の RMSF を解析する。
   * RMSF はゆらぎのパラメータ
-  * ''cpptraj'' に同梱のパッケージで解析ができる
+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>
 $ cpptraj
@@ 行 32: / 行 33: @@
 ===== 水素結合解析 =====
-  * 指定された対象の水素結合数の経時変化を解析する
+  * 指定された対象の水素結合数の経時変化を解析する。
-  * ''cpptraj'' に同梱のパッケージで解析ができる
+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>
 $ cpptraj
@@ 行 87: / 行 88: @@
 ===== 水和水分析 =====
-  * 指定された対象の周囲の水分子数をカウントする
+  * 指定された対象の周囲の水分子数をカウントする。
-  * ''cpptraj'' に同梱のパッケージで解析ができる
+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>
 $ cpptraj
@@ 行 103: / 行 104: @@
 ===== 距離解析 =====
-  * 指定された対象間の距離の経時変化を解析する
+  * 指定された対象間の距離の経時変化を解析する。
-  * ''cpptraj'' に同梱のパッケージで解析ができる
+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要
-  * 対象1 ''TARGET1_MASK'' と対象2 ''TARGET2_MASK'' の重心 (中心) 間の距離を測定する\\ <code bash>
+  * 対象1 ''TARGET1_MASK'' と対象2 ''TARGET2_MASK'' の重心 (中心) 間の距離を測定する。\\ <code bash>
 $ cpptraj
 > parm test.prmtop
@@ 行 117: / 行 118: @@
     * ''OUTPUT'': 出力ファイル
     * ''geom'': 幾何学的中心で測定 (指定しなければ、重心で測定)
+===== 水分子等を一定距離で切り出した構造の抽出 =====
+  * 生体分子から一定距離にある水分子のみを残した構造をトラジェクトリから抽出する。\\ <code bash>
+$ cpptraj
+> parm test.prmtop
+> trajin test.inpcrd
+> mask MASK_BIOMOL<:DIST maskpdb OUTPUT.pdb
+> go
+</code>
+    * ''MASK_BIOMOL'': 生体分子のマスク
+    * ''DIST'': 生体分子からの距離 (例: '':5.0'' = 残基単位で切り出す; ''@5.0'' = 原子単位で切り出す)
+    * ''OUTPUT.pdb'': 出力するファイル (''OUTPUT.pdb.1'' のように、ファイル名の後にフレーム番号が付けられ、フレーム毎に出力される)
+    * カウンターイオンが含まれており、水分子のみを切り出したい場合は、''(:1-20)|((:1-20<:5.0)&WAT)'' のように、論理演算子を使う。
+    * 距離の算出には trajin のそれぞれのフレームが用いられるが、別のリファレンス構造を用いたい場合は、''reference'' と ''activeref REFERENCE_INDEX'' を指定する。\\ <code bash>
+$ cpptraj
+> parm test.prmtop
+> trajin test.inpcrd
+> reference test2.inpcrd
+> activeref 0
+> mask MASK_BIOMOL<:DIST maskpdb OUTPUT.pdb
+> go
+</code>
+===== カーブフィッティング =====
+  * 生体分子に関わらず、与えられたデータをカーブフィッティングする。\\ <code bash>
+curvefit DATASET_NAME EQUATION INITIAL_VAL out OUTPUT.dat tol TOLERANCE maxit ITER resultsout PARAM.dat
+</code>
+    * ''DATASET_NAME'': データセット名 (''readdata'' で読み込んだファイルのデータセット名)
+    * ''EQUATION'': フィッティングする関数
+    * ''INITIAL_VAL'': フィティング関数に初めに与える初期値
+    * ''OUTPUT.dat'': フィッティング後のデータセットを出力するファイル
+    * ''TOLERANCE'': フィッティングの閾値
+    * ''ITER'': フィッティングのループ回数
+    * ''PARAM.dat'': フィッティング結果を出力するファイル (フィッティング関数で決定した係数や誤差など)
+  * 例: 一次関数でフィッティングする場合\\ <code bash>
+$ cpptraj
+> readdata "data.dat"
+> runanalysis curvefit "data.dat" \
+	"FitY = A0 * X + A1" \
+	A0=10 A1=20 \
+	out "curve.dat" \
+	tol 0.001 \
+	maxit 50 \
+	resultsout "param.dat"
+> go
+</code>
+  * 参考サイト:
+    * [[https://amberhub.chpc.utah.edu/curvefit/ | curvefit – AMBER-hub]]
 {{tag>Linux 分子シミュレーション アプリケーション}}