More than 3 years have passed since last update.

PLUMED入門その5 仮想原子を使ってCVを指定する

Posted at 2021-05-07

はじめに

本シリーズは下記チュートリアルを実施したメモです。

本記事はPLUMED入門の第４回です

過去回はこちら↓

Exercise 3 : 仮想原子

ときおり、CVを計算するときに、たくさんの原子の位置を使えないことがある。
その代わりに、他の原子の位置に基き生成された仮想原子を定義したいかもしれない。
例えば、原子集団の重心（center of mass; COM）や幾何中心（geometric center; CENTER）を使いたいかもしれない。

このエクササイズでは、仮想原子を指定してCVを定義する。
以下のインプットファイルを見てみよう。

# Geometric center of first residue
first: CENTER ATOMS=1,2,3,4,5,6,7,8
# Geometric center of last residue
last: CENTER ATOMS=427-436

# Distance between centers of first and last residues, with PBCs
d1: DISTANCE ATOMS=first,last
# Distance between centers of first and last residues, without PBCs
d2: DISTANCE ATOMS=first,last NOPBC

# Print the two distances on COLVAR file every step
PRINT ARG=__FILL__ FILE=COLVAR STRIDE=__FILL__

このインプットは、

1~8の原子の幾何中心を仮想原子"first"として定義する
427〜436の原子の幾何中心を仮想原子"last"として定義する
2つの仮想原子間の距離を計算し、変数d1として保存する
2つの仮想原子間の距離を計算し、変数d2として保存する。ただし、PBCは無視する
COLVARにd1,d2を出力する

firstとlastの原子の指定が異なる方法で行われている。

以下が私が作成したインプットファイルです。

# Geometric center of first residue
first: CENTER ATOMS=1,2,3,4,5,6,7,8
# Geometric center of last residue
last: CENTER ATOMS=427-436

# Distance between centers of first and last residues, with PBCs
d1: DISTANCE ATOMS=first,last
# Distance between centers of first and last residues, without PBCs
d2: DISTANCE ATOMS=first,last NOPBC

# Print the two distances on COLVAR file every step
PRINT ARG=d1,d2 FILE=COLVAR STRIDE=1

PLUMEDを実行する。

plumed driver --mf_xtc traj-broken.xtc --plumed plumed.dat

仮想原子は437,438番目の原子として追加されたよう。

PLUMED: Number of atoms: 436
...

PLUMED: Action CENTER
PLUMED:   with label first
PLUMED:   serial associated to this virtual atom is 437
PLUMED:   of atoms:
PLUMED:  1 2 3 4 5 6 7 8
...
PLUMED: Action CENTER
PLUMED:   with label last
PLUMED:   serial associated to this virtual atom is 438
PLUMED:   of atoms:
PLUMED:  427 428 429 430 431 432 433 434 435 436

COLVARを見ると、

# ! FIELDS time d1 d2
 0.000000 5.841553 5.841553
 1.000000 5.295298 8.311418
 2.000000 5.658301 8.467723
 3.000000 5.192408 7.894694
 4.000000 5.130006 7.779987
 5.000000 4.187504 7.114014

となっていて、d1とd2は異なる値となることがある。

whole.xtcで計算した場合でも

# ! FIELDS time d1 d2
 0.000000 5.841553 5.841553
 1.000000 5.295298 8.311418
 2.000000 5.658301 8.467723
 3.000000 5.192408 7.894694
 4.000000 5.130006 7.779987
 5.000000 4.187504 7.114014

であり、先程と同様にd1とd2は異る。

gnuplotでグラフを書くと、