A theoretical study of the excited states of CrH: potential energies, transition moments, and lifetimes

Ab initio calculations of low-lying electronic states of CrH are presented, including potential energies, dipole and transition dipole moment (TDM) functions, and radiative lifetimes for X (6)Sigma(+), A (6)Sigma(+), 3 (6)Sigma(+), 1 (6)Pi, 2 (6)Pi, 3 (6)Pi, and (6)Delta. Calculation of dynamic correlation effects was performed using the multistate complete active space second-order perturbation method, based on state-averaged complete active space self-consistent-field reference wave functions obtained with seven active electrons in an active space of 16 molecular orbitals. A relativistic atomic natural orbital-type basis set from the MOLCAS library was used for Cr. Good agreement is found between the current calculations and experiment for the lowest two (6)Sigma(+) states, the only states for which spectroscopic data are available. Potential curves for the 3 (6)Sigma(+) and 2 (6)Pi states are complicated by avoided crossings with higher states of the same symmetry, thus resulting in double-well structures for these two states. The measured bandhead T(0)=27 181 cm(-1), previously assigned to a (6)Pi<--X (6)Sigma(+) transition, is close to our value of T(0)=28 434 cm(-1) for the 2 (6)Pi state. We tentatively assign the ultraviolet band found experimentally at 30 386 cm(-1) to the 3 (6)Pi<--X (6)Sigma(+) transition for which the computed value is 29 660 cm(-1). The A (6)Sigma(+)<--X (6)Sigma(+) TDM and A (6)Sigma(+) lifetimes are found to be in reasonable agreement with previous calculations.     

Coupled-cluster with active space selected higher amplitudes: performance of seminatural orbitals for ground and excited state calculations

The active space approach for coupled-cluster models is generalized using the general active space concept and implemented in a string-based general coupled-cluster code. Particular attention is devoted to the choice of orbitals on which the subspace division is based. Seminatural orbitals are proposed for that purpose. These orbitals are obtained by diagonalizing only the hole-hole and particle-particle block of the one-electron density of a lower-order method. The seminatural orbitals are shown to be a good replacement for complete active space self-consistent field orbitals and avoid the ambiguities with respect to the reference determinant introduced by the latter orbitals. The seminatural orbitals also perform well in excited state calculations, including excited states with strong double excitation contributions, which usually are difficult to describe with standard coupled-cluster methods. A set of vertical excitation energies is obtained and benchmarked against full configuration interaction calculations, and alternative hierarchies of active space coupled-cluster models are proposed. As a simple application the spectroscopic constants of the C(2) B (1)Delta(g) and B(') (1)Sigma(g) (+) states are calculated using active space coupled-cluster methods and basis sets up to quadruple-zeta quality in connection with extrapolation and additivity schemes.     

Ab initio study of the low lying electronic states of ZnF and ZnF-

Highly correlated ab initio calculations have been performed for an accurate determination of the electronic structure and of the spectroscopy of the low lying electronic states of the ZnF system. Using effective core pseudopotentials and aug-cc-pVQZ basis sets for both atoms, the potential curves, the dipole moment functions, and the transition dipole moments between relevant electronic states have been calculated at the multireference-configuration-interaction level. The spectroscopic constants calculated for the X(2)Sigma(+) ground state are in good agreement with the most recent theoretical and experimental values. It is shown that, besides the X(2)Sigma(+) ground state, the B(2)Sigma(+), the C(2)Pi, and the D(2)Sigma(+) states are bound. The A(2)Pi state, which has been mentioned in previous works, is not bound but its potential presents a shoulder in the Franck-Condon region of the X(2)Sigma(+) ground state. All of the low lying quartet states are found to be repulsive. The absorption transitions from the v=0 level of the X(2)Sigma(+) ground state toward the three bound states have been evaluated and the spectra are presented. The potential energy of the ZnF(-) molecular anion has been determined in the vicinity of its equilibrium geometry and the electronic affinity of ZnF (EA=1.843 eV with the zero energy point correction) has been calculated in agreement with the photoelectron spectroscopy experiments.     

Theoretical investigation of the SO(2+) dication and the photo-double ionization spectrum of SO

Highly correlated ab initio methods were used in order to generate the potential energy curves of the electronic states of the SO(2+) dication and of the electronic ground state of the neutral SO molecule. These curves were used to predict the spectroscopic properties of this dication and to perform forward calculations of the double photoionization spectrum of SO. In light of spin-orbit calculations, the metastability of this doubly charged ion is discussed: for instance, the rovibrational levels of the X (1)Sigma(+) and A (3)Sigma(+) states are found to present relatively long lifetimes. In contrast, the other electronic excited states should predissociate to form S(+) and O(+) in their electronic ground states. The simulated spectrum shows structures due to transitions between the v=0 vibrational level of SO (X (3)Sigma(-)) and the vibrational levels below the barrier maximum of 11 of the calculated electronic states. The 2 (1)Sigma(+) electronic state of SO(2+) received further treatment: in addition to vibrational bands due to the below barrier energy levels of this electronic state, at least nine continuum resonances were predicted which are responsible for the special shape of the spectrum in this energy region. This work is predictive in nature and should stimulate future experimental investigations dealing with this dication.     

A theoretical study of the fine and hyperfine interactions in the NCO and CNO radicals

The geometries, the harmonic vibrational frequencies, and the Renner-Teller parameter have been reported for the NCO(+)(X (3)Sigma(-)), NCO(X (2)Pi,A (2)Sigma(+),B (2)Pi,2 (2)Sigma(+)), NCO(-)(X (1)Sigma(+)), CNO(+)(X), CNO(X (2)Pi,A (2)Sigma(+),B (2)Pi,2 (2)Sigma(+)), and CNO(-)(X (1)Sigma(+)) systems at the full valence-complete active space self-consistent-field (fv-CASSCF) level of theory. The (2)Pi electronic states of the NCO and CNO radicals have two distinct real vibrational frequencies for the bending modes and these states are subject to the type A Renner-Teller effect. The total energy of CNO(+) without zero point energy correction of the linear geometry is approximately 31 cm(-1) higher than the bent geometry at the fv-CASSCF level and the inversion barrier vanishes after the zero point energy correction; therefore, the ground state of the CNO(+) may possess a quasilinear geometry. The spin-orbit coupling constants estimated using atomic mean field Hamiltonian at the fv-CASSCF level of theory are in better agreement with the experimental values. The excitation energies, the electron affinity, and the ionization potential have been computed at the complete active space second order perturbation theory (CASPT2) and the multireference singles and doubles configuration (MRSD-CI) levels of theory. The computed values of the electric hyperfine coupling constants for the (14)N atom in the ground state of the NCO radical agree well with the experimental data. The magnetic hyperfine coupling constants (HFCC's) have been estimated employing the configuration selected MRSD-CI and the multireference singles configuration interaction (MRS-CI) methods using iterative natural orbitals (ino) as one particle basis. Sufficiently accurate value of the isotropic contribution to the HFCC's can be obtained using an MRS-CI-ino procedure.     

Extremely narrow peaks in predissociation of sodium dimer due to rovibronic coupling

In sodium dimer the 2 (3)Pi(g), 3 (3)Pi(g), and 4 (3)Sigma(g) (+) electronic states are coupled; the coupling of the two (3)Pi(g) states is due to vibrational motion while the nonadiabatic interaction between the (3)Sigma(g) (+) and the (3)Pi(g) states-in particular, the 3 (3)Pi(g) state-is mediated by rotational interaction. The resulting vibronic problem is studied in some detail. The bound vibrational states of the 3 (3)Pi(g) and 4 (3)Pi(g) (+) states lie in the dissociation continuum of the 2 (3)Pi(g) state and become resonances due to the prevailing nonadiabatic coupling. The resonances are calculated using the complex scaling method and the available ab initio adiabatic potential energy curves. It is demonstrated that the resonances associated with rotational nonadiabatic coupling are narrower by several orders of magnitude than those that emerge from the vibrational nonadiabatic coupling. The predissociation cross section is computed and compared with experiment.     

Quantum Monte Carlo study of the ground state and low-lying excited states of the scandium dimer

A large set of electronic states of scandium dimer has been calculated using high-level theoretical methods such as quantum diffusion Monte Carlo (DMC), complete active space perturbation theory as implemented in GAMESS-US, coupled-cluster singles, doubles, and triples, and density functional theory (DFT). The 3 Sigma u and 5 Sigma u states are calculated to be close in energy in all cases, but whereas DFT predicts the 5 Sigma u state to be the ground state by 0.08 eV, DMC and CASPT2 calculations predict the 3 Sigma u to be more stable by 0.17 and 0.16 eV, respectively. The experimental data available are in agreement with the calculated frequencies and dissociation energies of both states, and therefore we conclude that the correct ground state of scandium dimer is the 3 Sigma u state, which breaks with the assumption of a 5 Sigma u ground state for scandium dimer, believed throughout the past decades.     

Ab initio potential-energy surfaces of O2(X3Sigmag -, a1Deltag, b1Sigmag +) +O2 (X3Sigmag -, a1Deltag, b1Sigmag +): mechanism of quenching of O2 (a 1Deltag)

Ab initio computational studies were carried out in order to explore the possible mechanisms of quenching of O(2)(a (1)Delta(g)) by O(2)(X (3)Sigma(g) (-)): the self-quenching of O(2)(a (1)Delta(g)) and other energy-transfer processes involving two O(2) molecules. All eighteen states arising from two O(2) molecules in the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) states are considered. After scans at the state-averaged complete active space self-consistent field method to identify possible regions of crossing between states belonging to different asymptotes, complete active state second-order perturbation theory high-symmetry optimization and low-symmetry scans established that four different minima on the seams of crossing (MSXs), arising between the a+a manifold and the X+b manifold and responsible for self-quenching: O(2)(a (1)Delta(g))+O(2)(a (1)Delta(g))-->O(2)(X (3)Sigma(g) (-))+O(2)(b (1)Sigma(g) (+)), have coplanar C(2h) or C(2v) symmetries and are only 0.45 eV barrier relative to the a+a asymptote and energetically easily accessible. The rate constant for this process was estimated based on the Landau-Zener formalism. The MSXs for quenching of O(2)(a (1)Delta(g)) by the ground state O(2)(X (3)Sigma(g) (-)):O(2)(a (1)Delta(g))+O(2)(X (3)Sigma(g) (-))-->O(2)(X (3)Sigma(g) (-))+O(2)(X (3)Sigma(g) (-)) require higher energies and the process is not likely to be important.     

Electronic states and spectroscopic properties of SiTe and SiTe+

Ab initio based configuration interaction calculations have been carried out to study the low-lying electronic states and spectroscopic properties of the heaviest nonradioactive silicon chalcogenide molecule and its monopositive ion. Spectroscopic constants and potential energy curves of states of both SiTe and SiTe+ within 5 eV are reported. The calculated dissociation energies of SiTe and SiTe+ are 4.41 and 3.52 eV, respectively. Effects of the spin-orbit coupling on the electronic spectrum of both the species are studied in detail. The spin-orbit splitting between the two components of the ground state of SiTe+ is estimated to be 1880 cm(-1). Transitions such as 0+ (II)-X1Sigma(+)0+, 0+ (III)-X1Sigma(+)0+, E1Sigma(+)0+ -X1Sigma(+)0+, and A1Pi1-X1Sigma(+)0+ are predicted to be strong in SiTe. The radiative lifetime of the A1Pi state is less than a microsecond. The X(2)2Pi(1/2)-X(1)2Pi(3/2) transition in SiTe+ is allowed due to spin-orbit mixing. However, it is weak in intensity with a partial lifetime for the X2 state of about 108 ms. The electric dipole moments of both SiTe and SiTe+ in their low-lying states are calculated. The vertical ionization energies for the ionization of the ground-state SiTe to different ionic states are also reported.     

Ab initio configuration interaction study of the low-lying 1Sigma+ electronic states of LiCl

Ab initio configuration interaction calculations have been performed for the X 1Sigma+ and B 1Sigma+ electronic states of LiCl. Potential energy curves, dipole moment functions, and dipole transition moments have been computed for internuclear distances between R = 2.5a0 and 50a0. Single- and double-excitation configuration interaction wave functions were constructed using molecular orbitals obtained from a two-state averaged multiconfiguration self-consistent-field calculation. This procedure yielded an accurate energy splitting between the covalent and ionic separated-atom limits. The calculated avoided crossing of the X and B state curves occurs at R = 16.2a0, in close agreement with previous calculations using a semiempirical covalent-ionic resonance model. X 1Sigma+ state spectroscopic constants are in excellent agreement with experimental values.     

Observation of L uncoupling in the 5 1Deltag Rydberg state of Na2

The phenomenon of electronic orbital angular momentum L uncoupled from its internuclear axis has been observed in the sodium dimer using high-resolution cw optical-optical double-resonance spectroscopy. When L uncoupling occurs, the degeneracy of Lambda doubling is removed. In our experiment, the intermediate B (1)Pi(u) state of Na(2) is excited from the thermally populated ground X (1)Sigma(g) (+) state by a single-line Ar(+) laser. Then, a single-mode dye laser is used to probe the Rydberg states from the intermediate state. The signals are detected by monitoring the UV fluorescence from the triplet gerade states back to the a (3)Sigma(u) (+) state via collision energy transfer. Under our experimental resolution, the splitting of Lambda doubling in the 5 (1)Delta(g) state of Na(2) can be measured. A total of 136 rovibronic levels with ef parities have been assigned to the 5 (1)Delta(g) state. The Lambda-splitting constants deduced from these data are q(0)=0.376(90)x10(-4) cm(-1), q(v)=0.114(6)x10(-4) cm(-1), and mu=0.76(33)x10(-8) cm(-1). In general, the Lambda splitting of the Delta states is considerably smaller than that of the Pi states. However, the first-order splitting constants q(0) and q(v) reported here are larger than those in the B (1)Pi(u) state. This is due to the L uncoupling of the Rydberg states.     

Vacuum ultraviolet photoionization of C3

Photoionization efficiency (PIE) curves for C(3) molecules produced by laser ablation are measured from 11.0 to 13.5 eV with tunable vacuum ultraviolet undulator radiation. A step in the PIE curve versus photon energy, obtained with N(2) as the carrier gas, supports the conclusion of very effective cooling of C(3) to its linear (1)Sigma(g)(+) ground state. The second step observed in the PIE curve versus photon energy could be the first experimental evidence of the C(3)(+)((2)Sigma(g)(+)) excited state. The experimental results, complemented by ab initio calculations, suggest a state-to-state vertical ionization energy of 11.70 +/- 0.05 eV between the C(3)(X(1)Sigma(g)(+)) and the C(3)(+)(X(2)Sigma(u)(+)) states. An ionization energy of 11.61 +/- 0.07 eV between the neutral and ionic ground states of C(3) is deduced using the data together with our calculations. Accurate ab initio calculations are performed for both linear and bent geometries on the lowest doublet electronic states of C(3)(+) using Configuration Interaction (CI) approaches and large basis sets. These calculations confirm that C(3)(+) is bent in its electronic ground state, which is separated by a small potential barrier from the (2)Sigma(u)(+) minimum. The gradual increase at the onset of the PIE curve suggests a geometry change between the ground neutral and cationic states. The energies between several doublet states of the ion are theoretically determined to be 0.81, 1.49, and 1.98 eV between the (2)Sigma(u)(+) and the (2)Sigma(g)(+),( 2)Pi(u), (2)Pi(g) excited states of C(3)(+), respectively.     

An ab initio multireference perturbation theory study on the manganese dimer

The potential energy curves of the ground state and of some excited states of the manganese dimer have been calculated over a wide range of internuclear distances using the second order n-electron valence state perturbation theory applied to a complete active space self-consistent field reference wave function. The ground state of Mn(2), for which also the third order NEVPT has been used, is calculated to be a singlet belonging to the Sigma(g) (+) symmetry, characterized by a large equilibrium internuclear distance R(e) of 3.7-3.8 A, by a low dissociation energy D(e) of 0.07-0.08 eV, and by a small harmonic frequency omega(e) of 43 cm(-1). The experimental evidence that Mn(2) is a van der Waals molecule is thus confirmed. Among the excited states, (11)Pi(u), which is usually indicated as the ground state by density functional theory studies, appears as a low-lying state with R(e)=2.50 A, D(e)=1.35 eV, and omega(e)=246 cm(-1).     

Broken-symmetry unrestricted hybrid density functional calculations on nickel dimer and nickel hydride

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unacceptably large error, underestimating the dissociation energy of Ni(2) by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions. We also find that for Ni(2), the spin projection for the broken-symmetry unrestricted singlet states changes the ordering of the states, but the splittings are less than 10 meV. All our calculations predict a deltadelta-hole ground state for Ni(2) and delta-hole ground state for NiH. Upon spin projection of the singlet state of Ni(2), almost all of our calculations: Becke98 and FSLYP both AE and ECP and B3LYP/AE predict (1)(d(A)(x(2)-y(2)d(B)(x(2)-y(2)) or (1)(d(A)(xy) (d)(B)(xy)) ground state, which is a mixture of (1)Sigma(g) (+) and (1)Gamma(g). B3LYP/ECP predicts a (3)(d(A)(x(2)-y(2))d(B)(xy) (mixture of (3)Sigma(g) (-) and (3)Gamma(u)) ground state virtually degenerate with the (1)(d(A)(x(2)-y(2)d(B)(x)(2)-y(2)/(1)(d(A)(xy)D(B)(xy) state. The doublet delta-hole ground state of NiH predicted by all our calculations is in agreement with the experimentally predicted (2)Delta ground state. For Ni(2), all our results are consistent with the experimentally predicted ground state of 0(g) (+) (a mixture of (1)Sigma(g) (+) and (3)Sigma(g) (-)) or 0(u) (-) (a mixture of (1)Sigma(u) (-) and (3)Sigma(u) (+)).     

The ionization energy of Be2, and spectroscopic characterization of the (1)3Sigmau+, (2)3Pig, and (3)3Pig states

Low lying electronic states of the beryllium dimer were investigated by laser induced fluorescence (LIF) and resonance enhanced multiphoton ionization (REMPI) techniques. Be(2) was formed by pulsed laser ablation of Be metal in the presence of helium carrier gas, followed by a free jet expansion into vacuum. Several previously unobserved states of the dimer were characterized. These included transitions of the triplet manifold (2)(3)Pi(g) <-- (1)(3)Sigma(u)+ and (3)(3)Pi(g) <-- (1)(3)Sigma(u)+, for which rotationally resolved bands were obtained. In addition, transitions to the v' = 10-18 vibrational levels of the A (1)Pi(u) state were recorded. Photoionization efficiency (PIE) measurements were used to determine an accurate ionization energy (IE) for Be(2) of 7.418(5) eV and the term energy for (1)(3)Sigma(u)+. Above the ionization threshold the PIE spectrum was found to be highly structured, consisting of overlapping Rydberg series that converged on excited vibrational levels of Be(2)+. Analysis of these series yielded a vibration frequency for the X(2)Sigma(u)+ state of 498(20) cm(-1). The bond dissociation energy for Be(2)+, deduced from the IE measurement, was 16 072(40) cm(-1). Multi-reference configuration interaction (MRCI) calculations were carried out for Be(2) and Be(2)+, yielding results that were in excellent agreement with the experimental observations.     

Full configuration interaction calculation of BeH adiabatic states

An all-electron full configuration interaction (FCI) calculation of the adiabatic potential energy curves of some of the lower states of BeH molecule is presented. A moderately large ANO basis set of atomic natural orbitals (ANO) augmented with Rydberg functions has been used in order to describe the valence and Rydberg states and their interactions. The Rydberg set of ANOs has been placed on the Be at all bond distances. So, the basis set can be described as 4s3p2d1f3s2p1d(BeH)+4s4p2d(Be). The dipole moments of several states and transition dipole strengths from the ground state are also reported as a function of the R(Be-H) distance. The position and the number of states involved in several avoided crossings present in this system have been discussed. Spectroscopic parameters have been calculated from a number of the vibrational states that result from the adiabatic curves except for some states in which this would be completely nonsense, as it is the case for the very distorted curves of the 3s and 3p (2)Sigma(+) states or the double-well potential of the 4p (2)Pi state. The so-called "D complex" at 54 050 cm(-1) (185.0 nm) is resolved into the three 3d substates ((2)Sigma(+),(2)Pi,(2)Delta). A diexcited valence state is calculated as the lowest state of (2)Sigma(-) symmetry and its spectroscopic parameters are reported, as well as those of the 2 (2)Delta (4d) state The adiabatic curve of the 4 (2)Sigma(+) state shows a swallow well at large distances (around 4.1 A) as a result of an avoided crossing with the 3 (2)Sigma(+) state. The probability that some vibrational levels of this well could be populated is discussed within an approached Landau-Zerner model and is found to be high. No evidence is found of the E(4ssigma) (2)Sigma(+) state in the region of the "D complex". Instead, the spectroscopic properties obtained from the (4ssigma) 6 (2)Sigma(+) adiabatic curve of the present work seem to agree with those of the experimental F(4psigma) (2)Sigma(+) state. The FCI calculations provide benchmark results for other correlation models for the open-shell BeH system and evidence both the limitations and capabilities of the basis set.     

High resolution vacuum ultraviolet emission spectrum of D(2) from 78 to 103 nm: The D (1)Pi(u)-->X (1)Sigma(g) (+) and D(') (1)Pi(u) (-)-->X (1)Sigma(g) (+) band systems

The emission spectrum of the D(2) molecule has been studied at high resolution in the vacuum ultraviolet region 78.5-102.7 nm. A detailed analysis of the two D (1)Pi(u)-->X (1)Sigma(g) (+) and D(') (1)Pi(u) (-)-->X (1)Sigma(g) (+) electronic band systems is reported. New and improved values of the level energies of the two upper states have been derived with the help of the program IDEN [V. I. Azarov, Phys. Scr. 44, 528 (1991); 48, 656 (1993)], originally developed for atomic spectral analysis. A detailed comparison is made between the observed energy levels and solutions of coupled equations using the newest ab initio potentials by Wolniewicz and co-workers [J. Chem. Phys. 103, 1792 (1995); 99, 1851 (1993); J. Mol. Spectros. 212, 208 (2002); 220, 45 (2003)] taking into account the nonadiabatic coupling terms for the D (1)Pi(u) state with the lowest electronic states B (1)Sigma(u) (+), C (1)Pi(u), and B(') (1)Sigma(u) (+). A satisfactory agreement has been found for most of the level energies belonging to the D and D(') states. The remaining differences between observation and theory are probably due to nonadiabatic couplings with other higher electronic states which were neglected in the calculations.