MBPT studies of van der Waals molecules

The many-body perturbation theory is employed for the calculation of the interaction potential for the F^- … He system in the framework of the supermolecule method. A particular attention is paid to the choice of the basis set functions for the two subsystems and the related basis set superposition effects. It has been found that the main features of the interaction potential are recovered in the SCF approximation. The SCF potential has its minimum at the distance R = 6·4 a ₀ with the interaction energy of 53 cm^-1. The complete fourth-order MBPT method gives the potential minimum position and depth equal to 6·5 a ₀ and 64 cm^-1, respectively. The basis set superposition effects estimated by using the counterpoise technique are negligibly small for the SCF interaction potential, while at the correlated level their magnitude is comparable to the value of the total correlation contribution to the interaction energy. The basis set superposition effect in calculations of the electron correlation contribution to weak intermolecular interactions is found to be the major factor limiting the reliability of the corresponding theoretical data.     

MCSCF study of the spin-orbit coupling in OH (X 2Π i )

Full-valence MCSCF calculations were performed for the X ²Π _i state of OH in that it was hoped that these would predict the spin-orbit coupling A(r) to within a few tenths of a wavenumber over the entire range of r. The functional form was correct but the actual values did not meet our expectations. We briefly considered higher-order relativistic corrections.

A universal even-tempered basis set was tested relative to an ‘extended’, optimized basis. The universal basis set did remarkably well and highlighted significant r-dependent and, presumably, state dependent problems with the smaller basis.     

Basis set limit binding energies of hydrogen bonded clusters

The utility of the recently developed extrapolation method to estimate the binding energies of weakly bound clusters at the basis set limit exploiting the similar basis set convergence behaviour of correlation energies of the monomer and cluster in correlated calculations (J. chem. Phys., 116, 5389 (2002)) was tested for small to medium (HF)_n and (H₂O)_n (n = 2–5) clusters using various correlation consistent cc-pVXZ (X= D,T,Q,5) basis sets containing different numbers of diffuse functions and 6–31G type basis sets at the MP2 and CCSD(T) level for which accurate basis set limits are available for comparison. It is shown that the basis set limit binding energies estimated by this extrapolation method with modest size of basis sets (cc-pVDZ/cc-pVTZ or 6–3 1 G(d,p)/6-3 IG(2df,2pd)) are much closer to the exact basis set limits than the estimates by commonly used X ^?3 extrapolation or counterpoise corrected binding energies, signifying the importance of this extrapolation method for the study of large weakly bound clusters. It is also shown that the inclusion of appropriate diffuse functions in the basis sets can significantly improve the accuracy of the estimated basis set limits by this extrapolation method. For (HF)_n clusters the MP2 and CCSD(T) basis set limits estimated by this extrapolation method with aug-cc-pVDZ and aug-cc-pVTZ basis sets are 18.4 (18.5) and 18.9 (18.9) for the dimer, 61.8 (62) and 63.2 (63) for the trimer, 113.5 and 114.7 (116) for the tetramer, and 155.2 and 156.3 (158) for the pentamer, respectively, with the values in parentheses representing the apparent basis set limits, with the numbers in units of kJ mol^?1. The corresponding results for (H₂O)_n clusters are 20.5 (20.5) and 20.6 (20.7) for (H₂O)₂, and 60.5 (61) and 60.1 (60) for (H₂O)₃, respectively.     

A new basis set for molecular wavefunctions

A new basis set is proposed for molecular self-consistent field and configuration interaction calculations. Expansion functions are proposed in the form sin (r _A · p) exp (-ar_A ²) and cos (r _A · p) exp (-ar _A ²) and it is suggested that they are used as an alternative to the gaussian basis set x_A ⁱy_A ^jz_A ^k exp (-ar _A ²).

It is shown how the new basis can have the same effects, for suitable vectors p, as the gaussian basis. It is further shown that all the one and two-electron integrals in the new basis are evaluable in terms of exponentials, square roots and the complex error function erf (z). First calculations, using the new basis, are presented in LiH, HF and H₂O.     

A comparison of finite basis set and finite difference Hartree—Fock calculations for the open-shell (X 2Σ+) molecules BaF and YbF

A comparison is made of the accuracy with which the total electronic energy can be calculated by using the finite basis set approach (the algebraic approximation) and the finite difference method in calculations employing the Hartree—Fock model for the open shell ground (X ²Σ⁺) states of the fluorides BaF and YbF. The convergence of the calculations carried out within the algebraic approximation is monitored by employing systematically constructed basis sets of increasing size. The difference between the finite basis set and finite difference Hartree—Fock energies is 2.6μE _h for BaF and 2.8μE _h for YbF. Dipole moments determined within the algebraic approximation are also compared with the corresponding finite difference expectation values.     

TiO基态 (X 3 Δr) 的势能函数与光谱常数

应用群论及原子分子反应静力学方法推导了TiO分子基态(X³Δ_r)的离解极限.采用不同的计算方法,包括密度泛函B3LYP,B3P86,BP86,B3PW91和MP2,MP4方法,结合不同基组计算了TiO分子基态的平衡核间距、能量和振动频率.研究表明,使用B3LYP方法,对O原子使用6-311+G基组,Ti原子使用6-311+ +G^**基组时计算得到的平衡几何结构、分子离解能和谐振频率与实验值符合得最好.使用优选出的方法和基组对T 关键词： TiO 势能函数 光谱常数 密度泛函理论     

Accurate ab initio anharmonic force field and heat of formation for silane

From large basis set coupled cluster calculations and a minor empirical adjustment, an anharmonic force field for silane has been derived that is consistently of spectroscopic quality (±1 cm^?1 on vibrational fundamentals) for all isotopomers of silane studied. Inner-shell polarization functions have an appreciable effect on computed properties and even on anharmonic corrections. From large basis set coupled cluster calculations and extrapolations to the infinite-basis set limit, we obtain TAE₀ = 303.80 ± 0.18 kcal mol^?1, which includes an anharmonic zero-point energy (19.59 kcal mol^?1), inner-shell correlation (—0.36 kcal mol^?1), scalar relativistic corrections (— 0.70 kcal mol^?1) and atomic spin-orbit corrections (—0.43 kcal mol^?1). In combination with the recently revised ΔH ^o _{f, o}[Si(g)], we obtain ΔH ^o _f.o[SiH₄(g)] = 9.9 ± 0.4 kcal mol^?1 in between the two established experimental values.     

Norm Kernels and the Closeness Relation for Pauli-Allowed Basis Functions

The norm kernel of the generator-coordinate method is shown to be a symmetric kernel of an integral equation with its eigenfunctions defined in the Fock-Bargmann space and forming a complete set of orthonormalized states (classified with the use of the SU(3) symmetry indices) satisfying the Pauli exclusion principle. This interpretation has allowed us to develop a method which, even in the presence of the SU(3) degeneracy, provides for a consistent way to introduce additional quantum numbers for the classification of the basis states. In order to set the asymptotic boundary conditions for the expansion coefficients of a wave function in the SU(3) basis, a complementary basis of functions with partial angular momenta as good quantum numbers is needed. Norm kernels of the binary systems ⁶He+p, ⁶He+n, ⁶He+⁴He, and ⁸He+⁴He are considered in detail.Received December 16, 2002; accepted February 27, 2003 Published online August 25, 2003     

The infrared spectrum of isoxazole in the range 600–1400 cm−l,including a high-resolution study of the v 7(A′) band at 1370.9 cm−1 and the v 16(A″) band at 764.9 cm−l,together with ab initio studies of the full spectrum

Abstract

The high resolution infrared gas phase spectrum of isoxazole in the range 600–1400 cm^?1 has been recorded and more precise centres obtained for a number of bands; analyses of the v ₇(A′) band at 1370.9cm^?1 and the v ₁₆(A″) band at 764.9cm^?1 have been performed. Use of the Watson Hamiltonian, A-reduction, III^r-representation and simultaneous analysis of the present IR and previous microwave data, has led to rotation constants and quartic centrifugal distortion constants Δ _J , Δ _JK , ΔK, δ _J and δ _K for the ground state and for the v ₁₆ vibrationally excited states. The equilibrium structures and the derived harmonic frequencies have been calculated by ab initio methods using triple zeta+polarization (TZVP) and cc-pVTZ basis sets, with MP2, MP4 and CCSD(T) methods. At each methodology, for closest numerical agreement between the calculated and observed rotation constants, the optimum basis set seems to be TZVP rather than cc-pVTZ basis sets. However, the order of the highest A″ and lowest A′ symmetry vibrations is only resolved by the cc-pVTZ basis set with the MP4 methodology, which does generate the experimental sequence (v ₁₄>v ₁₃) The CCSD(T) methodology does not lead to significant difference over either MP2 or MP4 with the TZVP basis set.     

Accurately solving the electronic Schrödinger equation of atoms and molecules using explicitly correlated (r 12-)MR-CI IV. The helium dimer (He 2)

The helium dimer interaction potential is computed using the recently proposed (explicitly correlated) r ₁₂-MR-ACPF (averaged coupled-pair functional) method and a [11s8p6d5f4g] basis set. With an MR-ACPF ansatz that contains 121 references we obtain interaction energies that are close to full CI. In a smaller reference space containing 9 functions, however, even by successively adding [3h] and [2i] functions to the basis set mentioned above, the basis set limit could not be reached. While convergence to the basis set limit is slow, it nevertheless is monotonic and therefore allows for extrapolation to the limit. We obtain basis set corrections at R = 4 a ₀ and 5.6 a ₀ which we further extend to all distances and which we apply to the potential energy curve mentioned above. From our calculations, we conclude that a very recent potential which has been calculated using the SAPT (symmetry adapted perturbation theory) method, and which previously was assumed to be the most accurate available, is insufficiently repulsive at short distances. We correct our extrapolated potential for retardation and finally calculate the expectation value of the interatomic distance (?R?) and dissociation energy (D ₀) by solving the Schrödinger equation of the vibrating ⁴He₂. Our results (?R? = 41 ± 13 Å and D ₀ = 2.2 ± 1.0 mK) are in acceptable agreement with very recent calculations in the literature, but they disagree with a recent experiment.     

Interaction potentials,spectroscopy and transport properties of C+(2PJ) and C+(4PJ) with helium

We calculate accurate interatomic potentials for the interaction of a singly charged carbon cation with a helium atom. We employ the RCCSD(T) method, and basis sets of quadruple-ζ and quintuple-ζ quality; each point is counterpoise corrected and extrapolated to the basis set limit. We consider the two lowest C⁺(²P) and C⁺(⁴P) electronic states of the carbon cation, and calculate the interatomic potentials for the terms that arise from these: ²Π and ²Σ⁺, and ⁴Π and ⁴Σ^?, respectively. We additionally calculate the interatomic potentials for the respective spin–orbit levels, and examine the effect on the spectroscopic parameters. Finally, we employ each set of potentials to calculate transport coefficients, and compare these to the available data. Critical comments are made in the cases where there are discrepancies between the calculated values and measured data.     

boldmath B3Σu-–X3Σg- transition in selenium dimer: ab initio multireference configuration interaction calculations

Theoretical investigation of low-lying electronic states and B³Σ_u^-–X³Σ_g^- transition properties of selenium dimer using size-extensivity singly and doubly excitation multireference configuration interaction theory with nonrelativistic all-electron basis set and relativistic effective core potential plus its split valence basis set is presented in this paper. The spectroscopic constants of ten low-lying Λ-S bound states have been obtained and compared with experiments. Spin-orbit calculations for coupling between B³Σ_u^- sates and repulsive ¹Π_u, ⁵Π_u states have been made to interpret the predissociation mechanisms of the B³Σ_u^- state. The lifetimes of B³Σ_u^- (v=0～6) have been calculated with scalar relativistic effects included or excluded, respectively, and reasonably agree with experimental values.     

Theoretical study of Si+(2PJ)–RG complexes and transport of Si+(2PJ) in RG (RG = He–Ar)

We calculate accurate interatomic potentials for the interaction of a singly charged silicon cation with a rare gas atom of helium, neon or argon. We employ the RCCSD(T) method, and basis sets of quadruple-ζ and quintuple-ζ quality; each point is counterpoise-corrected and extrapolated to the basis set limit. We consider the lowest electronic state of the silicon atomic cation, Si⁺(²P), and calculate the interatomic potentials for the terms that arise from this: ²Π and ²Σ⁺. We additionally calculate the interatomic potentials for the respective spin-orbit levels, and examine the effect on the spectroscopic parameters; we also derive effective ionic radii for C⁺ and Si⁺. Finally, we employ each set of potentials to calculate transport coefficients, and compare these to available data for Si⁺ in He.     

A comparison of finite basis set and finite difference Hartree-Fock calculations for the open- shell (X2Σ+) molecules BeF,MgF, CaF and SrF

A comparison is made of the accuracy with which the total electronic energy can be calculated by using either the finite basis set approach (the algebraic approximation) or finite difference methods in calculations employing the restricted Hartree-Fock model for the open-shell ground (X²Σ+) states of the Group Ha fluorides BeF, MgF, CaF and SrF. The convergence of the calculations carried out within the algebraic approximation is monitored by employing systematically constructed basis sets of increasing size. By using two different grids, the accuracy of the finite difference calculations has been estimated to be of the order of 10^?2 μE _h. The average difference between the finite basis set and finite difference total Hartree-Fock energies is 2.75 μE_h. Dipole moments determined within the algebraic approximation are also compared with the corresponding finite difference expectation values.     

Investigations of spectroscopic parameters and molecular constants for X1Σg+, w3Δu, and W1Δu electronic states of P2 molecule

The potential energy curves (PECs) of three low-lying electronic states (X¹Σ_g⁺, w³Δ_u, and W¹Δ_u)of P₂ molecule are investigated using the full valence complete active space self-consistent field (CASSCF) method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach in conjunction with the correlation-consistent basis set in the valence range. The PECs of the electronic states involved are modified by the Davidson correction and extrapolated to the complete basis set (CBS) limit. With these PECs, the spectroscopic parameters of the three electronic states are determined and compared in detail with the experimental data. The comparison shows that excellent agreement exists between the present results and the available experimental data. The complete vibrational states are computed for the w³Δ_u, and W¹Δ_u electronic states when the rotational quantum number J equals zero and the vibrational level G(v), the inertial rotation constant B_v, and the centrifugal distortion constant D_v of the first 30 vibrational states are reported, which accord well with the experimental data. The present results show that the two-point extrapolation scheme can obviously improve the quality of spectroscopic parameters and molecular constants.     

The theoretical study on the potential energy curve for X~3△ state of TiO molecule

This paper applies the density functional theory method to optimise the structure for X ³Δ state of TiO molecule with the basis sets 6-31G, 6-31++G and 6-311G^**. Comparing the attained results with the experiments, it obtains the conclusion that the basis set 6-31++G is most suitable for the optimal structure calculations of X ³Δ state of TiO molecule. The whole potential energy curve for the electronic state is further scanned by using B3P86/6-31++G method for the ground state, then it uses a least square fitted to Murrell--Sorbie functions, at last it calculates the spectroscopic constants and force constants, which are in better agreement with the experimental data.     

密度泛函方法对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 势能函数 光谱常数     

Potential energy curves for neutral and multiply charged carbon monoxide

Potential energy curves of various electronic states of CO ⁿ⁺ (0 ≤ n ≤ 6) are generated at MRCI/CASSCF level using cc-pvQZ basis set and the results are compared with available experimental and theoretical data.