【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

MOLPRO学术成果

引言

氟作为卤素中最活跃的元素,在化学工业中,氟广泛应用于半导体、新能源等领域,这对未来的发展尤为关键。因此,F3的三维势能面也受到了很多关注。目前量子化学计算能够为该系统提供足够可靠的数据,势能面(PESs)可以帮助我们了解可能的反应机理和估计反应速率。采用经典轨道法和直接动态电子结构计算相结合的方法研究了F3的势能面。为了准确描述该系统,吉林大学黄旭日/王德权教授团队在本文中推导出了精确三维PESs。

计算方法

本文利用大基组(aVQZ),用从头算法计算了F3体系在MCSCF/MRCI水平下的最低三个绝热势能面,本文所有的计算都在MOLPRO 2012软件包中完成。

本文采用Jacobi坐标(r, R,ǿ),该坐标在我们之前关于三原子系统的工作[3,19 – 23]中可以很好地表示三维PESs。r表示两个氟原子之间的距离,R表示两个氟原子的折算质量到第三个氟原子的距离,ǿ是这两个矢量之间的夹角 (见图1)。

【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

结果与讨论

本文推导了F3(2A’’)体系最低能态的三维绝热势能,为了更清楚地了解该体系的特性,下面给出了二维PES。为了更清楚地讨论这些PESs,将F + F2的能量移至0.00 eV。
F3(2A’’)能量面以不同角度显示了二维最低能态绝热势能面(2D- PESs)。不同的四个角度(ǿ =0o 30 o 60 o 90 o)用来描述PES的生动变化。每张图的上半部分为2D图像PES,每张图的下半部分为2D-PES对应的等高线图。F3(2A’’)线性对称的2D-PES如图2所示。从图中可以看出,势能随着r距离的减小而增加,图中没有最小值。
【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面
如图所示,PES的特点为ǿ = 30o。在这个角度有一个最小值,最小值的几何形状是r = 1.40 Ǻ, R = 3.65 Ǻ。这个最小值的能量为-0.025 eV,略低于反应物的能量。这两个能量的差是如此之小,以至于最小值是不稳定的。如果你想让反应物一起转化成生成物ǿ= 30 o,系统通过TS应该超过0.88 eV能量,TS的几何形状为r = 1.69 Ǻ和R =2.37 Ǻ。
【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面
PES 角度等于 60o 时的特点(见图4)类似于ǿ = 30o。首先,这两个数字都有一个浅的最小值;能量的最小为ǿ = 60o-0.03 eV和这个最小值的几何是r = 1.40 Ǻ和R= 3.11 Ǻ。其次,反应物要通过这两个角度的反应途径转变为生成物,就需要给予系统额外的能量来克服过渡态。TS为ǿ = 60o 的能量是1.255 eV,高于ǿ = 30oTS的几何形状为r = 1.63 Ǻ和r =1.84 Ǻ。

【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

本文也给出了ǿ = 90o 的图表,如图5。当F原子攻击F2分子时ǿ = 90o,对于相同的r,当R>3.5 Ǻ时,能量几乎不随r的函数变化;当R减小到2.95A时, r =1.40 Ǻ时,系统达到浅井复合体,该复合体的能量为-0.026 eV。经过这一复杂过程后,随着R距离的进一步减小,系统的能量越来越高。

【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

F3 (2A’’) 为具有r(Fa Fb)和r(Fb Fc)两个键长的函数的2D-PES,在这两个键之间,固定在α= 180o如图6所示。在同一幅图中把重要区域放大了。蓝色虚线表示最容易的反应途径。同时,在图中绘制了这条虚线的示意图,表明反应势垒为104.715 kJ/mol,与Duggan和Grice的工作(108 kJ/mol)相比,结果几乎一致。反应路径的过渡态几何形状为r(FaFb) = 1.66 Ǻ, r(FbFc) = 1.65Ǻ。

 

总结

通过MOLPRO量子化学高能级从头算(MCSCF/MRCI)方法,得到了F3最低能态的电子绝热态PESs Jacobi坐标。在构型空间的较大区域上,利用aVQZ基组计算势能。通过精确的B-spline拟合方法,计算了27690个势能点,得到了全局精确的PESs。当ǿ = 0.0o时,PES中没有最小值,但是有一个过渡状态。沿角度ǿ = 0.0oFa + FbFc到FaFb + Fc应该给予额外的能量1.089 eV。但对于ǿ = 30o,该能量降低到0.894 eV;当角度增大到60 o时,反能量达到1.255 eV。当ǿ= 30 o 60 o90 o时, PES有浅复合体存在,能量是-0.025 eV, -0.031 eV和-0.026 eV。

【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

是在MOLPRO量化程序包下顺利完成,北京泰科博思科技有限公司是MOLPRO官方指定代理商,更多关于软件咨询或者技术支持,请咨询北京泰科。

作者简介

黄旭日   教授

【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

1981-1991年,本硕博毕业于吉林大学,现任职于吉林大学。研究方向:化学反应的绝热与非绝热势能面的研究;化学反应机理及微观反应动力学的理论研究;药物设计的理论化学计算及计算机模拟。获吉林省科技进步一等奖两次,获吉林省高级专家称号。
 
王德权   博士

【MOLPRO应用实例】吉林大学黄旭日/王德权教授团队成果展示:F3的三维势能面

2013年获德国锡根大学理论化学专业博士学位,现任职于吉林大学。研究方向:小分子反应动力学程序开发与研究。

公司简介

北京泰科博思科技有限公司(Beijing Tech-Box S&T Co. Ltd.)成立于2007年,是国内领先的分子模拟及虚拟仿真综合解决方案提供商。

北京泰科博思科技有限公司与国际领先的模拟软件厂商、开发团队深入合作,为高校、科研院所和企业在材料、化工、药物、生命科学、环境、人工智能及数据挖掘、虚拟仿真教学等领域提供专业的整体解决方案。用户根据需要在我们的平台上高效的进行各种模拟实验,指导实际的生产设计。

北京泰科博思科技有限公司拥有一支一流的技术服务团队和资深的专家咨询团队,以客户真正需求出发,服务客户,为客户创造价值。我们秉承“职业、敬业、担当、拼搏、合作”的企业精神,致力于用国际领先的软件产品和专业全面的技术支持服务,成为客户可信赖的合作伙伴。

评论