scholarly journals Benchmarking quantum mechanical methods for calculating reaction energies of reactions catalyzed by enzymes

2020 ◽  
Vol 2 ◽  
pp. e8 ◽  
Author(s):  
Jitnapa Sirirak ◽  
Narin Lawan ◽  
Marc W. Van der Kamp ◽  
Jeremy N. Harvey ◽  
Adrian J. Mulholland
Keyword(s):  
Density Functional ◽  
Good Choice ◽  
Basis Set ◽  
Basis Sets ◽  
Quantum Mechanical ◽  
Spin Component ◽  
Energy Calculation ◽  
Mp2 Calculations ◽  
Mechanical Methods ◽  
Reaction Energies

To assess the accuracy of different quantum mechanical methods for biochemical modeling, the reaction energies of 20 small model reactions (chosen to represent chemical steps catalyzed by commonly studied enzymes) were calculated. The methods tested included several popular Density Functional Theory (DFT) functionals, second-order Møller Plesset perturbation theory (MP2) and its spin-component scaled variant (SCS-MP2), and coupled cluster singles and doubles and perturbative triples (CCSD(T)). Different basis sets were tested. CCSD(T)/aug-cc-pVTZ results for all 20 reactions were used to benchmark the other methods. It was found that MP2 and SCS-MP2 reaction energy calculation results are similar in quality to CCSD(T) (mean absolute error (MAE) of 1.2 and 1.3 kcal mol−1, respectively). MP2 calculations gave a large error in one case, and are more subject to basis set effects, so in general SCS-MP2 calculations are a good choice when CCSD(T) calculations are not feasible. Results with different DFT functionals were of reasonably good quality (MAEs of 2.5–5.1 kcal mol−1), whereas popular semi-empirical methods (AM1, PM3, SCC-DFTB) gave much larger errors (MAEs of 11.6–14.6 kcal mol−1). These results should be useful in guiding methodological choices and assessing the accuracy of QM/MM calculations on enzyme-catalyzed reactions.

DFT calculation of crystallographic properties of dioctahedral 2:1 phyllosilicates

Clay Minerals ◽  
2008 ◽  
Vol 43 (3) ◽  
pp. 351-361 ◽  
Author(s):  
J. Ortega-Castro ◽  
N. Hernández-Haro ◽  
A. Hernández-Laguna ◽  
C. I. Sainz-Díaz
Keyword(s):  
Density Functional ◽  
Basis Sets ◽  
Quantum Mechanical ◽  
Functional Theory ◽  
Oh Groups ◽  
Low Degree ◽  
Mechanical Methods ◽  
Weak Interlayer ◽  
Low Charge ◽  
Cation Substitutions

AbstractThe low-charge dioctahedral 2:1 phyllosilicates are an important group of clay minerals that have a low degree of cation substitution and very weak interlayer interatomic interactions which are difficult to reproduce with quantum mechanical calculations. In order to study the crystallographic properties of these compounds with density functional theory (DFT) quantum-mechanical methods, an optimization of norm-conserving pseudopotentials of Al, Si, O, H and Na atoms has been carried out, and an optimization of the cutoff radii of the basis sets has been accomplished. Crystallographic properties and vibrational stretching frequencies of the OH groups, ν(OH), have been calculated, being consistent with previous computational and experimental results. All frequencies can be related to the different molecular environment of the OH groups. The effect of octahedral Fe3+ substitution on the ν(OH) frequency is reproduced. Several configurations of cation substitutions and interlayer cation (IC) positions are studied in low-charge dioctahedral 2:1 phyllosilicates, such as Al4(Si7–xAlx)O20(OH)4Nax, with x = 0.25, 0.50 and 0.75, indicating that the IVAl3+ is highly dispersed and the IC tends to be in the substituted ditrigonal hole. For the Al4(Si7Al)O20(OH)4Na composition, the trans-vacant form is more stable than the cis-vacant one.

Benchmark of Simplified Time-Dependent Density Functional Theory for UV-Vis Spectral Properties of Porphyrinoids

2019 ◽  
Author(s):  
Kamal Batra ◽  
Stefan Zahn ◽  
Thomas Heine
Keyword(s):  
Density Functional Theory ◽  
Density Functional ◽  
Large Scale ◽  
Absolute Error ◽  
Time Dependent ◽  
Basis Set ◽  
Basis Sets ◽  
Functional Theory ◽  
Framework Materials ◽  
Dependent Density

<p>We thoroughly benchmark time-dependent density- functional theory for the predictive calculation of UV/Vis spectra of porphyrin derivatives. With the aim to provide an approach that is computationally feasible for large-scale applications such as biological systems or molecular framework materials, albeit performing with high accuracy for the Q-bands, we compare the results given by various computational protocols, including basis sets, density-functionals (including gradient corrected local functionals, hybrids, double hybrids and range-separated functionals), and various variants of time-dependent density-functional theory, including the simplified Tamm-Dancoff approximation. An excellent choice for these calculations is the range-separated functional CAM-B3LYP in combination with the simplified Tamm-Dancoff approximation and a basis set of double-ζ quality def2-SVP (mean absolute error [MAE] of ~0.05 eV). This is not surpassed by more expensive approaches, not even by double hybrid functionals, and solely systematic excitation energy scaling slightly improves the results (MAE ~0.04 eV). </p>

AlN Nanotubes: A DFT Study of Al-27 and N-14 Electric Field Gradient Tensors

2007 ◽  
Vol 62 (12) ◽  
pp. 711-715 ◽  
Author(s):  
Ahmad Seif ◽  
Mahmoud Mirzaei ◽  
Mehran Aghaie ◽  
Asadollah Boshra
Keyword(s):  
Electric Field ◽  
Electric Field Gradient ◽  
Field Gradient ◽  
Density Functional ◽  
Quadrupole Coupling Constant ◽  
Basis Set ◽  
Basis Sets ◽  
B3lyp Method ◽  
Package Program ◽  
Nqr Parameters

Density functional theory (DFT) calculations were performed to calculate the electric field gradient (EFG) tensors at the sites of aliminium (27Al) and nitrogen (14N) nuclei in an 1 nm of length (6,0) single-walled aliminium nitride nanotube (AlNNT) in three forms of the tubes, i. e. hydrogencapped, aliminium-terminated and nitrogen-terminated as representatives of zigzag AlNNTs. At first, each form was optimized at the level of the Becke3,Lee-Yang-Parr (B3LYP) method, 6-311G∗∗ basis set. After, the EFG tensors were calculated at the level of the B3LYP method, 6-311++G∗∗ and individual gauge for localized orbitals (IGLO-II and IGLO-III) types of basis sets in each of the three optimized forms and were converted to experimentally measurable nuclear quadrupole resonance (NQR) parameters, i. e. quadrupole coupling constant (qcc) and asymmetry parameter (ηQ). The evaluated NQR parameters revealed that the considered model of AlNNT can be divided into four equivalent layers with similar electrostatic properties.With the exception of Al-1, all of the three other Al layers have almost the same properties, however, N layers show significant differences in the magnitudes of the NQR parameters in the length of the nanotube. Furthermore, the evaluated NQR parameters of Al-1 in the Al-terminated form and N-1 in the N-terminated form revealed the different roles of Al (base agent) and of N (acid agent) in AlNNT. All the calculations were carried out using the GAUSSIAN 98 package program.

REEXAMINATION OF LOW ENERGY STRUCTURES OF ${\rm Au}_{4}^{-}$ AND Au4

2010 ◽  
Vol 09 (supp01) ◽  
pp. 1-7 ◽  
Author(s):  
YI GAO ◽  
YU ZHAO ◽  
X. C. ZENG
Keyword(s):  
Free Energy ◽  
Density Functional ◽  
Photoelectron Spectrum ◽  
Global Minimum ◽  
Low Energy ◽  
Basis Sets ◽  
Cluster Method ◽  
Energy Calculation ◽  
Couple Cluster ◽  
B3lyp Method

Low energy isomers of [Formula: see text] and Au4 were reexamined using the hybrid density functional B3LYP method and the couple-cluster method with the aug-cc-pVDZ-PP and aug-cc-pVTZ-PP basis sets. For [Formula: see text], the B3LYP method favors the zigzag isomer and the second order Moller–Plesset perturbation (MP2) total energy calculation favors the D2h rhombus isomer, whereas the couple-cluster singles and doubles with perturbative triples [CCSD(T)] level of theory favors the Y-shaped C2v isomer. The pyramid isomer is much higher in energy and could be easily excluded. The Gibbs free energy correction based on harmonic approximation suggests that the zigzag isomer is lower in free energy than the D2h rhombus isomer at 298.15 K. These results confirm that the Y-shaped C2v isomer is the global minimum at both 0 K and room temperature and is thus the major isomer to account for the experimental photoelectron spectrum. The zigzag isomer is suggested, as a minor isomer, to account for the weak second peak at 3.40 eV in the experimental photoelectron spectrum. For neutral Au4 , the zigzag isomer is more stable than D2h rhombus isomer at the B3LYP level and the D2h rhombus isomer is the global minimum on basis of all post Hartree–Fock levels of theory.

scholarly journals Theoretical Study of Solvent Effects on the Electronic and Thermodynamic Properties of Tetrathiafulvalene (TTF) Molecule Based on DFT

2021 ◽  
pp. 42-54
Author(s):  
Rabiu Nuhu Muhammad ◽  
N. M. Mahraz ◽  
A. S Gidado ◽  
A. Musa
Keyword(s):  
Thermodynamic Properties ◽  
Bond Length ◽  
Density Functional ◽  
Theoretical Study ◽  
Energy Gap ◽  
Basis Set ◽  
Basis Sets ◽  
Frontier Molecular Orbital ◽  
Orbital Energy ◽  
Chemical Hardness

Tetrathiafulvalene () is an organosulfur compound used in the production of molecular devices such as switches, sensors, nonlinear optical devices and rectifiers. In this work, a theoretical study on the effects of solvent on TTF molecule was investigated and reported based on Density Functional Theory (DFT) as implemented in Gaussian 03 package using B3LYP/6-31++G(d,p) basis set. Different solvents were introduced as a bridge to investigate their effects on the electronic structure. The HUMO, LUMO, energy gap, global chemical index, thermodynamic properties, NLO and DOS analysis of the TTF molecule in order to determine the reactivity and stability of the molecule were obtained. The results obtained showed that the solvents have effects on the electronic and non-linear-optical properties of the molecule. The optimized bond length revealed that the molecule has strong bond in gas phase with smallest bond length of about 1.0834Å than in the rest of the solvents. It was observed that the molecule is more stable in acetonitrile with HOMO-LUMO gap and chemical hardness of 3.6373eV and 1.8187eV respectively. This indicates that the energy gap and chemical hardness of TTF molecule increases with the increase in polarity and dielectric constant of the solvents. The computed results agreed with the results in the literature. The thermodynamics and NLO properties calculation also indicated that TTF molecule has highest value of specific heat capacity (Cv), total dipole moment () and first order hyperpolarizability () in acetonitrile, while acetone has the highest value of entropy and toluene has a slightly higher value of zero point vibrational energy (ZPVE) than the rest of the solvents. The results show that careful selection of the solvents and basis sets can tune the frontier molecular orbital energy gap of the molecule and can be used for molecular device applications.

A density functional approach toward structural features and properties of C20…N2X2 (X = H, F, Cl, Br, Me) molecules

2014 ◽  
Vol 13 (04) ◽  
pp. 1450023 ◽  
Author(s):  
Reza Ghiasi ◽  
Morteza Zaman Fashami ◽  
Amir Hossein Hakimioun
Keyword(s):  
Density Functional ◽  
Chemical Shifts ◽  
Geometry Optimization ◽  
Structural Features ◽  
Nbo Analysis ◽  
Basis Set ◽  
Basis Sets ◽  
Basis Set Superposition Error ◽  
Homo And Lumo ◽  
Lowest Unoccupied Molecular Orbitals

In this work, the interaction of C 20 with N 2 X 2 ( X = H , F , Cl , Br , Me ) molecules has been explored using the B3LYP, M062x methods and 6-311G(d,p) and 6-311+G(d,p) basis sets. The interaction energies (IEs) obtained with standard method were corrected by basis set superposition error (BSSE) during the geometry optimization for all molecules at the same levels of theory. It was found C 20… N 2 H 2 interaction is stronger than the interaction of other N 2 X 2 ( X = F , Cl , Br , Me ) with C 20. Highest occupied and lowest unoccupied molecular orbitals (HOMO and LUMO, respectively) levels are illustrated by density of states spectra (DOS). The nucleus-independent chemical shifts (NICSs) confirm that C 20… N 2 X 2 molecules exhibit aromatic characteristics. Geometries obtained from DFT calculations were used to perform NBO analysis. Also, 14 N NQR parameters of the C 20… N 2 X 2 molecules are predicted.

scholarly journals Double-Hybrid DFT Functionals for the Condensed Phase: Gaussian and Plane Waves Implementation and Evaluation

Molecules ◽  
2020 ◽  
Vol 25 (21) ◽  
pp. 5174
Author(s):  
Frederick Stein ◽  
Jürg Hutter ◽  
Vladimir V. Rybkin
Keyword(s):  
Condensed Phase ◽  
Density Functional ◽  
Plane Waves ◽  
Basis Set ◽  
Basis Sets ◽  
Density Functionals ◽  
Basis Set Superposition Error ◽  
Super Cell ◽  
Long Range Interactions ◽  
Function Correlation

Intermolecular interactions play an important role for the understanding of catalysis, biochemistry and pharmacy. Double-hybrid density functionals (DHDFs) combine the proper treatment of short-range interactions of common density functionals with the correct description of long-range interactions of wave-function correlation methods. Up to now, there are only a few benchmark studies available examining the performance of DHDFs in condensed phase. We studied the performance of a small but diverse selection of DHDFs implemented within Gaussian and plane waves formalism on cohesive energies of four representative dispersion interaction dominated crystal structures. We found that the PWRB95 and ωB97X-2 functionals provide an excellent description of long-ranged interactions in solids. In addition, we identified numerical issues due to the extreme grid dependence of the underlying density functional for PWRB95. The basis set superposition error (BSSE) and convergence with respect to the super cell size are discussed for two different large basis sets.

scholarly journals Basis Set Effects in the Description of the Cl-O Bond in ClO and XClO/ClOX Isomers (X = H, O, and Cl) Using DFT and CCSD(T) Methods

2019 ◽  
Vol 2019 ◽  
pp. 1-23 ◽  
Author(s):  
Kenneth Irving ◽  
Martina Kieninger ◽  
Oscar N. Ventura
Keyword(s):  
Density Functional ◽  
Basis Set ◽  
Basis Sets ◽  
Enthalpies Of Formation ◽  
Isodesmic Reactions ◽  
Dft Methods ◽  
Functional Methods ◽  
The One ◽  
Isomeric Structures ◽  
Different Sources

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.

THEORETICAL STUDIES ON BOND DISSOCIATION ENERGIES FOR SOME THIOL COMPOUNDS BY DENSITRY FUNCTIONAL THEORY AND CBS-Q METHOD

2008 ◽  
Vol 07 (05) ◽  
pp. 943-951 ◽  
Author(s):  
XIAO-HONG LI ◽  
ZHENG-XIN TANG ◽  
ABRAHAM F. JALBOUT ◽  
XIAN-ZHOU ZHANG ◽  
XIN-LU CHENG
Keyword(s):  
Density Functional ◽  
Computational Cost ◽  
Basis Set ◽  
Basis Sets ◽  
Q Method ◽  
Bond Dissociation Energies ◽  
Functional Theory ◽  
Thiol Compounds ◽  
Bond Dissociation ◽  
Dissociation Energies

Quantum chemical calculations are used to estimate the bond dissociation energies (BDEs) for 15 thiol compounds. These compounds are studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86, PBE0) methods and the complete basis set (CBS-Q) method together with 6-311G** basis set. It is demonstrated that B3P86 and CBS-Q methods are accurate for computing the reliable BDEs for thiol compounds. In order to test whether the non-local BLYP method suggested by Fu et al.19 is general for our study and whether B3P86 method has a low basis set sensitivity, the BDEs for seven thiol compounds are also calculated using BLYP/6-31+G* and B3P86 method with 6-31+G*, 6-31+G**, and 6-311+G** basis sets for comparison. The obtained results are compared with the available experimental results. It is noted that B3P86 method is not sensitive to the basis set. Considering the inevitable computational cost of CBS-Q method and the reliability of the B3P86 calculations, B3P86 method with a moderate or a larger basis set may be more suitable to calculate the BDEs of the C–SH bond for thiol compounds.

Sign in / Sign up

Export Citation Format

Share Document