蒋彬[1]2012年在《从气相到金属表面的小分子势能面与量子动力学的理论研究》文中研究指明分子反应动力学是化学研究的前沿领域之一,也是联系微观动力学和宏观动力学的桥梁。原子和分子的微观运动遵循着量子力学的规律,原则上,我们可首先分离反应体系中电子和核的运动,然后运用量子化学方法计算不同构型下的势能,构建原子核运动的势能面,最后求解核运动的Schrodinger方程得到量子反应动力学的全部信息。势能面是精确研究反应动力学的基石,也是理论与计算化学的核心课题之一。势能面的构建在研究分子的性质和反应等方面起着关键作用,是在原子水平和量子态层次深入认识和揭示化学反应规律的基础。对反应体系建立精确的势能面以后,就可以进行量子动力学计算。量子态分辨的动力学在最基本的原子与分子的层次上对化学反应进行探索,从而为化学反应机理的研究提供了最详细的信息,使人们可以深入理解化学反应的本质。本文从气相中典型的基元反应出发,分别研究了ClH2,BrH2勺非绝热势能面,H20的非绝热势能面与态态光解动力学,并延伸到重要的金属表面催化基元反应,研究了H20在Cu(111)面的解离吸附动力学过程,获得了相关反应过程的细致的动力学机制。气相中Cl+H(≒)H+HCl和Br+H2≒H+HBr反应都是重要的基元反应。Cl(2P3/2,2P1/2)+H2反应中非绝热效应引起的激发态反应性一直存在着争论。多年前的实验结果发现Cl(2P1/2)+H2的反应截面比Cl(2P3/2)+H2反应截面更大,与后来的理论研究结果相悖。另一方面,基于较早构建的势能面,量子力学计算的H+HBr抽取反应速率常数比实验值大了两倍以上；另外Br(2P1/2)+H2中非绝热效应的大小也存在着争论。对这两个体系,我们利用高精度从头算方法和很大的基组,精确地计算了能量较低的三个电子态的势能面和旋轨耦合作用,并转换到非绝热表象中构建了包含旋轨耦合项的6×6非绝热势能矩阵。对于Cl(2P3/2,2P1/2)+H2反应,进一步的量子动力学计算验证了最新的实验结果,表明了Non-Born-Oppenheimer效应只在很低碰撞能下主导反应,导致Cl(2P1/2)有更大的反应性,而较高碰撞能下,Born-Oppenheimer近似就可以有效地描述反应动力学。对于H+HBr抽取反应,我们的绝热势能面在较大的温度区域很好的重现了实验测得的速率常数,比以前的势能面有很大的改进。同时我们的高精度非绝热势能面相对于之前的半经验拟合势能面,可以对Br(2P3/2,2P1/2)+H2可能发生的激发态淬灭和非绝热反应进行更好的描述。研究表明,Br(2P1/2)+H2(v)可以淬灭为Br(2P3/2)+H2(v+1),而非绝热跃迁效应则很小,即发生非绝热反应生成HBr+H的几率几乎可以忽略。进一步的量子动力学研究已经展开,可以激励更多对这个体系的实验研究。H20分子通过激发态B态的光解离也是一个典型的非绝热过程。最近,杨学明小组在B态吸收光谱中5个最强的吸收峰位置进行了高精度的实验动力学研究,得到了最详细的微分截面,但是实验中发现的现象并不能得到清楚的解释。我们在高水平从头算基础上构建了新的B-X非绝热势能面,发展了基于Chebyshev实数波包传播的计算态-态光解动力学的方法和程序,应用到B态非绝热解离过程中,第一次实现了三原子光解的态态微分截面的量子力学计算,并与实验值吻合得较好。我们的计算结果解释了实验上发现的OH(X)和OH(A)产物态态微分截面的明显区别,并清晰的表明两个产物通道经历不同的解离路径对解离图像有着重要影响。我们的计算结果对今后H2O(B)光解动力学的实验和理论研究有重要的启示意义。H20分子在过渡金属表面的解离吸附过程是非常重要的气相-表面反应范例,也是工业上众多多相催化过程必不可少的步骤。为了最终控制和提高催化过程的效率,对这个反应动力学更好的理解是至关重要的。但是,到目前为止,几乎没有实验或理论工作者对这个动力学过程进行研究。作为第一次的尝试,我们考察了H20的振动模式对反应性的影响。我们从一个包括了H20分子所有内振动态的减维模型出发,首先构造了六维的全域势能面,发展了基于Chebyshev波包传播的过渡态波包计算方法和程序,得到了H20分子不同初始振动态的解离几率。研究表明,H20分子与CH4分子在金属表面的解离有很大不同,CH4只有对称伸缩振动激发比相同能量平动能的增加更能增大反应性能,而H2O的对称伸缩、弯曲、反对称伸缩三个振动模式都能显著的增大反应。我们的理论预测将激励更多的态解析的动力学实验对这个重要的体系进行研究。
权伟龙[2]2006年在《含时量子波包法研究反应性散射F+HBr和Br+H_2》文中研究指明量子含时波包法,是研究原子分子碰撞与反应动力学最重要的理论方法之一,已经在单原子+双原子分子、双原子分子+双原子分子以及离子+双原子分子碰撞反应中得到了广泛应用,并取得了许多重要的研究结果。本文应用含时量子波包法研究了F+HBr和Br+H2反应性散射动力学。利用Persky等人构建的LEPS型势能面首先研究了F+HBr反应。研究结果表明,在碰撞能量很低时(几乎为零),反应过程就能够发生,且反应几率随能量急剧增大。随着碰撞能量的升高,反应几率略有下降,符合无阈能反应的一般特征。该反应的初始态转动-振动激发效应的结果表明, HBr分子的振动激发,会引起反应几率的降低,在低能碰撞时,下降尤为明显,从而不利于F+HBr反应的进行;HBr分子的转动激发也会引起反应几率的下降,但主要降低高能部分的反应几率。计算得到的反应截面在碰撞能量很低时下降很快,然后变得比较平坦,其变化趋势与实验测量的结果相吻合,呈现出强放热反应的基本特征。对速率常数的计算发现,在所考虑的整个温度区间内,含时量子波包法所得结果要比经典轨线法计算的结果略大一些,但更接近实验值,表明隧道效应在F+HBr反应中比较重要。另外,还利用含时量子波包法在e-LEPS和ab initio PES两个势能面上研究反应Br+H2。计算表明,在ab initio势能面上,Br+H2的反应阈能要比e-LEPS势能面上的反应阈能低一些,而且ab initio势能面上的反应几率、反应截面和速率常数都要大许多。并结合势能面的特征分析了这种差别的原因,探讨了势能面对反应特性的影响。我们还研究了H2分子的振动转动激发对反应特性的影响,结果显示,无论是在e-LEPS还是在从头算势能面上,H2的振动激发都会通过降低反应阈能强烈的促进反应进行,转动激发虽然对反应也有促进作用,但相对于振动激发而言,其作用要弱得多。把两个势能面上的速率常数和实验数据以及前人的计算结果进行比较,发现ab initio势能面上的速率常数更接近实验值,从而ab initio势能面要比e-LEPS势能面更适合于进行准确地量子动力学计算。
梁景娟[3]2012年在《多原子体系反应势能面及动力学》文中进行了进一步梳理本文主要包含两方面的内容，一部分是应用ab initio理论构建了He+H_2~+和He_2+H+离子分子反应的从头算势能面并计算了KH分子Χ1Σ+与Β1Π的解析势能函数。另一部分是利用准经典轨线方法研究了He+H_2~+、He_2+H+及H+HBr反应的标量及矢量性质。文章主体内容共分为六章，前两章介绍构建势能面及进行分子反应动力学的基本理论，后四章详细介绍了本人对He+H_2~+、He_2+H+、H2Br及KH不同体系所开展的研究工作。在第三章，我们运用多参考组态相互作用（MRCI）方法及aug-cc-pV5Z基组构造了HeH_2~+体系基态（12A'）从头算势能面。拟合时，我们共选择8840个能量点，并且采用多体展式及Aguado-Paniagua（AP）函数来描述原子间的相互作用，拟合误差为0.0677Kcal/mol。大规模从头算前，为选择合适的活性轨道，我们计算了(RHeH~+=2 bohr,RHH+=2 bohr,θ=180o）在不同活性空间下的能量。分析发现，当活性轨道数选为某值时，我们再增加轨道数，能量变化可忽略不计，因此我们最终选取该活性轨道进行从头算。为检测新势能面，我们绘出了势能面的等高线图，曲线非常光滑；我们采用准经典轨线方法计算了He+H_2~+→HeH~++H反应在碰撞能范围0-2 eV，反应物分子v=0-3四个不同振动态下的积分反应截面，计算结果比其它理论结果更接近实验结果。这些充分说明了我们构建的势能面精确可信，可作为动力学计算的基础。提取反应HeH~++H→He+H_2~+与交换反应H′+HeH~+→H++ HeH′都是典型的分子离子反应，两个反应享有同样的反应物，但是却有不同的质量组合。提取反应的质量因子为0.3994，属于light-light-heavy (LLH)模型；交换反应的质量因子为0.0405，属于LHL模型。在以前的研究中，人们已发现平动能对提取反应HeH~++H→He+H_2~+的标量和矢量性质都有显著的影响且平动能有利于该反应的进行。但是，反应物的振动激发效应还没有人研究过。因此，我们基于Ramachandran等人构建的HeH_2~+体系基态（11A/）势能面采用准经典轨线方法，研究了反应物HeH~+振动激发对两个反应通道标量及矢量性质的影响。首先，我们计算了当反应物的振动量子数分别为v=0，1，2时，两个反应通道的反应物的振动能、转动能、平动能所占的比例，揭示了反应物的振动能对产物的能量分配的影响。接下来，我们详尽研究了反应物处于不同振动态时，产物的振动、转动分布。为了比较产物振动激发对两个通道的作用，我们将反应物处于不同振动态时，产物相同振动态的转动分布进行了对比。最后，我们研究了振动激发对两个反应立体动力学性质的影响，我们计算了反应物不同振动态时的PDDCSS (2πσ)(dσ00 /dωt)，(2πσ)(dσ22 +/dωt)，(2πσ)(dσ21 /dωt)，P(θr)和P (φr )分布，发现提取反应HeH~++H→He+H_2~+的散射、取向程度比交换反应H′+HeH~+→H++ HeH′明显，但反应物的振动激发更有利于交换反应的进行。我们又研究了HeH~+H+/HeD+H+/HeT+H+提取反应的同位素效应对立体动力学性质的影响。当反应物的质量增加时，反应物的取向、定向及前向散射效应变弱，反应表现出明显的同位素效应。在第四章中，我们构建了HeH~+(X1∑+) + He交换反应势能面。在从头算过程中，我们采用d-aug-cc-pV5Z基组及MRCI方法共计算了16789个11A'态点。拟合时，我们采用非线性最小二乘法，拟合15682个能量低于10000cm-1的点。为了验证势能面的精确度，我们采用多种方法进行检测。我们基于两体项的势能曲线计算了光谱常数并与实验结果进行了对比，理论计算与实验的结果符合很好；计算了11A'解离极限H++He+He与H+He++He的能级间隔，结果与现有理论值一致；我们通过势能面计算出了稳定点的几何坐标和能量，并与从头算优化的结果及他人理论计算的结果进行了对比，三者符合得很好；绘出了势能面的等势线图及最小反应路径图，并与Panda等人绘制的结果进行了比较，结果表明两个势能面的基本特征相似，但是我们提供的反应路径在解离区域更加平滑；我们进一步计算了He' + H+He (v = 0, 1, 2, j = 0)→He'H++ He的积分反应截面，结果与量子结果符合很好。以上这些都充分说明了我们构建的势能面的精确、可信。He'+ H+He→He'H++ He是典型的离子分子反应，属于heavy-light-heavy（HLH）质量组合。但是关于该反应的动力学研究较少。在该反应的最小反应路径上有一个约为0.577eV的势阱，易生成长寿命的中间络合物，可是Bhattacharya等人[111]并未发现长寿命的中间络合物，取而代之的是发现了由直接反应机制产生的短寿命的中间体。Bhattacharya等人并未对该反应机制进行进一步的探究，因此，在第四章中我们采用准经典轨线方法对HeH~++ He'→He + He'H+交换反应的动力学性质进行了全面地研究。首先，我们研究了碰撞能对反应的影响：我们计算了Ec= 0.05 eV、0.1eV、0.2 eV、0.3 eV、0.5 eV时，反应时间；反应物平动能、振动能及转动能占总能量的比例；原子间距随时间变化图；反应物的振动、转动分布。研究发现，反应机制与碰撞能有关：当碰撞能较小时，反应以间接反应为主，反应时间较长，振动能所占比例较大；当碰撞能增加时，直接反应机制比例增加，间接反应机制比例减小，反应时间减小，平动能所占比例增加。低能时，反应轨线易被势阱所俘，形成长寿命的中间物，反应以间接反应为主；高能时，轨线较难俘获，反应一般为直接反应机制所主导，不存在长寿命的中间体。在计算了产物的振动及转动分布时，我们发现在所研究的碰撞能范围内，产物主要分布在振动基态，很少处于振动激发态。这种现象是由于振动基态与第一激发态的能级间隔较大，分子很难激发造成的。产物最可几转动态为j=1。振动、转动相关的积分反应截面都随碰撞能的增加而减小。我们还研究了碰撞能对产物分子散射方向及转动角动量的取向及定向效应的影响。我们发现，产物分子倾向于前向散射，转动角动量j′在与k平行的方向上有明显的取向，似乎分子倾向于在与散射平面平行的平面内转动，且当碰撞能增加时，该趋势减弱。接下来，我们采用准经典轨线方法研究了振动能、转动能及同位素取代对交换反应立体动力学性质的影响。计算时选取碰撞能为0.2eV。研究表明He'H+产物表现为明显的前向散射；产物分子转动角动量j′沿y轴方向有明显的取向且定向于y轴的正方向。当反应物分子振动能或是转动能增加，或因同位素取代反应物分子质量增加时，前向散射效应增强；产物分子转动角动量j′在与k平行的方向上取向减弱，在与k垂直方向上取向增强；沿y轴方向取向及定向效应增强。总体说来，振动能对反应的影响较为显著，转动能对反应的影响不明显。在第五章，我们研究了碰撞能对H' + HBr (ν= 0, j= 0)→H + H'Br交换反应的影响并计算了H + HBr (ν= 0, j = 0)→H_2~+ Br反应的热速率常数及积分反应截面。XH_2(X=卤族元素)的H提取和交换反应一直以来是理解和发展基础化学的动力学和反应动力学的基础和重要手段。由于Br原子比F、Cl原子电子数量多，原子质量大，因此想要构建BrH2体系精确的从头算势能面很难。Kurosaki等人采用MRCI方法及aug-cc-pVTZ基组构建了MB3势能面，此后很多理论计算都基于该势能面展开。但是并没有关于H+HBr交换反应的立体动力学性质的研究。因此，为了更全面的了解H+HBr交换反应的特点，我们基于MB3势能面进行了准经典轨线计算。首先，我们计算了碰撞能对该反应标量性质的影响。1.我们计算了碰撞能在0.5 - 2.0eV范围内的积分反应截面，结果与基于该势能面的量子结果符合的很好。因为该势能面有一个高度约为0.42 eV的势垒，因此碰撞能为0.5eV时的积分反应截面很接近于零，一旦碰撞能克服愈能的阻碍作用，积分反应截面随碰撞能增加迅速增加，表明交换反应存在大的acceptance cone。2.我们计算了碰撞能在0.5 -2.0 eV范围内对H + HBr (ν0= 0, j0= 0)交换反应产物能量分支比的影响，即振动能、转动能、平动能在总能量中所占的比例。当碰撞能增加时，平动能的比例显著减小，转动能的比例明显增大，振动能比例变化不大大。尽管平动能比例一直减小，碰撞能始终占据能量中的最大比例，显示了势能面的排斥性。3.我们计算了碰撞能为0.5、0.7、1.0、1.5、2.0eV时，反应的振动分布即产物H'Br不同振动态ν'的积分反应截面。产物的最可几振动态量子数始终为ν'= 0，并且该态相关的积分反应截面随碰撞能增大而减小。为了探究反应机制的特点，我们绘出了碰撞能为0.5, 0.7, 1.0及2.0 eV时，所有反应轨线的H'Br, HBr, and H'H原子间距随反应时间的变化。我们只发现了直接反应机制，这是由该反应路径上有一个高为12.8Kcal/mol的势垒造成的。在弄清楚反应机制特点之后，我们研究了碰撞能对立体动力学性质的影响。我们计算了极化微分反应截面(2π/σ)(dσ_(00) /dwt)，( 2πσ)(dσ_(22) +/dωt)，( 2π/σ)(dσ21 /dwt)，及P (θ_r)和P (φr)分布。从(2π/σ)(dσ_(00) /dwt)看出，产物分子主要是后向散射，当碰撞能增加时，向后散射的趋势变弱。正如书中[101]所描述的，在排斥势能面上，产物分子倾向于后向散射，如果是LHL模型，既使是在吸引型势能面上，产物分子也易表现出后向散射的趋势。所以作为典型的LHL质量组合，H' + BrH交换反应的产物H'Br主要表现为后向散射，这与势能面的排斥性质也有关系，且当碰撞能增加时，向后散射的趋势变弱。(2πσ)(dσ_(22) +/dωt)显示产物在沿y轴的方向上有明显取向效应，当碰撞能增大时，取向变弱。(2π/σ)(dσ_(21) /dwt)表明产物在沿x+z方向上有取向；在θ_ot=160时，反应物转动角动量存在强烈的极化，当碰撞能增大时，极化变弱，分布趋于各向同性。从P (θr)分布看出，取向分布关于k呈柱形对称分布，且j′在与k垂直方向取向明显。当碰撞能增加时，j′在与k垂直方向取向程度先减小后增加。取向因子<P_2(j',k)>随碰撞能变化趋势一致，与韩克利等人预测的LHL质量组合相一致。从P (φr)分布可以看出在产物H'Br倾向于在平行与散射平面k-k'的平面中转动，且转动角动量沿y轴负方向定向。当碰撞能增加时，定向效应变强且产物分子倾向于由“平面外机制”变为“平面内机制”。在本章最后，我们尝试采用准经典轨线计算了H + HBr (ν0= 0, j0= 0)→H_2~+ Br提取反应的速率常数，并与多组实验及理论结果进行了对比。本次计算结果与基于MB3势能面的含时量子波包计算的结果符合很好。尽管KH分子Χ~1Σ~+态在实验和理论方面都得到了广泛研究，但是到目前为止，解析的势能函数尚未发现。与Χ~1Σ~+态被广泛关注相比，Β~1Π态的有关研究较少。2009年，Lee等人第一次在实验测量了振转能级及光谱常数，这些数据可以作为理论计算的参考。目前已有的理论值与实验值相差较大，这些都说明尚有对该态进行精确的理论计算的必要。因此，我们采用多组态相互作用及大基组计算了Χ~1Σ~+及Β~1Π态势能曲线。从头算的过程中，我们计算时分别考虑了K原子三组不同电子相关（1、考虑3s~23p~64s~1，2、考虑3p~64s~1，3、只包含价电子4s1）和两组不同的活性空间（1、包含3d轨道，2、不包含3d轨道）。为方便进一步研究，势能曲线被拟合成解析的势能函数，函数形式为Murrel–Sorbie（MS）。为了检测势能函数的质量，我们基于势能函数计算了Χ~1Σ~+及Β~1Π振转能级及光谱常数，并与实验值进行了对比，结果与实验值符合很好。为了更加清楚地验证本次计算的势能函数，我们将转动能级T(ν,J),拟合为T (0 ,J)=A+A(J(J+1) 1)+A(J(J+1) 1)2000102形式，从而得到了三个光谱常数。两次计算的ωe相差较大，通过分析可知，实验处理数据方式较为粗糙，结果精确度稍低。为了进一步验证势能函数，我们还计算了Β~1Π态ν=0与Χ~1Σ~+态ν=0, 1之间P、R、Q支跃迁频率，并与实验值进行了比较，误差范围仅为0.1-0.95cm-1。
王华阳[4]2006年在《重—轻—重反应体系的动力学研究》文中认为近些年来，散射共振态的研究一直是化学动力学研究的一个重要的前沿课题。散射共振态在大多数化学反应中起决定性的作用，它控制了化学反应的分支比，产物的能量分布和角分布等重要的特征性质。对于探索各类化学反应的详细机理，提高反应体系中不同能量形式(平动能，转动能，振动能和电子激发能等等)的利用率，以及如何控制化学反应的方向及速率等各类研究，均需要对化学中的散射共振态进行详细的考察和分析。化学反应的散射共振态是一种综合的量子效应，主要包含量子化学理论和反应散射理论两方面的综合研究。其中，化学反应的各类势能面是研究散射共振态的最基本也是最重要的工具，通过对势能面的计算和分析，可以获得散射共振态的各类特征性质，如共振能，共振宽度或共振寿命等等。最近对于散射共振态的实验研究的飞速发展，更是大大地推动了理论化学家对于散射共振态的研究。其中，D M Neumark等人的光分离光谱实验研究和K P Liu等人的交叉分子束实验更是取得了重大的成就，对此学科的发展意义非凡。而对于散射共振态的理论计算研究还有很多问题尚未解决。 本论文对若干重-轻-重反应体系的散射共振态做了动力学的理论研究，主要内容包括以下几个方面： 第一章的绪论中，简要介绍了散射共振态研究的历史背景，共振态的不同类型，以及目前对散射共振态的实验研究和理论研究的发展现状。而后，关注了该领域发展的最新动态，以及对散射共振态研究的发展前景。 第二章中，详细介绍了目前研究散射共振态的理论方法。理论研究主要包括势能面的构建和量子反应散射理论对散射共振态的计算。然后介绍了偏分势能面的理论基础和构建方法，和如何利用偏分势能面研究化学反应的散射共振态，主要包括一维方势阱模型和构建散射矩阵两个方面。 第三章，一维方势阱模型被应用于对目标体系的共振态寿命的计算研究。首先介绍了该模型的理论依据，和目标体系的理论研究状况，然后分别对两个模型化学反应，I+HI(v)→IH(v')+I，Br+HBr(v)→BrH(v')+Br做了相关
参考文献：
[1]. 从气相到金属表面的小分子势能面与量子动力学的理论研究[D]. 蒋彬. 南京大学. 2012
[2]. 含时量子波包法研究反应性散射F+HBr和Br+H_2[D]. 权伟龙. 湘潭大学. 2006
[3]. 多原子体系反应势能面及动力学[D]. 梁景娟. 山东师范大学. 2012
[4]. 重—轻—重反应体系的动力学研究[D]. 王华阳. 山东大学. 2006
标签：化学论文; 动力学论文; 量子论文; 反应动力学论文; 量子效应论文; 分子和原子论文; 碰撞理论论文; 共振效应论文; 力学论文; 分子势能论文; 振动频率论文; 共振频率论文; 碰撞实验论文; 截面数据论文; 科学论文; 原子论文; 科普论文;