物理問(wèn)題的計(jì)算機(jī)模擬方法(1)-分子動(dòng)力學(xué)

物理問(wèn)題的計(jì)算機(jī)模擬方法(1)—分子動(dòng)力學(xué)編制僅供參考審核批準(zhǔn)生效日期地址：電話：傳真:郵編:碩士研究生課程《物理問(wèn)題的計(jì)算機(jī)模擬方法》講義適用專(zhuān)業(yè):凝聚態(tài)物理、材料物理與化學(xué)、理論物理、光學(xué)工程學(xué)時(shí)：30—40學(xué)時(shí)參考教材：1．[德]D.W.Heermann著，秦克誠(chéng)譯，理論物理中的計(jì)算機(jī)模擬方法，北京大學(xué)出版社，1996。2．[荷]Frenkel&Smit著，汪文川等譯，分子模擬—從算法到應(yīng)用，化學(xué)工業(yè)出版社，2002。3．M.P.AllenandD.J.Tildesley,ComputerSimulationofLiquids,ClarendonPress,Oxford,1989.4.A.R.Leach,MolecularModelling:PrinciplesandApplications,AddisonWesleyLongman,England,1996.5.[德]D.羅伯著，計(jì)算材料學(xué)，化學(xué)工業(yè)出版社，2002。6.[英]B.Chopard&MichelDroz著，物理系統(tǒng)的元胞自動(dòng)機(jī)模擬，祝玉學(xué)，趙學(xué)龍譯，清華大學(xué)出版社，2003。

目錄計(jì)算機(jī)模擬方法概論1.1序言1.2熱力學(xué)系統(tǒng)物理量的統(tǒng)計(jì)平均1.3分子動(dòng)力學(xué)方法模擬的基本思想1.4蒙特卡羅方法模擬的基本思想1.5元胞自動(dòng)機(jī)模擬的基本思想1.5.1簡(jiǎn)要的發(fā)展歷程1.5.2簡(jiǎn)單元胞自動(dòng)機(jī)：奇偶規(guī)則1.5.3元胞自動(dòng)機(jī)的一般定義確定性模擬方法—分子動(dòng)力學(xué)方法（MD）2.1分子動(dòng)力學(xué)方法2.2微正則系綜分子動(dòng)力學(xué)方法2.3正則系綜分子動(dòng)力學(xué)方法2.4等溫等壓系綜分子動(dòng)力學(xué)方法隨機(jī)性模擬方法—蒙特卡羅方法（MC）3.1預(yù)備知識(shí)3.2布朗動(dòng)力學(xué)（BD）3.3蒙特卡羅方法3.4微正則系綜蒙特卡羅方法3.5正則系綜蒙特卡羅方法3.6等溫等壓系綜蒙特卡羅方法3.7巨正則系綜蒙特卡羅方法離散性模擬方法—原胞自動(dòng)機(jī)（CA）4.1引言4.2元胞自動(dòng)機(jī)模擬*4.3元胞自動(dòng)機(jī)模擬的應(yīng)用

第一章計(jì)算機(jī)模擬方法概論§1.1序言什么是計(jì)算機(jī)模擬？SimulationModelling2．為什么要進(jìn)行計(jì)算機(jī)模擬？3．常用的計(jì)算機(jī)模擬方法確定性模擬方法：MD模擬MolecularDynamics隨機(jī)性模擬方法：MC模擬MonteCarlo離散性模擬方法：CA模擬CellularAutomata§1.2熱力學(xué)系統(tǒng)物理量的統(tǒng)計(jì)平均描述系統(tǒng)的坐標(biāo)（自由度）：x(t)={x1(t),x2(t),…xN(t)}系統(tǒng)的物理量：A(x(t))1．時(shí)間平均←分子動(dòng)力學(xué)（MD）模擬(1-1)2．系綜平均←蒙特卡羅(MC)模擬(1-2)—分布函數(shù)（幾率密度函數(shù)）(1-3)—配分函數(shù)(1-4)Ω—相空間H(x)—系統(tǒng)的哈密頓函數(shù)對(duì)于處于平衡態(tài)的系統(tǒng)，可以證明：對(duì)于實(shí)際的有限時(shí)間內(nèi)的平均，則有實(shí)際模擬的系統(tǒng)大小也是有限的：有限的粒子數(shù)N或有限的系統(tǒng)限度L對(duì)統(tǒng)計(jì)平均結(jié)果有影響?！?.3分子動(dòng)力學(xué)(MD)方法模擬的基本思想基本原理系統(tǒng)：N個(gè)粒子，體積V，粒子質(zhì)量為m描述一個(gè)粒子運(yùn)動(dòng)狀態(tài)的自由度：（ri,pi）（pi=mvi）相空間：6N維，相空間中的一點(diǎn)的坐標(biāo)XN=[rN,(mvN)]rN=(r1,r2,…,rN),vN=(v1,v2,…,vN)粒子間的相互作用勢(shì)：U(rN)=U(r1,r2,…,rN)=決定系統(tǒng)相軌跡XN（t）的運(yùn)動(dòng)方程：(1-5)物理量A的宏觀值，由A(XN)的時(shí)間平均獲得，即（離散情況：）對(duì)于平衡態(tài)：實(shí)際模擬時(shí)間總是有限的，模擬時(shí)間的長(zhǎng)短可通過(guò)判斷時(shí)間的增加對(duì)平均值的影響來(lái)確定，當(dāng)繼續(xù)增加時(shí)間帶來(lái)的平均值得變化在允許的誤差范圍之內(nèi)時(shí)，即可認(rèn)為模擬足夠長(zhǎng)了。計(jì)算步驟運(yùn)動(dòng)方程：即(1-6)或(1-7)數(shù)值求解：用差分近似表示微分（采用不同的差分格式，可得到不同的算法）。用顯示中心差分格式，將（7）式寫(xiě)為(1-8)由（7）和（8）式可得：(1-9)第一步：由（9）式計(jì)算第i個(gè)粒子在t+Δt時(shí)刻的位置坐標(biāo)。要啟動(dòng)計(jì)算，我們必須要知道最初兩點(diǎn)ri(0)和ri(Δt)第二步：對(duì)不同時(shí)刻t=Δt,2Δt,3Δt,……,LΔt(t0=0)計(jì)算物理量A(r1(lΔt),r2(lΔt),……,rN(lΔt))(l=1,2,……,N)第三步：計(jì)算物理量A的平均值L的大小由繼續(xù)增大L而<A>不變（或變化在誤差范圍內(nèi)）來(lái)確定?！?.4蒙特卡羅（MC）方法模擬的基本思想基本原理以正則系綜（T,V,N）為例正則分布：正則配分函數(shù)：系統(tǒng)能量：物理量：A(rN)=A(r1,r1,…,rN)系綜平均： （1-10）（位形積分）（1-11）用MC方法計(jì)算上述多維積分。計(jì)算步驟劃分原胞N個(gè)粒子—3N個(gè)空間自由度，3N維空間劃分成s個(gè)相等的原胞（s>>1）注意：由于積分中不含動(dòng)量，所以我們只需要在位置空間積分，而不需要在相空間中積分。當(dāng)系統(tǒng)的代表點(diǎn)落入第i個(gè)原胞時(shí)，則認(rèn)為系統(tǒng)處在狀態(tài)i，因此，s為系統(tǒng)可能的微觀狀態(tài)數(shù)目。于是，積分（10）和（11）可近似表述為(1-12)(1-13)建立馬爾可夫(Maρkoв)過(guò)程（鏈）將s個(gè)狀態(tài)看作一組隨機(jī)事件馬爾可夫鏈：從狀態(tài)i→j狀態(tài)j(i→j)的概率pij，只與i和j有關(guān)。，i=1,2,…,s若i經(jīng)歷n步到達(dá)j，其概率表示為，存在極限概率(j=1,2,…,s)uj為系統(tǒng)處在狀態(tài)j的概率。于是，沿?zé)o限長(zhǎng)的馬爾可夫鏈，物理量A的平均值可寫(xiě)為(1-14)選取，則（14）式為A的正則系綜平均值。（3）抽樣方法采用怎樣的抽樣方法所構(gòu)成的馬爾可夫鏈能得到上述平均值？粒子位置坐標(biāo)：粒子編號(hào)：=1,2,…N坐標(biāo)的三個(gè)方向：=1,2,3系統(tǒng)狀態(tài)：i=1,2,…,s給定粒子位置坐標(biāo)的變化量(小于系統(tǒng)體積的限度)給定系統(tǒng)的初態(tài)i，隨機(jī)選定4個(gè)隨機(jī)數(shù)，其中三個(gè)（=1，2，3），且–11，一個(gè)表示粒子編號(hào)=1,2,…N，由此隨機(jī)確定粒子位置的變化：（確保）若，則運(yùn)動(dòng)到新的位置，即系統(tǒng)由狀態(tài)i過(guò)渡到狀態(tài)j；若，則再選一個(gè)隨機(jī)數(shù)4（041），若，則粒子保留在原位置，不發(fā)生ij的躍遷；若，則發(fā)生ij的躍遷。由此進(jìn)行下去，則形成一個(gè)馬爾可夫鏈（或過(guò)程），此鏈的長(zhǎng)度L(即粒子行走的步數(shù)，遠(yuǎn)大于s)，由所計(jì)算的物理量的平均值(1-15)不再隨鏈的加長(zhǎng)而改變來(lái)確定。由此得到的平均值即正則平均值。一般來(lái)說(shuō)，L與N,V,T有關(guān)，比如，N=32~108,L=3000~5000。歸納起來(lái)，計(jì)算系統(tǒng)物理量的正則系綜平均值的具體步驟如下：第一步：給定系統(tǒng)的初始狀態(tài)(粒子的初始位置)ri和每一步的改變量；第二步：選擇四個(gè)隨機(jī)數(shù)，其中一個(gè)代表粒子的編號(hào)i(1iN)；另外三個(gè)表示粒子空間坐標(biāo)的改變x,y,z(-,=1,2,3);第三步：計(jì)算粒子i的新位置第四步：計(jì)算粒子在新舊兩個(gè)位置系統(tǒng)的能量之差第五步：由的大小判斷粒子i是否從ri運(yùn)動(dòng)到:若，則ri；若，則再選一個(gè)隨機(jī)數(shù)R(0<R<1),如果，則ri；如果，則ri不變，返回第二步。第六步：計(jì)算第七步：重復(fù)上述各個(gè)步驟，直到完成L步為止，最后利用公式（15）計(jì)算A的平均值。粒子間相互作用勢(shì)模型的選取最簡(jiǎn)單的兩種模型：（1）硬球模型(為硬球的直徑)（2）L—J勢(shì)§1.5元胞自動(dòng)機(jī)(CA)模擬的基本思想元胞自動(dòng)機(jī)：時(shí)間和空間都離散、物理參量只取有限數(shù)值集的物理系統(tǒng)的理想化模型cellularautomata或cellularautomaton—CA1.5.1簡(jiǎn)要的發(fā)展歷程1．自繁殖系統(tǒng)20世紀(jì)40年代，VonNeumann,構(gòu)造能解決非常復(fù)雜問(wèn)題的計(jì)算機(jī)，設(shè)想模仿人腦的行為——尋求與生物過(guò)程無(wú)關(guān)的情況下自繁殖機(jī)理的邏輯抽象。根據(jù)S.Ulam的建議，VonNeumann在由元胞構(gòu)成的完全離域的框架下處理這個(gè)問(wèn)題，構(gòu)造了一個(gè)完全離散的動(dòng)力學(xué)系統(tǒng)——元胞自動(dòng)機(jī)。第一個(gè)自復(fù)制元胞自動(dòng)機(jī)——由二維方形網(wǎng)絡(luò)組成，有數(shù)千個(gè)基本元胞構(gòu)成的自繁殖結(jié)構(gòu)。（1）一般機(jī)器只能構(gòu)造比自己簡(jiǎn)單的客體，而采用自復(fù)制元胞機(jī)，可獲得一種能產(chǎn)生新的、具有同樣復(fù)雜性和功能的“機(jī)器”；（2）VonNeumann的元胞自動(dòng)機(jī)規(guī)則具有所謂通用計(jì)算的性質(zhì)，這意味著，存在一種元胞自動(dòng)機(jī)的初始構(gòu)型，該元胞自動(dòng)機(jī)能產(chǎn)生任何計(jì)算機(jī)算法的解。通用計(jì)算的性質(zhì)指：用元胞自動(dòng)機(jī)演化規(guī)則能夠模擬任何計(jì)算機(jī)流程（邏輯選擇器開(kāi)關(guān)）。2．生命游戲機(jī)1970年，數(shù)學(xué)家JohnConway生命游戲機(jī)的概念，尋找能導(dǎo)致復(fù)雜行為的簡(jiǎn)單規(guī)則。設(shè)想一個(gè)類(lèi)似于棋盤(pán)的二維方形網(wǎng)格，每個(gè)元胞可能的狀態(tài)是活（狀態(tài)1）或死（狀態(tài)0），其更新規(guī)則是：有三個(gè)活元胞包圍的一個(gè)死元胞恢復(fù)為活元胞；由兩個(gè)以下或三個(gè)以上活元胞包圍的活元胞因孤立或擁擠而死亡。結(jié)果表明，生命游戲機(jī)有出乎意料的豐富行為，從原“湯”中顯示出來(lái)的復(fù)雜結(jié)構(gòu)，演變發(fā)展成為某些特殊的技藝，例如，可能形成所謂的滑翔機(jī)——緊鄰元胞的特殊排列，這些元胞具有沿直線彈道穿越空間運(yùn)動(dòng)的特性。生命游戲機(jī)也是具有計(jì)算通用性的元胞自動(dòng)機(jī)。3．模擬物理系統(tǒng)（1）20世紀(jì)70年代，Hardy,Pomeau和Pazzis建立了所謂的HPP格子氣體模型，用以在質(zhì)量和動(dòng)量守恒的情況下在方形網(wǎng)格上模擬粒子的碰撞行為。（2）1986年，F(xiàn)risch,Hasslacher和Pomeau提出了著名的FHP模型，這是在六邊形網(wǎng)格上模擬二維流體動(dòng)力學(xué)的第一個(gè)嚴(yán)格模型——全離散計(jì)算機(jī)模型替代風(fēng)洞試驗(yàn)。HPP和FHP通常稱(chēng)之為格子氣自動(dòng)機(jī)（LGA—LatticeGasAutomata）（3）Ising自旋動(dòng)力學(xué)模型，20世紀(jì)80年代末，Vichniac提出Q2R規(guī)則。（4）格子Boltzmann方法與多粒子模型格子Boltzmann方法或模型（LBM）：網(wǎng)格上定義的物理模型，在這個(gè)網(wǎng)格上與每個(gè)格位相關(guān)聯(lián)的變量是平均粒子數(shù)，或具有一定速度粒子出現(xiàn)的概率。該模型可以用均化或因子分解方法由元胞自動(dòng)機(jī)動(dòng)力學(xué)推導(dǎo)出來(lái)，或者自定義，而與特定的實(shí)現(xiàn)無(wú)關(guān)。格子Boltzmann模型保持了元胞自動(dòng)機(jī)方法的微觀水平解釋?zhuān)雎粤硕囿w的相關(guān)函數(shù)，但這種方法已經(jīng)成為目前模擬物理系統(tǒng)中最有前途的方法之一。在嚴(yán)格的元胞自動(dòng)機(jī)方法與較靈活的格子Boltzmann模型之間，有一種目前處于發(fā)展中的模型——多粒子模型。這種模型保留量化狀態(tài)的概念，但接受無(wú)限數(shù)值集，因此既保證了數(shù)值穩(wěn)定性（與LBM相反），又考慮了多體相關(guān)性。在模擬物理系統(tǒng)時(shí)，大量的可能狀態(tài)提供更多的靈活性，并產(chǎn)生小的統(tǒng)計(jì)噪聲。但多粒子動(dòng)力學(xué)更難設(shè)計(jì)，且在數(shù)值計(jì)算上比格子Boltzmann方法更慢。1.5.2簡(jiǎn)單元胞自動(dòng)機(jī)：奇偶規(guī)則簡(jiǎn)單元胞自動(dòng)機(jī)的演化規(guī)則（20世紀(jì)70年代，EdwardFredkin提出，定義在二維方形網(wǎng)上）：網(wǎng)格的每一個(gè)格位是一個(gè)元胞，以其位置r=(i,j)來(lái)標(biāo)記，其中i和j為行和列的標(biāo)號(hào)。函數(shù)描述每個(gè)元胞在時(shí)間t的狀態(tài)，其值為0或1。從時(shí)間t=0的初始條件及網(wǎng)格上給定的構(gòu)形值開(kāi)始，t=1時(shí)狀態(tài)按下列步驟求得：（1）對(duì)于每個(gè)格位r，都計(jì)算出其位于東、南、西、北4個(gè)最近鄰格位的值之和。應(yīng)使系統(tǒng)在i和j兩個(gè)方向循環(huán)（如同在環(huán)面上），從而確定出所有格位的計(jì)算值。（2）如果這個(gè)和值為偶數(shù)，則新?tīng)顟B(tài)為0（白色），否則為1（黑色）。重復(fù)上述步驟，得到t=2,3,4,…的狀態(tài)。這個(gè)元胞自動(dòng)機(jī)的奇偶規(guī)則可表示為：式中，符號(hào)代表“異或”邏輯運(yùn)算，也即模2和：11=00=0，10=01=1。反復(fù)迭代這個(gè)規(guī)則時(shí)，可得到非常精美的幾何圖形。1.5.3元胞自動(dòng)機(jī)的一般定義定義：規(guī)整的元胞網(wǎng)格覆蓋d維空間的一部分；歸屬于網(wǎng)格的每個(gè)格位r的一組布爾變量(r,t)={1(r,t),2(r,t),…,m(r,t)}給出每個(gè)元胞在時(shí)間t=0,1,2,…的局部狀態(tài)；演化規(guī)則R={R1,R2,…,Rm}按下列方式指定狀態(tài)(r,t)的時(shí)間演化過(guò)程：j(r,t+1)=Rj[(r,t),(r+1,t),(r+2,t),…,(r+q,t)]式中r+q指定從屬于元胞r的給定鄰居。鄰居：二維元胞自動(dòng)機(jī)的兩種鄰居：（1）VonNeumann鄰居，有一個(gè)中心元胞（要演化的元胞）和四個(gè)位于其近鄰東西南北方位的元胞組成；（2）Moore鄰居，除了前面涉及的最近鄰元胞外，還包括次近鄰的4個(gè)元胞，共九個(gè)元胞。還有一種有用的鄰居稱(chēng)為Margolus鄰居，將空間劃分成22元胞的鄰接單元塊，這個(gè)規(guī)則對(duì)位于單元塊內(nèi)的位置即左上、右上和左下、右下很敏感。三種鄰居如下圖所示。邊界條件：周期性邊界條件；固定邊界條件；絕熱邊界條件；映射邊界條件。備注：確定型元胞自動(dòng)機(jī)：演化規(guī)則確定，給定的初始狀態(tài)將始終演化出同樣的式樣。概率型自動(dòng)元胞機(jī)：演化規(guī)則包含一定的隨機(jī)性，給定的初始狀態(tài)可能演化出不同的式樣。

第二章確定性模擬方法—分子動(dòng)力學(xué)方法（MD）§2.1分子動(dòng)力學(xué)方法用分子動(dòng)力學(xué)方法模擬氣體、液體或固體系統(tǒng)1．分子動(dòng)力學(xué)的描述形式：哈密頓描述拉格朗日描述牛頓運(yùn)動(dòng)方程描述2．劃分MD元胞選取適當(dāng)大小的V=L3，太大耗費(fèi)計(jì)算時(shí)間，太小不能準(zhǔn)確反映系統(tǒng)的性質(zhì)。目的：維持系統(tǒng)恒定的密度。對(duì)平衡態(tài)系統(tǒng)，液體和氣體原胞的形狀無(wú)關(guān)緊要，但對(duì)晶體則不然。周期性邊界條件目的：減少表面效應(yīng)相互作用勢(shì)、最小影像約定和切斷距離n=(n1,n2,n3)元胞內(nèi)的粒子間元胞內(nèi)的粒子與影像的相互作用粒子的相互作用最小影像約定(minimumimageconvention):最小影像元胞中的一個(gè)粒子只與元胞中的另外N-1個(gè)粒子中的每一個(gè)或其最近鄰影像發(fā)生相互作用。切斷距離（cutoffdistance）：rcL/2（切斷半徑）上述約定相當(dāng)于用rc來(lái)切斷位勢(shì)，即rij>rc的相互作用忽略不計(jì)，代價(jià)是忽略了背景。積分格式于是，運(yùn)動(dòng)方程可寫(xiě)為(2-1)（1）遞推公式在第一章中求解運(yùn)動(dòng)方程時(shí)，我們是直接求解的關(guān)于粒子位置坐標(biāo)的二階微分方程，得到的遞推公式需要知道最初的兩點(diǎn)位置才能啟動(dòng)計(jì)算，但在實(shí)際計(jì)算中，我們常常是給出最初的位置和速度，于是，我們可通過(guò)選取一定的差分格式，有運(yùn)動(dòng)方程(16)得到關(guān)于位置和速度的遞推公式?；蛞阎f推公式的具體形式取決于差分各式的選取。（2）時(shí)間步長(zhǎng)選擇時(shí)間步長(zhǎng)的原則：在保證計(jì)算精度的前提下，盡量節(jié)省計(jì)算時(shí)間。實(shí)例：Ar原子系統(tǒng)，L—J勢(shì)時(shí)間步長(zhǎng)取為(=100ps=10fs)（3）約束條件保證能量、動(dòng)量或角動(dòng)量守恒。（4）減小計(jì)算誤差的技巧數(shù)值計(jì)算不可避免有誤差，與誤差有關(guān)的因素主要有：差分格式時(shí)間步長(zhǎng)切斷半徑最小影像約定等等6.計(jì)算熱力學(xué)量若系統(tǒng)的NVE恒定，則（微正則系綜平均—microcannonicalensembleaverage）若系統(tǒng)的NVT恒定，則（正則系綜平均—cannonicalensembleaverage）在實(shí)際計(jì)算中，往往還需要要求系統(tǒng)的總動(dòng)量P=0或恒定，因整個(gè)系統(tǒng)不受外界的作用。溫度與能量均分定理N—總粒子數(shù)，Nc—約束條件個(gè)數(shù),一般情況下，N>>NcP=0Nc=3位勢(shì)切斷誤差與修正g(r)—對(duì)關(guān)聯(lián)函數(shù)（paircorrelationfunction）—粒子數(shù)密度—當(dāng)原點(diǎn)位置上一個(gè)粒子時(shí)，在r附近dr內(nèi)發(fā)現(xiàn)有一個(gè)粒子的幾率。對(duì)氣體或液體，（各向同性）平均每個(gè)粒子的能量切斷誤差修正為：計(jì)算機(jī)模擬的組織初始化——給定粒子的初始位置、初始速度等；趨衡——從初始轉(zhuǎn)臺(tái)開(kāi)始，按運(yùn)動(dòng)方程要求的規(guī)律，從非平衡態(tài)趨向平衡態(tài)；投產(chǎn)——計(jì)算物理量的統(tǒng)計(jì)平均值?！?.2微正則系綜分子動(dòng)力學(xué)方法微正則系綜：NVE恒定，且P=0（分子動(dòng)力學(xué)方法的自然選擇）1．Verlet算法運(yùn)動(dòng)方程的顯示中心差分格式由可得到如下遞推公式（2-2）此即Verlet算法（A2）,其特點(diǎn)是：要啟動(dòng)此算法，必須已知兩點(diǎn)，所以又稱(chēng)為二步法；和不能同步算出，第n+1步算出的第n步的速度。2．Verlet算法的速度形式由此可得如下遞推公式：（2-3）此即Verlet算法的速度形式（A3）,其特點(diǎn)是：?jiǎn)?dòng)此算法，必須已知兩點(diǎn)；（2）和同步計(jì)算。（的計(jì)算需要知道）此算法比A2算法更穩(wěn)定。3．趨恒階段的能量調(diào)整由于系統(tǒng)初始狀態(tài)的能量離我們要模擬系統(tǒng)的能量（或溫度）有一定差異，這就需要在系統(tǒng)趨于平衡的階段進(jìn)行調(diào)整。最常用的方法之一就是進(jìn)行速度標(biāo)度：標(biāo)度因子：（2-4）Tref為設(shè)定的系統(tǒng)的溫度值。能量設(shè)定值：能量參考值：要注意的問(wèn)題：一般不需要每一步都進(jìn)行標(biāo)度，可每隔若干步（比如50步）調(diào)整一次；當(dāng)能量達(dá)到所設(shè)定的值后，停止標(biāo)度，在以下的模擬時(shí)間內(nèi)能量保持不變。實(shí)例氬原子系統(tǒng)，N=256

L作用力：取時(shí)間單位：取長(zhǎng)度單位：=0.3405nm取質(zhì)量單位：m=6.6338210-26kg使用上述單位，可得到約化單位下的作用力表達(dá)式：標(biāo)度因子：約化溫度：取時(shí)間步長(zhǎng)h=210-14s=20fs,相當(dāng)于約化時(shí)間0.064約化溫度：T*=2.53303K,T*=0.72286.5K約化數(shù)密度：*=N/L3,N=256*=0.636L=7.83,*=0.83134當(dāng)N=64時(shí)，*=0.83134L=4.25,所以，取rc=2.5是錯(cuò)的，因?yàn)橐笥?jì)算源程序見(jiàn)p129,程序PL1§2.3正則系綜分子動(dòng)力學(xué)方法正則系綜：NVT恒定，且P=0如何保持溫度T恒定？1．速度標(biāo)度方法算法A4:規(guī)定初始位置和初始速度；計(jì)算；計(jì)算；對(duì)速度進(jìn)行標(biāo)度：，返回（2）。說(shuō)明：在前面的微正則系綜的模擬中也對(duì)速度進(jìn)行了標(biāo)度，但其目的是使總能量達(dá)到所需要得值，而且不必每步進(jìn)行速度標(biāo)度，可隔若干步標(biāo)度一次，一旦總能量達(dá)到所設(shè)定的值，即可停止標(biāo)度；而在正則系綜的模擬中，每步都必須進(jìn)行標(biāo)度，以始終保持溫度恒定。通過(guò)一個(gè)廣義位勢(shì)引進(jìn)能量漲落，連同細(xì)致耦合的一種特殊選擇，會(huì)導(dǎo)致速度標(biāo)度機(jī)制。（見(jiàn)書(shū)中的證明，公式3.59有誤）由于控制溫度的反饋環(huán)內(nèi)的時(shí)間延遲，溫度將在一定程度上漲落。2．隨機(jī)方法（Anderson熱浴）體系與一強(qiáng)加了溫度的熱浴相耦合，與熱浴的耦合由偶爾作用于隨