scholarly journals DFT calculation of the Renner Teller Effect in NCO: Preliminary assessment of Exact Exchange Energy on the accuracy of the X2Π Renner Coefficient

Author(s):  
David Kashinski ◽  
Tyler Radziewicz ◽  
Matthew Suarez ◽  
Constantine Stephens ◽  
Edward Byrd

Assessment of DFT methods through analysis of the Renner-Teller Effect (RTE) in the XΠ state of the NCO radical was completed. Our results suggest that the amount of exact exchange at long range is important for an accurate description of the RTE in NCO. DFT functionals from the B3LYP, PBE, TPSS, M06, and M11 families with standard Correlation Consistent, 6-311G split valence family, as well as Sadlej, and Sapporo polarized triple-ζ basis sets were assessed. Our Renner coefficients are compared with previously published theoretical and experimental results to characterize the overall accuracy of various functional/basis set combinations in determining the RTE splitting in the Π (bending) modes of NCO(XΠ). We suggest that this method of analysis can be extended to other systems, serve as an accuracy metric when selecting a functional, and provide a means to create training sets for machine learning in computational molecular physics applications.

2018 ◽  
Author(s):  
Danilo Carmona ◽  
David Contreras ◽  
Oscar A. Douglas-Gallardo ◽  
Stefan Vogt-Geisse ◽  
Pablo Jaque ◽  
...  

The Fenton reaction plays a central role in many chemical and biological processes and has various applications as e.g. water remediation. The reaction consists of the iron-catalyzed homolytic cleavage of the oxygen-oxygen bond in the hydrogen peroxide molecule and the reduction of the hydroxyl radical. Here, we study these two elementary steps with high-level ab-initio calculations at the complete basis set limit and address the performance of different DFT methods following a specific classification based on the Jacob´s ladder in combination with various Pople's basis sets. Ab-initio calculations at the complete basis set limit are in agreement to experimental reference data and identified a significant contribution of the electron correlation energy to the bond dissociation energy (BDE) of the oxygen-oxygen bond in hydrogen peroxide and the electron affinity (EA) of the hydroxyl radical. The studied DFT methods were able to reproduce the ab-initio reference values, although no functional was particularly better for both reactions. The inclusion of HF exchange in the DFT functionals lead in most cases to larger deviations, which might be related to the poor description of the two reactions by the HF method. Considering the computational cost, DFT methods provide better BDE and EA values than HF and post--HF methods with an almost MP2 or CCSD level of accuracy. However, no systematic general prediction of the error based on the employed functional could be established and no systematic improvement with increasing the size in the Pople's basis set was found, although for BDE values certain systematic basis set dependence was observed. Moreover, the quality of the hydrogen peroxide, hydroxyl radical and hydroxyl anion structures obtained from these functionals was compared to experimental reference data. In general, bond lengths were well reproduced and the error in the angles were between one and two degrees with some systematic trend with the basis sets. From our results we conclude that DFT methods present a computationally less expensive alternative to describe the two elementary steps of the Fenton reaction. However, choice of approximated functionals and basis sets must be carefully done and the provided benchmark allows a systematic validation of the electronic structure method to be employed


1991 ◽  
Vol 69 (6) ◽  
pp. 927-933 ◽  
Author(s):  
Ted Schaefer ◽  
Rudy Sebastian ◽  
Christian Beaulieu

The inversion potentials, obtained from STO-3G, STO-3G(*), 3-21G, 3-21G(*), and 4-31G basis sets, are reported for thianthrene and thioxanthene, molecules in which both or only one of the methylene groups have been replaced by sulfur in 9,10-dihydroanthracene. Comparison with the available experimental data suggests that the split-valence bases lead to an overestimate, possibly by about 10 kJ/mol, of the inversion barrier in the crystal, whereas the STO-3G and STO-3G* basis sets underestimate this barrier. It appears that the inversion barrier for thianthrene is much lower in solution than in the crystal. The long-range coupling constants between the methylene and ring protons for thioxanthene in solution are consistent with an inversion barrier somewhat smaller than those obtained with the split-valence bases but rather larger than those predicted with the STO-3G basis set. The bond lengths and angles in the equilibrium structures of the two molecules, as computed with the 3-21G(*) basis, agree reasonably well with those in their crystals, except that the theoretical folding angles are smaller than measured. These discrepancies become less marked when expectation values are calculated from the theoretical inversion potentials at finite temperatures. Key words: MO calculations, inversion potentials of thianthrene and thioxanthene; 1H NMR, thioxanthene; spin–spin coupling constants, long range, in thioxanthene.


2013 ◽  
Vol 91 (7) ◽  
pp. 559-572 ◽  
Author(s):  
Jennifer L. Kellie ◽  
Stacey D. Wetmore

When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.


2019 ◽  
Vol 2019 ◽  
pp. 1-23 ◽  
Author(s):  
Kenneth Irving ◽  
Martina Kieninger ◽  
Oscar N. Ventura

The performance of a group of density functional methods of progressive complexity for the description of the ClO bond in a series of chlorine oxides was investigated. The simplest ClO radical species and the two isomeric structures XClO/ClOX for each X = H, Cl, and O were studied using the PW91, TPSS, B3LYP, PBE0, M06, M06-2X, BMK, and B2PLYP functionals. Geometry optimizations and reaction enthalpies and enthalpies of formation for each species were calculated using Pople basis sets and the (aug)-cc-pVnZ Dunning sets, with n = D, T, Q, 5, and 6. For the calculation of enthalpies of formation, atomization and isodesmic reactions were employed. Both the precision of the methods with respect to the increase of the basis sets, as well as their accuracy, were gauged by comparing the results with the more accurate CCSD(T) calculations, performed using the same basis sets as for the DFT methods. The results obtained employing composite chemical methods (G4, CBS-QB3, and W1BD) were also used for the comparisons, as well as the experimental results when they are available. The results obtained show that error compensation is the key for successful description of molecular properties (geometries and energies) by carefully selecting the method and basis sets. In general, expansion of the one-electron basis set to the limit of completeness does not improve results at the DFT level, but just the opposite. The enthalpies of formation calculated at the CCSD(T)/aug-cc-pV6Z for the species considered are generally in agreement with experimental determinations and the most accurate theoretical values. Different sources of error in the calculations are discussed in detail.


1976 ◽  
Vol 29 (8) ◽  
pp. 1635 ◽  
Author(s):  
L Radom

Ab initio molecular orbital theory with the minimal STO-3G and split-valence 4-31G basis sets is used to obtain geometries of 18 anions:OH-, NH2-, HF2-, BH4-, BF4-, C22-, CN-, NCN2-, N3-, NO2-, NO3-, 0CCO2-, CO32-, HCOO-, CH3COO-, C2O42-, C4O42- and C(CN)3-. The theoretical results are compared with experimental results from the literature. The STO-3G basis set performs somewhat worse for anions than for neutral molecules. On the other hand, the 4-31G basis set gives good results and predicts bond lengths to within 0.02� for all the molecules considered. Limited information on bond angle predictions suggests that these are of comparable quality to those for neutral molecules. The tricyanomethanide ion is predicted to be planar.


2018 ◽  
Author(s):  
Danilo Carmona ◽  
David Contreras ◽  
Oscar A. Douglas-Gallardo ◽  
Stefan Vogt-Geisse ◽  
Pablo Jaque ◽  
...  

The Fenton reaction plays a central role in many chemical and biological processes and has various applications as e.g. water remediation. The reaction consists of the iron-catalyzed homolytic cleavage of the oxygen-oxygen bond in the hydrogen peroxide molecule and the reduction of the hydroxyl radical. Here, we study these two elementary steps with high-level ab-initio calculations at the complete basis set limit and address the performance of different DFT methods following a specific classification based on the Jacob´s ladder in combination with various Pople's basis sets. Ab-initio calculations at the complete basis set limit are in agreement to experimental reference data and identified a significant contribution of the electron correlation energy to the bond dissociation energy (BDE) of the oxygen-oxygen bond in hydrogen peroxide and the electron affinity (EA) of the hydroxyl radical. The studied DFT methods were able to reproduce the ab-initio reference values, although no functional was particularly better for both reactions. The inclusion of HF exchange in the DFT functionals lead in most cases to larger deviations, which might be related to the poor description of the two reactions by the HF method. Considering the computational cost, DFT methods provide better BDE and EA values than HF and post--HF methods with an almost MP2 or CCSD level of accuracy. However, no systematic general prediction of the error based on the employed functional could be established and no systematic improvement with increasing the size in the Pople's basis set was found, although for BDE values certain systematic basis set dependence was observed. Moreover, the quality of the hydrogen peroxide, hydroxyl radical and hydroxyl anion structures obtained from these functionals was compared to experimental reference data. In general, bond lengths were well reproduced and the error in the angles were between one and two degrees with some systematic trend with the basis sets. From our results we conclude that DFT methods present a computationally less expensive alternative to describe the two elementary steps of the Fenton reaction. However, choice of approximated functionals and basis sets must be carefully done and the provided benchmark allows a systematic validation of the electronic structure method to be employed


2018 ◽  
Author(s):  
Oscar Ventura ◽  
Kenneth Irving ◽  
Martina Kieninger

<p>The performance of a group of density functional methods of progressive complexity for the description of the ClO bond in a series of chlorine oxides was investigated. The simplest ClO radical species as well as the two isomeric structures XClO/ClOX for each X=H, Cl and O were studied using the PW91, TPSS, B3LYP, PBE0, M06, M06-2X, BMK and B2PLYP functionals. Geometry optimizations as well as reaction enthalpies and enthalpies of formation for each species were calculated using Pople basis sets and the (aug)-cc-pVnZ Dunning sets, with n=2-6. For the calculation of enthalpies of formation, atomization as well as isodesmic reactions were employed. Both the precision of the methods with respect to the increase of the basis sets, as well as their accuracy, were gauged by comparing the results with the more accurate CCSD(T) calculations, performed using the same basis sets as for the DFT methods. The results obtained employing composite chemical methods (G4, CBS-QB3 and W1BD) were also used for the comparisons, as well as the experimental results when they are available. The results obtained show that error compensation is the key for successful description of molecular properties (geometries and energies) by carefully selecting method and basis sets. In general, expansion of the one-electron basis set to the limit of completeness does not improve results at the DFT level, but just the opposite. The enthalpies of formation calculated at the CCSD(T)/aug-cc-pV6Z for the species considered are generally in agreement with experimental determinations, and the most accurate derived theoretically up to present. Different sources of error in the calculations are discussed in detail.</p>


2017 ◽  
Author(s):  
Saurav Dutta ◽  
Bhabani S. Mallik

<div> <table> <tr> <td> <p>Knowledge of the transition state is crucial in determining the mechanism in order to diversify the applicability of the reaction. The computational method is the most convenient way to locate the transition state in the absence any efficient experimental technique. We have applied the method of the transition state search on the Diels-Alder reaction computationally by means of combined linear synchronous transit and quadratic synchronous transit methods. Here we have shown that, of various methods adopted, BOP functional with numerical basis set provides a computationally economical alternative to the widely used B3LYP functional with higher Gaussian basis sets in the transition state search. It can reproduce the experimental parameters like activation energy of the Diels-Alder reaction, and the calculations are much faster than the corresponding other functional based calculations.</p> </td> </tr> </table> </div>


Sign in / Sign up

Export Citation Format

Share Document