確定起始構型
進行分子動力學模擬的第一步是確定起始構型,一個能量較低的起始構型是進行分子模擬的基礎,一般分子的起始構型主要來自實驗數據或量子化學計算。 在確定起始構型之后要賦予構成分子的各個原子速度,這一速度是根據波爾茲曼分布隨機生成的,由于速度的分布符合波爾茲曼統(tǒng)計,因此在這個階段,體系的溫度是恒定的。另外,在隨機生成各個原子的運動速度之后須進行調整,使得體系總體在各個方向上的動量之和為零,即保證體系沒有平動位移。
進入平衡相
由上一步確定的分子組建平衡相,在構建平衡相的時候會對構型、溫度等參數加以監(jiān)控。
進入生產相
進入生產相之后體系中的分子和分子中的原子開始根據初始速度運動,可以想象其間會發(fā)生吸引、排斥乃至碰
分子動力學
撞,這時就根據牛頓力學和預先給定的粒子間相互作用勢來對各個粒子的運動軌跡進行計算,在這個過程中,體系總能量不變,但分子內部勢能和動能不斷相互轉化,從而體系的溫度也不斷變化,在整個過程中,體系會遍歷勢能面上的各個點,計算的樣本正是在這個過程中抽取的。步驟
以下是做模擬的一般性步驟,具體的步驟和過程依賴于確定的系統(tǒng)或者是軟件,但這不影響我們把它當成一個入門指南:
1)首先我們需要對我們所要模擬的系統(tǒng)做一個簡單的評估, 三個問題是我們必須要明確的:
做什么(what to do)為什么做(why to do)怎么做(how to do)
2)選擇合適的模擬工具,大前提是它能夠實現你所感興趣的目標,這需要你非常謹慎的查閱文獻,看看別人用這個工具都做了些什么,有沒有和你相關的,千萬不要做到一半才發(fā)現原來這個工具根本就不能實現你所感興趣的idea,切記! 考慮1:軟件的選擇,這通常和軟件主流使用的力場有關,而軟件本身就具體一定的偏向性,比如說,做蛋白體系,Gromacs,Amber,Namd均可;做DNA, RNA體系,首選肯定是Amber;做界面體系,Dl_POLY比較強大,另外做材料體系,Lammps會是一個不錯的選擇
考慮2:力場的選擇。力場是來描述體系中最小單元間的相互作用的,是用量化等方法計算擬合后生成的經驗式,有人會嫌它粗糙,但是它確確實實給我們模擬大系統(tǒng)提供了可能,只能說關注的切入點不同罷了。常見的有三類力場:全原子力場,聯合力場,粗?;?;當然還有所謂第一代,第二代,第三代力場的說法,這里就不一一列舉了。
再次提醒注意:必須選擇適合于我們所關注體系和我們所感興趣的性質及現象的力場。 3)通過實驗數據或者是某些工具得到體系內的每一個分子的初始結構坐標文件,之后,我們需要按我們的想法把這些分子按照一定的規(guī)則或是隨機的排列在一起,從而得到整個系統(tǒng)的初始結構,這也是我們模擬的輸入文件。
4)結構輸入文件得到了,我們還需要力場參數輸入文件,也就是針對我們系統(tǒng)的力場文件,這通常由所選用的力場決定,比如鍵參數和非鍵參數等勢能函數的輸入參數。 分子動力學
5)體系的大小通常由你所選用的box大小決定,我們必須對可行性與合理性做出評估,從而確定體系的大小,這依賴于具體的體系,這里不細說了。6)由于初始構象可能會存在兩個原子挨的太近的情況(稱之為bad contact),所以需要在正式模擬開始的第一步進行體系能量最小化,比較常用的能量最小化有兩種,最速下降法和共軛梯度法,最速下降法是快速移除體系內應力的好方法,但是接近能量極小點時收斂比較慢,而共軛梯度法在能量極小點附近收斂相對效率高一些,所有我們一般做能量最小化都是在最速下降法優(yōu)化完之后再用共軛梯度法優(yōu)化,這樣做能有效的保證后續(xù)模擬的進行。 7)以平衡態(tài)模擬為例,你需要設置適當的模擬參數,并且保證這些參數設置和力場的產生相一致,舉個簡單的例子,gromos力場是用的范德華勢雙截斷來定范德華參數的,若你也用gromos力場的話也應該用雙截斷來處理范德華相互作用。常見的模擬思路是,先在NVT下約束住你的溶質(劑)做限制性模擬,這是一個升溫的過程,當溫度達到你的設定后, 接著做NPT模擬,此過程將調整體系的壓強進而使體系密度收斂。 經過一段時間的平衡模擬,在確定系統(tǒng)弛豫已經完全消除之后,就可以開始取數據了。如何判斷體系達到平衡,這個問題是比較技術性的問題,簡單的講可以通過以下幾種方式,一,看能量(勢能,動能和總能)是否收斂;二,看系統(tǒng)的壓強,密度等等是否收斂;三看系統(tǒng)的RMSD是否達到你能接受的范圍,等等。
8)運行足夠長時間的模擬以確定我們所感興趣的現象或是性質能夠被觀測到,并且務必確保此現象出現的可重復性。
9)數據拿到手后,很容易通過一些可視化軟件得到軌跡動畫,但這并不能拿來發(fā)文章。真正的工作才剛剛開始——分析數據,你所感興趣的現象或性質只是表面,隱含在它們之中的機理才是文章中的主題。