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