Rozen's epoxidation reagent, CH3CN.HOF: a theoretical study of its structure, vibrational spectroscopy, and reaction mechanism

Rozen's epoxidation reagent, CH(3)CN.HOF, and a prototype epoxidation reaction employing it, have been subjected to an extensive ab initio and density functional study. Its anharmonic force field reveals a very strong red shift for the OH stretch and a strong blue shift for the HOF bend, in semiquantitative agreement with experiment. The very strong hydrogen bond (8.20 kcal/mol at the W1 level) not only serves to stabilize the reactant but also considerably lowers the barrier height for epoxidation of ethylene. Moreover, the reaction byproduct HF is found to act autocatalytically. The OH moiety acquires HO(+) character in the transition state. Our W1 benchmark data for the reaction profile allow the performance of various DFT functionals to be assessed. In general, "kinetics" functionals overestimate barrier heights, the BMK functional less so than the others. The B1B95 and TPSS33B95 meta-GGA functionals both perform very well, whereas general-purpose hybrid GGAs underestimate barrier heights. The simple PBE0 functional does reasonably well.     

Testing the performance of density functionals for the calculation of energetic properties of complex-forming radical-molecule reactions

The performance of some popular and some more recent density functional methods for the calculation of energies of stationary points on the potential surfaces of radical-molecule reactions was examined. The functionals studied are B3-LYP, BH&H, BH&H-LYP, MPW1K, MPWB1K, TPSS, TPSSh, BB1K, M05 and M05-2X, in conjunction with nine different AO basis sets. The reaction energies, barrier heights and the relative energies of the pre-and post-reaction complexes were compared with those obtained at the CCSD(T)/CBS limit for the reactions of OH radicals with HOOH and CH₃OOH. Very poor barrier heights are provided by the B3-LYP, TPSS and TPSSh functionals. The best overall performance was obtained with the BB1K, MPW1K and MPWB1K functionals. In these reactions all of the studied functionals provide converged results only if they are used with large basis sets like aug-cc-pVTZ and def2-TZVP. The data show that before relying on a functional for a specific reaction, it is desirable to make some test calculations on the performance. The same functional can predict some relative energies very well and some others very poorly even in systems including chemically similar reactants.     

Double-hybrid functionals for thermochemical kinetics

We propose two new double-hybrid functionals, denoted B2K-PLYP and mPW2K-PLYP, which yield thermochemical performance comparable to existing double-hybrid functionals but offer superior performance for barrier heights of various kinds. We show that the new functionals yield excellent performance for all of the following: (a) main-group thermochemistry; (b) main-group thermochemical kinetics; (c) late transition metal reactions. In addition, B2K-PLYP performs well for weak interactions.     

Comparative assessment of density functional methods for 3d transition-metal chemistry

In the present study, we comparatively assessed the newly developed M05 functional against a data set of reaction energies for transition-metal chemistry. The functionals to which we compare are BLYP, B3LYP, B97-2, MPWLYP1M, TPSS, and TPSSh. We draw the following conclusions: (1) TPSS gives the best performance for calculating the binding energies of three transition-metal dimers (Sc(2), Ni(2), and V(2)) that have severe multireference character, (2) B97-2 gives the best performance for calculating the binding energies of the nine metal-ligand diatomics (three monohydrides, three monoxide, and three monofluorides), and (3) M05 gives the overall best performance for all 18 data in the assessment, and it has a mean unsigned error 55% lower than the popular B3LYP functional. Since the M05 functional also gives good performance for main-group thermochemistry, for noncovalent chemistry, and for calculating barrier heights, M05 can be applied to a wide range of problems where nonhybrid functionals or functionals designed for kinetics fail.     

DFT and ab initio calculations on two reactions between hydrogen atoms and the fire suppressants 2-H heptafluoropropane and CF3Br

Reaction enthalpies and barrier heights of the reactions CF3Br+H-->CF3+HBr {reaction (1)} and CF3CHFCF3+H-->CF3CFCF3+H2 {reaction (2)} have been calculated at the near state-of-the-art ab initio level, and also by employing the B3LYP, BH&HLYP, BB1K, MPW1K, MPWB1K and TPSS1KCIS functionals. In addition, the integrated molecular orbital+molecular orbital (IMOMO) method has been used to study reaction (2). The ab initio benchmark values of the reaction enthalpy (298 K) and barrier height (0 K) of reaction (2) are reported for the first time {-(0.7+/-0.7) and 13.3+/-0.5 kcal/mole respectively}. When density functional theory (DFT) results are compared with ab initio benchmarks for both reactions (1) and (2), the MPWB1K functional is found to have the best performance of the six functionals used. The IMOMO method with the RCCSD/aug-cc-pVTZ and/or RCCSD(T)/aug-cc-pVTZ levels, as the high levels of calculation on the model system, gives reaction enthalpies and barrier heights of reaction (2), which agree with ab initio benchmark values to within 1 kcal/mole. Computed key geometrical parameters and imaginary vibrational frequencies of the transition state structures of reactions (1) and (2) obtained at different levels of calculation are compared. The magnitudes of the computed imaginary vibrational frequencies of the transition states of both reactions considered are found to be very sensitive to the levels of calculation used to obtain them. The heat of formation (298 K) of CF3CFCF3 calculated at the near state-of-the-art level has a value of -(318+/-3) kcal/mole.     

Performance of density functionals for calculation of reductive ring-opening reaction energies of Li<Superscript>+</Superscript>-EC and Li<Superscript>+</Superscript>-VC

Reaction energies were determined for reductive ring-opening reactions of Li⁺-coordinated ethylene carbonate (EC) and vinylene carbonate (VC) by using various density functional theory (DFT) and ab-initio methods applying the basis sets up to Dunnings aug-cc-pVQZ. The methods examined include the local density functional (SVWN), the pure gradient-corrected density functionals (BLYP and BPW91), and the hybrid density functionals (B3LYP, B1LYP, B3PW91, and mPW1PW91). Comparison of the DFT results with ab-initio results indicates that the mPW1PW91 approach introduced by Adamo and Barone, is superior to all the other DFT methods (including B3LYP). The performance of more cost-effective Pople-type basis sets ranging from 6-31G(d,p) to 6-311++G(3df,3pd) was assessed at DFT levels of theory by calibrating them with the aug-cc-pVQZ results     

Design of density functionals that are broadly accurate for thermochemistry, thermochemical kinetics, and nonbonded interactions

This paper develops two new hybrid meta exchange-correlation functionals for thermochemistry, thermochemical kinetics, and nonbonded interactions. The new functionals are called PW6B95 (6-parameter functional based on Perdew-Wang-91 exchange and Becke-95 correlation) and PWB6K (6-parameter functional for kinetics based on Perdew-Wang-91 exchange and Becke-95 correlation). The resulting methods were comparatively assessed against the MGAE109/3 main group atomization energy database, against the IP13/3 ionization potential database, against the EA13/3 electron affinity database, against the HTBH38/4 and NHTBH38/04 hydrogen-transfer and non-hydrogen-transfer barrier height databases, against the HB6/04 hydrogen bonding database, against the CT7/04 charge-transfer complex database, against the DI6/04 dipole interaction database, against the WI7/05 weak interaction database, and against the new PPS5/05 pi-pi stacking interaction database. From the assessment and comparison of methods, we draw the following conclusions, based on an analysis of mean unsigned errors: (i) The PW6B95, MPW1B95, B98, B97-1, and TPSS1KCIS methods give the best results for a combination of thermochemistry and nonbonded interactions. (ii) PWB6K, MPWB1K, BB1K, MPW1K, and MPW1B95 give the best results for a combination of thermochemical kinetics and nonbonded interactions. (iii) PWB6K outperforms the MP2 method for nonbonded interactions. (iv) PW6B95 gives errors for main group covalent bond energies that are only 0.41 kcal (as measured by mean unsigned error per bond (MUEPB) for the MGAE109 database), as compared to 0.56 kcal/mol for the second best method and 0.92 kcal/mol for B3LYP.     

Benchmark database of barrier heights for heavy atom transfer, nucleophilic substitution, association, and unimolecular reactions and its use to test theoretical methods

A benchmark database of forward and reverse barrier heights for 19 non-hydrogen-transfer reactions has been developed by using Weizmann 1 calculations, and 29 DFT methods and 6 ab initio wave-function theory (WFT) methods have been tested against the new database as well as against an older database for hydrogen atom transfer reactions. Among the tested hybrid DFT methods without kinetic energy density, MPW1K is the most accurate model for calculations of barrier heights. Among the tested hybrid meta DFT methods, BB1K and MPWB1K are the two most accurate models for the calculations of barrier heights. Overall, the results show that BB1K and MPWB1K are the two best DFT methods for calculating barrier heights, followed in order by MPW1K, MPWKCIS1K, B1B95, MPW1B95, BH and HLYP, B97-2, mPW1PW91, and B98. The popular B3LYP method has a mean unsigned error four times larger than that of BB1K. Of the methods tested, QCISD(T) is the best ab initio WFT method for barrier height calculations, and QCISD is second best, but QCISD is outperformed by the BB1K, MPWB1K, MPWKCIS1K, and MPW1K methods.     

一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

DFT studies of the phosphinidene derivative HPNaF and its insertion reaction with R–H(R=F,OH, NH<Subscript>2</Subscript>, CH<Subscript>3</Subscript>)

The formations of the phosphinidene derivative HPNaF and its insertion reactions with R–H (R=F, OH, NH₂, CH₃) have been systematically investigated employing the density functional theory (DFT), such as the B3LYP and MPW1PW91 methods. A comparison with the results of MP2 calculations shows that MPW1PW91 underestimates the barrier heights for the four reactions considered. Similarly, the same is also true for the B3LYP method depending on the selected reactions, but by much less than MPW1PW91, where the barrier heights of the four reactions are 25.2, 85.7, 119.0, and 142.4 kJ/mol at the B3LYP/6-311+G* level of theory, respectively. All the mechanisms of the four reactions are identical to each other, i.e., an intermediate has been located during the insertion reaction. Then, the intermediate could dissociate to substituted phosphinidane(H₂RP) and NaF with a barrier corresponding to their respective dissociation energies. Correspondingly, the reaction energies for the four reactions are −92.2, −68.1, −57.2, and −44.3 kJ/mol at the B3LYP/6-311+G* level of theory, respectively, where both the B3LYP and MPW1PW91 methods underestimate the reaction energies compared with the MP2 results. The linear correlations between the calculated barrier heights and the reaction energies have also been observed. As a result, the relative reactivity among the four insertion reactions should be as follows: H–F > H–OH > H–NH₂ > H–CH₃.     

Assessment of theoretical procedures for calculating barrier heights for a diverse set of water-catalyzed proton-transfer reactions

Accurate electronic barrier heights are obtained for a set of nine proton-transfer tautomerization reactions, which are either (i) uncatalyzed, (ii) catalyzed by one water molecule, or (iii) catalyzed by two water molecules. The barrier heights for reactions (i) and (ii) are obtained by means of the high-level ab initio W2.2 thermochemical protocol, while those for reaction (iii) are obtained using the W1 protocol. These three sets of benchmark barrier heights allow an assessment of the performance of more approximate theoretical procedures for the calculation of barrier heights of uncatalyzed and water-catalyzed reactions. We evaluate initially the performance of the composite G4 procedure and variants thereof (e.g., G4(MP2) and G4(MP2)-6X), as well as that of standard ab initio procedures (e.g., MP2, SCS-MP2, and MP4). We find that the performance of the G4(MP2)-type thermochemical procedures deteriorates with the number of water molecules involved in the catalysis. This behavior is linked to deficiencies in the MP2-based basis-set-correction term in the G4(MP2)-type procedures. This is remedied in the MP4-based G4 procedure, which shows good performance for both the uncatalyzed and the water-catalyzed reactions, with mean absolute deviations (MADs) from the benchmark values lying below the threshold of "chemical accuracy" (arbitrarily defined as 1 kcal mol(-1) ≈ 4.2 kJ mol(-1)). We also examine the performance of a large number of density functional theory (DFT) and double-hybrid DFT (DHDFT) procedures. We find that, with few exceptions (most notably PW6-B95 and B97-2), the performance of the DFT procedures that give good results for the uncatalyzed reactions deteriorates with the number of water molecules involved in the catalysis. The DHDFT procedures, on the other hand, show excellent performance for both the uncatalyzed and catalyzed reactions. Specifically, almost all of them afford MADs below the "chemical accuracy" threshold, with ROB2-PLYP and B2K-PLYP showing the best overall performance.     

Benchmark Study of the Performance of Density Functional Theory for Bond Activations with (Ni,Pd)-Based Transition-Metal Catalysts

The performance of 23 density functionals, including one LDA, four GGAs, three meta-GGAs, three hybrid GGAs, eight hybrid meta-GGAs, and ten double-hybrid functionals, was investigated for the computation of activation energies of various covalent main-group single bonds by four catalysts: Pd, PdCl⁻, PdCl₂, and Ni (all in the singlet state). A reactant complex, the barrier, and reaction energy were considered, leading to 164 energy data points for statistical analysis. Extended Gaussian AO basis sets were used in all calculations. The best functional for the complete benchmark set relative to estimated CCSD(T)/CBS reference data is PBE0-D3, with an MAD value of 1.1 kcal mol⁻¹ followed by PW6B95-D3, the double hybrid PWPB95-D3, and B3LYP-D3 (1.9 kcal mol⁻¹ each). The other tested hybrid meta-GGAs perform less well (M06-HF: 7.0 kcal mol⁻¹; M06-2X: 6.3 kcal mol⁻¹; M06: 4.9 kcal mol⁻¹) for the investigated reactions. In the Ni case, some double hybrids show larger errors due to partial breakdown of the perturbative treatment for the correlation energy in cases with difficult electronic structures (partial multi-reference character). Only double hybrids either with very low amounts of perturbative correlation (e.g., PBE0-DH) or that use the opposite-spin correlation component only (e.g., PWPB95) seem to be more robust. We also investigated the effect of the D3 dispersion correction. While the barriers are not affected by this correction, significant and mostly positive results were observed for reaction energies. Furthermore, six very recently proposed double-hybrid functionals were analyzed regarding the influence of the amount of Fock exchange as well as the type of perturbative correlation treatment. According to these results, double hybrids with <50–60 % of exact exchange and ∼30 % perturbative correlation perform best.     

Benchmark of density functional theory methods on the prediction of bond energies and bond distances of noble-gas containing molecules

We have tested three pure density functional theory (DFT) functionals, BLYP, MPWPW91, MPWB95, and ten hybrid DFT functionals, B3LYP, B3P86, B98, MPW1B95, MPW1PW91, BMK, M05-2X, M06-2X, B2GP-PLYP, and DSD-BLYP with a series of commonly used basis sets on the performance of predicting the bond energies and bond distances of 31 small neutral noble-gas containing molecules. The reference structures were obtained using the CCSD(T)∕aug-cc-pVTZ theory and the reference energies were based on the calculation at the CCSD(T)∕CBS level. While in general the hybrid functionals performed significantly better than the pure functionals, our tests showed a range of performance by these hybrid functionals. For the bond energies, the MPW1B95∕6-311+G(2df,2pd), BMK∕aug-cc-pVTZ, B2GP-PLYP∕aug-cc-pVTZ, and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 2.0-2.3 kcal∕mol per molecule. For the bond distances, the MPW1B95∕6-311+G(2df,2pd), MPW1PW91∕6-311+G(2df,2pd), and B3P86∕6-311+G(2df,2pd), DSD-BLYP∕6-311+G(2df,2pd), and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 0.008-0.013 A? per bond. The current study showed that a careful selection of DFT functionals is very important in the study of noble-gas chemistry, and the most recommended methods are MPW1B95∕6-311+G(2df,2pd) and DSD-BLYP∕aug-cc-pVTZ.     

甲烷与HO2反应的微观动力学研究

应用从头算方法和变分过渡态理论,在B3LYP/6-311+G^**方法下和300～2000K温度范围内研究甲烷与HO₂反应的微观动力学特性,得到由过渡态向反应物方向、向产物方向的能垒分别是11.83和102.90kJ/mol,理论计算正向反应速率常数与实验值之比为1.08～2.85,用此方法还可以预测没有实验数据的温度点反应的速率常数.     

Hybrid density functional theory with specific reaction parameter: hydrogen abstraction reaction of fluoromethane by the hydroxyl radical

Three potential energy surfaces with specific reaction parameters are developed and tested for the OH + CH(3)F --> H(2)O + CH(2)F reaction. The goal of this work is to determine surfaces that provide calculated reaction rate constants that are comparable to the experimental data. The potential energy surfaces are constructed using hybrid and hybrid meta density functional theory methods, and the levels of electronic structure theory used in this study are mPW1PW91, B1B95, and mPW1B95 in conjunction with the 6-31+G(d,p) basis set. The reaction rate constants are calculated over the range 200-1,500 K using variational transition state theory with multidimensional tunneling contributions. New specific-reaction-parameter Hartree-Fock contributions are determined, and the hybrid density functional theory methods with these new contributions (35.5 +/- 1.2% for mPW1PW91, 36.6 +/- 1.2% for B1B95, and 40.7 +/- 1.2% for mPW1B95, respectively) reproduce accurate rate constants over an extended temperature range. On these potential energy surfaces, the classical barrier height for the hydrogen abstraction reaction is determined to be 3.4-3.8 kcal/mol, with a best estimate value of 3.6 kcal/mol.     

Highly accurate first-principles benchmark data sets for the parametrization and validation of density functional and other approximate methods. Derivation of a robust, generally applicable, double-hybrid functional for thermochemistry and thermochemical kinetics

We present a number of near-exact, nonrelativistic, Born-Oppenheimer reference data sets for the parametrization of more approximate methods (such as DFT functionals). The data were obtained by means of the W4 ab initio computational thermochemistry protocol, which has a 95% confidence interval well below 1 kJ/mol. Our data sets include W4-08, which are total atomization energies of over 100 small molecules that cover varying degrees of nondynamical correlations, and DBH24-W4, which are W4 theory values for Truhlar's set of 24 representative barrier heights. The usual procedure of comparing calculated DFT values with experimental atomization energies is hampered by comparatively large experimental uncertainties in many experimental values and compounds errors due to deficiencies in the DFT functional with those resulting from neglect of relativity and finite nuclear mass. Comparison with accurate, explicitly nonrelativistic, ab initio data avoids these issues. We then proceed to explore the performance of B2x-PLYP-type double hybrid functionals for atomization energies and barrier heights. We find that the optimum hybrids for hydrogen-transfer reactions, heavy-atoms transfers, nucleophilic substitutions, and unimolecular and recombination reactions are quite different from one another: out of these subsets, the heavy-atom transfer reactions are by far the most sensitive to the percentages of Hartree-Fock-type exchange y and MP2-type correlation x in an (x, y) double hybrid. The (42,72) hybrid B2K-PLYP, as reported in a preliminary communication, represents the best compromise between thermochemistry and hydrogen-transfer barriers, while also yielding excellent performance for nucleophilic substitutions. By optimizing for best overall performance on both thermochemistry and the DBH24-W4 data set, however, we find a new (36,65) hybrid which we term B2GP-PLYP. At a slight expense in performance for hydrogen-transfer barrier heights and nucleophilic substitutions, we obtain substantially better performance for the other reaction types. Although both B2K-PLYP and B2GP-PLYP are capable of 2 kcal/mol quality thermochemistry, B2GP-PLYP appears to be the more robust toward nondynamical correlation and strongly polar character. We additionally find that double-hybrid functionals display excellent performance for such problems as hydrogen bonding, prototype late transition metal reactions, pericyclic reactions, prototype cumulene-polyacetylene system, and weak interactions.     

Is the uniform electron gas limit important for small Ag clusters? Assessment of different density functionals for Ag(n) (n < or = 4)

Twenty-three density functional theory (DFT) methods, including the second- and the third-generation functionals, are tested in conjunction with two basis sets (LANL2DZ and SDD) for studying the properties of neutral and ionic silver clusters. We find that DFT methods incorporating the uniform electron gas limit in the correlation functional, namely, those with Perdew's correlation functionals (PW91, PBE, P86, and TPSS), Becke's B95, and the Van Voorhis-Scuseria functional VSXC, generally perform better than the other group of functionals, e.g., those incorporating the LYP correlation functional and variations of the B97 functional. Strikingly, these two groups of functionals can produce qualitatively different results for the Ag3 and Ag4 clusters. The energetic properties and vibrational frequencies of Ag(n) are also evaluated by the different functionals. The present study shows that the choice of DFT methods for heavy metals may be critical. It is found that the exact-exchange-incorporated PBE functional (PBE1PBE) is among the best for predicting the range of properties.     

Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09

A benchmark study on all possible density functional theory (DFT) methods in Gaussian09 is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (E_int) for small non‐covalently interacting molecular dimers in “dispersion‐dominated” (class 1), “dipole‐induced dipole” (class 2), and “dipole‐dipole” (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and E_int was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good E_int, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional “recommended” for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the “recommended” list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and E_int. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in Gaussian09 to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.     

Ab initio benchmark calculations on Ca(II) complexes and assessment of density functional theory methodologies

A set of benchmark results for the geometries, binding energies, and protonation affinities of 24 complexes of small organic ligands with Ca(II) is provided. The chosen level of theory is CCSD(T)/CBS obtained by means of a composite procedure. The performance of four density functionals, namely, PW91, PBE, B3LYP, and TPSS and several Pople-type basis sets, namely, 6-31G(d), 6-31+G(d), 6-31+G(2d,p) and 6-311+G(d) have been assessed. Additionally, the nature of the metal ligand bonding has been analyzed by means of the Symmetry Adapted Perturbation Theory (SAPT). We have found that the B3LYP hybrid functional, in conjunction with either the polarized double-ζ 6-31+G(2d,p) or the triple-ζ 6-311+G(d) basis sets, yields the closest results compared to the benchmark data. The SAPT analysis stresses the importance of induction effects in the binding of these complexes and suggests that consideration of classical electrostatic contributions alone may not be reliable enough for the prediction of relative binding energies for Ca(II) complexes.