概述
分子動力學法(MD,MolecularDynamicsmethod),是用計算機模擬的一種,是調查物質諸性質時候使用的手法之一。根據在計算機中每時每刻的追蹤全部的粒子的運動規律,導出物質全體的性質,這就是分子動力學法。
分子動力學法是要嚴格求解每個粒子的運動方程,通過分析系統來確定粒子的運動狀態。MD方法一般認為粒子服從牛頓運動規律,粒子所受到的作用力通過粒子間相互作用勢計算,因此粒子的初始條件和運動方程中的受力狀況一旦確定,系統就可被精確求解。
經典MD模擬,其系統規模在一般的計算機上也可達到數萬個原子,模擬時間為納秒量級。2006年進行了三千二百億個原子的模擬(IBMlueGene/L)。
分子動力學
分子動力學是一套分子模擬方法,該方法主要是依靠計算機來模擬分子、原子體系的運動,是一種多體模擬方法。通過對分子、原子在一定時間內運動狀態的模擬,從而以動態觀點考察系統隨時間演化的行為。通常,分子、原子的軌跡是通過數值求解牛頓運動方程得到,勢能(或其對笛卡爾坐標的一階偏導數,即力)通常可以由分子間相互作用勢能函式、分子力學力場、從頭算給出。對於考慮分子本身的量子效應的體系,往往採用波包近似處理或採用量子力學的費恩曼路徑積分表述方式處理。分子動力學也常常被採用作為研究複雜體系熱力學性質的採樣方法。以在由分子體系的不同狀態構成的系綜中抽取樣本,從而計算體系的構型積分,並以構型積分的結果為基礎進一步計算體系的熱力學量和其他巨觀性質。分子動力學最早在20世紀50年代由物理學家提出,如今廣泛套用與物理、化學、生物體系的理論研究中。
基本步驟
確定起始構型
進行分子動力學模擬的第一步是確定起始構型,一個能量較低的起始構型是進行分子模擬的基礎,一般分子的起始構型主要來自實驗數據或量子化學計算。
在確定起始構型之後要賦予構成分子的各個原子速度,這一速度是根據玻爾茲曼分布隨機生成的,由於速度的分布符合玻爾茲曼統計,因此在這個階段,體系的溫度是恆定的。另外,在隨機生成各個原子的運動速度之後須進行調整,使得體系總體在各個方向上的動量之和為零,即保證體系沒有平動位移。
進入平衡相
由上一步確定的分子組建平衡相,在構建平衡相的時候會對構型、溫度等參數加以監控。
進入生產相
進入生產相之後體系中的分子和分子中的原子開始根據初始速度運動,可以想像其間會發生吸引、排斥乃至碰撞,這時就根據牛頓力學和預先給定的粒子間相互作用勢來對各個粒子的運動軌跡進行計算,在這個過程中,體系總能量不變,但分子內部勢能和動能不斷相互轉化,從而體系的溫度也不斷變化,在整個過程中,體系會遍歷勢能面上的各個點(理論上,如果模擬時間無限)。計算分析所用樣本正是從這個過程中抽取的。
計算結果
用抽樣所得體系的各個狀態計算當時體系的勢能,進而計算構型積分。
作用勢與動力學計算
作用勢的選擇與動力學計算的關係極為密切,選擇不同的作用勢,體系的勢能面會有不同的形狀,動力學計算所得的分子運動和分子內部運動的軌跡也會不同,進而影響到抽樣的結果和抽樣結果的勢能計算,最初的分子動力學計算採用比較簡單的剛球勢,現在更多地採用蘭納-瓊斯勢,後者能夠更好的與粒子間相互作用擬合。
時間步長與約束動力學
分子動力學計算的基本思想是賦予分子體系初始運動狀態之後利用分子的自然運動在相空間中抽取樣本進行統計計算,時間步長就是抽樣的間隔,因而時間步長的選取對動力學模擬非常重要。太長的時間步長會造成分子間的激烈碰撞,體系數據溢出;太短的時間步長會降低模擬過程搜尋相空間的能力,因此一般選取的時間步長為體系各個自由度中最短運動周期的十分之一。
但是通常情況下,體系各自由度中運動周期最短的是各個化學鍵的振動,而這種運動對計算某些巨觀性質並不產生影響,因此就產生了禁止分子內部振動或其他無關運動的約束動力學,約束動力學可以有效地增長分子動力學模擬的時間步長,提高搜尋相空間的能力。
套用
分子動力學的計算過程給定了體系的總能量,因此適用於對微正則系綜的模擬計算,另外由於分子動力學計算過程始終是時間的函式,因此一些與時間有關的巨觀量如擴散係數的模擬必須套用分子動力學。
另外,在實際套用中,經常把分子動力學方法和蒙特·卡羅方法聯合使用。在近年,多尺度模擬計算已經得到了很多學者的關注。在多尺度模擬計算中,分子動力學方法研究納米級現象,蒙特·卡羅方法研究微觀形態,有限元方法套用於巨觀領域。通過多種尺度的多種模擬計算方法的聯合套用,令納觀與巨觀聯結起來。