Theoretical study of spectroscopic and molecular properties of several low‐lying electronic states of CO molecule

The potential energy curves (PECs) of eight low‐lying electronic states (X¹Σ⁺, a³Π, a′³Σ⁺, d³Δ, e³Σ^?, A¹Π, I¹Σ^?, and D¹Δ) of the carbon monoxide molecule have been studied by an ab initio quantum chemical method. The calculations have been performed using the complete active space self‐consistent field method, which is followed by the valence internally contracted multireference configuration interaction (MRCI) approach in combination with the correlation‐consistent aug‐cc‐pV5Z basis set. The effects on the PECs by the core‐valence correlation and relativistic corrections are included. The way to consider the relativistic corrections is to use the third‐order Douglas–Kroll Hamiltonian approximation at the level of a cc‐pV5Z basis set. Core‐valence correlation corrections are performed using the cc‐pCVQZ basis set. To obtain more reliable results, the PECs determined by the MRCI calculations are corrected for size‐extensivity errors by means of the Davidson modification (MRCI+Q). The spectroscopic parameters (D_e, T_e, R_e, ω_e, ω_ex_e, ω_ey_e, B_e, α_e, and γ_e) of these electronic states are calculated using these PECs. The spectroscopic parameters are compared with those reported in the literature. Using the Breit–Pauli operator, the spin–orbit coupling effect on the spectroscopic parameters is discussed for the a³Π electronic state. With the PECs obtained by the MRCI+Q/aug‐cc‐pV5Z+CV+DK calculations, the complete vibrational states of each electronic state have been determined. The vibrational manifolds have been calculated for each vibrational state of each electronic state. The vibrational level G(ν), inertial rotation constant B_ν, and centrifugal distortion constant D_ν of the first 20 vibrational states when the rotational quantum number J equals zero are reported and compared with the experimental data. Comparison with the measurements demonstrates that the present spectroscopic parameters and molecular constants determined by the MRCI+Q/aug‐cc‐pV5Z+CV+DK calculations are both reliable and accurate. © 2012 Wiley Periodicals, Inc.     

Spectroscopic constants and molecular properties of A3Σu+, B3Πg,W3Δu,and B′3Σu− electronic states of the N2 molecule

The potential energy curves (PECs) of A³Σ, B³Π_g, W³Δ_u, and B^′3Σ electronic states of the N₂ molecule have been studied for internuclear separations from 0.05 to 2.0 nm using the full valence complete active space self‐consistent‐field method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach in conjunction with the correlation‐consistent basis sets. Effects on the PECs by the core–valence correlation and relativistic corrections are taken into account. The way to consider the relativistic correction is to use the second‐order Douglas‐Kroll Hamiltonian approximation. The core–valence correlation correction is made with the cc‐pCV5Z basis set. And the relativistic correction is performed at the level of cc‐pV5Z basis set. To obtain more reliable results, the PECs determined by the MRCI calculations are also corrected for size‐extensivity errors by the Davidson modification (MRCI+Q). These PECs are extrapolated to the complete basis set limit. The spectroscopic parameters of ¹⁴N₂, ¹⁴N¹⁵N, and ¹⁵N₂ isotopologs have been evaluated and compared with those reported in the literature. Excellent agreement has been found between the present results and the Rydberg‐Klein‐Rees (RKR) data. With the PECs obtained by the MRCI+Q/CV+DK+56 calculations, the first 30 vibrational states for three species are computed for each electronic state. And for each electronic state of each species, the vibrational level G(ν), inertial rotation constant B_ν, and centrifugal distortion constant D_ν have been determined, which agree well with the RKR data. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Spectroscopic Parameters of X3∑-, a1△, and A'3△ Electronic States of SO Radical

The potential energy curves (PECs) of three low-lying electronic states (X³∑^-, a¹△, and A^'3△) of SO radical have been studied by ab initio quantum chemical method. The calcula-tions were carried out with the full valence complete active space self-consistent field method followed by the highly accurate valence internally contracted multireference configuration in-teraction (MRCI) approach in combination with correlation-consistent basis sets. Effects of the core-valence correlation and relativistic corrections on the PECs are taken into account. The core-valence correlation correction is carried out with the cc-pCVDZ basis set. The way to consider the relativistic correction is to use the second-order Douglas-Kroll Hamiltonian approximation, and the correction is performed at the level of cc-pV5Z basis set. To obtain more reliable results, the PECs determined by the MRCI calculations are also corrected for size-extensivity errors by means of the Davidson modification (MRCI+Q). These PECs are extrapolated to the complete basis set limit by the two-point energy extrapolation scheme. With these PECs, the spectroscopic parameters are determined.     

Accurate spectroscopic constants of the lowest three electronic states in halonitrenes with multireference configuration interaction

Internally contracted multireference configuration interaction (icMRCI) calculations of the ground state (X³Σ⁻), the first excited state (a¹Δ) as well as the second excited state (b¹Σ⁺) have been performed for a series of halogenated nitrenes NXs (X = Cl, Br, and I). Accurate spectroscopic constants of these lowest three electronic states of each NX were obtained in this work using MRCI methods with aug‐cc‐pVXZ (X = T, Q, 5) basis sets and complete basis set (CBS) limit. In addition, various corrections, including the Davidson correction, scalar relativistic effect, core‐valence correlation, and spin‐orbit coupling effect, have been studied in calculating spectroscopic constants, especially for heavy‐atom nitrenes. Comparisons have been made with previous computational and experimental results where available. The icMRCI + Q calculations presented in this work provide a comprehensive series of results at a consistent high level of theory for all of the halogenated nitrenes.     

Spectroscopic parameter and molecular constant investigations for low‐lying electronic states of P ion

The potential energy curves (PECs) of three low‐lying electronic states of P ion, X²Π_u, A²Σ, and B²Σ, have been studied using the full valence complete active space self‐consistent field method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach and MRCI with Davidson correction (+Q). The correlation‐consistent basis sets, aug‐cc‐pV5Z and aug‐cc‐pV6Z, are used and the total energies are extrapolated to the complete basis set limit. Using these PECs obtained with the MRCI+Q/56‐extrapolation, the spectroscopic parameters for these electronic states are determined and compared in detail with experimental data and those of previous studies reported in the literature. The comparison shows that excellent agreement exists between the present results and the available experiments. The first 40 vibrational states for the three electronic states are also computed when the rotational quantum number J equals zero. For each vibrational state, the vibrational level G(υ), inertial rotation constant B_υ, and centrifugal distortion constant D_υ are determined when J = 0, which are in good accord with the available measurements. © 2012 Wiley Periodicals, Inc.     

Theoretical determination of the X 2Σ+ and A 2Π potentials of CsO using relativistic effective core potentials

Theoretical potential energy curves are computed for the X ²Σ⁺ and A ²Π states of CsO using a relativistic effective core potential and a large valence Gaussian basis set. Seventeen electrons are correlated by a CI (SD ) calculation from each HF reference. We find the X ²Σ⁺ state lower by 497 and 726 cm^?1 at the HF and CI(SD) levels. Our calculated ω_e of 312 cm^?1 for the X ²Σ⁺ state agrees well with experimental values deduced from studies in matrices.     

One‐electron pseudopotential study of the alkali hydride cation NaH+: Structure,spectroscopy, transition dipole moments,and radiative lifetimes

The structure and spectroscopic properties of the ground and the lowest excited electronic states of the alkali hydride cation NaH⁺ have been investigated using an ab initio approach. In this approach, a nonempirical pseudopotential for the Na⁺ core has been used and a core–core and a core‐valence correlation corrections have been added. The adiabatic potential energy curves and the molecular spectroscopic constants for numerous electronic states of ²Σ⁺, ²Π, and ²Δ symmetries, dissociating up to Na (4d) + H⁺ and Na⁺ + H (3d), have been calculated. As no experimental data are available, we discuss our results by comparing with the available theoretical calculations. A satisfying agreement has been found for the ground state with previous works. However, a clear disagreement between this study and the model potential work of Magnier (Magnier, J. Phys. Chem. A 2005, 109, 5411) has been observed for several excited states. Numerous avoided crossings between electronic states of ²Σ⁺ and ²Π symmetries have been found and analysed. They are related to the interaction between the potential energy curves and to the charge transfer process between the two ionic systems Na⁺H and NaH⁺. Furthermore, we provide an extensive set of data concerning the transition dipole moments from X²Σ⁺ and the 2²Σ⁺ states to higher excited states of ²Σ⁺ and ²Π symmetries. Finally, the adiabatic potential energy curves of the ground (X²Σ⁺) and the first (2²Σ⁺) excited states and the transition dipole moments between these states are used to evaluate the radiative lifetimes for the vibrational levels of the 2²∑⁺ state for the first time. In addition to the bound–bound contribution, the bound‐free term has been evaluated and added to the total radiative lifetime. © 2012 Wiley Periodicals, Inc.     

CASPT2 study on low‐lying states of HBN and HNB radicals

In this study, some low‐lying states of the HBN and HNB radicals have been studied using multiconfiguration second‐order perturbation theory. The geometries of all stationary points along the potential energy surfaces (PESs) were optimized at the CASPT2/cc‐pVQZ level. The ground and the first excited states of HBN were predicted to be X²Π and A²Σ⁺ states, respectively. It was predicted that the ground state of HNB is X²Σ⁺ state. The A²Π state of HNB has unique imaginary frequency, which was different from the previously published results. A bending local minimum M1 was found for the first time along the 1²A″ PES, and the A²Π state of HNB should be the transition state of the isomerization reactions for M1 ? M1. The CASPT2/ANO potential energy curves (PECs) of isomerization reactions for HBN ? HNB were calculated as a function of HBN bond angle. By comparing the CASPT2 and CASSCF calculated results, we concluded that the influence of the dynamic electron correlation on HBN ? HNB system is not large. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Transition state theory thermal rate constants and RRKM‐based branching ratios for the N(2D) + CH4 reaction based on multi‐sState and multi‐reference Ab Initio calculations of interest for the titan's chemistry

Multireference single and double configuration interaction (MRCI) calculations including Davidson (+Q) or Pople (+P) corrections have been conducted in this work for the reactants, products, and extrema of the doublet ground state potential energy surface involved in the N( ² D) + CH₄ reaction. Such highly correlated ab initio calculations are then compared with previous PMP4, CCSD(T), W1, and DFT/B3LYP studies. Large relative differences are observed in particular for the transition state in the entrance channel resolving the disagreement between previous ab initio calculations. We confirm the existence of a small but positive potential barrier (3.86 ± 0.84 kJ mol^?1 (MR‐AQCC) and 3.89 kJ mol^?1 (MRCI+P)) in the entrance channel of the title reaction. The correlation is seen to change significantly the energetic position of the two minima and five saddle points of this system together with the dissociation channels but not their relative order. The influence of the electronic correlation into the energetic of the system is clearly demonstrated by the thermal rate constant evaluation and it temperature dependance by means of the transition state theory. Indeed, only MRCI values are able to reproduce the experimental rate constant of the title reaction and its behavior with temperature. Similarly, product branching ratios, evaluated by means of unimolecular RRKM theory, confirm the NH production of Umemoto et al., whereas previous works based on less accurate ab initio calculations failed. We confirm the previous findings that the N( ² D) + CH₄ reaction proceeds via an insertion–dissociation mechanism and that the dominant product channels are CH₂NH + H and CH₃ + NH. © 2012 Wiley Periodicals, Inc.     

Configuration-Interaction study of lower excited states of O2: Valence and rydberg characters of the two lowest 3Σu − states

Configuration-interaction calculations, with an extended basis, are carried out on the ground and lower excited states of O₂ and O₂⁺ at and near the equilibrium internuclear distance (R = 2.3 a.u.) of the ground state of O₂. Particular attention has been paid to the two lowest ³Σ_u^? states, and the mixing of the valence and Rydberg characters in these states are studied. The lowest ³Σ_u^? state is a Rydberg-type state for R < 2.3 a.u., but becomes valence-type for R ? 2.3 a.u. The second ³Σ_u^? state, which is 1.6 eV above the lowest ³Σ_u^? at R = 2.3 a.u., changes its character from Rydberg to valence, valence to Rydberg, and then to valence again when R increases from 1.9 to 3.1 a.u. Satisfactory agreement between the calculated and experimental vertical excitation energies is obtained.     

Theoretical prediction on HAlS+ and HSAl+ cations using multiconfiguration second‐order perturbation theory

Some low‐lying states of the HAlS⁺ and HSAl⁺ cations have been studied for the first time by large‐scale theoretical calculations using three methods: complete active space self‐consistent field (CASSCF), complete active second‐order perturbation theory (CASPT2), and density functional theory Becke's three‐parameter hybrid function with the nonlocal correlation of Lee–Yang–Parr (B3LYP) with the contracted atomic natural orbital (ANO‐L) and cc‐pVTZ basis sets. The geometries of all stationary points along the potential energy surfaces (PESs) were optimized at the CASSCF/ANO‐L and B3LYP/cc‐pVTZ levels. The ground and the first excited states of linear HAlS⁺ are predicted to be X²Π and A²Σ⁺ states, respectively. For the linear HSAl⁺ structure, the first excited state is A²Σ⁺. The X²Π state of linear HSAl⁺ is a second‐order saddle point, because it has two imaginary frequencies. Two bent global minima M1 and M2 were found along the 1²A′ and 1²A″ PESs, respectively. The CASPT2/ANO‐L potential energy curves of isomerization reactions were calculated as a function of HAlS bond angle. According to our calculations, the ground‐state HAlS⁺ is linear, whereas the ground‐state HSAl⁺ is bent. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical prediction on HBeN− and HNBe− anions using multiconfiguration second‐order perturbation theory

The HBeN^? and HNBe^? anions have been investigated for the first time using the CASSCF, CASPT2, and DFT/B3LYP methods with the contracted atomic natural orbital (ANO) and cc‐pVTZ basis sets. The geometries of all stationary points along the potential energy surfaces were optimized at the CASSCF/ANO and B3LYP/cc‐pVTZ levels. The ground and the first excited states of HBeN^? are predicted to be X²Π and A²Σ⁺ states, respectively. It was predicted that the ground state of HNBe^? is X²Σ⁺ state. The A²Π state of HNBe^? has unique imaginary frequency. A bend local minimum M1 was found along the 1²A″ potential energy surface and the A²Π state of HNBe^? should be the transition state of the isomerization reactions for M1 ? M1. The CASPT2/ANO potential energy curves of isomerization reactions were calculated as a function of HBeN bond angle. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

O‐loss photodissociation of the OCS+ ion in the low‐lying electronic states studied using multiconfiguration second‐order perturbation theory

The CASPT2 potential energy curves (PECs) for O‐loss dissociation from the X²Π, A²Π, B²Σ⁺, C²Σ⁺, 1⁴Σ^?, 1²Σ^?, and 1⁴Π states of the OCS⁺ ion were calculated. The PEC calculations indicate that X²Π, 1⁴Σ^?, 1²Σ^?, and 1⁴Π correlate with CS⁺(X²Σ⁺) + O(³P_g); A²Π and B²Σ⁺ correlate with CS⁺(A²Π) + O(³P_g); and C²Σ⁺ probably correlates with CS⁺(X²Σ⁺) + O(¹D_g). The CASSCF minimum energy crossing point (MECP) calculations were performed for the C²Σ⁺/1⁴Σ^?, C²Σ⁺/1⁴Π, A²Π/1⁴Σ^?, A²Π/1²Σ^?, A²Π/1⁴Π, and B²Σ⁺/1²Σ^? state pairs and the spin‐obit couplings were calculated at the located MECPs. A conical intersection point between the B²Σ⁺ and C²Σ⁺ potential energy surfaces was found at the CASSCF level. Based on our calculations, seven O‐loss predissociation processes of the C²Σ⁺ state are suggested and an appearance potential value of 7.13 eV for the CS⁺ + O product group is predicted. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Laser Spectroscopy on the MgH A2Π - X2Σ+ Band System

The A²Π - X²Σ⁺ electronic transition of MgH has been studied by the laser excitation spectroscopy. Some new transitions have been observed for the first time. Rotational parameters of the X and A state have been derived and compared with other experimental values.     

Geometries,vibrational frequencies,and excitation energies of a series of fluorine‐substituted carbenes,FCX (X = H,F, Cl,Br, and I): A high‐level multireference configuration interaction study

High‐level calculations using internally contracted multireference configuration interaction including Davidson correction (icMRCI+Q) method have been carried out for the ground singlet states, the first excited states, and the lowest triplet states of a series of fluorine‐substituted carbenes FCX (X = H, F, Cl, Br, and I). Equilibrium geometries and vibrational frequencies of the three electronic states, adiabatic transition energy of the first excited singlet state, as well as the ground singlet—lowest triplet energy gap (S‐T gap) of each of FCX carbenes have been obtained. Effects of the basis set of icMRCI+Q calculation on the geometries and energies have been investigated. In addition, various corrections, including the scalar relativistic effect, spin‐orbit coupling, and core‐valence correlation, have been studied in calculating the transition energies and the S‐T gaps, especially for heavy‐atom carbenes. This results have been compared with previous calculations using a variety of methods. Our icMRCI+Q results are in very good agreement with the high‐resolution laser‐based spectroscopic results where available. Some structure and spectroscopic constants of the fluorine‐substituted carbenes which are void in the literature have been provided with consistent high‐level calculations. © 2013 Wiley Periodicals, Inc.     

Oscillator strengths and radiative lifetimes for C2: Swan, Ballik-Ramsay, Phillips, and d 3Pig<--c 3Sigmau+ systems

High level ab initio calculations, using multireference configuration interaction (MRCI) techniques, have been carried out to investigate the spectroscopic properties of the singlet A 1Piu<--X 1Sigmag+ Phillips, the triplet d 3Pig<--a 3Sigmau Swan, the b 3Sigmag-<--a 3Piu Ballik-Ramsay, and the d 3Pig<--c 3Sigmau+ transitions of C2. The MRCI expansions are based on full-valence complete active space self-consistent-field reference states and utilize the aug-cc-pV6Z basis set to resolve valence electron correlation. Core and core-valence correlations and scalar relativistic energy corrections were also incorporated in the computed potential energy surfaces. Nonadiabatic and spin-orbit effects were explored and found to be of negligible importance in the calculations. Harmonic frequencies and rotational constants are typically within 0.1% of experiment. The calculated radiative lifetimes compare very well with the available experimental data. Oscillator strengths are reported for all systems: fv'v", where 0相似文献   

High Resolution Laser Excitation Spectra and Franck-Condon Factors of A2Π-X2Σ+ Electronic Transition of MgF

Magnesium monofluoride (MgF) is proposed as an ideal candidate radical for direct laser cooling. Here, the rotationally resolved laser spectra of MgF for the A²Π-X²Σ⁺ electronic transition system were recorded by using laser induced fluorescence technique. The MgF radicals were produced by discharging SF₆/Ar gas mixtures between the tips of two magnesium needles in a supersonic jet expansion. We recorded a total of 19 vibrational bands belonging to three sequences of Δv=0, ±1 in the region of 348-370 nm. Accurate spectroscopic constants for both X²Σ⁺ and A²Π states are determined from rotational analysis of the experimental spectra. Spectroscopic parameters, including the Franck-Condon factors (FCFs), are determined from the experimental results and the Rydberg-Klein-Rees (RKR) calculations. Significant discrepancies between the experimentally measured and RKR-calculated FCFs are found, indicating that the FCFs are nearly independent of the spin-orbit coupling in the A²Π state. Potential energy curves (PECs) and FCFs determined here provide necessary data for the theoretical simulation of the laser-cooling scheme of MgF.     

Accurate ab initio potential energy surfaces for the 3A' and 3A' electronic states of the O(3P)+HBr system

In this work, we report the construction of potential energy surfaces for the (3)A(') and (3)A(') states of the system O((3)P) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O((3)P) + HBr → OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A(') electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A(') surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A(') and 4.16 kcal/mol for the (3)A(') state.     

Global X2A′ potential energy surface of Li2H and quantum dynamics of H + Li2 (X1Σg+) → Li + LiH (X1Σ+) reaction

A global potential energy surface (PES) for the electronic ground state of Li₂H system is constructed over a large configuration space. About 30 000 ab initio energy points have been calculated by MRCI‐F12 method with aug‐cc‐pVTZ basis set. The neural network method is applied to fit the PES and the root mean square error of the current PES is only 1.296 meV. The reaction dynamics of the title reaction has been carried out by employing time‐dependent wave packet approach with second order split operator on the new PES. The reaction probability, integral cross section and thermal rate constant are obtained from the dynamics calculation. In most of the collision energy regions, the integral cross sections are in well agreement with the results reported by Gao et al. The rate constant calculated from the new PES increases in the temperature range of present investigation.     

Nonempirical calculations of the one‐bond 29Si–13C spin–spin coupling constants taking into account relativistic and solvent corrections

The computational study of the one‐bond ²⁹Si–¹³C spin–spin coupling constants has been performed at the second‐order polarization propagator approximation (SOPPA) level in the series of 60 diverse silanes with a special focus on the main factors affecting the accuracy of the calculation including the level of theory, the quality of the basis set, and the contribution of solvent and relativistic effects. Among three SOPPA‐based methods, SOPPA(MP2), SOPPA(CC2), and SOPPA(CCSD), the best result was achieved with SOPPA(CCSD) when used in combination with Sauer's basis set aug‐cc‐pVTZ‐J characterized by the mean absolute error of calculated coupling constants against the experiment of ca 2 Hz in the range of ca 200 Hz. The SOPPA(CCSD)/aug‐cc‐pVTZ‐J method is recommended as the most accurate and effective computational scheme for the calculation of ¹J(Si,C). The slightly less accurate but essentially more economical SOPPA(MP2)/aug‐cc‐pVTZ‐J and/or SOPPA(CC2)/aug‐cc‐pVTZ‐J methods are recommended for larger molecular systems. It was shown that solvent and relativistic corrections do not play a major role in the computation of the total values of ¹J(Si,C); however, taking them into account noticeably improves agreement with the experiment. The rovibrational corrections are estimated to be of about 1 Hz or 1–1.5% of the total value of ¹J(Si,C). Copyright © 2014 John Wiley & Sons, Ltd.