2005年第29卷石油大學(xué)學(xué)報(自然科Vol. 29 No 4第4期Journal of the University of Petroleum, ChinaAug.2005文章編號:1000-5870(2005)04014304綜述t從頭計算分子動力學(xué)方法及其應(yīng)用形1藍(lán)建慧,盧貴武,黃喬松,李英峰,朱閣(中國石油大學(xué)物理科學(xué)與技術(shù)學(xué)院,山東東營257061)摘要:從頭計算分子動力學(xué)方法把密度泛函理論和分子動力學(xué)方法有機(jī)地結(jié)合起來,使電子的極化效應(yīng)及化學(xué)鍵的本質(zhì)均可用計算機(jī)分子模擬方法進(jìn)行研究,是目前計算機(jī)模擬實驗中最先進(jìn)、最重要的方法之一。文章簡述了從頭計算分子動力學(xué)方法的基本原理,介紹了該方法在水、水溶液及其他氫鍵液體的結(jié)構(gòu)與動力學(xué)研究中的應(yīng)用關(guān)鍵詞:從頭計算;密度泛函理論;分子動力學(xué);計算機(jī)分子模擬中圖分類號:O35文獻(xiàn)標(biāo)識碼:AMethod of ab initio molecular dynamics and its applicationsLAN Jian-hui, LU Gui-wu, HUANG Qiaosong, LI Ying-feng, ZHU GeCollege of Physics Science and Technology in China University of Petroleum, Dongying 257061, China)Abstract: The ab initio molecular dynamics method, which combines the density functional theory with the molecular dynamics methodology, made it convenient to study the electronic polarization effects and the nature of the chemical bonds interm of the computer molecular simulation. The method is one of the most important and advanced computer simulation experiment methods. The basic principle of the ab initio molecular dynamics method and its applications in structure and dynamics research of liquid water, aqueous solutions and other hydrogen-bond liquids were introducedKey words: ab initio density functional theory: molecular dynamics; computer molecular simulation現(xiàn)代凝聚態(tài)理論研究應(yīng)用最普遍的方法之一是同反應(yīng)設(shè)置不同的參數(shù),模擬結(jié)果仍可能偏離反應(yīng)分子動力學(xué)(MD)方法。MD方法把原子的運動與路徑{2。為了解決上述問題,從頭計算分子動力學(xué)特定的軌道聯(lián)系在一起,通過求解原子的牛頓運動方法( ab initio molecular dynamics,AIMD)被提出。方程得到體系的熱力學(xué)性質(zhì)。MD計算的核心是選本文中重點介紹從頭計算分子動力學(xué)方法的基本原擇合適的力場,即在系統(tǒng)中引入簡單數(shù)學(xué)模型來描理及其在水、水溶液和其他氫鍵液體研究中的應(yīng)用。述原子間的結(jié)合、彎曲和二面角勢及原子間的范德華力和靜電作用并預(yù)設(shè)實驗數(shù)據(jù)或模型參數(shù)進(jìn)行1從頭計算分子動力學(xué)方法計算。利用該方法既可得到原子的運動軌跡,也可從頭計算分子動力學(xué)(AIMD)方法主要基于以在計算過程中對平衡或非平衡系統(tǒng)的微觀細(xì)節(jié)進(jìn)行下3個假設(shè):(1)忽略系統(tǒng)的核量子效應(yīng);(2)認(rèn)為系研究,因此被廣泛應(yīng)用于研究各種物理、化學(xué)問題,統(tǒng)滿足軌道近似(即單電子近似);(3)認(rèn)為系統(tǒng)滿足在處理凝聚態(tài)體系時獲得了極大的成功。但是絕熱近似。其中電子基態(tài)本征函數(shù)和本征值的計算MD方法是基于力場的,因此存在嚴(yán)重的缺陷。首是AIMD的核心內(nèi)容。先它忽略了電子極化效應(yīng)。后來有人提出了校正該電子基態(tài)計算屬于復(fù)雜的量子多體問題,需引缺陷的極化模型,但極化模型只適用于特定的問題,人密度泛函理論( density functional theory,DFT)以目前還沒有得到廣泛的應(yīng)用。此外,MD方法無法簡化計算量,把復(fù)雜的多體問題轉(zhuǎn)化為一組自洽的描述化學(xué)鍵形成或斷裂的本質(zhì)問題,若采用經(jīng)驗價單鍵方法或其他的半經(jīng)驗方法作近似處理即使對不和原中國煤化T1方程并根據(jù)電子的影響程度,對交CNMHG收稿日期:2004-12-21基金項目:山東省自然科學(xué)基金資助項目(Y2003A01)和石油科技中青年創(chuàng)新基金(04E7038)作者簡介:藍(lán)建慧(1979-),女(漢族),山東即墨人,碩士研究生,專業(yè)方向為計算物理。144石油大學(xué)學(xué)報(自然科學(xué)版)2005年8月?lián)Q勢采用局域密度近似( local density approxima-后設(shè)定虛擬系統(tǒng)廣義拉格朗日-歐拉方程中的折合tion,LDA)或廣義梯度近似( general gradient ap-質(zhì)量等參量,使電子的運動標(biāo)度遠(yuǎn)小于離子的運動proximation,GGA),從而方程可解。這也是目前凝標(biāo)度,從而保證電子隨著離子位形的改變盡量趨于聚態(tài)物理中計算電子結(jié)構(gòu)普遍采用的方法。電子基基態(tài)(即電子在等能面E附近作簡諧振動),這樣無態(tài)計算不僅得到原子間相互作用勢,也為分子動力論高溫還是低溫條件下,整個耦合系統(tǒng)就處在BO學(xué)研究提供了精確的力場,這種將密度函數(shù)理論和( Born-Oppenheimer)面上,然后開始電子和離子耦分子動力學(xué)結(jié)合起來的方法,是目前從頭計算分子合系統(tǒng)的動力學(xué)模擬。根據(jù)第一過程的計算結(jié)果即動力學(xué)中普遍采用的算法之一(即 Car-Parrinello,可計算離子運動方程中的“虛擬力”,求解離子運動CP算法)2。該算法的基本步驟見圖1。由圖可方程并得到t=t+△t時刻的離子坐標(biāo)。見,AIM模擬可歸結(jié)為電子系統(tǒng)動力學(xué)計算、電子和離子耦合系統(tǒng)動力學(xué)計算兩個過程。2從頭計算分子動力學(xué)方法的應(yīng)用始坐標(biāo)從頭計算分子動力學(xué)方法在凝聚態(tài)體系的微觀YES幾何結(jié)構(gòu)和電子結(jié)構(gòu)以及反應(yīng)動力學(xué)的研究中應(yīng)用計算基態(tài)單電子波函數(shù)廣泛,本文中主要介紹該方法在液態(tài)水及其他氫鍵液體的結(jié)構(gòu)與動力學(xué)研究中的應(yīng)用。計算離子間相互作用勢V及電子系統(tǒng)總能量E2.1液態(tài)水的從頭計算分子動力學(xué)模擬計算離子虛擬力計算電子的虛擬力由于300K水的特殊性,對該溫度下水的結(jié)構(gòu)的AMD研究13-7比較多。文獻(xiàn)[8]中研究了由64求解離子運動方程求解電子運動方程個水分子構(gòu)成的小系統(tǒng),模擬布里淵區(qū)中心點采用 Troullier- Martins型原子贗勢,電子相關(guān)效應(yīng)新的離子坐標(biāo)新的體系波函數(shù)用BLYP方法處理。運行周期為11ps,時間步長為0.17fs,得到的徑向分布函數(shù)與實驗結(jié)果一致。文t=t+△t獻(xiàn)[7]采用 Perdew交換關(guān)聯(lián)函數(shù),當(dāng)運行周期75ps時,得到的O-O徑向分布函數(shù)與實驗吻合。圖1從頭計算分子動力學(xué)模擬流程AMD模擬還被用于計算水的紅外光譜和平均第一過程主要有以下幾個步驟:①根據(jù)電子系偶極矩。文獻(xiàn)[4]和[9]中分別對含32個水分子和統(tǒng)動力學(xué)計算原子初始構(gòu)型的基態(tài)單電子軌道波函64個水分子的體系進(jìn)行了研究采用BLYP函數(shù)且數(shù)。通常采用嘗試法,先用少量的平面波基函數(shù)將運行周期為10~12p時,水的紅外光譜和平均偶分子軌道展開,按矩陣對角化技術(shù)得到分子軌道本極矩的計算結(jié)果與實驗結(jié)果吻合較好。另外,將路征波函數(shù),再以按照擴(kuò)展平面波基函數(shù)展開的分子徑取樣法10和AMD技術(shù)結(jié)合起來,可以闡述水軌道作嘗試波函數(shù)進(jìn)行虛擬動力學(xué)模擬。②根據(jù)密的自身離解以及H3OOH離子對的形成問題。度泛函理論計算電子體系總能量E和離子間的相應(yīng)當(dāng)指出,由于AIMD模擬中系統(tǒng)和基組大小互作用勢v。體系的基態(tài)能量對應(yīng)于E的極小值,的影響3-以及直接從中子和x射線衍射數(shù)據(jù)中是電子基態(tài)密度的泛函。離子間的相互作用勢v提取這些量尚存在困難通過AIMm方法得到300(即電子和離子耦合系統(tǒng)的勢能)以及V對離子坐K水的準(zhǔn)確徑向分布函數(shù)仍是一個具有挑戰(zhàn)性的課標(biāo)的微分根據(jù)密度泛函理論和 Hellman- Feynman定題。理計算。③計算電子波函數(shù)正交化產(chǎn)生的“虛擬2.2液態(tài)氨的從頭計算分子動力學(xué)模擬力”,求解運動方程,得到新的電子基態(tài)波函數(shù)。如液氨化學(xué)類似于水化學(xué),水中的反應(yīng)類型在液果運動方程中未包含由波函數(shù)正交化產(chǎn)生的“虛擬氨中基本上都可以見到。液氨分子中氮原子的電負(fù)力”需要對新的波函數(shù)正交化。④將新的波函數(shù)作性比才V山中國煤化工建電子更容易被利為t=t+△時刻的輸人,進(jìn)行下一步的電子動力用故CNMH(有機(jī)反應(yīng)和與金學(xué)模擬直至t>tm(tm即AIMD模擬的最大設(shè)屬反應(yīng)的溶液中。其結(jié)構(gòu)因子和徑向分布函數(shù)已于定時間)。最近通過中子衍射技術(shù)在實驗中得到確認(rèn),通過第二過程算法步驟如下:得到初始本征波函數(shù)AIMD研究結(jié)果和純液態(tài)氨實驗數(shù)據(jù)的比較,可以第29卷第4期藍(lán)建慧等:從頭計算分子動力學(xué)方法及其應(yīng)用145對AIMD方法做出評價。為了減小輕氫原子核量子效應(yīng)的影響,用KOD代目前,基于CP運動方程的AMD模擬被應(yīng)用替KOH。交換關(guān)聯(lián)采用BLYP方法處理,核心電子于含32個氨分子、箱長為1127A且滿足周期性邊用 Troullier-Martins腰勢和一個鉀的半核贗勢處界條件的樣品,實驗溫度分別為260K和273理。模擬周期為10ps,采用600au虛擬電子群時K12.13]。交換關(guān)聯(lián)函數(shù)用 BLYP GGA方法4.1處間步長為5ps此模擬長度可容納2ps的平衡態(tài)。理,核心電子用 Troullier Martins型贗勢處理。系模擬計算了偶極矩并給出紅外光譜。由于模擬長度統(tǒng)平衡態(tài)可保持的最長時間為22p,運行周期為相對較短,采用最大熵技術(shù)4提取紅外光譜。比較60ps,時間步長為0.12fs。260K液氨的結(jié)構(gòu)因計算得到的KOD水溶液的紅外光譜和14moL子和徑向分布函數(shù)的計算結(jié)果與實驗結(jié)果1)吻合KOH溶液的實驗測量結(jié)果2發(fā)現(xiàn),前者再現(xiàn)了14極好。另外根據(jù)計算1確定的自擴(kuò)散常數(shù)為1.1 mol/L Koh溶液的實驗所得光譜的所有特征×10-4cm2·s-1,比較接近實驗值1.0×10-4cm23結(jié)束語2.3300K液態(tài)甲醇的從頭計算分子動力學(xué)模擬介紹了結(jié)合有限溫度分子動力學(xué)(MD)與虛擬另一種重要的氫鍵液體是甲醇(CH3OH)。和電子結(jié)構(gòu)計算的從頭計算分子動力學(xué)模擬方法,以液氨一樣,甲醇作為溶劑被廣泛地應(yīng)用于許多普通及液態(tài)水、液態(tài)氨、液態(tài)甲醇和KOD水溶液中的的有機(jī)反應(yīng)中,尤其在新興的燃料電池技術(shù)中發(fā)揮AM模擬應(yīng)用。此外,AMD方法在液體結(jié)構(gòu)、酸著重要作用,是一種重要的工業(yè)用液體。利用中子化學(xué)、工業(yè)和生物催化劑以及材料學(xué)等研究領(lǐng)域也衍射技術(shù),液體甲醇的結(jié)構(gòu)最近得到確認(rèn)18.9。與得到廣泛應(yīng)用揭示了許多經(jīng)驗?zāi)Ｐ蜔o法解釋的現(xiàn)從頭計算相比,利用該技術(shù)更容易得到結(jié)構(gòu)因子和象,同時啟發(fā)人們對一些實驗數(shù)據(jù)做出新的解釋,設(shè)徑向分布函數(shù)。 Morrone等201對液體甲醇中質(zhì)子計新的實驗研究方案的輸運特性進(jìn)行了AIMD研究,模擬采用含32個應(yīng)當(dāng)指出,從頭計算分子動力學(xué)方法使基于密甲醇分子、周期性箱的尺寸為12.93A的樣品,交換度泛函理論的第一性原理計算直接應(yīng)用于統(tǒng)計力學(xué)與關(guān)聯(lián)用 BIYP GGA方法處理。分別采用兩種不模擬成為可能極大地擴(kuò)展了計算機(jī)分子模擬的廣同的模擬方案進(jìn)行研究。第一方案:選用25個平面度和深度,成為計算機(jī)分子模擬最先進(jìn)和最重要的波并采用超軟贗勢法(21,運行周期為20ps;第二方方法之一。但與以力場為基礎(chǔ)的傳統(tǒng)計算機(jī)模擬方案:采用量子力學(xué)/分子力學(xué)(QMMM)方法2),法相比,AIMD方法的計算量顯著增加。傳統(tǒng)計算OH基在QM級別上用BLYP函數(shù)處理14.1,CH3機(jī)分子模擬可應(yīng)用于時間標(biāo)度在幾十個納秒數(shù)量基在MM級別上用 AMBER力場處理。研究結(jié)果級原子數(shù)約為104~106的大分子體系,AMD計發(fā)現(xiàn),對于結(jié)構(gòu)因子,QMMM模型超軟贗勢法與算則主要應(yīng)用于時間標(biāo)度在幾十個皮秒數(shù)量級而原實驗三者吻合較好,而對于OO,OH和CC的徑向子數(shù)在幾十到幾百的小分子體系。這些特點在實際分布函數(shù),QM^MM模型比超軟贗勢法好得多應(yīng)用中需格外注意。2.4KOD溶液的從頭計算分子動力學(xué)模擬了解水環(huán)境中由于H或OH的增加而產(chǎn)生的參考文獻(xiàn)電荷轉(zhuǎn)移行為在酸堿化學(xué)、生物工程燃料電池、皂[1]趙宇軍姜明,曹培林從頭計算分子動力學(xué)J物理化和工業(yè)催化劑等領(lǐng)域至關(guān)重要。含水合氫離子學(xué)進(jìn)展,1998,18(1):47H3O和氫氧根離子OH的水溶液是良好的質(zhì)子導(dǎo)ZHAO Yu-jun, JIANG Ming, CAO Pei-lin. Ab initio體,其輸運過程的微觀細(xì)節(jié)目前仍存在爭議。部分molecular dynamics[J]. Progress in Physics, 1998, 18爭議在于難以找到能證明其特殊機(jī)制的明確的實驗(1):47-75證據(jù)。高水平AMD模擬在揭開這些機(jī)制的細(xì)節(jié)2] MARK E Tuckerman. Ab initio molecular dynamics:tend novel applications [J].J方面可以發(fā)揮特殊作用。在文獻(xiàn)[23]的研究中,1.5中國煤化工5mL的溶液含32個水分子和一個孤立的KOD分[3C N M H GRRINELLO M, et al子,置于尺寸為10.25A的周期性立方箱中;13Ab initio liquid water [J]. J Chem Phys, 1993, 99(11)mo/L的溶液由27個水分子和8個孤立的KOD分9080-9089子構(gòu)成,置于尺寸為10.15A的周期性立方箱中。[4] SILVESTRELLI P L, BERNASCONI M and PAR石油大學(xué)學(xué)報(自然科學(xué)版)2005年8月RINELLO M. Ab initio infrared spectrum of liquidSalvetti correlation-energy formula into a functional ofter[J. Chem Phys Lett,1997,277(5-6):478-482the electron density[j]. Phys Rev B, 1988, 37(2): 785[5] TROUT B L and PARRINELLO M. X-ray spectroscopicnd quantum-chemical study of carbon tubes produced in [16] RICCI M A, NARDONE M, RICCI F P, et al.Miarc-discharge [J]. Chem Phys Lett, 1998, 288 (3-4)croscopic structure of low temperature liquid ammonia: aneutron diffraction experiment [J]. J Chem Phys 1995[6 SPRIK M. Computation of the pK of liquid water using102(19):7650-7655coordination constraints [J]. Chem Phys, 2000, 258(2-[17]O'REILLY D E, PETERSON E M and SCHEIE CE.Self-diffusion in liquid ammonia and deuteroammonia[7] SCHWEGLER E, GALLI G, GYGI F, et al. Dissocia[J]. J Chem Phys,1973,58(10):4072-4075tion of water under pressure [J]. Phys Rev Lett, 2001, [18] YANAGUCHI T, HIDAKA K and SOPER AK.The87(26):265-501structure of liquid methanol revisited: a neutron diffrac[ 8] IZVEKOV S and VOTH G. A Car- Parrinello moleculartion expenment at-80℃and+25℃[J].Moldynamics simulation of liquid water: new results [J].JPhys,1999,96(8):1159-1168Chem Phys,202,116(23):10372-10376[19] YANAGUCHI T, HIDAKA K and SOPER AK. Erra9] SILVESTRELLI P L and PARRINDLLO M. Watum: the structure of liquid methanol revisited a neu-[J]ron diffraction experiment at-80℃and+25℃[J]hys Rev Lett,99,82(16):3308-3311Mol Phys,1999,97(4):603-605[10] DELLAGO C, BOLHUIS P G, CSAJKA F S, et al. [20] MORRONE J A and TUKERMANME. Ab initio mlecTransition path sampling and the calculation of rate con-ular dynamics study of proton bility in liquid methanol[J]stants[J]. J Chem Phys,1998,108(5):1964-1977J Chem Phys,2002,117(9:4403-4413[11] DELLAGE C, BOLHUIS P G and CHANDLER D On [21] LAASONEN K, PASQUARELLO A, CAR R, et althe calculation of reaction rate constants in the transitionCar-Parrinello molecular dynamics with Vanderbilt ulpath ensemble [J]. J Chem Phys, 1999, 110(14): 6617rasoft pseudopotentials [J]. Phys Rev B, 1993, 47(16):10142-1015[12] DIRAISON M, MARTYNA G J and TUCKERMAN [22] TSUCHIDA E, KANADA Y and TSUKADA M.Den-ME. Simulation studies of liquid ammoniasity-functional study of liquid methanol[J]. Chem Phab initio, classical, and path-integral molecular dynam-tt,1999,311(3):236-240ics[J]. j Chem Phys,199,01(3):1096-1103[23] ZHU Z and TUCKERMAN M E. Ab initio molecular[13] LIU Y and TUCKERMAN M E. Protonic defectsdynamics investigation of the concentration dependencehydrogen bonded liquids: structure and dynamics in am-ofmonia and comparison with water[J]. J Phys Chem B,tion of the infrared spectrum [J]. J Phys Chen B2001,105(28):6598-66102002,106(33):8009-8018[14] BECKE A D. Density-functional exchange-energy ap- [24] NAILONG W. The Maximum Entropy Method [M]proximation with correct asymptotic behavior[J]. PhysNew York: Springer, 1997RevA,1988,38(6):3098-3100(編輯陳淑嫻」[15] LEE W Y C and PARR R C. Development of the Colle-油田開發(fā)過程中油層保護(hù)及解堵技術(shù)研究”達(dá)國際領(lǐng)先水平由中國石油大學(xué)石油工程學(xué)院蔣官澄博士為首的課題組與中國石化勝利油田有限公司孤東采油廠共同完成的油田開發(fā)過程中油層保護(hù)及解堵技術(shù)研究”,近日通過V凵中國煤化工家鑒定。鑒定委員會一致認(rèn)為,該項研究達(dá)到了同類研究的國際領(lǐng)先水平。CNMHG目前,該項研究成果已在孤東油田進(jìn)行了1877井次的現(xiàn)場試驗,提高了施工成功率和有效期,保護(hù)了油層,共創(chuàng)直接經(jīng)濟(jì)效益29664.67萬元,總的投入產(chǎn)出比達(dá)到1:2.7。(摘自中國石油大學(xué)校園網(wǎng))