Theoretical investigations of spectroscopic parameters and molecular constants for electronic ground state of Cl2 and its isotopes

The potential energy curve of the Cl 2(X1Σg+) is investigated by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach in combination with the largest correlation-consistent basis set, aug-cc-pV6Z, in the valence range. The theoretical spectroscopic parameters and the molecular constants of three isotopes, 35Cl2, 35Cl37Cl and 37Cl2, are studied. For the 35Cl2 (X1Σg+), the values of D0 , De , Re , ωe , ωeχe , αe and Be are obtained to be 2.3921 eV, 2.4264 eV, 0.19939 nm, 555.13 cm-1 , 2.6772 cm-1 , 0.001481 cm-1 and 0.24225 cm-1 , respectively. For the 35Cl37Cl(X1Σg+), the values of D0 , De , Re , ωe , ωeχe , αe and Be are calculated to be 2.3918 eV, 2.4257 eV, 0.19939 nm, 547.68 cm-1 , 2.6234 cm-1 , 0.00140 cm-1 and 0.23572 cm-1 , respectively. And for the 37C2 (X1Σg+), the values of D0 , De , Re , ωe , ωeχe , αe and Be are obtained to be 2.3923 eV, 2.4257 eV, 0.19939 nm, 540.06 cm-1 , 2.5556 cm-1 , 0.00139 cm-1 and 0.22919 cm-1 , respectively. These spectroscopic results are in good agreement with the available experimental data. With the potential of Cl2 molecule determined at the MRCI/aug-cc-pV6Z level of theory, the total of 59 vibrational states is predicted for each isotope when the rotational quantum number J equals zero (J = 0). The theoretical vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants are determined when J = 0, which are in excellent accordance with the available experimental findings.     

Elastic scattering of two H（^2Sg） and N（^4Su） atoms at low temperatures and accurate spectroscopic parameters of NH （X^3∑^-） radical

This paper reports that the interaction potential for the X3Z- state of NH radical is constructed at the CCSD（T）/ cc-PV6Z level of theory. Using this potential, this paper calculates the spectroscopic parameters （De, Re, ωe, ωeχe, αe and Be） and their values are of 3.578eV, 0.10368nm, 3286.833cm^-1, 78.433cm^-1, 0.6469cm^-1 and 16.6735cm^-1 respectively, which are in excellent agreement with the experiments. Then the total of 14 vibrational states has been found when J=0 by solving the radial Schrodinger equation of nuclear motion. For each vibrational state, the vibrational manifolds are reported for the first time. And last, the total cross sections, s-wave, p-wave and d-wave cross sections are computed for the elastic collisions between two ground-state atoms （hydrogen and nitrogen） at low temperatures. It finds that the total elastic cross sections are dominated by s-wave scattering when the collision energy is below 10^-6a.u. The pronounced shape resonance is found at energy of 6.1 × 10^-6a.u. Calculations have shown that the shape resonance comes from the p-wave contributions.     

Investigation of analytical harmonic frequency and potential energy function,vibrational levels for the X^2∑^＋ and A^2Л states of CN radical

This paper calculates the equilibrium structure and the potential energy functions of the ground state （X^2∑^＋） and the low lying excited electronic state （A^2Л） of CN radical are calculated by using CASSCF method. The potential energy curves are obtained by a least square fitting to the modified Murrell-Sorbie function. On the basis of physical theory of potential energy function, harmonic frequency （ωe） and other spectroscopic constants （ωeχe, βe and αe） are calculated by employing the Rydberg-Klei-Rees method. The theoretical calculation results are in excellent agreement with the experimental and other complicated theoretical calculation data. In addition, the eigenvalues of vibrational levels have been calculated by solving the radial one-dimensional SchrSdinger equation of nuclear motion using the algebraic method based on the analytical potential energy function.     

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(B e,α e and ω e χ e) of these states.For the X 2 Σ + state,the obtained values of D e,B e,α e,ω e χ e,R e 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 A 2Ⅱ state,the present values of D e,B e,α e,ω e χ e,R e and ω e 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 B 2 Σ + state,the present values of D e,B e,α e,ω e χ e,R e and ω e 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 Schro¨dinger 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.     

Ab initio calculations on the spectroscopic constants, vibrational levels and classical turning points for the 21∏u state of dimer 7Li2

The accurate dissociation energy and harmonic frequency for the highly excited 2^1Пu state of dimer ^7Li2 have been calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space. The calculated results are in excellent agreement with experimental measurements. The potential energy curves at numerous basis sets for this state are obtained over a wide internuclear separation range from about 2.4a0 to 37.0a0. And the conclusion is gained that the basis set 6-311＋＋G（d,p） is a most suitable one. The calculated spectroscopic constants De, Re, ωe, ωeχe, ae and Be at 6-311＋＋G（d,p） are 0.9670 eV, 0.3125 nm, 238.6 cm^-1, 1.3705 cm^-1, 0.0039 cm^-1 and 0.4921 cm^-1, respectively. The vibrational levels are calculated by solving the radial SchrSdinger equation of nuclear motion. A total of 53 vibrational levels are found and reported for the first time. The classical turning points have been computed. Comparing with the measurements, in which only the first nine vibrational levels have been obtained so far, the present calculations are very encouraging. A careful comparison of the present results of the parameters De and We with those obtained from previous theories clearly shows that the present calculations are much closer to the measurements than previous theoretical results, thus representing an improvement on the accuracy of the ab initio calculations of the potentials for this state.     

Ab initio calculation on accurate analytic potential energy functions and harmonic frequencies of c^3∑g^＋ and B^1-Пu states of dimer 7Li2

The comparison between single-point energy scanning （SPES） and geometry optimization （OPT） in determining the equilibrium geometries of c^3∑g^＋ and B^1-Пu states of dimer 7Li2 is made at numerous basis sets by using a symmetryadapted-cluster configuration-interaztion （SAC-CI） method in the Gaussian 03 program package. In this paper the difference of the equilibrium geometries obtained by SPES and by OPT is reported. The results obtained by SPES are found to be more reasonable than those obtained by OPT in full active space at the present SAC-CI level of theory. And the conclusion is attained that the cc-PVTZ is a most suitable basis set for these states. The calculated dissociation energies and equilibrium geometries are 0.8818 eV and 0.3090 nm for c^3∑g^＋ state, and 0.3668 eV and 0.2932 nm for B^1-Пu state respectively. The potential energy curves are calculated over a wide internuclear distance range from about 2.5α0 to 37α0 and have a least-squares fit into the Murrell-Sorbie function. According to the calculated analytic potential energy functions, the harmonic frequencies （We） and other spectroscopic data （ωeXe, Be and αe） are calculated. Comparison of the theoretical determinations at present work with the experiments and other theories clearly shows that the present work is the most complete effort and thus represents an improvement over previous theoretical results.     

Elastic collisions of sulfur and hydrogen in their ground states at low temperatures and spectroscopic parameters of SH（X^2∏） radical

This paper constructs the interaction potential of the SH（X^2∏） radical by using the coupled-cluster singlesdoubles-approximate-triples theory combining the correlation-consistent quintuple basis set augmented with the diffuse functions, aug-cc-pV5Z, in the valence range. Employing the potential, it accurately determines the spectroscopic parameters. The present De, Re, ωe, ωeχe, ae and Be values are of 3.7767eV, 0.13424nm, 2699.846 cm^-1, 47.7055 cm^-1, 0.2639cm^-1 and 9.4414 cm^-1, respectively, which are in excellent agreement with those obtained from the measure- ments. A total of 19 vibrational states has been found when J = 0 by solving the radial SchrSdinger equation of nuclear motion. The complete vibrational levels, classical turning points, initial rotation and centrifugal distortion constants when J = 0 are reported for the first time, which are in good accord with the experimental results. The total and various partial-wave cross sections are computed for the elastic collisions of sulfur and hydrogen in their ground states at low temperatures when two atoms approach each other along the SH（X^2∏） potential energy curve. Over the impact energy range from 1.0×10^-11 to 1.0×10^-4 a.u., eight shape resonances have been found in the total elastic cross sections. For each shape resonance, the resonant energy is accurately calculated. Careful investigations have pointed out that these resonances result from the 1 = 0, 1, 2, 3, 4, 6, 7, 8 partial-wave contributions.     

Structural, Spectroscopic and Vibrational Properties for Low-Lying Electronic States of LiCl

A multireference configuration interaction(MRCI) study has been carried out on the LiCl molecule. The potential energy has been calculated over a wide range of internuclear separation for the 21 low-lying electronic states of the LiCl molecule dissociating into Li(2S,2P,3S)+Cl(2P). The(4)1Σ+,(3)1Π, 1-33Σ+, 1-33Π,1,3,1,3Σ,(5)1Σ+,(4)3Σ+,(4)1Π,(4)3Π excited states are studied for the first time in theory. Molecular spectroscopic constants(Re, De,ωe, ωeχe, Be and αe) have been derived for the 9 bound states(X1Σ+,(3)1Σ+,(2)3Σ+,1,3,1,3Σ,(4)1Π,(4)3Π) with a regular shape, and the spectroscopic constants of ground states X1Σ+are in good agreement with available experimental and theoretical values. The relative diferences between experimental values and our values for Re, De, ωe, ωeχe, Be andαe are 1.02%, 0.60%, 1.72%, 9.46%, 2.0%, and 0.75%, respectively. Moreover, vibrational levels of 9 bound states, which have not been investigated experimentally, are computed.     

Elastic collisions between Si and D atoms at low temperatures and accurate analytic potential energy function and molecular constants of the SiD（χ^2П） radical

Interaction potential of the SiD(X2Π) radical is constructed by using the CCSD(T) theory in combination with the largest correlation-consistent quintuple basis set augmented with the diffuse functions in the valence range. Using the interaction potential, the spectroscopic parameters are accurately determined. The present D0, De, Re, ωe, αe and Be values are of 3.0956 eV, 3.1863 eV, 0.15223 nm, 1472.894 cm-1, 0.07799 cm-1 and 3.8717 cm-1, respectively, which are in excellent agreement with the measurements. A total of 26 vibrational states is predicted when J=0 by solving the radial Schro¨dinger equation of nuclear motion. The complete vibrational levels, classical turning points, initial rotation and centrifugal distortion constants when J=0 are reported for the first time, which are in good accord with the available experiments. The total and various partial-wave cross sections are calculated for the elastic collisions between Si and D atoms in their ground states at 1.0×10-11–1.0×10-3 a.u. when the two atoms approach each other along the SiD(X2Π) potential energy curve. Four shape resonances are found in the total elastic cross sections, and their resonant energies are of 1.73×10-5, 4.0×10-5, 6.45×10-5 and 5.5×10-4 a.u., respectively. Each shape resonance in the total elastic cross sections is carefully investigated. The results show that the shape of the total elastic cross sections is mainly dominated by the s partial wave at very low temperatures. Because of the weakness of the shape resonances coming from the higher partial waves, most of them are passed into oblivion by the strong s partial-wave elastic cross sections.     

Potential energy curves and analytical potential energy functions of the metastable states of B2^＋＋

The multi-reference configuration interaction method and aug-cc-pvqz （AVQZ） have been used to calculate potential energy curves （PECs） of the singlet and triplet states of the riu and rig symmetry of B2＋＋. All of the four states （^l∏u, ^1∏g, ^3∏u and ^3∏g） are found to be metastable states, though the potential well of ^3∏u symmetry is very shallow. Based on the PECs, the analytical potential energy functions （APEFs） of these states have been fitted using the least square fitting method and two models of function. The spectroscopic parameters of each state are also calculated, and are compared with other investigations in the literature. The credibility and veracity of the two functions are evaluated. Some ideas to improve the fitting accuracy are presented. Also the vibrational levels for each state are predicted by solving the SchrSdinger equation of nuclear motion.     

Molecular constants of LiCI（X^1∑＋） and elastic collisions of two ground-state Cl and Li atoms at low and ultralow temperatures

Interaction potentials for LiCI（X^1∑＋） are constructed by the highly accurate valence internally contracted multireference configuration interaction in combination with a number of large correlation-consistent basis sets, which are used to determine the spectroscopic parameters （D0, De, Re, ωe, ωeχe, Be and αe）. The potentials obtained at the basis sets, i.e., aug-cc-pV5Z-JKFI for Cl and cc-pV5Z for Li, are selected to study the elastic collision properties of Li and Cl atoms at the impact energies from 1.0 ×10^-12 to 1.0× 10-4 a.u. The derived total elastic cross sections are very large and almost constant at ultralow temperatures, and their shapes are mainly dominated by the s-partial wave at very low impact energies. Only one shape resonance can be found in the total elastic cross sections over the present collision energy regime, which is rather strong and obviously broadened by the overlap contributions of the abundant resonances coming from various partial waves. Abundant resonances exist for the elastic partial-wave cross sections until l= 22 partial waves. The vibrational manifolds of the LiCI（X^1∑＋） molecule, which are predicted at the present level of theory and the basis sets cc-pV5Z for Li and the aug-cc-pV5Z-JKFI for Cl, should achieve much high accuracy due to the employment of the large correlation-consistent basis sets.     

Spectroscopic investigations on NO＋（X1∑＋, a3∑＋, A1II） ion using multi-reference configuration interaction method and correlation-consistent sextuple basis set augmented with diffuse functions

Three low-lying electronic states （X1∑, a3∑＋, and A1II） of NO＋ ion are studied using the complete active space self-consistent-field （CASSCF） method followed by highly accurate valence internally contracted multi-reference configuration interaction （MRCI） approach in combination of the correlation-consistent sextuple basis set augmented with diffuse functions, aug-cc-pV6Z. The potential energy curves （PECs） of the NO＋（X1∑＋, a3∑＋, A1II） are calculated. Based on the PECs, the spectroscopic parameters Re, De, We, WeXe, ae, Be, and D0 are reproduced, which are in excellent agreement with the available measurements. By numerically solving the radial SchrSdinger equation of nuclear motion using the Numerov method, the first 20 vibrational levels, inertial rotation and centrifugal distortion constants of NO＋（X1∑＋, a3∑＋, A1II） ion are derived when the rotational quantum number J is equal to zero （J = 0） for the first time, which accord well with the available measurements. Finally, the analytical potential energy functions of these states are fitted, which are used to accurately derive the first 20 classical turning points when J = 0. These results are compared in detail with those of previous investigations reported in the literature.     

Ab initio calculations on the α^3∑u^＋ state properties of dimer ^7Li2

The comparison between single-point energy scanning （SPES） and geometry optimization （OPT） in determining the equilibrium geometry of the α^3∑u^＋ state for ^7Li2 is made at numerous basis sets such as 6-311＋＋G（2df）, cc-PVTZ, 6-311＋＋G（2df, p）, 6-311G（3df,3pd）, 6-311＋＋G（2df,2pd）, D95（3df,3pd）, 6-311＋＋G, DGDZVP, 6-311＋＋G（3df,2pd）, 6-311G（2df,2pd）, D95V＋＋, CEP-121G, 6-311＋＋G（d,p）, 6-311＋＋G（2df, pd） and 6-311＋＋G（3df,3pd） in full active space using a symmetry-adapted-cluster/ symmetry-adapted-cluster configuration-interaction （SAC/SAC=CI） method presented in Gaussian03 program package. The difference of the equilibrium geometries obtained by SPES and by OPT is reported. Analyses show that the results obtained by SPES are more reasonable than those obtained by OPT. We have calculated the complete potential energy curves at those sets over a wide internuclear distance range from about 3.0α0 to 37.0α0, and the conclusion is that the basis set cc-PVTZ is the most suitable one. With the potential obtained at ccopVTZ, the spectroscopic data （Te, De, D0, ωe,ωeХe, αe and Be） are computed and they are 1.006 eV, 338.71 cm^-1, 307.12 cm^-1, 64.88 cm^-1, 3.41 cm^-1, 0.0187 cm^-1 and 0.279 cm^-1, respectively, which are in good agreement with recent measurements. The total 11 vibrational states are found at J=0. Their corresponding vibrational levels and classical turning points are computed and compared with available RKR data, and good agreement is found. One inertial rotation constant （By） and six centrifugal distortion constants （Dr Hv, Lv, My, Nv, and Ov） are calculated. The scattering length is calculated to be -27.138α0, which is in good accord with the experimental data.     

Investigations on molecular constants of the CD（X^2∏） radical and elastic collisions between ground-state C and D atoms at low temperatures

The potential energy curve of the CD(X~2Π) radical is obtained using the coupled-cluster singles-doublesapproximate-triples [CCSD(T)] theory in combination with the correlation-consistent quintuple basis set augmented with diffuse functions,aug-cc-pV5Z.The potential energy curve is fitted to the Murrell-Sorbie function,which is used to determine the spectroscopic parameters.The obtained D0,De,Re,ωe,ωeχe,αe and Be values are 3.4971 eV,3.6261 eV,0.11197 nm,2097.661 cm 1,34.6963 cm 1,0.2083 cm 1 and 7.7962 cm 1,respectively,which conform almost perfectly to the available measurements.With the potential obtained at the UCCSD(T)/aug-cc-pV5Z level of theory,a total of 24 vibrational states have been predicted for the first time when J = 0 by solving the radial Schr¨odinger equation of nuclear motion.The complete vibrational levels,the classical turning points,the inertial rotation constants and centrifugal distortion constants are reproduced from the CD(X~2Π) potential when J = 0,and are in excellent agreement with the available measurements.The total and the various partial-wave cross sections are calculated for the elastic collisions between the ground-state C and D atoms at energies from 1.0×10 11 to 1.0×10 4 a.u.when the two atoms approach each other along the CD(X~2Π) potential energy curve.Only one shape resonance is found in the total elastic cross sections,and the resonant energy is 8.36×10 6 a.u.The results show that the shape of the total elastic cross section is mainly dominated by the s partial wave at very low temperatures.Because of the weak shape resonances coming from higher partial waves,most of them are passed into oblivion by the strong total elastic cross sections.     

Accurate Analytic Potential Energy Function and Spectroscopic Study for G^1Ⅱg State of Dimer ^7Li2

The reasonable dissociation limit for the G^1Ⅱg, state of dimer ^7Li2 is determined. The equilibrium internuclear distance, dissociation energy, harmonic frequency, vibrational zero energy, and adiabatic excitation energy are calculated using a symmetry-adapted-cluster configuration-interactlon method in complete active space in Gaussian03 program package at such numerous basis sets as 6-311 ＋＋G, 6-311 ＋＋G（2df,2pd）, 6-311 ＋＋G（2df, p）, cc-PVTZ, 6- 311＋＋G（3df,3pd）, CEP-121G, 6-311＋＋G（2df, pd）, 6-311＋＋G（d,p）,6-311G（3df,3pd）, D95（3df,3pd）, 6-311＋＋G（3df,2p）, 6-311＋＋G（2df）, 6-311＋＋G（df, pd） D95V＋＋, and DGDZVP. The complete potential energy curves are obtained at these sets over a wide internuclear distance range and have least squares fitted to Murrell-Sorbie function. The conclnsion shows that the basis set 6-311＋＋G（2df, p） is a most suitable one for the G^1Ⅱg state. At this basis set, the calculated spectroscopic constants Te, De, Eo, Re, ωe, ωeXe, ae, and Be are of 3.9523 eV, 0.813 06 eV, 113.56 cm^-1, 0.320 15 nm, 227.96 cm^-1, 1.6928 cm^-1, 0.004 436 cm^-1, and 0.4689 cm^-1, respectively, which are in good agreement with measurements whenever available. The total 50 vibrational levels and corresponding inertial rotation constants are for the first time calculated and compared with available RKR data. And good agreement with measurements is obtained.     

Structure and analytical potential energy function for the ground state of the BCx (x=0, －1)

In this paper, the electronic states of the ground states and dissociation limits of BC and BC- are correctly determined based on group theory and atomic and molecular reaction statics. The equilibrium geometries, harmonic frequencies and dissociation energies of the ground state of BC and BC- are calculated by using density function theory and quadratic CI method including single and double substitutions. The analytical potential energy functions of these states have been fitted with Murrell-Sorbie potential energy function from our ab initio calculation results. The spectroscopic data （αe, ωe and ωeχe） of each state is calculated via the relation between analytical potential energy function and spectroscopic data. All the calculations are in good agreement with the experimental data.     

Elastic Collisions Between two Ground-State P and D Atoms at Low and Ultralow Temperatures

The PD（X^3∑^-） interaction potential is constructed using the CCSD（T） theory and the basis set, augcc-pV5Z. Using this potential, the spectroscopic parameters are accurately determined. The present Do, De, Re, ωe, ωeχe, αe, and Be are of 3.056 99 eV, 3.161 75 eV, 0.142 39 nm, 1701.558 cm^-1, 23.6583 cm^-1, 0.085 99 cm^-1, and 4.3963 cm^-1, respectively, which almost perfectly conform with the measurements. A total of 26 vibrational states is predicted when J = 0 by solving the radial Sehrodinger equation of nuclear motion. The complete vibrational levels, classical turning points, initial rotation and centrifugal distortion constants when J = 0 are reported for the first time, which favorably agree with the experiments. The total and various partial-wave cross sections are calculated for the elastic impact between two ground-state P and D atoms at 1.0 × 10^-12 - 1.0 × 10^-4 a.u. when they approach each other along the PD（X^3∑^-） potential. No shape resonances exist in the total elastic cross sections, though the peaks can be found for each partial wave until l=6. The shape of the total elastic cross sections is dominated by the s partial wave at very low temperatures. Due to the weakness of the shape resonances of each partial wave, they are all passed into oblivion by the strong total elastic cross sections.     

Accurate studies on the full vibrational energy spectra and molecular dissociation energies for some electronic states of N_2 molecule

1 IntroductionAccurate results of high-lying vibrational energies and molecular dissociation energy De of diatomic electronic states are very important for thermodynamics, molecular spec- troscopy, reactive scatterings, and the collision physics of ultracold atoms. For example, the binding energy of the highest bound vibrational state to the ground vibrational state determines the s-wave scattering length. This in turn determines the low- energy (pre- dominantly s-wave) atomic elastic-scatteri…     

MRCI potential energy curves and analytical potential energy functions for the X^2∑^＋ and 2^∏ states of BO molecule

The potential energy curves （PECs） of BO molecule, including ∑^＋and ∏ symmetries with doublet spin multiplicities, are obtained employing multi-reference configuration interaction （MRCI） method and Dunning＇s correlation consistent basis sets. The analytical potential energy functions （APEFs） are fitted using the Murrell-Sorbie （MS） function and the least square method. Based on the PECs, the spectroscopic constants of the states have been determined and compared with the theoretical and experimental results available to affirm the accuracy and liability of the calculations. The root-mean-square （RMS） errors between the fitted results and the ab initio values are too little in comparison with the chemical accuracy （349.755 cm^-1）. It is shown that the present APEFs are accurate and can display the interaction between the atoms well. The present APEFs can be used to construct more complicated APEF or do some dynamic investigations.     

The theoretical study on the potential energy curves for X^1 ∑^＋, A^1 ⅡI and C^1 ∑^- states of SiO molecule

This paper applies the symmetry-aziapted-cluster/symmetry-adapted-cluster configuration-interaction （SAC/SACCI） method to optimize the structures for X^1∑^＋, A^1 Ⅱ and C^1 ∑^- states of SiO molecule with the basis sets D95＋＋, 6-311＋＋G and 6-311＋＋G^＊＊. Comparing the obtained results with the experiments, it gets the conclusion that the basis set 6-311＋＋G^＊＊ is most suitable for the optimal structure calculations of X^1.∑^＋, A^Ⅱ and C^1∑^- states of SiO molecule. The whole potential energy curves for these electronic states are further scanned by using SAC/6-311＋＋G^＊＊ method for the ground state and SAC-CI/6-311＋＋G^＊＊ method for the excited states, then use a least square method to fit Murrell~Sorbie functions, at last the spectroscopic constants and force constants are calculated, which are in good agreement with the experimental data.