MPW1K Performs Much Better than B3LYP in DFT Calculations on Reactions that Proceed by Proton-Coupled Electron Transfer (PCET)

DFT calculations have been performed with the B3LYP and MPW1K functional on the hydrogen atom abstraction reactions of ethenoxyl with ethenol and of phenoxyl with both phenol and alpha-naphthol. Comparison with the results of G3 calculations shows that B3LYP seriously underestimates the barrier heights for the reaction of ethenoxyl with ethenol by both proton-coupled electron transfer (PCET) and hydrogen atom transfer (HAT) mechanisms. The MPW1K functional also underestimates the barrier heights, but by much less than B3LYP. Similarly, comparison with the results of experiments on the reaction of phenoxyl radical with alpha-naphthol indicates that the barrier height for the preferred PCET mechanism is calculated more accurately by MPW1K than by B3LYP. These findings indicate that the MPW1K functional is much better suited than B3LYP for calculations on hydrogen abstraction reactions by both HAT and PCET mechanisms.     

IR laser-induced dehydrohalogenation of halogenated ethanes

Pulsed IR radiation from a CO₂ TEA laser has been shown to cause dehydrohalogenation of CF₃CH₃, CF₂ClCH₃ and C₂H₅Cl under collisionless conditions. IR emission has been observed and is attributed to HF² and CF₂CH₂² from CF₃CH₃, and to HF² and HCl² from CF₂ClCH₃.     

Are 1,5-disubstituted semibullvalenes that have C2v equilibrium geometries necessarily bishomoaromatic?

Unrestricted density functional theory (UB3LYP), CASSCF, and CASPT2 calculations have been employed to compute the relative energies of the C(s) and C(2v) geometries of several 1,5-disubstituted semibullvalenes. Substitution at these positions with R = F, -CH(2)-, or -O- affords semibullvalenes that are predicted to have C(2v) equilibrium geometries. Calculated singlet-triplet energy splittings and the energies of isodesmic reactions are used to assess the amount of bishomoaromatic character at these geometries. The results of these calculations show that employing strain to destabilize the C(s) geometries of semibullvalenes can lead to a significant decrease in the amount of bishomoaromatic stabilization of the C(2v) geometries, due to reduced through-space interaction between the two allyl groups. However, the C(2v) equilibrium geometries of the 1,5-disubstituted semibullvalenes with R = F and -RR- = -O- do benefit from stabilizing through-bond interactions between the two allyl groups. These interactions involve mixing of the bisallyl HOMO with the low-lying C-F or C-O sigma orbital combinations of the same symmetry. In contrast, for -RR- = -CH(2)-, through-bond interactions destabilize the bisallyl HOMO and are predicted to make the ground state of this semibullvalene a triplet.     

Ligand-assisted reduction of osmium tetroxide with molecular hydrogen via a [3+2] mechanism

Osmium tetroxide is reduced by molecular hydrogen in the presence of ligands in both polar and nonpolar solvents. In CHCl3 containing pyridine (py) or 1,10-phenanthroline (phen), OsO4 is reduced by H2 to the known Os(VI) dimers L2Os(O)2(mu-O)2Os(O)2L2 (L2 = py2, phen). However, in the absence of ligands in CHCl3 and other nonpolar solvents, OsO4 is unreactive toward H2 over a week at ambient temperatures. In basic aqueous media, H2 reduces OsO4(OH)n(n-) (n = 0, 1, 2) to the isolable Os(VI) complex, OsO2(OH)4(2-), at rates close to that found in py/CHCl3. Depending on the pH, the aqueous reactions are exergonic by deltaG = -20 to -27 kcal mol(-1), based on electrochemical data. The second-order rate constants for the aqueous reactions are larger as the number of coordinated hydroxide ligands increases, k(OsO4) = 1.6(2) x 10(-2) M(-1) s(-1) < k(OsO4(OH)-) = 3.8(4) x 10(-2) M(-1) s(-1) < k(OsO4(OH)2(2-)) = 3.8(4) x 10(-1) M(-1) s(-1). The observation of primary deuterium kinetic isotope effects, k(H2)/k(D2) = 3.1(3) for OsO4 and 3.6(4) for OsO4(OH)-, indicates that the rate-determining step in each case involves H-H bond cleavage. Density functional calculations and thermochemical arguments favor a concerted [3+2] addition of H2 across two oxo groups of OsO4(L)n and argue against H* or H- abstraction from H2 or [2+2] addition of H2 across one Os=O bond. The [3+2] mechanism is analogous to that of alkene addition to OsO4(L)n to form diolates, for which acceleration by added ligands has been extensively documented. The observation that ligands also accelerate H2 addition to OsO4(L)n highlights the analogy between these two reactions.     

Controlling the DNA binding specificity of bHLH proteins through intramolecular interactions

Reversible control of the conformation of proteins was employed to probe the relationship between flexibility and specificity of the basic helix-loop-helix protein MyoD. A fusion protein (apaMyoD) was designed where the basic DNA binding helix of MyoD was stablized by an amino-terminal extension with a sequence derived from the bee venom peptide apamin. The disulfide-stabilized helix from apamin served as a nucleus for a helix that extended for a further ten residues, thereby holding apaMyoD's DNA recognition helix in a predominantly alpha-helical conformation. The thermal stability of the DNA complexes of apaMyoD was increased by 13 degrees C relative to MyoD-bHLH. Measurements of the fluorescence anisotropy change on DNA binding indicated that apaMyoD bound to E-box-containing DNA sequences with enhanced affinity relative to MyoD-bHLH. Consequently, the DNA binding specificity of apaMyoD was increased 10-fold.     

Do deviations from bond enthalpy additivity define the thermodynamic stabilities of diradicals?

Deviations from bond enthalpy additivity (DeltaBEA) are frequently used to assess the thermodyamic stabilities of diradicals. (U)B3LYP/6-31G calculations have been performed in order to determine how well DeltaBEA values actually do reflect the thermodynamic stabilities of the triplet states of diradicals in which one or both nonbonding electrons occupy a delocalized pi orbital. The calculations find that different pathways for forming sigma,pi-diradicals, such as alpha,2- and alpha,4-dehydrotoluene (4 and 6), give DeltaBEA values that differ by ca. 1 kcal/mol. The path dependency of the DeltaBEA values is computed to be one order of magnitude larger for non-Kekulé hydrocarbon diradicals, such as m-benzoquinodimethane (12) and 1,3-dimethylenecyclobutane-2,4-diyl (15), than for sigma,pi-diradicals. Since the DeltaBEA values for forming 4, 6, 12, and 15 are all path dependent, we conclude that DeltaBEA values for diradicals with one or two delocalized, nonbonding pi electrons do not, in general, uniquely define the thermodynamic stabilities of the diradicals. Hence, DeltaBEA values should not be used for this purpose, especially for non-Kekulé hydrocarbon diradicals.     

Deactivation of vibrationally excited ozone by O(3P) atoms

Using laser-induced fluorescence of ozone (to measure the rate of disappearance of O₃²) and NO₂ titration (to determine O atom concentrations), we have determined bimolecular rate constants for the deactivation by O(³P atoms) of ozone in excited stretching and bending modes. These experiments do not distinguish between deactivation by (a) the exchange of vibrational and translational energy or (b) the chemical reaction O₃ + O → 2O₂. If the non-reactive pathway (a) is assumed to dominate, then O(³P) is 150 times more effective than O₂ in deactivating O²₃. If chemical reaction (b) is dominant, the bimolecular rate constant for O²₃ + O(³P) is larger by a factor of 150–1500 than that for ground-state ozone.     

The search for bishomoaromatic semibullvalenes and barbaralanes: computational evidence of their identification by UV/Vis and IR spectroscopy and prediction of the existence of a blue bishomoaromatic semibullvalene

Time-dependent B3LYP/6-31G calculations have been performed at the optimized C(2) or C(2v) geometries of several substituted semibullvalenes (1(deloc)) and barbaralanes (2(deloc)), to compare the computed vertical electronic excitation energies with the temperature-dependent, long-wavelength absorptions that have been observed in the UV/vis spectra of some of these compounds by Quast and co-workers. The excellent agreement between the calculated vertical excitation energies and the observed values of lambda(max) provides strong support for the identification of the bishomoaromatic species 1(deloc) and 2(deloc) as the source of these absorptions. Furthermore, the CN stretching frequencies, computed for the C(2) geometry of 1,5-dimethyl-2,6-dicyano-4,8-diphenylsemibullvalene (1f(deloc)), fit the low-frequency absorptions seen in the IR spectrum of 1f, thus furnishing independent evidence that bishomoaromatic geometries of semibullvalenes have, in fact, been observed spectroscopically. B3LYP/6-31G calculations predict that 2,6-dicyano-4,8-diphenylsemibullvalene 1c has a C(2) equilibrium geometry (1c(deloc)) and that the long-wavelength UV/vis absorption (lambda(max) = 585 nm) and CN stretching frequencies (2192 and 2194 cm(-1)) computed for 1c(deloc) should serve to identify this bishomoaromatic semibullvalene when it is synthesized.     

Computational studies of the thermal fragmentation of P-arylphosphiranes: have arylphosphinidenes been generated by this method?

CASSCF, CASPT2, CCSD(T), and (U)B3LYP electronic structure calculations have been performed in order to investigate the thermal fragmentation of P-phenylphosphirane (1) to phenylphosphinidene (PhP) and ethylene. The calculations show that generation of triplet PhP via a stepwise pathway is 21 kcal mol(-1) less endothermic and has a 12 kcal mol(-1) lower barrier height than concerted fragmentation of 1 to give singlet PhP. The formation of singlet PhP via a concerted pathway is predicted to be stereospecific, whereas formation of triplet PhP is predicted to occur with complete loss of stereochemistry. However, calculations on fragmentation of anti-cis-2,3-dimethyl-P-mesitylphosphirane (cis-1Me) to triplet mesitylphosphinidene (MesP) indicate that this reaction should be more stereospecific, in agreement with the experimental results of Li and Gaspar. Nevertheless, with a predicted free energy of activation of 42 kcal mol(-1), the formation of MesP from cis-1Me is not likely to have occurred in an uncatalyzed reaction at the temperatures at which this phosphirane has been pyrolyzed.     

Like (CO)4, Do (CS)4 and (CSe)4 Have a Triplet Ground State?

Cyclobutane‐1,2,3,4‐tetraone, (CO)₄, was computationally predicted and, subsequently, experimentally confirmed to have a triplet ground state, in which a b_2g σ MO and an a_2u π MO were each singly occupied. In contrast, the (U)CCSD(T) calculations reported herein found that cyclobutane‐1,2,3,4‐tetrathione, (CS)₄, and cyclobutane‐1,2,3,4‐tetraselenone, (CSe)₄, both had singlet ground states, in which the b_2g σ MO was doubly occupied and the a_2u π MO was empty. Our calculations showed that both the longer C?X distances and smaller coefficients on the carbon atoms in the b_2g and a_2u MOs of (CS)₄ and (CSe)₄ contributed to the difference between the ground states of these two molecules and the ground state of (CO)₄. An experimental test of the prediction of a singlet ground state for (CS)₄ is proposed.