Ab initio chemical kinetics for reactions of ClO with Cl2O2 isomers

The mechanisms for the reactions of ClO with ClOClO, ClOOCl, and ClClO(2) have been investigated at the CCSD(T)/6-311+G(3df)//PW91PW91∕6-311+G(3df) level of theory. The rate constants for their low energy channels have been calculated by statistical theory. The results show that the main products for the reaction of ClO with ClOClO are ClOCl + ClOO, which can be produced readily by ClO abstracting the terminal O atom from ClOClO. This process occurs without an intrinsic barrier, with the predicted rate constant: k (ClO + ClOClO) = 7.26 × 10(-10) T(-0.15) × exp (-40/T) cm(3)molecule(-1)s(-1) for 200-1500 K. For the reactions of ClO + ClOOCl and ClClO(2), the lowest abstraction barriers are 7.2 and 7.3 kcal/mol, respectively, suggesting that these two reactions are kinetically unimportant in the Earth's stratosphere as their rate constants are less than 10(-14) cm(3)molecule(-1)s(-1) below 700 K. At T = 200-1500 K, the computed rate constants can be represented by k (ClO+ ClOOCl) = 1.11 × 10 (-14) T (0.87) exp (-3576/T) and k (ClO+ ClClO(2)) = 4.61 × 10(-14) T(0.53) exp (-3588/T) cm(3)molecule(-1)s(-1). For these systems, no experimental or theoretical kinetic data are available for comparison.     

Density functional study of chlorine–oxygen compounds related to the ClO self-reaction

Geometrical parameters, vibrational frequencies, relative stabilities, and dissociation energies of the three stable Cl₂O₂ isomers and the OClO and ClOO radicals were investigated by density functional theory (DFT). The present analysis shows that DFT using hybrid functionals is capable of describing these systems to at least the same degree of accuracy as ab initio methods. The average absolute bond-length deviation of ClClO₂, ClOOCl, and ClO₂ from experimental results is 0.024/0.027 Å, with a maximum deviation for the dichlorine peroxide O(SINGLE BOND)O bond equal to 0.072/0.063 Å, for the B3PW91 and B3LYP functionals, respectively. The average absolute bond-angle deviation for the hybrid functionals is 0.8°. Harmonic vibrational frequencies calculated with DFT give for all Cl(SINGLE BOND)O compounds good agreement with experiments. The dissociation energies of ClOOCl, OClO, and ClOO were found to be in good agreement with experiments, the average error being less than 1.2 kcal/mol. The two isomers chloryl chloride (ClClO₂) and dichlorine peroxide (ClOOCl) were found to be approximately 9 kcal/mol more stable than the chlorine chlorite (ClOClO) isomer. The ClOO isomer is predicted to be 3.0 kcal/mol more stable than OClO, in accordance with the experimental value of 4 kcal/mol. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 203–217, 1998     

Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.     

Kinetics of the OH + ClOOCl and OH + Cl2O reactions: experiment and theory

The rate coefficients for the reactions OH + ClOOCl --> HOCl + ClOO (eq 5) and OH + Cl2O --> HOCl + ClO (eq 6) were measured using a fast flow reactor coupled with molecular beam quadrupole mass spectrometry. OH was detected using resonance fluorescence at 309 nm. The measured Arrhenius expressions for these reactions are k5 = (6.0 +/- 3.5) x 10(-13) exp((670 +/- 230)/T) cm(3) molecule(-1) s(-1) and k6 = (5.1 +/- 1.5) x 10(-12) exp((100 +/- 92)/T) cm(3) molecule(-1) s(-1), respectively, where the uncertainties are reported at the 2sigma level. Investigation of the OH + ClOOCl potential energy surface using high level ab initio calculations indicates that the reaction occurs via a chlorine abstraction from ClOOCl by the OH radical. The lowest energy pathway is calculated to proceed through a weak ClOOCl-OH prereactive complex that is bound by 2.6 kcal mol(-1) and leads to ClOO and HOCl products. The transition state to product formation is calculated to be 0.59 kcal mol(-1) above the reactant energy level. Inclusion of the OH + ClOOCl rate data into an atmospheric model indicates that this reaction contributes more than 15% to ClOOCl loss during twilight conditions in the Arctic stratosphere. Reducing the rate of ClOOCl photolysis, as indicated by a recent re-examination of the ClOOCl UV absorption spectrum, increases the contribution of the OH + ClOOCl reaction to polar stratospheric loss of ClOOCl.     

Accurate heats of formation of the "Arduengo-type" carbene and various adducts including H2 from ab initio molecular orbital theory

The heats of formation of saturated and unsaturated diaminocarbenes (imadazol(in)-2-ylidenes) have been calculated by using high levels of ab initio electronic structure theory. The calculations were done at the coupled cluster level through noniterative triple excitations with augmented correlation consistent basis sets up through quadruple. In addition, four other corrections were applied to the frozen core atomization energies: (1) a zero point vibrational correction; (2) a core/valence correlation correction; (3) a scalar relativistic correction; (4) a first-order atomic spin-orbit correction. The value of DeltaHf( 298) for the unsaturated carbene 1 is calculated to be 56.4 kcal/mol. The value of DeltaHf( 298) for the unsaturated triplet carbene (3)1 is calculated to be 142.8 kcal/mol, giving a singlet-triplet splitting of 86.4 kcal/mol. Addition of a proton to 1 forms 3 with DeltaHf( 298)(3) = 171.6 kcal/mol with a proton affinity for 1 of 250.5 kcal/mol at 298 K. Addition of a hydrogen atom to 1 forms 4 with DeltaHf( 298)(4) = 72.7 kcal/mol and a C-H bond energy of 35.8 kcal/mol at 298 K. Addition of H- to 1 gives 5 with DeltaHf( 298)(5) = 81.2 kcal/mol and 5 is not stable with respect to loss of an electron to form 4. Addition of H2 to the carbene center forms 6 with DeltaHf( 298)(6) = 41.5 kcal/mol and a heat of hydrogenation at 298 K of -14.9 kcal/mol. The value of DeltaHf( 298) for the saturated carbene 7 (obtained by adding H2 to the C=C bond of 1) is 47.4 kcal/mol. Hydrogenation of 7 to form the fully saturated imidazolidine, 8, gives DeltaHf( 298)(8) = 14.8 kcal/mol and a heat of hydrogenation at 298 K of -32.6 kcal/mol. The estimated error bars for the calculated heats of formation are +/-1.0 kcal/mol.     

Thermochemical parameters of CHFO and CF2O

Thermochemical properties of CHFO and CF 2O and their derivatives were calculated by using coupled-cluster theory (U)CCSD(T) calculations with the aug-cc-pV nZ ( n = D, T, Q, 5) basis sets extrapolated to the complete basis set limit with additional corrections. The predicted properties include the following. Enthalpies of formation (298 K, kcal/mol): Delta H f (CF 2O) = -144.7, Delta H f(CHFO) = -91.1, Delta H f (CFO (*)) = -41.6. Bond dissociation energy (0 K, kcal/mol): BDE(CFO-F) = 120.7, BDE(CHO-F) = 119.1, BDE(CFO-H) = 100.2. Ionization potential (eV): IP 1(CF 2O) = 13.04, IP 2(CF 2O) = 14.09, IP 1(CHFO) = 12.41, IP 2(CHFO) = 13.99, IP 1(CFO (*)) = 9.34. Proton affinity (298 K, kcal/mol), PA O(CF 2O) = 148.8, PA O(CHFO) = 156.7, PA F(CHFO) = 154.5 kcal/mol. Electron affinity: EA(CFO (*)) = 2.38 eV. Triplet-singlet separation gap (eV): Delta E T1-S0(CF 2O) = 4.47, Delta E T1-S0(CHFO) = 4.36. Triplet-triplet transition energy (eV): Delta E T2-T1(CF 2O) = 0.44. The new calculated values contribute to solving some persistent discrepancies in the literature. The effects of F-atoms on thermochemical parameters are not linearly additive, and the changes are largely dominated by the first F-substitution. On the basis of the calculated proton affinities of CF 2O and CF 3OH, the nucleophilicities of the oxygen atoms are, within computational errors, the same in both compounds.     

An ab initio study of the electronic structure of diimide

Ab initio calculations on the structure and geometry of the three isomers of N₂H₂ (trans-diimide, cis-diimide, and 1,1-dihydrodiazine) were performed both on HF and CI level using gaussian basis sets with polarization functions. The trans and cis isomers have singlet ground states; the trans isomer is found to be lower in energy than the cis isomer by 6.9 kcal/mol (HF) and 5.8 kcal/mol (CI), respectively. The barrier for the trans-cis isomerization is predicted to be 56 (HF) and 55 (CI) kcal/mol. H₂ N=N has a triplet ground state with a non-planar equilibrium geometry and a rather long NN bond of 1.34 Å. Its lowest singlet state, however, is planar with an NN double bond of 1.22 Å; it is found to lie about 3 kcal/mol above the triplet and 26 kcal/mol above the singlet ground state of trans-diimide.     

Is the 1,3,5-tridehydrobenzene triradical a cyclopropenyl radical analogue?

The gas-phase heat of formation (DeltaH(f,298)) of the 1,3,5-tridehydrobenzene triradical has been determined by using a negative ion thermochemical cycle. The first three measurements carried out were of the gas-phase acidity of 3,5-dichlorobenzoic acid, the enthalpy for decarboxylation of 3,5-dichlorobenzoate, and the enthalpy for chloride loss from 3,5,-dichlorophenide and constitute the measurement of the heat of formation for 5-chloro-m-benzyne. The last two measurements, the electron affinity of 5-chloro-m-benzyne, and the threshold for chloride loss from 5-chloro-m-benzyne, when combined with DeltaH(f,298) of 5-chloro-m-benzyne, give the heat of formation of the triradical. The 5-chloro-m-benzyne heat of formation is 116.2 +/- 3.7 kcal/mol. The heat of formation of the 1,3,5-tridehydrobenzene triradical measured in this work is 179.1 +/- 4.6 kcal/mol. This heat of formation was used to derive the bond dissociation energy (BDE) at the 5-position of m-benzyne, a third BDE in benzene. The BDE, at 109.2 +/- 5.6 kcal/mol, is ca. 4 kcal/mol lower than the first BDE in benzene (112.9 kcal/mol) and significantly higher than the BDE of phenyl radical at the meta position. The agreement between the first and third BDEs implies that the triradical is best described as a phenyl radical that interacts little with a m-benzyne moiety. The experimentally measured BDE is in good agreement with multireference configuration interaction calculations, which predict a (2)A(1) ground state for the Jahn-Teller distorted triradical. The trends in the first, second, and third BDEs of benzene are similar to those found for cyclopropane, suggesting a cyclopropenyl-like electronic structure within the six-membered ring of the 1,3,5-benzene triradical.     

Peculiar structure of the HOOO(-) anion

The HOOO(-) anion (1) can adopt a triplet state (T-1) or a singlet state (S-1), where the former is 9.8 kcal/mol (DeltaH(298) = 10.3 kcal/mol) more stable than the latter. S-1 possesses a strong O-OOH bond with some double bond character and a weakly covalent OO-OH bond (1.80 A) according to CCSD(T)/6-311++G(3df,3pd) calculations (the longest O-O bond ever found for a peroxide). In aqueous solution, S-1 adopts a geometry closely related to that of HOOOH (OO(O), 1.388 A; (O)OO(H), 1.509 A; tau(OOOH), 78.3 degrees ), justifying that S-1 is considered the anion of HOOOH. Dissociation into HO anion and O(2)((1)Delta(g)) requires 15.4 (DeltaH(298) = 14.3; DeltaG(298) = 8.9) kcal/mol. Structure T-1 corresponds to a van der Waals complex between HO anion and O(2)((3)Sigma(g)(-)) having a binding energy of 2.7 (DeltaH(298) = 2.1) kcal/mol. Modes of generating S-1 in aqueous solution are discussed, and it is shown that S-1 represents an important intermediate in ozonation reactions.     

Photolysis of (3-methyl-2H-azirin-2-yl)-phenylmethanone: direct detection of a triplet vinylnitrene intermediate

The photoreactivity of (3-methyl-2H-azirin-2-yl)-phenylmethanone, 1, is wavelength-dependent (Singh et al. J. Am. Chem. Soc. 1972, 94, 1199-1206). Irradiation at short wavelengths yields 2P, whereas longer wavelengths produce 3P. Laser flash photolysis of 1 in acetonitrile using a 355 nm laser forms its triplet ketone (T(1K), broad absorption with λ(max) ~ 390-410 nm, τ ~ 90 ns), which cleaves and yields triplet vinylnitrene 3 (broad absorption with λ(max) ~ 380-400 nm, τ = 2 μs). Calculations (B3LYP/6-31+G(d)) reveal that T(1K) of 1 is located 67 kcal/mol above its ground state (S(0)) and has a long C-N bond (1.58 ?), and the calculated transition state to form 3 is only 1 kcal/mol higher in energy than T(1K) of 1. The calculations show that 3 has significant 1,3-carbon iminyl biradical character, which explains why 3 reacts efficiently with oxygen and decays by intersystem crossing to the singlet surface. Photolysis of 1 in argon matrixes at 14 K produced ketene imine 7, which presumably is formed from 3 intersystem crossing to 7. In comparison, photolysis of 1 in methanol with a 266 nm laser produces mainly ylide 2 (λ(max) ~ 380 nm, τ ~ 6 μs, acetonitrile), which decays to form 2P. Ylide 2 is formed via singlet reactivity of 1, and calculations show that the first singlet excited state of the azirine chromophore (S(1A)) is located 113 kcal/mol above its S(0) and that the singlet excited state of the ketone (S(1K)) is 85 kcal/mol. Furthermore, the transition state for cleaving the C-C bond in 1 to form 2 is located 49 kcal/mol above the S(0) of 1. Thus, we theorize that internal conversion of S(1A) to a vibrationally hot S(0) of 1 forms 2, whereas intersystem crossing from S(1K) to T(1K) results in 3.     

Ab initio chemical kinetic study on Cl + ClO and related reverse processes

The reaction of ClO with Cl and its related reverse processes have been studied theoretically by ab initio quantum chemical and statistical mechanical calculations. The geometric parameters of the reactants, products, and transition states are optimized by both UMPW1PW91 and unrestricted coupled-cluster single and double excitation (UCCSD) methods with the 6-311+G(3df) basis set. The potential energy surface has been further refined (with triple excitations, T) at the UCCSD(T)/6-311+G(3df) level of theory. The results show that Cl(2) and O ((3)P) can be produced by chlorine atom abstraction via a tight transition state, while ClOCl ((1)A(1)) and ClClO ((1)A') can be formed by barrierless association processes with exothermicities of 31.8 and 16.0 kcal/mol, respectively. In principle the O ((1)D) atom can be generated with a large endothermicity of 56.9 kcal/mol; on the other hand, its barrierless reaction with Cl(2) can readily form ClClO ((1)A'), which fragments rapidly to give ClO + Cl. The rate constants of both forward and reverse processes have been predicted at 150-2000 K by the microcanonical variational transition state theory (VTST)/Rice-Ramsperger-Kassel-Marcus (RRKM) theory. The predicted rate constants are in good agreement with available experimental data within reported errors.     

New selective haloform-type reaction yielding 3-hydroxy-2,2-difluoroacids: theoretical study of the mechanism

Experimental results of an unprecedented haloform-type reaction in which 4-alkyl-4-hydroxy-3,3-difluoromethyl trifluoromethyl ketones undergo base-promoted selective cleavage of the CO-CF(3) bond, yielding 3-hydroxy-2,2-difluoroacids and fluoroform, are rationalized using DFT (B3LYP) calculations. The gas-phase addition of hydroxide ion to 1,1,1,3,3-pentafluoro-4-hydroxypentan-2-one (R) is found to be a barrierless process, yielding a tetrahedral intermediate (INT), involving a DeltaG(r)(298 K) of -61.4 kcal/mol. The CO-CF(3) bond cleavage in INT leads to a hydrogen-bonded [CH(3)CHOHCF(2)CO(2)H...CF(3)](-) complex by passage through a transition structure (TS1) with a DeltaG()(298 K) of 20.8 kcal/mol and a DeltaG(r)(298 K) of 9.8 kcal/mol. This complex undergoes a proton transfer between its components, yielding a hydrogen-bonded [CH(3)CHOHCF(2)CO(2)...CHF(3)](-) complex. This process has associated with it a DeltaG()(298 K) of only 3.1 kcal/mol and a DeltaG(r)(298 K) of -43.3 kcal/mol. The CO-CF(2) bond cleavage in INT leads to a hydrogen-bonded [CH(3)CHOHCF(2)...CF(3)CO(2)H](-) complex by passage through a transition structure (TS3) with a DeltaG()(298 K) of 29.2 kcal/mol and a DeltaG(r)(298 K) of 25.1 kcal/mol. The lower energy barrier found for CO-CF(3) bond cleavage in INT is ascribed to the larger number of fluorine atoms stabilizing the negative charge accumulated on the CF(3) moiety of TS1, as compared to the number of fluorine atoms stabilizing the negative charge on the CH(3)CHOHCF(2) moiety of TS3. The solvent-induced effects on the two pathways, introduced within the SCRF formalism through PCM calculations, do not reverse the predicted preference of the CO-CF(3) over the CO-CF(2) bond cleavage of R in the gas phase.     

The very low-pressure pyrolysis of phenyl ethyl ether,phenyl allyl ether,and benzyl methyl ether and the enthalpy of formation of the phenoxy radical

A value of the enthalpy of formation of the phenoxy radical in the gas phase, ΔH°_,298K (?O·, g) = 11.4 ± 2.0 kcal/mol, has been obtained from the kinetic study of the unimolecular decompositions of phenyl ethyl ether, phenyl allyl ether, and benzyl methyl ether

1 Trivial names for ethoxy benzene, 2-propenoxy (allyloxy) benzene, and α-methoxytoluene, respectively

at very low pressures. Bond fission, producing phenoxy or benzyl radicals, respectively, is the only mode of decomposition in each case. The present value leads to a bond dissociation energy BDE(?O—H) = 86.5 ± 2 kcal/mol,

2 1 kcal = 4.18674 kJ (absolute)

in good agreement with recent estimates made on the basis of competitive oxidation steps in the liquid phase. A comparison with bond dissociation energies of aliphatic alcohols, BDE(RO—H) = 104 kcal/mol, reveals that the stabilization energy of the phenoxy radical (17.5 kcal/mol) is considerably greater than the one observed for the isoelectronic benzyl radical (13.2 kcal/mol). Decomposition of phenoxy radicals into cyclopentadienyl radicals and CO has been observed at temperatures above 1000°K, and a mechanism for this reaction is proposed.     

The very low-pressure pyrolysis of phenyl methyl sulfide and benzyl methyl sulfide. The enthalpy of formation of the methylthio and phenylthio radicals

The kinetics and mechanisms of the unimolecular decompositions of phenyl methyl sulfide (PhSCH₃) and benzyl methyl sulfide (PhCH₂SCH₃) have been studied at very low pressures (VLPP). Both reactions essentially proceed by simple carbon-sulfur bond fission into the stabilized phenylthio (PhS·) and benzyl (PhCH₂·) radicals, respectively. The bond dissociation energies BDE(PhS-CH₃) = 67.5 ± 2.0 kcal/mol and BDE(PhCH₂-SCH₃) = 59.4 ± 2 kcal/mol, and the enthalpies of formation of the phenylthio and methylthio radicals ΔH° _,298K(PhS·, g) = 56.8 ± 2.0 kcal/mol and ΔH°_{f, 298K}(CH₃S·, g) = 34.2 ± 2.0 kcal/mol have been derived from the kinetic data, and the results are compared with earlier work on the same systems. The present values reveal that the stabilization energy of the phenylthio radical (9.6 kcal/mol) is considerably smaller than that observed for the related benzyl (13.2 kcal/mol) and phenoxy (17.5 kcal/mol) radicals.     

Bond dissociation energies in second-row compounds

Heats of formation at 0 and 298 K are predicted for PF3, PF5, PF3O, SF2, SF4, SF6, SF2O, SF2O2, and SF4O as well as a number of radicals derived from these stable compounds on the basis of coupled cluster theory [CCSD(T)] calculations extrapolated to the complete basis set limit. In order to achieve near chemical accuracy (+/-1 kcal/mol), additional corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, a correction for first-order atomic spin-orbit effects, and vibrational zero-point energies. The calculated values substantially reduce the error limits for these species. A detailed comparison of adiabatic and diabatic bond dissociation energies (BDEs) is made and used to explain trends in the BDEs. Because the adiabatic BDEs of polyatomic molecules represent not only the energy required for breaking a specific bond but also contain any reorganization energies of the bonds in the resulting products, these BDEs can be quite different for each step in the stepwise loss of ligands in binary compounds. For example, the adiabatic BDE for the removal of one fluorine ligand from the very stable closed-shell SF6 molecule to give the unstable SF5 radical is 2.8 times the BDE needed for the removal of one fluorine ligand from the unstable SF5 radical to give the stable closed-shell SF4 molecule. Similarly, the BDE for the removal of one fluorine ligand from the stable closed-shell PF3O molecule to give the unstable PF2O radical is higher than the BDE needed to remove the oxygen atom to give the stable closed-shell PF3 molecule. The same principles govern the BDEs of the phosphorus fluorides and the sulfur oxofluorides. In polyatomic molecules, care must be exercised not to equate BDEs with the bond strengths of given bonds. The measurement of the bond strength or stiffness of a given bond represented by its force constant involves only a small displacement of the atoms near equilibrium and, therefore, does not involve any reorganization energies, i.e., it may be more appropriate to correlate with the diabatic product states.     

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.     

A semiempirical study of the optimized ground and excited state potential energy surfaces of retinal and its protonated Schiff base.

The electronic ground and first excited states of retinal and its Schiff base are optimized for the first time using the semiempirical AM1 Hamiltonian. The barrier for rotation about the C(11)-C(12) double bond is characterized by variation of both the twist angle delta(C(10)-C(11)-C(12)-C(13)) and the bond length d(C(11)-C(12)). The potential energy surface is obtained by varying these two parameters. The calculated ground state rotational barrier is equal to 15.6 kcal/mol for retinal and 20.5 kcal/mol for its Schiff base. The all-trans conformation is more stable by 3.7 kcal/mol than the 11-cis geometry. For the first excited state, S(1,) the 90 degrees twisted geometry represents a saddle point for retinal with the rotational barrier of 14.6 kcal/mol. In contrast, this conformation is an energy minimum for the Schiff base. It can be easily reached at room temperature from the planar minima since it is separated from them by a barrier of only 0.6 kcal/mol. The 90 degrees minimum conformation is more stable than the all-trans by 8.6 kcal/mol. We are thus able to present a reaction path on the S(1) surface of the retinal Schiff base with an almost barrier-less geometrical relaxation into a twisted minimum geometry, as observed experimentally. The character of the ground and first excited singlet states underscores the need for the inclusion of double excitations in the calculations.     

Reaction of 4-oxo-2,3,5,6-tetrafluorocyclohexa-2,5-dienylidene with acetylene: a carbene to carbene reaction.

The EPR spectrum of triplet 4-oxo-2,3,5,6-tetrafluorocyclohexa-2,5-dienylidene 1 was recorded in solid argon at 15 K. Carbene 1 reacts with acetylene under the conditions of matrix isolation yielding triplet vinylmethylene 4, which was characterized by its IR, UV-vis, and EPR spectrum. Carbene 4 is photolabile and is converted to spiro compound 5 on irradiation with lambda > 515 nm. The reaction of triplet carbene 1 with acetylene to produce triplet carbene 4 is predicted to be exothermic by 55 kcal mol(-1) at the B3LYP/6-31G(d,p) level of theory. The cis isomer is calculated to be only 0.4 kcal mol(-1) less stable than trans-4 at this level of theory. According to our calculations, singlet carbene S-4 is not a minimum on the C(8)F(4)H(2)O potential energy surface; however, at the T-4 geometry, the lowest lying singlet state is predicted to be 20.7 kcal mol(-1) higher in energy. The subsequent photochemical cyclization of T-4 yielding spiro compound 5 is exothermic by 10.3 kcal mol(-1) relative to T-4 and by 31.1 kcal mol(-1) relative to S-4. 4-Ethinyl-2,3,5,6-tetrafluorocyclohexa-2,5-dienone 9, the C-H insertion product of 1 and acetylene, was not observed experimentally, although it is favored energetically by 4.3 kcal mol(-1) over 5.     

Computational study of bond dissociation enthalpies for substituted β-O-4 lignin model compounds

The biopolymer lignin is a potential source of valuable chemicals. Phenethyl phenyl ether (PPE) is representative of the dominant β-O-4 ether linkage. DFT is used to calculate the Boltzmann-weighted carbon-oxygen and carbon-carbon bond dissociation enthalpies (BDEs) of substituted PPE. These values are important for understanding lignin decomposition. Exclusion of all conformers that have distributions of less than 5% at 298 K impacts the BDE by less than 1 kcal mol(-1). We find that aliphatic hydroxyl/methylhydroxyl substituents introduce only small changes to the BDEs (0-3 kcal mol(-1)). Substitution on the phenyl ring at the ortho position substantially lowers the C-O BDE, except in combination with the hydroxyl/methylhydroxyl substituents, for which the effect of methoxy substitution is reduced by hydrogen bonding. Hydrogen bonding between the aliphatic substituents and the ether oxygen in the PPE derivatives has a significant influence on the BDE. CCSD(T)-calculated BDEs and hydrogen-bond strengths of ortho-substituted anisoles, when compared with M06-2X values, confirm that the latter method is sufficient to describe the molecules studied and provide an important benchmark for lignin model compounds.     

Reactions of 1,3-cyclohexadiene with singlet oxygen. A theoretical study.

A thorough study of the reaction of singlet oxygen with 1,3-cyclohexadiene has been made at the B3LYP/6-31G(d) and CASPT2(12e,10o) levels. The initial addition reaction follows a stepwise diradical pathway to form cyclohexadiene endoperoxide with an activation barrier of 6.5 kcal/mol (standard level = CASPT2(12e,10o)/6-31G(d); geometries and zero-point corrections at B3LYP/6-31G(d)), which is consistent with an experimental value of 5.5 kcal/mol. However, as the enthalpy of the transition structure for the second step is lower than the diradical intermediate, the reaction might also be viewed as a nonsynchronous concerted reaction. In fact, the concertedness of the reaction is temperature dependent since entropy differences create a free energy barrier for the second step of 1.8 kcal/mol at 298 K. There are two ene reactions; one is a concerted mechanism (DeltaH(double dagger) = 8.8 kcal/mol) to 1-hydroperoxy-2,5-cyclohexadiene (5), while the other, which forms 1-hydroperoxy-2,4-cyclohexadiene (18), passes through the same diradical intermediate (9) as found on the pathway to endoperoxide. The major pathway from the endoperoxide is O-O bond cleavage (22.0 kcal/mol barrier) to form a 1,4-diradical (25), which is 13.9 kcal/mol less stable than the endoperoxide. From the diradical, two low-energy pathways exist, one to epoxyketone (29) and the other to the diepoxide (27), where both products are known to be formed experimentally with a product ratio sensitive to the nature of substitutents. A significantly higher activation barrier leads to C-C bond cleavage and direct formation of maleic aldehyde plus ethylene.