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.     

2-Pyridinethiol/2-pyridinethione tautomeric equilibrium. A comparative experimental and computational study

The gas phase and solvent dependent preference of the tautomerization between 2-pyridinethiol (2SH) and 2-pyridinethione (2S) has been assessed using variable temperature Fourier transform infrared (FTIR) experiments, as well as ab initio and density functional theory computations. No spectroscopic evidence (nu(S)(-)(H) stretch) for 2SH was observed in toluene, C(6)D(6), heptane, or methylene chloride solutions. Although, C(s)() 2SH is 2.61 kcal/mol more stable than C(s)() 2S (CCSD(T)/cc-pVTZ//B3LYP/6-311+G(3df,2p)+ZPE), cyclohexane solvent-field relative energies (IPCM-MP2/6-311+G(3df,2p)) favor 2S by 1.96 kcal/mol. This is in accord with the FTIR observations and in quantitative agreement with the -2.6 kcal/mol solution (toluene or C(6)D(6)) calorimetric enthalpy for the 2S/2SH tautomerization favoring the thione. As the intramolecular transition state for the 2S, 2SH tautomerization (2TS) lies 25 (CBS-Q) to 30 kcal/mol (CCSD/cc-pVTZ) higher in energy than either tautomer, tautomerization probably occurs in the hydrogen bonded dimer. The B3LYP/6-311+G(3df,2p) optimized C(2) 2SH dimer is 10.23 kcal/mol + ZPE higher in energy than the C(2)(h)() 2S dimer and is only 2.95 kcal/mol + ZPE lower in energy than the C(2) 2TS dimer transition state. Dimerization equilibrium measurements (FTIR, C(6)D(6)) over the temperature range 22-63 degrees C agree: K(eq)(298) = 165 +/- 40 M(-)(1), DeltaH = -7.0 +/- 0.7 kcal/mol, and DeltaS = -13.4 +/- 3.0 cal/(mol deg). The difference between experimental and B3LYP/6-311+G(3df,2p) [-34.62 cal/(mol deg)] entropy changes is due to solvent effects. The B3LYP/6-311+G(3df,2p) nucleus independent chemical shifts (NICS) are -8.8 and -3.5 ppm 1 A above the 2SH and 2S ring centers, respectively, and the thiol is aromatic. Although the thione is not aromatic, it is stabilized by the thioamide resonance. In solvent, the large 2S dipole, 2-3 times greater than 2SH, favors the thione tautomer and, in conclusion, 2S is thermodynamically more stable than 2SH in solution.     

Unimolecular reactions of chemically activated CF2BrCF2CH3 and CF2BrCF2CD3: evidence for 1,2-FBr interchange

Vibrationally excited CF2BrCF2CH3 and CF2BrCF2CD3 molecules were prepared with 96 kcal mol-1 energy at room temperature by the recombination of CF2BrCF2 and CH3 (CD3) radicals. The observed unimolecular reactions are 1,2-BrF interchange to give CF3CFBrCH3 (CD3) molecules and 2,3-FH (FD) elimination; the rate constants are 2.2 x 10(5) (1.5 x 10(5)) s(-1) and 2.0 x 105 (0.75 x 10(5)) s(-1), respectively. The CF3CFBrCH3 (CD3) molecules rapidly, relative to the reverse reaction, eliminate HBr or DBr to give the observed product CF3CF=CH2 (CD2). Density functional theory at the B3PW91/6-311+G(2d,p) level was used to obtain vibrational frequencies and moments of inertia of the molecule and transition states for subsequent calculations of statistical rate constants for CF2BrCF2CH3 and CF2BrCF2CD3. Matching experimental and calculated rate constants gave threshold energies of 62 and 66 kcal mol-1 for 1,2-BrF interchange and 2,3-FH elimination, respectively. The BrF interchange reaction is compared to ClF interchange from CF2ClCF2CH3 and CF2ClCHFCH3.     

Effective monkey saddle points and berry and lever mechanisms in the topomerization of SF(4) and related tetracoordinated AX(4) species

The topomerization mechanisms of the SF(4) and SCl(2)F(2) sulfuranes, as well as their higher (SeF(4), TeF(4)) and isoelectronic analogues PF(4)(-), AsF(4)(-), SbF(4)(-), SbCl(4)(-), ClF(4)(+), BrF(4)(+), BrCl(2)F(2)(+), and IF(4)(+)), have been computed at B3LYP/6-31+G and at B3LYP/6-311+G. All species have trigonal bipyramidal (TBP) C(2)(v)() ground states. In such four-coordinated molecules, Berry rotation exchanges both axial with two equatorial ligands simultaneously while the alternative "lever" mechanism exchanges only one axial ligand with one equatorial ligand. While the barrier for the lever exchange in SF(4) (18.8 kcal mol(-1)) is much higher than that for the Berry process (8.1 kcal mol(-1)), both mechanisms are needed for complete ligand exchange. The F(ax)F(ax) and F(eq)F(eq) isomers of SF(2)Cl(2) have nearly the same energy and readily interconvert by BPR with a barrier of 7.6 kcal mol(-1). The enantiomerization of the F(ax)F(eq) chiral isomer can occur by either the Berry process (transition state barrier 8.3 kcal mol(-1)) or the "lever" mechanism via either of two C(s)() transition states, based on the TBP geometry: Cl(ax) <--> Cl(eq) or F(ax) <--> F(eq) exchanges with barriers of 6.3 and 15.7 kcal mol(-1), respectively. Full scrambling of all ligand sites is possible only by inclusion of the lever mechanism. Planar, "tetrahedral", and triplet forms are much higher in energy. The TBP C(3)(v) structures of AX(4) either have two imaginary frequencies (NIMAG = 2) for the X = F, Cl species or are minima (NIMAG = 0) for the X = Br, I compounds. These "effective monkey saddle points" have degenerate modes with two small frequencies, imaginary or real. Although a strictly defined "monkey saddle" (with degenerate frequencies exactly zero) is not allowed, the flat C(3)(v) symmetry region serves as a "transition state" for trifurcation of the pathways. The BPR mechanism also is preferred over the alternative lever process in the topomerization of the selenurane SeF(4) (barriers 5.9 vs. 12.1 kcal mol(-1)), the tellurane TeF(4) (2.1 vs. 6.4), and the interhalogen cations ClF(4)(+) (2.5 vs 14.8), BrF(4)(+) (4.7 vs. 11.3), BrF(2)Cl(2)(+) (14.6 vs. 17.4), and IF(4)(+) (1.4 vs. 6.0), as well as for the series PF(4)(-) (7.0 vs. 9.0), AsF(4)(-) (9.3 vs. 17.2), and SbF(4)(-) (3.8 vs. 5.3 kcal mol(-1)), all computed at B3LYP/6-311+G with the inclusion of quasirelativistic pseudopotentials for Te, I, and Sb. The heavier halogens increasingly favor the lever process, where the barrier (2.6 kcal mol(-1)) pertaining to the effective monkey saddle point (C(3)(v) minimum for SbCl(4)(-)) is less than that for the Berry process (8.2 kcal mol(-1)).     

Halogenation of cubane under phase-transfer conditions: single and double C-H-bond substitution with conservation of the cage structure.

The first highly selective C-H chlorination, bromination, and iodination of cubane (1) utilizing polyhalomethanes as halogen sources under phase-transfer (PT) conditions is described. Isomeric dihalocubanes with all possible combinations of chlorine, bromine, and iodine in ortho, meta, and para positions were also prepared by this method; m-dihalo products form preferentially. Ab initio and density functional theory (DFT) computations were used to rationalize the pronounced differences in the reactions of 1 with halogen (Hal(*)) vs carbon-centered trihalomethyl (Hal(3)C(*)) radicals (Hal = Cl, Br). For Hal(3)C radicals the C-H abstraction pathway is less unfavorable (DeltaG(double dagger)(298) = 21.6 kcal/mol for Cl(3)C(*) and 19.4 kcal/mol for Br(3)C(*) at B3LYP/6-311+G//B3LYP/6-31G) than the fragmentation of the cubane skeleton via S(H)2-attack on one of the carbon atoms of 1 (DeltaG(double dagger)(298) = 33.8 and 35.1 kcal/mol, respectively). In stark contrast, the reaction of 1 with halogen atoms preferentially follows the fragmentation pathway (DeltaG(double dagger)(298) = 2.1 and 7.5 kcal/mol) and C-H abstraction is more unfavorable (DeltaG(double dagger)(298) = 4.6 and 12.0 kcal/mol). Our computational results nicely agree with the behavior of 1 under PT halogenation conditions (where Hal(3)C(*) is involved in the activation step) and under free-radical photohalogenation with Hal(2) (Della, E. W., et al. J. Am. Chem. Soc. 1992, 114, 10730). The incorporation of a second halogen atom preferentially in the meta position of halocubanes demonstrates the control of the regioselectivity by molecular orbital symmetry.     

Cheletropic decomposition of cyclic nitrosoamines revisited: the nature of the transition states and a critical role of the ring strain

The cheletropic decompositions of 1-nitrosoaziridine (1), 1-nitroso-Delta(3)-pyrroline (2), 7-nitroso-7-azabicyclo[2.2. 1]hepta-2,5-diene (3), and 6-nitroso-6-azabicyclo[2.1.1]hexa-4-ene (4) have been studied theoretically using high level ab initio computations. Activation parameters of the decomposition of nitrosoaziridine 1 were obtained experimentally in heptane (DeltaH()(298) = 18.6 kcal mol(-)(1), DeltaS()(298) = -7.6 cal mol(-)(1) K(-)(1)) and methanol (20.3 kcal mol(-)(1), 0.3 cal mol(-)(1) K(-)(1)). Among employed theoretical methods (B3LYP, MP2, CCD, CCSD(T)//CCD), the B3LYP method in conjunction with 6-31+G, 6-311+G, and 6-311++G(3df,2pd) basis sets gives the best agreement with experimental data. It was found that typical N-nitrosoheterocycles 2-4 which have high N-N bond rotation barriers (>16 kcal mol(-)(1)) extrude nitrous oxide via a highly asynchronous transition state with a planar ring nitrogen atom. Nitrosoaziridine 1, with a low rotation barrier (<9 kcal mol(-)(1)) represents a special case. This compound can eliminate N(2)O via a low energy linear synperiplanar transition state (DeltaH()(298) = 20.6 kcal mol(-)(1), DeltaS()(298) = 2.5 cal mol(-)(1) K(-)(1)). Two higher energy transition states are also available. The B3LYP activation barriers of the cheletropic fragmentation of nitrosoheterocycles 2-4 decrease in the series: 2 (58 kcal mol(-)(1)) > 3 (18 kcal mol(-)(1)) > 4 (12) kcal mol(-)(1). The relative strain energies increase in the same order: 2 (0 kcal mol(-)(1)) < 3 (39 kcal mol(-)(1)) < 4 (52 kcal mol(-)(1)). Comparison of the relative energies of 2-4 and their transition states on a common scale where the energy of nitrosopyrroline 2 is assumed as reference indicates that the thermal stability of the cyclic nitrosoamines toward cheletropic decomposition is almost entirely determined by the ring strain.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

o-Quinone methide as alkylating agent of nitrogen, oxygen, and sulfur nucleophiles. The role of H-bonding and solvent effects on the reactivity through a DFT computational study

The reactivity of the alkylating agent o-quinone methide (o-QM) toward NH(3), H(2)O, and H(2)S, prototypes of nitrogen-, oxygen-, and sulfur-centered nucleophiles, has been studied by quantum chemical methods in the frame of DF theory (B3LYP) in reactions modeling its reactivity in water with biological nucleophiles. The computational analysis explores the reaction of NH(3), H(2)O, and H(2)S with o-QM, both free and H-bonded to a discrete water molecule, with the aim to rationalize the specific and general effect of the solvent on o-QM reactivity. Optimizations of stationary points were done at the B3LYP level using several basis sets [6-31G(d), 6-311+G(d,p), adding d and f functions to the S atom, 6-311+G(d,p),S(2df), and AUG-cc-pVTZ]. The activation energies calculated for the addition reactions were found to be reduced by the assistance of a water molecule, which makes easier the proton-transfer process in these alkylation reactions by at least 12.9, 10.5, and 6.0 kcal mol(-1) [at the B3LYP/AUG-cc-pVTZ//B3LYP/6-311+G(d,p) level], for ammonia, water, and hydrogen sulfide, respectively. A proper comparison of an uncatalyzed with a water-catalyzed reaction mechanism has been made on the basis of activation Gibbs free energies. In gas-phase alkylation of ammonia and water by o-QM, reactions assisted by an additional water molecule H-bonded to o-QM (water-catalyzed mechanism) are favored over their uncatalyzed counterparts by 5.6 and 4.0 kcal mol(-1) [at the B3LYP/6-311+G(d,p) level], respectively. In contrast, the hydrogen sulfide alkylation reaction in the gas phase shows a slight preference for a direct alkylation without water assistance, even though the free energy difference (DeltaDeltaG(#)) between the two reaction mechanisms is very small (by 1.0 kcal mol(-1) at the B3LYP/6-311+G(d,p),S(2df) level of theory). The bulk solvent effect, evaluated by the C-PCM model, significantly modifies the relative importance of the uncatalyzed and water-assisted alkylation mechanism by o-QM in comparison to the case in the gas phase. Unexpectedly, the uncatalyzed mechanism becomes highly favored over the catalyzed one in the alkylation reaction of ammonia (by 7.0 kcal mol(-1)) and hydrogen sulfide (by 4.0 kcal mol(-1)). In contrast, activation induced by water complexation still plays an important role in the o-QM hydration reaction in water as solvent.     

Absolute rate constants for reactions of tributylstannyl radicals with bromoalkanes, episulfides, and alpha-halomethyl-episulfides, -cyclopropanes, and -oxiranes: new rate expressions for sulfur and bromine atom abstraction

Arrhenius rate expressions were determined for the abstraction of bromine atom from 2-phenethyl bromide by tri-n-butylstannyl radical (Bu(3)Sn(*)) in benzene using transient absorption spectroscopy, (log(k(abs,Br)/M(-1) s(-1)) = (9.21 +/- 0.20) - (2.23 +/- 0.28)/theta, theta = 2.3RT kcal/mol, errors are 2sigma) and for the abstraction of sulfur atom from propylene sulfide to form propylene, (log(k(s)/M(-1) s(-1)) = (8.75 +/- 0.91) - (2.35 +/-1.33)/theta). Rate constants for reactions of organic bromides, RBr, with Bu(3)Sn(*) were found to vary as R = benzyl (15.6) > thiiranylmethyl (6.2) > oxiranylmethyl (3.1) > cyclopropylmethyl (1.3) > 2-phenethyl (1.0), with k(abs,Br) = 6.8 x 10(7) M(-1) s(-1) at 353 K for 2-phenethyl bromide. Bromine abstraction from alpha-bromomethylthiirane is about 7-fold faster than sulfur atom abstraction and is comparable to the reactivity of a secondary alkyl bromide. The potential surface for the vinylthiomethyl --> allylthiyl radical rearrangement at UB3LYP/6-31G(d) and UB3LYP/6-311+G(2d,2p) levels of theory suggests that the thiiranylmethyl radical is produced about 9 kcal/mol above the allylthiyl radical on the rearrangement surface, consistent with the observed enhancement of the Br atom abstraction from the thiirane and with synchronous C-S bond scission of the thiirane ring. The selectivities reported in this work for S vs Cl and Br abstraction provide applications for radical-based synthesis and new competition basis rate expressions for trialkylstannyl radicals.     

Conformational energies for 2-substituted butanes

The conformational free energies for some 2-substituted butanes where X = F, Cl, CN, and CCH were calculated using G3-B3, CBS-QB3, and CCSD(T)/6-311++G(2d,p) as well as other theoretical levels. The above methods gave consistent results with free energies relative to the trans conformers as follows: X = CCH, g+ = 0.77 +/- 0.05 kcal/mol. g- = 0.88 +/- 0.05 kcal/mol; X = CN, g+ = 0.85 +/- 0.05 kcal/mol, g- = 0.75 +/- 0.05 kcal/mol; X = Cl, g+ = 0.70 +/- 0.05 kcal/ml, g- = 0.80 +/- 0.05 kcal/mol; and X = F, g+ = 0.53 +/- 0.05 kcal/mol, g- = 0.83 +/- 0.05 kcal/mol. The conformational free energies also were estimated using the observed liquid phase IR spectra and intensities calculated using B3LYP/6-311++G** and MP2/6-311++G**. The rotational free energy profiles for all of the compounds were estimated at the G3-B3 level.     

Can the ebselen derivatives catalyze the isomerization of peroxynitrite to nitrate?

The reaction of ebselen and its derivatives (1-7) with peroxynitrite anion (ONOO(-); PN) has been studied in gas phase and in aqueous, dichloromethane, benzene, and cyclohexane solutions using B3LYP/6-311+G(d,p)//B3LYP/6-311G(d,p) and PCM-B3LYP/6-311+G(d,p)//B3LYP/6-311G(d,p) approaches, respectively. It was shown that the reaction of 2 (R=H) with PN proceeds via 2 + PN --> 2-PN --> 2-TS1 (O-O activation) --> 2-O(NO(2)(-)()) --> 2-SeO + NO(2)(-) pathway with a rate-determining barrier of 25.3 (14.8) kcal/mol at the NO(2)(-) dissociation step (numbers presented without parentheses are enthalpies, and those in parentheses are Gibbs free energies). The NO(3)(-) formation process, starting from the complex 2-O(NO(2)(-)()), requires by (7.9) kcal/mol more energy than the NO(2)(-) dissociation process and is unlikely to compete with the latter. Thus, in the gas phase, the peroxynitrite --> nitrate isomerization catalyzed by complex 2 is unlikely to occur. It is shown that the NO(3)(-) formation process is slightly more favorably than the NO(2)(-) dissociation process for complex 4, with a strongest electron-withdrawing ligand R=CF(3). Therefore, complex 4 (as well as complex 6 with R=OH) is predicted to be a good catalyst for peroxynitrite <--> nitrite isomerization in the gas phase. Solvent effects (a) change the rate-determining step of the reaction 2 + PN from NO(2)(-) dissociation in the gas phase to O-O activation, which occurs with barriers of (13.9), (8.4), (8.4), and (8.2) kcal/mol in water, dichloromethane, benzene, and cyclohexane, respectively, and (b) significantly reduce the NO(2)(-) dissociation energy, while only slightly destabilizing the NO(3)(-) formation barrier, and make the peroxynitrite <--> nitrate isomerization process practically impossible, even for complex 4.     

Falloff curves for the recombination reaction Cl + FC(O)O + M --> FC(O)OCl + M

The pressure dependence of the recombination reaction Cl + FC(O)O + M --> FC(O)OCl + M has been investigated at 296 K. FC(O)O radicals and Cl atoms were generated by laser flash photodissociation of FC(O)OO(O)CF at 193 nm in mixtures with Cl2 and He or SF6 over the total pressure range 8-645 Torr. The measured FC(O)O radical and F atom yields in the photolysis are 0.33 +/- 0.06 and 0.67 +/- 0.06. The reaction lies in the falloff range approaching the high-pressure limit. The extrapolations toward the limiting low- and high-pressure ranges were carried out using a reduced falloff curves formalism, which includes a recent implementation for the strong-collision broadening factors. The resulting values for the low-pressure rate coefficients are (2.2 +/- 0.4) x 10(-28)[He], (4.9 +/- 0.9) x 10(-28)[SF6], (1.9 +/- 0.3) x 10(-28)[Cl2] and (5.9 +/- 1.1) x 10(-28)[FC(O)OO(O)CF] cm3 molecule(-1) s(-1). The derived high-pressure rate coefficient is (4.4 +/- 0.8) x 10(-11) cm3 molecule(-1) s(-1). For the reaction Cl + FC(O)OCl --> Cl2 + FC(O)O a rate coefficient of (1.6 +/- 0.3) x 10(-11) cm3 molecule(-1) s(-1) was determined. The high-pressure rate coefficient was theoretically interpreted using SACM/CT calculations on an ab initio electronic potential computed at the G3S level of theory. Standard heat of formation values of -99.9 and -102.5 kcal mol(-1) were computed at the G3//B3LYP/6-311++G(3df,3pd) level of theory for cis-FC(O)OCl and trans-FC(O)OCl, respectively. The computed electronic barrier for the conversion between the trans and cis conformers is 8.9 kcal mol(-1). On the basis of the present results, the above reactions are expected to have a negligible impact on stratospheric ozone levels.     

Gas-phase identity nucleophilic substitution reactions of cyclopropenyl halides

The gas-phase identity nucleophilic substitution reactions of halide anions (X = F, Cl, and Br) with cyclopropenyl halides, X(-) + (CH)(3)X <= => X(CH)(3) + X(-), are investigated theoretically at four levels of theory, B3LYP/6-311+G**, MP2/6-311+G**, G2(+)MP2//MP2/6-311+G**, and G2(+)//MP2/6-311+G**. Four types of reaction paths, the sigma-attack S(N)2, pi-attack S(N)2'-syn, and S(N)2'-anti and sigmatropic 1,2-shift, are possible for all the halides. In the fluoride anion reactions, two types of stable adducts, syn- and anti-1,2-difluorocyclopropyl anions, can exist on the triple-well-type potential energy surface of the identity substitution reactions with rearrangement of double bond (C=C), S(N)2'-syn, and S(N)2'-anti processes. The TSs for the sigma-attack S(N)2 paths have "open" (loose) structures so that the ring positive charges are high rendering strong aromatic cyclopropenyl (delocalized) cation-like character. In contrast, in the pi-attack S(N)2' paths, a lone pair is formed at the unsubstituted carbon (C3), which stabilizes the 1,2-dihalocyclopropyl (delocalized) anion-like TS by two strong n(C)-sigma*(C-F) vicinal charge-transfer delocalization interactions. The barrier height increases in the order S(N)2'-anti < sigma-attack S(N)2 < S(N)2'-syn for X = Cl and Br, whereas for X = F the order is changed to S(N)2'-anti < S(N)2'-syn < sigma-attack S(N)2 due to the stable difluoro adduct formation. The sigmatropic 1,2-shift (circumambulatory) reactions have high activation barriers and cannot interfere with the substitution reactions.     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Modified Gaussian-2 level investigation of the identity ion-pair SN2 reactions of lithium halide and methyl halide with inversion and retention mechanisms

Identity ion-pair S(N)2 reactions LiX + CH(3)X --> XCH(3) + LiX (X = F, Cl, Br, and I) have been investigated in the gas phase and in solution at the level of the modified Gaussian-2 theory. Two possible reaction mechanisms, inversion and retention, are discussed. The reaction barriers relative to the complexes for the inversion mechanism [DeltaH(cent) ( not equal )(inv)] are found to be much higher than the corresponding values for the gas phase anionic S(N)2 reactions, decreasing in the following order: F (263.6 kJ mol(-1)) > Cl (203.3 kJ mol(-1)) > Br (174.7 kJ mol(-1)) > I (150.7 kJ mol(-1)). The barrier gaps between the two mechanisms [DeltaH(cent) ( not equal ) (ret) - DeltaH(cent) ( not equal ) (inv)] increase in the order F (-62.7 kJ mol(-1)) < Cl (4.4 kJ mol(-1)) < Br (24.9 kJ mol(-1)) < I (45.1 kJ mol(-1)). Thus, the retention mechanism is energetically favorable for fluorine and the inversion mechanism is favored for other halogens, in contrast to the anionic S(N)2 reactions at carbon where the inversion reaction channel is much more favorable for all of the halogens. The stabilization energies for the dipole-dipole complexes CH(3)X. LiX (DeltaH(comp)) are found to be similar for the entire set of systems with X = F, Cl, Br, and I, ranging from 53.4 kJ mol(-1) for I up to 58.9 kJ mol(-1) for F. The polarizable continuum model (PCM) has been used to evaluate the direct solvent effects on the energetics of the anionic and ion-pair S(N)2 reactions. The energetic profiles are found to be still double-well shaped for most of the ion-pair S(N)2 reactions in the solution, but the potential profile for reaction LiI + CH(3)I is predicted to be unimodal in the protic solvent. Good correlations between central barriers [DeltaH(cent) ( not equal ) (inv)] with the geometric looseness of the inversion transition state %C-X( not equal ), the dissociation energies of the C-X bond (D(C-X)) and Li-X bond (D(Li-X)) are observed, respectively.     

Corannulene as a Lewis base: computational modeling of protonation and lithium cation binding.

A computational modeling of the protonation of corannulene at B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p) and of the binding of lithium cations to corannulene at B3LYP/6-311G(d,p)//B3LYP/6-31G(d,p) has been performed. A proton attaches preferentially to one carbon atom, forming a sigma-complex. The isomer protonated at the innermost (hub) carbon has the best total energy. Protonation at the outermost (rim) carbon and at the intermediate (bridgehead rim) carbon is less favorable by ca. 2 and 14 kcal mol(-)(1), respectively. Hydrogen-bridged isomers are transition states between the sigma-complexes; the corresponding activation energies vary from 10 to 26 kcal mol(-)(1). With an empirical correction obtained from calculations on benzene, naphthalene, and azulene, the best estimate for the proton affinity of corannulene is 203 kcal mol(-)(1). The lithium cation positions itself preferentially over a ring. There is a small energetic preference for the 6-ring over the 5-ring binding (up to 2 kcal mol(-)(1)) and of the convex face over the concave face (3-5 kcal mol(-)(1)). The Li-bridged complexes are transition states between the pi-face complexes. Movement of the Li(+) cation over either face is facile, and the activation energy does not exceed 6 kcal mol(-)(1) on the convex face and 2.2 kcal mol(-)(1) on the concave face. In contrast, the transition of Li(+) around the corannulene edge involves a high activation barrier (24 kcal mol(-)(1) with respect to the lowest energy pi-face complex). An easier concave/convex transformation and vice versa is the bowl-to-bowl inversion with an activation energy of 7-12 kcal mol(-)(1). The computed binding energy of Li(+) to corannulene is 44 kcal mol(-)(1). Calculations of the (7)Li NMR chemical shifts and nuclear independent chemical shifts (NICS) have been performed to analyze the aromaticity of the corannulene rings and its changes upon protonation.     

Theoretical study of the equilibrium structure, vibrational spectrum, and thermochemistry of the peroxynitrate CF2BrCFBrOONO2

The results of a theoretical study of the molecular structure and conformational mobilities of the peroxynitrate CF(2)BrCFBrOONO(2) and its radical decomposition product CF(2)BrCFBrOO are reported in this paper. The most stable structures were calculated from ab initio G3(MP2)B3 and G4(MP2) methods and from density functional theory at the B3LYP/6-311+G(d) and B3LYP/6-311+G(3df) levels of theory. The equilibrium conformation of CF(2)BrCFBrOONO(2) indicates that the bromine atoms lie in position anti to each other and possess a COON dihedral angle of 114°. A quantum statistical analysis shows that about 40% of the internal rotors can freely rotate at room temperature. Our best values for the standard enthalpies of formation of CF(2)BrCFBrOONO(2) and CF(2)BrCFBrOO at 298 K obtained from isodesmic reactions at the G3(MP2)//B3LYP/6-311+G(3df) level of theory are -144.7 and -127.0 kcal mol(-1). From these values and the enthalpy of formation of the NO(2) radical, a CF(2)BrCFBrOO-NO(2) bond dissociation enthalpy of 26.0 ± 2 kcal mol(-1) was estimated.     

The X1 method for accurate and efficient prediction of heats of formation

We propose the X1 method which combines the density functional theory method with a neural network (NN) correction for an accurate yet efficient prediction of heats of formation. It calculates the final energy by using B3LYP6-311+G(3df,2p) at the B3LYP6-311+G(d,p) optimized geometry to obtain the B3LYP standard heats of formation at 298 K with the unscaled zero-point energy and thermal corrections at the latter basis set. The NN parameters cover 15 elements of H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, P, S, and Cl. The performance of X1 is close to the Gn theories, giving a mean absolute deviation of 1.43 kcalmol for the G399 set of 223 molecules up to 10 nonhydrogen atoms and 1.48 kcal/mol for the X107 set of 393 molecules up to 32 nonhydrogen atoms.     

Molecule-induced alkane homolysis with dioxiranes.

The mechanisms of C-H and C-C bond activations with dimethyldioxirane (DMD) were studied experimentally and computationally at the B3LYP/6-311+G**//B3LYP/6-31G* density functional theory level for the propellanes 3,6-dehydrohomoadamantane (2) and 1,3-dehydroadamantane (3). The sigma(C-C) activation of 3 with DMD (Delta G(*) = 23.9 kcal mol(-1) and Delta G(r) = -5.4 kcal mol(-1)) is the first example of a molecule-induced homolytic C-C bond cleavage. The C-H bond hydroxylation observed for 2 is highly exergonic (Delta G(r) = -74.4 kcal mol(-1)) and follows a concerted pathway (Delta G(*) = 34.8 kcal mol(-1)), in contrast to its endergonic molecule-induced homolysis (Delta G(*) = 28.8 kcal mol(-1) and Delta G(r) = +9.2 kcal mol(-1)). The reactivities of 2 and 3 with CrO(2)Cl(2), which follow a molecule-induced homolytic activation mechanism, parallel the DMD results only for highly reactive 3, but differ considerably for more stable propellanes such as 4-phenyl-3,6-dehydrohomoadamantane (1) and 2.     

Carbene Rearrangements Unsurpassed: Details of the C(7)H(6) Potential Energy Surface Revealed

The rearrangement of phenylcarbene (1) to 1,2,4,6-cycloheptatetraene (3) has been studied theoretically, using SCF, CASSCF, CASPT2N, DFT (B3LYP), CISD, CCSD, and CCSD(T) methods in conjunction with the 6-31G, 6-311+G, 6-311G(2d,p), cc-pVDZ, and DZd basis sets. Stationary points were characterized by vibrational frequency analyses at CASSCF(8,8)/6-31G and B3LYP/6-31G. Phenylcarbene (1) has a triplet ground state ((3)A") with a singlet-triplet separation (DeltaE(ST)) of 3-5 kcal mol(-)(1). In agreement with experiment, chiral 3 is the lowest lying structure on this part of the C(7)H(6) potential energy surface. Bicyclo[4.1.0]hepta-2,4,6-triene (2) is an intermediate in the rearrangement of 1 into 3, but it is unlikely to be observable experimentally due to a barrier height of only 1-2 kcal mol(-)(1). The enantiomers of 3 interconvert via the (1)A(2) state of cycloheptatrienylidene (4) with an activation energy of 20 kcal mol(-)(1). The "aromatic" (1)A(1) state, previously believed to be the lowest singlet state of 4, is roughly 10 kcal mol(-)(1) higher in energy than the (1)A(2) state, and, in violation of Hund's rule, (3)A(2) is also calculated to lie above (1)A(2) in energy. Thus, even if (3)A(2) were populated, it is likely to undergo rapid intersystem crossing to (1)A(2). We suggest (3)B(1)-4 is the metastable triplet observed by EPR.