第六章分子動力學(xué)模擬課件

第六章分子動力學(xué)模擬課件一、系綜理論本章主要內(nèi)容二、分子動力學(xué)方法分子動力學(xué)模擬三、模擬細節(jié)四、參量的計算五、液態(tài)水的MD模擬六、誤差分析七、分子動力學(xué)模擬方法的應(yīng)用

分子動力學(xué)模擬（moleculardynamicssimulation，簡稱MD）方法首先是由Alder和Wainwright提出的，現(xiàn)已逐漸成為預(yù)測系統(tǒng)特性、驗證理論和改進模型的計算工具。

MD方法的基本思想是把物質(zhì)看成由原子和分子組成的粒子系統(tǒng)，從該體系的某一假定的位能模型出發(fā)，并假定體系粒子的運動遵循經(jīng)典力學(xué)或量子力學(xué)描述的規(guī)律，若已知粒子的所有受力作用，則可以求解出運動方程而得到系統(tǒng)中全體粒子在相空間中的軌道，然后統(tǒng)計得到系統(tǒng)的熱力學(xué)參數(shù)、結(jié)構(gòu)和輸運特性等。也就是由體系的微觀性質(zhì)來求算其宏觀性質(zhì)。屬于微觀尺度的模擬技術(shù)。一、系綜理論分子動力學(xué)模擬

相空間或稱相宇，是經(jīng)典的系綜理論中為了描述系統(tǒng)的微觀狀態(tài)而引入的一個基本概念。1.1Phasespace

空間中的每一點，代表系統(tǒng)一種可能的運動狀態(tài)，系統(tǒng)的不同運動狀態(tài)，可以用相空間中不同的點來代表。從經(jīng)典力學(xué)可知，系統(tǒng)的運動狀態(tài)隨時間變化，因此相空間中代表系統(tǒng)運動狀態(tài)的點將在相空間中移動，形成軌道。分子動力學(xué)模擬

處在一定宏觀條件下，大量性質(zhì)和結(jié)構(gòu)完全相同的系統(tǒng)的集合。又稱統(tǒng)計系綜。系綜是采用統(tǒng)計方法描述系統(tǒng)的統(tǒng)計規(guī)律性時引入的一個最基本的概念。1.2系綜(Ensemble)

相空間中的不同點代表了系統(tǒng)的不同運動狀態(tài)，也可以把它們看作是性質(zhì)和結(jié)構(gòu)相同而初始條件不同的系統(tǒng)所處的狀態(tài)。這些系統(tǒng)的集合就是系綜。統(tǒng)計系綜的各個系統(tǒng)反映了實際系統(tǒng)在不同時刻的面貌。分子動力學(xué)模擬等溫等壓系綜（NPT）3．水的擴散系數(shù)和粘滯系數(shù)如果單時間步長選擇太短，小于體系分子位置的演變，也會造成誤差的增加，所以模擬時要結(jié)合所選取的求解分子運動方程的算法，合理選擇步長以減小誤差。TheinitialdistributionoftheMoleculardynamicssimulationisgeneratedinarandomdistribution.表3水的擴散系數(shù)模擬結(jié)果與文獻值和實驗值的比較表1水的勢能模擬結(jié)果與文獻值和實驗值的比較統(tǒng)計系綜的各個系統(tǒng)反映了實際系統(tǒng)在不同時刻的面貌。③、Verlet速度相加形式微正則系綜（NVE）：它是孤立的、保守的系統(tǒng)的統(tǒng)計系綜。指平衡狀態(tài)下理想氣體分子速度分布的統(tǒng)計規(guī)律。水溶液中A(7)、C(11)、E(15)和H(21)分子在Fe和FeCO3表面的平衡吸附構(gòu)型系綜的方法可用于由彼此存在相互作用的大數(shù)量粒子所組成的系統(tǒng)，按實際系統(tǒng)所處的不同宏觀條件，應(yīng)當(dāng)采用不同的系綜（微正則系綜、正則系綜、巨正則系綜）及其對應(yīng)的不同形式的分布函數(shù)。模擬體系的內(nèi)能由兩部分組成2Periodicboundaryconditions3Calculationofinteractions1.3準(zhǔn)各態(tài)歷經(jīng)假說

一個力學(xué)體系在長時間的運動中，它的代表點可以無限接近相空間中的任何點，系綜平均等于長時間的時間平均。

系綜的方法可用于由彼此存在相互作用的大數(shù)量粒子所組成的系統(tǒng)，按實際系統(tǒng)所處的不同宏觀條件，應(yīng)當(dāng)采用不同的系綜（微正則系綜、正則系綜、巨正則系綜）及其對應(yīng)的不同形式的分布函數(shù)。

統(tǒng)計方法的特點是由微觀量求宏觀量。分子動力學(xué)模擬1.4常用到的系綜正則系綜（NVT）：一個粒子數(shù)為N,體積為V,溫度為T和總動量為守恒量的系綜,在這個系綜中系統(tǒng)的粒子數(shù)(N),體積(V)和溫度(T)都保持不變,并且總動量為零。微正則系綜（NVE）：它是孤立的、保守的系統(tǒng)的統(tǒng)計系綜。在這種系綜中，體系與外界不交換能量，體系的粒子數(shù)守恒，體系的體積也不發(fā)生變化。巨正則系綜：由與溫度恒定的大熱源和化學(xué)勢恒定的大粒子源接觸，具有確定體積的系統(tǒng)構(gòu)成的統(tǒng)計系綜。等溫等壓系綜（NPT）分子動力學(xué)模擬1.5統(tǒng)計系綜（Statisticalensembles）④、(T,PN)constantpressure（等溫等壓系綜）①、(E,V,N)microcanonical（微正則系綜）,②、(T,V,N)canonical,constantvolume（正則系綜）③、(T,V,)grandcanonical分子動力學(xué)模擬2.1Newtonianmechanics

IntheMDmethod,everynewdistributionisderivedfromthepreviousonebyusingtheinteractionsbetweentheparticles.Theinteractionsdependsonthepositionoftheparticles.InthatpotentialtheparticlefeelsaforceAccordingtoNewton’ssecondlaw二、分子動力學(xué)方法分子動力學(xué)模擬2.2Potentialenergyfunctions對一個由N個原子構(gòu)成的簡單系統(tǒng)，其勢能項由下式給出

式中右端第一項是外場（如電場、磁場、聲場等）對系統(tǒng)的作用；第二項是兩體勢即系統(tǒng)中每兩個粒子間的相互作用；第三項是三體勢，表示系統(tǒng)中每三個粒子間的相互作用……有效兩體勢它包含多體效應(yīng)，可很好地反映系統(tǒng)粒子間的相互作用。分子動力學(xué)模擬簡單粒子系統(tǒng)采用Lennard-Joans勢、硬球勢、軟球勢、方阱勢等模型復(fù)雜粒子系統(tǒng)

采用多中心的位置-位置（site-site）相互作用模型，但其中心間的相互作用仍然采用簡單的位能模型。帶電分子庫侖（Coulomb）相互作用。處于外場中靜電相互作用是長程相互作用分子動力學(xué)模擬1、Lennard-Joans勢下面僅對簡單系統(tǒng)的相互作用模型給予簡介分子動力學(xué)模擬2.硬球勢（Hard-Sphere）3.軟球勢（Soft-Sphere）通常，v是為整數(shù)的參數(shù)。4.方阱勢（Square-Well）分子動力學(xué)模擬2.3Calculationsofforce,velocity,position

TheinitialdistributionoftheMoleculardynamicssimulationisgeneratedinarandomdistribution.Insimulation:Eachparticleisalsoassignedaninitialvelocityvi

分子動力學(xué)模擬TheforcecausesanaccelerationWhichinturnmodifiestheinitialvelocityviasAndmodifiestheinitialpositionriasTheMDsimulationcandescribesystemsthatevolveintime.ThenewpositionsarederivedfromtheNewtonianlawofmotionsandthereforedeterministic.分子動力學(xué)模擬2.4Equationsofmotion

為了在計算機上解運動方程，必須為微分方程建立一個有限差分格式，從差分方程中再導(dǎo)出位置和速度的遞推關(guān)系式。這些算法是一步一步執(zhí)行的，先算t時刻的位置和速度，然后在此基礎(chǔ)上計算t+1時刻的位置和速度。

微分方程最為直接的離散化格式來自泰勒展開:1、有限差分方法-預(yù)測校正法分子動力學(xué)模擬2、有限差分方法-Verlet算法①、Verlet算法的一般形式為了用數(shù)值方法求解微分方程，采用有限差分方法離散化得到差分格式兩式相加得分子動力學(xué)模擬

上式提供了一個方法，從粒子在前兩步（t和t-h）時刻的位置以及t時刻的作用力來得到粒子在t+h時刻的位置。

速度對軌道計算沒有關(guān)系，但對估算動能及速度相關(guān)函數(shù)（用來研究粒子的輸運特性）非常有用。速度形式為分子動力學(xué)模擬②、Verlet蛙跳格式（leap-frog）Thevelocityvariessoonehastochooseareasonableaveragevaluetobeused.Thevelocityatthemiddleofthestepoughttobeagoodcompromise,分子動力學(xué)模擬③、Verlet速度相加形式分子動力學(xué)模擬分子動力學(xué)模擬3.1初始化位型

模擬開始時首先要初始化位型分布，既首先要給定微元中分子的初始位置和初始速度。①、簡單立方晶格②、體心立方晶格③、面心立方晶格abc三、模擬細節(jié)分子動力學(xué)模擬abcabcabcabcabcabcabcabcabc分子動力學(xué)模擬初始速度模擬時，各粒子的初始速度按麥克斯韋速度分布取樣。Maxwell’sdistributionlawofvelocity

指平衡狀態(tài)下理想氣體分子速度分布的統(tǒng)計規(guī)律。在平衡狀態(tài)下，分布在任一速率區(qū)間內(nèi)的分子數(shù)與總分子數(shù)的比率為分子動力學(xué)模擬3.2Periodicboundaryconditions

通常選取小數(shù)體系（幾十個到數(shù)千個分子）作為研究對象，但由于位于表面和內(nèi)部的分子受力差別較大，將會產(chǎn)生表面效應(yīng)。

為了消除此效應(yīng)，同時建造出一個準(zhǔn)無窮大體積，使其可以更精確地代表宏觀系統(tǒng)，必須引入周期性邊界條件。

把小數(shù)體系看作一個中心元胞，此中心元胞被與中心元胞相同的其它元胞所包圍，當(dāng)某分子離開中心元胞時，該分子將在整個點格上移動以致它從中心元胞對面重新進入這個元胞。周期性邊界條件：分子動力學(xué)模擬D基本元胞ABFEHGC分子動力學(xué)模擬3.3Calculationof

interactionsi,j兩個離子的相互作用勢體系總勢能可表示為：

在計算勢能時必須考慮相互作用的力程，對粒子數(shù)為N的模擬系統(tǒng)，原則上任何兩個粒子間都有相互作用，在計算體系勢能時須進行N（N-1）/2次運算。

為了提高計算效率，實際模擬過程中通常在一個方便的力程上將位勢截斷，用以減少計算勢能所消耗的時間。分子動力學(xué)模擬

對勢能的最大貢獻來自于粒子的近鄰區(qū)域，位勢截斷常用的方法是球形截斷，截斷半徑一般取2.5σ或3.6σ，對截斷距離之外分子間相互作用能按平均密度近似的方法進行校正。分子動力學(xué)模擬★Cellmethod★Neighborlistmethod1、近鄰表-Verlet近鄰表2、近鄰表-網(wǎng)格近鄰表分子動力學(xué)模擬Thevelocityvariessoonehastochooseareasonableaveragevaluetobeused.ThenewpositionsarederivedfromtheNewtonianlawofmotionsandthereforedeterministic.對勢能的最大貢獻來自于粒子的近鄰區(qū)域，位勢截斷常用的方法是球形截斷，截斷半徑一般取2.Ru-Al合金斷裂過程動態(tài)模擬所以實際模擬過程中，必須認真分析個方面的情況，在精度和誤差之間作出平衡，統(tǒng)籌兼顧，以達到最佳的模擬效果。3．水的擴散系數(shù)和粘滯系數(shù)如果單時間步長選擇太短，小于體系分子位置的演變，也會造成誤差的增加，所以模擬時要結(jié)合所選取的求解分子運動方程的算法，合理選擇步長以減小誤差。式中右端第一項是外場（如電場、磁場、聲場等）對系統(tǒng)的作用；模擬體系分子數(shù)擴大一倍，總模擬時間并不僅增加一倍，而是呈級數(shù)形式上升。本文將水系統(tǒng)視為正則系綜進行MD模擬，模擬所用分子數(shù)為256，水分子初始位型按面心立方晶格分布，晶格所占空間體積由水的密度導(dǎo)出。Andmodifiestheinitialpositionrias模擬開始時首先要初始化位型分布，既首先要給定微元中分子的初始位置和初始速度。由于受計算機硬件的限制，在模擬過程中不可能對時間積分至無窮。同時受計算機硬件的限制，也不可能選擇更接近實際的復(fù)雜模型，如多體相互作用模型，從而導(dǎo)致誤差出現(xiàn)。3Calculationsofforce,velocity,position3.4標(biāo)度與趨衡

正則系綜的MD模擬，模擬過程中必須保證體系溫度恒定，所以在每一時間步要對分子的速度進行標(biāo)度使其滿足正則系綜的條件，標(biāo)度因子為，其中T是指定的體系溫度，T*是每一時間步的瞬時動能溫度。標(biāo)度：分子動力學(xué)模擬

一種按上述方法建立的系統(tǒng)不會具有所要的能量，而且很可能這個狀態(tài)不對應(yīng)一個平衡態(tài)，而所期望的體系參量是在平衡態(tài)通過統(tǒng)計平均得到的，所以為使系統(tǒng)達到平衡需要一個趨衡階段。趨衡：

對于正則系綜，體系趨衡是通過增減體系能量來實現(xiàn)的。如果在進行足夠多的運算步數(shù)以后，系統(tǒng)能持續(xù)給出確定的平均動能和平均勢能的數(shù)值，那么就認為平衡已經(jīng)建立。分子動力學(xué)模擬3.5TheMDunits

這里我們以能量（焦耳，J）和粒子間距離（米，m）為基本單位來對其它量進行約化，使其成為無量綱量。僅舉例如下：粒子數(shù)密度：溫度：能量：壓強：時間：力：分子動力學(xué)模擬Asimulationrunproducestherawdatainformofaverylargediskfile,Thefilecontainsthecompletestateofthesystemateachstep.Thediskfileisprocessedafterthesimulationisfinished.Itcontainsatleastallthepositionsandvelocitiesofallparticles.Thisinformationissufficienttocalculateallthepropertiesofthesystem.However,itismoreeconomicaltocalculatepropertiesduringthesimulationandstoretheminthefileratherthanreadingthefileandcalculatingthemafterwards.Thepotentialsandforcescertainlybelongtothisclassbecausetheyareneededateachtimestepanywayinordertocalculatethenewpositionandvelocitiesoftheparticles.四、參量的計算分子動力學(xué)模擬1．內(nèi)能

模擬體系的內(nèi)能由兩部分組成2．比熱

漲落是體系中熱力學(xué)量在不同時刻出現(xiàn)的波動。系統(tǒng)比熱可以通過能量的漲落導(dǎo)出，表示如下分子動力學(xué)模擬3.壓強

體系的壓強可通過維里定理導(dǎo)出：分子動力學(xué)模擬4．?dāng)U散系數(shù)

在模擬過程中，體系的遷移系數(shù)可以通過平衡相關(guān)函數(shù)導(dǎo)出。由Green-Kubo方法，計算擴散系數(shù)的公式如下：

由于受計算機硬件的限制，在模擬過程中不可能對時間積分至無窮。為了選擇合適的積分上限，在模擬過程中應(yīng)首先計算各水分子的速度自相關(guān)函數(shù)，從水體系歸一化的速度自相關(guān)函數(shù)可以看出在什么長度的時間間隔時，速度自相關(guān)函數(shù)逐漸趨零，也就是說在什么時間范圍內(nèi)速度基本上不再相關(guān)，那么求擴散系數(shù)所選取的積分區(qū)間就可以確定。為速度自相關(guān)函數(shù)（VACF）分子動力學(xué)模擬5、粘滯系數(shù)

利用平衡分子動力學(xué)模擬方法計算平衡系統(tǒng)的粘滯系數(shù)的基礎(chǔ)是Green-Kobo理論。根據(jù)該理論，粘滯系數(shù)由下式給定分子動力學(xué)模擬6、徑向分布函數(shù)------g(r)

徑向分布函數(shù)是與時間無關(guān)的關(guān)聯(lián)性的量度，精確地說就是與中心粒子間距為r~r+dr的體積元內(nèi)分子的平均數(shù)密度。在實際模擬過程中我們以元胞邊長的1/2為半徑作圓球，將其分割為厚度相等的100個球殼，計算分子間距R滿足（i<j，且i,j依次取值）的分子對個數(shù)N(i)，i=1,2,...100.分子動力學(xué)模擬★分子動力學(xué)模擬方法的基本步驟1.確定研究對象，如確定所研究的體系是屬于微正則系綜、正則系綜、巨正則系綜呢還是屬于等溫等壓系綜等等；2.建立適當(dāng)?shù)奈荒苣Ｐ?，如假定所研究的體系位勢形式滿足硬球勢、軟球勢或勒納德—瓊斯勢等等；3.初始化位形和速度；4.建立體系粒子的運動方程并進行求解，求解時需對位勢進行截斷，并應(yīng)用周期性邊界條件和最小影像數(shù)約定，每一時間步對速度進行標(biāo)度以使體系達到趨衡；5.通過統(tǒng)計平均來求解一些我們感興趣的熱力學(xué)參量等。五、液態(tài)水的MD模擬分子動力學(xué)模擬1．模擬參數(shù)

本文將水系統(tǒng)視為正則系綜進行MD模擬，模擬所用分子數(shù)為256，水分子初始位型按面心立方晶格分布，晶格所占空間體積由水的密度導(dǎo)出。各分子起始速度按Maxwell分布取樣。模擬溫度分別取=273K、293K、313K、333K、353K、373K。模擬過程中為保證體系溫度恒定，在每一時間步對分子的速度進行標(biāo)度；計算體系位能時采用球形截斷法，截斷半徑為2.5σ，對截斷距離之外分子間相互作用能按平均密度近似的方法進行校正；運動方程的求解采用VelocityVerletAlgorithm算法；時間步長為2×10-15s，每次模擬的總步數(shù)為105步，其中前50000步使體系趨衡，后50000步用來統(tǒng)計各種平衡性質(zhì)和結(jié)構(gòu)性質(zhì)；模擬過程中各參量均采用約化單位。分子動力學(xué)模擬

計算體系的內(nèi)能有著重要的意義，不僅是可以通過能量求出比熱、表面張力等熱力學(xué)量，更為重要的是內(nèi)能（或勢能）的演變是我們在模擬過程中判斷體系平衡的依據(jù)。體系是否達到平衡直接關(guān)系到熱力學(xué)量的真實性、穩(wěn)定性和準(zhǔn)確性。2．水的內(nèi)能和比熱分子動力學(xué)模擬表1水的勢能模擬結(jié)果與文獻值和實驗值的比較方法實驗值[34]本文MDMC[35]MD[24]勢能（kJ/mol）-44.77*-41.99-36.09-42.20模型×L-JTIPSPC系綜×NVTNVTEVN表2水的比熱的模擬結(jié)果和實驗結(jié)果的比較溫度（K）273293313333353373模擬值（J/gK）4.01773.97853.95923.94133.93293.9254實驗值（J/gK）4.20744.18014.17884.18554.1985×分子動力學(xué)模擬3．水的擴散系數(shù)和粘滯系數(shù)表3水的擴散系數(shù)模擬結(jié)果與文獻值和實驗值的比較方法實驗值[34]本文MDMD[35]MD[36]擴散系數(shù)(×10-5cm2/s)2.32.9622.952.8模型×L-JSPCMCY系綜×NVTNVTNVT分子動力學(xué)模擬4．水的徑向分布函數(shù)

徑向分布函數(shù)是對物質(zhì)結(jié)構(gòu)的反映，在研究物質(zhì)特別是液態(tài)時有著重要的作用。分子動力學(xué)模擬1、位能模型

位能模型是對體系中分子或離子之間相互作用勢的反映，位能模型選擇的正確與否決定了模擬的成敗。但是許多情況下由于研究對象的復(fù)雜性，使對問題的認識還不夠深入，導(dǎo)致所選擇的位能模型和實際情況存在一定差距；同時受計算機硬件的限制，也不可能選擇更接近實際的復(fù)雜模型，如多體相互作用模型，從而導(dǎo)致誤差出現(xiàn)。模型的名稱ST2TIPSCITIPS2SPCBFPPC模型的類型經(jīng)驗公式經(jīng)驗公式從頭計算經(jīng)驗公式經(jīng)驗公式半經(jīng)驗公式偶極子及交換作用表1關(guān)于純水結(jié)構(gòu)的理論模型六、誤差分析分子動力學(xué)模擬2、模擬分子數(shù)

分子動力學(xué)模擬的主要任務(wù)就是要從微觀來求算體系的宏觀性質(zhì)，所以模擬過程中參與模擬的分子數(shù)越多，模擬結(jié)果就越真實，越接近體系的宏觀屬性。例如Tanaka等人采用1.36×104個水分子對水體系進行了模擬，其擴散系數(shù)更加接近實驗值。但也應(yīng)考慮到：模擬過程中由于分子受力和勢能的計算需要涉及微元內(nèi)所有的分子對，所以計算量很大，非常消耗機時，一般占總計算量的80%左右。模擬體系分子數(shù)擴大一倍，總模擬時間并不僅增加一倍，而是呈級數(shù)形式上升。

要想解決大數(shù)粒子體系的模擬問題，還須引入一些模擬技巧，如VerletNeighbourList和LinkList等方法。分子動力學(xué)模擬粒子數(shù)本文模擬LCVerlettableLLC2560.68sec×××1372×0.45sec0.072sec0.037sec2048×0.57sec0.114sec0.058sec4000×××0.094sec6912×1.59sec0.529sec0.164sec表2不同方法每時間步所消耗CPU時間的比較分子動力學(xué)模擬3、求解運動方程的算法

求解分子運動方程的方法有兩種，一種是VerletAlgorithm，它是對二階微分運動方程直接求解的一種方法；另一種則是Gear預(yù)測—校正算法（GearPredictor-CorrectorAlgorithm），它是通過泰勒級數(shù)展開的方法實現(xiàn)的。

在時間步長較小的情況下，通過高階Gear算法得到的結(jié)果更為精確，但是模擬者對大步長更感興趣，因為相同的總模擬步數(shù)，單步長越大，模擬體系演化的時間越長，就越接近實際情況，模擬結(jié)果就越真實，在這種情況下VerletAlgorithm算法則更具吸引力。分子動力學(xué)模擬4、總時間步、時間步長

一般來說模擬總時間步越多，單步長越大，模擬結(jié)果就接近實際情況，但是由于受計算機硬件的限制，不可能對一個足夠長的時間過程進行模擬，同時單步長較大，對體系演變反映的跨度超長，計算誤差也就越大；如果單時間步長選擇太短，小于體系分子位置的演變，也會造成誤差的增加，所以模擬時要結(jié)合所選取的求解分子運動方程的算法，合理選擇步長以減小誤差。

求解中的另一個重要問題是時間步長的選擇：步長h連同實現(xiàn)的MD步數(shù)，決定了相空間有多大一部分被抽樣，看起來越大越好，這樣可以對更大的部分抽樣，但是由h所決定的時間標(biāo)尺和實際體系發(fā)生變化的時間標(biāo)尺并不一定完全相同，所以時間步長的選取應(yīng)根據(jù)研究對象的具體情況而定分子動力學(xué)模擬5、位能截斷

為提高計算效率，實際模擬過程應(yīng)進行位勢截斷，對截斷距離之外分子間相互作用能按平均密度近似的方法進行校正，雖然進行了校正，但難以避免誤差的出現(xiàn)，特別是在分子間存在長程相互作用時，進行位勢截斷會導(dǎo)致嚴重的失真，在這種情況下則必須采用Ewald求和的方法。

上述的誤差分析可以看出，誤差產(chǎn)生的原因具有多樣性。為了減少誤差、提高精度，單純從某一個方面入手是盲目和片面的，在誤差產(chǎn)生的過程中它們是相互影響、相互制約的。所以實際模擬過程中，必須認真分析個方面的情況，在精度和誤差之間作出平衡，統(tǒng)籌兼顧，以達到最佳的模擬效果。分子動力學(xué)模擬Adsorptionofmethaneonzeolite

分子篩是具有均勻的微孔、其孔徑與一般分子大小相當(dāng)為一類吸附劑或薄膜類物質(zhì)。在化學(xué)化工中經(jīng)常被用作催化劑或者催化劑的載體。七、分子動力學(xué)模擬方法的應(yīng)用分子動力學(xué)模擬Layer-CellModelereasilycreatessimulationofwater-benzeneinterface

在體系內(nèi)部物理性質(zhì)和化學(xué)性質(zhì)完全均一的一部分稱為“相”（Phase)。相與相之間在指定的條件下有明顯的界面，在界面上，從宏觀的角度看，性質(zhì)的改變是飛躍式的。分子動力學(xué)模擬MgOSubstrateSputteredbyArPlasma分子動力學(xué)模擬AqueousLiBrsubjecttoexternalelectricfield分子動力學(xué)模擬分子動力學(xué)模擬0ps200ps1000ps2250ps2760ps2910ps3500ps石墨烯/硅納米線復(fù)合材料自卷曲過程平行放置相互靠近完成第一層包覆剩余石墨烯完成包覆形成穩(wěn)定結(jié)構(gòu)分子動力學(xué)模擬

實驗測定與

MD

模擬結(jié)果的比較W.Linert,andF.Renz,J.Chem.Inf.Comput.Sci.,