1、自旋污染的基本概念以及對計算結果的影響關于Hartree-Fock計算（通常采用Rootaan的自洽場（SCF）方式）介紹性的描述集中在單 重態體系，在這樣體系中所有電子自旋都已經配對。通過假設計算是限制在每個占據軌道有兩 個電子計算，計算可以相對地容易進行。這通常是指自旋限制的Hartree-Fock計算或RHF。對于自旋多重度不為一的體系，是不可能象自旋多重度為一的體系那樣用RHF方式的。在 一個非限制性的計算中，有兩組完全 的軌道軌道，一組是alpha電子的，一組是beta電子的。 通常這兩組軌道用同樣的基函數但是不同的分子軌道參數。非限制性計算的優點是他們計算的效率非常高。缺點是波函數

2、不再是總自旋的本征函數， 因此計算中可能會導入一些誤差。這種誤差成為自旋污染。自選污染如何影響計算結果：自旋污染導致在我們看來是期望的波函數中混合進了一點其他的自旋態。這種情況下由于 有了較多的變分自由度，有時候會導致輕微的降低計算的總能量。更多的情況是因為較高能量態 的混入導致總能量的輕微升高。然而這些變化是不正確的波函數 的結果.因為這不是一種系統 誤差，不同態間的能量差可能受到相反的影響。高自旋的污染可以影響稽核和布局數分析，和 顯著地影響自旋密度。作為自旋污染出現的一種檢查，大多ab initio程序會打印出期望的總自旋值，如果沒有自 旋污染，它會等于s（s+1），這里s等于1/2的未

3、成對電子數倍。由有機分子計算經驗得到的 一個簡單規律是，當和s（s+1）的差別小于10%時，自旋污染的影響可以忽略。雖然這個規則 提供了一種快速檢驗的方法，通常用實證據或更嚴格的計算來檢查你的結果是明智的。自旋污染通常在非限制的Hartree-Fock （UHF）和非限制性的Mller-Plesset（UMP2, UMP3, UMP4）計算中出現。在DFT計算中很少會發現顯著的自旋污染，甚至在采用非 限制性的Kohn-Sham軌道的情況.（完整word版）自旋污染的基本概念以及對計算結果的影響 非限制性計算通常結合一步自旋湮滅步驟，這個步驟從計算的波函數中除去很大百分比的自旋污染。它有助于最小

4、化自旋污染但是無法完全阻止自旋污染。的最后值是對出現的自旋污 染的最好的檢查。在Gaussian中，選項” iop（5/14=2）” 告訴程序用湮滅后的波函數產生布 局數分析.我不知道是否有程序在執行幾何優化過程中使用湮滅的波函數.限制性開殼成計算有 可能進行自旋限制性開殼層計算（ROHF）.這種方法的好處是沒有自旋污染。缺點是會需要額外 的CPU時間在對單占據的和雙占據的軌道都進行正確處理以及它們之間的相互作用上。作為采 用的數學算法的一個結果，ROHF計算給出好的總能量和波函數但是單占據的軌道不嚴格服從 Koopman 定理.當自旋污染導致的誤差不可接受時，先執行開殼層計算是給出可靠波函數

5、的最好方法.在Gaussian程序中，限制性開殼層計算可以執行HartreeFock，密度泛函理論，MP2和 一些半經驗波函數計算。但是ROMP2方式不知吃解析梯度，因此最快的方式是在其他方式優 化的幾何上進行單點能計算。如果幾何優化一定要在這樣的理論水平下進行，必須采用一種 非基于梯度的方式，如FletcherPowell優化。（注意在G94手冊中對所有情況都還沒有這種 功能）自旋投影方式另一種近似是先運行一種非限制性計算然后在得到的波函數上投影掉自旋污染 PUHF,PUMP2）。自旋投影結果并不給出用限制性開殼層計算得到的能量。這是因為非限制性的軌道被優化 來描述被污染的自旋態而不是被優化

6、來描述自旋投影的態。用自旋強制的UHF方式（SUHF）可以得到類似的效果。在這種方式中，通過一個拉格朗日乘 因子法，UHF波函數中的自旋污染誤 差被限制。當乘的因子趨于無窮時，這種方法完全地除去 了自旋污染。實際上小的正值就可以除去大多數的自旋污染.(完整word版)自旋污染的基本概念以及對計算結果的影響 半電子近似對自由基計算，半經驗程序常用半電子近似。半電子方式是RHF計算中處理一個單占據軌道的數學技術。它在擁有一個近似的波函數和軌道能量的情況下能得到一致性的 總能量。因為采用的是單行列式計算，因此沒有自旋污染。一致性的總能量使得這種方法可以計算單態一三態能隙：對單重態用RHF，對三重態用

7、半電 子近似。半電子近似計算不服從 Koopman定理。也不能得到自旋密度。Mulliken布居分析通 常是相當合理的.結合到具體的計算方法，使用CASSCF,CASPT2，MRMP2，MRCI等multireference方法 可以比較有效的消除spin contamination的問題.根據個人體會使用PUMP2(projected UMP2) 似乎也是不錯的選擇(在Gaussian中的實現為在route section直接使用UMP2(對開殼層體 系)關鍵詞，然后在結果中選用projected后的能量即可)，而且速度比基于multireference 的方法更快(CASSCF除外).ab

8、 initio的計算主要基于e rror cancelling比如反應熱，反應 勢壘等，沒有人對絕 對能量感興趣。anyway，雖然UHF有spin contamination的問題，但是 使用UHF計算的ground state的能量一般比ROHF更低，因為前者中LCAO的自由度更多.補充信息一些討論和結果在：W. Jo Hehre， Lo Radom， Po v。R. Schleyer, J. Ao Pople Ab Initio Molecular Orbita l Theory” Wiley (1986)An article that compares unrestricted， re

9、stricted and projected results isM. W. Wong, Lo Radom J。Phys。Chemo 99, 8582 (1995)Some specific examples and a discussion of the half electron method are given in To Clark A Handbook of Computational Chemistry” Wiley (1985)A more mathematical treatment can be found in the paperJ. S. Andrews， Do Jaya

10、tilake， R. G。Ao Bone，No C. Handy， R. D. Amos Chemo Physo Lett. 183， 423 （1991）SUHF results are examined inPo Ko Nandi, T. Kar， Ao B. Sannigrahi Journal of Molecular Structure （Theochem） 362, 69 （1996）誰能解釋一下什么是自旋污染？就我的理解：自旋污染表現在計算結果中操作S2 （2指平方）不等于其本征值：s （s+1）例如有一個單電子的體系，s2=0.75。但是計算結果可能要大于0o 75，若比0.7

11、5大 10%,那么計算結果就不可靠了。這也是對結果的影響.只有開殼層才有自旋污染，閉殼層沒有此問題.直觀上來看，造成自旋污染的原因是 混入高自旋的態，所以總自旋變大.從計算原理上說，開殼層的a，b軌道能級不一樣，會 自旋極化；若極化過多，則會產生自旋污染（呵呵，這話我也不太懂，不知道是否有錯誤 。哪位大蝦指點一下）。根據文獻報導，UMP2是自旋污染最大的。若有未成對電子，需要用PMP2 （即投影 MP2的能量）其他的，UHF，UMP3，UMP4的自旋污染比較大要避免的話，HF就選擇ROHF方法.不過此法有不少缺點，目前基本不用。第二個方法 就是采用自旋投影方法，gaussian基本使用此法。（

12、我沒用過，哪位用過說說怎么用的） 第三個方法就是半電子近似的方法，這個用于半經驗方法，象半經驗程序mopac中有.第（完整word版）自旋污染的基本概念以及對計算結果的影響 四個辦法，就是用casscf方法了，它沒有自旋污染問題。還有.。.應該還有其他辦法：）其實我沒有處理過這個問題，但是我對此也感興趣。不知我說的對不對，與大家討論討論。自旋污染看 “before annihilation” 一行MP2能量看“EUMP2行有ROMP2算法的，可能能解決自旋污染問題。LZ可以試試看.怎樣消除自旋污染Re： QC強制收斂般是沒有辦法了再強制，所用時間很長.conconver=8是默認 可以調整=5

13、降低收斂標準IOP里面有IOP（5/13=1）就是防止自選污染 還可以算下去Re：用DFT方法進行優化自旋污染會小一些。個人意見，僅供參考Re：我用的就是DFT方法。但是也挺嚴重的。因為是高自旋，沒有辦法。Re:如果是高自旋的體系是不是用MCSCF更為合理一些。CASSCF 優化，CASPT2校正能量。Re： RO方法沒有自旋污染，如RODFT（ROB3LYP，ROOLYP。.），ROMP2，.。但是可能會不收斂guass03出錯信息及解決辦法整理（待續）收斂問題的調整如果SCF計算收斂失敗，可以采用以下技巧：1。在Guess關鍵字中使用Core,Huckel或Mix選項，試驗不同的初始猜測。

14、2。對開殼層體系，嘗試收斂到同一分子的閉殼層離子，接下來用作開殼層計算的初始 猜測。添加電子可以給出更合理的虛軌道，但是作為普遍的經驗規則，陽離子比陰離子更容易 收斂。選項Guess=Read定義初始猜測從Gaussian計算生成的checkpoint文件中讀取.3。另一個初始猜測方法是首先用小基組進行計算，由前一個波函得到用于大基組計算 的初始猜測（Guess二Read自動進行）。4。嘗試能級移動（SCF=Vshift）。5。如果接近SCF但未達到，收斂標準就會放松或者忽略收斂標準。這通常用于不是在 初始猜測而是在平衡結構收斂的幾何優化。SCF二Sleazy放松收斂標準，Conver選項給出

15、更多的 控制。6。一些程序通過減小積分精度加速SCF。對于使用彌散函數，長程作用或者低能量激 發態的體系，必須使用高積分精度：SCF二NoVarAcc.7。嘗試改變結構。首先略微減小鍵長，接下來略微增加鍵長，接下來再對結構作一 點改變。8。考慮使用不同的基組。9。考慮使用不同理論級別的計算.這并不總是實用的，但除此之外，增加迭代數量總 是使得計算時間和使用更高理論級別差不多.10。關閉DIIS外推（SCF二NoDIIS）。同時進行更多的迭代（SCF=（MaxCycle二N）.11。更多的SCF迭代（SCF（MaxCycle二N），其中N是迭代數）。這很少有幫助，但值得一試。12。使用強制的收斂

16、方法。SCF二QC通常最佳，但在極少數情況下SCF=DM更快.不要忘記給 計算額外增加一千個左右的迭代.應當測試這個方法獲得的波函，保證它最小，并且正好不是穩 定點（使用Stable關鍵字）。(完整word版)自旋污染的基本概念以及對計算結果的影響13。 試著改用DIIS之外其它方法(SCF二SD或SCF=SSD)。自旋污染自旋污染是在非限制從頭計算中的誤差，通常在總能量中輕微出現。這是由于基態波函和更高 能態發生混合。高的自旋污染會影響結構和布居數分析，并極大地影響自旋密度.它還會影響收 斂，特別是MPn任務.可以通過比較總自旋，的預期值和計算值驗證自旋污染。S應當等于s(s+1)。 對于有

17、機分子，偏差不超過正確值10%的計算通常是可接受的。未成對電子數s預期S值0 0 0 1 0.5 0。75 2 1。0 2。00 3 1。5 3。75 4 2。0 6.00 5 2.5 8.75。自旋污染經常在非限制HartreeFock (UHF)計算以及非限制M ller-Plesset (UMP2, UMP3, UMP4)計算中遇到。而在DFT, CI或CC計算中幾乎見不到任何顯著的自旋污染。在包含 金屬原子的體系中和在過渡結構中，出現高自旋污染是很平常的.非限制計算經常與自旋猝滅步驟結合，它能從某些點計算的波函中刪除相當比例的自旋污 染。這使得自旋污染最小，但不能徹底防止。在Gauss

18、ian中，選項iop (5/14=2)”通知程序 使用猝滅波函產生布居數分析.可以進行自旋限制開殼層計算(ROHF, ROMP2)。它的優勢是沒有自旋污染。缺點是需要 額外的CPU時間，用于正確控制單占據軌道和雙占據軌道以及它們的相互作用。 這個數學方法 的結果是，ROHF計算給出很好的總能量和波函，但是單占據軌道能量不嚴格服從Koopman理論。 前綴RO可以用在HF，所有的DFT方法，AM1, MINDO3, MNDO,和MP2波函中。ROMP2方法尚不支 持分析梯度。另一個方法是運行非限制計算，在獲得波函(即Gaussian輸出文件中打印的PUHF和PMP2 值)之后，接著投影除去自旋污

19、染.這給出的是對能量而不是對波函的修正。自旋投影結果不會 給出使用限制開殼層計算得到的能量.這是因為優化的非限制軌道描述污染態而不是自旋投影 態。半經驗程序經常對原子團計算使用半電子方法。半電子方法是一種數學技巧，用于處理(完整word版)自旋污染的基本概念以及對計算結果的影響 RHF計算中的單占據軌道.這個結果依靠近似波函和軌道能量，與總能量一致。由于使用了單重 行列式計算，所以沒有自旋污染。總能量的一致使得它可以對單重態用RHF，對三重態用半電子 計算單重-三重裂距。半電子計算不服從Koopman定理.同樣，也不能得到自旋密度.而Mulliken 布居數分析通常還算是合理的。如果自選污染小

20、，對自旋猝滅波函和自旋投影能量更適宜使用非限制方法.不要對DFT使用 自旋投影.其中自旋污染的總和對使用限制開殼層方法更加重要。參見：Jo Baker，Ao Scheiner，J. Andzelm， Chem。Phys. Lett. 216， 380 (1993)。M. Wo Wong, Lo Radom, J。Physo Chemo 99， 8582 (1995) o H. B。Schlegel， Encyclopedia of Computational Chemistry 4，2665， John Wiley & Sons， New York (1998)使用贗勢從頭計算不能僅僅忽略核電

21、子，但它們可以被哈密頓量中的勢能項代替.這稱為核勢，有效核勢 (ECP)或相對論有效核勢(RECP)。核勢必須連同相應的價電子基組一起使用。這不僅減少了計 算時間，核勢還可以包含相對論質量虧損和對重原子近核區非常重要的自旋耦合項的影響。這 對Rb和更高的重元素是首選。有些使用贗勢的計算結果不同于全電子計算的結果，輸出的總能量是價電子能量。這個能量的大小依賴于核勢模擬多少個電子殼層，其中電子明確描述 為價電子和價電子-1殼層，計算一般更準確.同樣，維里定理的檢查是沒有意義的。怎么解決S2的問題# freq ub3lyp/gen scf=(maxcyc=120) geom二check guess二

22、read這個的計算命令但是最后算出 HF=-316。6140566S2=1.702056按道理這個S2=0.75 2重態的？想請教下這個怎么解決？這個算的東西還可靠，嗎？如何消？高斯使用手冊 沒有這方面的 提示就是什么投影方法？這個有誰用過嗎？能不能給點實例 介紹 學習！下!盼高手解決！ ！!問題是這樣的我算的這個體系作為初始步4重態能算出來而且結果很好但是2重態卻這么也算不出 來，好不容易計算個S2=1.70幾 自旋污染嚴重，我看別人文獻里面說他們說這種結構不存在？！ ! ！我是該放棄尋找還是消除自旋污染繼續找？迷茫中？望高手解決中！我覺得自選污染顯然就是四重態混雜引起的你算出來的2重態能量

23、是不是明顯比四重態高，如果是，似乎沒有什么必要進一步研究了現在是2重態能量 比四重態高，問題就是四重態于2重態在這里會有交叉這個時候怎么辦？用QC IOP命令強制算出來的東西 都是自旋污染的（S2=1。7）的 這個時候我怎么辦？掃描這方法我也想了下 但是覺得不實際 而且太費時我該怎么辦？ 這個點算出來 該不改放棄？用 rob3lyp可能不如rdft/udft收斂快，甚至根本不收斂方法已定 這個時候 有什么辦法沒找這個這個點？或者跟別人樣說這個點不存在換方法，用bp86換方法能解決自旋污染？如果是自選污染 我可不可以認為這里有自旋交叉？放棄尋找這個點？消除自旋污染只能用多參考態方法Density

24、 Functional TheoryIn density functional theory (DFT) the energy of a system is given as a sum of six components:DFT NNTv coulexch corrtheory。The definitions for the nuclearnuclear repulsion 知，the nuclearelectron attraction E, and the classical electron一electron Coulomb repulsion E l energies are the

25、 same as those used in HartreeFock theory o The kinetic energy of the electrons Et as well as the non-classical electron-electron exchange energies E h are, however, different from those used in HartreeFock theoryo The last term E describes the correlated movement of electrons of different spin and

26、is not accounted for in Hartree-Fock theory. Due to these differences, the exchange energies calculated exactly in Hartreetheory。in DFT methods. or the electronVarious approaches exist to calculate the exchange and correlation energy terms These approaches differ in using either only the electron de

27、nsity ( local methods) density as well as its gradients (gradient corrected methods or generalized gradient approximation, GGA). Aside from these pure DFT methods, another group of hybrid functionals exists, in which mixtures of DFT and Hartree-Fock exchange energies are used.in DFT methods. or the

28、electronLocal methodswith the localApproximation (LDA) methodoalso referred to as the Local Spin Density Approximation (LSDA)o with either the LSDA or SVWN keyword 。 with the localApproximation (LDA) methodoalso referred to as the Local Spin Density Approximation (LSDA)o with either the LSDA or SVWN

29、 keyword 。 A sample input for theFor open shell systems, using unrestricted wavefunctions, this is This method is used in Gaussian frequency calculation for themethanol-1yl radical using the SVWN functional and the 6-31G (d) basis set is as follows:#SVWN/631G(d)freq#SVWN/631G(d)freqSVWN/631G (d)耳 3S

30、VWN/631G (d)耳 3H5h*0 2C102 1 r2H3 2 r3 1 a3H4 1 r4 2 a4 3 d4H5 1 r5 2 a5 4 d5r2=1.35454381 r3=0。 97675734r4=1。 09713107r5=1。 09208721a3=108.95266528a4=119.30335589a5=113。 0329959d4=-26。 58949552d5=-149。 77421707Please observe that the acrynom VWN is used to indicate different functionals in differen

31、t software packages. While in Gaussianthis keyword refers to functional III described in the VWNpaper (S. Jo Vosko， Lo Wilk，M。Nusair， Can。J。Phys。1980, 58，1200)， most other packages refer to functional V in the same paper. This latter functional can also be used inGaussian with the VWN5 keywordoGradi

32、ent corrected methodsThe gradient一corrected exchange functionals available in Gaussian include (given with their abbreviations in parenthesis):Becke88 (B)Perdew-Wang (PW91)Modified Perdew-Wang by Barone and Adamo (MPW)Gill96 (G96)Perdew-Burke-Ernzerhof (PBE)and the gradientcorrected correlation func

33、tionals areLYP by Lee, Yang, and Parr (LYP)Perdew-Wang (PW91)Perdew 86 (P86)Becke96 (B96)Perdew-Burke-Ernzerhof (PBE)In all cases， the names of these functionals refer to their respective authors and the year of publication。 All combinations of exchange and correlation functionals are possible， the

34、keywords being composed of the acronyms for the two functionals。 The frequently used BLYP method, for example, combines Beckes 1988 exchange functional with the correlation functional by Lee, Yang, and Parr. An input file for the BLYP frequency calculation of the methanol-1-yl radical is:#P BLYP/6-3

35、1G(d) freqBLYP/6-31G (d)opt min methanol radicalBLYP/6-31G (d)opt min methanol radicalC1021r2H32r3 1 a3H41r4 2 a4 3 d4H51r5 2 a5 4 d5r2=1.38575119r3=0.98034945r4=1.09658704r5=1.09123991a3=108。 13230593a4=118。 45265645a5=112。 14316572d4=-29。 40173022d5=145.83450795Another frequently used GGA function

36、al is BP86 composed ofhe Becke 1988 exchange functional and the Perdew 86 correlation functional. The PW91 functional combines exchange and correlation functionals developed by the same authors in 1991. The keywords used Gaussian for a particular GGA functional are combinations of the acronyms for e

37、xchange and correlation functionals. Use of the PW91 exchange and correlation functionals thus requires the keyword PW91PW91oHybrid functionalsThe basic idea behind the hybrid functionals is to mix exchange energies calculated in an exact (Hartree-Focklike) manner with those obtained from DFT method

38、s in order to improve performance.Frequently used methods are: Becke-Halfand-HalfLYP (BHandHLYP) uses a 1: 1 mixture of DFT and exact exchange energies:XCE = 0.5大E (HF) + Go 5*E (B88) + E (LYP)XCBecke一3-LYP (B3LYP) uses a different mixing scheme involving three mixing parameters Exc = 0.2*Ex(HF) + G

39、o 8*Ex(LSDA) + 0。72大DEx (B88) + 0。81大EC(LYP) + 0。19大EC (VWN)In this latter case, the B88 gradient correction to the local LSDA exchange energies carries its own scaling factor of 0。72。 The three scaling factors have been derived through fitting the parameters to a set of thermochemical data (G1 set)

40、。PBE0 uses a 1:3 mixture of DFT and exact exchange energies: Exc = 0.25*Ex (HF) + 0.75大Ex (PBE) + EC (mPW91)B98 (Becke98) is based on a 10-parameter equation scaling components of exact exchange, GGA exchange and GGA correlation。 All parameters have been optimized simultaneously to fit thermochemica

41、 data collected in the extended G2 data set:Exc = 0.2198*Ex (HF) + Ex + ECThe Becke-halfandhalf-LYP method is used with the BHandHLYP keyword, the very popular Becke3-LYP method is used with the B3LYP keyword, and the PBE0 method with the PBE1PBE keyword。 A sample input file for the B3LYP frequency

42、calculation on the methanol-1-yl radical is:#P B3LYP/631G(d) freqBecke3LYP/631G(d) opt min methanol radicalC1r202 r2=1.37007444 r3=0.96903862 r4=1.08886325r2H3r3a3.coH4r4a4 3 d4H5r5H3r3a3.coH4r4a4 3 d4H5r5a5 4 d5r5=1。 08381987 a3=108.87657926a4=118。 49962766a5=112.61039767d4=29。19863769d5=-146。 7345

43、0089It is important to note that implementations of the B3LYP and BHLYP functionals in different programs (Gaussiar, Turbomole, Jaguar, Q-Chem) vary somewhat. Other hybrid functionals available inGaussian are： B3P86 uses the P86 correlation functional instead of LYP, but retains the three parameters

44、 derived for B3LYP; B3PW91 uses the PW91 correlation functional instead of LYP, but retains the three parameters derived for B3LYP; B1B96, a one parameter functional developed by Becke using the Becke 96 correlation functional; MPW1PW91 developed by Barone and Adamo using a modified version of the P

45、W91 exchange functional in combination with the original PW91 correlation functional and a mixing ratio of exact and DFT exchange of 0。 25 : 0。 75。Userdefined modelsCalculation of reaction and activation energies in either radical reactions or transition metal mediated reactions often shows that the

46、 mixing ratio of H and DFT-exchange is not optimal for a given purpose. The definition of user-defined models tailored to a certain purpose is therefore sometimes necessary o This can be achieved by first choosing gradient corrected exchange and correlation functionals from the list of available cho

47、ices and then specifying the mixing ratios of the various exchange and correlation energy terms. For this latter step Gaussianuses the general formula:E = P E hf + P (P E Slater + PdE non-local) + PE local + PdE non-local XC 2 X14 X3 X6 C5 CThe values for the parameters P1 - P6 are provided through

48、the IOP switches IOP(5/45二xxxxyyyy), IOP (5/46二xxxxyyyy), and IOP(5/47=xxxxyyyy) with P 二 xxxx/1000 from IOP (5/45) etc. The use of this facility will be illustrated with the MPW1K functional by Truhlar and coworkers (B。 J。 Lynch, Po L. Fast, Mo Harris, D. G. Truhlar, J。Phys. Chem. A 2000, 104, 4811

49、) developed to optimize reaction and activation energies of free radical reactions This method is very similar to the MPW1PW91 hybrid functional, but uses a mixing ratio of exact and DFT exchange energy terms of 0.428 :0 o 572. An input file for the methanol-1 一yl radical frequency calculation at th

50、eMPW1K/631G(d) level is:#P MPWPW91/631G(d) freqIOP(5/45=10000428) IOP(5/46=05720572) IOP(5/47=10001000)MPW1K/6-31G(d)opt min methanol radical(5/47=10001000)MPW1K/6-31G(d)opt min methanol radical0 2C1 TOC o 1-5 h z 021r2a3a4 3 d4a3a4 3 d4a5 4 d5H41r42H51r52 r2=1.35355619r3=0。 95690342r4=1.08170148r5=

51、1.07685321a3=109。 38001654a4=118.56662303a5=112.99211737d4=-29.64928765d5=-147.42512309The gradient corrected MPWPW91 functional is chosen here as the basis of the hybrid method and the amount of exact exchange is specified through parameter P2O The composition of the functional is also described in

52、 the output file in a compact format:IExCor= 908 DFT=T Ex=mPW+HF Corr=PW91 ScaHFX= 0。 4280ScaDFX= 0o 5720 0.5720 1。 0000 1.0000The value given after ScaHFX corresponds to P while the four values listed after ScaDFX are P3-P6. Please observe that theses lines are only printed when using the extended

53、output option of Gaussiar。 Before using this facility with a particular version of Gaussiar, please make sure that your IOPdirectives are read correctly。 In some versions ofGaussian, the options P 一 P are set with IOP(3/76) IOP(3/78) instead of IOP (5/45) IOP(5/47). In this latter case, the values o

54、f P1 - P6 are given in five digit format o The keyword line for the PBE0 frequency calculation thus becomes:#P PBEPBE/6-31G(d) freqIOP(3/76=1000002500) IOP(3/77=0750007500) IOP (3/78=1000010000)Absolute and relative CPUtimes for DFT- and HartreFockcalculations depend dramatically on the size of the

55、system under investigation, the computer hardware and the operating system usedo The following table contains CPU times ( in seconds) for all frequency calculations on the methanol1yl radical described above with Gaussian03, Rev。 B。03. The platforms used are the CIP-room Pentium 4/2 GHz computers un

56、der LINUX. In all cases the frequency calculation has been performed on the structure optimized at the given theoretical levelo The 631G (d) basis set has been used in all cases.methodCPU (P4/2GHz)secvib(1)cm1vib (2) cm1vib(3) cm1 USVWN/6-31G (d)31369014551036UBLYP/631G (d)69359914631042UB3LYP/631G

57、(d)71376115081072UMPW1K/6-31G (d)73395015521102UHF/631G(d)11412316271156obsvd.-365014591056Despite all platformdependend variations, it is clear that hybrid DFT calculations are similarly laborious as GGA calculations , while local DFT methods are somewhat and HF-calculations much cheaper o In match

58、ing calculated (harmonic) gas phase vibrational frequencies with experimentally observed anharmonic values measured in a matrix (Ro D。 Johnson III, J. W。 Hudjens, J. Phys Chemo , 1996,100, 19874), please remember that uniform scaling is often performed to improvecorrespondence (see our earlier discu

59、ssion of this point) o The values reported here are unscaledo A general trend visible in this example is the increase of vibrational frequency wavenumbers with an increase in the amount of HartreeFock exchange contributions (varying between 0 and 100% from the second entry to the bottom in Table 1).

60、Further technical considerationsEvaluation of the exchangecorrelation energies of all density functional methods implemented in Gaussian involves a gridbased numerical integration stepo The computational effort required for this step strongly depends on the selected grid size. The larger the number