Ab initio investigation of titanium hydroxide isomers and their cations, TiOH(0,+) and HTiO(0,+)

We studied the electronic and geometrical structure of the [Ti, O, H](0,+) species, using large basis sets and both single-reference coupled cluster and multireference configuration interaction methodologies. The electronic structure of HTiO(0,+) is interpreted qualitatively in terms of a hydrogen atom bonding to TiO(0,+), while the structure of TiOH(0,+) is interpreted in terms of Ti(+,2+) bonding to OH(-). Potential energy profiles are reported as functions of the Ti-OH and H-TiO bond lengths, and of the H-Ti-O angle. For a total of 33 stationary points on the potential energy surfaces, we report absolute energies, geometries, and harmonic frequencies. For the neutral species, dipole moments are also given.     

A density functional study of the structures and ionization energies of some scandium compounds with hydrogen and oxygen

The geometries, frequencies and energies of various Sc species with hydrogen and oxygen which might exist in a flame doped with Sc were studied with density functional methods. Two functionals, B3LYP and BP86, equipped with two basis sets, TZVP and 6-311++G(3df,3pd), were employed. In predicting geometries and frequencies, the BP86 functional performs consistently better than does the B3LYP functional, but the B3LYP functional performs better for energy calculations. For most species, the effects of using different basis-set effects are minor. The calculated ionization energies of 6.47 eV for Sc and 6.6 eV for ScO are in good agreement with the experimental values of 6.561 for Sc and 6.6±0.3 eV for ScO. For other species, suggested values of the calculated ionization energies are 6.3 eV for ScH, 6.1 for ScOH, 9.3 for HScO, 6.5 for HScOH, 9.5 for OScOH, 6.2 for Sc(OH)₂, 9.6 for HSc(OH)₂ and 10.1 eV for Sc(OH)₃, all with an approximate error of 0.2–0.3 eV. These high values indicate that the formation of Sc-containing ions in flames by thermal (collisional) ionization is unlikely.     

Ab initio and quantum-defect calculations for the Rydberg states of ArH

Potential energy curves were evaluated for the ground and thirteen low-lying excited electronic states of the ArH molecule over a wide range of internuclear distances by the multi-reference averaged quadratic coupled cluster method. The ab initio energy differences and transition dipole moments were used to estimate Einstein emission coefficients, absorption oscillator strengths and radiative lifetimes. Diagonal and off-diagonal quantum defects, as functions of internuclear distance, were extracted from ab initio potentials of the lowest Rydberg states of the neutral ArH molecule by taking account of configuration interaction between Rydberg series converging to the ground and two electronic excited states of the ArH(+) cation. The derived quantum-defect functions were used to generate manifolds of higher excited Rydberg states. The agreement between experimental and calculated energies and radiative transition probabilities was found to be as good as or better than that obtained by earlier calculations.     

Theoretical investigation of excited states of C(3)

In this work, we present ab initio calculations for the potential energy surfaces of C(3) in different electronic configurations, including the singlet ground state [X (1)Sigma(g) (+),((1)A(1))], the triplet ground state [a (3)Pi(u),((3)B(1), (3)A(1))], and some higher excited states. The geometries studied include triangular shapes with two identical bond lengths, but different bond angles between them. For the singlet and triplet ground states in the linear geometry, the total energies resulting from the mixed density functional--Hartree-Fock and quadratic configuration interaction methods reproduce the experimental values, i.e., the triplet occurs 2.1 eV above the singlet. In the geometry of an equilateral triangle, we find a low-lying triplet state with an energy of only 0.8 eV above the energy of the singlet in the linear configuration, so that the triangular geometry yields the lowest excited state of C(3). For the higher excited states up to about 8 eV above the ground state, we apply time-dependent density functional theory. Even though the systematic error produced by this approach is of the order of 0.4 eV, the results give different prospective to insight into the potential energy landscape for higher excitation energies.     

Interaction of the early 3d transition metals Sc, Ti, V, and Cr with N2: an ab initio study

The interaction of the early 3d transition elements M=Sc, Ti, V, and Cr with N2(X 1Sigmag+) has been studied by coupled-cluster and multiconfigurational techniques in conjunction with quantitative basis sets. We investigated both triatomic (MN2) and tetratomic (M2N2) species but focused mainly on high-spin linear and T-shaped triatomics. The lowest bound states of ScN2(4B1),TiN2(5Delta), and VN2(6Sigma+) correlate to the first excited state of the M atom, with M-N2 binding energies (De) of 24, 14, and 8 kcal/mol, respectively. In CrN2, the first bound state (7) product operator correlates to the sixth excited state of the Cr atom (7P) with De = 27 kcal/mol. The M-N2-M bond strength of high-spin linear tetratomics is twice as large the binding energy of the corresponding M-N2 linear triatomics, M = Sc, Ti, V, and Cr.     

The electronic spectrum of protonated adenine: theory and experiment

In this work we present the results of a combined experimental and theoretical study concerned with the question how a proton changes the electronic spectrum and dynamics of adenine. In the experimental part, isolated adenine ions have been formed by electro-spray ionisation, stored, mass-selected and cooled in a Paul trap and dissociated by resonant photoexcitation with ns UV laser pulses. The S(0)-S1 spectrum of protonated adenine recorded by fragment ion detection lies in a similar energy range as the first pipi* transition of neutral 9H-adenine. It shows a flat onset with a broad substructure, indicating a large S(0)-S1 geometry shift and an ultra-short lifetime. In the theoretical part, relative energies of the ground and the excited states of the most important tautomers have been calculated by means of a combined density functional theory and multi-reference configuration interaction approach. Protonation at the nitrogen in position 1 of the neutral 9H-adenine tautomer yields the most stable protonated adenine species, 1H-9H-A+. The 3H-7H-A+ and the 3H-9H-A+ tautomers, formed by protonation of 7H- and 9H-adenine in 3-position, are higher in energy by 162 cm(-1) and 688 cm(-1), respectively. Other tautomers lie at considerably higher energies. Calculated vertical absorption spectra are reported for all investigated tautomers whereas geometry optimisations of excited states have been carried out only for the most interesting ones. The S1 state energies and geometries are found to depend on the protonation site. The theoretical data match best with the experimental onset of the spectrum for the 1H-9H-A+ tautomer although we cannot definitely exclude contributions to the experimental spectrum from the 3H-7H-A+ tautomer at higher energies. The vertical S(0)--> S1 excitation energy is similar to the one in neutral 9H-adenine. As for the neutral adenine, we find a conical intersection of the S1 of protonated adenine with the ground state in an out-of-plane coordinate but at lower energies and accessible without barrier.     

Electronic structure of 2,2'-bipyridine organotransition-metal complexes. Establishing the ligand oxidation level by density functional theoretical calculations

A density functional theoretical (DFT) study (B3LYP) has been carried out on 20 organometallic complexes containing η(5)- and/or η(3)-coordinated cyclopentadienyl anions (Cp(-)) and 2,2'-bipyridine (bpy) ligand(s) at varying oxidation levels, i.e., as the neutral ligand (bpy(0)), as the π-radical monoanion (bpy(?-))(-), or as the diamagnetic dianion (bpy(2-))(2-). The molecular and electronic structures of these species in their ground states and, in some cases, their first excited states have been calculated using broken-symmetry methodology. The results are compared with experimental structural and spectroscopic data (where available) in order to validate the DFT computational approach. The following electron-transfer series and complexes have been studied: [(Cp)(2)V(bpy)](0,+,2+) (1-3), [(Cp)(2)Ti(bpy)](-,0,+,2+) (4-7), [(Cp)(2)Ti(biquinoline)](0,+) (8 and 9), [(Cp*)(2)Ti(bpy)](0) (10) (Cp* = pentamethylcyclopentadienyl anion), [Cp*Co(bpy)](0,+) (11 and 12), [Cp*Co(bpy)Cl](+,0) (13 and 14), [Fe(toluene)(bpy)](0) (15), [Cp*Ru(bpy)](-) (16), [(Cp)(2)Zr(bpy)](0) (17), and [Mn(CO)(3)(bpy)](-) (18). In order to test the predictive power of our computations, we have also calculated the molecular and electronic structures of two complexes, A and B, namely, the diamagnetic dimer [Cp*Sc(bpy)(μ-Cl)](2) (A) and the paramagnetic (at 25 °C) mononuclear species [(η(5)-C(5)H(4)(CH(2))(2)N(CH(3))(2))Sc((m)bpy)(2)] (B). The crystallographically observed intramolecular π-π interaction of two N,N'-coordinated π-radical anions in A leading to an S = 0 ground state is reliably reproduced. Similarly, the small singlet-triplet gap of ~600 cm(-1) between two antiferromagnetically coupled (bpy(?-))(-) ligands in B, two ferromagnetically coupled radical anions in the triplet excited state of B, and the structures of A and B is reproduced. Therefore, we are confident that we can present computationally obtained, detailed electronic structures for complexes 1-18. We show that N,N'-coordinated neutral bpy(0) ligands behave as very weak π acceptors (if at all), whereas the (bpy(2-))(2-) dianions are strong π-donor ligands.     

Theoretical study of the electronic structure of MCH(2)(+)(M=Fe,Co,Ni)

State of the art coupled cluster (CC) methods are applied to accurately characterize the ground state electronic structure and photoelectron spectra of transition metal carbene ions MCH(2) (+) (M=Fe, Co, and Ni). The geometries and energies of the lowest energy quartet, triplet, and doublet electronic states as well as several low-lying vertical excitation energies of FeCH(2) (+), CoCH(2) (+), and NiCH(2) (+) are reported. The excitation energies are computed using the equation-of-motion CC and for states of different symmetries, by the energy differences of single reference ground and excited states (Delta-CC). The latter employ several reference states; the unrestricted Hartree-Fock, restricted open shell Hartree-Fock, and unrestricted Kohn-Sham. We conclude that the (2)A(1) electronic ground state of NiCH(2) (+) is separated by about 30.0 kJ/mol from the next highest state, and the lowest (4)B(1) and (4)B(2) states of FeCH(2) (+) as well as the (3)A(2) and (3)A(1) states of CoCH(2) (+) are nearly degenerate. The presence of metal-pi*(MCH(2)) charge transfer states with significant oscillator strengths in the visible/near-UV energy domain of the theoretical spectra of FeCH(2) (+) and CoCH(2) (+) are at the origin of the photofragmentation of these compounds observed after irradiation between 310 and 360 nm.     

The photodissociation of NO(2) by visible and ultraviolet light

We present velocity map images of the NO, O((3)P(J)) and O((1)S(0)) photofragments from NO(2) excited in the range 7.6 to 9.0 eV. The molecule was initially pumped with a visible photon between 2.82-2.95 eV (440-420 nm), below the first dissociation threshold. A second ultraviolet laser with photon energies between 4.77 and 6.05 eV (260-205 nm) was used to pump high-lying excited states of neutral NO(2) and/or probe neutral photoproducts. Analysis of the kinetic energy release spectra revealed that the NO photofragments were predominantly formed in their ground electronic state with little kinetic energy. The O((3)P(J)) and O((1)S(0)) kinetic energy distributions were also dominated by kinetically 'cold' fragments. We discuss the possible excitation schemes and conclude that the unstable photoexcited states probed in the experiment were Rydberg states coupled to dissociative valence states. We compare our results with recent time-resolved studies using similar excitation and probe photon energies.     

Accurate theoretical study of PS(q) (q = 0,+1,-1) in the gas phase

Highly correlated ab initio methods were used in order to generate the potential energy curves and spin-orbit couplings of electronic ground and excited states of PS and PS(+). We also computed those of the bound parts of the electronic states of the PS(-) anion. We used standard coupled cluster CCSD(T) level with augmented correlation-consistent basis sets, internally contacted multi-reference configuration interaction, and the newly developed CCSD(T)-F12 methods in connection with the explicitly correlated basis sets. Core-valence correction and scalar relativistic effects were examined. Our data consist of a set of spectroscopic parameters (equilibrium geometries, harmonic vibrational frequencies, rotational constants, spin-orbit, and spin-spin constants), adiabatic ionization energies, and electron affinities. For the low laying electronic states, our calculations are consistent with previous works whereas the high excited states present rather different shapes. Based on these new computations, the earlier ultraviolet bands of PS and PS(+) were reassigned. For PS(-) and in addition to the already known anionic three bound electronic states (i.e., X(3)Σ(-), (1)Δ, and 1(1)Σ(+)), our calculations show that the (1)Σ(-), (3)Σ(+), and the (3)Δ states are energetically below their quartet parent neutral state (a(4)Π). The depletion of the J = 3 component of PS(-)((3)Δ) will mainly occur via weak interactions with the electron continuum wave.     

Theoretical studies of the excited doublet states of SF and SCl and singlet states of SF2, SFCl, and SCl2

In previous work, we reported that the lowest-lying excited states of SF, SCl, SF(2), SFCl, and SCl(2) have recoupled pair bonds. In this study, we examine the analogous low-spin states--the (2)Σ(-) and (2)Δ states of SF and SCl and the excited singlet states of SF(2), SCl(2), and SFCl--which also possess recoupled pair bonds. In contrast to the excited states treated previously, the states studied in the present work have the same spin multiplicities as their respective ground states and are thus potentially observable via electronic excitation. Of particular interest are the minima on the (1)A″ potential energy surface of SFCl corresponding to bond-stretch isomers analogous to those found on the (3)A″ surface. In addition, we discovered that the first two excited states ((1)A″) accessible via vertical excitations from the ground state of SFCl have the electronic structure of the bond-stretch isomers. Thus, electronic excitation spectroscopic studies of SFCl could reveal a signature of the bond-stretch isomers. We will also present limited data on the lowest singlet Rydberg states of the triatomic species. Calculations were performed at the MRCI+Q/aug-cc-pV(Q+d,5+d)Z levels of theory.     

Structure and dynamics of the aniline-argon complex as derived from its potential energy surface

The structure and dynamics of the van der Waals (vdW) complex of aniline (An) with argon (Ar) are studied using ab initio methods. The inversion potential of the aniline-argon (AnAr) complex perturbed by the weak vdW interaction is calculated taking into account subtle corrections from the zero-point energy of the vdW modes and from the frequency shifts of the An normal modes modified by the complexation. The intermolecular potential energy surface (PES) of the AnAr complex is determined by performing a large-scale computation of the interaction energy and the fitting of the analytical many-body expansion to the set of single-point interaction energies. The PES determined shows two deep local minima corresponding to the anti and syn AnAr conformers. The difference in the energies of these two minima is only 15 cm-1, but it is sufficient to localize the inversion wave functions and to form the two conformers. In the conformer anti (syn) of lower (higher) energy, Ar is bound to the An ring opposite (adjacent) the amino-hydrogens. In the additional local minima higher in energy, Ar approaches the aniline ring between the C-H bonds near its plane. An additional local minimum is located opposite of nitrogen between the two N-H bonds. The high-energy minima are, however, too flat to form stable conformers. The perturbation of the interaction of Ar with the phenyl ring by the NH2 group is described by the vdW hole, which is responsible for unusually strong intermode mixing in the excited intermolecular vibrational states. The analysis of these states calculated for the ground (S0) as well as the first excited electronic state (S1) resolves difficulties faced earlier with the assignment of the observed vibronic bands of AnAr.     

Spectroscopic properties of novel aromatic metal clusters: NaM4 (M=Al,Ga,In) and their cations and anions

The ground- and several excited states of metal aromatic clusters, namely NaM(4) and NaM(4) (+/-) (M=Al,Ga,In) clusters have been investigated by employing complete active-space self-consistent-field followed by multireference singles and doubles configuration interaction computations that included up to 10 million configurations and other methods. The ground states NaM(4) (-) of aromatic anions are found to be symmetric C(4nu) ((1)A(1)) electronic states with ideal square pyramid geometries. While the ground state of NaIn(4) is also predicted to be a symmetric C(4nu) ((2)A(1)) square pyramid, the ground state of the NaAl(4) cluster is found to have a C(2nu) ((2)A(1)) pyramid with a rhombus base, and the ground state of NaGa(4) possesses a C(2nu) ((2)A(1)) pyramid with a rectangle base. In general, these structures exhibit two competing geometries, viz., an ideal C(4nu) structure and a distorted rhomboidal or rectangular pyramid structure (C(2nu)). All of the ground states of the NaM(4) (+) (M=Al,Ga,In) cations are computed to be C(2nu) ((3)A(2)) pyramids with rhombus bases. The equilibrium geometries, vibrational frequencies, dissociation energies, adiabatic ionization potentials, adiabatic electron affinities for the electronic states of NaM(4) (M=Al,Ga,In), and their ions are computed and compared with experimental results and other theoretical calculations. On the basis of our computed excited states energy separations, we have tentatively suggested assignments to the observed X and A states in the anion photoelectron spectra of Al(4)Na(-) reported by Li et al. [X. Li, A. E. Kuznetov, H. F. Zheng, A. I. Boldyrev, and L. S. Wang, Science 291, 859 (2001)]. The X state can be assigned to a C(2nu) ((2)A(1)) rhomboidal pyramid. The A state observed in the anion spectrum is assigned to the first excited state ((2)B(1)) of the neutral NaAl(4) with the C(4nu) symmetry. The assignments of the excited states are consistent with the experimental excitation energies and the previous Green's function-based methods for the vertical transition energy separations between the X and A bands.     

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.     

Ab initio spin-orbit configuration interaction calculations for high-lying states of the HeNe quasimolecule

Multireference configuration interaction (MRD-CI) calculations are reported for a large series of electronic states of the HeNe quasimolecule up to 170000 cm(-1) excitation energy, including those that dissociate to the 3S1 and 2 1S0 excited states of the He atom. Spin-orbit coupling is included through the use of relativistic effective core potentials (RECPs). Good agreement is obtained with experimental spectroscopic data for the respective atomic levels, although there is a tendency to systematically underestimate the energies of the Ne atom by 1000-1500 cm(-1) because of differences in the correlation effects associated with its ground and Rydberg excited states. Potential curves are calculated for each of these states, and a number of relatively deep minima are found. The CI Omega-state wave functions are sufficiently diabatic until r = 4-5 a0 to allow for a clear identification of the He 1s-2s excited states. Electric dipole transition moments are computed between these states and the HeNe X 0+ ground state up to r = 4.0 a0, and it is found that the 2 (1)S0 - X maximum value is over an order of magnitude larger than that for the corresponding (3)S1 - X excitation process.     

Electronic structure and thermochemical properties of silicon-doped lithium clusters Li(n)Si0/+, n = 1-8: new insights on their stability

A theoretical investigation on small silicon-doped lithium clusters Li(n)Si with n = 1-8, in both neutral and cationic states is performed using the high accuracy CCSD(T)/complete basis set (CBS) method. Location of the global minima is carried out using a stochastic search method and the growth pattern of the clusters emerges as follows: (i) the species Li(n)Si with n ≤ 6 are formed by directly binding one Li to a Si of the smaller cluster Li(n-1)Si, (ii) the structures tend to have an as high as possible symmetry and to maximize the coordination number of silicon. The first three-dimensional global minimum is found for Li(4)Si, and (iii) for Li(7)Si and Li(8)Si, the global minima are formed by capping Li atoms on triangular faces of Li(6)Si (O(h)). A maximum coordination number of silicon is found to be 6 for the global minima, and structures with higher coordination of silicon exist but are less stable. Heats of formation at 0 K (Δ(f)H(0)) and 298 K (Δ(f)H(298)), average binding energies (E(b)), adiabatic (AIE) and vertical (VIE) ionization energies, dissociation energies (D(e)), and second-order difference in total energy (Δ(2)E) of the clusters in both neutral and cationic states are calculated from the CCSD(T)/CBS energies and used to evaluate the relative stability of clusters. The species Li(4)Si, Li(6)Si, and Li(5)Si(+) are the more stable systems with large HOMO-LUMO gaps, E(b), and Δ(2)E. Their enhanced stability can be rationalized using a modified phenomenological shell model, which includes the effects of additional factors such as geometrical symmetry and coordination number of the dopant. The new model is subsequently applied with consistency to other impure clusters Li(n)X with X = B, Al, C, Si, Ge, and Sn.     

Dissociative recombination of H(3)(+) in the ground and excited vibrational states

The article presents calculated dissociative recombination (DR) rate coefficients for H(3) (+). The previous theoretical work on H(3) (+) was performed using the adiabatic hyperspherical approximation to calculate the target ion vibrational states and it considered just a limited number of ionic rotational states. In this study, we use accurate vibrational wave functions and a larger number of possible rotational states of the H(3) (+) ground vibrational level. The DR rate coefficient obtained is found to agree better with the experimental data from storage ring experiments than the previous theoretical calculation. We present evidence that excited rotational states could be playing an important role in those experiments for collision energies above 10 meV. The DR rate coefficients calculated separately for ortho- and para-H(3) (+) are predicted to differ significantly at low energy, a result consistent with a recent experiment. We also present DR rate coefficients for vibrationally excited initial states of H(3) (+), which are found to be somewhat larger than the rate coefficient for the ground vibrational level.     

Photoelectron spectroscopy of anions at 118.2 nm: observation of high electron binding energies in superhalogens MCl4- (M=Sc, Y, La)

High energy photon is needed for photoelectron spectroscopy (PES) of anions with high electron binding energies, such as superhalogens and O-rich metal oxide clusters. The highest energy photon used for anion PES in the laboratory has been 157 nm (7.866 eV) from F2 eximer lasers. Here, we report an anion PES experiment using coherent vacuum ultraviolet radiation at 118.2 nm (10.488 eV) by tripling the third harmonic output (355 nm) of a Nd:YAG laser in a XeAr cell. Our study focuses on a set of superhalogen species, MCl(4) (-) (M=Sc, Y, La), which were expected to possess very high electron binding energies. While the 157 nm photon can only access the ground state detachment features for these species, more transitions to the excited states at binding energies higher than 8 eV are observed at 118.2 nm. The adiabatic detachment energies are shown to be, 6.84, 7.02, and 7.03 eV for ScCl(4) (-), YCl(4) (-), and LaCl(4) (-) eV, respectively, whereas their corresponding vertical detachment energies are measured to be 7.14, 7.31, and 7.38 eV.