On the Reaction of FNO2 with CH3, tert-Butyl, and C13H21

Solvent Effect on the Potential Energy Surfaces of the F- + CH3CH2Br Reaction

Although substantial work has been undertaken on reaction pathways involved in base-promoted elimination reactions and bimolecular nucleophilic substitution reaction of F^- on CH₃CH₂X (X = Cl, Br, I), the effect of solvents with varying dielectric constants on the stereochemistry of each of the reaction species involved across the reaction profile have not yet been clearly understood. The present investigation reports the effect of solvents on the potential energy surfaces (PES) and structures of the species appearing in the reaction pathway of F^- with bromoethane. The PESs in the gas phase have been computed at MP2 level and CCSD(T) level. The performance of several hybrid density functional, such as B3LYP, M06, M06L, BHandH, X3LYP, M05, M05-2X, and M06-2X have also been investigated toward describing the elimination and nucleophilic substitution reactions. With respect to MAE values and to make the computation cost-effective, we have explored the implicit continuum solvent model, CPCM in solvents like cyclohexane, methanol, acetonitrile, dimethyl sulfoxide and water. The reactant complexes proceed through the subsequent steps to produce fluoroethane as the substitution product and ethylene as one of the elimination products. For elimination reaction both syn and anti elimination have been explored. The calculated relatives energies values, which are negative in the gas phase, are found to be positive in polar solvents since the point charge in the separated reactants are more stabilized than the dispersed charge in the transient complex, which has also been analyzed through NBO analysis.

E2 and SN2 Reactions of X(-) + CH3CH2X (X = F, Cl); an ab Initio and DFT Benchmark Study

We have computed consistent benchmark potential energy surfaces (PESs) for the anti-E2, syn-E2, and SN2 pathways of X(-) + CH3CH2X with X = F and Cl. This benchmark has been used to evaluate the performance of 31 popular density functionals, covering local-density approximation, generalized gradient approximation (GGA), meta-GGA, and hybrid density-functional theory (DFT). The ab initio benchmark has been obtained by exploring the PESs using a hierarchical series of ab initio methods [up to CCSD(T)] in combination with a hierarchical series of Gaussian-type basis sets (up to aug-cc-pVQZ). Our best CCSD(T) estimates show that the overall barriers for the various pathways increase in the order anti-E2 (X = F) < SN2 (X = F) < SN2 (X = Cl) ∼ syn-E2 (X = F) < anti-E2 (X = Cl) < syn-E2 (X = Cl). Thus, anti-E2 dominates for F(-) + CH3CH2F, and SN2 dominates for Cl(-) + CH3CH2Cl, while syn-E2 is in all cases the least favorable pathway. Best overall agreement with our ab initio benchmark is obtained by representatives from each of the three categories of functionals, GGA, meta-GGA, and hybrid DFT, with mean absolute errors in, for example, central barriers of 4.3 (OPBE), 2.2 (M06-L), and 2.0 kcal/mol (M06), respectively. Importantly, the hybrid functional BHandH and the meta-GGA M06-L yield incorrect trends and qualitative features of the PESs (in particular, an erroneous preference for SN2 over the anti-E2 in the case of F(-) + CH3CH2F) even though they are among the best functionals as measured by their small mean absolute errors of 3.3 and 2.2 kcal/mol in reaction barriers. OLYP and B3LYP have somewhat higher mean absolute errors in central barriers (5.6 and 4.8 kcal/mol, respectively), but the error distribution is somewhat more uniform, and as a consequence, the correct trends are reproduced.

Oxidative Addition of the Fluoromethane C−F Bond to Pd. An ab Initio Benchmark and DFT Validation Study

We have computed a state-of-the-art benchmark potential energy surface (PES) for two reaction pathways (oxidative insertion, OxIn, and S(N)2) for oxidative addition of the fluoromethane C-F bond to the palladium atom and have used this to evaluate the performance of 26 popular density functionals, covering LDA, GGA, meta-GGA, and hybrid density functionals, for describing these reactions. The ab initio benchmark is obtained by exploring the PES using a hierarchical series of ab initio methods (HF, MP2, CCSD, CCSD(T)) in combination with a hierarchical series of seven Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account through a full four-component all-electron approach. Our best estimate of kinetic and thermodynamic parameters is -5.3 (-6.1) kcal/mol for the formation of the reactant complex, 27.8 (25.4) kcal/mol for the activation energy for oxidative insertion (OxIn) relative to the separate reactants, 37.5 (31.8) kcal/mol for the activation energy for the alternative S(N)2 pathway, and -6.4 (-7.8) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). Our work highlights the importance of sufficient higher angular momentum polarization functions for correctly describing metal-d-electron correlation. Best overall agreement with our ab initio benchmark is obtained by functionals from all three categories, GGA, meta-GGA, and hybrid DFT, with mean absolute errors of 1.4-2.7 kcal/mol and errors in activation energies ranging from 0.3 to 2.8 kcal/mol. The B3LYP functional compares very well with a slight underestimation of the overall barrier for OxIn by -0.9 kcal/mol. For comparison, the well-known BLYP functional underestimates the overall barrier by -10.1 kcal/mol. The relative performance of these two functionals is inverted with respect to previous findings for the insertion of Pd into the C-H and C-C bonds. However, all major functionals yield correct trends and qualitative features of the PES, in particular, a clear preference for the OxIn over the alternative S(N)2 pathway.

Oxidative addition of the ethane CC bond to Pd. Anab initiobenchmark and DFT validation study

We have computed a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the ethane C-C bond to the palladium atom and have used this to evaluate the performance of 24 popular density functionals, covering LDA, GGA, meta-GGA, and hybrid density functionals, for describing this reaction. The ab initio benchmark is obtained by exploring the PES using a hierarchical series of ab initio methods [HF, MP2, CCSD, CCSD(T)] in combination with a hierarchical series of five Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Our best estimate of kinetic and thermodynamic parameters is -10.8 (-11.3) kcal/mol for the formation of the reactant complex, 19.4 (17.1) kcal/mol for the activation energy relative to the separate reactants, and -4.5 (-6.8) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). Our work highlights the importance of sufficient higher angular momentum polarization functions for correctly describing metal-d-electron correlation. Best overall agreement with our ab initio benchmark is obtained by functionals from all three categories, GGA, meta-GGA, and hybrid DFT, with mean absolute errors of 1.5 to 2.5 kcal/mol and errors in activation energies ranging from -0.2 to -3.2 kcal/mol. Interestingly, the well-known BLYP functional compares very reasonably with a slight underestimation of the overall barrier by -0.9 kcal/mol. For comparison, with B3LYP we arrive at an overestimation of the overall barrier by 5.8 kcal/mol. On the other hand, B3LYP performs excellently for the central barrier (i.e., relative to the reactant complex) which it underestimates by only -0.1 kcal/mol.

Ab initio and DFT benchmark study for nucleophilic substitution at carbon (SN2@C) and silicon (SN2@Si)

To obtain a set of consistent benchmark potential energy surfaces (PES) for the two archetypal nucleophilic substitution reactions of the chloride anion at carbon in chloromethane (S(N)2@C) and at silicon in chlorosilane (S(N)2@Si), we have explored these PESes using a hierarchical series of ab initio methods [HF, MP2, MP4SDQ, CCSD, CCSD(T)] in combination with a hierarchical series of six Gaussian-type basis sets, up to g polarization. Relative energies of stationary points are converged to within 0.01 to 0.56 kcal/mol as a function of the basis-set size. Our best estimate, at CCSD(T)/aug-cc-pVQZ, for the relative energies of the [Cl(-), CH(3)Cl] reactant complex, the [Cl-CH(3)-Cl](-) transition state and the stable [Cl-SiH(3)-Cl](-) transition complex is -10.42, +2.52, and -27.10 kcal/mol, respectively. Furthermore, we have investigated the performance for these reactions of four popular density functionals, namely, BP86, BLYP, B3LYP, and OLYP, in combination with a large doubly polarized Slater-type basis set of triple-zeta quality (TZ2P). Best overall agreement with our CCSD(T)/aug-cc-pVQZ benchmark is obtained with OLYP and B3LYP. However, OLYP performs better for the S(N)2@C overall and central barriers, which it underestimates by 2.65 and 4.05 kcal/mol, respectively. The other DFT approaches underestimate these barriers by some 4.8 (B3LYP) to 9.0 kcal/mol (BLYP).

Computational exploration of rearrangements related to the vitamin B12-dependent ethanolamine ammonia lyase catalyzed transformation

DFT (B3LYP/6-31G) and ab initio molecular orbital theory (QCISD/cc-pVDZ) are used to investigate several possible mechanisms involving free radical intermediates as well as their protonated forms for processes related to the coenzyme B(12)-dependent rearrangement catalyzed by ethanolamine ammonia lyase. Two major types of rearrangements are discussed in detail, intramolecular migration and dissociation of the amine/ammonia groups, for both of which several scenarios are considered. According to the calculations, the complete dissociation of the migrating group and its subsequent association constitute an unlikely route for both the protonated and the unprotonated reactant because of the high-energy barriers (more than 23 kcal/mol) involved in these steps. Direct migration of the protonated amine group is far more favorable (10.4 kcal/mol) and therefore presents the most likely candidate for the actual enzymatic reaction. The calculations further imply that the direct loss of an ammonium cation (10.6 kcal/mol) represents a feasible pathway as well. Comparing the rearrangements for the aminoethanol radical and its protonated counterpart, in line with previous findings reported by Golding, Radom, and co-workers, we find that the migration of a protonated group is in general associated with lower energy barriers, suggesting that the actual enzyme substrate quite likely corresponds to (partially) protonated aminoethanol. As the extent of the substrate protonation/deprotonation by the active site of the enzyme may vary, the actual energy barriers are expected to range between the values calculated for the two extreme cases of a substrate, that is, the aminoethanol radical 2 and its fully protonated form 6.