量子化学是以量子力学为理论基础

1、量子化学综述量子化学简介量子化学是以量子力学为理论基础， 以计算机为主要计算工具来研究物质的微观结构与宏观性质的关系科学， 用以解释物质和化学反应所具有的特性的内在本质及其规律性。1926年，薛定诗成功地解决了量子态也（r , t） 是如何随时 间演化及各种情况下求出波函数的问题，提出了著名的薛定谔方程。 1927 年， 化学家heitler和london等人成功地利用量子力学理论解释了h分子的形成，开辟了用量子力学方法研究分子中电子行为的广阔领域，标志着量子化学的开 始。 随着量子化学基础理论和计算方法的发展， 不但使其成为解释化学现象微观 本质的强有力工具，而且，使通过量子化学计算来预测化

2、合物性能成为可能。量子化学是理论化学的一个分支学科， 是将量子力学的原理和方法应用到化学问题研究中而产生的一门学科， 经过化学家们的努力， 量子化学理论和计算方法在近几十年来取得了很大的发展， 在定性和定量地阐明许多分子、 原子和电子尺度级问题上已经受到足够的重视。 目前， 量子化学已被广泛应用于化学的各个分支以及生物、 医药、 材料、 环境、 能源、 军事等领域， 取得了丰富的理论成果， 并对实际工作起到了很好的指导作用。量子化学可分基础研究和应用研究两大类， 基础研究主要是寻求量子化学中的自身规律， 应用研究是利用量子化学方法处理化学问题， 用量子化学的结果解 释化学现象。二 现代价键理论

3、计算方法1 价键波函数的计算在价键方法中一个重要的性质是使用了非正交轨道，所有的n!项都对矩阵元有贡献，至今为止仍然没有高效的算法来计算hamiltonian 和重叠矩阵元，这就是价键理论中著名的"n!"困难。计算hamiltonian矩阵元和重叠矩阵元有两种 方法:一种是经典的行列式展开方法，另一种则是对不变式方法。a经典的行列式展开一个hlsp波函数可以表达2n个slater行列式的线性组合："仃(1一k 0(5).i一n为共价键的个数，d(qk)为d(qk)的slater行列式，pi是第i个共价键中交换该 键的两个成键电子的算符。hamiltonian矩阵元

4、可以表达为：rtsr<u <f其中frs和grs,ut为单电子和双电子积分，d(srs)和d(sstru)为两个价键行列式重叠 矩阵的一级和二级余子式。b.无自旋价键理论的对不变式方法对不变式(pairedpermaneni) 1定义为:ppd")= 2需(尸)/与5_/_ 产”其中a=aj ,i,j=1,2,.n(n为偶数)是n*n的方阵，为了简单起见，这里只考虑自旋为0的情况。指标i=2k-l(2k)的配对指标定义为2k(2k-l),记为i。 则对不变式为：ppd(4)= e4,ppd4；ppd(/)k<i其中当l=工时，dkl=1;当/.时，dkl=1/2,p

5、pd看为2阶的对不变式，为从a中去掉第k, l行和第i , i列得到的(n-2)阶对不变式，其列指标保持原顺序，行指标为：如果l=不，原顺序不变；如果1手支:土，配对，其余不变。注意到价键波函数k是一个对不变式，可以表示为：m二岫,n)两个价键波函数k和l的重叠矩阵可以表示为对不变式形式ml ppd（s）其中s是价键轨道重叠矩阵。哈密顿矩阵元可以表达为：=££内琛ppd&s煞）+ ££工力/“需ppq&6微） / k<ji<-j k </其中：f； = fiks-f +力为+儿与+玉& + g存/+益我哈= &#

6、163;n2所叩2价键理论的从头算计算方法hamiltonian矩阵元和重叠矩阵元求得之后，便可通过求得本征值方程计算 得到能量，这个本征方程的维数与价键结构的数目有关，而上述两个矩阵元和价键轨道系数有关。由此便有两类价键计算方法 ：第一类是属于价键自恰场方法通 过当前价键波函数获得hamiltonian矩阵元和重叠矩阵元，算出能量后在这基础 上进行变分，从而使价键波函数和能量同时得到优化， 这类方法有vbsc和bovb 等，它们主要考虑静态电子相关（bovb考虑了部分电子相关）；第二类方法不属 于自恰场方法，通过产生激发价键结构，扩大本征方程的维数来达到改善能量计 算结果的目的。这类方法有

7、vbci等，主要考虑动态电子相关。这些从头算价键 方法都采用非正交价键轨道并基于 hlsp函数和无自旋价键方法。balint-kurti 和van lenthe 等人建立了价键自洽场（vbscf方法，vbscf 方法是一种非正交的多组态（价键结构）scf方法。vbscf介键波函数甲是一系列 价键结构的线性组合，对价键结构的线性组合系数和价键轨道同时优化，使得体系的总能量最低。vbscft法由于将体系波函数表示为价键结构的线性组合，因 此，vbscff法描述了体系的静态电子相关效应。但是，由于其没有包括动态相 关效应，vbscft法的定量计算结果不能让人满意。hiberty等提出了 bovbj法

8、2。bovbj法允许不同价键结构的价键轨道不完全相同。由于各个轨道可以根据各自所处价键结构的特点进行调节，bovb方法包括了部分动态相关能，而这一部分相关能是描述化学键形成和断裂所需要 的。然而由于bovbj法允许在vbscft法中相同的轨道不完全相同，这些轨道 只有微小的差别，因此，在实际应用中，bovb方法存在收敛困难等问题。从本质来讲，bovb1扩展了的vbscfj法。随着计算机技术的发展，从头价键计算的计算将成为可能。然而，由于vbscf 方法的计算结果并不能令人满意，因此，发展高精度的价键理论计算方法将是今 后价键理论研究的一个重要课题。近年来分子轨道理论的多重参考技术己成为高 精度

9、量子化学计算必需的技术，定域轨道的应用也己成为相关能计算的有效途 径，而这两点是价键理论方法所具有的特点。可以相信，随着计算方法的不断完善，价键方法将改变目前这种只作为定性讨论问题的工具，成为分子轨道理论方法之外的另一个精确量子化学计算方法。3无自旋价键计算程序xmvbxmvb是一个无自旋价键从头计算价键程序。它建立在无自旋价键方法以及对不变行 列式算法基础的从头算价键计算程序。可以利用这个程序包来使用任何形式的vb轨道做任何类型的vb计算。它可以实现上述 vbscf、bovb、和vbci等计算，也可以将价键方法 和一些分子轨道方法联合，比如vbmpz , vb-dft等等。这里简单介绍一下这

10、个程序包的算法和功能。a.矩阵元的计算根据不同的要求可采用行列式或者对不变式来计算矩阵元。b.价键轨道的优化在vbsciw bovb勺计算中，价键轨道是需要优化的。xiamen1序中价键轨 道可以自由定义和优化，可以根据研究的需要使用重叠加强基组(overiapenhaneed orbital , oeos)纯定域的杂化原子轨道 (hybrid atom orbital , haos和成键变形基组(bdos)。xmv则序采用了两种轨道优化的策略，第一种 是基于广义布里渊定理 (generalized brillouln theorem):优|方-耳力)=0 其中vij是将基态的占据轨道i替换成

11、未占据虚轨道j而后形成的单重态激发态。 在vbscf计算中虚轨道是任意选择的，可以通过投影的方法来得到。另外一个策略就是利用数值方法计算能量梯度。向前差分近似计算中能量梯 度向量可以表达为：3e e- e3ae其中勺2是指在能量表达式中系数cij替换成cij +所得到的能量值。程序 也可以使用中心差商的方法来计算梯度。梯度计算出来以后，程序使用 dfp(davidon-fleteher-powell)方法来优化价键轨道。三量子化学的应用80年代以来，随着计算机和大型工作站的出现以及有关量子化学计算软件的 大量推出(例如gsmess系歹1，hondo系列j, gaussian系列等程序包的推出)

12、， 使量子化学计算方法成为化学计算理论研究的主流。量子化学计算方法可用于解 释或预测原子，分子和晶体的各种性质，包括分子的能量和结构，化学反应路径， 过渡态的能量和结构，分子轨道，原子电荷，成健和化学反应能量，极化率和超 极化率，各向异性超精细耦合常数，静电势以及电子密度等。其应用范围也不断 扩大，涉及到化学各分支学科，并取得了可喜的成果。1在基础化学研究中的应用量子化学在化学反应机理与动力学研究，分子结构与价键分析等基础研究中 的应用一直是化学工作者研究的热点。在研究反应机理，动力学，电子结构，能 带分析以及价键分析等方面应用较多。junjing gu等网人分别用价键从头计算法，vbscf以

13、及vbci等方法计算研究 多烯类抗菌素c2nh2n+2 (n=2 8)和多烯类自由基c2n-1h2n+1 (n=28)的基态和激发 态的相对能量和电子结构。价键从头计算法算得的能量值与现有的实验值一致， 同时与mrci, caspt3以及以前开发的价键半经验法(vbdft(s)所算的结果一致。他们又在价键理论基础上发展了一种新的价键半经验法(vbdft(s), 在化学价键结构方面， 用以研究讨论多烯类聚合物的共价激发态的电子结构， 所算的结 果与所有的从头计算法所得结果一致。2 在有机金属材料中的应用金属有机材料， 是一种典型的纳米多孔材料， 在化学吸附和化学分离技术等方面比传统的纳米材料存在

14、许多潜在的优点。量子化学计算方法应用于这一领域，用以研究mofs 的结构和电子性质4 。近年来引起了的广泛注意，得到了迅速的发展。d. f. bahr和c. a. bauer等根据密度泛函理论(dft)对立方竣酸锌单晶 的弹性性能进行计算，分别采用了局域密度近似(lda) 计算方法和广义梯度近似(gga)计算方法。结果算的密度泛函理论值(lda-gga平均值)为e = 21.6 ±0.3 gpa,cii = 0.28 仕01 gpa, ci2= 0.11 0.01 gpa以及c44= 0.03 0.02 gpa lda 计算结果证实了较小值的模量c44预示着irmof-1结构的不稳定

15、，同时，gga预 测出较大的晶格常数，将导致其较小的弹性常数。将密度泛函理论( dft )计算 出的弹性系数与通过纳米压膜测量的金属有机材料的弹性塑性响应结果相比较， 杨氏模值较密度泛函理论(dft )计算值偏低，同时各向异性弹性常数也较dft计算值偏低。samant铸6使用dft-lda和更为精确的paw方法进行计算，得到显著不同 的结果，尤其是他们发现了一个小值的模量c44,预示着这种mofs结构的不稳定 性。其后，zhou和yildirim等4人报道了使用dft-lda计算得出的irmof-1的弹 性常数值及其通过计算与测量得出的振动特性。他们证实了samant繇发现的小值模量c44o他

16、们的测量和计算中子非弹性散射结果基本上支持了通过dft-lda计算得出的弹性性能。 计算得出的点阵参数也十分接近通过中子衍射实验所得出 的测量值。3 在生物大分子体系研究中的应用生物大分子体系的量子化学计算一直是一个具有挑战性的研究领域， 尤其是生物大分子体系的理论研究具有重要意义。电子水平上对体系进行精细的理论研究，是其它理论研究方法所难以替代的。因此要深入理解有关酶的催化作用、基因的复制与突变、药物与受体之间的识别与结合过程 及作用方式等，都很有必要运用量子化学的方法对这些生物大分子体系进行研 究。毫无疑问，这种研究可以帮助人们有目的地调控酶的催化作用，甚至可以有目的地修饰酶的结构、设计并

17、合成人工酶；可以揭示遗传与变异的奥秘，进而调控基因的复制与突变，使之造福于人类；可以根据药物与受体的结合过程和作 用特点设计高效低毒的新药等等，可见运用量子化学的手段来研究生命现象是十 分有意义的。4在锂离子电池研究中的应用锂离子二次电池因为具有电容量大、工作电压高、循环寿命长、安全可靠、无记忆效应、重量轻等优点，被人们称之为“最有前途的化学电源”，被广泛应用于便携式电器等小型设备，并已开始向电动汽车、军用潜水艇、飞机、航空等 领域发展。锂离子电池又称摇椅型电池，电池的工作过程实际上是 li离子在正负两电 极之间来回嵌入和脱嵌的过程。因此，深入锂的嵌入 -脱嵌机理对进一步改善锂 离子电池的性能

18、至关重要。ago等用半经验分子轨道法以。2 h14作为模型碳结构研究了锂原子在碳层问 的插入反应。认为锂最有可能掺杂在碳环中心的上方位置。ago等 用abinitio分子轨道法对掺锂的芳香族碳化合物的研究表明，随着锂含量的增加，锂的离子 性减少，预示在较高的掺锂状态下有可能存在一种li - c和具有共价性的li -li的混合物。satoru等用分子轨道计算法，对低结晶度的炭素材料的掺锂反应 进行了研究，研究表明，锂优先插入到石墨层间反应，然后掺杂在石墨层中不同 部位里。随着人们对材料晶体结构的进一步认识和计算机水平的更高发展，相信量子化学原理在锂离子电池中的应用领域会更广泛、更深入、更具指导性

19、。5在非线性光学材料中的应用近年来，非线性光学材料（nlo）的研究工作获得了迅猛发展，在很多方面得 到了的应用，但对于其结构与性质还需进一步研究。paresh c. ra野7对推拉系列卟啉类非线性光学材料进行了研究， 对中心金属离子对其第一极化率和双光子吸收(tpa)的影响进行了量子化学分析。通过 b3lyp/6-31g(d,p)水平优化，获得 了该口卜咻分子的几何构型；采用zindo/cv方法对该吓咻和tpa的物化性质进行 了分析计算，结果发现金属离子d轨道和吓咻大环冗轨道形成d-碘辄体系，对该金 属吓咻的nlo、tpa性质有很重要的作用。四 结论随着科学技术以及软件技术的不断发展和开发，

20、量子化学的计算方法也在不断更新， 其所应用的领域也越来越广。 量子化学计算的研究越来越多地与实验相结合， 其研究方法和手段呈现多样化。 量子化学计算的发展给科研提供了许多可以借鉴的经验和方法，与之相互结合、取长补短，可使科研走向一个新的高度。参考文献1 wu, w.; wu, a.; mo, y e. a. int. j. quantum chem. 1998, 67, 287.2 hiberty, p. c.; humbel, s.; van lenthe, j. h.; byrman, c. p. j. chem. phys. 1994,101, 5969.3 junjing gu, yonghui lin et al. covalent excit