非绝热-连化学物理研究所课件.ppt
非绝热反应动力学,韩克利 中国科学院大连化学物理研究所,1.多原子分子的哈密顿量和波函数 H=Te+TN+V e e+V e N+VN N H(Ra,ri)=E(Ra,ri)2.B-O近似下核运动和电子运动的分离 a.核不动时电子运动的薛方程 电子哈密顿量：He=Te+V e e+V e N(Te+V e e+V e N)e(Ra,ri)=(U-VN N)e(Ra,ri)E e(Ra)=U-VN N(Ra)b.核运动薛方程 分子总波函数(Ra,r i)=e(Ra,r i)N(Ra)TNe(Ra,r i)N(Ra)+(TE+V E E+V E N+VN N)e(Ra,r i)N(Ra)=Ee(Ra,r i)N(Ra),研究分子的电子结构的理论称为量子化学，而处理原子核之间的相对运动的理论称为动力学。通过量子化学计算可以得到原子或分子间的相互作用势，即势能面。在此基础上即可做动力学研究。不同电子态的势能面间相互作用很小（没有交叉）时，其反应过程是绝热的。但当势能面间的相互作用很大或发生交叉，即可能会发生非绝热动力学过程时，进行动力学研究时就必须考虑电子的运动。现在的量子化学理论处理势能面间的相互作用基本是采用微扰理论。,The effect of the spin-orbital coupling on the reactions F(2P3/2,2P1/2),Cl(2P3/2,2P1/2)+H2,HD,意义,实际应用:作为Cl2+H2链式反应的重要基元反应，Cl+H2及其同位素反应在大气化学和光化学空气污染中占有非常重要的地位。理论：最典型的三原子反应体系，检验双分子反应速率理论的样板。动力学理论的研究：非绝热动力学能使人们更加深入细致地了解化学反应的本质。,刘国平从实验上发现自旋轨道激发态反应活性更大。对三原子分子体系绝热理论研究相当成熟，对非绝热研究还比较少。,Cl(2P3/2)、Cl(2P1/2)+H2,参考文献,参考文献,哈密顿,展开基函数,K=k+,基函数可以定义为：,势能矩阵元,三个非绝热的电子态的对角元,绝热表象和非绝热表象的变换,绝热态和非绝热态可以通过一个正交变换相联系，这个变换可以用变换角来表示：变换角的物理意义，Px与Pz的转动，可以通过轨道角动量来求。,非绝热势能,非绝热基组下的能量为,自旋轨道耦合,使用的MOLPRO量化计算程序。,对基态和激发态使用多体展开基态双体项其中三体项,激发态双体项,激发态的HCl是排斥态，选用排斥函数为保证核间距很小时,仍是排斥势最大误差低于1.0Kcal/Mol,V2的函数形式角度函数形式的选择考虑到了对称性的问题。,自旋轨道相互作用 A,其中damp是,自旋轨道相互作用 B,以上的函数形式除了基态以及激发态的三体项外是参考别的势能面自己创造的。,Cl+H2共线情况的基态势能面,Cl+H2共线情况下的激发态的势能面,Cl+H2的计算结果几率 J=0.5,Cl+H2的计算结果几率 J=15.5,基态几率差别不大，激发态有很大的差别。在我们的势能面上计算的结果比CW势能面计算的结果大的多。在低能处，由于具有更高的内能，激发态几率大，随能量的增加，基态几率迅速增加，并超过激发态。随J增大，几率交叉点高能移动。,Cl+H2的计算结果 截面,Cl+H2替换比较,Cl+H2替换比较,替换V1以及V2 基本没有影响替换A 基本没有影响替换B 减小，而且与CW结果类似，所以B的影响最大。,Cl+H2 势能矩阵元B 比较,Cl+H2 势能矩阵元B 比较,Cl+H2 势能矩阵元B 比较,Cl+D2的计算结果,自旋轨道相互作用B对Cl（2P3/2)+H2反应的影响比较小。自旋轨道相互作用B对Cl（2P1/2)+H2反应有很大的影响。我们势能面中的自旋轨道相互作用更合理些。,回顾和展望,发展了三维量子散射理论处理非绝热光解动力学、非绝热反应动力学、传能动力学以及光电离（三个势能面以上）的理论计算方法。进一步工作包括发展处理更多原子的非绝热过程的理论计算方法。,致谢,主要合作者：解廷献、张岩。感谢张增辉教授和刘国平教授。感谢国家自然科学基金和973的支持。Comments from Kopin Liu：From my first reading last night,it is indeed an outstanding work,and will be like dropping a bomb in the community when it is published.Congratulation!,液相飞秒化学,染料分子在溶液中的超快（飞秒）振动弛豫动力学。叶绿素a分子电子态之间的能量转化。,Femtochemistry,Gas PhaseMolecular Beam,Condensedphase,Mesoscopicphase,ClustersDropletsNanocavitiesNanoparticles metallic colloidalNanorods,Dense FluidsSuper FluidsLiquidsSolidsSurfacesInterfacesPolymers,Femtobiology,StructuresUED X-rays,Control,The Schematic Diagram of FS TR SEP FD,The Schematic Diagram of the Experimental Setup,叶绿素a电子态间的能量转移,飞秒激光出现后，立即被用于研究光合作用中的超快过程。例如：LH I 和 LH II以及反应中心。叶绿素a和b构成了这些分子。对叶绿素a的研究亦非常多，但很少有人研究电子态之间的能量转移。,致谢,国家自然科学基金委东北师范大学,电子相关能,自洽场(SCF)方法中，考虑了电子间平均作用,但没考虑电子间瞬时相关,没有考虑“Coulomb孔”,降低了其它电子出现的几率。过高地估计了电子的排斥能，使得计算得到的体系总能量比实际值要高。电子相关能是指该Hamilton量在该状态的精确本征值和它的限制的HF极限期望值之差。,电子相关能,电子相关能在体系总能量中占的比例为0.31%。但在研究电子激发、反应途径(势能面)、分子解离等过程时,由于相关能的数值与一般化学过程中反应热或活化能具有相同的数量级,所以必须在HF基础上考虑电子相关能。处理相关能组态相互作用理论(CI),耦合簇理论(CC)和微扰理论(MP)等。,电子相关能,从一组完备的单电子基函数k(x)出发,可构造出一个完备的行列式函数集合k k=(N!)-1/2|k1(x1)k2(x2)kN(xN)|在组态相互作用(CI)方法中,将多电子波函数近似展开为有限个行列式波函数的线性组合(CI展开),电子相关能,multireference configuration interaction(MRCI)波函数求的能量,是通过complete active space self-consistent(CASSCF)波函数计算的能量。是scaling factor,基函数,Slater型基函数(STO)和Gauss型基函数(GTO)。STO在描述电子云分布方面较其它基函数具有一定的优越性。GTO计算包含Slater原子轨道的三中心和四中心积分时比较容易。可以用Slater函数向Gauss函数集合展开的方法来选择基函数。分裂价基:每个价层轨道用若干个新的具有不同指数的STO-KG来逼近,Cl+H2能量示意图,势能面的意义,势能面是反应速率和反应历程分析与计算的基础，势能面的形状反应出整个化学反应过程的全貌以及反应的始终态、中间体和过渡态的基本态势。可以说势能面提供了反应过程的舞台，包含了反应过程的信息库。势能面是动力学的基础和前提条件，势能面的好坏对计算结果影响很大。,势能面的来源,基于光谱学的、气体动力学和分子束散射等实验测定，是基本的方法，但是只能提供势能面的局域特征。量子化学理论计算的方法得到，理论计算原则上可以给出整个势能面的信息。,拟合势能面的原因,拟合势能面主要用于做动力学计算，进行一次动力学计算，需要成千上万个点的能量，对每一个点的能量都进行计算，在理论上是可行的，但是现在的硬件水平决定了这是很困难的。拟合势能面就是通过一些已知的分立点的能量来计算任意一点的能量。,动力学计算对势能面的要求,对拟合的势能面的要求：给出决定分子体系几何构型的参数通过简单计算就能得到体系的能量。因此拟合的势能面有许多种形式，可以给出解析式，例如多体展开，也可以是非解析式，例如泰勒展开,多体展开基本步骤,一、首先计算一些分子构型下的能量二、选取合适的含有待定系数的基函数三、确定待定系数四、编写拟合势能面的程序,计算某些分子构型下的能量,选取待计算的能量点。一般情况下，首先寻找该体系的势垒和势阱的位置，在这些位置附近要多计算一些点。这个过程主要由量化计算程序来完成，例如用的非常广泛的高斯程序。,选取合适的基函数,函数形式的选取主要是凭经验。要求函数尽可能的接近真实的势能面形状，考虑势能的渐近形式，例如bw2势能面就含有键长的倒数项，这样就能保证当核间距很小时，分子势能变的很大。现在有很多人使用不同的函数来拟合势能面。,确定待定参数,由于选取的基函数中一般含有非线性的待定参数，因此这是一个非线性的拟合过程。最常用的、最广泛的拟合方法就是最小二乘法。,最小二乘法,由于Q取极值，因此对m个待定系数的导数等于零：这样就得到m个方程，求解这m个方程组成的方程组可以得到待定系数的值。这些方程一般是非线性的，有许多方法来求解，例如牛顿迭代法、梯度法等。,编写计算势能面的程序,根据计算出来的系数，编写计算势能面的程序，基本上就是各个量的加和，这一步最容易的。,
