2016-03-28 20 views
0

私は水の箱にミセルの分子動力学の軌道を持っています。各フレームについて、ミセルの重心(COM)を求め、COMからウォーターボックスの側面までの水分子の密度を計算したいと考えています。誰かがそれを行うための最善の方法を知っていますか?誰かがこのような問題を以前に処理したかどうか疑問に思っていました。MD軌道上の放射密度分布

おかげ Sadegh

+0

こんにちは、Sadegh!あなたは写真で何を意味しているのかを教えてください。そしてCOMから水箱の側面までの水分子の密度を計算してください。脂質は何を使用していますか? (例で正しい名前を書く必要がある) – XuMuK

答えて

0

まず - 便利なコマンドの多くとその使用方法はhereに書かれています。

あなたはタイプによって選択するよりも、あなたの脂質を(私はリン脂質から構築されたミセルを想定)を選択する必要があり

- 脂質:

set sel1 [atomselect top resname DPC] 

そしてより:

center_of_mass $sel1 
set sel1 [atomselect top "lipids"] 

またはresnameはによって

水分子の密度を計算するCOMから 水ボックス

の両側にSは私はあなたが話している内容を正確に把握していない - ミセル内のいくつかの水の分子があるのですか?または水槽はミセルから一定の距離内の水分子ですか?いくつかの情報を提供してください(説明のあるシステムの方がいいです)。

+0

こんにちはXuMuK!ミセルは水分子によって囲まれた立方体の水箱の中にある。また、ミセル内部には水分子が浸透している。ミセルの外側にある水分子の数やミセルの外側にある水分子の数を調べるには、ミセルCOMに対する水分子のラジカル密度関数を計算する必要があります。あなたは、ミセルCOMの周りにいくつかの球状層を想像することができます。私はこれらの層のそれぞれの中にある水分子の数をrの関数として求めたいと思います。しかし、LOOSパッケージにはプラグインがあることがわかりました。 – s110faramarzi

+0

@ s110faramarziしたがって、resname DPPC "の中の' atomselect top "を使うべきです" http://www.ks.uiuc.edu/Research/vmd/vmd-1.7.1/ug/node108.htmlまた:http:///www.ks.uiuc.edu/Research/vmd/allversions/vmd_faq.htmlこれを行うには、 'sqr(x-3)+ sqr(y-4)+ sqr(z + 5)< (x、y、z)=(3,4、-5)の2Å以内のすべての原子を選択する。私たちはいつかこの選択を簡単にするために点x y zの距離内に選択項目を追加することがあります。 – XuMuK

関連する問題