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.     

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.     

Alkali-metal ion catalysis and inhibition in nucleophilic displacement reactions at carbon,phosphorus and sulfur centres. IX. p-Nitrophenyl diphenyl phosphate

We report on the catalytic effects by alkali-metal ions in the ethanolysis of p-nitrophenyl diphenyl phosphate, in continuation of our studies on alkali-metal ion catalysis and inhibition in nucleophilic displacement reactions at carbon, phosphorus and sulfur centres. The following selectivity order of catalytic reactivity was observed for nucleophilic displacement at the phosphorus center with p-nitrophenoxide as leaving group: Li+ > Na+ > K+ > Cs+. A minor reaction pathway with phenoxide leaving was also found. The kobs data have been dissected into reaction pathways by free ions (kEtO) and by ion pairs (kMOEt), with the latter being dominant, in a 4-membered transition-state. Further analysis is given in terms of initial-state and transition-state stabilization by the alkali-metal ions in terms of the Eisenman model (electrostatic interaction vs. desolvation). Results of ab-initio MO calculations are presented based on interaction between M+ and a model bipyramidal phosphorane intermediate and compared with the sulfurane analogue.     

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.     

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.     

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.     

Structure of Amido Pyridinium Betaines: Persistent Intermolecular C−H⋅⋅⋅N Hydrogen Bonding in Solution

A hydrogen bond of the type C?H???X (X=O or N) is known to influence the structure and function of chemical and biological systems in solution. C?H???O hydrogen bonding in solution has been extensively studied, both experimentally and computationally, whereas the equivalent thermodynamic parameters have not been enumerated experimentally for C?H???N hydrogen bonds. This is, in part, due to the lack of systems that exhibit persistent C?H???N hydrogen bonds in solution. Herein, a class of molecule based on a biologically active norharman motif that exhibits unsupported intermolecular C?H???N hydrogen bonds in solution has been described. A pairwise interaction leads to dimerisation to give bond strengths of about 7 kJ mol^?1 per hydrogen bond, which is similar to chemically and biologically relevant C?H???O hydrogen bonding. The experimental data is supported by computational work, which provides additional insight into the hydrogen bonding by consideration of electrostatic and orbital interactions and allowed a comparison between calculated and extrapolated NMR chemical shifts.     

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.     

Heavy‐Atom Tunneling in the Ring Opening of a Strained Cyclopropene at Very Low Temperatures

The highly strained 1H‐bicyclo[3.1.0]‐hexa‐3,5‐dien‐2‐one 1 is metastable, and rearranges to 4‐oxacyclohexa‐2,5‐dienylidene 2 in inert gas matrices (neon, argon, krypton, xenon, and nitrogen) at temperatures as low as 3 K. The kinetics for this rearrangement show pronounced matrix effects, but in a given matrix, the reaction rate is independent of temperature between 3 and 20 K. This temperature independence means that the activation energy is zero in this temperature range, indicating that the reaction proceeds through quantum mechanical tunneling from the lowest vibrational level of the reactant. At temperatures above 20 K, the rate increases, resulting in curved Arrhenius plots that are also indicative of thermally activated tunneling. These experimental findings are supported by calculations performed at the CASSCF and CASPT2 levels by using the small‐curvature tunneling (SCT) approximation.