GMX中的特殊分組

Gromacs中的角度中剩,距離等分析工具的使用都涉及到特殊的分組方法舶胀,能不能以合適的格式分組是計算成功的關鍵。

選取分組主要有兩種方法唤锉,一是直接在命令中使用Select語法進行分組冬骚,二是建立相應的索引文件椅贱。

第一種方法更加靈活并支持動態(tài)選區(qū),可以在命令中直接用''括起來放入其中只冻。當然前提是整個體系中原子命名等信息不能有重復庇麦。例如在一條碳鏈中,每個碳最好分別命名為"C1""C2"...而不能統(tǒng)一命名為"C"喜德,這樣將很難通過Select語法來甄別山橄。例如,要是要計算一條碳鏈與z軸之間的夾角舍悯,則可以使用命令gmx gangle -f ..xtc -s ..tpr -n ..ndx -g1 vector -group1 'name C1 C5 permute 1 2' -g2 z -oav -oall -oh航棱。此處的'name C1 C5 permute 1 2'意為從名為C1的原子指向名為C5原子的向量,1 2表明向量的方向從第一組指向第二組萌衬。有關Select詳細的語法和使用實例在李老師的博文中有詳細的介紹饮醇,感興趣的可以訪問這里GROMACS選區(qū)(selection)語法及用法

若是沒有明智地提前將原子名稱區(qū)分好秕豫,或是即使用Select語法也無法選定特定的組朴艰,就只好通過建立相應的索引文件了。在distance與gangle等命令中所用的組往往需要選定兩個位置,這兩個位置應當在索引文件中成對出現(xiàn)呵晚。若是作為向量,則是從第一個向量指向第二個向量沫屡。同理饵隙,對于需要三個位置的組,則應三個為一行沮脖。下面是個簡單的例子:
[ Vector1 ] 1 5 2 6 3 7 ... [ Plane1 ] 1 3 5 2 4 6 ...

這樣金矛,在做好相應的索引文件后,就可以直接在命令執(zhí)行時選擇相應的組就可以了勺届。這時驶俊,上述的例子就可以寫成: gmx gangle -f ..xtc -s ..tpr -n ..ndx -g1 vector -group1 -g2 z -oav -oall -oh,在執(zhí)行命令后會出現(xiàn)選擇分組的提示免姿,根據(jù)提示選擇相應的組就可以了饼酿。

另外,在建立特定索引組的過程中胚膊,可以尋求出序號之間的規(guī)律故俐,寫一個小程序直接生成疮蹦。索引文件的格式很簡單窖剑,保證組名和數(shù)字就好。并不是非得使用gmx make_ndx工具的白魂。

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
  • 序言:七十年代末喻犁,一起剝皮案震驚了整個濱河市槽片,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌肢础,老刑警劉巖还栓,帶你破解...
    沈念sama閱讀 222,378評論 6 516
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異传轰,居然都是意外死亡蝙云,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,970評論 3 399
  • 文/潘曉璐 我一進店門路召,熙熙樓的掌柜王于貴愁眉苦臉地迎上來勃刨,“玉大人,你說我怎么就攤上這事股淡∩硪” “怎么了?”我有些...
    開封第一講書人閱讀 168,983評論 0 362
  • 文/不壞的土叔 我叫張陵唯灵,是天一觀的道長贾铝。 經(jīng)常有香客問我,道長,這世上最難降的妖魔是什么垢揩? 我笑而不...
    開封第一講書人閱讀 59,938評論 1 299
  • 正文 為了忘掉前任玖绿,我火速辦了婚禮,結果婚禮上叁巨,老公的妹妹穿的比我還像新娘斑匪。我一直安慰自己,他們只是感情好锋勺,可當我...
    茶點故事閱讀 68,955評論 6 398
  • 文/花漫 我一把揭開白布蚀瘸。 她就那樣靜靜地躺著,像睡著了一般庶橱。 火紅的嫁衣襯著肌膚如雪贮勃。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 52,549評論 1 312
  • 那天苏章,我揣著相機與錄音寂嘉,去河邊找鬼。 笑死枫绅,一個胖子當著我的面吹牛垫释,可吹牛的內容都是我干的。 我是一名探鬼主播撑瞧,決...
    沈念sama閱讀 41,063評論 3 422
  • 文/蒼蘭香墨 我猛地睜開眼棵譬,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了预伺?” 一聲冷哼從身側響起订咸,我...
    開封第一講書人閱讀 39,991評論 0 277
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎酬诀,沒想到半個月后脏嚷,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 46,522評論 1 319
  • 正文 獨居荒郊野嶺守林人離奇死亡瞒御,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 38,604評論 3 342
  • 正文 我和宋清朗相戀三年父叙,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片肴裙。...
    茶點故事閱讀 40,742評論 1 353
  • 序言:一個原本活蹦亂跳的男人離奇死亡趾唱,死狀恐怖,靈堂內的尸體忽然破棺而出蜻懦,到底是詐尸還是另有隱情甜癞,我是刑警寧澤,帶...
    沈念sama閱讀 36,413評論 5 351
  • 正文 年R本政府宣布宛乃,位于F島的核電站悠咱,受9級特大地震影響蒸辆,放射性物質發(fā)生泄漏。R本人自食惡果不足惜析既,卻給世界環(huán)境...
    茶點故事閱讀 42,094評論 3 335
  • 文/蒙蒙 一躬贡、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧眼坏,春花似錦拂玻、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,572評論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽擂仍。三九已至囤屹,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間逢渔,已是汗流浹背肋坚。 一陣腳步聲響...
    開封第一講書人閱讀 33,671評論 1 274
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留肃廓,地道東北人智厌。 一個月前我還...
    沈念sama閱讀 49,159評論 3 378
  • 正文 我出身青樓,卻偏偏與公主長得像盲赊,于是被迫代替她去往敵國和親铣鹏。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 45,747評論 2 361

推薦閱讀更多精彩內容

  • 最近在寫個性化推薦的論文哀蘑,經(jīng)常用到Python來處理數(shù)據(jù)诚卸,被pandas和numpy中的數(shù)據(jù)選取和索引問題繞的比較...
    shuhanrainbow閱讀 4,568評論 6 19
  • 一. Java基礎部分.................................................
    wy_sure閱讀 3,814評論 0 11
  • 1. Java基礎部分 基礎部分的順序:基本語法,類相關的語法绘迁,內部類的語法合溺,繼承相關的語法,異常的語法缀台,線程的語...
    子非魚_t_閱讀 31,664評論 18 399
  • 文/宇坤字數(shù)1.6K 微信運營手冊V1.0系統(tǒng)整理棠赛,暫定7文,本月底前發(fā)文完畢膛腐,大致計劃是: 1睛约、框架:名稱、頭像...
    大坤筆記閱讀 1,131評論 3 14
  • 夜色那么凄美哲身,猶如轉瞬即逝的人生痰腮,苦短卻韻味十足。不知還要像斷線的風箏一樣漂泊多久律罢,我才能徹底的掙脫那根寂寥的線...
    竹鴻初閱讀 194評論 0 0