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.     

Theoretical study of the electronic nonadiabatic transitions in the photoelectron spectroscopy of F2O

The photoelectron spectrum of F2O pertaining to ionizations to the ground (X2B1) and low-lying excited electronic states (A2B2, B2A1, and C2A2) of F2O+ is investigated theoretically. The near equilibrium potential energy surfaces of the ground electronic state (X2B1) of F2O and the mentioned ground and excited electronic states of F2O+ reported by Wang et al. ( J. Chem. Phys. 2001, 114, 10682) for the C2v configuration are extended for the Cs geometry assuming a harmonic vibration along the asymmetric stretching mode. The vibronic interactions between the A2B2 and B2A1 electronic states of F2O+ are treated within a linear coupling approach, and the strength of the vibronic coupling parameter is calculated by an ab initio method. The nuclear dynamics is simulated by both time-independent quantum mechanical and time-dependent wave packet approaches. Although the first photoelectron band exhibits resolved vibrational progression along the symmetric stretching mode, the second one is highly overlapping. The latter is attributed to the nonadiabatic interactions among the energetically close A2B2, B2A1, and C2A2 electronic states of F2O+. The theoretical findings are in good accord with the available experimental results.     

Global potential energy surfaces for the Al+(1S) + H2 system

Global, three-dimensional multireference ab initio potential energy surfaces have been calculated for the AlH2+ system for the two lowest energy singlet states and the lowest energy triplet state. These surfaces were calculated using the multireference configuration interaction level of theory with a large basis set. The accuracy of the surfaces were checked against available experimental data and previous theoretical investigations. The areas of surface crossings between the ground state singlet surface and the lowest energy triplet surface and the first excited singlet surface have been thoroughly investigated in all three dimensions and found to give rise to two regions of surface crossings--an "early" crossing (reduced H2 distance) and a "late" crossing (enlarged H2 distance). It is anticipated that both of these crossings will be important in modeling the dynamics of the system. Each of the global potential energy surfaces were fit by interpolation methodology to obtain analytic representations of the surfaces. A representative classical simulation on the ground state singlet surface was performed and discussed.     

Theoretical study of the photoinduced transfer among the ground state and two metastable states in [Fe(CN)5NO]2-

Ab initio calculations were performed to investigate photoinduced transfers among the ground state (GS) and two metastable states (MS1 and MS2) of [Fe(CN)5NO]2-. We obtained the global potential energy surface of the electronic ground state by a scheme of multireference singly and doubly excited configuration interaction followed by a Davidson-type quadruple correction (MRSDCI+Q). The ground state surface has three local minima corresponding to GS, MS1, and MS2. The character of bond between Fe and the nitrosyl group are discussed. We carried out calculations of the lower five electronic excited states by MRSDCI+Q. The main configurations of these lower five excited states were represented by the dFe-->pi*NO transition accompanied by considerable back-donation. The potential energy surfaces of the six states, including the ground state, were obtained by state averaged complete active space self-consistent field calculations. The surfaces have several conical intersections and avoided crossings in the reaction pathway. The photoinduced transfers among GS, MS1, and MS2 are caused by the nonadiabatic effect near these crossings.     

Multireference coupled-cluster study of the symmetry breaking in the C2B radical

The potential energy surfaces (PESs) for both the ground and the excited electronic states of the C(2)B radical are investigated using various multireference (MR) coupled-cluster (CC) approaches. In the ground state case we employ the reduced MR (RMR) CC approach with singles (S) and doubles (D), the RMR CCSD method, as well as its RMR CCSD(T) version corrected for secondary triples, relying on various model spaces and basis sets. The reliability of this approach is also tested against the benchmark full configuration interaction results obtained for a small Dunning-Hay (DH) basis set. The results imply a clear preference for a cyclic structure which, however, breaks the C(2v) symmetry. This symmetry breaking manifests itself strongly at the level of the independent particle model, as represented by the restricted open-shell Hartree-Fock approximation, but the tendency toward symmetry breaking diminishes with the increasing size of the basis set employed as well as with the enhanced account of the correlation effects. It is likely to disappear in the complete basis set limit. The general model space CCSD method is then used to compute vertical excitation energies for a number of excited states as well as the cuts of the PES as the boron atom moves around the C(2) fragment. These results also explain why no symmetry breaking is found when relying on a spin contaminated unrestricted Hartree-Fock reference, as in the UMP2 method.     

QM/MM connection atoms for the multistate treatment of organic and biological molecules

We present an extension of our semiempirical floating occupation MO-CI approach for the determination of ground and excited state potential energy surfaces of interest in photochemistry. The QM/MM variant of the method, which allows for electrostatic and van der Waals interactions between the QM and MM subsystems, is supplemented with a treatment of covalent interactions based on Antes and Thiels connection atom approach. We concentrate on the correct treatment of electrostatic interactions concerning the connection atom, on the specific requirements for the representation of excited states, and on the transferability of the optimal parameters. We show the viability of the method with four examples of connection atoms: S in a thioether bridge, acylic C, aliphatic C, and N in a peptide. The results obtained with the QM/MM treatment compare well with all-QM results of the same level.     

The predicted spectrum of the hypermetallic molecule MgOMg

The present study of MgOMg is a continuation of our theoretical work on Group 2 M(2)O hypermetallic oxides. Previous ab initio calculations have shown that MgOMg has a linear (1)Σ(g)+ ground electronic state and a very low lying first excited triplet electronic state that is also linear; the triplet state has (3)Σ(u)+ symmetry. No gas phase spectrum of this molecule has been assigned, and here we simulate the infrared absorption spectrum for both states. We calculate the three-dimensional potential energy surface, and the electric dipole moment surfaces, of each of the two states using a multireference configuration interaction (MRCISD) approach based on full-valence complete active space self-consistent field (FV-CASSCF) wavefunctions with a cc-pCVQZ basis set. A variational MORBID calculation using our potential energy and dipole moment surfaces is performed to determine rovibrational term values and to simulate the infrared absorption spectrum of the two states. We also calculate the dipole polarizability of both states at their equilibrium geometry in order to assist in the interpretation of future beam deflection experiments. Finally, in order to assist in the analysis of the electronic spectrum, we calculate the vertical excitation energies, and electric dipole transition matrix elements, for six excited singlet states and five excited triplet states using the state-average full valence CASSCF-MRCISD/aug-cc-pCVQZ procedure.     

An Anderson impurity model for efficient sampling of adiabatic potential energy surfaces of transition metal complexes

We present a model intended for rapid sampling of ground and excited state potential energy surfaces for first-row transition metal active sites. The method is computationally inexpensive and is suited for dynamics simulations where (1) adiabatic states are required "on-the-fly" and (2) the primary source of the electronic coupling between the diabatic states is the perturbative spin-orbit interaction among the 3d electrons. The model Hamiltonian we develop is a variant of the Anderson impurity model and achieves efficiency through a physically motivated basis set reduction based on the large value of the d-d Coulomb interaction U(d) and a Lanczos matrix diagonalization routine to solve for eigenvalues. The model parameters are constrained by fits to the partial density of states obtained from ab initio density functional theory calculations. For a particular application of our model we focus on electron transfer occurring between cobalt ions solvated by ammonium, incorporating configuration interaction between multiplet states for both metal ions. We demonstrate the capability of the method to efficiently calculate adiabatic potential energy surfaces and the electronic coupling factor we have calculated compares well to previous calculations and experiment. (     

Singlet excitation energies of thiophene derivatives of fluorene: TD‐DFT study

Fluorene‐thiophene (FT)‐based oligomers and polymers and their derivatives are good candidates for organic blue light‐emitting diodes. In this work, the intrinsic properties of the ground and excited states of FT monomer and its derivatives are studied. The ground‐state optimized structures and energies are obtained using molecular orbital theory and density functional theory (DFT). The ground‐state potential energy curves or surfaces of FT and its derivatives are also obtained. All derivatives are nonplanar in their electronic ground states. The character and energy of the first 20 singlet–singlet electronic transitions are investigated by applying the time‐dependent density functional theory (TD‐DFT) approximations to the correspondingly optimized ground‐state geometries. The lowest singlet state is studied with the configuration interaction (singles) approach (CIS). Excitation energies are red shifted when the FT unit or its derivatives are extended longitudinally. CIS results suggest geometry relaxation in the first singlet excited state. When available, a comparison is made with experimental results. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Predictions of the geometries and fluorescence emission energies of oxyluciferins

The complete active space self-consistent field (CASSCF) method and multiconfigurational second-order perturbation theory (CASPT2) have been used to study the structures and spectra of oxyluciferins (OxyLH2). The ground and lowest-lying singlet excited states geometries have been optimized using CASSCF. CASPT2 has been used to predict relaxed emission energies. The focus is on the lowest-lying singlet excited states of the anionic keto and enol forms of OxyLH2(-1) at the optimized excited-state geometries. The planar keto and enol forms of OxyLH2(-1) are minima on both the S0 and the S1 potential energy surfaces. The twisted keto and enol forms of OxyLH2(-1) are transition states on the S0 and S1 potential energy surfaces. The S1 --> S0 fluorescence emission energies are in the range of 54.2-58.4 kcal/mol for the anionic planar keto forms of OxyLH2, and in the range of 55.7-63.2 kcal/mol for the anionic enol forms of OxyLH2. S0 and S1 potential energy surfaces and thus are not implicated in the emission spectra in the gas phase.     

Theoretical analysis of reaction kinetics with singlet oxygen molecules

A comparative analysis of predictive ability of three approaches to estimate the rate constants of reactions of H(2), H, H(2)O and CH(4) with electronically excited O(2)(a(1)Δ(g)) and O(2)(b(1)Σ(g)(+)) molecules is conducted. The first approach is based on a detailed ab initio study of potential energy surfaces. The second one is known as the "bond energy-bond order" method, and the third approach is a modification of the updated method of vibronic terms that makes it possible to evaluate the activation energy of reactions involving electronically excited species. The comparison showed that the estimates of the energy barrier by the updated method of vibronic terms for some reactions can be in good agreement with ab initio calculations and available experimental data. It was revealed that reactions of O(2)(b(1)Σ(g)(+)) molecules with H(2), H(2)O and CH(4) molecules and with the H atom result in the formation of electronically excited species. The reactivity of O(2)(b(1)Σ(g)(+)) molecules is smaller than that of O(2)(a(1)Δ(g)) ones, but much higher as compared to the reactivity of ground state O(2) molecules. For each reaction under study involving oxygen molecules in the excited electronic states O(2)(a(1)Δ(g)) and O(2)(b(1)Σ(g)(+)) the recommended temperature-dependent rate constants are presented.     

Theoretical studies on the photochemical reaction mechanisms of o-xylylene. Effects of phenyl group for electrocyclic reaction mechanism

The potential energy surfaces (PESs) of the electrocyclic reactions of o-xylylene at the ground and the lowest excited states are calculated by CASSCF molecular orbital and MRMP2 methods. The lowest excited state geometry of o-xylylene has C(2v) symmetry and is about 65 kcal mol(-1) in energy above the ground state. The PESs in the vicinity of the conical intersection are different from those of the electrocyclic reaction of cis-butadiene. In the vicinity of the conical intersection, the transition state at the ground state relating to methylene-cycloheptadienyl carbene is located. The transition state is only 4.3 kcal mol(-1) lower in energy than the conical intersection at the CASSCF(10,10)/6-31G(d) level and 0.5 kcal mol(-1) lower at the MRMP2/6-311+G(d,p) level. The transition state corresponding to benzocyclobutene does not locate in the vicinity of the conical intersection because of the resonance energy between benzene ring and methylene group.     

Ab initio calculations and vibrational energy level fits for the lower singlet potential-energy surfaces of C3

Ab initio multireference configuration interaction potential energy surfaces are computed for the eight lowest singlet surfaces of C(3). These reveal several important features, including several conical intersections in linear, nonlinear, and equilateral triangle geometries. These intersections are important because, particularly for the excited A (1)Pi(u) state, reasonable ab initio results could only be obtained by including nearby, near degenerate, (1)Sigma(u) (-) and (1)Delta(u) states that cross the A (1)Pi(u) state around 4500 cm(-1) above the equilibrium geometry, and a (1)Pi(g) state whose potential in turn crosses the other states about 2000 cm(-1) further up. These states are probably responsible for the complexity of the shorter wavelength UV absorption spectrum of C(3). The computed potential energy surface for the ground, X (1)Sigma(g) (+), state and for the lowest two excited singlet surfaces (which both correlate with the A (1)Pi(u) state in a collinear geometry) are fitted to analytic functional forms. Vibrational energy levels are calculated for both states, taking account of the Renner-Teller coupling in the excited A (1)Pi(u) state. The potential parameters for both states are then least-squares fitted to experimental data. The ground-state fit covers a range of approximately 8500 cm(-1) above the lowest level, and reproduces 100 observed vibrational levels with an average error of 2.8 cm(-1). The A (1)Pi(u) state surfaces cover a range of 3250 cm(-1) above the zero-point level, and reproduce the 44 observed levels in this range with an average error of 2.8 cm(-1).     

Ground and excited state isomeric forms of the anthracene-diethylaniline molecular exciplex

The potential energy surfaces for the molecular complex formed between anthracene (the electron acceptor) andN ,N-diethylaniline (DEA) (the electron donor) were computed as the quasi-adiabatic states resulting from the configuration interaction between the ground (AD), locally excited (A^*D) and charge-transfer (A⁻D⁺) excited electronic configurations. The results clearly indicate the existence of three geometrically and energetically different isomeric forms of the complex in the ground state. In the excited state, the potential energy surfaces reveal the existence of five well-defined equilibrium configurations separated by energy barriers and characterized by different admixtures of the (A^*D) and (A⁻D⁺) electronic configurations. Such a variety of equilibrium configurations in the ground and excited states is, in part, accounted for by the existence of two different conformational forms of DEA that can form complexes with anthracene, and are characterized by different balances between steric effects and interactions of electronic charge distributions in the complex components. The energies of transitions between the relevant ground and excited state equilibrium configurations were calculated and compared with spectroscopic data of a jet-cooled complex obtained in supersonic beam experiments. These transitions were successfully assigned to the observed resonance-like and exciplex-like spectra, and this enabled interpretation of observed changes in the fluorescence excitation and fluorescence spectra of the complex upon excess excitation energy.     

The potential energy surfaces and the radiationless dynamics of excited states of benzene and pyrazine

The potential energy surfaces of the lowest excited states of benzene and pyrazine are investigated as a function of some of the symmetry-adapted internal coordinates by means of the INDO/S method. A large stabilization of the T₂ (ππ*) state of pyrazine (≈ 0.5 eV) along the S_8b vibrational coordinate is found. The calculated potential energy in some excited states (T₁ in benzene, T₂ and S₂ in pyrazine) is a very flat function of the S_16b vibrational coordinate, leading to a crossing with the potential energy of the ground state at relatively small excess of vibrational energy (≈ 1 eV). Thus the ν_16b vibrational mode is postulated to play an important role in the radiationless relaxation to the ground states of these systems. No such crossing has been found near the “channel three” threshold of benzene.     

Hückel－Hubbard参数化法及氮叶立德［2，3］和氢［1，3］σ键迁移的研究

本文建议一种Ｈüｃｋｅｌ－Ｈｕｂｂａｒｄ参数化法，并用Ｈüｃｋｅｌ－Ｈｕｂｂａｒｄ理论首次计算了氮叶立德［２，３］和氢［１，３］σ键迁移反应的基态和低激发态势能面，根据计算得到的势能面，对相应的基态和激发态反应途径进行了讨论，得到有价值的结论。     

Relaxation of excited Franck–Condon states of molecular electron donor–acceptor complexes in liquid solution: the role of orientational isomers

The direct consequences of the presence of ground state orientational isomers of molecular complexes are discussed in terms of the adiabatic potential energy surfaces calculated for the ground and excited states of electron donor–acceptor complexes of tetracyanobenzene with toluene and with mesitylene. Some earlier experimental results that confirm the presence of orientational isomers are also recalled and reviewed, together with the recent results for molecular exciplexes under supersonic molecular beam conditions. Exploration of potential energy surfaces shows that the relaxation pathways of excited Franck–Condon states of the ground state isomers may differ considerably and in liquid solution may be sensitive to physical conditions, which in fact is observed in time-resolved fluorescence spectra of the electron donor–acceptor systems under consideration, upon excitation of high-energy Van der Waals orientational isomers. It is concluded that, in weak electron donor–acceptor complexes in liquid solutions, the role of such isomers may be limited, but it may become crucial for the kinetics and dynamics of excited states if the system is simultaneously capable of forming an exciplex.     

On the nature of the lowest triplet excited state of the [Rh 2(1,3-diisocyanopropane)(4)]2+ ion

The nature of the ground state and the lowest triplet excited state of the [Rh(2)(1,3-diisocyanopropane)(4)](2+) ion have been investigated by the density functional theory. Two locally stable geometrical conformations are found on the potential energy surfaces of both the ground and excited states, corresponding to the eclipsed and twisted conformations, the eclipsed conformation being more stable and having the shorter Rh-Rh bond length. While the Rh-Rh distances of the two conformations differ by approximately 0.4 A, they shorten to the same value upon excitation ( approximately 3.1 A). The excited state originates from the d(z)()()2 (metal antibonding) to p(z)() (ligand-metal bonding) electronic transition. The Mayer Rh-Rh bond order increases from approximately 0.2 to more than 0.8 upon excitation, while the Rh-C(N) bond order shows a slight decrease. A topological bond path between the Rh atoms is found in both the ground and excited states, while the electron localization function (ELF) indicates weak Rh-Rh covalent bonding for the excited state only.     

Nonadiabatic effects in C-Br bond scission in the photodissociation of bromoacetyl chloride

Bromoacetyl chloride photodissociation has been interpreted as a paradigmatic example of a process in which nonadiabatic effects play a major role. In molecular beam experiments by Butler and co-workers [J. Chem. Phys. 95, 3848 (1991); J. Chem. Phys. 97, 355 (1992)], BrCH2C(O)Cl was prepared in its ground electronic state (S0) and excited with a laser at 248 nm to its first excited singlet state (S1). The two main ensuing photoreactions are the ruptures of the C-Cl bond and of the C-Br bond. A nonadiabatic model was proposed in which the C-Br scission is strongly suppressed due to nonadiabatic recrossing at the barrier formed by the avoided crossing between the S1 and S2 states. Recent reduced-dimensional dynamical studies lend support to this model. However, another interpretation that has been given for the experimental results is that the reduced probability of C-Br scission is a consequence of incomplete intramolecular energy redistribution. To provide further insight into this problem, we have studied the energetically lowest six singlet electronic states of bromoacetyl chloride by using an ab initio multiconfigurational perturbative electronic structure method. Stationary points (minima and saddle points) and minimum energy paths have been characterized on the S0 and S1 potential energy surfaces. The fourfold way diabatization method has been applied to transform five adiabatic excited electronic states to a diabatic representation. The diabatic potential energy matrix of the first five excited singlet states has been constructed along several cuts of the potential energy hypersurfaces. The thermochemistry of the photodissociation reactions and a comparison with experimental translational energy distributions strongly suggest that nonadiabatic effects dominate the C-Br scission, but that the reaction proceeds along the energetically allowed diabatic pathway to excited-state products instead of being nonadiabatically suppressed. This conclusion is also supported by the low values of the diabatic couplings on the C-Br scission reaction path. The methodology established in the present study will be used for the construction of global potential energy surfaces suitable for multidimensional dynamics simulations to test these preliminary interpretations.     

Theoretical investigation on the electronic and geometric structure of GaN2+ and GaN4+

The electronic and geometric structures of gallium dinitride cation, GaN2+ and gallium tetranitride cation, GaN4+ were systematically studied by employing density functional theory (DFT-B3LYP) and perturbation theory (MP2, MP4) in conjunction with large basis sets, (aug-)cc-pVxZ, x = T, Q. A total of 7 structures for GaN2+ and 24 for GaN4+ were identified, corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces, and bonding mechanisms for some low-lying electronic states. The calculated dissociation energy (De) of the ground state of GaN2+, X1Sigma+, is 5.6 kcal/mol with respect to Ga+(1S) + N2(X1Sigmag+) and that of the excited state, ?3Pi, is 24.8 kcal/mol with respect to Ga+(3P) + N2(X1Sigmag+). The ground state and the first excited minimum of GaN4+ are of 1A1(C2v) and 3B1(C2v) symmetry with corresponding De of 11.0 and 43.7 kcal/mol with respect to Ga+(1S) + 2N2(X1Sigmag+) for X1A1 and Ga+(3P) + 2N2(X1Sigmag+) for 3B1.