Theoretical study of the reaction of ethane with oxygen molecules in the ground triplet and singlet delta States

Quantum chemical calculations are carried out to study the reaction of ethane with molecular oxygen in the ground triplet and singlet delta states. Transition states, intermediates, and possible products of the reaction on the triplet and singlet potential energy surfaces are identified on the basis of the coupled-cluster method. The basis set dependence of coupled-cluster energy values is estimated by the second-order perturbation theory. The values of energy barriers are also refined by using the compound CBS-Q and G3 techniques. It was found that the C(2)H(6) + O(2)(X(3)Σ(g)(-)) reaction leads to the formation of C(2)H(5) and HO(2) products, whereas the C(2)H(6) + O(2)(a(1)Δ(g)) process produces C(2)H(4) and H(2)O(2) molecules. The appropriate rate constants of these reaction paths are estimated on the basis of variational and nonvariational transition-state theories assuming tunneling and possible nonadiabatic transitions in the temperature range 500-4000 K. The calculations showed that the rate constant of the C(2)H(6) + O(2)(a(1)Δ(g)) reaction path is much greater than that of the C(2)H(6) + O(2)(X(3)Σ(g)(-)) one. At the same time, the singlet and triplet potential surface intersection is detected that leads to the appearance of the nonadiabatic quenching channel O(2)(a(1)Δ(g)) + C(2)H(6) → O(2)(X (3)Σ(g)(-)) + C(2)H(6). The rate constant of this process is estimated with the use of the Landau-Zener model. It is demonstrated that, in the case of the existence of thermal equilibrium in the distribution of molecules over the electronic states, at low temperatures (T < 1200 K) the main products of the reaction of C(2)H(6) with O(2) are C(2)H(4) and H(2)O(2), rather than C(2)H(5) and HO(2). At higher temperature (T > 1200 K) the situation is inverted.     

The Ni + O2 reaction: a combined IR matrix isolation and theoretical study of the formation and structure of NiO2

The reaction of Ni atoms with molecular oxygen has been reinvestigated experimentally in neon matrices and theoretically at the DFT PW91PW91/6311G(3df) level. Experimental results show that i) the nature of the ground electronic state of the superoxide metastable product is the same in neon and argon matrices, ii) two different photochemical pathways exist for the conversion of the superoxide to the dioxide ground state (involving 1.6 or 4 eV photons) and iii) an important matrix effect exists in the Ni + O(2)--> Ni(O(2)) or ONiO branching ratios. Theoretical results confirm that the electronic ground state of the metastable superoxide corresponds to the singlet state, in agreement with former CCSD(T) calculations, but in contradiction with other recent works. Our results show that the ground electronic state of the dioxide is (1)Sigma(+)(g) with the lowest triplet and quintet states at slightly higher energy, consistent with the observation of weak vibronic transitions in the near infrared. The potential energy profiles are modelled for the ground state and nine electronic excited states and a pathway for the Ni(triplet) + O(2)(triplet) --> Ni(O(2)) or ONiO (singlet) reaction is proposed, as well as for the Ni(O(2)) --> ONiO photochemical reaction, accounting for the experimental observations.     

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.].     

Experimental and theoretical investigations of the inelastic and reactive scattering dynamics of O(3p) + D2

This paper presents a combined experimental and theoretical study of the dynamics of O((3)P) + D(2) collisions, with emphasis on a center-of-mass (c.m.) collision energy of 25 kcal mol(-1). The experiments were conducted with a crossed-molecular-beams apparatus, employing a laser detonation source to produce hyperthermal atomic oxygen and mass spectrometric detection to measure the product angular and time-of-flight distributions. The novel beam source, which enabled these experiments to be conducted, contributed unique challenges to the experiments and to the analysis, so the experimental methods and approach to the analysis are discussed in detail. Three different levels of theory were used: (1) quasiclassical trajectories (QCT), (2) time-independent quantum scattering calculations based on high-quality potential surfaces for the two lower-energy triplet states, and (3) trajectory-surface-hopping (TSH) studies that couple the triplet surfaces with the lowest singlet surface using a spin-orbit Hamiltonian derived from ab-initio calculations. The latter calculations explore the importance of intersystem crossing in the dynamics. Both experiment and theory show that inelastically scattered O atoms scatter almost exclusively in the forward direction, with little or no loss of translational energy. For the reaction, O((3)P) + D(2) --> OD + D, the experiment shows that, on average, approximately 50% of the available energy goes into product translation and that the OD product angular distributions are largely backward-peaked. These results may be interpreted in light of the QCT and TSH calculations, leading to the conclusion that the reaction occurs mainly on triplet potential energy surfaces with, at most, minor intersystem crossing to a singlet surface. Reaction on either of the two low-lying reactive triplet surfaces proceeds through a rebound mechanism in which the angular distributions are backward-peaked and the product OD is both vibrationally and rotationally excited. The quantum scattering results are in good agreement with QCT calculations, indicating that quantum effects are relatively small for this reaction at a collision energy of 25 kcal mol(-1).     

Experimental and theoretical studies of the C2F4 + O reaction: nonadiabatic reaction mechanism

In this work, the C(2)F(4)(X(1)A(g)) + O((3)P) reaction was investigated experimentally using molecular beam-threshold ionization mass spectrometry (MB-TIMS). The major primary products were observed to be CF(2)O (+ CF(2)) and CF(3) (+ CFO), with measured approximate yields of % versus %, respectively, neglecting minor products. Furthermore, the lowest-lying triplet and singlet potential energy surfaces for this reaction were constructed theoretically using B3LYP, G2M(UCC, MP2), CBS-QB3, and G3 methods in combination with various basis sets such as 6-31G(d), 6-311+G(3df), and cc-pVDZ. The primary product distribution for the multiwell multichannel reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis. It was found that the observed production of CF(3) (+ CFO) can only occur on the singlet surface, in parallel with formation of ca. 5 times more CF(2)O(X) + CF(2)(X(1)A(1)). This requires fast intersystem crossing (ISC) from the triplet to the singlet surface at a rate of ca. 4 x 10(12) s(-1). The theoretical calculations combined with the experimental results thus indicate that the yield of triplet CF(2)(?(3)B(1)) + CF(2)O formed on the triplet surface prior to ISC is < or =35%, whereas singlet CF(2)(X(1)A(1)) + CF(2)O is produced with yield > or =60%, after ISC. In addition, the thermal rate coefficients k(O + C(2)F(4)) in the T = 150-1500 K range were computed using multistate transition state theory and can be expressed as k(T) = 1.67 x 10(-16) x T(1.48) cm(3) molecule(-1) s(-1); they are in agreement with the available experimental results in the T = 298-500 K range.     

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.     

Computational study of carbon atom (3P and 1D) reaction with CH2O. Theoretical evaluation of 1B1 methylene production by C (1D)

Singlet and triplet free energy surfaces for the reactions of C atom ((3)P and (1)D) with CH(2)O are studied computationally to evaluate the excited singlet ((1)B(1)) methylene formation from deoxygenation of CH(2)O by C ((1)D) atom as suggested by Shevlin et al. Carbon atoms can react by addition to the oxygen lone pair or to the C=O double bond on both the triplet and singlet surfaces. Triplet C ((3)P) atoms will deoxygenate to give CO plus CH(2) ((3)B(1)) as the major products, while singlet C ((1)D) reactions will form ketene and CO plus CH(2) ((1)A(1)). No definitive evidence of the formation of excited singlet ((1)B(1)) methylene was found on the singlet free energy surface. A conical intersection between the (1)A' and (1)A' ' surfaces located near an exit channel may play a role in product formation. The suggested (1)B(1) state of methylene may form via the (1)A' ' surface only if dynamic effects are important. In an effort to interpret experimental observation of products trapped by (Z)-2-butene, formation of cis- and trans-1,2-dimethylcyclopropane is studied computationally. The results suggests that "hot" ketene may react with (Z)-2-butene nonstereospecifically.     

O(3P) + CO2 collisions at hyperthermal energies: dynamics of nonreactive scattering, oxygen isotope exchange, and oxygen-atom abstraction

The dynamics of O((3)P) + CO(2) collisions at hyperthermal energies were investigated experimentally and theoretically. Crossed-molecular-beams experiments at = 98.8 kcal mol(-1) were performed with isotopically labeled (12)C(18)O(2) to distinguish products of nonreactive scattering from those of reactive scattering. The following product channels were observed: elastic and inelastic scattering ((16)O((3)P) + (12)C(18)O(2)), isotope exchange ((18)O + (16)O(12)C(18)O), and oxygen-atom abstraction ((18)O(16)O + (12)C(18)O). Stationary points on the two lowest triplet potential energy surfaces of the O((3)P) + CO(2) system were characterized at the CCSD(T)/aug-cc-pVTZ level of theory and by means of W4 theory, which represents an approximation to the relativistic basis set limit, full-configuration-interaction (FCI) energy. The calculations predict a planar CO(3)(C(2v), (3)A') intermediate that lies 16.3 kcal mol(-1) (W4 FCI excluding zero point energy) above reactants and is approached by a C(2v) transition state with energy 24.08 kcal mol(-1). Quasi-classical trajectory (QCT) calculations with collision energies in the range 23-150 kcal mol(-1) were performed at the B3LYP/6-311G(d) and BMK/6-311G(d) levels. Both reactive channels observed in the experiment were predicted by these calculations. In the isotope exchange reaction, the experimental center-of-mass (c.m.) angular distribution, T(θ(c.m.)), of the (16)O(12)C(18)O products peaked along the initial CO(2) direction (backward relative to the direction of the reagent O atoms), with a smaller isotropic component. The product translational energy distribution, P(E(T)), had a relatively low average of = 35 kcal mol(-1), indicating that the (16)O(12)C(18)O products were formed with substantial internal energy. The QCT calculations give c.m. P(E(T)) and T(θ(c.m.)) distributions and a relative product yield that agree qualitatively with the experimental results, and the trajectories indicate that exchange occurs through a short-lived CO(3)* intermediate. A low yield for the abstraction reaction was seen in both the experiment and the theory. Experimentally, a fast and weak (16)O(18)O product signal from an abstraction reaction was observed, which could only be detected in the forward direction. A small number of QCT trajectories leading to abstraction were observed to occur primarily via a transient CO(3) intermediate, albeit only at high collision energies (149 kcal mol(-1)). The oxygen isotope exchange mechanism for CO(2) in collisions with ground state O atoms is a newly discovered pathway through which oxygen isotopes may be cycled in the upper atmosphere, where O((3)P) atoms with hyperthermal translational energies can be generated by photodissociation of O(3) and O(2).     

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.     

Potential energy surfaces, product distributions and thermal rate coefficients of the reaction of O(3P) with C2H4(X1Ag): a comprehensive theoretical study

The potential energy surface for the O((3)P) + C(2)H(4) reaction, which plays an important role in C(2)H(4)/O(2) flames and in hydrocarbon combustion in general, was theoretically reinvestigated using various quantum chemical methods, including G3, CBS-QB3, G2M(CC,MP2), and MRCI. The energy surfaces of both the lowest-lying triplet and singlet electronic states were constructed. The primary product distribution for the multiwell multichannel reaction was then determined by RRKM statistical rate theory and weak-collision master equation analysis using the exact stochastic simulation method. Intersystem crossing of the "hot" CH(2)CH(2)O triplet adduct to the singlet surface, shown to account for about half of the products, was estimated to proceed at a rate of approximately 1.5 x 10(11) s(-1). In addition, the thermal rate coefficients k(O + C(2)H(4)) in the T = 200-2000 K range were computed using multistate transition state theory and fitted by a modified Arrhenius expression as k(T) = 1.69 x 10(-16) x T(1.66) x exp(-331 K/T) . Our computed rates and product distributions agree well with the available experimental results. Product yields are found to show a monotonic dependence on temperature. The major products (with predicted yields at T = 300 K/2000 K) are: CH(3) + CHO (48/37%), H + CH(2)CHO (40/19%), and CH(2)(X(3)B(1)) + H(2)CO (5/29%), whereas H + CH(3)CO, H(2) + H(2)CCO, and CH(4) + CO are all minor (< or =5%).     

A new mechanism for the production of highly vibrationally excited OH in the mesosphere: an ab initio study of the reactions of O2(A 3Sigmau+ and A' 3Deltau)+H

In an attempt to explain the observed nightglow emission from OH(v=10) in the mesosphere that has the energy greater than the exothermicity of the H+O(3) reaction, potential energy surfaces were calculated for reactions of high lying electronic states of O(2)(A (3)Sigma(u) (+) and A' (3)Delta(u)) with atomic hydrogen H((2)S) to produce the ground state products OH((2)Pi)+O((3)P). From collinear two-dimensional scans, several adiabatic and nonadiabatic pathways have been identified. Multiconfigurational single and double excitation configuration interaction calculations show that the adiabatic pathways on a (4)Delta potential surface from O(2)(A' (3)Delta)+H and a (4)Sigma(+) potential surface from O(2)(A (3)Sigma(u) (+))+H are the most favorable, with the zero-point corrected barrier heights of as low as 0.191 and 0.182 eV, respectively, and the reactions are fast. The transition states for these pathways are collinear and early, and the reaction coordinate suggests that the potential energy release of ca. 3.8 eV (larger than the energy required to excite OH to v=10) is likely to favor high vibrational excitation.     

Oxidation of CO by SO2: a theoretical study

The elementary reaction SO(2) + CO --> CO(2) + SO((3)Sigma) (1) and the subsequent reaction SO((3)Sigma) + CO --> CO(2) + S((3)P) (2) have been studied by the application of the Gaussian-3//B3LYP quantum chemical approach to characterize the potential energy surfaces and transition state kinetic analysis to derive rate coefficients. Reaction 1 is found to take place via two transition states (TS), a cis-OSOCO TS and a trans-OSOCO TS. Reaction via the cis-TS is concerted and takes place on a singlet surface. Intersystem crossing to the final products occurs after passage through the barrier on the singlet surface. The trans-TS leads to a very weakly bound singlet OSOCO intermediate that then passes through a second TS (on the triplet surface) to form the products. Reaction 2 takes place on triplet surfaces. There is a concerted reaction through a cis-SOCO TS and a weakly bound trans-SOCO has also been identified. Reaction 2 is analogous to the reaction CO + O(2)((3)Sigma) --> CO(2) + O((3)P) (3), and this reaction has been reinvestigated at a similar level of theory and the rate coefficient derived by quantum chemistry is compared with experiment. The sensitive effects of trace impurities such as H(2), H(2)O, and hydrocarbons on the accurate experimental determination of the rate coefficient of reaction 3 is discussed. Using rate coefficients for reactions 1 and 2 obtained via quantum chemical calculations, we have been unable to model the extent of decomposition of SO(2) measured in a shock tube study of reaction between SO(2) and CO [Bauer, S. H.; Jeffers, P.; Lifshitz, A.; Yadava, B. P. Proc. Combust. Inst. 1971, 13, 417]. In light of the known sensitivity of reaction 3 to trace impurities, we have incorporated trace amounts of H(2), CH(4), or H(2)O, together with our rate coefficients for (1) and (2), in a kinetic model of Alzueta et al. [Combust. Flame 2001, 127, 2234], which is then shown to be able to substantially model the SO(2) data of Bauer et al. In the course of this modeling study we also computed heats of formation for a number of sulfur-containing small molecules: HS, HSO, HSOH, HOSO, HS(2), HSO(2), HOSO(2), HOSOH, and HOSHO.     

Photoelectron spectroscopic study of the oxyallyl diradical

The photoelectron spectrum of the oxyallyl (OXA) radical anion has been measured. The radical anion has been generated in the reaction of the atomic oxygen radical anion (O(?-)) with acetone. Three low-lying electronic states of OXA have been observed in the spectrum. Electronic structure calculations have been performed for the triplet states ((3)B(2) and (3)B(1)) of OXA and the ground doublet state ((2)A(2)) of the radical anion using density functional theory (DFT). Spectral simulations have been carried out for the triplet states based on the results of the DFT calculations. The simulation identifies a vibrational progression of the CCC bending mode of the (3)B(2) state of OXA in the lower electron binding energy (eBE) portion of the spectrum. On top of the (3)B(2) feature, however, the experimental spectrum exhibits additional photoelectron peaks whose angular distribution is distinct from that for the vibronic peaks of the (3)B(2) state. Complete active space self-consistent field (CASSCF) method and second-order perturbation theory based on the CASSCF wave function (CASPT2) have been employed to study the lowest singlet state ((1)A(1)) of OXA. The simulation based on the results of these electronic structure calculations establishes that the overlapping peaks represent the vibrational ground level of the (1)A(1) state and its vibrational progression of the CO stretching mode. The (1)A(1) state is the lowest electronic state of OXA, and the electron affinity (EA) of OXA is 1.940 ± 0.010 eV. The (3)B(2) state is the first excited state with an electronic term energy of 55 ± 2 meV. The widths of the vibronic peaks of the X? (1)A(1) state are much broader than those of the a? (3)B(2) state, implying that the (1)A(1) state is indeed a transition state. The CASSCF and CASPT2 calculations suggest that the (1)A(1) state is at a potential maximum along the nuclear coordinate representing disrotatory motion of the two methylene groups, which leads to three-membered-ring formation, i.e., cyclopropanone. The simulation of b? (3)B(1) OXA reproduces the higher eBE portion of the spectrum very well. The term energy of the (3)B(1) state is 0.883 ± 0.012 eV. Photoelectron spectroscopic measurements have also been conducted for the other ion products of the O(?-) reaction with acetone. The photoelectron imaging spectrum of the acetylcarbene (AC) radical anion exhibits a broad, structureless feature, which is assigned to the X? (3)A' state of AC. The ground ((2)A') and first excited ((2)A') states of the 1-methylvinoxy (1-MVO) radical have been observed in the photoelectron spectrum of the 1-MVO ion, and their vibronic structure has been analyzed.     

Reaction mechanism of the CCN radical with nitric oxide

To investigate the possibility of the carbyne radical CCN in removal of nitric oxide, a detailed computational study is performed at the Gaussian-3//B3LYP/6-31G(d) level on the CCN + NO reaction by constructing the singlet and triplet electronic state [C(2)N(2)O] potential energy surfaces (PESs). The barrierless formation of the chain-like isomers NCCNO (singlet at -106.5, triplet cis at -48.2 and triplet trans at -47.6 kcal/mol) is the most favorable entrance attack on both singlet and triplet PESs. Subsequently, the singlet NCCNO takes an O-transfer to form the branched intermediate singlet NCC(O)N (-85.6), which can lead to the fragments CN + NCO (-51.2) via the intermediate singlet NCOCN (-120.3). The simpler evolution of the triplet NCCNO is the direct N-O rupture to form the weakly bound complex triplet NCCN...O (-56.2) before the final fragmentation to NCCN + (3)O (-53.5). However, the lower lying products (3)NCN + CO (-105.6) and (3)CNN + CO (-74.6) are kinetically much less competitive. All the involved transition states for generation of CN + NCO and NCCN + (3)O lie much lower than the reactants. Thus, the novel reaction CCN + NO can proceed effectively even at low temperatures and is expected to play a role in both combustion and interstellar processes. Significant differences are found on the singlet PES between the CCN + NO and CH + NO reaction mechanisms.     

Theoretical study of the reaction of Cu with OCS

The reactivity of Cu⁺ with OCS on both singlet and triplet potential energy surfaces (PES) has been investigated at the UB3LYP/6-311+G(d) level. The object of this investigation was the elucidation of the reaction mechanism. The calculated results indicated that both the C–S and C–O bond activations proceed via an insertion–elimination mechanism. Intersystem crossing between the singlet and triplet surfaces may occur along both the C–S and C–O bond activation branches. The ground states of CuS⁺ and CuO⁺ were found to be triplets, whereas CuCO⁺ and CuCS⁺ have singlet ground states. The C–S bond activation is energetically much more favorable than the C–O bond activation. All theoretical results are in line with early experiments.     

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.     

The almost bottleable triplet carbene: 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene.

Computations on 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene (1) using ab initio and density functional theory methods underscore the unusual stability of the triplet over the singlet state. At the B3LYP/6-311G(d,p) level, the triplet state had a slightly bent central C-C-C bond angle of 167 degrees, whereas this angle in the singlet was 134 degrees. The B3LYP singlet-triplet splitting (12.2 kcal/mol) was larger than that of the parent molecule (5.8 kcal/mol), diphenylcarbene (2), which also has a triplet ground state. The energy of a suitable isodesmic reaction showed the triplet and singlet states of (1) to be destabilized, by 6.3 and 12.5 kcal/mol, respectively, due to the combined effects of the CF3, Br, and alkyl substituents. The linear-coplanar form of (3)(1), which might facilitate dimerization or electrophilic attack at the more exposed diradical center, was prohibitively (35.9 kcal/mol) higher in energy. Our results confirm Tomioka's conclusion that the triplet diarylcarbene, ortho-substituted with bulky CF3 and Br substituents, is persistent due to steric protection of the diradical center. Dimerization and other possible reaction pathways are inhibited, not only by the bulky ortho substituents but also by the para alkyl groups. The increase in stability of the triplet ((3)(1)) state relative to the singlet ((1)(1)) state does not influence the reactivity directly.     

Global potential energy surfaces for O((3)P) + H2O((1)A1) collisions

Global analytic potential energy surfaces for O((3)P) + H(2)O((1)A(1)) collisions, including the OH + OH hydrogen abstraction and H + OOH hydrogen elimination channels, are presented. Ab initio electronic structure calculations were performed at the CASSCF + MP2 level with an O(4s3p2d1f)/H(3s2p) one electron basis set. Approximately 10(5) geometries were used to fit the three lowest triplet adiabatic states corresponding to the triply degenerate O((3)P) + H(2)O((1)A(1)) reactants. Transition state theory rate constant and total cross section calculations using classical trajectories to collision energies up to 120?kcal mol(-1) (～11?km s(-1) collision velocity) were performed and show good agreement with experimental data. Flux-velocity contour maps are presented at selected energies for H(2)O collisional excitation, OH + OH, and H + OOH channels to further investigate the dynamics, especially the competition and distinct dynamics of the two reactive channels. There are large differences in the contributions of each of the triplet surfaces to the reactive channels, especially at higher energies. The present surfaces should support quantitative modeling of O((3)P) + H(2)O((1)A(1)) collision processes up to ～150?kcal mol(-1).     

Unusual enlargment of magnetic interactions in the lowest excited states of nitroxide radical chromium(III) complexes as revealed by magnetic circular dichroism and luminescence

Magnetic circular dichroism and NIR luminescence of nitroxide radical complexes, [Cr(III)(beta-diketonato)(2)(NIT2py or IM2py)]PF(6), demonstrate that the energy gaps between the singlet ((1)L(D)) and triplet ((3)L(D)) spin coupled levels in the lowest excited (2)E(g),(2)T(1g) states are much larger than those in the ground state. This is the first observation of magnetic interactions in the excited states of radical complexes, which could be elucidated in terms of the exchange mechanism.     

Donor stabilized borylnitrene: a highly reactive BN analogue of vinylidene

A donor atom stabilized borylnitrene, 2-nitreno-1,3,2-benzodioxaborole 4c, is characterized by matrix isolation IR, UV, and ESR spectroscopy as well as multiconfiguration SCF and CI computations. UV irradiation (lambda = 254 nm) of the corresponding azide 6c, isolated in solid argon at 10 K, produces 4c in high yield. The oxygen donor atoms in 4c result in a triplet ground state (|D/hc| = 1.492 cm(-)(1), |E/hc| = 0.004 cm(-)(1)) for the borylnitrene. The lowest energy singlet state ((1)A(1)) is 33 kcal mol(-)(1) higher in energy and closely related to the ground state of vinylidene. Under the conditions of matrix isolation, triplet 4c is photochemically and thermally stable toward rearrangement to the corresponding cyclic iminoborane. Photochemical irradiation (lambda > 550 nm) of 4c rather causes an efficient reaction with molecular nitrogen, lying in matrix sites nearby, to give 6c. Similarly, photochemical, but not thermal, trapping of 4c with CO is possible and results in the corresponding isocycanate 9c. Thermal reaction of 4c with O(2) in doped argon matrixes at 35 K could be observed by IR spectroscopy to result in borylnitroso-O-oxide 17c as shown by (18)O(2) labeling experiments and DFT computations. The diradical 17c is very photolabile and quickly rearranges to the nitritoborane 16c upon irradiation.