非线性三原子分子高激发振动态的理论研究: Lie代数方法

利用Ｌｉｅ代数方法研究了弯曲三原子分子的振动高激发态能谱，并以ＳＯ_２为例， 拟合 ３０条光谱能级得到的 ＲＭＳ误差是 １．６６ ｃｍ~－１．结果表明，所得到的分子 Ｈａｍｉｌｔｏｎｉａｎ 的代数展开式可以很好地再现实验能级，它预测了ＳＯ_２分子振动总量子数达１０的全部 振动能级．从得到的该代数Ｈａｍｉｌｔｏｎｉａｎ还可以得到势能面等．     

Ne-CO2的从头算势能面及微波光谱

采用三重激发校正的耦合簇[CCSD(T)]方法和大基组计算了范德华复合物Ne-CO2的分子间势能面. 分子间相互作用能的计算采用考虑了基组重叠误差修正的超分子方法. 计算结果表明, 该势能面有两个极小值点, 分别对应T形构型和线性Ne-OCO构型. 采用离散变量表象(DVR)方法及Lanczos算法计算了Ne-CO2的振转能级. 计算结果表明, 体系势能面支持22个振动束缚态. 计算得到的微波光谱的跃迁频率与实验值吻合得很好.     

He-N2O的从头算势能面及振转能级

采用超分子MP4方法和较大的基组计算得到了He-N₂O体系的分子间势能面,发现该势能面有3个极小值点,分别对应T形构型及两个线性He-ONN和He-NNO构型.同时采用离散变量表象方法预测了体系的振转能级,计算结果表明,MP4势能面支持5个振动束缚态.     

复合物Kr-CS2的一个新的四维的势能面及其红外光谱的预测

用高精度的量子化学从头算方法构建了Kr-CS₂体系的精确的四维势能面.对该势能面不仅考虑了分子间的振动方式,而且考虑了单体CS₂分子内的ν₁对称伸缩和ν₃反对称伸缩.采用了量子化学超分子方法在[CCSD（T）]-F12水平上构建了从头算势能面.将得到的四维势能面作积分,便得到CS₂处于振动基态和激发态的复合物平均势能面,这两个态的势能面均有一个T型的全局最小值和两个相等的线性极小值.还通过径向部分采用离散变量表象法（DVR）和角度部分采用有限基组表象法（FBR）结合Lanczos循环算法对Kr-CS₂复合物的振转能级进行了计算,对光谱常数进行了预测.另外,也预测了Kr-CS₂复合物的带心位移（-1.2357 cm^-1）.     

基态UC2分子的结构和势能函数

采用密度泛函理论(DFT)的B3LYP方法和相对论有效原子实势理论模型(RECP),对UC2分子可能的结构进行优化计算,得到UC2分子稳定构型为角形C-U-C(C2v);由微观可逆性原理,判断了UC2分子的离解极限;并且导出了基态UC2分子(X 5B1)的多体项展式势能函数,其势能面等值图展现了C-U-C(C2v)稳定结构;根据势能面等值图,讨论了C+UC(X 3П)反应和U+C2(X 1∑+g)反应的势能面静态特征.     

Xe-N2O复合物的分子间势能面和振转光谱的理论研究

采用CCSD(T)方法研究了范德华分子体系Xe—N2O复合物的势能面和振转光谱性质,研究表明,该势能面有两个极小点,分别对应T构型和线性Xe—ONN构型,采用离散变量表象和Lanczos算法计算了体系的振转能级,计算结果表明,CCSD(T)势能面支持97个振动束缚态,并对能级进行了指认,计算得到的Xe—N2O转动跃迁频率与实验值吻合得很好。     

锂(100)面上H_2H+H反应动力学的理论研究

本文在H-Li(100)面吸附扩散的ab initio SCF势能面基础上构造了(H_2H+H)/Li(100)面相互作用的推广LEPS势能面,并用QCT方法研究了该体系的反应动力学行为。分析势能面特征得到:H_2在Li(100)面上的吸附无需活化能,H_2在Li(100)面上的解离吸附与吸附位及吸附模式密切相关,H_2的卧式解离比立式解离要容易得多。分析各种碰撞轨迹得到:低覆盖度下双氢原子的表面复合几率很小,H_2的表面解离几率受到H_2振动量子数的控制。本文构造了一种适合于动力学研究的气体-金属表面相互作用势能面,并且,动力学QCT计算结果能够对H_2表面活化的分子束实验作出合理的解释。     

PuC和PuC2的分子结构与势能函数

采用密度泛函B3LYP方法和相对论有效原子实理论模型优化出PuC和PuC2分子稳定构型,其电子状态分别为X5Σ－和X5A2.PuC2分子为C2v构型,其∠CPuC=147.67°,平衡核间距Re=0.22819 nm, 离解能De=5.543 eV, 并计算出谐振动频率:ν1=61.736 cm－1、ν2=229.894 cm－1、ν3=305.582 cm－1.在此基础上,运用多体项展式理论方法,导出了基态PuC2分子的分析势能函数,该势能面准确地再现了PuC2分子的稳定结构,并根据势能面等值图讨论了PuC＋C反应和Pu＋C2反应的势能面静态特征.     

基于神经网络的分子体系势能面构造

本文介绍了近年来使用神经网络构造分子体系势能面的进展,在以前的研究基础上完善了一套完整的分子构型选择方案,在严格的势能面评价标准下,完全解决了势能面构造过程中分子构型选择的难题.基于此方法,采用一系列拟合技巧,结合大量反应动力学计算,成功构造了一系列重要体系的高精度从头算势能面,并能得到可靠的化学动力学结果.     

Ne-HCl势能面和振转光谱的理论研究

利用量子化学计算方法CCSD(T)和大基组aug-cc-pVTZ加键函数3_s3_p2_d对Ne-HCl体系的分子间势能面进行了理论研究.结果表明,势能面上有两个势阱,分别对应于线性Ne-ClH和Ne-HCl构型.通过精确求解核运动方程发现,该从头算势能面分别支持5个(对Ne-HCl)和7个(Ne-DCl)振动束缚态.计算得到的振转跃迁频率与实值值吻合.     

Quasiclassical trajectory study of the O(3P) + CH4 --> OH + CH3 reaction with a specific reaction parameters semiempirical Hamiltonian

We present a theoretical study of the O(3P) + CH4 --> OH + CH3 reaction using electronic structure, kinetics, and dynamics calculations. We calculate a grid of ab initio points at the PMP2/AUG-cc-pVDZ level to characterize the potential energy surface in regions of up to 1.3 eV above reagents. This grid of ab initio points is used to derive a set of specific reaction parameters (SRP) for the MSINDO semiempirical Hamiltonian. The resulting SRP-MSINDO Hamiltonian improves the quality of the standard Hamiltonian, particularly in regions of the potential energy surface beyond the minimum-energy reaction path. Quasiclassical-trajectory calculations are used to study the reaction dynamics with the original and the improved MSINDO semiempirical Hamiltonians, and a prior surface. The SRP-MSINDO semiempirical Hamiltonian yields OH rotational distributions in agreement with experimental results, improving over the results of the other surfaces. Thermal rate constants estimated with Variational Transition State Theory using the SRP-MSINDO Hamiltonian are also in agreement with experiments. Our results indicate that reparametrized semiempirical Hamiltonians are a good alternative to generating potential energy surfaces for accurate dynamics studies of polyatomic reactions.     

Potential energy surface and highly excited vibrational lines of NO2 via algebraic approach

The algebraic Hamiltonian of NO₂ is optimized using U(4) algebra via fitting to 102 observed vibrational lines. The RMS error of the fitting is 2.39 cm^?1. We calculated highly excited vibrational energy levels using this optimized Hamiltonian, and then obtained the potential energy surface for the electronic ground state by using the classical limit of the U(4) algebraic Hamiltonian. We also calculated the dissociation energies, the force constants etc. Our results are in good agreement with the other theoretical results. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Complete nuclear motion Hamiltonian in the irreducible normal mode tensor operator formalism for the methane molecule

A rovibrational model based on the normal-mode complete nuclear Hamiltonian is applied to methane using our recent potential energy surface [A. V. Nikitin, M. Rey, and Vl. G. Tyuterev, Chem. Phys. Lett. 501, 179 (2011)]. The kinetic energy operator and the potential energy function are expanded in power series to which a new truncation-reduction technique is applied. The vibration-rotation Hamiltonian is transformed systematically to a full symmetrized form using irreducible tensor operators. Each term of the Hamiltonian expansion can be thus cast in the tensor form whatever the order of the development. This allows to take full advantage of the symmetry properties for doubly and triply degenerate vibrations and vibration-rotation states. We apply this model to variational computations of energy levels for (12)CH(4), (13)CH(4), and (12)CD(4).     

Theoretical study of highly vibrational states of nonlinear triatomic molecules using Lie algebraic approach

The vibrational excitations of bent triatomic molecules are studied by using Lie algebra. The RMS error of fitting 30 spectroscopic data is 1.66 cm-1 for SO2. The results show that the expansion of a molecular algebraic Hamiltonian can well describe the experimental data. And the total vibrational levels can be calculated using this Hamiltonian. At the same time, the potential energy surface can also be obtained with the algebraic Hamiltonian.     

Classical trajectory study of the dynamics of the reaction of Cl atoms with ethane

We present an electronic-structure and dynamics study of the Cl + C2H6 --> HCl + C2H5 reaction. The stationary points of the ground-state potential energy surface have been characterized using various electronic-structure methods and basis sets. Our best calculations, CCSD(T) extrapolated to the complete basis limit, using geometries and harmonic frequencies obtained at the MP2/aug-cc-pVTZ level, are in agreement with the experimental reaction energy. Ab initio information has been used to reparameterize a semiempirical Hamiltonian so that the predictions of the improved Hamiltonian agree with the higher-level calculations in key regions of the potential energy surface. The improved semiempirical Hamiltonian is then used to propagate quasiclassical trajectories. Computed kinetic energy release and scattering angle distributions at a collision energy of approximately 5.5 kcal mol(-1) are in reasonable agreement with experiments, but no evidence was found for the low translational energy HCl products scattered in the backward hemisphere reported in recent experiments.     

Semiempirical QM/MM potential with simple valence bond (SVB) for enzyme reactions. Application to the nucleophilic addition reaction in haloalkane dehalogenase

 We present a method for the correction of errors in combined QM/MM calculations using a semiempirical Hamiltonian for enzyme reactions. Since semiempirical models can provide a reasonable representation of the general shape of the potential energy surface for chemical reactions, we introduce a simple valence bond-like (SVB) term to correct the energies at critical points on the potential energy surface. The present SVB term is not a stand-alone potential energy function, but it is used purely for introducing small energy corrections to the semiempirical Hamiltonian to achieve the accuracy needed for modeling enzymatic reactions. We show that the present coupled QM-SVB/MM approach can be parameterized to reproduce experimental and ab initio results for model reactions, and have applied the PM3-SVB/MM potential to the nucleophilic addition reaction in haloalkane dehalogenase. In a preliminary energy minimization study, the PM3-SVB/MM results are reasonable, suggesting that it may be used in free energy simulations to assess enzymatic reaction mechanism. Received: 1 November 2001 / Accepted: 6 September 2002 / Published online: 19 February 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 222nd National Meeting of the American Chemical Society, 2001 Correspondence to: Lakshmi S. Devi-Kesavan e-mail: kesavan@chem.umn.edu Acknowledgments. The work is partially supported by the NIH and the NSF.     

Quasi-classical trajectory study of the dynamics of the Cl + CH4→ HCl + CH3 reaction

We present an on-the-fly classical trajectory study of the Cl + CH(4)→ HCl + CH(3) reaction using a specific reaction parameter (SRP) AM1 Hamiltonian that was previously optimized for the Cl + ethane reaction [S. J. Greaves et al., J. Phys Chem A, 2008, 112, 9387]. The SRP-AM1 Hamiltonian is shown to be a good model for the potential energy surface of the title reaction. Calculated differential cross sections, obtained from trajectories propagated with the SRP-AM1 Hamiltonian compare favourably with experimental results for this system. Analysis of the vibrational modes of the methyl radical shows different scattering distributions for ground and vibrationally excited products.     

A reaction surface Hamiltonian study of malonaldehyde

We report calculations using a reaction surface Hamiltonian for which the vibrations of a molecule are represented by 3N-8 normal coordinates, Q, and two large amplitude motions, s(1) and s(2). The exact form of the kinetic energy operator is derived in these coordinates. The potential surface is first represented as a quadratic in Q, the coefficients of which depend upon the values of s(1),s(2) and then extended to include up to Q(6) diagonal anharmonic terms. The vibrational energy levels are evaluated by solving the variational secular equations, using a basis of products of Hermite polynomials and appropriate functions of s(1),s(2). Our selected example is malonaldehyde (N=9) and we choose as surface parameters two OH distances of the migrating H in the internal hydrogen transfer. The reaction surface Hamiltonian is ideally suited to the study of the kind of tunneling dynamics present in malonaldehyde. Our results are in good agreement with previous calculations of the zero point tunneling splitting and in general agreement with observed data. Interpretation of our two-dimensional reaction surface states suggests that the OH stretching fundamental is incorrectly assigned in the infrared spectrum. This mode appears at a much lower frequency in our calculations due to substantial transition state character.     

Independent particle theory with electron correlation

We formulate an effective independent particle model where the effective Hamiltonian is composed of the Fock operator and a correlation potential. Within the model the kinetic energy and the exchange energy can be expressed exactly leaving the correlation energy functional as the remaining unknown. Our efforts concentrate on finding a correlation potential such that exact ionization potentials and electron affinities can be reproduced as orbital energies. The equation-of-motion coupled-cluster approach enables us to define an effective Hamiltonian from which a correlation potential can be extracted. We also make the connection to electron propagator theory. The disadvantage of the latter is the inherit energy dependence of the potential resulting in a different Hamiltonian for each orbital. Alternatively, the Fock space coupled-cluster approach employs an effective Hamiltonian which is energy independent and universal for all orbitals. A correlation potential is extracted which yields the exact ionization potentials and electron affinities and a set of associated molecular orbitals. We also describe the close relationship to Brueckner theory.     

New method for calculating bound states: the A(1) states of Li(3) on the spin-aligned Li(3)(1 (4)A(')) potential energy surface

In this paper, we present a calculation for the bound states of A(1) symmetry on the spin-aligned Li(3)(1 (4)A(')) potential energy surface. We apply a mixture of discrete variable representation and distributed approximating functional methods to discretize the Hamiltonian. We also introduce a new method that significantly reduces the computational effort needed to determine the lowest eigenvalues and eigenvectors (bound state energies and wave functions of the full Hamiltonian). In our study, we have found the lowest 150 energy bound states converged to less than 0.005% error, and most of the excited energy bound states converged to less than 2.0% error. Furthermore, we have estimated the total number of the A(1) bound states of Li(3) on the spin-aligned Li(3)(1 (4)A(')) potential surface to be 601.