†Corresponding author. E-mail: chengxl@scu.edu.cn
*Project supported by the National Natural Science Foundation of China (Grant Nos. 11374217 and 11474207).
An accurate electric dipole moment function (EDMF) is obtained for the carbon monoxide (CO) molecule (X1Σ+) by fitting the experimental rovibrational transitional moments. Additionally, an accurate ab initio EDMF is found using the highly accurate, multi-reference averaged coupled-pair functional (ACPF) approach with the basis set, aug-cc-pV6Z, and a finite-field with ±0.005 a.u. (The unit a.u. is the abbreviation of atomic unit). This ab initio EDMF is very consistent with the fitted ones. The vibrational transition matrix moments and the Herman–Wallis factors, calculated with the Rydberg–Klein–Rees (RKR) potential and the fitted and ab initio EDMFs, are compared with experimental measurements. The consistency of these line intensities with the high-resolution transmission (HITRAN) molecular database demonstrates the improved accuracy of the fitted and ab initio EDMFs derived in this work.
Carbon monoxide (CO) is a stable and important diatomic molecule and has been detected in interstellar clouds as well as the Earth’ s atmosphere.[1] The spectral parameters of CO are included in well-known spectroscopic databases such as the high-resolution transmission (HITRAN) molecular database, [2] the high-temperature version of HITRAN called the HITEMP database, [3] and others. For CO and other molecules, calculations of line parameters are closely related to their electric dipole moment functions (EDMFs). With the progress of experimental technology, further investigation of the EDMF for CO (X1Σ + ) based on new experimental data is possible.
Several band systems of CO have been observed. Numerous spectral line parameter measurements of CO have been collected for the fundamental band 1– 0, [4– 7] the first overtone 2– 0, [7– 11] and the second overtone 3– 0.[8, 12– 17] Many studies in the past decade have focused on the 3– 0 band by using a variety of experimental designs. Chackerian and Valero[18] measured the absolute rovibrational intensities of the 4– 0 band of carbon monoxide, and more recently Ogilvie et al.[19] have evaluated the electric dipolar moments based on new experimental measurements of line strengths including the 4– 0 band of 12C16O. Furthermore, Chung et al.[20] adapted a cavity ringdown spectral technique in the near-infrared region to measure the 5– 0 band for the first time.
Several EDMFs for the ground state of CO have been deduced with various methods based on experimental data with Δ υ ≤ 4.[8, 12, 19, 21– 26] For instance, based on the inversion of experimental intensity information, Chackerian and Tipping[22] derived a hybrid EDMF in a Padé approximant. It was used to calculate the rovibrational intensities of the ground state CO (X1Σ + ) by Goorvitch et al.[27, 28, 29] In consideration of these pure vibrational matrix elements of the bands υ ′ – 0 with 0 ≤ υ ′ ≤ 5, Chung et al. derived an EDMF in a polynomial form.[20] A considerable number of quantum chemical calculations of EDMFs for CO have been presented in various papers.[19, 30– 37] However, we find that the vibrational transition matrix elements (or rotationless transition matrix elements) calculated using the quantum chemical EDMFs along with the Rydberg– Klein– Rees (RKR) potential, or ab initio potential, actually disagree with the experimental results. Therefore, more accurate theoretical calculations are needed.
In this work, we adopt a direct fit approach to determine an EDMF for the CO X1Σ + state by utilizing individual, high-quality measurements for line intensities including the 5– 0 band. We also perform a series of theoretical calculations to determine the EDMF of the electronic state of CO (X1Σ + ). Based on the RKR potential, the eigenfunctions, Ψ υ , J(r), are determined by solving the radial Schrö dinger equation of nuclear motion. Eigenfunctions, Ψ υ , J(r), and EDMFs are used to determine the vibrational matrix elements and the Herman– Wallis factors for several bands of CO (X1Σ + ). Finally, we calculate the rovibrational line intensities of some bands of CO (X1Σ + ) by using the present fitted and ab initio EDMFs and compare them with the published experimental data.
The method of obtaining the EDMF experimentally, as developed by Tipping and Herman, is to assume a power series expansion
and to fit the coefficient, Mi, to the measured rotationless transition matrix elements 〈 υ ′ , 0| M(r)| υ ″ , 0〉 .[38] The single and double primes correspond to the upper and lower states, respectively. In Eq. (1), x = (r − re)/re, where re is the equilibrium separation and M0 represents the permanent dipole moment at the minimum potential where r = re. In a harmonic potential, M1x is related to the Δ υ = 1 transitions, while M1x2 is related to the Δ υ = 2 transitions and so on.
Using Eq. (1), the rovibrational matrix elements can be written as
Thus, the inclusion of rotational components allows a better determination of the Mi coefficients by fitting to a greater number of points.[39, 40] In order to derive the new EDMF from Eq. (2), we identify the best available experimental results for each band. We use the results of Zou et al.[7] for the fundamental band, Brulat et al.’ s[41] for the first overtone band, Sung et al.’ s[17] for the 3– 0 band, Chackerian et al.’ s[18] for the 4– 0 band, and Chung et al.’ s[20] for the 5– 0 band. Moreover, in order to obtain the accurate coefficient M0, the value of the matrix moment μ υ ′ = 0, υ ″ = 0 = 0.10982 Debye as reported by Chackerian et al.[6] is also included in the fitting work. The weighted and unweighted fits of all data are performed.
The ground state carbon monoxide molecule is a closed-shell molecule, and the electric configuration of the CO X1Σ + is
The EDMF calculations are performed by using the MOLPRO 2009.1 program package, [44] with and without a finite-field approach with field values of ± 0.005 a.u. We separately use the aug-cc-pV6Z (aV6Z) and spdfghi aug-cc-pV6Z (spdfghi aV6Z) basis sets of Dunning et al. for both C and O atoms.[45– 47] These basis sets are sufficiently diffuse to describe long-range interactions so that a more accurate EDMF can be determined at larger internuclear separations. In this work, we only consider the electric correlation energy of valence electrons and only choose the full valence active spaces consisting of eight valence orbitals, which is denoted by (4220) in C2V point group symmetry.
With the aid of Le Roy’ s LEVEL 8.0 program, [48] the piecewise polynomial interpolation is performed by using numerical input potential energies and electric dipole moments. The rovibrational matrix elements are obtained through the use of the fitted EDMF and the theoretically calculated EDMF, together with an RKR potential energy function. The RKR potential curves for CO are constructed by using the purely numerical RKR method with Dunham coefficients from Velichko et al. who provided accurate molecular parameters for higher-lying vibrational levels (up to υ = 41).[49, 50]
The RKR potential is employed to calculate the transition wavenumbers, lower-state energies, and the expectation values 〈 υ ′ , J′ | xi| 0, J″ 〉 by using Le Roy’ s LEVEL 8.0 program.[48] The rovibrational matrix elements are then evaluated from the measured line intensities according to
where S is the line intensity for pure 12C16O in cm− 1/mol; h, c, and k are Planck constant, velocity of light, and Boltzmann constant, respectively; σ 0 is the transition wavenumber in cm− 1; T is the temperature; Q(T) is the total internal partition; m = 1/2[J′ (J′ + 1) − J″ (J″ + 1)] is a running number of values J″ + 1 for a line in branch R or − J″ for a line in branch P. A sequence of signs of the rovibrational matrix moments (+ − + − + + ) is adopted in a similar way to that in a previous study.[20] Table 1 lists the rovibrational matrix elements derived from the experimental intensities and their experimental uncertainties. Finally, a fit of all listed rovibrational matrix elements to the system of linear Eq. (2) is carried out by using a least-squares fit program.[51] Six coefficients of the EDMF, M0, M1, M2, … , M5, are determined, and the polynomial equation of the weighted EDMF is written as
whereas the unweighted EDMF equation is
Both equations implement an equilibrium separation of re = 1.12822 Å obtained from the constructed RKR potential. For comparison, the determined fitted EDMFs and those reported by Chung et al. are plotted against internuclear distance in Fig. 1.[20] Meanwhile, the theoretical EDMFs calculated through using the ACPF method at the aV6Z set level with the field values of ± 0.005 a.u., those by Langhoff and Bauschlicher[36] through using the aug-cc-pVQZ ACPF-field method, and those by Ogilvie et al.[19] through using the MCSCF method are all plotted in Fig. 1.
![]() | Table 1. Rovibrational transitional moments Mυ 0(m) of CO in units of Debye. The subscripts υ and 0 of Mυ 0(m) represent the υ -th upper vibrational level and the low vibrational level, respectively. |
Within the respective valid range, the three fitted EDMFs are in close agreement with each other. According to the RKR potential, the classical turning points of the vibrational level υ = 5 are found to be approximately 0.967 Å and 1.322 Å . Therefore, the fitted EDMF is most reliable in the region between these two limits, which are slightly different from the results shown in Chung et al. (0.99 Å and 1.33 Å ).[20] The valid range of the fitted EDMF is controlled by the largest overtones which are very sensitive to determining the higher coefficients of Mi in Eq. (1). Thus, it is important to obtain experimental measurements for the largest possible values of υ ′ and Δ υ . The weighted and unweighted fitted EDMFs overlap together at 0.85 Å ≤ r ≤ 1.45 Å in the present work. In the valid range, all the fitted EDMFs agree well with the theoretical EDMF calculated using the ACPF method with a finite-field ± 0.005 a.u at the aV6Z basis set level (aV6Z ACPF-field). They also agree well with the theoretical results of Langhoff and Bauschlicher.[36] The three lines outside the valid range are not consistent with the two theoretical EDMFs at large internuclear distances. Actually, the degree of agreement between the fitted EDMF of Chung et al. and the two theoretical EDMFs become continually worse when r ≥ 1.33 Å .
According to geometrical optimization of the CO molecule, the equilibrium separation and the minimum energy are 1.1324 Å and − 113.197964 a.u. for the aV6Z basis set, and 1.1323 Å and − 113.198196 a.u for the spdfghi aV6Z basis set, respectively. These two values of re agree well with the experimental values of 1.1283 Å and 1.1282 Å obtained from the constructed RKR potential in the present work.[52] The calculated energies and dipole moments each as a function of r are listed in Table 2 at the ACPF (4220) level using the aV6Z and spdfghi aV6Z basis sets. The two EDMFs almost converge both with and without using the finite-field approach. Here, we only compare the aV6Z ACPF-field EDMF with the three fitted and previous theoretical ones in Fig. 1.
![]() | Table 2. Energies and electric dipole moments (EDMs) (a.u.)a for the X1Σ + state CO using the aug– cc– pV6Z and spdfgh aug– cc– pV6Z basis sets. |
The aV6Z ACPF-field EDMF agrees well with the three fitted ones in their valid ranges. Compared with the theoretical results given by Langhoff and Bauschlicher, [36] the values of the dipole moments are smaller at small r and slightly larger at large r. The latest work of Oglivie et al.[19] disagrees with the three fitted and the two theoretical ones.
By combining the RKR potential, the EDMFs from the present weighted fits, and the aV6Z ACPF-field calculations, we calculate the rovibrational matrix moments that are used to compute the line intensities for the υ ′ – 0 (υ ′ = 1– 5) bands, in units of cm/mol at 296 K. The line intensities are then compared with the most current HITRAN (2012) database and measurements of the fundamental band and overtone bands up to υ ′ = 5.[2] The results of the comparisons are plotted in Figs. 3– 7.
In Fig. 3, the calculated line intensities for the 1– 0 band of CO (X1Σ + ) are compared using measurements from Zou et al.[7] and the HITRAN database. The line intensities calculated using the present unweighted fitted EDMF are in better agreement with the experimental results of Zou et al.[7] than the intensities in HITRAN. At − 23 ≤ m ≤ 25, the discrepancies in line intensities between the fitted EDMF and the experimental results of Zou et al. are estimated to be − 2.36%∼ 2.65%, while the discrepancies from the HITRAN comparison with the experimental results range from − 0.38 to − 5.38. The intensities in the HITRAN database are typically lower than the measurements given by Zou et al. The line intensities calculated using the aV6Z ACPF-field EDMF are somewhat larger than the ones of this fitted calculation, the HITRAN database, and Zou et al. Their discrepancies with the measurements are about 1.12%∼ 6.38, which indicates that our ab initio results are also acceptable. For the 1– 0 band, the band strengths calculated using the unweighted fitted and aV6Z ACPF-field EDMFs are 1.027× 10− 17 cm/mol and 1.061× 10− 17 cm/mol, respectively. These values agree well with the value of 1.022± 0.036× 10− 17 cm/mol derived by Varanasi and Sarangi and other data.[4, 6, 7, 10]
![]() | Fig. 3. Comparisons of the present calculated line intensities with the measurements and the HITRAN database for the band 1– 0 of the CO. |
![]() | Fig. 4. Comparisons of the present calculated line intensities with the measurements and the HITRAN database for the band 2– 0 of the CO. |
Figure 4 illustrates the comparison for the 2– 0 band. The results show that the line intensities of the unweighted fit agree well with the HITRAN database and the experimental line intensities of Zou et al. and Devi et al.[2, 7, 10, 11] The line intensities from the aV6Z field EDMF are slightly larger than the measurements of Zou et al. and Devi et al., [7, 10, 11] and their discrepancies are typically within ± 3. The band strengths that are obtained from the unweighted fitted and aV6Z ACPF-field EDMFs are 7.649× 10− 20 cm/mol and 7.768× 10− 20 cm/mol. Both agree well with the value of 7.466± 0.002× 10− 20 cm/mol computed using a fitted vibrational matrix moment and the experimental value of 7.532± 0.066× 10− 20 cm/mol.[7, 41]
The measurements for the 3– 0 band given in Refs. [2], [13]– [15], and [17] are compared with the measurements made in this work and with those in the HITRAN database as shown in Fig. 5. The line intensities from the present study and HITRAN database[2] are consistent with the 3– 0 band measurements presented in the literature. While the results of the unweighted fitted EDMFs and the aV6Z ACPF-field EDMFs are within 1% in agreement with the values provided by Chackerian et al.[15] and Sung et al., [17] and the line intensities in the HITRAN [2] database correspond more closely to those of Henningsen et al.[13] and Picqué et al.[14]
![]() | Fig. 5. Comparisons of the present calculated line intensities with the measurements and the HITRAN database for the band 3– 0 of the CO. |
Though Ogilvie et al. measured the strengths of lines for vibration– rotational bands of 12C16O, including the 4– 0 band, the parameters of the line intensities are not listed in their paper.[19] Therefore, only the early experimental data of Chackerian and Valero[18] are used for fits and comparison. These measurements as shown in Fig. 6, are compared with the calculated line intensities of the 4– 0 band and the corresponding values in the HITRAN database.[2] The line intensities of the P and R branches as determined with the unweighted fit, are consistent with the referenced experimental results. Furthermore, the line intensities of the aV6Z field EDMF of the P branch are in better agreement with the experimental results than those of the R branch. In contrast, the line intensities from the HITRAN database match better with the measurements in the R branch.[2] However, the values from the unweighted fitted EDMF, the aV6Z ACPF-field EDMF, and HITRAN do not agree with the experimental results when considering the intensities of the P(15), R(20), and R(23) lines. The discrepancy of the R(23) line exceeds 20%, and it may arise from the large error of the measurement. Still, values from the two EDMFs and HITRAN[2] are consistent with each other.
The line intensity measurements for the 5– 0 overtone band are only performed by Chung et al. who merely presented 16 lines of the P branch with an overly high reported uncertainty.[20] The minimum uncertainty of the line strength is likely to be 15%. The comparisons of line intensities for the 5– 0 band appear in Fig. 7. There are rather large discrepancies between the measurements and calculations of Chung et al., i.e., the line intensities of the unweighted fit and of the aV6Z ACPF-field EDMF.[20] The unweighted fit line intensities agree better with the results using the aV6Z ACPF-field EDMF. The calculated results again show that the line intensities in branch P are truly greater than those in branch R. For the 5– 0 band, the band strengths obtained from the weighted fitted and aV6Z ACPF-field EDMFs are 5.585× 10− 27 cm/mol and 5.674× 10− 27 cm/mol, which are in close agreement with the experimental results of 5.1± 1.3× 10− 27 cm/mol.[20] Considering the rather large experimental uncertainties, we argue that the present intensities of the unweighted fitted line and the aV6Z ACPF-field EDMFs are still reliable and reflect the basic qualitative characteristics of spectral lines.
![]() | Fig. 6. Comparisons of the present calculated line intensities with the measurements and the HITRAN database for the band 4– 0 of the CO. |
![]() | Fig. 7. Comparisons of the line intensities of the CO 5– 0 band with the measurements and the calculations of Chung et al.[20] |
The comparisons of line intensity, presented here, between the values from the HITRAN database and experimental measurements clearly show that the EDMF can be directly constructed from the experimental rovibrational transition matrix elements and can provide accurate values for CO (X1Σ + ). Although the line intensities calculated using the aV6Z ACPF-field EDMF are greater than the measurements at the υ ′ – 0 (υ ′ = 1– 5) band, the discrepancies are within the admissible range of deviation between theoretical calculation and measurement.
By using a weighted and unweighted approach to fit the best available data of individual rovibrational transitions of CO (X1Σ + ), we obtain an accurate EDMF within the valid range of 0.967 Å ≤ r ≤ 1.322 Å , which agrees well with that of Chung et al.[20] Meanwhile, an accurate ab initio EDMF is also calculated at the ACPF/aug-cc-pV6Z level with a finite-field of ± 0.005 a.u. The fitted EDMF in the valid range shows close agreement with the ab initio EDMF, and for small values of r, the electric dipole moment values are slightly smaller than the previous theoretical results of Langhoff and Bauschlicher.[36] Combining the RKR potential curve and these EDMFs, the transition matrix moments are computed and compared with the experimental ones. The vibrational transition matrix moments μ υ ′ , υ ″ calculated by using our unweighted fitted EDMF are in excellent agreement with the experimental values. The Herman– Wallis factors calculated using the present aV6Z ACPF-field EDMF also agree well with the experimental measurements. Comparisons are made among line intensities calculated here, reported in the HITRAN[2] database and obtained from experimental measurements. The results clearly show that the present unweighted fitted and aV6Z ACPF-field EDMFs are accurate and useful in predicting the spectral properties of the CO molecule.
1 |
|
2 |
|
3 |
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|
17 |
|
18 |
|
19 |
|
20 |
|
21 |
|
22 |
|
23 |
|
24 |
|
25 |
|
26 |
|
27 |
|
28 |
|
29 |
|
30 |
|
31 |
|
32 |
|
33 |
|
34 |
|
35 |
|
36 |
|
37 |
|
38 |
|
39 |
|
40 |
|
41 |
|
42 |
|
43 |
|
44 |
|
45 |
|
46 |
|
47 |
|
48 | [Cited within:2] |
49 | [Cited within:1] |
50 |
|
51 |
|
52 |
|
53 |
|