用神经网络方法拟合的反应体系H＋CH4←→H2＋CH3的一个全域势能面

利用神经网络力法,基于47783个高精度从头算能量点构建了反应体系H＋CH4←→H2＋CH3的一个全域势能面．通过大最的准经典轨线以及量子散射计算测试了势能面的收敛性质．这个势能面对于拟合过程以及从头算点的数目都已经完全收敛,拟合误差很小县比Shepard插值的势能面计算速度更快,代表了此标志性多原子反应体系最好的势能面．     

OCS电子基态势能面与振动光谱的理论研究

本文采用键长-键角内标系下的自洽场-组态相互作用方法精确计算了OCS分子的振动高激发态能级,并结合实验观测到的振动能级利用非线性最小二乘法优化电子基态势能函数中的势能参数。由优化所得的势能面计算出的振动激发态能级与50个实验观测到的振动能级比较,标准偏差为0.08cm^-^1。此外,还用该势能面计算了OCS同位素分子的振动能级,计算结果与实验值也十分吻合。     

势能面特征与过渡态熵及A因子

反应体系的势能面,对了解反应的微观过程起着重要的作用,它的特征决定了化学反应的机理.原则上,由反应体系的Schrodinger方程的解,可得到体系能量随核间距变化的函数,从而获得势能面.除少数简单反应外,几乎无法精确求得复杂反应体系的势能面.因而,除从头算法外,人们先后发展了计算势能面的一些半经验方法.对某一反应,文献中可能记载好几个势能面,因此,在分析反应或计算反应的各物理量时,应当说明所应用的是何种势能面.我们曾指出,过渡态熵的可靠性,有赖于提供过渡态参数的势能面.本文从下述基元反应     

构造SO~2分子势能面的李代数方法

把李代数方法得到的SO~2分子的代数Hamiltonian,利用相干态基经典化并找到一个新的变换,将分子的键角引入,而得到SO~2分子的势能面。由该势能面计算的解离能,所给出的势能面的立体图和相应的等高线以及力常数与其他方法给出的相一致。该方法可以推广到多原子分子及反应体系。     

势能面拓扑结构失稳与A+B2反应途径分叉条件

研究势能面拓扑结构对认识势能面的复杂性，分析势能面的关键区域，简化量化计算以及了解反应机理等都具有极其重要的意义＊．近些年来，关于势能面反应途径分又研究表明：势能面拓扑结构失稳会导致反应途径分叉，从而使得分子体系出现一系列新的现象和规律［’－’］．根本上来说，势能面拓扑结构最终由临界点的拓扑性质以及其数目所决定［‘1．因而研究势能面临界点的种类及其变化有至关重要的作用·本文用微分动力学系统定性方法分析了势能面上临界点的拓扑结构类型，并给出平面线性线系统势能面拓扑结构失稳的条件；针对三原子A＋BZ反…     

两态反应Ni2+与c-C6H12的机理及自旋-轨道耦合研究

在密度泛函理论的B3LYP方法下, 对两态反应Ni₂⁺与环己烷体系进行了较为系统的研究. 结果表明, 反应分别在第一个氢迁移(⁴IM1→⁴TS1/2), 第三个氢迁移(⁴TS15/16→⁴IM15)以及 翻转过程(⁴IM5→⁴TS5/6, ²TS11/12→²IM12)发生了二、四重态势能面的交叉, 本文运用内禀反应坐标单点垂直激发计算的方法得到势能面大致的交叉点(CP), 进一步利用Crossing2004程序获得精确的最低能量交叉点(MECP). 对MECP附近的自旋轨道耦合(SOC)常数进行了计算. MECP1～MECP4处的SOC值分别为318.01, 396.89, 268.74和306.67 cm^-1. 较大的SOC值说明不同势能面间发生了有效地跃迁并使反应沿着最低反应通道进行. 对反应通道的研究发现, 反应中同面脱氢是主要反应通道. 异面脱氢由于翻转过程中决速步骤势垒为33 kcal/mol(吸热3 kcal/mol), 只生成少量的异面脱氢产物. 计算结果解释了实验现象.     

H3反应势能面的构建——计算化学实验设计

介绍一个针对高年级本科生的物理化学探索性实验。通过计算化学手段构建H3反应势能面,使学生初步掌握Gaussian03W,Gaussview5．0和Origin软件的使用,深入理解反应过渡态理论,并进一步了解目前势能面的研究动态。     

多原子反应体系的高精度拟合势能面

本文介绍了近几年来我们组构建多原子反应体系的高精度拟合势能面的进展。我们基于神经网络(NN)方法,成功构建了多原子气相分子体系和气相分子在金属表面解离的一系列势能面。这些势能面的拟合精度相当高,并且经过了严格的量子动力学测试,能广泛应用到动力学研究中。我们还提出了一种新的置换不变势能面的拟合方法,即基本不变量神经网络方法(FI-NN)。基本不变量的使用极大地减少了神经网络输入层多项式的个数,有效提高了势能面的计算速度。     

偏分势能面及其应用(Ⅰ):概念、构造方法及几个简例

提出了偏分势能面的概念,偏分势能面可由完全势能面抽取出来,也可采用ab initio方法进行构造．作为范例。给出了F H2→FH H,H H2→H2 H,I HI→IH I及Na I2→Na^ I2^-等体系中几种偏分势能面的构造和应用．可以看到,应用偏分势能面对于分析反应机理及散射共振态的生成方面显示突出优点．     

偏分势能面及其应用(I):概念、构造方法及几个简例

提出了偏分势能面的概念,偏分势能面可由完全势能面抽取出来,也可采用abinitio方法进行构造.作为范例,给出了F+H2→FH+H,H+H2→H2+H,I+HI→IH+I及Na+I2→Na++I-等体系中几种偏分势能面的构造和应用.2可以看到,应用偏分势能面对于分析反应机理及散射共振态的生成方面显示突出优点.     

A full-dimensional coupled-surface study of the photodissociation dynamics of ammonia using the multiconfiguration time-dependent Hartree method

Full-dimensional quantum mechanical computations are carried out to investigate the photodissociation dynamics of A? state NH(3) and ND(3) using the multiconfiguration time-dependent Hartree (MCTDH) method with recently developed coupled ab initio potential energy surfaces (PESs) [Z. H. Li, R. Valero, and D. G. Truhlar, Theor. Chim. Acc. 118, 9 (2007)]. To use the MCTDH method efficiently the PESs are represented as based on the high-dimensional model representation. The A? ← X? absorption spectra for both isotopomers were calculated for the zeroth vibrational state of the ground electronic state. With a view to treating larger systems, Jacobi coordinates are used. Computations on the coupled PES are carried out for two-, three-, five-, and six-dimensional model systems to understand the validity of reduced-dimensional calculations. In addition to the fully coupled calculations, the effect of nonadiabatic coupling on absorption spectra is shown by propagating the initial wavepacket only in the A? electronic state. The calculated absorption spectra are shown to be in good agreement with available theoretical and experimental observations. Comparisons with calculations using Radau and valence coordinates show the effect of including the symmetry of the system explicitly. Finally, branching ratios for loss of a hydrogen atom via the two available channels are calculated. These predict that the nonadiabatic product increases with the dimension of the calculations and confirm the importance of the full-dimensional calculations.     

Theoretical studies on the photochemical reaction mechanisms of o-xylylene. Effects of phenyl group for electrocyclic reaction mechanism

The potential energy surfaces (PESs) of the electrocyclic reactions of o-xylylene at the ground and the lowest excited states are calculated by CASSCF molecular orbital and MRMP2 methods. The lowest excited state geometry of o-xylylene has C(2v) symmetry and is about 65 kcal mol(-1) in energy above the ground state. The PESs in the vicinity of the conical intersection are different from those of the electrocyclic reaction of cis-butadiene. In the vicinity of the conical intersection, the transition state at the ground state relating to methylene-cycloheptadienyl carbene is located. The transition state is only 4.3 kcal mol(-1) lower in energy than the conical intersection at the CASSCF(10,10)/6-31G(d) level and 0.5 kcal mol(-1) lower at the MRMP2/6-311+G(d,p) level. The transition state corresponding to benzocyclobutene does not locate in the vicinity of the conical intersection because of the resonance energy between benzene ring and methylene group.     

High-dimensional ab initio potential energy surfaces for reaction dynamics calculations

There has been great progress in the development of potential energy surfaces (PESs) for reaction dynamics that are fits to ab initio energies. The fitting techniques described here explicitly represent the invariance of the PES with respect to all permutations of like atoms. A review of a subset of dynamics calculations using such PESs (currently 16 such PESs exist) is then given. Bimolecular reactions of current interest to the community, namely, H + CH(4) and F + CH(4), are focused on. Unimolecular reactions are then reviewed, with a focus on the photodissociation dynamics of H(2)CO and CH(3)CHO, where so-called "roaming" pathways have been discovered. The challenges for electronically non-adiabatic reactions, and associated PESs, are presented with a focus on the OH* + H(2) reaction. Finally, some thoughts on future directions and challenges are given.     

Benchmarking the Polyatomic Reaction Dynamics of X+Methane

With recent developments of sophisticated experimental techniques and advanced theoretical methods/computations, the field of chemical dynamics has reached the point that theoryexperiment comparisons can be made at a quantitative level in very fine details for a prototypical A+BC system. As the system becomes larger, more degrees of freedom are involved and the complexity increases exponentially. At the same time, the multifaceted nature of polyatomic systems also opens up the possibilities for observing many new chemistry and novel phenomena|a land of opportunities. For the past 15 years or so my laboratory has delved into the reaction dynamics of methane+X (X: F, Cl, O(³P), and OH). This effort shifts the paradigm in the field of reaction dynamics by making the title reaction a benchmark polyatomic system. In this account, I shall disclose my thinking behind some of the key concepts and methods we introduced and how the unexpectedly discovered phenomena led to other uncharted territories. Those ndings not only enrich our understanding of the specific reactions we studied at the most fundamental level and inspire the theoretical developments, but also shape our thinking and lay the foundation for future explorations of different aspects of the multifaceted nature of polyatomic reactivity.     

用交叉束激光诱导荧光技术研究氟原子与碘化物的反应

在570-620 nm波长范围内用交叉束激光诱导荧光技术研究了F原子与C_2H_5I, C_2F_5I和C_3H_5I的反应, 获得了反应产物IF的振功和转功激发谱; 测得了反应产物IF在v=3, 4和5能级中的初生相对振动布居密度, 和部分的相对细致振动速率常数, 转动温度和振功能级中转动能量的平均部分数; 并对反应机制进行了讨论。     

Spin-orbit coupled potential energy surfaces and properties using effective relativistic coupling by asymptotic representation

A new method has been reported recently [H. Ndome, R. Welsch, and W. Eisfeld, J. Chem. Phys. 136, 034103 (2012)] that allows the efficient generation of fully coupled potential energy surfaces (PESs) including derivative and spin-orbit (SO) coupling. The method is based on the diabatic asymptotic representation of the molecular fine structure states and an effective relativistic coupling operator and therefore is called effective relativistic coupling by asymptotic representation (ERCAR). The resulting diabatic spin-orbit coupling matrix is constant and the geometry dependence of the coupling between the eigenstates is accounted for by the diabatization. This approach allows to generate an analytical model for the fully coupled PESs without performing any ab initio SO calculations (except perhaps for the atoms) and thus is very efficient. In the present work, we study the performance of this new method for the example of hydrogen iodide as a well-established test case. Details of the diabatization and the accuracy of the results are investigated in comparison to reference ab initio calculations. The energies of the adiabatic fine structure states are reproduced in excellent agreement with reference ab initio data. It is shown that the accuracy of the ERCAR approach mainly depends on the quality of the underlying ab initio data. This is also the case for dissociation and vibrational level energies, which are influenced by the SO coupling. A method is presented how one-electron operators and the corresponding properties can be evaluated in the framework of the ERCAR approach. This allows the computation of dipole and transition moments of the fine structure states in good agreement with ab initio data. The new method is shown to be very promising for the construction of fully coupled PESs for more complex polyatomic systems to be used in quantum dynamics studies.     

Ab initio ground and excited state potential energy surfaces for NO-Kr complex and dynamics of Kr solids with NO impurity

The intermolecular potentials for the NO(X 2Pi)-Kr and NO(A 2Sigma+)-Kr systems have been calculated using highly accurate ab initio calculations. The spin-restricted coupled cluster method for the ground 1 2A' state [NO(X 2Pi)-Kr] and the multireference singles and doubles configuration interaction method for the excited 2 2A' state [NO(A 2Sigma+)-Kr], respectively, were used. The potential energy surfaces (PESs) show two linear wells and one that is almost in the perpendicular position. An analytical representation of the PESs has been constructed for the triatomic systems and used to carry out molecular dynamics (MD) simulations of the NO-doped krypton matrix response after excitation of NO. MD results are shown comparatively for three sets of potentials: (1) anisotropic ab initio potentials [NO molecule direction fixed during the dynamics and considered as a point (its center of mass)], (2) isotropic ab initio potentials (isotropic part in a Legendre polynomial expansion of the PESs), and (3) fitted Kr-NO potentials to the spectroscopic data. An important finding of this work is that the anisotropic and isotropic ab initio potentials calculated for the Kr-NO triatomic system are not suitable for describing the dynamics of structural relaxation upon Rydberg excitation of a NO impurity in the crystal. However, the isotropic ab initio potential in the ground state almost overlaps the published experimental potential, being almost independent of the angle asymmetry. This fact is also manifested in the radial distribution function around NO. However, in the case of the excited state the isotropic ab initio potential differs from the fitted potentials, which indicates that the Kr-NO interaction in the matrix is quite different because of the presence of the surrounding Kr atoms acting on the NO molecule. MD simulations for isotropic potentials reasonably reproduce the experimental observables for the femtosecond response and the bubble size but do not match spectroscopic results. A general overall view of the results suggests that, when the Kr-NO interaction takes place inside the matrix, potentials are rather symmetric and less repulsive than those for the triatomic system.