Electronic and geometric structure of the 3d-transition metal monocarbonyls MCO, M=Sc, Ti, V, and Cr

The electronic and geometric structure of the 3d-transition metal monocarbonyls MCO, M=Sc, Ti, V, and Cr was investigated through coupled cluster (CC) and multireference variational methods (MRCI) combined with large basis sets. For the ground and a few low-lying excited states complete potential energy profiles were constructed at the CC-level of theory. The M-CO dissociation energies of the ground states X 4Sigma-,X 5Delta,X 6Sigma+, and X 7A' are calculated to be 36, 27, 18, and 2 kcal/mol for ScCO, TiCO, VCO, and CrCO, with respect to Sc(4F),Ti(5F),V(6D),Cr(7S)+CO(X 1Sigma+). The bonding is rather complicated and could be attributed mainly to pi-conjugation effects between the M and CO pi-electrons, along with weak sigma-charge transfer from CO to M atoms. Almost in all cases the metal atoms appear to be slightly positively charged, at least according to the direction of the dipole moment vectors and the MRCI population densities.     

Electronic structure and bonding of the 3d transition metal borides, MB, M=Sc, Ti, V, Cr, Mn, Fe, Co, Ni, and Cu through all electron ab initio calculations

The electronic structure and bonding of the ground and some low-lying states of all first row transition metal borides (MB), ScB, TiB, VB, CrB, MnB, FeB, CoB, NiB, and CuB have been studied by multireference configuration interaction (MRCI) methods employing a correlation consistent basis set of quintuple cardinality (5Z). It should be stressed that for all the above nine molecules, experimental results are essentially absent, whereas with the exception of ScB and CuB the remaining seven species are studied theoretically for the first time. We have constructed full potential energy curves at the MRCI/5Z level for a total of 27 low-lying states, subsequently used to extract binding energies, spectroscopic parameters, and bonding schemes. In addition, some 20 or more states for every MB species have been examined at the MRCI/4Z level of theory. The ground state symmetries and corresponding binding energies (in kcal/mol) are 5Sigma-(ScB), 76; 6Delta(TiB), 65; 7Sigma+(VB), 55; 6Sigma+(CrB), 31; 5Pi(MnB), 20; 4Sigma-(FeB), 54; 3Delta(CoB), 66; 2Sigma+(NiB), 79; and 1Sigma+(CuB), 49.     

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相似文献   

The electronic structure of transition metal dihelide dications

Multi-reference configuration interaction (MRCI) calculations have been employed to characterize the low-lying states of first-row transition metal dihelide dications, He(2)TM(2+) (TM = Sc-Cu). The most important state-ordering principles were determined to be the occupation of the 4s orbital and orientation of the occupied 3d orbital. The ground states of all species are predicted to be of D(infinityh) symmetry arising from a 3d(n+1) electronic configuration. For excited states with singly occupied 4s or doubly occupied 3d(sigma) orbitals, bending to C(2v) symmetry typically lowers the energy and shortens the He-TM bond length. Coupled cluster singles and doubles with a perturbative treatment of triple excitations (CCSD(T)) results for ground state spectroscopic properties are in agreement with the MRCI predicted trends.     

Ab initio investigation of the electronic and geometric structure of magnesium diboride, MgB2

Employing multireference variational (MRCI) and coupled cluster (CC) methods combined with quadruple-zeta quality correlation-consistent basis set, we have studied 36 states of the magnesium diboride (MgB(2)) molecule as well as 17 states of the experimentally unknown diatomic MgB. For all states of MgB(2), we report geometries, atomization energies, and dipole moments, while for the first 5 states, potential energy profiles have been also constructed. The ground state is formally of (1)A(1) V-shaped symmetry with an atomization energy of 108.1(109) kcal/mol at the MRCI(MRCI + Davidson correction) level. The first excited state ((3)B(1)) is less than 1 kcal/mol above the X(1)A(1) state, with the next state of linear Mg-B-B geometry (b(3)Sigma(-)) located 10 kcal/mol higher. In all states, bent or linear, the bonding is complicated and unconventional because of the extraordinary bonding agility of the boron atom(s).     

Approximately size extensive local Multireference Singles and Doubles Configuration Interaction

Multi-reference Configuration Interaction (MRCI) is often used to predict the electronic structures and reaction energetics of small molecules with very high accuracy. Unfortunately, MRCI is inapplicable to large or even medium-sized molecules for two reasons: its computational cost scales poorly with molecule size and MRCI methods are not size extensive, leading to large energy errors. We have developed a new local (L) and approximately size extensive MRCI method that addresses both shortcomings. Truncating long-range electron correlation in a local orbital basis as well as efficient processing of two-electron integrals via Cholesky decomposition (CD) and integral screening reduce the computational cost to O(N(3)) with a small prefactor. A priori and a posteriori size extensivity corrections are both considered. The Multi-reference Averaged Coupled-Pair Functional (MRACPF) provides approximate size extensivity by modifying the energy functional. The very inexpensive Davidson-Silver and Pople a posteriori schemes also produce quite accurate corrections over a large range of molecular size. Employing CD-LMRACPF is slightly more expensive than using a Davidson-type correction, but the former gives superior results. Molecules with up to 50 heavy atoms can be treated with our CD-LMRACPF method thus far.     

All-electron first principles calculations of the ground and some low-lying excited states of BaI

The electronic structure of the heavy diatomic molecule BaI has been examined for the first time by ab initio multiconfigurational configuration interaction (MRCI) and coupled cluster (RCCSD(T)) methods. The effects of special relativity have been taken into account through the second-order Douglas-Kroll-Hess approximation. The construction of Omega(omega,omega) potential energy curves allows for the estimation of "experimental" dissociation energies (De) of the first few excited states by exploiting the accurately known De experimental value of the X2Sigma+ ground state. All states examined are of ionic character with a Mulliken charge transfer of 0.5 e- from Ba to I, and this is reflected to large dipole moments ranging from 6 to 11 D. Despite the inherent difficulties of a heavy system like BaI, our results are encouraging. With the exception of bond distances that on the average are calculated 0.05 A longer than the experimental ones, common spectroscopic parameters are in fair agreement with experiment, whereas De values are on the average 10 kcal/mol smaller.     

On the ground and some low-lying excited states of ScB: a multiconfigurational study

The electronic structure of a series of low-lying excited triplet and quintet states of scandium boride (ScB) was examined using multireference configuration interaction (including Davidson's correction for quadruple excitations) and single-reference coupled cluster (CC) methods with averaged natural orbital (ANO) basis sets. The CC approach was used only for the lowest quintet state. The authors have analyzed eight low-lying triplets 3Sigma-(2), 3Sigma+, 3Pi(3), and 3Delta(2) dissociating to Sc(2D)/B(2P) atoms and eight low-lying quintet states 5Sigma-, 5Sigma+, 5Pi(2), 5Phi, and 5Delta(3) dissociating to Sc(4F)/B(2P) atoms. They report the potential energy curves and spectroscopic parameters of ScB obtained with the multireference configuration interaction (MRCI) technique including all singly and doubly excited configurations obtained with the ANO-S basis set. For the two lowest states they obtained also improved ANO-L spectroscopic constants, dipole and quadrupole moments as well as scalar relativistic effects based on the Douglas-Kroll-Hess Hamiltonian. They provide the analysis of the bonding based on Mulliken populations and occupation numbers. Since the two lowest states, 3Sigma- and 5Sigma-, lie energetically very close, their principal goal was to resolve the nature of the ground state of ScB. Their nonrelativistic MRCI(Q) (including Davidson correction) results indicate that the quintet is more stable than the triplet by about 800 cm(-1). Inclusion of scalar relativistic effects reduces this difference to about 240 cm(-1). The dissociation energies for 5Sigma- ScB range from 3.20 to 3.30 eV while those for the 3Sigma- range from 1.70 to 1.80 eV.     

Electronic structure of cobalt carbide, CoC

The ground and 18 low lying excited states of the diatomic molecule cobalt carbide, CoC, have been examined by multireference variational methods (MRCI) combined with quantitative basis sets. All calculated states are bound and correlate adiabatically to the ground-state atoms, Co(a4F) + C(3P). We report complete potential energy curves, equilibrium bond distances, dissociation energies (De), spectroscopic constants, electric dipole moments and spin-orbit splittings. The bonding character of certain states is also discussed with the help of Mulliken distributions and valence-bond-Lewis diagrams. We are practically certain that the ground state is of 2Sigma+ symmetry with a state of 2Delta symmetry lying less than 3 kcal/mol higher, in agreement with the relevant experimental findings. Our best estimate of the X 2Sigma+ dissociation energy is De(D0) = 83(82) kcal/mol at r(e) = 1.541 A, 0.02 A shorter than the experimental bond length.     

Multireference configuration interaction studies on higher valence and Rydberg states of OClO, ionization potentials, and electron detachment energies

MRCI results are reported for the vertical excitation energies (VEE) and oscillator strengths f of doublet states of OClO up to 11 eV, including 3b(1) → 4s, 4p, 3d, 5s, 5p, 4d, and most 1a(2), 8a(1), 5b(2) → 4s and 4p Rydberg states. The lowest Rydberg states 3b(1) → 4s and 3b(1) → 4p(x) have mixed valence-Rydberg character. The observed spectral bands were reassigned to include valence states which have generally higher oscillator strengths. The well-known valence state 1(2)A(2) has a VEE of 3.63 eV, and a relatively high f of 0.042. Overall, the calculated oscillator strengths are in good agreement with measured values. The lowest quartet state, 1(4)B(2), lies at 6.95 eV. Quartet Rydberg states start with 1a(2) → 4s at 9.28 eV. According to calculated vertical ionization potentials (VIP) of OClO, the second VIP at 12.59 eV is reassigned from 1(3)B(1) to 1(3)B(2) (ionization from 1a(2), rather than 8a(1)), and the third VIP at 12.63 eV from 1(1)B(1) to 1(3)B(1) (ionization from 8a(1)). Vertical electron detachment energies of OClO(-) have been calculated up to 8.9 eV. There is good agreement with experimental values.     

A comprehensive theoretical investigation of the electronic states of Ca2 up to the Ca(4s(2) 1S) + Ca(4s5p 1P) dissociation limit

A theoretical survey of the electronic structure of Ca(2) is presented using two-electron pseudopotentials complemented by core-polarization operators on Ca atoms and multireference configuration interaction/quasidegenerate perturbation theory (MRCI/QDPT) treatment of molecular excited states. The spectroscopic constants of 70 electronic states up to 30,000?cm(-1) above the ground state are determined. This implies all Ca(2) states dissociating up to the Ca(4s(2) (1)S) + Ca(4s5p (3,1)P) dissociation limits. All spin states (singlet, triplet, and quintet) are investigated. The work emphasizes the variety of interactions implying singly valence and lowest Rydberg excited states, doubly excited states generated by atom pairs (3)P(4s4p) + (3)P(4s4p), or (3)P(4s4p) + (3)D(4s3d), 4p3d double excitations asymptotically localized on a single-atom. Zwitterionic Ca(+) + Ca(-) configurations are evidenced and shown to induce specific electronic patterns in (1)Σ(g)(+), (3)Σ(g)(+), (1)Σ(u)(+), (3)Σ(u)(+), (1)Π(g), (3)Π(g), (1)Π(u), and (3)Π(u) symmetry manifolds. They also provide insight for qualitative features (barriers) found for the lower electronic states already investigated in previous publications by other authors.     

Ab initio study on the oxidation of NCN by O (3P): prediction of the total rate constant and product branching ratios

The reaction of NCN with O is relevant to the formation of prompt NO according to the new mechanism, CH+N2-->cyclic-C(H)NN- -->HNCN-->H+NCN. The reaction has been investigated by ab initio molecular orbital and transition state theory calculations. The mechanisms for formation of possible product channels involved in the singlet and triplet potential energy surfaces have been predicted at the highest level of the modified GAUSSIAN-2 (G2M) method, G2M (CC1). The barrierless association/dissociation processes on the singlet surface were also examined with the third-order Rayleigh-Schr?dinger perturbation (CASPT3) and the multireference configuration interaction methods including Davidson's correction for higher excitations (MRCI+Q) at the CASPT3(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) and MRCI+Q(6,6)/6-311+G(3df)//UB3LYP/6-311G(d) levels. The rate constants for the low-energy channels producing CO+N2, CN+NO, and N(4S)+NCO have been calculated in the temperature range of 200-3000 K. The results show that the formation of CN+NO is dominant and its branching ratio is over 99% in the whole temperature range; no pressure dependence was noted at pressures below 100 atm. The total rate constant can be expressed by: kt=4.23x10(-11) T0.15 exp(17/T) cm3 molecule(-1) s(-1).     

Quantum chemical investigation of the electronic spectra of the keto, enol, and keto-imine tautomers of cytosine

The low-lying excited singlet states of the keto, enol, and keto-imine tautomers of cytosine have been investigated employing a combined density functional/multireference configuration interaction (DFT/MRCI) method. Unconstrained geometry optimizations have yielded out-of-plain distorted structures of the pi --> pi and n --> pi excited states of all cytosine forms. For the keto tautomer, the DFT/MRCI adiabatic excitation energy of the pi --> pi state (4.06 eV including zero-point vibrational energy corrections) supports the resonant two-photon ionization (R2PI) spectrum (Nir et al. Phys. Chem. Chem. Phys. 2002, 5, 4780). On its S1 potential energy surface, a conical intersection between the 1pipi state and the electronic ground state has been identified. The barrier height of the reaction along a constrained minimum energy path amounts to merely 0.2 eV above the origin and explains the break-off of the R2PI spectrum. The 1pipi minimum of the enol tautomer is found at considerably higher excitation energies (4.50 eV). Because of significant geometry shifts with respect to the ground state, long vibrational progressions are expected, in accord with experimental observations. For the keto-imine tautomer, a crossing of the 1pipi potential energy surface with the ground-state surface has been found, too. Its n --> pi minimum (3.27 eV) is located well below the conical intersection between the pi --> pi and S0 states, but it will be difficult to observe because of its small transition moment. The identified conical intersections of the pi --> pi excited states of the keto cytosine tautomers are made responsible for the ultrafast decay to the electronic ground states and thus may explain their subpicoseconds lifetimes.     

Photostability via a sloped conical intersection: a CASSCF and RASSCF study of pyracylene

An extensive ab initio study of the ground- and excited-state potential energy surfaces of pyracylene is presented in this work. CASSCF calculations show that there is an accessible sloped S0/S1 conical intersection, which leads to ultrafast internal conversion and explains the observed photostability. RASSCF calculations (using a well-defined subset of the CASSCF configurations) are shown to be able to reproduce CASSCF results satisfactorily and will therefore be useful for larger systems where CASSCF is currently too expensive. MRCI and MRPT2 energy corrections are computed to assess the ionic character of the excited states. Finally, MMVB calculations are also benchmarked against CASSCF, to assess the reliability of this parametrized method for excited states of large conjugated polycyclic aromatic hydrocarbons.     

Ab initio study of KN

The potential energy curves for the lowest (3)Sigma(-), (3)Pi, and (5)Sigma(-) states of the KN molecule have been calculated by the multireference singles and doubles configuration interaction method, including Davidson's corrections for quadruple excitations [MRCI(+Q)]. It is shown that the former two are bound, while the last one is repulsive. The electronic ground state of KN is predicted as (3)Sigma(-) state, although the term energy of the (3)Pi state is very small, 177.3 cm(-1). The binding energy for the (3)Sigma(-) state is evaluated as 0.838 eV, the rotational constant B(0) as 0.250 63 cm(-1), and harmonic frequency as 324.4 cm(-1). The spin-orbit coupling effects between the (3)Sigma(-) and (3)Pi states of KN are evaluated and discussed. The same MRCI(+Q) computational procedures are applied to the isovalent LiN, KC, KO, and KCl to confirm the accuracy of present calculations. Theoretical spectroscopic constants presented here will inspire experimental studies of KN.     

Theoretical electronic structure of the lowest-lying states of the YI molecule

CAS-SCF/MRCI calculations have been performed for 15 molecular states in the representation ^2S+1Λ^(+/−) (neglecting spin–orbit effects) for the molecule YI. The corresponding 33 molecular states in the representation Ω^(+/−) (including spin–orbit effects) have been calculated using a semi-empirical spin–orbit pseudopotential built up for yttrium. Calculated potential energy curves and spectroscopic constants are reported, to the best of our knowledge they are the first ones from ab initio methods for this molecule. Present results are compared to experimental accurate data available for the ground X¹Σ⁺ and 3 excited states (1)¹Π, (2)¹Σ⁺ and (2)¹Π.     

Intramolecular hydrogen bonding plays a crucial role in the photophysics and photochemistry of the GFP chromophore

In commonly studied GFP chromophore analogues such as 4-(4-hydroxybenzylidene)-1,2-dimethyl-1H-imidazol-5(4H)-one (PHBDI), the dominant photoinduced processes are cis-trans isomerization and subsequent S(1) → S(0) decay via a conical intersection characterized by a highly twisted double bond. The recently synthesized 2-hydroxy-substituted isomer (OHBDI) shows an entirely different photochemical behavior experimentally, since it mainly undergoes ultrafast intramolecular excited-state proton transfer, followed by S(1) → S(0) decay and ground-state reverse hydrogen transfer. We have chosen 4-(2-hydroxybenzylidene)-1H-imidazol-5(4H)-one (OHBI) to model the gas-phase photodynamics of such 2-hydroxy-substituted chromophores. We first use various electronic structure methods (DFT, TDDFT, CC2, DFT/MRCI, OM2/MRCI) to explore the S(0) and S(1) potential energy surfaces of OHBI and to locate the relevant minima, transition state, and minimum-energy conical intersection. These static calculations suggest the following decay mechanism: upon photoexcitation to the S(1) state, an ultrafast adiabatic charge-transfer induced excited-state intramolecular proton transfer (ESIPT) occurs that leads to the S(1) minimum-energy structure. Nearby, there is a S(1)/S(0) minimum-energy conical intersection that allows for an efficient nonadiabatic S(1) → S(0) internal conversion, which is followed by a fast ground-state reverse hydrogen transfer (GSHT). This mechanism is verified by semiempirical OM2/MRCI surface-hopping dynamics simulations, in which the successive ESIPT-GSTH processes are observed, but without cis-trans isomerization (which is a minor path experimentally with less than 5% yield). These gas-phase simulations of OHBI give an estimated first-order decay time of 476 fs for the S(1) state, which is larger but of the same order as the experimental values measured for OHBDI in solution: 270 fs in CH(3)CN and 230 fs in CH(2)Cl(2). The differences between the photoinduced processes of the 2- and 4-hydroxy-substituted chromophores are attributed to the presence or absence of intramolecular hydrogen bonding between the two rings.     

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.     

叠氮化氰的光解离机理

采用多参考态方法,在MRCI+Q//CAS(10,9)/6-311+G(2df)水平上对叠氮化氰(N3CN)的光解离机理进行理论研究.优化得到基态(S0)和低激发态(S1、S2、T1)势能面上的极小点、过渡态、内转换交叉点(IC-S1/S0)和隙间窜跃交叉点(ISC-S1/T1)的结构和能量,构建反应势能面.在MRCI+Q//CAS(10,9)水平上计算N3CN的垂直激发能,并和实验值进行对比.结果表明,在S0、S1、S2和T1态势能面上,N—N键断裂生成N2+NCN是主要解离途径,而C—N键断裂通道是次要通道.实验观测到220 nm处的吸收峰对应分子由S0态到S1态的激发,对应主要光解离产物为NCN[a1△g];而在275 nm处的吸收峰则对应分子被激发到T1态,然后直接生成基态产物NCN[X3Σg-].我们的理论结果与实验测量符合得很好.     

Electronic states and spin-forbidden cooling transitions of AlH and AlF

The feasibility of laser cooling AlH and AlF is investigated using ab initio quantum chemistry. All the electronic states corresponding to the ground and lowest two excited states of the Al atom are calculated using multi-reference configuration interaction (MRCI) and the large AV6Z basis set for AlH. The smaller AVQZ basis set is used to calculate the valence electronic states of AlF. Theoretical Franck-Condon factors are determined for the A(1)Π→ X(1)Σ(+) transitions in both radicals and found to agree with the highly diagonal factors found experimentally, suggesting computational chemistry is an effective method for screening suitable laser cooling candidates. AlH does not appear to have a transition quite as diagonal as that in SrF (which has been laser cooled) but the A(1)Π→ X(1)Σ(+) transition transition of AlF is a strong candidate for cooling with just a single laser, though the cooling frequency is deep in the UV. Furthermore, the a(3)Π→ X(1)Σ(+) transitions are also strongly diagonal and in AlF is a practical method for obtaining very low final temperatures around 3 μK.