MOLPRO 分子專家程序包是一套由各種進行分 子體系電子結(jié)構計算的量子化學從頭算程序組 成的完整體系。它主要由Hans-Joachim Werner 和Peter Knowles 設計維護,同時MOLPRO 分 子專家程序包之中也包含很多其他作者的貢 獻。與其他常用的量子化學軟件包不同的是, MOLPRO 分子專家程序包側(cè)重于高精度的計 算,這一點主要體現(xiàn)在對電子相關廣泛的處理 上,使用的方法包括多參考態(tài)的組態(tài)相互作用 方法(MRCI)、耦合簇方法(CC)以及其他 相關方法。使用發(fā)展的密度擬合近似技術 (Density Fitting),結(jié)合局域相關方法(Local correlation methods),可以異常顯著的降低隨 著分子尺度的增加而增加的計算量,顯著的改 善計算量對于分子尺度的標度(scaling)。這 樣,相對于其他軟件包,MOLPRO 分子專家程 序包可以對更大的分子體系進行精確的量子化 學從頭計算。更進一步的,新發(fā)展的直接相關 MP2 和直接相關耦合簇方法可以得到非常接近 基組極限(basis set limit)的結(jié)果。
Molpro的功能與特點
★ 能量計算
*生成一般收縮對稱化高斯基的積分(spdfghi)
*有效核勢(ECP)和核極化勢(CPP)
*閉殼層和開殼層(包括自旋約束和非約束)自 洽場方法(RHF, UHF)
*Kohn-Sham 框架下的密度泛函理論(DFT),支持大量一般梯度近似(GGA)交換泛函和相關泛函
*單參考態(tài)電子相關方法(閉殼層方法有: Moller-Plesset 微擾理論(MP2-MP4),耦合簇方 法(CCSD(T)),二次組態(tài)相互作用方法
(QCISD(T)以及Brueckner 耦合簇方法 (BCCD(T))。開殼層方法有:開殼層自旋非 約束和自旋約束的MP2 方法(UHF-UMP2, ROHF-RMP2)和(部分約束的)開殼層耦合簇 方法(RHF-UCCSD(T), RHF-RCCSD(T)))
*計算量相對于分子尺度增長為線性標度的局域 電子相關方法(LMP2, LQCISD(T), LCCSD(T))
*開殼層和閉殼層的直接相關(Explicitly correlated) 方法(MP2-F12, LMP2-F12, CCSD(T)-F12)
*可以直接計算激發(fā)能的CIS, CC2 和EOM-CCSD 方法
*多組態(tài)自洽場(MCSCF)方法,包括CASSCF 和RASSCF,并可以處理一般態(tài)平均(general state-averaged)的能量泛函
*內(nèi)部收縮的多參考態(tài)組態(tài)相互作用(MRCI), 包括多種變形,例如:MR-ACPF,MR-AQCC
*一般的多參考態(tài)二階和三階微擾理論(MRPT2/ CASPT2, MP-PT3/CASPT3),多態(tài)(multistate) CASPT2
*全組態(tài)相互作用方法(FCI)
*使用直接積分技術(Integral-direct)的Hartree- Fock,DFT,MCSCF 和電子對相關方法(MPn, CCSD, MRCI)。
*使用密度擬合技術(Density Fitting)的Hartree- Fock,DFT,MP2,LMP2,LCCSD(T),MP2- F12,LMP2-F12 方法。
★ 梯度計算和幾何結(jié)構優(yōu)化
*可以計算解析能量梯度的方法有:SCF, DFT, MP2, QCISD(T), LMP2, LQCISD, SA-MCSCF/ CASSCF, MR-PT2/CASPT2。密度擬合方法已應 用于SCF,DFT 和LMP2 的梯度計算。
*閉殼層HF 和(單態(tài),single state)MCSCF/ CASSCF 的解析能量二階導數(shù)(Hessians)
*在所有方法中都能計算數(shù)值能量梯度和數(shù)值能量二階導數(shù)(Hessians)
*自動進行幾何構型優(yōu)化,過渡態(tài)尋找,反應路徑跟蹤。
*MCSCF 的解析非絕熱耦合(non-adiabatic coupling)矩陣元計算,MRCI 的數(shù)值非絕熱耦合(non-adiabatic coupling)矩陣元計算。
*自動計算簡諧和非簡諧振動頻率,振動強度和 熱力學性質(zhì)。
*點電荷陣列的能量梯度(為做QM/MM優(yōu)化)
★ 性質(zhì)計算和波函數(shù)分析
*以下方法可以計算多種單電子性質(zhì):HF, DFT, MCSCF/CASSCF, MRCI, EOM-CCSD。使用以 下方法計算多種躍遷性質(zhì)(Transition
Properties):MCSCF/CASSCF, MRCI, EOMCCSD( CASSCF 和MRCI 也可計算不同的軌道 組成的波函數(shù)之間的躍遷)。MCSCF/CASSCF 能計算一些雙電子性質(zhì)(例如角動量算符)
*CASSCF 和MRCI 的自旋-軌道耦合(Spin-orbit coupling)。
*使用有效核勢方法或使用Breit- Pauli 形式。
*Mulliken 電荷分析和多極矩分布分析。
*NBO 和自然布居數(shù)(natural population)分析
*CASSCF 波函數(shù)的價鍵(Valence-Bond)分析和 能量化的價鍵波函數(shù)計算。
*軌道局域化(Boys,Pipek-Mezey)
★ 接口
*可視化程序:MOLDEN,XMOL,…
*QM/MM 程序:Chemshell, QoMMMa, QMPOT
*WebMO (http://www.webmo.net)
*MRCC 程序(Kallay & Gauss)
*以XML 格式輸出CML 幾何構型和所有結(jié)果,基 組,軌道,可以在瀏覽器上看到其中包含的三 維幾何構型