應(yīng)用多中心分區(qū)方法構(gòu)建H3分子反應(yīng)勢(shì)能面

戴 偉 陳柳楊 鄭利敏 楊明暉,*

(1湖北第二師范學(xué)院物理與機(jī)電工程學(xué)院, 武漢 430205;2中國(guó)科學(xué)院武漢物理與數(shù)學(xué)研究所, 武漢 430071)

應(yīng)用多中心分區(qū)方法構(gòu)建H3分子反應(yīng)勢(shì)能面

戴 偉1,2陳柳楊2鄭利敏2楊明暉2,*

(1湖北第二師范學(xué)院物理與機(jī)電工程學(xué)院, 武漢 430205;2中國(guó)科學(xué)院武漢物理與數(shù)學(xué)研究所, 武漢 430071)

勢(shì)能面在分子反應(yīng)動(dòng)力學(xué)的研究中起著非常重要的作用. 本文提出了一種新的勢(shì)能面構(gòu)建方法——多中心分區(qū)法. 通過(guò)比較London-Eyring-Polanyi-Sato (LEPS)勢(shì)能、多體展開勢(shì)能、置換對(duì)稱不變多項(xiàng)式三種方法, 確定了H3分子的最佳勢(shì)能函數(shù)表達(dá)形式, 并應(yīng)用準(zhǔn)經(jīng)典軌線方法分析了勢(shì)能面的合理性, 結(jié)果表明置換對(duì)稱不變多項(xiàng)式能很好地描述H3分子的勢(shì)能面特征. 結(jié)合置換對(duì)稱不變多項(xiàng)式和本文提出的多中心分區(qū)方法,可以有效改善H3分子勢(shì)能面的精度并可能推廣到高維反應(yīng)勢(shì)能面.

勢(shì)能面擬合; 準(zhǔn)經(jīng)典軌線法; 多中心分區(qū)

1 引 言

隨著實(shí)驗(yàn)技術(shù)的快速發(fā)展, 分子反應(yīng)動(dòng)力學(xué)實(shí)驗(yàn)方法1取得了長(zhǎng)足進(jìn)步, 理論研究能夠解釋并預(yù)測(cè)基元化學(xué)反應(yīng)的動(dòng)態(tài)結(jié)構(gòu)、反應(yīng)過(guò)程和反應(yīng)機(jī)理,是實(shí)驗(yàn)研究的重要補(bǔ)充手段. 勢(shì)能面是理論研究中的一個(gè)重要概念, 它描述反應(yīng)體系中多個(gè)原子間的相互作用, 因此需要構(gòu)造準(zhǔn)確的勢(shì)能面以盡可能真實(shí)表達(dá)分子反應(yīng)的勢(shì)能面,2通常構(gòu)造勢(shì)能面的方法可以分為兩類, 一種是插值方法,3–23Collins等3–15在高維插值勢(shì)能面方面做了很多系統(tǒng)的工作, 還有應(yīng)用三次樣條插值、16,17Shepard 插值、18,19Modified Shepard 插值20–22來(lái)構(gòu)建分子反應(yīng)勢(shì)能面的研究工作,通過(guò)高維插值構(gòu)建勢(shì)能面, 需要較多的從頭算數(shù)據(jù)點(diǎn). 另一種是函數(shù)擬合,23–30Corchado等24通過(guò)解析表達(dá)式擬合構(gòu)建了H+CH4勢(shì)能面, Varandas等25通過(guò)雙多體展開(DMBE2001)構(gòu)建了H3勢(shì)能面, Bowman等26通過(guò)選擇滿足置換對(duì)稱性的多項(xiàng)式作為勢(shì)能函數(shù), 構(gòu)建了置換對(duì)稱不變多項(xiàng)式(PIP)勢(shì)能面, 以及最近得到廣泛應(yīng)用的的神經(jīng)網(wǎng)絡(luò)方法.27–30選擇合適的勢(shì)能函數(shù)形式能夠有效地提高勢(shì)能面精度.

在勢(shì)能面構(gòu)建過(guò)程中, 基于誤差是按區(qū)域分布這一特點(diǎn), 我們認(rèn)為如果按區(qū)域多中心分區(qū)擬合,將會(huì)有利于改善勢(shì)能面的精度. 基于這一思想, 本文發(fā)展了多中心分區(qū)擬合方法并應(yīng)用到H3分子反應(yīng)勢(shì)能面的構(gòu)建中, 相比高維插值方法, 節(jié)約計(jì)算資源, 相比傳統(tǒng)的函數(shù)擬合方法, 提高勢(shì)能面的精度,該項(xiàng)工作可為研究較大體系的反應(yīng)動(dòng)力學(xué)提供一個(gè)新的思路.

2 H3分子勢(shì)能面的解析表達(dá)式

2.1 訓(xùn)練集的數(shù)據(jù)來(lái)源及要求

本文H3結(jié)構(gòu)與能量的訓(xùn)練集數(shù)據(jù)由DMBE2001勢(shì)能面25產(chǎn)生. 其中R1, R2, R3取值區(qū)間為0.05–0.5 nm,步長(zhǎng)為0.01 nm, 坐標(biāo)的選取如圖1所示. 剔除能量高于40000 cm–1的能量點(diǎn), 同時(shí)剔除不能構(gòu)成幾何形狀的數(shù)據(jù)點(diǎn).

2.2 勢(shì)能函數(shù)形式的選取

改善勢(shì)能面擬合精度的關(guān)鍵在于選擇合適的勢(shì)能函數(shù)形式. Connor31指出, 理想的勢(shì)能函數(shù)應(yīng)滿足以下條件:

(1) 能準(zhǔn)確描述反應(yīng)物和產(chǎn)物通道的漸進(jìn)性質(zhì);

(2) 能正確顯示出體系的對(duì)稱性;

圖1 H3體系坐標(biāo)Fig.1 Coordinates for H3system

(3) 在有實(shí)驗(yàn)數(shù)據(jù)或者(非經(jīng)驗(yàn)的)理論數(shù)據(jù)的反應(yīng)區(qū)域能準(zhǔn)確描述體系的能量;

(4) 在沒(méi)有實(shí)驗(yàn)數(shù)據(jù)和理論數(shù)據(jù)的反應(yīng)區(qū)域能合理描述體系的行為;

(5) 保證各個(gè)區(qū)域光滑連接;

(6) 勢(shì)能函數(shù)及其微分的代數(shù)形式盡可能簡(jiǎn)單; (7) 盡可能少的數(shù)據(jù)點(diǎn)就能達(dá)到擬合的精度.本文討論以下三種勢(shì)能函數(shù)在擬合H3反應(yīng)勢(shì)能面的精度.

2.2.1 LEPS勢(shì)能

Eyring和Polanyi32利用London對(duì)三原子體系的量子力學(xué)勢(shì)能近似構(gòu)建了London-Eyring-Polanyi (LEP)勢(shì)能面, Sato33又在他們的基礎(chǔ)上進(jìn)行了修正,去掉了勢(shì)壘頂端不合理的勢(shì)阱, 得到了London-Eyring-Polanyi-Sato (LEPS)勢(shì)能面. LEPS也適合于多原子分子, 由于形式比較簡(jiǎn)單, 參數(shù)少, 能夠反映出體系的主要特征, 在分子動(dòng)力學(xué)發(fā)展的早期應(yīng)用特別廣泛. 三原子體系的能量采用London形式表示:

其中, Q為庫(kù)能積分, J為交換積分,

式中1Ei為雙原子Morse函數(shù),3Ei為反Morse函數(shù), 表達(dá)式如下:

式中Di、和βi分別表示雙原子分子基態(tài)離解能、平衡核間距和光譜常數(shù), 可通過(guò)從頭計(jì)算確定, Si是針對(duì)不同體系引入的可調(diào)參數(shù), βi為光譜常數(shù), 可表示為

c為光速, ωe為振動(dòng)光譜常數(shù), μ為折合質(zhì)量.

為確定LEPS勢(shì)能中的可調(diào)參數(shù)Si, 掃描了一條均方根偏差(RMSD)隨Si變化的曲線, 由圖2可以看出, Si取0.10時(shí), 均方根偏差最小. 構(gòu)建的LEPS勢(shì)能面如圖3所示, 勢(shì)能中所用參數(shù)如表1所示.

2.2.2 多體展開勢(shì)

多體展開法34是構(gòu)建多原子分子勢(shì)能面的常見(jiàn)方法之一, 它由單體項(xiàng)、二體項(xiàng)和多體項(xiàng)構(gòu)成. 表達(dá)形式如下:

其中,

其中RAB為雙原子體系A(chǔ)B的鍵長(zhǎng), DeAB是離解能, ReAB是平衡核間距, αAB是Morse參數(shù). 同理可得AC、BC雙原子系統(tǒng)的Morse勢(shì)表達(dá)式. H3體系的離解能為0.293526 hartree, 平衡核間距為1.390 Bohr, Morse參數(shù)為8.10036.

式中的為非線性系數(shù), ρBC和 ρAC的表達(dá)式與此類似, d是擬合參數(shù), M是展開的階數(shù), 本文中M = 5.

圖2 均方根偏差(RMSD)隨Si參數(shù)的變化關(guān)系Fig.2 Relationship between the root mean square deviation (RMSD) and the Siparameter

圖3 H3分子30°的勢(shì)能等值圖Fig.3 Contour plots of potential energy surface for H3at 30°

表1 LEPS勢(shì)能參數(shù)數(shù)值Table1 Values of LEPS potential parameters

2.2.3 置換對(duì)稱不變多項(xiàng)式(PIP)方法

美國(guó)Emory大學(xué)化學(xué)系Joel M. Bowman教授在構(gòu)建多原子分子勢(shì)能面時(shí), 選擇滿足置換對(duì)稱性的多項(xiàng)式作為勢(shì)能函數(shù), 對(duì)勢(shì)能面精度有很大改善,其擬合表達(dá)式27如下:

式中a為可調(diào)參數(shù), 在本文中取1.8 bohr, n1, n2, n3= 1, 2,, 10且n1+ n2+ n3≤ 18, N為擬合表達(dá)式中待定參數(shù)總和, 文中取970項(xiàng).

2.2.4 計(jì)算結(jié)果

為了進(jìn)一步研究勢(shì)能函數(shù)帶來(lái)數(shù)值誤差的分布狀況, 本文統(tǒng)計(jì)了三種不同勢(shì)能函數(shù)在不同散射角度下的均方根偏差(見(jiàn)表2). 結(jié)果表明: (1) PIP方法能準(zhǔn)確描述H3分子勢(shì)能面的結(jié)構(gòu)特征, 統(tǒng)計(jì)誤差最小; (2) 所有勢(shì)能函數(shù)在60°散射時(shí), 誤差最大, 分析原因可能是DMBE2001勢(shì)能面25提供的數(shù)據(jù)集在這一角度本身不能很好描述H3體系的相互作用, 需要通過(guò)從頭計(jì)算做進(jìn)一步的驗(yàn)證.

圖4、圖5、圖6是三種擬合表達(dá)式60°時(shí)的勢(shì)能面等值圖, 由圖可以看出, 勢(shì)能面整體趨勢(shì)區(qū)別不大, 第一種擬合表達(dá)式近程排斥區(qū)域最大, 第二種擬合表達(dá)式近程排斥區(qū)域最小, 在三體解離區(qū), 也存在一定差別, 其它區(qū)域三種形式基本吻合.

表2 不同散射角下勢(shì)能面的RMSDTable2 RMSD of potential energy surface under different scattering angles

圖7是三種擬合勢(shì)能面及LEPS勢(shì)能面的準(zhǔn)經(jīng)典軌線計(jì)算結(jié)果, 總體來(lái)說(shuō), 與文獻(xiàn)提供的勢(shì)能面在

圖4 公式1擬合下60°時(shí)的勢(shì)能面等值圖Fig.4 Contour plots of potential energy surface for the Function 1 at 60°

圖5 公式2擬合下60°時(shí)的勢(shì)能面等值圖Fig.5 Contour plots of potential energy surface for the Function 2 at 60°

圖6 公式3擬合下60°時(shí)的勢(shì)能面等值圖Fig.6 Contour plots of potential energy surface for the Function 3 at 60°

圖7 不同勢(shì)能函數(shù)表達(dá)式的準(zhǔn)經(jīng)典軌線結(jié)果比較Fig.7 Comparison of the quasi-classical trajectory results between different potential energy functions

基本趨勢(shì)上是一致的, LEPS勢(shì)能面參數(shù)少, 使用方便, 但高能區(qū)域偏差較大, 第一種擬合表達(dá)勢(shì)能面在低能區(qū)域偏差較大, 不適宜用來(lái)描述H3分子相互作用, 第二種和第三種勢(shì)能面在低能區(qū)與文獻(xiàn)值符合很好, 高能區(qū)域, 第三種勢(shì)能面符合更好, 綜合動(dòng)力學(xué)結(jié)果分析, PIP方法能更好地描述H3體系的相互作用.

3 應(yīng)用多中心分區(qū)擬合構(gòu)建H3分子勢(shì)能面

3.1 多中心分區(qū)擬合的思想由來(lái)

圖8是本文構(gòu)建勢(shì)能面與數(shù)據(jù)測(cè)試集的差分圖(散射角60°), 圖中表示的是采用PIP方法構(gòu)建的勢(shì)能面與測(cè)試數(shù)據(jù)的差值, 用來(lái)描述勢(shì)能面的回歸特性. 由圖可以看出, 誤差呈區(qū)域分布, 大部分區(qū)域吻合較好, 誤差不大, 只有在鞍點(diǎn)區(qū)域到三體離解區(qū)域過(guò)渡時(shí), 兩種勢(shì)能差別較大. 擬合勢(shì)能面與測(cè)試數(shù)據(jù)集的吻合度按區(qū)域分布, 有些地方高于測(cè)試值,有些地方低于測(cè)試值, 鑒于誤差是按區(qū)域分布這一特點(diǎn), 我們預(yù)測(cè)如果在構(gòu)建勢(shì)能面時(shí), 按區(qū)域多中心分區(qū)擬合, 對(duì)于提高勢(shì)能面的精度會(huì)有改善.

圖8 置換對(duì)稱不變多項(xiàng)式勢(shì)能面與測(cè)試數(shù)據(jù)的差值等高圖Fig.8 Contour plots of the difference between the permutation-invariant polynomial method and the test data

3.2 多中心展開構(gòu)建H3勢(shì)能面

3.2.1 坐標(biāo)的選取

為了描述H + H2反應(yīng)體系, 選取H2為中心原子,通過(guò)兩個(gè)鍵長(zhǎng)(R1, R2)和一個(gè)鍵角(θ)表示, 坐標(biāo)選取如圖1所示.

3.2.2 計(jì)算方法

3.2.2.1 訓(xùn)練集數(shù)據(jù)的選取

采用DMBE2001勢(shì)能面25作為測(cè)試集, 取R1, R2∈ [0.05 nm, 0.5 nm]均勻100 個(gè)點(diǎn); 均勻取θ ∈[0°, 180°]30個(gè)點(diǎn), 得到各數(shù)據(jù)點(diǎn)的能量值.

3.2.2.2 分區(qū)原則

在H + H2→ H2+ H的最小反應(yīng)路徑上取7個(gè)(R1, R2)(單位: nm), 作為七個(gè)中心點(diǎn), 計(jì)算訓(xùn)練集中的數(shù)據(jù)點(diǎn)(R1, R2)與七個(gè)中心的距離, 依據(jù)最近及次近距離將數(shù)據(jù)集分成7個(gè)區(qū)域, 當(dāng)R1= R2時(shí), 存在兩個(gè)次小數(shù)據(jù)點(diǎn)與各中心點(diǎn)間距值, 將此點(diǎn)納入相應(yīng)的兩個(gè)分區(qū), 擬合數(shù)據(jù)的分布與誤差統(tǒng)計(jì)如表3所示.

3.2.2.3 邊界處理

合并擬合分區(qū), 將各數(shù)據(jù)點(diǎn)根據(jù)公式(18)計(jì)算相應(yīng)的能量值: centq為此數(shù)據(jù)點(diǎn)所在數(shù)據(jù)域的中心點(diǎn), q = 1, 2, 3, , 7. wk為分區(qū)能量在總能量中的權(quán)重, 表示如下:

表3 擬合數(shù)據(jù)分布及誤差統(tǒng)計(jì)Table3 Distribution of the fitting data and the error statistics

圖9 多中心分區(qū)法計(jì)算準(zhǔn)經(jīng)典軌線的結(jié)果Fig.9 Results of quasiclassical trajectory by “multi-center partition” method

當(dāng)此數(shù)據(jù)點(diǎn)僅存在于兩個(gè)分區(qū)時(shí), w3= 0, vfitj3= 0, vfitjk表示處在第j個(gè)分區(qū)時(shí)，第k次近中心點(diǎn)對(duì)總勢(shì)能面的貢獻(xiàn). 則多中心分區(qū)能量(單位為cm–1)為:

3.3 多中心展開構(gòu)建H3勢(shì)能面

采用多中心分區(qū)擬合后, 構(gòu)建勢(shì)能面的準(zhǔn)經(jīng)典軌跡計(jì)算結(jié)果與文獻(xiàn)完全一致, 如圖9所示. 由圖9可以看出, 多中心分區(qū)擬合后, 構(gòu)建勢(shì)能面與原勢(shì)能面回歸性很好, 相比全域擬合(圖7), 精度進(jìn)一步提高, 無(wú)論是在低能區(qū)域, 還是在高能區(qū)域, 都與文獻(xiàn)值完全吻合.

4 結(jié) 論

構(gòu)建了H3體系的LEPS勢(shì)能面, 采用三種不同勢(shì)能函數(shù), 構(gòu)建了H3分子勢(shì)能面, 為多中心分區(qū)方法確定了最佳擬合函數(shù)形式, 應(yīng)用準(zhǔn)經(jīng)典軌跡法驗(yàn)證了勢(shì)能面的合理性, 最后采用多中心分區(qū)擬合, 構(gòu)建了H3分子勢(shì)能面. 研究結(jié)果表明:

(1) LEPS勢(shì)能面參數(shù)少, 除Si參數(shù)需調(diào)節(jié)確定,其余參數(shù)都可通過(guò)從頭計(jì)算得到, 勢(shì)能函數(shù)形式簡(jiǎn)單, 使用方便, 在定性計(jì)算中可以廣泛使用, 但不適宜描述高精度反應(yīng)勢(shì)能面, 動(dòng)力學(xué)驗(yàn)證誤差較大.

(2) 比較了三種不同的勢(shì)能函數(shù)形式. 訓(xùn)練集采用了213200個(gè)構(gòu)型的能量信息, 截?cái)嗄転?0000 cm–1,數(shù)據(jù)量大, 能域?qū)? PIP方法統(tǒng)計(jì)誤差最小, 動(dòng)力學(xué)信息最吻合, 可以用來(lái)描述H3分子勢(shì)能面的能量特征.

(3) 擬合勢(shì)能面的精度與表達(dá)式展開的項(xiàng)數(shù)、數(shù)據(jù)集的數(shù)據(jù)量大小密切相關(guān). 多中心分區(qū)之后,每一區(qū)域數(shù)據(jù)集變少, 改變了數(shù)據(jù)集和擬合參數(shù)的相對(duì)比, 可以降低擬合表達(dá)式展開項(xiàng)的要求, 同時(shí)有效改善擬合精度. 從計(jì)算結(jié)果可以看出: 多中心分區(qū)擬合法相比全域擬合, 勢(shì)能面回歸性很好, 無(wú)論是在低能區(qū)域, 還是在高能區(qū)域, 都與文獻(xiàn)值吻合很好.

(1)Sun, Z. F.; Gao, Z.; Wu, X. K.; Tang, G. Q.; Zhou, X. G.; Liu, S. L. Acta Phys. -Chim. Sin. 2015, 31, 829. [孫中發(fā), 高 治, 吳向坤, 唐國(guó)強(qiáng), 周曉國(guó), 劉世林. 物理化學(xué)學(xué)報(bào), 2015, 31, 829.] doi: 10.3866/PKU.WHXB201503041

(2)Li, W.; Qu, J. Y.; Zhao, X. S. Acta Phys. -Chim. Sin. 2003, 19 (8), 751. [李 巍, 屈軍艷, 趙新生. 物理化學(xué)學(xué)報(bào), 2003, 19 (8), 751.] doi: 10.3866/PKU.WHXB20030816

(3)Le, H. A.; Frankcombe, T. J.; Collins, M. A. J. Phys. Chem. A 2010, 114 (40), 10783. doi: 10.1021/jp1060182

(4)Ramazani, S.; Frankcombe, T. J.; Andersson, S.; Collins, M. A. J. Chem. Phys. 2009, 130 (24), 244302. doi: 10.1063/1.3156805

(5)Moyano, G. E.; Jones, S. A.; Collins, M. A. J. Chem. Phys. 2006, 124 (12), 124318. doi: 10.1063/1.2181571

(6)Moyano, G. E.; Collins, M. A. Theor. Chem. Acc. 2005, 113 (4), 225. doi: 10.1007/s00214-004-0626-8

(7)Evenhuis, C. R.; Collins, M. A.; Lin, X.; Zhang, O. H. Amer. Chem. Soc. 2004, 227, 259.

(8)Moyano, G. E.; Pearson, D.; Collins, M. A. J. Chem. Phys. 2004, 121 (24), 12396. doi: 10.1063/1.1810479

(9)Castillo, J. F.; Aoiz, F. J.; Ba?ares, L.; Collins, M. A. J. Phys. Chem. A 2004, 108 (32), 6611. doi: 10.1021/jp048366b

(10)Moyano, G. E.; Collins, M. A. J. Chem. Phys. 2003, 119 (11), 5510. doi: 10.1063/1.1599339

(11)Fuller, R. O.; Bettens, R. P. A.; Collins, M. A. J. Chem. Phys. 2001, 114 (24), 10711. doi: 10.1063/1.1377602

(12)Song, K.; Collins, M. A. Chem. Phys. Lett. 2001, 335 (5), 481.

(13)Collins, M. A.; Zhang, D. H. J. Chem. Phys. 1999, 111 (22), 9924. doi: 10.1063/1.480344

(14)Bettens, R. P. A.; Hansen, T. A.; Collins, M. A. J. Chem. Phys. 1999, 111 (14), 6322. doi: 10.1063/1.479937

(15)Jordan, M. J. T.; Collins, M. A. J. Chem. Phys. 1996, 104 (12), 4600. doi: 10.1063/1.471207

(16)Li, H.; Le, R. R. J. J. Chem. Phys. 2006, 125 (4), 044307. doi: 10.1063/1.2212933

(17)Li, A. Y.; Xie, D. Q. J. Chem. Phys. 2010, 133 (14), 144306. doi: 10.1063/1.3490642

(18)Wu, T.; Manthe, U. J. Chem. Phys. 2003, 119 (1), 14. doi: 10.1063/1.1577328

(19)Ishida, T.; Schatz, G. C. J. Chem. Phys. 1997, 107 (9), 3558. doi: 10.1063/1.474695

(20)Takata, T.; Taketsugu, T.; Hirao, K.; Gordon, M. S. J. Chem. Phys. 1998, 109 (11), 4281. doi: 10.1063/1.477032

(21)Crespos, C.; Collins, M. A. J. Chem. Phys. 2004, 120 (5), 2392. doi: 10.1063/1.1637337

(22)Wang, M.; Sun, X.; Bian, W.; Cai, Z. J. Chem. Phys. 2006, 124 (23), 234311. doi: 10.1063/1.2203610

(23)Dawes, R.; Thompson, D. L.; Guo, Y.; Wagner, A. F.; Minkoff, M. J. Chem. Phys. 2007, 126 (18), 184108. doi: 10.1063/1.2730798

(24)Corchado, J. C.; Bravo, J. L.; Espinosa-Garcia, J. J. Chem. Phys. 2009, 130 (18), 184314. doi: 10.1063/1.3132223

(25)Varandas, A. J. C.; Brown, F. B.; Mead, C. A. J. Chem. Phys. 1987, 86, 6258. doi: 10.1063/1.452463

(26)Bowman, J. M.; Czako, G.; Fu, B. Phys. Chem. Chem. Phys. 2011, 13 (18), 8094. doi: 10.1039/c0cp02722g

(27)Pukrittayakamee, A.; Malshe, M.; Hagan, M.; Raff, L. M.; Narulkar, R.; Bukkapatnum, S.; Komanduri, R. J. Chem. Phys. 2009, 130 (13), 134101. doi: 10.1063/1.3095491

(28)Le, H. M.; Raff, L. M. J. Phys. Chem. A 2010, 114 (1), 45. doi: 10.1021/jp907507z

(29)Sumpter, B. G.; Noid, D. W. Chem. Phys. Lett. 1992, 192 (5), 455.

(30)Blank, T. B.; Brown, S. D.; Calhoun, A. W.; Doren, D. J. J. Chem. Phys. 1995, 103 (10), 4129. doi: 10.1063/1.469597

(31)Connor, J. N. L. Comput. Phys. Commun. 1979, 17, 117. doi: 10.1016/0010-4655(79)90075-4

(32)Eyring, H.; Polanyi, M. Phys. Chem. Abt. B 1931, 12, 279.

(33)Sato, S. J. Chem. Phys. 1955, 23 (12), 2465.

(34)Aguado, A.; Paniagua, M. J. Chem. Phys. 1992, 96 (2), 1265. doi: 10.1063/1.462163

Application of the Multi-Center Partition Method to Construct the Potential Energy Surface of H3

DAI Wei1,2CHEN Liu-Yang2ZHENG Li-Min2YANG Ming-Hui2,*
(1School of Physics and Mechanical & Electrical Engineering, Hubei University of Education, Wuhan 430205, P. R. China;2Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, P. R. China)

The potential energy surface plays an important role in studying molecular reaction dynamics. In this work, a new method, namely the “multi-center partition” method, is proposed to construct the potential energy surface of H3. The optimized function is first determined by comparing the London-Eyring-Polanyi-Sato (LEPS) potential, the many-body expansion potential, and the permutation-invariant polynomial potential. This comparison shows that the permutation-invariant polynomial fitting proposed by Bowman is the most efficient method for describing the topology of the H3system. The quasi-classical trajectory method is used to analyze the rationality of those potential energy surfaces. By combining the multi-center partition method with the permutation-invariant polynomial method, the accuracy of the H3molecular potential energy surface is greatly improved and could possibly be used in the fitting of potential energy surfaces in other systems.

Potential energy surface fitting; Quasi-classical trajectory method; Multi-center partition

O641.3

10.3866/PKU.WHXB201509143

Received: July 6, 2015; Revised: September 14, 2015; Published on Web: September 14, 2015.

*Corresponding author. Email: yangmh@wipm.ac.cn; Tel: +86-27-87197783.

This project was supported by the China Scholarship Council (201408420174), Science and Technology Research Program of the Education Department, Hubei Province, China (Q20133005), and Natural Science Foundation of Hubei Province, China (2014CFB428, 2015CFB502).

國(guó)家留學(xué)基金委項(xiàng)目(201408420174), 湖北省教育廳科學(xué)技術(shù)項(xiàng)目(Q20133005)和湖北省自然科學(xué)基金(2014CFB428, 2015CFB502)資助

?Editorial office of Acta Physico-Chimica Sinica