三原子分子振转激发态成簇现象理论研究

在只考虑弯曲振动与总角动量的耦合,而冻结伸缩振动的模型下,采用Jacobi多项式作为弯曲振动的基函数,用严格的变分法研究了H_2O分子的振转激发态的成簇现象.本文计算了H_2O分子的振转能级和函数,研究了振动激发态下转动高激发态光谱中出现的成簇态.     

三原子分子振转激发态成簇现象理论研究

在只考虑弯曲振动与总角动量的耦合,而冻结伸缩振动的模型下,采用Jacobi多项式作为弯曲振动的基函数,用严格的变分法研究了H~2O分子的振转激发态的成簇现象。本文计算了H~2O分子的振转能级和函数,研究了振动激发态下转动高激发态中出现的成簇态。     

CO2势能面和振动激发态的理论研究

采用振动自洽场-组态相互作用(SCF-CI)方法通过实验振动光谱优化了CO₂分子的势能函数,由该势能函数计算得到的纯振动光谱数据与实验值相比,所有能级的误差均在4cm^-1以内,均方根偏差为1.50cm^-1,所预测的Π态振转光谱也与实验值很接近.     

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

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

CH_3F-Ar体系的高精度势能面及预测的振转光谱

本文在从头算CCSD(T)/AVQZ水平上计算了CH_3F-Ar体系包含CH_3F分子Q_3振动的四维势能面.以摩斯长程势(MLR)为模型,分析的三维高精度分子间势能面通过非线性拟合振动平均的基态(v_3=0)和激发态(v_3=1)的能量而得到.其中对2038个格点拟合的方差仅为0.08 cm-1.我们发现,该体系相互作用能的主要贡献来源于短程的交换排斥能和长程的色散能,势能面有3个极小值点,在θ方向表现出强烈的各向异性,而在φ方向各向异性比较弱.与CH_3F-He体系相比,分子间相互作用较强,分子间相互作用对单个分子的束缚较明显,CH_3F和Ar表现得更像紧束缚的超分子.随后,采用径向离散变量表象(DVR)和角度有限基组表象(FBR)相结合的杂化基函数方法进行展开,用Lanzcos迭代的方法计算了该体系的振转能级,并首次报道了该体系的微波、红外光谱.值得指出的是,由于相互作用比较强,形成的势阱较深,该体系在两个不同的势阱中分别形成了分子间振动基态和激发态,这与CH_3F-He体系中全空间离域的能级态形成鲜明的对比.     

确定三原子分子势能面的SCF-CI方法

建议根据振动高激发态的实验能级确定三原子分子势能面的SCF-CI方法.此法中使用键长-键角内坐标系下的SCF-CI方法来精确地计算振动高激发态的能级及其对势能参数的一阶微分,并使用LMF算法来优化势能参数.为验证此方法,优化了水分子的势能面,计算出水分子的振动高激发态能级与70个观测到的振动能级相比较,标准偏差为1.15cm~(-1).     

复合物C6H5CH3…Ar分子间的外部振动频率

利用双光子共振电离光谱和飞行时间质谱技术在超声分子束中观察到C6H5CH3…Ar的振动光谱.借助同位素光谱效应、内转动能级和分子间振动能级的理论计算,合理地归属了涉及CH3转动和Ar原子振动的光谱,并由此获得复合物分子间各种模式的振动频率.     

三原子分子振转激发态的理论研究

三原子分子振转激发态的理论研究谢代前,鄢国森,田安民（四川大学化学系,成都,６１００６４）关键词振转激发态,三原子分子,变分法能级较高的振转激发态通常包含大振幅运动,其波函数分布于很广的势能面区域内,传统的正则模理论已不适合于解决这类问题．近年来,Ｈ...     

双原子分子振转激发态

发展了一种较简便的方法来计算双原子分子的振转激发态。原则上,该方法可以给出同数值积分方法一样精确的结果。作为例子, 给出了H_2的162个振转激发态的准确能级和波函数。     

甲基氟苯-氩复合物的分子间振动

在超声分子束装置上,利用多光子共振光电离和飞行时间质谱技术,观测到了三种甲基氟苯-氩复合物S1←S0电子态跃迁的两光子共振电离光谱。通过对光谱的分析和理论计算,获得了三种复合物内部分子音质伸缩振动和弯曲振动的频率。结果显示,F原子和CH3在苯环上取代H原子后,它们的相对位置对伸缩振影响很小,而与弯曲振动有较大的相关。另外,在复合物内部,Ar原子对甲基氟苯中的CH3内转动有明显的阻碍作用。     

An SCF-CI method for determining the potential energy surface of a triatomic molecule

A self-consistent-field (SCF)-configuration interaction (CI) (SCF-CI) method for determining the potential energy surface of a triatomic molecule from the observed vibrational band origins has been suggested. By this method, the SCF-CI procedure in the internal coordinates is used to calculate the vibrational bond origins and their first derivatives with respect to parameters in the potential energy function using the exact vibrational Hamiltonian, and the optimizer LMF in the nonlinear-squares problem is employed to optimize parameters in the potential energy function. This approach is used to optimize the potential energy function of the water molecule. The standard deviation of this fitting to the 70 observed band origins is 1.154cm-1.     

A five-dimensional potential energy surface and predicted infrared spectra for the N2O-hydrogen complexes

We present a five-dimensional potential energy surface for the N(2)O-hydrogen complex using supermolecular approach with the full counterpoise correction at the coupled-cluster singles and doubles with noniterative inclusion of connected triple level. The normal mode Q(3) for the nu(3) antisymmetric stretching vibration of the N(2)O molecule was included in the calculations of the potential energies. The radial discrete variable representation/angular finite basis representation method and Lanczos algorithm were employed to calculate the rovibrational energy levels for four species of N(2)O-hydrogen complexes (N(2)O-para-H(2), -ortho-H(2), -ortho-D(2), and -para-D(2)) without separating the inter- and intramolecular vibrations. The calculated band origins are all blueshifted relative to the isolated N(2)O molecule and in good agreement with the experimental values. The calculated rotational spectroscopic constants and molecular structures agree well with the available experimental results. The frequencies and line intensities of the rovibrational transitions in the nu(3) region of N(2)O for the van der Waals ground vibrational state were calculated and compared with the observed spectra. The predicted infrared spectra are consistent with the observed spectra and show that the N(2)O-H(2) complexes are mostly a-type transitions while both a-type and b-type transitions are significant for the N(2)O-D(2) complexes.     

Fully quantum rovibrational calculation of the He(H2) bound and resonance states

In a recent paper [J. Chem. Phys. 2005, 122, 124318], a full-dimensional quantum method, designed to efficiently compute the rovibrational states of triatomic systems with long-range interactions, was applied to the benchmark Li-(H2) ion-molecule system. The method incorporates several key features in order to accurately represent the rovibrational Hamiltonian using only modestly sized basis sets: (1) exact analytical treatment of Coriolis coupling; (2) a single bend-angle basis for all rotational states; (3) phase space optimization of the vibrational basis; (4) G(4) symmetry adaptation of the rovibrational basis. In this paper, the same methodology is applied for the first time to a van der Waals complex system, He(H2). As in the Li-(H2) study, all of the rovibrational bound states, and a number of resonance states, are computed to very high accuracy (1/10,000 of a wavenumber or better). Three different isotopologues are considered, all of which are found to have a single bound state with a very low binding energy. Several extremely long-lived Feshbach resonances are also reported.     

Free Rotor Model or Rigid Rotor Model for CH3F-Ne Complex and Comparison with Other CH3F-Rare Gas Systems

The assignment of the rovibrational spectra of molecule-Ne complexes is always a challenge to study van der Waals systems, since they usually exhibit behavior intermediate between free rotor and rigid rotor. In this paper, the microwave and infrared spectra of CH₃F-Ne, a model system for symmetric-top-atom dimer, were firstly predicted and analyzed based on the four-dimensional ab initio intermolecular potential energy surfaces(PESs), which explicitly incorporate the v₃(C-F) stretch normal model coordinate of the CH₃F monomer. Analytic three-dimensional PESs were obtained by least-squares fitting vibrationally averaged interaction energies for v₃(CH₃F)=0 and 1 to the Morse/long-range(MLR) potential function for symmetry top impurity with atom model. These PESs fitting to 2340 points have root-mean-square(RMS) deviations of 0.07 cm^-1, and require only 167 parameters. Based on the analytical vibrationally averaged PESs, the rovibrational energy levels were calculated by employing Lanczos algorithm, with combined radial discrete variable representation and parity-adapted angular finite basis representation. Based on the wavefunction analysis and comparison of CH₃F-Ne with CH₃F-He and CH₃F-Ar complexes, the bound states were assigned. Spectral parameters for CH₃F-Rg(Rg:rare gas, Rg=He, Ne, Ar) complexes were fitted and discussed. Temperature dependent transition intensities for CH₃F-Ne were also reported and analyzed. The complete microwave and infrared spectra information for CH₃F-Ne made it possible to provide important guidance for future experimental spectroscopic assignments.     

Symmetry assignment in the distributed Gaussian functions method to study homonuclear rotating trimers

An approximate method based on the use of distributed Gaussian functions (DGF) to describe the interparticle distances is employed to study the rovibrational spectrum of trimers. Rotational energy levels are obtained by assuming that vibration and rotation are separated. Thus, eigenstates of the Hamiltonian for the zero total angular momentum, J = 0, are used as basis set to solve the rotational Hamiltonian. A procedure to identify the corresponding symmetry character for the rovibrational bound states is proposed. The DGF approach is applied to the case of the rotating Ar₃ trimer. The reliability of the method is tested by comparison with results from an exact hyperspherical coordinate calculation for J = 0, 1 and 6.     

Beryllium dimer,a critical test case of MBPT and CI methods

The ground-state potential curve for the beryllium dimer is calculated as a critical test case for methods based on many-body perturbation theory (MBPT ) and configuration interaction (CI ). In particular, the recently proposed double excitation (DE ) MBPT method is compared to the standard SCF-CI method including single and double excitations from a single reference determinant. The SCF-CI method is shown to give surprisingly accurate results compared to more complete CI calculations including a larger configuration space, whereas the DE-MBPT method breaks down more or less completely, particularly for larger basis sets. The results thus demonstrate the importance of including the renormalization terms in this case. Finally, Davidson's correction and related methods lead to an even more severe breakdown than the DE-MBPT method.     

Möbius and hückel systems in the SCF and CI approximations

Over the last 15 years since the introduction of Möbius-Hückel theory, a number of varied questions has accumulated. The most interesting of these deals with the question of whether or not the Möbius-Hückel theory is valid in the SCF and SCF-CI approximation. This paper presents a treatmet which shows that the repulsion and exchange contributions are independent of the Möbius vs Hückel nature of the orbital array. Also it is shown that the one-electron terms are symmetry determined and derive from SCF coefficients. An analytical SCF-CI treatment if given. Several further unanswered questions are also considered.     

Convergence of the density functional one-centre expansion for the molecular continuum: N2 and (CH3)3N

A large-scale one-centre expansion with a radial B-spline basis set is implemented for bound and continuum states. A Kohn-Sham hamiltonian is employed with Hartree and exchange-correlation potentials calculated from the SCF electron density taken from a previous LCAO calculation. An inverse iteration method is used to obtain the continuum wavefunction, from which the cross section and asymmetry parameter are calculated. The convergence with respect to angular momentum and cut-off radius is analysed for N₂. The relevance of multipolar contributions even at large distances is shown and suggestions for further improvements are given. In order to show that the present method is suitable to treat systems of moderate size, the (CH₃)₃N molecule has also been calculated and the results are compared with experiment. Received: 19 May 1998 / Accepted: 20 August 1998 / Published online: 7 December 1998     

Accurate ab initio pair potentials between helium and the heavier group 2 elements

Interactions between the heavier Group 2 metals (Ca, Sr, and Ba) and helium were studied using the well-tempered model core potential method. Accurate pair potentials, calculated at the coupled-cluster level of theory with very large basis sets, were used in bound state calculations. Three bound rovibrational states were found for each complex. The pair-potential parameters were used to predict how each of the metal atoms would be solvated by a helium nanodroplet. The Ca atom is not fully solvated by the droplet and the interaction between the helium and the metal decreases from Ca to Ba. This agrees with the experimental observation that the spectra of these atoms in a nanodroplet are intermediate between the spectra of the free atoms and the spectra in liquid helium.     

A new ab initio intermolecular potential energy surface and predicted rotational spectra of the Ar-H2S complex

We report a reliable three-dimensional ab initio intermolecular potential energy surface for the Ar-H(2)S complex with H(2)S monomer fixed at its experimental average structure. The potential energies were evaluated using the supermolecular approach at the coupled-cluster level with a large basis set including bond functions. The full counterpoise procedure was used to correct the basis set superposition error. The potential has a planar T-shaped global minimum with a well depth of 177.48 cm(-1) at the intermolecular distance of 3.72 ?. An additional planar local minimum is also found and is separated from the global minimum with an energy barrier with a height of 47.46 cm(-1). The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to calculate the rovibrational energy levels for three isotopic species of Ar-H(2)S complexes (Ar-H(2)(32)S, Ar-H(2)(33)S, and Ar-H(2)(34)S). The rotational transition frequencies and structural parameters for the three isotopomers were also determined for the ground and the first excited states, which are all in good agreement with the available experimental values.