A comparative examination of density functional performance against the ISOL24/11 isomerization energy benchmark

2016 ◽  
Vol 1090 ◽  
pp. 147-152 ◽  
Author(s):  
Sierra Rayne ◽  
Kaya Forest
2016 ◽  
Vol 94 (10) ◽  
pp. 865-869 ◽  
Author(s):  
Sierra Rayne ◽  
Kaya Forest

The photoisomerization equilibrium between the model tannins (-)-catechin and (-)-epicatechin in aqueous solution was investigated at the density functional level of theory to gain insights into the action of these compounds as natural sunscreens in aquatic systems. Increasing water temperature, as might be expected on seasonal and diurnal bases, is predicted to shift the equilibrium further in favor of catechin. The isomerization energy between catechin and epicatechin was also considered in a range of polar protic, polar aprotic, apolar protic, and apolar aprotic solvents using the solvation model based on density (SMD) and integral equation formalism polarizable continuum model (IEFPCM). The IEFPCM yielded a modest range in isomerization energies depending on solvent polarity or proticity, whereas a substantial variation was observed with the SMD. The SMD results suggest that the solvation environment around catechin and epicatechin will play a major role on the photoisomerization equilibrium between these two compounds. As the freely dissolved monomer in aquatic systems, the catechin–epicatechin photoisomerization equilibrium will be in the range of 11:1 to 14:1. In the less polar environments of associations with dissolved organic matter or within a larger tannin structural framework, the theoretical modeling efforts indicate that the catechin–epicatechin photoisomerization equilibrium could be as low as 3:1.


2019 ◽  
Author(s):  
Ernest Awoonor-Williams ◽  
William Isley ◽  
Stephen Dale ◽  
Erin Johnson ◽  
Haibo Yu ◽  
...  

Targeted covalent inhibitor drugs require computational methods that go beyond simple molecular-mechanical force fields in order to model the chemical reactions that occur when they bind to their targets. Here, several semi-empirical and density-functional theory (DFT) methods are assessed for their ability to describe the potential energy surface and reaction energies of the covalent modification of a thiol by an electrophile. Functionals such as PBE and B3LYP fail to predict a stable enolate intermediate. This is largely due to delocalization error, which spuriously stabilizes the pre-reaction complex, in which excess electron density is transferred from the thiolate to the electrophile. Functionals with a high-exact exchange component, range-separated DFT functionals, and variationally-optimized exact exchange (i.e., the LC-B05minV functional) correct this issue to various degrees. The large gradient behaviour of the exchange enhancement factor is also found to significantly affect the results, leading to the improved performance of PBE0. While ωB97X-D and M06-2X were easonably accurate, no method provided quantitative accuracy for all three electrophiles, making this a very strenuous test of functional performance. Additionally, one drawback of M06-2X was that MD simulations using this functional were only stable if a fine integration grid was used. The low-cost semi-empirical methods, PM3, AM1, and PM7, provide a qualitatively correct description of the reaction mechanism, although the energetics are not quantitatively reliable. As a proof of concept, the potential of mean force for the addition of methylthiolate to MVK was calculated using QM/MM MD in an explicit polarizable aqueous solvent.<br>


2020 ◽  
Vol 153 (21) ◽  
pp. 214116
Author(s):  
Subrata Jana ◽  
Abhilash Patra ◽  
Szymon Śmiga ◽  
Lucian A. Constantin ◽  
Prasanjit Samal

Author(s):  
Ernest Awoonor-Williams ◽  
William Isley ◽  
Stephen Dale ◽  
Erin Johnson ◽  
Haibo Yu ◽  
...  

Targeted covalent inhibitor drugs require computational methods that go beyond simple molecular-mechanical force fields in order to model the chemical reactions that occur when they bind to their targets. Here, several semi-empirical and density-functional theory (DFT) methods are assessed for their ability to describe the potential energy surface and reaction energies of the covalent modification of a thiol by an electrophile. Functionals such as PBE and B3LYP fail to predict a stable enolate intermediate. This is largely due to delocalization error, which spuriously stabilizes the pre-reaction complex, in which excess electron density is transferred from the thiolate to the electrophile. Functionals with a high-exact exchange component, range-separated DFT functionals, and variationally-optimized exact exchange (i.e., the LC-B05minV functional) correct this issue to various degrees. The large gradient behaviour of the exchange enhancement factor is also found to significantly affect the results, leading to the improved performance of PBE0. While ωB97X-D and M06-2X were easonably accurate, no method provided quantitative accuracy for all three electrophiles, making this a very strenuous test of functional performance. Additionally, one drawback of M06-2X was that MD simulations using this functional were only stable if a fine integration grid was used. The low-cost semi-empirical methods, PM3, AM1, and PM7, provide a qualitatively correct description of the reaction mechanism, although the energetics are not quantitatively reliable. As a proof of concept, the potential of mean force for the addition of methylthiolate to MVK was calculated using QM/MM MD in an explicit polarizable aqueous solvent.<br>


2011 ◽  
Vol 130 (4-6) ◽  
pp. 851-857 ◽  
Author(s):  
Jong-Won Song ◽  
Takao Tsuneda ◽  
Takeshi Sato ◽  
Kimihiko Hirao

2021 ◽  
Author(s):  
Sebastian Spicher ◽  
Eike Caldeweyher ◽  
Andreas Hansen ◽  
Stefan Grimme

<p>The strongly attractive noncovalent interactions of charged atoms or molecules with p-systems are important binding motifs in many chemical and biological systems. These so-called ion-pi interactions play a major role in enzymes, molecular recognition, and for the structure of proteins. In this work, a molecular test set termed IONPI19 is compiled for inter- and intramolecular ion-pi interactions, which is well balanced between anionic and cationic systems. The IONPI19 set includes interaction energies of significantly larger molecules (up to 133 atoms) than in other ion-pi test sets and covers a broad range of binding motifs. Accurate (local) coupled cluster values are provided as reference. Overall, 18 density functional approximations, including seven (meta-)GGAs, seven hybrid functionals, and four double hybrid functionals combined with three different London dispersion corrections, are benchmarked for interaction energies. DFT results are further compared to wave function based methods such as MP2 and dispersion corrected Hartree-Fock. Also the performance</p><p>of semiempirical QM methods such as the GFNn-xTB and PMx family of methods is tested. It is shown that dispersion-uncorrected DFT underestimates ion-pi interactions significantly, even though electrostatic interactions dominate the overall binding. Accordingly, the new charge dependent D4 dispersion model is found to be consistently better than the standard D3 correction. Furthermore, the functional performance trend along Jacob’s ladder is generally obeyed and the reduction of the self-interaction error leads to an improvement of (double) hybrid functionals over (meta-)GGAs, even though the effect of the SIE is smaller than expected. Overall, the double hybrids PWPB95-D4/QZ and revDSD-PBEP86-D4/QZ turned out to be the most reliable among all assessed methods in predicting ion-pi interactions, which opens up new perspectives for systems where coupled cluster calculations are no longer computationally feasible.</p>


2021 ◽  
Author(s):  
Sebastian Spicher ◽  
Eike Caldeweyher ◽  
Andreas Hansen ◽  
Stefan Grimme

<p>The strongly attractive noncovalent interactions of charged atoms or molecules with p-systems are important binding motifs in many chemical and biological systems. These so-called ion-pi interactions play a major role in enzymes, molecular recognition, and for the structure of proteins. In this work, a molecular test set termed IONPI19 is compiled for inter- and intramolecular ion-pi interactions, which is well balanced between anionic and cationic systems. The IONPI19 set includes interaction energies of significantly larger molecules (up to 133 atoms) than in other ion-pi test sets and covers a broad range of binding motifs. Accurate (local) coupled cluster values are provided as reference. Overall, 18 density functional approximations, including seven (meta-)GGAs, seven hybrid functionals, and four double hybrid functionals combined with three different London dispersion corrections, are benchmarked for interaction energies. DFT results are further compared to wave function based methods such as MP2 and dispersion corrected Hartree-Fock. Also the performance</p><p>of semiempirical QM methods such as the GFNn-xTB and PMx family of methods is tested. It is shown that dispersion-uncorrected DFT underestimates ion-pi interactions significantly, even though electrostatic interactions dominate the overall binding. Accordingly, the new charge dependent D4 dispersion model is found to be consistently better than the standard D3 correction. Furthermore, the functional performance trend along Jacob’s ladder is generally obeyed and the reduction of the self-interaction error leads to an improvement of (double) hybrid functionals over (meta-)GGAs, even though the effect of the SIE is smaller than expected. Overall, the double hybrids PWPB95-D4/QZ and revDSD-PBEP86-D4/QZ turned out to be the most reliable among all assessed methods in predicting ion-pi interactions, which opens up new perspectives for systems where coupled cluster calculations are no longer computationally feasible.</p>


2019 ◽  
Author(s):  
Ernest Awoonor-Williams ◽  
William Isley ◽  
Stephen Dale ◽  
Erin Johnson ◽  
Haibo Yu ◽  
...  

Targeted covalent inhibitor drugs require computational methods that go beyond simple molecular-mechanical force fields in order to model the chemical reactions that occur when they bind to their targets. Here, several semi-empirical and density-functional theory (DFT) methods are assessed for their ability to describe the potential energy surface and reaction energies of the covalent modification of a thiol by an electrophile. Functionals such as PBE and B3LYP fail to predict a stable enolate intermediate. This is largely due to delocalization error, which spuriously stabilizes the pre-reaction complex, in which excess electron density is transferred from the thiolate to the electrophile. Functionals with a high-exact exchange component, range-separated DFT functionals, and variationally-optimized exact exchange (i.e., the LC-B05minV functional) correct this issue to various degrees. The large gradient behaviour of the exchange enhancement factor is also found to significantly affect the results, leading to the improved performance of PBE0. While ωB97X-D and M06-2X were easonably accurate, no method provided quantitative accuracy for all three electrophiles, making this a very strenuous test of functional performance. Additionally, one drawback of M06-2X was that MD simulations using this functional were only stable if a fine integration grid was used. The low-cost semi-empirical methods, PM3, AM1, and PM7, provide a qualitatively correct description of the reaction mechanism, although the energetics are not quantitatively reliable. As a proof of concept, the potential of mean force for the addition of methylthiolate to MVK was calculated using QM/MM MD in an explicit polarizable aqueous solvent.<br>


2019 ◽  
Vol 21 (44) ◽  
pp. 24478-24488 ◽  
Author(s):  
Martin Gleditzsch ◽  
Marc Jäger ◽  
Lukáš F. Pašteka ◽  
Armin Shayeghi ◽  
Rolf Schäfer

In depth analysis of doping effects on the geometric and electronic structure of tin clusters via electric beam deflection, numerical trajectory simulations and density functional theory.


2000 ◽  
Vol 98 (20) ◽  
pp. 1639-1658 ◽  
Author(s):  
Yuan He, Jurgen Grafenstein, Elfi Kraka,

Sign in / Sign up

Export Citation Format

Share Document