耦合簇理论的方法综述

耦合簇理论的方法综述1耦合簇理论的方法综述厦门大学化学系量化专业苏培锋耦合簇方法于上世纪六十年代由Cízek和Paldus【1】首先引入在量子化学研究中，它很快成为最可靠，最准确的量子化学计算方法之一。经过几十年发展，这种方法日趋成熟，并正在蓬勃发展，它从单参考态耦合簇方法（SRCC）开始，发展到多重参考态耦合簇方法（MRCC），SRCC只能解决闭壳层的多电子体系，而MRCC已经能够计算开壳层，高度近似简并的多电子体系，不仅如此，它还与其他方法互相结合，渗透，产生新的计算方法，比如价键与耦合簇的结合（VBCC）。当然耦合簇目前也存在一些困难和问题有待解决（后文将提及），但是瑕不掩瑜，应该说它是一种比较好而且有着较大发展潜力的理论方法。一、耦合簇方法的发展过程作为一个方法，耦合簇较晚才被理论化学所接受，这是因为早期的发展者采用了优美但却又深奥的数学方法，比如Feynman图、二次量子化等[1-4]。但是耦合簇方法后来所表现出来的优点很快使它成为量子化学家的至爱。首先发展的是单参考态耦合簇方法[1-4]，作为计算电子相关的方法，它比传统的非完全CI优越。必须指出的是完全的CC与完全的CI是等价的，与CI相似，实际计算的CC也采用了非完全的形式（CCD、CCSD、CCSDT、CCSDTQ），本文准备着重介绍CCSD。由于SRCC不可克服的弱点（只能用于闭壳层体系），后来发展了多参考态耦合簇方法，一般公认的MRCC有三种：1.FOCKSPACEMRCC[5-9]、2HILBERTSPACEMRCC[10-13]、3STATESSELECTEDMRCC[14-16]，由于篇幅所限，只作一般介绍。二、单参考态耦合簇方法（SRCC）1.基本要素及其二次量子化形式SRCC方法的基本要素有：参考态（波函数），指数簇算符，哈密顿的新形式HN。SRCC是一种基于以下拟设（anastz）的方法：一个确定的多粒子波函数可以写成指数簇算符作用于一个独立粒子波函数上[1-4]：-（1）这里的波函数和指数簇T都采用二次量子化形式，其中a+是产生算符；a是湮灭算符上式是指在有相互作用的多粒子体系中，N个粒子体系的状态波函数可用粒子在单粒子态上的占据数来描述，而单粒子态占据又可用一系列产生算符作用在真空态|上来得到。耦合簇理论的方法综述2一般的，我们有：（2）T1，T2分别对应单电子相关，双电子相关；Tn则是多电子相关，它可以分为连接相关图簇（linkedcluster）和非连接相关簇（unlikedcluster），前者表示多电子直接相关，后者表示同时分别在不同区域发生的较小电子簇相关。比如四电子相关，它可以是真正四个电子相关，也可以是两对双电子相关，还可以是单电子、三电子相关，等等多种形式。我们把（2）带入（1），并展开，与CI的展开形式作比较发现耦合簇算符和CI展开系数有以下关系：C1=T1C2=T2+(1/2)T12C3=T3+T1T2+(1/3!)T13………………………………从以上的关系式就可以看出，对应于CI的三激发，耦合簇用了三项来描述，考虑了非连接相关簇，这使得它具有大小一致性，这是它比CI优越的地方之一。要解方程就必须提到哈密顿算符的二次量子化形式：对产生湮灭算符的乘积采用WICK定理，将编时乘积改为正规乘积，可得新的哈密顿算符的形式如下：H＝（3）仔细分析，可以发现，该式子的第一二项是Fock算符，最后两项是HF的能量EHF经改写，（11）可有如下形式：（4）或（5）耦合簇理论的方法综述3这里定义一个新的算符：HN这个算符将对以后的计算有极大的帮助。2.CCSD的能量方程、系数方程虽然原则上所有电子相关簇都有存在的可能性，但从两点考虑，双电子相关是最重要的。首先哈密顿算符只包含双电子算符，电子间相互作用可分解为各对电子作用之和；其次，pauli原理指出三电子以上直接相关簇可以忽略。故本文只考虑CC的近似：CCSD。在CCSD中，求能量用能量方程，求波函数用系数方程（求得系数，则可知簇算符的形式，从而波函数可知）。本文侧重于方程是如何得到的，及其引入过程，方程的解是用场论方法得到的，在此不祥述。我们注意到：exp（T）可以做级数展开：（6）若Tn取到T2，即只考虑到二级激发，采取这种近似后，即可得到CCSD中的参考态波函数。将（6）式带入（1）式，经简化，有：（7）得到新的波函数形式，我们知道量化计算的关键是解以下形式的薛定厄方程：（8）把得到的波函数形式带入薛定厄方程中，有（9）把它投影到零级波函数上，得到关于能量的表达式：（10）注意这里用了一个假设：中间归一化。同样的，若把（6）式投影到n级激发态波函数上，可以得到：（11）实际上，如果将（6）式同时左乘exp（-T），再投影到零级n级波函数上，有（12）耦合簇理论的方法综述4（13）以后可以看到，9,10两式比7,8两式更为有用。这是因为可以用“TheHausdorffExpansion”来展开，此式是以后计算能量和t系数的基础：=（14)在CCSD中，该式的右边简化为：（15）在（15）式中，将HN＝FN＋VN，并将FN、VN、、T1、T2的二次量子化形式带入式中的对易式中，同时考虑产生湮灭算符的WICK定理，（14）式可再次简化为：＝（16）该式下缀C表示connected，这是说在展开式中由于对易而产生的众多代数项之中，只有HN与至少每个簇算符有一个连接（wick定理中称为缩并）的这些项才不为零。这就大大减少了项的数目。有文献也称之为“linkedclustertheorm”。将（12），（14），（16）结合起来，我们就得到CCSD的能量方程：（17）将的具体形式（16）式带入，解得耦合簇理论的方法综述5（18）结合（13）（14）（16）式，可得求CCSD的t的系数方程：（19）（20）t系数的求解方程及其麻烦，在此只列出T1的，T2省略。=0除了代数方法，也可用图形方法求解，具体可见文献。对得出的能量，系数方程，采取一系列近似和计算数学的处理，加以程序化，就可以用计算机求解了。从上面的解析式可知，CCSD数学形式复杂，求解繁琐，而且方程是非线性的，所以它必然要消耗较大的计算资源，一般的，它只能计算小分子。本节结尾总结一下SRCCSD与CISD的比较：从形式上看，耦合簇展开与CI展开无太大区别，都是用激发态组态函数展开，但两种展开的思想是不同的，CI展开只是把组态函数作为基组，简单将其按激发等级分类，故而出现大小不一致的问题；而耦合簇展开是以连接相关簇为基组，并考虑了不同相关簇的作用。所以很好解决大小不一致的问题。同时作CI展开时，最多只能展到四电子激发，再高激发就被忽略，然而高激发项系数有低激发项的贡献，CI将其忽略了，不仅产生较大误差，而且不易收敛，而CC考虑的多体算符里有不相连相关簇，极大减少了误差。三、SRCC的缺点和多参考态耦合簇MRCC的发展SRCC成功解决了闭壳层非简并态，和大小一致性问题，但对于开壳层和高度近似简并的情形的处理是失败的，这是因为在这种情况下CC方法里用的参考态不能再用单一的波函数来描述；必须采用多重参考态，也就是多个单一波函数的叠加。一般的，根据参考态波函数的展开方式，可以分为三种MRCC1．FOCKSPACEMRCC[5-9]这种方法考虑了所有可能产生的离子，即多电子体系的所有价电子的激发，它的耦合簇算符比较简单。它可以较好的计算离子化能和电子亲和力，但对势能面的描述比较差。2．HILBERTSPACEMRCC[10-13]这种方法中，每一个参考态对应一个簇算符，单个簇算14符作用于单个参考态上，叠耦合簇理论的方法综述6加后组成多重参考态。它在几何构型转变时对基态，激发态描述比较精确，而且与MRCI有紧密联系，可以利用MRCI的一些结论和校正，但它需要计算大量波函数系数，在数学处理上极为繁琐。这以上两种方法统称genuineOSCC他们的共同点是基于同样的BLOCH方程，并一样定义了作用于选定的模型空间的有效哈密顿算符。所以它们的式子的基本形式是一样的。这两种方法在以下四点上存在困难：1.模型空间（参考态）的选择2.对不同状态的截取有多种可能3.易于产生侵入态（intruderstates）4.数学解析式太繁琐，不易计算。以上诸多缺点，给了SSMRCC发展的空间和“土壤”。3.STATES-SELECTMRCC[14-16]从数学形式上看SSMRCC不够严密，但在实际计算上要比前两种方法来得简单。SSMRCC的主要是基于以下思想：在许多MRCC情形中，大多数高次项都不是很重要，只有一些是关键的。所以我们可以只选择最重要的簇来描述、组成多重参考态的波函数。说不准确点，它是利用两步激发，第一步用内算符作用于一个选定的参考态上，产生多重参考态，再用外算符作用在多重参考态上来描述我们需要的解析式。SSMRCC利用了SRCC的算符表达式，所以数学形式较为简单。实际的计算也较可行。四、耦合簇方法今后的发展与改进尽管耦合簇方法在量化计算中取得巨大成功，但必须指出，它还远远不能满足实际研究的需要，今后研究者认为它的发展方向有：1.大分子体系[17-18]对于CCSD而言，它计算所需的时间和资源与基函数数目的七次方成正比，这使得它计算大分子几乎是不可能的，多年来许多研究者在这方面作出很大的努力。2.与其他方法的结合[19-20]比如它与价键方法、组态相互作用、多体微绕的相互结合，极大拓展耦合簇方法的研究范围和适用性。3.高精度计算[21-23]众所周知，耦合簇方法的计算精度是较高的，在计算解离能时用大基组其精确度可以达到“化学精度”（±1kal/mol）水平。近年来一些研究者提出能否使其精度提高到“光谱水平”（±1cm-1）并为此做了大量工作。耦合簇理论的方法综述7参考文献：1.J.Cízek,J.Chem.Phys.,45,4256(1966).2.J.Cízek,Adv.Chem.Phys.,14,35(19693.CízekandJ.Paldus,Int.J.QuantumChem.,5,359(1971)4.A.C.Hurley,ElectronCorrelationinSmallMolecules,AcademicPress,London,1976.5.A.HaqueandD.Mukherjee.J.Chem.Phys.80,5058（1984）6.S.KochandD.Mukherjee,Chem.Phys.Lett,145,321(1988)7.K.B.GhoseandS.Pal,J.Chem.Phys,97,3863(1988)8.W.Kutzelnigg,J.Chem.Phys80,822(1984)9.K.Jankowski,J.Paldus,I.Grabowski,J.Chem.Phys97,7600(1992)10.B.JeziorskiandH.J.Monkhorst,Phys.ReV.A24,1668(1981)11B,JeziorskiandJ.PaldusJ.Chem.Phys885673(1988)12W.D.Laidig,P.SaxeandR.J.Bartiett,J.Chem.Phys86887(1987)13W.D.Laidig,andR.J.Bartiett,Chem.Phys.Lett104,424(1984)14.K.B.Ghose,P.Piecuch,S.Pal,L.Adamowicz,J.Chem.Phys10417(1996)15U.SMahapatra,B.Datta,andD.Mukherjee.J.Chem.Phys.1106171(1999)16P.Piecuch,N.Oliphant,andL.adamowicz,J.Chem.Phys.993(1993)17P.Pulay,Chem.Phys.Lett.,100,151(1983)18S.SæbøandP.Pulay,Chem.Phys.Lett.,113,13(1985)19J.Paldus,J.PlanellesTheoCh