量子化学计算方法.docx
- 文档编号:7719888
- 上传时间:2023-01-26
- 格式:DOCX
- 页数:9
- 大小:67.53KB
量子化学计算方法.docx
《量子化学计算方法.docx》由会员分享,可在线阅读,更多相关《量子化学计算方法.docx(9页珍藏版)》请在冰豆网上搜索。
量子化学计算方法
量子化学计算方法-正文
主要分为:
①分子轨道法(简称MO法,见分子轨道理论);②价键法(简称VB法,见价键理论)。
以下只介绍分子轨道法,它是原子轨道对分子的推广,即在物理模型中,假定分子中的每个电子在所有原子核和电子所产生的平均势场中运动,即每个电子可由一个单电子函数(电子的坐标的函数)来表示它的运动状态,并称这个单电子函数为分子轨道,而整个分子的运动状态则由分子所有的电子的分子轨道组成(乘积的线性组合),这就是分子轨道法名称的由来。
简史 分子轨道法的核心是哈特里-福克-罗特汉方程,简称HFR方程,它是以三个在分子轨道法发展过程中做出卓著贡献的人的姓命名的方程。
1928年D.R.哈特里提出了一个将N个电子体系中的每一个电子都看成是在由其余的N-1个电子所提供的平均势场中运动的假设。
这样对于体系中的每一个电子都得到了一个单电子方程(表示这个电子运动状态的量子力学方程),称为哈特里方程。
使用自洽场迭代方式求解这个方程(见自洽场分子轨道法),就可得到体系的电子结构和性质。
哈特里方程未考虑由于电子自旋而需要遵守的泡利原理。
1930年,B.A.福克和J.C.斯莱特分别提出了考虑泡利原理的自洽场迭代方程,称为哈特里-福克方程。
它将单电子轨函数(即分子轨道)取为自旋轨函数(即电子的空间函数与自旋函数的乘积)。
泡利原理要求,体系的总电子波函数要满足反对称化要求,即对于体系的任何两个粒子的坐标的交换都使总电子波函数改变正负号,而斯莱特行列式波函数正是满足反对称化要求的波函数。
将哈特里-福克方程用于计算多原子分子,会遇到计算上的困难。
C.C.J.罗特汉提出将分子轨道向组成分子的原子轨道(简称AO)展开,这样的分子轨道称为原子轨道的线性组合(简称LCAO)。
使用LCAO-MO,原来积分微分形式的哈特里-福克方程就变为易于求解的代数方程,称为哈特里-福克-罗特汉方程,简称HFR方程。
RHF方程 闭壳层体系是指体系中所有的电子均按自旋相反的方式配对充满某些壳层(壳层指一个分子能级或能量相同的即简并的两个分子能级)。
这种体系的特点,是可用单斯莱特行列式表示多电子波函数(分子的状态),描述这种体系的HFR方程称为限制性的HFR方程,所谓限制性,是要求每一对自旋相反的电子具有相同的空间函数。
限制性的HFR方程简称RHF方程。
UHF方程 开壳层体系是指体系中有未成对的电子(即有的壳层未充满)。
描述开壳层体系的波函数一般应取斯莱特行列式的线性组合,这样,计算方案就将很复杂。
然而对于开壳层体系的对应极大多重度(所谓多重度,指一个分子因总自旋角动量的不同而具有几个能量相重的状态)的状态(即自旋角动量最大的状态)来说,可以保持波函数的单斯莱特行列式形式(近似方法)。
描述这类体系的最常用的方法是假设自旋向上的电子(α自旋)和自旋向下的电子(β自旋)所处的分子轨道不同,即不限制自旋相反的同一对电子填入相同的分子轨道。
这样得到的HFR方程称为非限制性的HFR方程,简称UHF方程。
从头计算法 原则上讲,有了HFR方程(不论是RHF方程或是UHF方程),就可以计算任何多原子体系的电子结构和性质。
真正严格的计算称之为从头计算法。
RHF方程的极限能量与非相对论薛定谔方程的严格解之差称为相关能。
对于某些目的,还需要考虑体系的相关能。
UHF方程考虑了相关能的一小部分,更精密的作法则须取多斯莱特行列式的线性组合形式的波函数,由变分法求得这些斯莱特行列式的组合系数。
这些由一个斯莱特行列式或数个斯莱特行列式按某种方式组合所描述的分子的电子结构称为组态,所以这种取多斯莱特行列式波函数的方法称为组态相互作用法(简称CI)。
甚至是对于单行列式的RHF方程和UHF方程的从头计算法,计算工作量也是惊人的,主要的困难是RHF方程和UHF方程中的电子间排斥积分的计算,这些排斥积分数目多而且计算繁。
这样就有了一些针对简化排斥积分计算的方案,称为近似计算方法。
近似的RHF方程和UHF方程一般得到的结果较差。
取一些经验参量代替近似RHF方程和UHF方程中的一些积分,同时忽略某些积分,则能得到令人满意的结果。
这种方法称为半经验近似计算方法。
在从头计算法里,分子轨道由组成体系的原子的全部原子轨道线性组合而成。
对于原子轨道(也称基组),有不同的选法。
斯莱特型轨道适于描写电子云的分布,但在计算一些积分时包含对无穷级数的积分,十分麻烦,所以在从头计算法里,常取高斯型函数做为基函数,取一个高斯型函数或数个高斯型函数的线性组合模拟一个原子轨道。
已经有了不少进行多原子体系的从头计算法的标准的计算机程序,如J.A.波普尔的研究集体推出的计算机程序系列:
高斯系列(因其将原子轨道的斯莱特型轨道展开为高斯型轨道,故名)。
该系列第一个公开的版本为高斯70(70指1970年),以后差不多每两年更新一次,现已出到高斯82,功能逐次增多,算法也日趋完善。
赝势价轨道从头算法 从直观的化学观念出发,可以想见,在原子形成分子时,仅仅原子的价层电子发生了较大的变形,而内层电子分布则改变较小。
为了节省计算时间,而又不失去计算精度,从20世纪70年代以来,开始出现只考虑原子的价电子,而把内层电子和原子核看成一个凝固的原子实,用一个模型势来代替内层电子与价电子的相互作用的从头计算法,称为赝势价轨道从头算法。
不同的方案对于赝势的取法不尽相同,且大都能得到与全电子从头算法相近的结果,因而大大节省了计算时间,特别是对于含有重原子的体系,恰好是全电子从头算难以处理的。
赝势价轨道从头算法将会在过渡金属络合催化的量子化学研究方面发挥重要作用。
组态相互作用法 比从头算法精确度更高的组态相互作用法有两个显著的优点:
一是不依赖于试探波函数的形式,就能在原则上提供非相对论薛定谔方程的精确解;二是原则上可用于原子或分子体系的任何稳定态。
组态相互作用的概念是相当简单的,只是将线性变分原理直接用于写成多行列式线性组合的试探函数,而这些行列式又是单电子轨道的反对称化乘积,就可以进行CI计算。
从原理上讲,这种做法对于开壳层几乎就象对于基态那么容易;对于激发态几乎就象对于基态那么容易;对于远离平衡的几何构型几乎就象在平衡点附近那么容易。
但实际上,除了迭代收敛上的困难外,CI计算还不象从头计算法那么容易掌握。
CI计算需要考虑原子基组、分子轨道和组态函数(简称CF,即满足自旋匹配,有时还要满足对称性匹配的斯莱特行列式或斯莱特行列式的线性组合)的选择。
由于对分子轨道和组态函数CF有不同的选择,所以有不同的组态相互作用等级。
同时优化分子轨道和组态函数的组合系数(分子轨道向原子轨道展开,体系总波函数向组态函数展开)的方法称为多组态自洽场法(简称MCSCF法)。
组态相互作用法和多组态自洽场法在原子光谱、电子-原子散射和其他原子体系以及范德瓦耳斯力的研究上是非常重要而有用的。
近似计算法 在量子化学半经验计算方法里,一般基组都只取价原子轨道,引入一些近似方法,忽略掉全部三中心积分、四中心积分和部分双中心积分。
由于没有多中心(三中心、四中心)积分的计算困难,所以半经验近似自洽场方法的原子轨道基组一般都取斯莱特型的原子轨道。
目前最常用的近似是零微分重叠(简称ZDO)近似,即假定某些原子基之间的重叠为零:
φμ(i)φv(i)=δφv
式中φμ和φv为原子轨道基;δφv为克罗内克符号,当μ等于v时,δφv等于1;μ不等于v时,δφv为零。
有不同等级的零微分重叠近似。
近似程度最高的是全略微分重叠(简称CNDO)近似,是J.A.波普尔等在1965年提出的。
在CNDO近似里,取所有不同原子轨道之间的重叠都为零。
这样,双电子排斥积分为:
式中r12为电子1和电子2之间的距离。
只有在μ=v、λ=σ时上式才不为零。
波普尔等还通过与从头算的计算结果的拟合来得到一些参量(这些参量的选取要考虑旋转不变性),以代替某些积分值,提出了两个对CNDO进行参量比的方案,分别称为CNDO/1和CNDO/2。
CNDO/2可计算含有从氢到第三周期元素的原子的分子,并且能得到与实验对应较好的几何构型,所以被广泛用于有机分子的计算。
以后又有人优化出过渡金属的CNDO参数,称为CNDO/M法。
还有人通过拟合光谱数据得到CNDO参数,这种参量化的CNDO方法称为CNDO/S法,可以较好地重复光谱实验的结果,被用于分子的光谱性质的研究。
在CNDO法里,忽略了不同原子轨道上电子的交换相关,因而对两个具有平行或反平行自旋的电子间实际存在的不同相互作用未予适当考虑,尤其当两个电子同属一个原子时。
为了弥补这个缺陷,波普尔等进一步提出了不仅保留单中心库仑积分,也保留单中心交换积分的方法,称为间略微分重叠(简称INDO)近似法。
对于分子的平衡几何构型,INDO能得到与CNDO相近的结果,而对于分子的波谱(核磁共振和电子自旋共振谱)的研究,由于考虑了单中心的交换作用,却能得到CNDO不能得到的结果。
波普尔参量化的CNDO和INDO所得到的能量值都较差。
为了使分子的全面性质都能得到好的结果,M.J.S.杜瓦等用一些有代表性的有机分子的实验生成热和平衡几何构型参数做标准,用最小二乘法优化出INDO的参量,并用一些函数来代替某些参量,对大部分有机分子,都能同时得到较好的几何参数和生成热数值。
杜瓦的这种方法称为改进的INDO法,简称MINDO。
杜瓦等前后推出三个MINDO版本,分别为MINDO/1、MINDO/2和MINDO/3,其中最成功的并且至今还在较广泛应用的是MINDO/3。
MINDO/3对有些有机分子的生成热计算甚至能达到化学精度(4~8×103焦/摩尔)。
MINDO/3对一些有机反应的研究也能得到好的结果。
在ZDO近似中,与从头算最接近的为忽略双原子微分重叠近似法,简称NDDO。
它只对属于不同原子的原子轨道φμ和φv才使用零微分重叠近似,而保留了全部双中心排斥积分。
由于剩下的单中心和双中心积分的数目还是相当多的,所以计算依然很繁,而得到的结果却较差,因此很少有人用NDDO进行计算。
1977年杜瓦等提出了一个改进的NDDO近似法,称为MNDO,它采用半经验模型处理双中心排斥积分,即把双原子的电荷分布间的相互作用看成是两堆电荷分布的多极矩相互作用之和,然后也用实验数据优化参量。
这样,计算时间只比MINDO/3增加百分之二十,而计算结果(与实验数据相比)却比MINDO/3好。
在近似计算方法中,还有休克尔分子轨道法和推广的休克尔分子轨道法。
Xα方法 在量子化学计算方法中,除了MO法和VB法外,还有一种半经验方法,叫做Xα方法,是斯莱特在1951年首先提出来的。
他是将电子交换作用能用一统计平均来近似(用总电子电荷密度乘以某个系数来代替交换能),代入哈特里-福克方程,就得到Xα方程。
Xα这个名词是在推导Xα方程时用的,X指交换能用某种方法近似,同时这个表示式中有一个参数α。
原子和晶体的Xα方程易于求解,由于分子没有空间周期性,在相当长一段时间内,Xα方程不能用于分子的计算。
从1965年以来,发展了多重散射波Xα方法(简称SCF-Xα-SW)。
用它可求解自洽场方程,首先是用松饼罐头近似,将分子分为三个区:
原子内区Ⅰ、原子间区Ⅱ和原子外区Ⅲ。
根据Xα方程,分别对三个区域求解,利用边界条件得到久期方程,从而求得分子轨道和能级。
多重散射波Xα方法对原子簇计算结果较好,也比从头计算法节省计算时间,但计算精度比从头计算法差。
对于三个区域的划分,需要计算者的经验和摸索,所以掌握起来比从头算困难,且应用对象有一定限制(即对于某些分子不大适用)。
自洽场分子轨道法-正文
利用自洽迭代过程处理分子轨道的方法。
分子轨道法认为在分子中存在着一系列单电子空间波函数,称为分子轨道,每一个分子轨道都有确定的能量与之对应,而整个分子的波函数,可以近似地用分子轨道的乘积来描述。
按照泡利原理,总波函数必须是反对称的,即交换任何两个电子的坐标后,波函数将只改变一个符号,斯莱特行列式波函数满足这个要求。
例如,对于一个有2N个电子的闭壳层分子,有一系列正交归一化的分子轨道ψ1、ψ2、…、ψN,按能量从低到高排列,考虑到电子可以有两种自旋状态,α或β,则对应一个分子轨道ψi有两个自旋空间轨道ψiα和ψiβ,在以下的行列式中简写作ψi和徰i,对于分子的基态,斯莱特行列式波函数形式为:
(1)
这种表达方式称为“轨道近似”。
总波函数ψ满足分子中电子运动的薛定谔方程:
(2)
(3)
式中H为分子中电子的总哈密顿算符;h(i)为第i个电子的哈密顿算符,它包含单电子的动能和这个电子与所有原子核间的相互吸引能。
把式
(1)代入式
(2),可以得到分子基态的电子总能量E的表达式:
(4)
(5)
式中r12为电子1和电子2间的距离;Jij和Kij分别为库仑积分和交换积分。
在分子处于基态时,E应为最小值,如何选择一组正交归一化的分子轨道ψ1,ψ2,…,ψN使E成为最小,可以用变分法来达到。
变分结果要求分子轨道满足的方程为:
(6)
这个方程称为哈特里-福克方程,是B.A.福克在1930年导出的,F
(1)称为哈特里-福克算符:
(7)
式中库仑算符Jj
(1)和交换算符Kj
(1)的定义为:
哈特里-福克方程虽然表面上具有通常本征值问题的形式,但却不能用通常的方法求解,这是因为算符F
(1)本身还包含着ψ1、ψ2、…、ψN,需要用迭代法求解。
先给定一组初始的ψi
(1),代入式(7)得到F
(1),然后求解式(6)可得新的一组ψi
(2);再用ψi
(2)重复上面的过程,如此循环,直至最后两次的结果符合到所规定的要求为止。
这个过程称为自洽迭代的过程,这种分子轨道的处理方法称为自洽场分子轨道法。
如果把分子轨道写成原子轨道线性组合的形式,则组合系数满足罗特汉方程。
这种自洽场分子轨道法在量子化学研究中广为应用,是很重要的理论方法。
分子轨道理论-正文
一种化学键理论,是原子轨道理论对分子的自然推广。
其基本观点是:
物理上存在单个电子的自身行为,只受分子中的原子核和其他电子平均场的作用,以及泡利不相容原理的制约;数学上则企图将难解的多电子运动方程简化为单电子方程处理。
因此,分子轨道理论是一种以单电子近似为基础的化学键理论。
描写单电子行为的波函数称轨道(或轨函),所对应的单电子能量称能级。
对于任何分子,如果求得了它的系列分子轨道和能级,就可以像讨论原子结构那样讨论分子结构,并联系到分子性质的系统解释。
有时,即便根据用粗糙的计算方案所得到的部分近似分子轨道和能级,也能分析出很有用处的定性结果。
分子轨道理论
分子轨道理论
氢分子离子的分子轨道 正如在原子轨道理论中,氢原子的严格解提供了进一步发展的理论模式,氢分子离子H娚中,单个电子在固定核间距R的双质子场中的波动方程解,是分子轨道理论进程中的基石。
H娚的分子轨道用符号σ、π、
、…表征,对应于精确解中的量子数m=0,±1,±2,…,它描述相对于核间距R的轨道对称行为。
此外,还需用g和u表征相对于分子中心反演的对称行为。
综合起来,H娚的分子轨道用σg、σu、πg、πu、…等符号表征,借助精确求解固定核间距R的波动方程获得。
图1给出两个最低轨道1σg和1σu的能量E随R的变化曲线。
1σg能级有一极小值-1.20Ry(里德伯能量),出现在R=2a0处(a0为玻尔半径,图2),代表基态;当R增大以至无穷时,1σg能量趋近-1.0Ry。
两者差值0.20Ry就是H娚的离解能。
1σu的行为不同,能量随R减小而单调上升,显示排斥态的本质。
1σg和1σu也被称作成键轨道和反键轨道。
1σg均取正值,1σu则在中心两端发生符号变化,但极值均出现在质子所在处,且伴随R变小;1σg在核间区数值增大,描写了电子在分子中的转移。
随着R的增大,1σg和1σu的函数值渐近于式
(1)~
(2):
(1)
(2)
式中Ng和Nu是归一化系数;
和
分别是两个组成氢原子上的1s原子轨道。
图3
给出R=2a0时,1σg和1σu轨道的精确值和按式
(1)与
(2)的近似值的比较,说明式
(1)与
(2)的近似程度是很好的。
式
(1)和
(2)表示,分子轨道可以近似地当作原子轨道的线性组合,简写为LCAO(见量子化学计算方法)。
当R很大时,结果是准确的,即使R达到分子核间距大小时,也给出不错的结果。
这一点启发了对复杂分子也可采用LCAO方法去寻找近似分子轨道。
这是因为在分子中,靠近一个核的电子主要受到该核的势场的作用;而受到其余核的联合作用,则小得多,因此在近核处,分子轨道必定近似于该核的原子轨道。
对于整个空间的任何一点,可以设想分子轨道由有关的原子轨道线性组合而成。
一般的形式是:
(3)
式中ψ代表分子轨道或轨函;φv是属于各组成原子的原子轨道;cv是待定系数,由变分法确定。
还应指出,LCAO是一种可行的近似方式,但不是唯一的近似方式。
任意双原子分子的分子轨道 用原子轨道线性组合法LCAO近似来讨论任意双原子分子中,分属两个原子的一对原子轨道形成分子轨道的最优条件。
这时,式(3)采取以下简单形式:
ψ=c1φa+c2φb(4)
代入波动方程Hψ=Eψ,得到近似能级E:
(5)
式中Haa和Hbb分别是原子轨道φa和φb的库仑积分,可看作φa和φb的能量,即Haa=Ea,Hbb=Eb;Hab=β称共振积分,与φa和φb的重叠情况有关,一般取负值,Sab称重叠积分,可当作零处理,而不影响定性结论。
最优条件由E取极值确定,即:
(6)
上式给出一组c1、c2的齐次方程组,由系数组成的久期方程(见休克尔分子轨道法)得到:
(7)
设Ea<Eb,则式(7)的含意可用图4表示。
分子轨道理论
由此看出,有效的成键作用决定于B值的大小,后者又与原子轨道能量差|Ea-Eb|以及重叠情况β有关,从而可归纳为三个条件:
①能量近似条件:
指|Ea-Eb|越小越好,当|Ea-Eb|=0时,B最大,等于|β|。
②最大重叠条件:
φa与φb的重叠越大(图5
),β的绝对值也可能越大。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 量子 化学 计算方法
![提示](https://static.bdocx.com/images/bang_tan.gif)