Three-dimensional Potential Energy Surface and Bound States of the Ar2-Ne Complex

The first three-dimensional interaction potential energy surface (PES) of the Ar₂-Ne complex is developed using the single and double excitation coupled cluster theory with noniterative treatment of triple excitations CCSD(T). The aug-cc-pVQZ basis sets are employed for all atoms, including an additional (3s3p2d2f1g) set of midpoint bond functions. The calculated single point energies are fitted to an analytic two-dimensional potential model at each of seven fixed r_Ar2 values. The seven model potentials are then used to construct the three-dimensional PES by interpolating along (r—r_e) using a sixth-order polynomial. The PES is used in the following rovibrational energy levels calculations. The comparisons of theoretical transition frequencies and spectroscopic constants with the experimental results are given.     

范德华复合物Kr-C2H2的势能面和振转光谱的理论研究

采用[CCSD（T）]-F12方法和aug-cc-pVTZ基组,同时引入中心键函数（3s3p2d1f1g）构建了Kr-C₂H₂体系的高精度四维势能面.在构建势能面时考虑了分子间的振动方式及C₂H₂单体内的ν₁对称伸缩和ν₃反对称伸缩振动.将计算得到的四维势能面在Q₁方向和Q₃方向分别做积分得到C₂H₂单体分别处于振动基态和（ν₁,ν₃）=（1,1）激发态的平均势能面.计算结果表明,这2个平均势能面均存在2个等价的T型全局极小值和2个等价线性极小值.全局极小值的几何构型位于R=0.41 nm,θ=65.6°/114.4°,势阱深度为151.88 cm^-1.对径向部分采用离散变量表象法（DVR）,角度部分采用有限基组表象法（FBR）,并结合Lanczos循环算法计算了Kr-C₂H₂的振转能级和束缚态.计算结果表明,复合物在（ν₁,ν₃）=（1,1）区域的带心位移为-1.48 cm^-1,表现为红移,与实验值-1.38 cm^-1很接近;计算得到的红外跃迁频率也与实验值相吻合,说明得到的从头算势能面具有高精度.     

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

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

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

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

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

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

An ab initio potential energy surface and vibrational energy levels of HXeBr

A three-dimensional global potential energy surface for the electronic ground state of HXeBr molecule is constructed from more than 4200 ab initio points. These points are generated using an internally contracted multi-reference configuration interaction method with the Davidson correction （icMRCI ＋ Q） and large basis sets. The stabilities and dissociation barriers are identified from the potential energy surfaces. The three-body dissociation channel is found to be the dominate dissociation channel for HXeBr. Based on the obtained potentials, low-lying vibrational energy levels of HXeBr calculated using the Lanczos algorithm is found to be in good agreement with the available experimental band origins.     

Product Vibrational State Distributions of F+CH3OH Reaction on Full-Dimensional Accurate Potential Energy Surface

The hydrogen abstraction reaction of methanol with fluorine atoms can produce HF and CH\begin{document}$ _3 $\end{document}O or CH\begin{document}$ _2 $\end{document}OH radicals, which are important in the environment, combustion, radiation, and interstellar chemistry. In this work, the dynamics of this typical reaction is investigated by the quasi-classical trajectory method based on a recently developed globally accurate full-dimensional potential energy surface. Particularly, the vibrational state distributions of the polyatomic products CH\begin{document}$ _3 $\end{document}O and CH\begin{document}$ _2 $\end{document}OH are determined by using the normal mode analysis method. It is found that CH\begin{document}$ _3 $\end{document}O and CH\begin{document}$ _2 $\end{document}OH are dominantly populated in the ground state when the reactants are at the ground ro-vibrational state. The OH stretching mode, torsional mode, H\begin{document}$ _2 $\end{document}CO out-of-plane bending mode and their combination bands in the CH\begin{document}$ _2 $\end{document}OH product can be effectively excited once the OH stretching mode of the reactant CH\begin{document}$ _3 $\end{document}OH is excited to the first vibrationally excited state. Most of the available energy flows into the HF vibrational energy and the translational energy in both channels, while the radical products, CH\begin{document}$ _3 $\end{document}O or CH\begin{document}$ _2 $\end{document}OH, receive a small amount of energy, consistent with experiment, which is an indication of its spectator nature.     

SO2基态势能面和振转能级的理论研究

在键长-键角坐标下,精确求解SO₂的核振动方程,并通过与实验数据比较来优化势能参数,由所得势能面计算得到38个振动能级,与实验值相比,均方根误差为0.93cm^-1.计算了³⁴SO₂的部分振动能级以及³²SO₂的J=6以下的部分振动能级,所得结果均与实验值较为吻合     

Ab initio electronic and rovibrational structure of

The electronic structure of the ground state of has been investigated using relativistically-corrected CCSD(T) in conjunction with ANO-RCC (Mg) and aug-cc-pVQZ (H) basis sets. The molecular potential energy surface possessed minima corresponding to both ¹A₁ and equilibrium structures (with a ¹Σ⁺ transition state). The ¹A₁ structure possessed R_e and θ_e values of 2.0297 Å and of 22.09°, respectively. The higher-energy structure exhibited an R_e value of 2.1658 Å. Property surfaces were constructed to calculate rovibrational energies and spectral line intensities for the ground states of , (¹A′) MgHD²⁺ and . For the vibration ground state of , the vibration-averaged R_e and θ_e values were calculated to be 2.0209 Å and 22.53°, respectively. The A, B and C rotational constants were calculated to be 58.0, 2.21 and 2.11 cm⁻¹, respectively.     

The Role of High Excitations in Constructing Sub-spectroscopic Accuracy Intermolecular Potential of He-HCN: Critically Examined by the High-Resolution Spectra with Resonance States

Interpreting high-resolution rovibrational spectra of weakly bound complexes commonly requires spectroscopic accuracy (<1 cm^-1) potential energy surfaces (PES). Constructing high-accuracy ab initio PES relies on the high-level electronic structure approaches and the accurate physical models to represent the potentials. The coupled cluster approaches including single and double excitations with a perturbational estimate of triple excitations (CCSD(T)) have been termed the "gold standard" of electronic structure theory, and widely used in generating intermolecular interaction energies for most van der Waals complexes. However, for HCN-He complex, the observed millimeter-wave spectroscopy with high-excited resonance states has not been assigned and interpreted even on the ab initio PES computed at CCSD(T) level of theory with the complete basis set (CBS) limit. In this work, an effective three-dimensional ab initio PES for HCN-He, which explicitly incorporates dependence on the Q₁ (C-H) normal-mode coordinate of the HCN monomer has been calculated at the CCSD(T)/CBS level. The post-CCSD(T) interaction energy has been examined and included in our PES. Analytic two-dimensional PESs are obtained by least-squares fitting vibrationally averaged interaction energies for v₁(C-H)=0, and 1 to the Morse/Long-Range potential function form with root-mean-square deviations (RMSD) smaller than 0.011 cm^-1. The role and significance of the post-CCSD(T) interaction energy contribution are clearly illustrated by comparison with the predicted rovibrational energy levels. With or without post-CCSD(T) corrections, the value of dissociation limit (D₀) is 8.919 or 9.403 cm^-1, respectively. The predicted millimeter-wave transitions and intensities from the PES with post-CCSD(T) excitation corrections are in good agreement with the available experimental data with RMS discrepancy of 0.072 cm^-1. Moreover, the infrared spectrum for HCN-He complex is predicted for the first time. These results will serve as a good starting point and provide reliable guidance for future infrared studies of HCN doped in (He)_n clusters.     

A Theoretical Study on the Potential Energy Surface of the C3 + NO Reaction

The reaction of linear form carbon cluster C₃ molecules in their ¹Σ ground state with NO (X²Π) radicals is explored theoretically to investigate the formation of hitherto undetected NCCCO molecules in the interstellar medium via a neutral-neutral reaction. The doublet potential energy surface is worked out by the ab initio MO calculations at the CCSD(T)/cc-pVTZ//B3LYP/6-311G(d,p)+ZPE level of theory. It is shown that the main pathway of the C₃(¹Σ)+NO(X²Π) reaction involves the N-atom of NO attacking the side C-atom of the ¹C₃ molecule first to form the adduct CCCNO, followed by the N-shift to give I6 CCNCO, and then to the main products P₁ (CCN+CO). Alternatively, I6 can be converted via the N-shift again to I9 (CN)CCO, and then it leads to the minor products P₂ (CNC+CO) and P₃ NCCCO. Since the three pathways have zero threshold energy and proceed via strongly bound energized complexes, they should possess the character of negative temperature dependence, and might be quite rapid even at very low temperature. The reaction represents facile neutral-neutral pathways to produce hitherto undetected CCN, CNC and NCCCO molecules in interstellar environments.     

A fully ab initio potential energy surface for ClH2 reactive system

Anab initio analytical potential energy surface called BW3 for the ClH₂ reactive system is presented. The fit of this surface is based on about 1 200ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H₂ and HCl, and the endothermicity of the Cl + H₂ abstraction reaction. For the Cl + H₂ abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

New ab initio Potential Energy Surfaces for Cl（^2P3/2,^2P1/2）＋H2 Reaction

New global three dimensional potential energy surfaces for the Cl＋H2 reactive system have been constructed using accurate multireference configuration interaction calculations with a large basis set. The three lowest adiabatic potential energy surfaces correlating asymptotically with Cl（^2p）＋H2 have been transformed to adiabatic representation, which leads to a fourth coupling potential for non-linear geometries. In addition, the spin-orbit coupling surfaces have also been computed using the Breit-Pauli Hamiltonian. Properties of the new potential are described. Reaction dynamics based on the new potential agrees with the recent experimental results quite well.     

次氟酸分子的势能函数及其振动光谱的同位素效应的理论研究

用三原子振动激发态的变分计算程序(TRIATOM)精确计算次氟酸分子H¹⁶OF的振动激发态的能级以及次氯酸分子中的H和O分别被D和¹⁸O取代后的H¹⁸OF,D¹⁶)OF和D¹⁸OF的同位素效应,理论计算值与已有的实验结果吻合较好。预测了一些尚未观测到的谱线频率及同位素效应,并确立了一个同位素位移的加和规则。     

Ab Initio Study on the Potential Energy Surfaces of B4 Cluster

The singlet and triplet potential energy surfaces (PES) for the isomerization and dissociation reactions of B₄ isomers have been investigated using ab initio methods. Ten B4 isomers have been identified and of these 10 species, 4 have not been reported previously. The singlet rhombic structure ¹1 is found to be the most stable on the B₄ surface, in agreement with the results of previous reports. Several isomerization and dissociation pathways have been found. On the singlet PES, the linear ¹3b can rearrange to rhombus ¹1 directly, while ¹3c rearranges to ¹1 through two‐step reactions involving a cyclic intermediate. On the triplet PES, the capped triangle structure ³2 undergoes ring opening to the linear isomer ³3b with a barrier of 34.8 kcal/mol and 44.9 kcal/mol, and the latter undergoes ring closure to the square structure ³1 with a barrier of 30.4 kcal/mol and 33.0 kcal/mol at the MP4/6–311+G(3df)//MP2/6–311G(d) and CCSD/aug‐cc‐pVTZ//MP2/6–311G(d) levels of theory, respectively. The direct decomposition of singlet B₄ yielding to B₃+B is shown to have a large endothermicity of 87.3 kcal/mol (CCSD), and that producing 2B₂ to have activation energy of 133.4 kcal/mol (CCSD).     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

Analysis of Potential Energy Surface for Butanone Isomerization

The potential energy surfaces for butanone isomerization have been investigated by density function theory calculation. Six main reaction pathways are confirmed using the intrinsic reaction coordinate method, and the corresponding isomerization products are 1-buten-2-ol, 2-buten-2-ol, butanal or 1-buten-l-ol, methyl 1-propenyl ether, methyl allyl ether, and ethyl vinyl ether, respectively. Among them, there are three pathways through butylene oxide, indicating butylene oxide is an important intermediate product during butanone isomer ization. The calculated vertical ionization energies of the reactant and its products are in a good agreement with the experimental values available. From the consideration for the relative energies Of transition states and the number of high-energy barriers we infer that the reaction pathway butanone-＊l-buten-2-ol---2-buten-2-oi is the most competitive. The obtained results are informative for future studies on isomerization of ketone molecules.     

乙烯酮自由基(HCCO&#183;)与氧气(O2)反应势能面的理论研究

在G2(B3LYP/MP2/CC)水平上对反应HCCO+O2进行了计算,得到了反应势能面,提出了3种可能的反应机理:(1)四元环反应机理得到产物P1(HCO+CO2);(2)三元环反应机理得到产物P2(CO+HCO2);(3)O—O键断裂反应机理得到产物P3(O+OCC(O)H)和P4(O+CO+HCO).由反应势能面推测产物P1(HCO+CO2)为主要产物,产物P2(CO+HCO2),P3(O+OCC(O)H)和P4(O+CO+HCO)为次要产物.     

Limiting high pressure rate coefficients for the HO+NO2&→HONO2 reaction on ab initio potential energy surfaces

Phase space limiting high pressure rate coefficients for the title reaction on ab initio potential energy surfaces have been calculated at 50-600 K. Calculated rigidity factors at different levels of theory are presented. The best limiting high pressure rate coefficient obtained at 300 K, 4.0x10^-11cm³ molecule^-1 s^-1, compares very well with the latest IUPAC recommended value. This revised version was published online in June 2006 with corrections to the Cover Date.     

1CH2＋N2O反应的势能面

利用密度泛函理论（B3LYP）计算了1CH2＋N2O反应的反应物、中间体、过渡态及产物 的几何构型.进而用从头算方法（QCISD（T））计算了单点能量.由此描绘了反应的势能面, 确定了反应的最终产物通道为N2＋H2CO和NO＋HCN＋H.后者比前者有更大的分支比.N2、H2CO 、NO、HCN的存在有待于实验检测.作者认为,反应在室温下是加成-消除机理,而在高温下 可以通过直接取代的机理获得N2＋H2CO.