Accurate ab initio-based DMBE potential energy surface for HLi2(X 2A′) via scaling of the external correlation

A globally accurate potential energy surface is reported for the electronic ground-state HLi₂ by fitting ab initio energies to double many-body expansion formalism. The total 3726 ab initio energies used to map the HLi₂ potential energy surface are calculated using the multi-reference configuration interaction method, with their dynamical correlation being semiempirically corrected by the double many-body expansion-scaled external correlation method. The current potential energy surface generates an excellent fit of the ab initio energies, showing a small root-mean squared derivation of 0.636 ? kcal ? mol^-1. The topographical features of the HLi₂ potential energy surface are examined in detail, which concludes that the H + Li₂(X ? ¹ Σ _g ) → Li + LiH(X ? ¹ Σ) reaction is essentially barrierless and the exothermicity is calculated to be 33.668 ? kcal ? mol^-1, thus corroborates the available experimental and theoretical results.     

Globally accurate ab initio based potential energy surface of H_2O~+(X~4A')

A globally accurate potential energy surface is reported for the electronic ground-state H2O+.The ab initio energies utilized to map the potential energy surface are calculated at the multireference configuration interaction method employing the aug-cc-p VQZ basis set and the full valence complete active space wave function as reference.In order to improve accuracy of the resulting raw ab initio energies,they are then extrapolated to the complete basis set limit and most importantly to the full configuration-interaction limit by semiempirically correcting the dynamical correlation using the double manybody expansion-scaled external correlation method.The topographical features of the current potential energy surface were examined in detail,which agree nicely with those of other theoretical work.     

CuTe,Cu2和Cu2Te的结构与势能函数

在Cu和Te的RECP(Relativistic Effective Core Potential)近似下,运用B3LYP方法,在LANL2DZ基组水平上对CuTe,Cu₂和Cu₂Te分子体系的结构进行优化计算.结果表明,CuTe和Cu₂分子的基电子状态分别为²Π和¹∑_g⁺,Cu₂Te分子的基态为单重态的C_2V构型,其电子状态为¹A₁.同时还计算了Cu₂Te分子基态的离解能、力常数和振动频率.采用最小二乘法拟合出CuTe和Cu₂分子Murrell-Sorbie势能函数参数.在此基础上,运用多体展式理论方法导出Cu₂Te分子基态势能函数的解析表达式,其势能面准确复现了平衡态的结构特征.     

UN、UN2分子的结构和势能函数

用相对论有效原子实势(RECP)和密度泛函(B3LYP)方法对UN和UN2分子的结构进行优化,得到了它们的平衡几何构型和谐振频率.采用最小二乘法拟合出UN 分子的Murrell-Sorbie势能函数,在此基础上推导出光谱数据和力常数;并通过多体展式理论导出UN2分子的势能函数,正确地反应了其平衡构型特征.     

AuY、Au2Y分子的结构和势能函数

用密度泛函DFT中的 B3LYP 方法,选择LANL2DZ基组,对AuY、Au2Y分子的结构进行优化,得到了它们的平衡几何构型和谐振频率.采用最小二乘法拟合出AuY分子的 Murrell-Sorbie势能函数,在此基础上推导出光谱数据和力常数;并通过多体展示理论导出Au2Y分子的势能函数,正确地反应了其平衡构型特征.     

Mg2Ni分子的结构和势能函数

采用密度泛函B3LYP方法在6-311+g(d,p)基组水平上对MgNi、Mg2和Mg2Ni分子的各种可能的结构进行优化,得到了它们的几何构型、平衡核间距、离解能和谐振频率.采用最小二乘法拟合MgNi和Mg2分子的Murrell-Sorbie势能函数,在此基础上导出光谱数据和力常数.通过多体项展式理论导出Mg2Ni分子的解析势能函数,并绘出了Mg2Ni分子的等值势能面,其等值势能面正确地反应了基态Mg2Ni分子的平衡构型特征.     

ScH~+和ScH_2~+分子离子的结构与势能函数(英文)

本文利用相对论有效原子实势(RECP)和密度泛函(B3PW91)的方法对ScH 和ScH2 分子离子的结构进行了优化,得到了它们的平衡几何构型和谐振频率.采用最小二乘法拟合出ScH 分子离子的Murrell-Sorbie势能函数,在此基础上,推导出光谱数据和力常数,并通过多体展式理论导出ScH2 分子离子的势能函数,其等值势能面图准确地再现了ScH<,2< 分子离子的结构特征和离解能,由此讨论了反映势能面的静态特征,并利用杂化轨道理论解释了ScH2 分子离子的结构.     

Ar-H2(D2, T2)碰撞体系的振转相互作用势及散射截面的理论计算

用量子力学从头算的耦合族CCSD(T)方法, 使用相关一致基组aug-cc-pV5Z并加3s3p2d1f1g高斯键函数计算了Ar原子与H₂分子的振转相互作用和电荷分布, 采用Boys和Bernardi提出的均衡法消除了基组重叠误差(BSSE). 然后用Tang-Toennies势能函数拟合得到Ar-H₂体系相互作用势的解析表达式. 在该相互作用势下, 用密耦方法计算了Ar原子入射能量为83 meV时, Ar-H₂(D₂, T₂)碰撞体系的散射截面. 计算Ar-D₂体系的微分截面与实验值比较符合很好. 计算结果及分析表明, 在长程吸引势的散射中, 色散能起主要作用; 在短程排斥势的散射中, 交换能起重要作用. 当碰撞参数在0.27 nm至0.47 nm的范围时, Ar-H₂(D₂, T₂)碰撞体系的径向偶极发生两次转向.     

TiO,O2 和TiO2的分析势能函数及光谱研究

用Gaussian09程序包的密度泛函理论DFT方法,在BP86/6-311++g(d,p)水平上对O₂, TiO和TiO₂ 分子进行了优化.得到该系列分子的基态电子态分别为：O₂(X³Σ_g), TiO(X³Π_g), TiO₂(X¹ A₁), TiO₂分子的稳定构型为C_2v构型. 用Murrell-Sorbie势能函数对TiO和O₂分子的扫描势能点进行拟合, 其扫描点都与四参数Murrell-Sorbie函数拟合曲线符合得很好,在此基础上推导出它们的光谱数据和力常数. 用多体项展开理论导出TiO₂分子的全空间解析势能函数,在固定键角∠OTiO=110.5° 的情况下, R_Ti-O = 0.1652 nm处存在一个深度为15.09 eV的势阱, 表明在该处易形成稳定的TiO₂分子. 关键词： TiO 2和TiO₂')" href="#">O₂和TiO₂ 密度泛函理论 势能函数     

Si_3分子基态(X~1A_1)的结构与分析势能函数

应用群论及原子分子反应静力学方法推导Si分子的电子态及其离解极限,在B3P86/CC-PVTZ水平上,对Si3分子基态进行优化计算,得出Si3基态的单重态能量最低,其稳定构性为的C2V构型,平衡核间距Re=0.2176nm、∠213=79.7°,能量为-869.2057a.u..同时计算出基态的简正振动频率:对称伸缩振动频率ν(B2)=547.6446cm-1,弯曲振动频率ν(A1)=185.6100cm-1和反对称伸缩振动频率ν(A1)=559.6090cm-1.在此基础上,使用多体项展式理论方法,导出了基态Si3分子的全空间解析势能函数,该势能函数准确再现了Si3(C2V)平衡结构.     

SeN_2自由基解析势能函数的耦合簇理论研究

利用单双迭代耦合簇理论CCSD结合相关一致四重基组cc-pVQZ对SeN2基态的平衡结构和谐振频率进行了优化计算.计算结果表明:基态SeN2自由基分子稳定态为C2v构型,基态电子组态为X1A1,平衡核间距RSe-N=0.1691 nm,RN-N=0.1970 nm,αN-Se-N=71.289?,离解能De=4.78 eV.基态简正振动频率分别为:ν1=326.9288 cm-1,ν2=808.0161 cm-1以及ν3=948.3430 cm-1.对SeN基态和N2基态采用上述相同方法进行几何构型与谐振频率的计算并进行单点能扫描,使用Murrell-Sorbie函数进行最小二乘拟合得到其势能函数和光谱常数,通过和其他理论值以及实验值做比较,显示本文的计算工作达到了很高的精度.应用多体项展式理论导出了基态SeN2的全空间解析势能函数,其势能函数等值势能图准确再现了SeN2分子的结构特征和能量变化.     

N2O异构体结构与解析势能函数

在B3P86/cc-PVTZ水平上,对N₂O异构体进行优化计算,得出N₂O基态的单重态能量最低,其稳定构型为C_∞v构型,平衡核间距R₁=0.1121nm,R₂=0.1177nm,α＝180°,能量为-185.1188a.u.同时计算出基态的简正振动频率ω₁(Π)=601.5010 cm^{关键词： 异构体 多体项展式理论 解析势能函数}     

GeO分子激发态的电子结构和跃迁性质的组态相互作用方法研究

应用多参考组态相互作用方法计算了GeO分子的第一解离极限(Ge(³P_g)+O(³P_g))对应的18个Λ-S电子态的电子结构. 计算中纳入了Ge原子的3d轨道电子的内壳层-价壳层电子关联效应、标量相对论效应和Davidson修正. 基于计算的电子态的电子结构, 通过求解径向Schrödinger方程获得了束缚电子态的光谱常数R_e, T_e, ω_e, ω_eχ_e, B_e, 理论计算给出的这些电子态的光谱常数与之前的实验结果符合得很好. 计算了电子态的电偶极矩随核间距的变化, 分析了电子态的组态成分的变化对电偶极矩的影响. 计算的势能曲线表明, 激发态A¹Π, 1¹Σ^-, D¹Δ, a³Π, a’³Σ⁺, d³Δ 和 e³Σ^-的绝热激发能密集地分布于26000-37000 cm^-1范围内, 这些密集分布的电子态之间的相互作用对振动波函数有明显扰动作用. 借助于激发态之间的自旋-轨道耦合矩阵元, 阐明了邻近的激发态对A¹Π和a³Π的扰动作用. 基于计算的A¹Π-X¹Σ⁺和A’¹Σ⁺-X¹Σ⁺跃迁的电偶极跃迁矩和Franck-Condon 因子, 给出了A¹Π 和A’¹Σ⁺态的最低的六个振动能级的辐射寿命.     

SiF$lt;sub$gt;2$lt;/sub$gt;($lt;sup$gt;1$lt;/sup$gt;A$lt;sub$gt;1$lt;/sub$gt;)自由基的从头算及势能函数

采用从头算的单双取代的二次组态相互作用方法及耦合簇理论对SiF₂自由基的基态进行结构优化, 发现用单双取代的二次组态相互作用方法配有基组6-311G(2df)计算得到的结构参数、谐振频率、离解能及力常数与实验值最接近并优于文献值. 借助多体项展式理论导出SiF₂自由基的势能函数并绘制了等值势能图. SiF₂自由基对称伸缩振动和旋转势能图显示: 在SiF+F →SiF₂反应通道上有鞍点出现, F原子需要越过4.38 eV的能垒才能生成稳定的SiF₂自由基; 只能通过SiF+F→SiF₂两个等价的通道才能生成稳定的SiF₂ 自由基, 并且该反应是有阈能的吸热反应. 关键词： 2')" href="#">SiF₂ 从头算 势能曲线     

Theoretical study of spin-forbidden cooling transitions of indium hydride using ab initio methods

The feasibility of spin-forbidden cooling of the In H molecule is investigated based on ab initio quantum chemistry calculations. The potential energy curves for the X~1Σ_(0~+)~+, a~3Π_~(0~-), a~3Π_(0~+), a~3Π_1, a~3Π_2, A~1Π_1, 1~3Σ_(0~-)~+, and 1~3Σ_1~+states of In H are obtained based on multi-reference configuration interaction plus the Davidson corrections method. The calculated spectroscopic constants are in good agreement with the available experimental data. In addition, the influences of the active space and spin–orbit coupling effects on the potential energy curves and spectroscopic constants are also studied. For Re of a~3Π_(0~-), a~3Π_(0~+), a~3Π_1, and a~3Π_2 states, the error from large active space is small. The potential energy curve of the A~1Π_1state is not smooth for small active space. The spin–orbit coupling effects have great influences on the potential well depth and equilibrium internuclear distance of the A~1Π state. The Franck–Condon factors and radiative lifetimes are obtained on the basis of the transition dipole moments of the a~3Π_(0~+) → X~1Σ_(0~+)~+, a~3Π_1 → X~1Σ_(0~+)~+, and A~1Π_1 → X~1Σ_(0~+)~+ transitions. Our calculation indicates that the a~3Π_1( ν'= 0) → X~1Σ_(0~+)~+(ν = 0) transition provides a highly diagonally distributed Franck–Condon factor and a short radiative lifetime for the a3Π1 state, which can ensure rapid and efficient laser cooling of In H.The proposed laser drives a~3Π_1 → X~1Σ_(0~+)~+ transitions by using three wavelengths.     

反应物NO的振动激发对反应C(3P)+NO(X2Π)→CO(X1Σ+)+N(2D)立体动力学性质的影响

采用准经典轨线方法研究碰撞能为0.23 eV时,反应物分子NO在不同初始振动态（v=0~3）下发生在两个电子态（²A″和²A'）势能面上反应C（³P）+NO（X²Π）→CO（X¹Σ⁺）+N（²D）的立体动力学性质。计算反应产物的转动角动量矢量分布（P（θ_r）和P（φ_r））以及微分散射截面（P₀₀（ω_t）,P₂₀（ω_t））。结果表明：反应物的振动激发对发生在²A'势能面上的反应立体动力学性质影响显著,而对发生在²A″势能面上的反应立体动力学性质影响较小,这可能与两个势能面的构型有关。在两个电子态势能面上发生的反应均受内平面反应机理支配。     

Energy and rotation-dependent stereodynamics of H(~2S) + NH(a~1?) → H_2(X~1Σ_g~+) + N(~2D) reaction

Quasi-classical trajectory calculations are performed to study the stereodynamics of the H(~2S) + NH(a~1?) →H_2(X~1Σ_g~+) + N(~2D) reaction based on the first excited state NH_2(1~2A') potential energy surface reported by Li et al.[Li Y Q and Varandas A J C 2010 J. Phys. Chem. A 114 9644] for the first time. We observe the changes of differential cross-sections at different collision energies and different initial reagent rotational excitations. The influence of collision energy on the k–k' distribution can be attributed to a purely impulsive effect. Initial reagent rotational excitation transforms the reaction mechanism from insertion to abstraction. The effect of initial reagent rotational excitations on k–k' distribution can be explained by the rotational excitation enlarging the rotational rate of reagent NH in the entrance channel to reduce the probability of collision between incidence H atom and H atom of target molecular. We also investigate the changes of vector correlations and find that the rotational angular momentum vector j' of the product H_2 is not only aligned, but also oriented along the y axis. The alignment parameter, the disposal of total angular momentum and the reaction mechanism are all analyzed carefully to explain the polarization behavior of the product rotational angular moment.     

Theoretical study of HOCl adsorption on ice surface

Adsorption of HOCl on ice surface was studied using the ab initio molecular orbtial theory. We applied Hartree–Fock (HF) self-consistent field and the second-order Møller–Plesset (MP2) level of theory to cluster models of the (0001) surface ice Ih to optimize adsorption structures and binding energies. In all stable binding configurations, HOCl acts as a proton donor in a hydrogen bond. The presence of neighboring water molecules can strengthen the interaction of HOCl with ice. In the HOCl·(H₂O)₄ system, interaction hydrogen bond length is about 1.85 Å, and binding energies are −10.063−11.149 kcal mol⁻¹. We also calculated the vibrational frequencies of HOCl affected by the ice surface.     

激光冷却OH分子的理论研究

采用高精度的多参考组态相互作用方法计算OH分子基态和第一激发态的势能曲线.为获得更精确的计算结果,在计算过程中考虑Davidson修正、标量相对论效应、核价相关效应和自旋轨道耦合效应.基于计算的Λ-S和?态的势能曲线,对一维径向薛定谔方程进行数值求解,得到各个电子态的光谱数据,与已报道的实验值和理论值相符合.获得OH分子的永久偶极矩、跃迁偶极矩、振动能级、Franck-Condon因子及辐射寿命,结果表明,A~2Σ~+→X~2Π跃迁具有高度对角化的Franck-Condon因子(0.9053)和短的辐射寿命(5.8363×10~(-7)s),符合激光冷却分子的条件.制定了激光冷却OH分子的具体方案,计算得到激光冷却跃迁A~2Σ~+→X~2Π所需的三束激光波长,主光束波长为307.1532 nm,两束重抽运激光波长为344.9163和349.7659 nm.计算结果为超冷OH分子的实验制备提供重要的理论依据.     

基于变分Monte Carno计算的H2分子势能曲线的可靠性

基于Alexander and Coldwell报道的变分Monte Carno计算的势能曲线,计算了H2分子的振转能级.利用Murrell-Sorbie函数和最小二乘法拟合出了解析势能函数,以此为基础计算出光谱常数.通过与实验比较发现:由于该势能曲线考虑了相对论效应和非绝热效应,振转能级和光谱常数计算结果与实验值符合良好,与目前最高水平从头计算的结果相当,说明这种计算方法准确可靠,可进一步推广应用于势能函数计算之中.