Relativistic potential energy surfaces of initial oxidations of Si(100) by atomic oxygen: the importance of surface dimer triplet state

The non-relativistic and relativistic potential energy surfaces (PESs) of the symmetric and asymmetric reaction paths of Si(100)-2×1 oxidations by atomic oxygen were theoretically explored. Although only the singlet PES turned out to exist as a major channel leading to "on-dimer" product, both the singlet and triplet PESs leading to "on-top" products are attractive. The singlet PESs leading to the two surface products were found to be the singlet combinations (open-shell singlet) of the low-lying triplet state of surface silicon dimer and the ground (3)P state of atomic oxygen. The triplet state of the "on-top" product can also be formed by the ground singlet state of the surface silicon dimer and the same (3)P oxygen. The attractive singlet PESs leading to the "on-dimer" and "on-top" products made neither the intersystem crossings from triplet to singlet PES nor high energy (1)D of atomic oxygen necessary. Rather, the low-lying triplet state of surface silicon dimer plays an important role in the initial oxidations of silicon surface.     

Theory of the photodissociation of ozone in the Hartley continuum: potential energy surfaces, conical intersections, and photodissociation dynamics

Ab initio potential energy and transition dipole moment surfaces are presented for the five lowest singlet even symmetry electronic states of ozone. The surfaces are calculated using the complete active space self consistent field method followed by contracted multireference configuration interaction (MRCI) calculations. A slightly reduced augmented correlation consistent valence triple-zeta orbital basis set is used. The ground and excited state energies of the molecule have been computed at 9282 separate nuclear geometries. Cuts through the potential energy surfaces, which pass through the geometry of the minimum of the ground electronic state, show several closely avoided crossings. Close examination, and higher level calculations, very strongly suggests that some of these seemingly avoided crossings are in fact associated with non-symmetry related conical intersections. Diabatic potential energy and transition dipole moment surfaces are created from the computed ab initio adiabatic MRCI energies and transition dipole moments. The transition dipole moment connecting the ground electronic state to the diabatic B state surface is by far the strongest. Vibrational-rotational wavefunctions and energies are computed using the ground electronic state. The energy level separations compare well with experimentally determined values. The ground vibrational state wavefunction is then used, together with the diabatic B<--X transition dipole moment surface, to form an initial wavepacket. The analysis of the time-dependent quantum dynamics of this wavepacket provides the total and partial photodissociation cross sections for the system. Both the total absorption cross section and the predicted product quantum state distributions compare well with experimental observations. A discussion is also given as to how the observed alternation in product diatom rotational state populations might be explained.     

Singlet and Triplet Excited States and Intersystem Crossing in Free‐Base Porphyrin: TDDFT and DFT/MRCI Study

Extensive time-dependent DFT (TDDFT) and DFT/multireference configuration interaction (MRCI) calculations are performed on the singlet and triplet excited states of free-base porphyrin, with emphasis on intersystem crossing processes. The equilibrium geometries, as well as the vertical and adiabatic excitation energies of the lowest singlet and triplet excited states are determined. Single and double proton-transfer reactions in the first excited singlet state are explored. Harmonic vibrational frequencies are calculated at the equilibrium geometries of the ground state and of the lowest singlet and triplet excited states. Furthermore, spin–orbit coupling matrix elements of the lowest singlet and triplet states and their numerical derivatives with respect to nuclear displacements are computed. It is shown that opening of an unprotonated pyrrole ring as well as excited-state single and double proton transfer inside the porphyrin cavity lead to crossings of the potential energy curves of the lowest singlet and triplet excited states. It is also found that displacements along out-of-plane normal modes of the first excited singlet state cause a significant increase of the 2|H_so|S₁>, 1|H_so|S₁>, and 1|H_so|S₀> spin–orbit coupling matrix elements. These phenomena lead to efficient radiationless deactivation of the lowest excited states of free-base porphyrin via intercombination conversion. In particular, the S₁→T₁ population transfer is found to proceed at a rate of ≈10⁷ s⁻¹ in the isolated molecule.     

Ab initio multireference investigation of disjoint diradicals: singlet versus triplet ground states

We present improved virtual orbital (IVO) complete active space (CAS) configuration interaction (IVO‐CASCI) and IVO‐CASCI‐based multireference Møller–Plesset perturbation theory (MRMPPT) calculations with an aim to elucidate the electronic structure of tetramethyleneethane (TME) in its lowest singlet and triplet state and to quantify their order and extent of splitting. The potential surfaces of singlet and triplet states for the twisting of TME are also studied. We found that the triplet state is higher in energy than the singlet one in the whole range of twisting angles with the energy gap minimum at a twisting angle of about 45°. Harmonic vibrational frequencies of TME have also been calculated for both the states. We also report the ground to first excited triplet state transition energies. Our results are analyzed with respect to the results available in the literature to illustrate the efficacy of our methods employed. We also demonstrate that the spin character of the ground state of disjoint, TME‐like diradicals can be manipulated by using appropriate selection of annulenic spacer to separate the allyl groups of TME.     

Singlet and triplet potential surfaces for the O2+C2H4 reaction

Electronic structure calculations at the CASSCF and UB3LYP levels of theory with the aug-cc-pVDZ basis set were used to characterize structures, vibrational frequencies, and energies for stationary points on the ground state triplet and singlet O(2)+C(2)H(4) potential energy surfaces (PESs). Spin-orbit couplings between the PESs were calculated using state averaged CASSCF wave functions. More accurate energies were obtained for the CASSCF structures with the MRMP2/aug-cc-pVDZ method. An important and necessary aspect of the calculations was the need to use different CASSCF active spaces for the different reaction paths on the investigated PESs. The CASSCF calculations focused on O(2)+C(2)H(4) addition to form the C(2)H(4)O(2) biradical on the triplet and singlet surfaces, and isomerization reaction paths ensuing from this biradical. The triplet and singlet C(2)H(4)O(2) biradicals are very similar in structure, primarily differing in their C-C-O-O dihedral angles. The MRMP2 values for the O(2)+C(2)H(4)→C(2)H(4)O(2) barrier to form the biradical are 33.8 and 6.1 kcal/mol, respectively, for the triplet and singlet surfaces. On the singlet surface, C(2)H(4)O(2) isomerizes to dioxetane and ethane-peroxide with MRMP2 barriers of 7.8 and 21.3 kcal/mol. A more exhaustive search of reaction paths was made for the singlet surface using the UB3LYP/aug-cc-pVDZ theory. The triplet and singlet surfaces cross between the structures for the O(2)+C(2)H(4) addition transition states and the biradical intermediates. Trapping in the triplet biradical intermediate, following (3)O(2)+C(2)H(4) addition, is expected to enhance triplet→singlet intersystem crossing.     

Study of acyl group migration by femtosecond transient absorption spectroscopy and computational chemistry

The primary photophysical and photochemical processes in the photochemistry of 1-acetoxy-2-methoxyanthraquinone (1a) were studied using femtosecond transient absorption spectroscopy. Excitation of 1a at 270 nm results in the population of a set of highly excited singlet states. Internal conversion to the lowest singlet npi* excited state, followed by an intramolecular vibrational energy redistribution (IVR) process, proceeds with a time constant of 150 +/- 90 fs. The 1npi* excited state undergoes very fast intersystem crossing (ISC, 11 +/- 1 ps) to form the lowest triplet pipi* excited state which contains excess vibrational energy. The vibrational cooling occurs somewhat faster (4 +/- 1 ps) than ISC. The primary photochemical process, migration of acetoxy group, proceeds on the triplet potential energy surface with a time constant of 220 +/- 30 ps. The transient absorption spectra of the lowest singlet and triplet excited states of 1a, as well as the triplet excited state of the product, 9-acetoxy-2-methoxy-1,10-anthraquinone (2a), were detected. The assignments of the transient absorption spectra were supported by time-dependent DFT calculations of the UV-vis spectra of the proposed intermediates. All of the stationary points for acyl group migration on the triplet and ground state singlet potential energy surfaces were localized, and the influence of the acyl group substitution on the rate constants of the photochemical and thermal processes was analyzed.     

Theoretical Investigation for Two-state Reactivity of CO_2 Hydrogenation Catalyzed by Ru in Gas Phase

Gas-phase CO_2 catalyzed activation hydrogenation by Ru atoms was studied with density functional theory. Based on the structure optimization of different potential energy surfaces,there are two crossing points between singlet and triplet potential energy surfaces and there is a crossing point between quintet and triplet potential energy surfaces in the whole catalytic cycle. Spin transition probabilities in the vicinity of the intersections have been calculated by the Landau-Zener model theory. There are three minimum energy crossing points(MECPs) with strong spin-orbital coupling effect and higher spin transition probability,and all spin inversion occurred in s orbital and different d orbitals of ruthenium,indicating this is a typical two-state reactivity(TSR) reaction. Finally,the lowest energy reaction path is ensured.     

Mechanism of the photochemical process of singlet oxygen production by phenalenone

Phenalenone (PN) is a very efficient singlet oxygen sensitiser in a wide range of solvents. This work uses ab initio quantum chemical calculations (CASSCF/CASPT2 protocol) to study the mechanism for populating the triplet state of PN responsible for this reaction, the (3)(π-π*) state. To describe in detail this reaction path, the singlet and triplet low-lying excited states of PN have been studied, the critical points of the potential energy surfaces corresponding to these states located and the vertical and adiabatic energies calculated. Our results show that, after the initial population of the S(2) excited state of (π-π*) character, the system undergoes an internal conversion to the (1)(n-π*) state. After populating the dark S(1) state, the system relaxes to the (1)(n-π*) minimum, but rapidly populates the triplet manifold through a very efficient intersystem crossing to the (3)(π-π*) state. Although the population of the minimum of this triplet state is strongly favoured, a conical intersection with the (3)(n-π*) surface opens an internal conversion channel to this state, a path accessible only at high temperatures. Radiationless deactivation processes are ruled out on the basis of the high-energy barriers found for the crossings between the excited states and the ground state. Our computational results satisfactorily explain the experimental findings and are in very good agreement with the experimental data available. In the case of the frequency of fluorescence, this is the first time that these data have been theoretically predicted in good agreement with the experimental results.     

Triplet-state formation along the ultrafast decay of excited singlet cytosine

We address the possibility of populating the lowest triplet state of cytosine by an "intrinsic"mechanism, namely, intersystem crossing (ISC) along the ultrafast internal conversion pathway of the electronically excited singlet species. For this purpose, we present a discussion of the ISC process and triplet-state reactivity based on theoretical analysis of the spin-orbit strength and the potential energy surfaces for the relevant singlet and triplet states of cytosine. High-level ab initio computations show that ISC is possible in wide regions of the singlet manifold along the reaction coordinate that controls the ultrafast internal conversion to the ground state. Thus, the ISC mechanism documented here provides a possibility to access the triplet state, which has a key role in the photochemistry of the nucleic acid bases.     

On the reversible O2 binding of the Fe-porphyrin complex

Electronic mechanism of the reversible O(2) binding by heme was studied by using Density Functional Theory calculations. The ground state of oxyheme was calculated to be open singlet state [Fe(S =1/2) + O(2)(S = 1/2)]. The potential energy surface for singlet state is associative, while that for triplet state is dissociative. Because the ground state of the O(2)+ deoxyheme system is triplet in the dissociation limit [Fe(S = 2) + O(2)(S = 1)], the O(2) binding process requires relativistic spin-orbit interaction to accomplish the intersystem crossing from triplet to singlet states. Owing to the singlet-triplet crossing, the activation energies for both O(2) binding and dissociation become moderate, and hence reversible. We also found that the deviation of the Fe atom from the porphyrin plane is also important reaction coordinate for O(2) binding. The potential surface is associative/dissociative when the Fe atom locates in-plane/out-of-plane.     

Density functional theory study of 1,2‐dioxetanone decomposition in condensed phase

The decomposition of 1,2‐dioxetanone into a CO₂ molecule and into an excited state formaldehyde molecule was studied in condensed phase, using a density functional theory approach. Singlet and triplet ground and excited states were all included in the calculations. The calculations revealed a novel mechanism for the chemiluminescence of this compound. The triplet excitation can be explained by two intersystem crossings (ISCs) with the ground state, while the singlet excitation can be accounted by an ISC with the triplet state. The experimentally verified small excitation yield can then be explained by the presence of an energy barrier present in the potential energy surface of the triplet excited state, which will govern both triplet and singlet excitation. It was also found that the triplet ground state interacts with both the triplet excited and singlet ground states. A MPWB1K/mPWKCIS approach provided results in agreement with the existent literature. © 2012 Wiley Periodicals, Inc.     

Model studies of hydrogen atom addition and abstraction processes involving ortho-, meta-, and para-benzynes.

H-atom addition and abstraction processes involving ortho-, meta-, and para-benzyne have been investigated by multiconfigurational self-consistent field methods. The H(A) + H(B)...H(C) reaction (where r(BC) is adjusted to mimic the appropriate singlet-triplet energy gap) is shown to effectively model H-atom addition to benzyne. The doublet multiconfiguration wave functions are shown to mix the "singlet" and "triplet" valence bond structures of H(B)...H(C) along the reaction coordinate; however, the extent of mixing is dependent on the singlet-triplet energy gap (DeltaE(ST)) of the H(B)...H(C) diradical. Early in the reaction, the ground-state wave function is essentially the "singlet" VB function, yet it gains significant "triplet" VB character along the reaction coordinate that allows H(A)-H(B) bond formation. Conversely, the wave function of the first excited state is predominantly the "triplet" VB configuration early in the reaction coordinate, but gains "singlet" VB character when the H-atom is close to a radical center. As a result, the potential energy surface (PES) for H-atom addition to triplet H(B)...H(C) diradical is repulsive! The H3 model predicts, in agreement with the actual calculations on benzyne, that the singlet diradical electrons are not coupled strongly enough to give rise to an activation barrier associated with C-H bond formation. Moreover, this model predicts that the PES for H-atom addition to triplet benzyne will be characterized by a repulsive curve early in the reaction coordinate, followed by a potential avoided crossing with the (pi)1(sigma*)1 state of the phenyl radical. In contrast to H-atom addition, large activation barriers characterize the abstraction process in both the singlet ground state and first triplet state. In the ground state, this barrier results from the weakly avoided crossing of the dominant VB configurations in the ground-state singlet (S0) and first excited singlet (S1) because of the large energy gap between S0 and S1 early in the reaction coordinate. Because the S1 state is best described as the combination of the triplet X-H bond and the triplet H(B)...H(C) spin couplings, the activation barrier along the S0 abstraction PES will have much less dependence on the DeltaE(ST) of H(B)...H(C) than previously speculated. For similar reasons, the T1 potential surface is quite comparable to the S0 PES.     

On the intrinsic population of the lowest triplet state of thymine

The population of the lowest triplet state of thymine after near-UV irradiation has been established, on the basis of CASPT2//CASSCF quantum chemical calculations, to take place via three distinct intersystem crossing mechanisms from the initially populated singlet bright 1pipi* state. Two singlet-triplet crossings have been found along the minimum-energy path for ultrafast decay of the singlet state at 4.8 and 4.0 eV, involving the lowest 3npi* and 3pipi* states, respectively. Large spin-orbit coupling elements predict efficient intersystem crossing processes in both cases. Another mechanism involving energy transfer from the lowest 1npi* state with much larger spin-orbit coupling terms can also be proposed. The wavelength dependence measured for the triplet quantum yield of pyrimidine nucleobases is explained by the location and accessibility of the singlet-triplet crossing regions.     

Quantum chemical study of the mechanism of reaction between NH (X 3sigma-) and H2, H2O, and CO2 under combustion conditions

Reactions of ground-state NH (3sigma-) radicals with H2, H2O, and CO2 have been investigated quantum chemically, whereby the stationary points of the appropriate reaction potential energy surfaces, that is, reactants, products, intermediates, and transition states, have been identified at the G3//B3LYP level of theory. Reaction between NH and H2 takes place via a simple abstraction transition state, and the rate coefficient for this reaction as derived from the quantum chemical calculations, k(NH + H2) = (1.1 x 10(14)) exp(-20.9 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K, is found to be in good agreement with experiment. For reaction between triplet NH and H2O, no stable intermediates were located on the triplet reaction surface although several stable species were found on the singlet surface. No intersystem crossing seam between triplet NH + H2O and singlet HNO + H2 (the products of lowest energy) was found; hence there is no evidence to support the existence of a low-energy pathway to these products. A rate coefficient of k(NH + H2O) = (6.1 x 10(13)) exp(-32.8 kcal mol(-1)/RT) cm3 mol(-1) s(-1) between 1000 and 2000 K for the reaction NH (3sigma-) + H2O --> NH2 (2B) + OH (2pi) was derived from the quantum chemical results. The reverse rate coefficient, calculated via the equilibrium constant, is in agreement with values used in modeling the thermal de-NO(x) process. For the reaction between triplet NH and CO2, several stable intermediates on both triplet and singlet reaction surfaces were located. Although a pathway from triplet NH + CO2 to singlet HNO + CO involving intersystem crossing in an HN-CO2 adduct was discovered, no pathway of sufficiently low activation energy was discovered to compare with that found in an earlier experiment [Rohrig, M.; Wagner, H. G. Proc. Combust. Inst. 1994, 25, 993.].     

Dissociative addition of water to a main group 13 metal cluster: computational study of the reaction of gallium dimer with H2O

We have investigated the lowest triplet and singlet potential energy surfaces (PESs) for the reaction of Ga(2) dimer with water. Under thermal conditions, we predict formation of the triplet ground state addition complex Ga(2)···OH(2)((3)B(1)) involving Ga···O···Ga bridge interaction. At the coupled cluster CCSD(T)/AE (CCSD(T)/ECP) computational levels, Ga(2)···OH(2)((3)B(1)) is bound by 5.5 (5.7) kcal/mol with respect to the ground state reactants Ga(2)((3)Π(u)) + H(2)O. Identification of the addition complex is in agreement with the experimental evidence from matrix isolation infrared (IR) spectroscopy reported recently by Macrae and Downs. The located minimum energy crossing points (MECPs) between the triplet and singlet energy surfaces on the entrance channel of Ga(2) + H(2)O are not expected to be energetically accessible under the matrix conditions, consistent with the lack of occurrence of Ga(2) insertion into the O-H bond under such conditions. The computed energies and harmonic and anharmonic vibrational frequencies for the triplet and singlet Ga(2)(H)(OH) insertion isomers indicate the singlet double-bridged Ga(μ-H)(μ-OH)Ga isomer to be the most stable and support the experimental IR identification of this species. The energy barrier for elimination of H(2) from the second most stable singlet HGa(μ-OH)Ga insertion isomer found to be 13.9 (12.9) kcal/mol is also consistent with the available experimental data.     

Molecular dynamics and metadynamics simulations of [2 + 2] photocycloaddition

Triplet state mechanism of [2 + 2] photocycloaddition forming a cyclobutane ring from two ethylenes is investigated in the context of photocatalysis. High‐level ab initio calculations are combined with ab initio adiabatic molecular dynamics and ab initio metadynamics for rare events modeling. In a photocatalytic scheme, a reactant reaches the triplet state either via intersystem crossing (ISC) or triplet sensitization. The model system adopts a biradical structure, which represents energy intersection with the ground state. The system either completes cyclization or undergoes fragmentation into two olefinic units. The potential and free energy surfaces of the cyclobutane/ethylenes system are mapped with multireference approaches describing possible reaction pathways. To obtain a full picture of a double bond photoreactivity, ab initio adiabatic dynamical calculations were used to estimate reaction yields and to model the effects of excess energy. The potential use of density functional theory based approaches for [2 + 2] photocycloaddition was investigated for future simulations and design of realistic photocatalytic systems. Dynamical aspects of [2 + 2] photocycloaddition via a triplet state manifold are investigated by combining ab initio multireference methods and ab initio molecular dynamics and metadynamics. The reaction pathways are studied for a model system of two ethylenes forming a cyclobutane ring to provide a basis for further studies on design of photocatalytic systems.     

High-level ab initio studies of unimolecular dissociation of the ground-state N3 radical

A comprehensive study of the unimolecular dissociation of the N(3) radical on the ground doublet and excited quartet potential energy surfaces has been carried out with multireference single and double excitation configuration interaction and second-order multireference perturbation methods. Two forms of the N(3) radical have been located in the linear and cyclic region of the lowest doublet potential energy surface with an isomerization barrier of 62.2 kcal/mol above the linear N(3). Three equivalent C(2v) minima of cyclic N(3) are connected by low barrier, meaning the molecule is free to undergo pseudorotation. The cyclic N(3) is metastable with respect to ground state products, N((4)S)+N(2), and dissociation must occur via intersystem crossing to a quartet potential energy surface. Minima on the seams of crossing between the doublet and quartet potential surfaces are found to lie substantially higher in energy than the cyclic N(3) minima. This strongly suggests that cyclic N(3) possesses a long collision-free lifetime even if formed with substantial internal excitation.     

Dynamical calculations of charge-transfer-to-solvent excited states of small I- (CH3CN)n clusters

Relaxation dynamics of photoexcited charge-transfer-to-solvent (CTTS) states for the I(-)(CH(3)CN)(n) (n = 2 and 3) clusters has been theoretically studied using electronic structure methods. First, we have calculated several lowest singlet and triplet potential energy surfaces using the multireference configuration interaction method. It was found that the character of the singlet CTTS excited-state potential surfaces is very similar to that of the triplet CTTS states. Due to a small singlet-triplet splitting, the lowest triplet potential energy surface was used as a good model to understand the dynamics of the photoexcited singlet CTTS states. We have carried out direct molecular dynamics simulations on the lowest triplet surface at the B3LYP level. When an I(-) anion is exteriorly solvated by CH(3)CN molecules, we found that the (CH(3)CN)(n)(-) anion cluster is effectively produced. In addition, when the I(-) anion is placed in the interior in I(-)(CH(3)CN)(n) clusters, photoexcitation gives an acetonitrile monomer anion plus neutral monomers. However, if the initial geometric configuration is distorted from the minimum structure, we also found that the (CH(3)CN)(2)(-) anion cluster, where an excess electron is internally trapped, is formed via I(-)(CH(3)CN)(2) + hnu --> I + (CH(3)CN)(2)(-) process.     

Ab initio quantum chemical studies on the reactions of CF3O2 with OH

The potential energy surface for the CF₃O₂ + OH reaction has been theoretically investigated using the DFT (B3LYP/6-311G(d,p)) level of theory. Both singlet and triplet potential energy surfaces are investigated. The reaction mechanism on the triplet surface is simple. However, the reaction mechanism on the singlet surface is more complicated. It is revealed that the formation of CF₃O + HO₂ is the dominant channel on the triplet surface. The potential energy surface (PES) for this reaction has been given according to the relative energies calculated at the DFT/B3LYP/6-311G(d,p) level. Because this reaction involves both triplet and singlet states, triplet–singlet intersystem crossing (ISC) crossing also have been investigated in this paper.     

Multireference configuration interaction studies on the ground and excited states of N2O2: the potential energy curves of N2O2 along N-N distance

In this paper, the ground and excited states of N2O2 were studied at the multireference configuration interaction (MRCI) level of theory with Dunning's [J. Chem. Phys. 90, 1007 (1985); 96, 6796 (1992)] correlation consistent basis sets augo-cc-pVDZ and aug-cc-pVTZ. The geometry optimizations were performed for the ground state of N2O2. The vertical excitation energies and transition moments were calculated for the low-lying singlet states of N2O2 including the lowest three 1A1 states, two 1B1 states, one 1B2 state, and two 1A2 states at the MRCI level of theory with Dunning's correlation consistent basis sets aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ. Furthermore, for the first time, the potential energy curves were calculated at the complete active space self-consistent-field and MRCI levels of theory for as many as 12 N2O2 singlet electronic states along the N-N distance. The dissociation asymptotes of these 12 N2O2 singlet electronic states were discussed.