The analytical potential energy function of flue gas SO2（X1A1）

The equilibrium structure of flue gas SO2 is optimized using the density functional theory （DFT）/B3P86 method and CC-PV5Z basis. The result shows that it has a bent （C2v, X1A1） ground state structure with an angle of 119.1184°. The vibronic frequencies and the force constants are also calculated. Based on the principles of atomic and molecular reaction statics （AMIIS）, the possible electronic states and reasonable dissociation limits for the ground state of SO2 molecule are determined. The potential functions of SO and 02 are fitted by the modified Murrell-Sorbie＋c6 （M-S＋c6） potential function and the fitted parameters, the force constants and the spectroscopic constants are obtained, which are all close to the experimental values. The analytic potential energy function of the SO2 （X1A1） molecule is derived using the many-body expansion theory. The contour liues are constructed, which show the static properties of SO2 （XIA1）, such as the equilibrium structure, the lowest energies, the most possible reaction channel, etc.     

Temperature-dependent ionization potential of sodium clusters

The ionization potential of sodium clusters () at a finite temperature is studied using density functional theory and ab initio molecular dynamics. The threshold regions of the photoionization efficiency curves are deduced from the integrated IP distributions, which are obtained from the energy eigenvalues of the highest occupied Kohn-Sham states during molecular dynamics by applying a theoretically well-defined shift. The calculated ionization potentials are directly compared to the experimental values. The energetically best geometry of Na₅₅ is found to be a slightly distorted icosahedron. Received 16 April 1999 and Received in final form 6 July 1999     

A computational study of C?X (X = H,C, F,Cl) bond dissociation enthalpies (BDEs) in polyhalogenated methanes and ethanes

Bond dissociation enthalpies (BDEs) play a significant role in the photolysis of Chlorofluorocarbons (CFCs) and hydrochlorofluorocarbons (HCFCs), which lead to the depletion of stratospheric ozone. In this work, we estimate the performance of Density Functional Theory (DFT) methods in calculating BDEs of CFCs and HCFCs, and find that DFTs are unreliable for this system. The reasons for the unreliability of DFT methods in this system are also concluded. Furthermore, composite ab initio methods G3 and G3B3 are demonstrated to accurately estimate BDEs of polyhalogenated lower alkanes. Eighty two experimental values from Comprehensive Handbook of Chemical Bond Energies (2007, 2nd edition) are re‐evaluated. Eight of them are doubted as having a deviation exceeding 20.0 kJ/mol between the theoretical and experimental values. We also systematically predict the BDEs in polyhalogenated methanes and ethanes. A further study is conducted on their relationships of structures and properties. Copyright © 2010 John Wiley & Sons, Ltd.     

XY(H,Li, Na)分子基态的结构与势能函数

运用群论及原子分子反应静力学方法，推导了XY(H，Li，Na)分子基态的电子态及相应的离解极限.并采用密度泛函方法(B3LYP)和二次组态相互作用方法(QCISD)优化计算了XY(H，Li，Na)分子基态的平衡结构、振动频率和离解能.使用QCISD/6-311++G(3df,3pd)方法，对XY(H，Li，Na)分子基态进行了单点能扫描计算，采用最小二乘法拟合Murrell-Sorbie函数得到了相应的势能函数和与该基态相对应的光谱常数(B_e,α_e,ω_e和ω_eχ_e)，计算结果与实验数据符合得相当好. 关键词： 组态相关方法 基态 势能函数     

7Li2分子23Πu激发态的解析势能函数、振动能级及其转动惯量

使用Gaussian 03程序包中的“对称性匹配簇-组态相互作用”方法、在0.13—2.0nm的核间距范围内利用6-311+〖KG-*3〗+G(d,p)基组对⁷Li₂(2³Π_u)分子的势能曲线进行了计算, 同时使用最小二乘法将计算结果拟合成了解析势能函数. 利用拟合出的解析势能函数并结合Rydberg-Klein-Rees方法, 计算了该态的谐振频率, 进而计算了该态的其他光谱常数, 分别为T关键词： 解析势能函数 谐振频率 振动能级 转动惯量     

Bound states and critical behavior of the Yukawa potential

We investigate the bound states of the Yukawa potential V (r)=−λexp(−αr)/r, using different algorithms: solving the Schr?dinger equation numerically and our Monte Carlo Hamiltonian approach. There is a critical α = α_C, above which no bound state exists. We study the relation between α_C and λ for various angular momentum quantum number l, and find in atomic units, α_C(l) = λ[A ₁ exp(−l/B ₁) + A ₂ exp(−l/B ₂)], with A ₁ = 1.020(18), B ₁ = 0.443(14), A ₂ = 0.170(17), and B ₂ = 2.490(180).     

He原子与N2分子相互作用势的理论研究

用量子力学从头算方法深入研究了He原子与N2分子的相互作用势,用不同的方法和基组计算了该体系的各项异性相互作用势的单点能数据,得出了空间势能面的分布情况,并选定CCSD(T)/6-311++G(3df,2pd)方法和基组,同时采用了Boys和Bernardi提出的Full Couterpoise方法,消除了计算中的基组重叠误差,得出了较为准确的He-N2体系相互作用势的解析表达式.通过计算得到的微分截面与实验值符合较好,同时得出了不同碰撞能量时He原子与N2分子碰撞的微分截面的规律.     

Accurate potential energy function and spectroscopic study of the X^2∑＋, A^2П and B^2∑＋ states of the CP radical

This paper calculates the equilibrium internuclear separations, the harmonic frequencies and the potential energy curves of the X^2∑＋, A^2П and B^2∑＋ states of the CP radical by the highly accurate valence internally contracted multireference configuration interaction method with correlation-consistent basis sets （aug-cc-pV6Z for C atom and aug-cc-pVQZ for P atom）. The potential energy curves are all fitted with the analytic potential energy function by the least-square fitting. Employing the analytic potential energy function, we determine the spectroscopic constants （Be, αe and ωeχe） of these states. For the X2∑＋ state, the obtained values of De, Be, αe, ωeχe, Re and ωe are 5.4831 eV, 0.792119 cm-1, 0.005521 cm-1, 6.89653 cm-1, 0.15683 nm, 12535.11 cm-1, respectively. For the A2H state, the present values of De, Be,αe, ωeχe, Re and We are 4.586 eV, 0.703333 cm-1, 0.005458 cm-1, 6.03398 cm-1, 0.16613 nm, 1057.89 cm-1, respectively. For the B2E＋ state, the present values of De, Be, αe, ωeχe, Re and We are 3.506 eV, 0.677561 cm-1, 0.00603298 cm-1, 5.68809 cm-1, 0.1696 nm, 822.554 cm-1, respectively. For these states, the vibrational states with the rotational quantum number J equals zero （J = 0） are studied by solving the radial nuclear Schr6dinger equation using the Numerov method. For each vibrational state, the vibrational level, the classical turning points, the rotational inertial and the centrifugal distortion constants are calculated. Comparison is made with recent theoretical and experimental results.     

密度泛函方法对SiO分子基态（X 1Σ+）势能函数的研究

利用密度泛函B3P86方法,分别选用STO-3G,D95^**,6-311G,6-311++G,6-311++G^**,cc-PVTZ基组对SiO分子基态(X¹Σ⁺)进行结构优化计算.通过比较得出,cc-PVTZ基组为对SiO分子基态(X¹Σ⁺)进行结构优化最优基组的结论.使用密度泛函B3P86方法,选用cc-PVTZ基组进 关键词： B3P86 SiO 势能函数 光谱常数     

High Pressure effect on fluorescence lifetime τ for magnetic dipole 5D0→5F1 transitions in YAG:Eu3+

The fluorescence lifetime for magnetic dipole ⁵D₀→⁷F₁ transition in yttrium aluminum garnet doped with Eu³⁺ (YAG:Eu³⁺) crystal was studied under the pressure of up to 10.4?GPa at room temperature. The fluorescence lifetime τ (⁵D₀→⁷F₁ transition) slowly decreased with pressure. The pressure effect on τ (⁵D₀→⁷F₁ transition) was explained with a model which considered pressure effect on line position: inter-ionic distance, ion volume, molecular volume, ion polarizability, molecular polarizability, sample refractive index, and surrounding hydrostatic medium refractive index. The fluorescence lifetime τ calculated by the presented model was in close correspondence with the experimental values.     

On the application of the theory of the brownian movement in a periodic potential to the study of inertial effects and dipole coupling

For a series of ten electron molecules (HF, H₂O, NH₃, CH₄) the molecular polarizability tensor and the derivatives with respect to the symmetry coordinates have been calculated from ab initio SCF wavefunctions using the finite field method as well as perturbation theory approaches. Raman intensities and degrees of depolarization derived from the finite field results agree well with the available experimental data.

The zeroth order bond polarizability model and the atom dipole interaction model have been analysed. Both models can be used to describe the computed static polarizabilities and the derivatives with respect to bond stretching, but fail for the derivatives with respect to the bending coordinates.     

利用WBEPM模型和稳定变分法计算氦原子的极化率和色散系数

以最弱受约束电子势模型(WBEPM)理论为基础,建立了计算氦原子多极动态极化率和相互作用色散系数的稳定变分方法 .导出了该方法中所涉及的矩阵元和线性方程组的解析表达式.作为应用,具体计算了基态氦原子的极化率和两体色散系数,将计算结果与用其他方法所得到的结果进行了比较,数据基本一致.     

利用WBEPM模型和稳定变分法计算氦原子的极化率和色散系数

以最弱受约束电子势模型(WBEPM)理论为基础,建立了计算氦原子多极动态极化率和相互作用色散系数的稳定变分方法.导出了该方法中所涉及的矩阵元和线性方程组的解析表达式.作为应用,具体计算了基态氦原子的极化率和两体色散系数,将计算结果与用其他方法所得到的结果进行了比较, 数据基本一致.     

B2C(1A1)和BC2(2A′)的结构与解析势能函数

采用单双取代的二次组态相互作用方法，分别选用6-311++G(d,p)和6-311G(df,pd)基组，对B₂C和BC₂分子的结构进行了优化，得到这两个分子的基态结构为C_2v和C_s，基态电子状态为¹A₁和²A′，同时还得到了它们的平衡几何结构、离解能、谐振频率和力常数. 关键词： 碳化硼 Murrell-Sorbie函数 谐振频率 势能函数     

The influence of intermolecular interactions on the Kerr effect in gases

Expressions are derived for the second and third Kerr virial coefficients, B _K and C _K, of spherical top molecules in terms of irreducible cluster polarizabilities, and values are calculated using the dipole-induced dipole model for argon, krypton, xenon, methane, tetrafluoromethane, neopentane and sulphur hexafluoride. For mixtures of rare gases it is shown that the collision-induced dipole moment makes a negligible contribution to B _K. The effect of the choice of intermolecular potential function on the calculated second Kerr virial coefficients is also demonstrated. It is found that the predominant contributions to C _K arise from the pair polarizability, and that the triplet polarizability is only of minor importance.     

CALCULATION OF THE FINE STRUCTURE OF OXYGEN-LIKE IONS USING THE POLARIZATION POTENTIAL FUNCTION

We have calculated the forbidden transition energies and magnetic dipole transition probabilities of 2s²2p⁴(³P₁－³P₂) and 2s²2p⁴ (³P₀－³P₁) of oxygen-like isoelectronic sequences (Z=10－32) by a method of polarization potential correction. The transition energies show good agreement with experiment and are much better than the calculations in the literature. These results also illustrate that it is feasible to use the dipole expansion of the polarization potential to deal with some dynamic and non-dynamic effects in the central field approach. The relation of polarizability and cut-off radius with atomic number is discussed. We also give the fitted formula between the polarizability α₁ and atomic number Z as α₁=0.73429-9.56644×10^-4Z+7.43016×10^-5Z^{2-2.53298×10-6Z3+2.08306×10-8Z4.}     

Calculations of ozone line shifting induced by N2 and O2 pressure

Ozone line shifts by nitrogen and oxygen pressure are computed for the ν₁ + ν₃, 2ν₁ and 2ν₃ bands of the 5 μm spectral region by a semiempirical approach. The calculated values agree with measurements better than 0.001 cm^?1 atm^?1 for 98% of O₃–N₂ lines and 87% of O₃–O₂ lines. In contrast with the water molecule case, the polarization components of the interaction potential are shown to contribute to the line shift more efficiently than the electrostatic interactions. As intermediate results, the mean dipole polarizability and the components of the polarizability tensor for the vibrational states (101), (200), and (002) of ozone molecule are determined by least-squares fitting of theoretical shifts to some experimental values. The temperature exponents for the ν₁ + ν₃ band lines are also estimated.     

OH, OCI, HOCI(1A')的从头算与势能曲线

采用Gassian09程序包中的多种方法对OH, OCI, HOCI分子的基态结构进行优化计算, 优选出QCISD/6-311G(2df), B3P86/6-311+G(2df)方法分别对OH(X²), OCI(X²)分子进行计算, 得到平衡核间距R_OH=0.09696 nm, R_OCI=0.1569 nm, 谐振频率ω(OH)=3745.37 cm^-1, ω(OCI)=892.046 cm^-1, 与实验结果非常符合. 用Murrell-Sorbie势能函数对OH和OCI分子的扫描势能点进行拟合, 其扫描点都与四参数Murrell-Sorbie函数拟合曲线符合得很好．优选出QCISD(T)/D95(df, pd)方法对HOCI分子进行计算, 得到基态为X¹A', 键长R_OH =0.0966 nm, 键角∠HOCI=102.3°, 谐振频率ω₁(a₁)=738.69 cm^-1, ω₂(b₂)=1260.25 cm^-1, 离解能D_e=2.24eV. 通过比较发现这些结果与实验值符合得很好,并优于文献报道的结果. 随后计算出了力常数, 在此基础上,推导出HOCI分子的多体展式势能函数.报道了HOCI分子对称伸缩振动势能图中在H+OCI →HOCI反应通道上有一鞍点, H原子需要越过1.74eV的能垒才能生成HOCI的稳定结构, 在Cl+OH→HOCI通道上不存在明显势垒, 容易形成稳定的HOCI分子.     

A quasi-classical trajectory study of the Cl + H<Subscript>2</Subscript> (D<Subscript>2</Subscript>) reaction on a new BW3 potential energy surface

Molecular reaction dynamics of Cl + H₂ (D₂) has been studied on the latest analytical potential energy surface called BW3 using the Monte Carlo quasi-classical trajectory method. Excitation functions, differential cross sections and angular distributions of HCl and DCl products have been calculated. The excitation functions of the Cl (²P_3/2) + n-H₂ and Cl(²P_3/2) + n-D₂ reactions are also studied. The results are compared with those of quasi-classical trajectory [M. Alagia et al.: Phys. Chem. Chem. Phys. 2 (2000); F. J. Aoiz et al.: J. Phys. Chem. 100 (1996)], quantum mechanical (QM) calculations [F. J. Aoiz et al.:J. Chem. Phys. 115 (2001)] and experimental data [S. H. Lee et al.: J. Chem. Phys. 110 (1999); F. Dong et al.: J. Chem. Phys. 115 (2001)]. Discussions are given to some new results.