1
|
Hancock AC, Giudici E, Goerigk L. How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models. J Comput Chem 2024; 45:1667-1681. [PMID: 38553847 DOI: 10.1002/jcc.27351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 03/07/2024] [Accepted: 03/12/2024] [Indexed: 06/04/2024]
Abstract
Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- ω PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.
Collapse
Affiliation(s)
- Amy C Hancock
- School of Chemistry, The University of Melbourne, Parkville, Victoria, Australia
| | - Erica Giudici
- School of Chemistry, The University of Melbourne, Parkville, Victoria, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria, Australia
| |
Collapse
|
2
|
Mun H, Lorpaiboon W, Ho J. In Search of the Best Low-Cost Methods for Efficient Screening of Conformers. J Phys Chem A 2024; 128:4391-4400. [PMID: 38754085 DOI: 10.1021/acs.jpca.4c01407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
Locating the lowest energy conformer is crucial for the accurate computation of equilibrium properties of molecular systems. This paper examines the performance of efficient low-cost methods in terms of the alignment and relative energies of their energy minima against the benchmark revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ potential energy surface. The low-cost methods considered include GFN-FF, GFN2-xTB, DFTB3, HF-3c, B97-3c, PBEh-3c, and r2SCAN-3c composite methods against a diverse test set of 20 compounds including alkanes, perfluoroalkyl molecules, peptides, open-shell radicals, and Zn(II) complexes of varying sizes. The "3c" composite methods are generally more accurate, but are at least 2-3 orders of magnitude more expensive than tight-binding methods which have energy minima that align well with the benchmark potential energy surface. The findings of this paper were further exploited to introduce a simple strategy involving Grimme's CENSO energy-sorting algorithm that resulted in up to an order of magnitude reduction in computational time for locating the lowest energy conformer on the revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ surface.
Collapse
Affiliation(s)
- Haedam Mun
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Wanutcha Lorpaiboon
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
3
|
Ferrero R, Pantaleone S, Gho CI, Hoti G, Trotta F, Brunella V, Corno M. Unveiling the synergy: a combined experimental and theoretical study of β-cyclodextrin with melatonin. J Mater Chem B 2024; 12:4004-4017. [PMID: 38568714 DOI: 10.1039/d3tb02795c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
Melatonin (MT) is a vital hormone controlling biorhythms, and optimizing its release in the human body is crucial. To address MT's unfavorable pharmacokinetics, we explored the inclusion complexes of MT with β-cyclodextrin (β-CD). Nano spray drying was applied to efficiently synthesize these complexes in three molar ratios (MT : β-CD = 1 : 1, 2 : 1, and 1 : 2), reducing reagent use and expediting inclusion. The complex powders were characterized through thermal analyses (TGA and DSC), Fourier transform infrared spectroscopy (FTIR), and in vitro MT release measurements via high-performance liquid chromatography (HPLC). In parallel, computational studies were conducted, examining the stability of MT : β-CD complexes by means of unbiased semi-empirical conformational searches refined by DFT, which produced a distribution of MT : β-CD binding enthalpies. Computational findings highlighted that these complexes are stabilized by specific hydrogen bonds and non-specific dispersive forces, with stronger binding in the 1 : 1 complex, which was corroborated by in vitro release data. Furthermore, the alignment between simulated and experimental FTIR spectra demonstrated the quality of both the structural model and computational methodology, which was crucial to enhance our comprehension of optimizing MT's release for therapeutic applications.
Collapse
Affiliation(s)
- Riccardo Ferrero
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| | - Stefano Pantaleone
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| | - Cecilia Irene Gho
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| | - Gjylije Hoti
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| | - Francesco Trotta
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| | - Valentina Brunella
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| | - Marta Corno
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, Via P. Giuria 7, 10125 Torino, Italy.
| |
Collapse
|
4
|
Khatun M, Paul S, Roy S, Dey S, Anoop A. Performance of Density Functionals and Semiempirical 3c Methods for Small Gold-Thiolate Clusters. J Phys Chem A 2023; 127:2242-2257. [PMID: 36877153 DOI: 10.1021/acs.jpca.2c07561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
In light of the recent surge in computational studies of gold thiolate clusters, we present a comparison of popular density functionals (DFAs) and three-part corrected methods (3c-methods) on their performance by taking a data set named as AuSR18 consisting of 18 isomers of Aun(SCH3)m (m ≤ n = 1-3). We have compared the efficiency and accuracy of the DFAs and 3c-methods in geometry optimization with RI-SCS-MP2 as the reference method. Similarly, the performance for accurate and efficient energy evaluation was compared with DLPNO-CCSD(T) as the reference method. The lowest energy structure among the isomers of the largest stoichiometry from our data set, AuSR18, i.e., Au3(SCH3)3, is considered to evaluate the computational time for SCF and gradient evaluations. Alongside this, the numbers of optimization steps to locate the most stable minima of Au3(SCH3)3 are compared to assess the efficiency of the methods. A comparison of relevant bond lengths with the reference geometries was made to estimate the accuracy in geometry optimization. Some methods, such as LC-BLYP, ωB97M-D3BJ, M06-2X, and PBEh-3c, could not locate many of the minima found by most of the other methods; thus, the versatility in locating various minima is also an important criterion in choosing a method for the given project. To determine the accuracy of the methods, we compared the relative energies of the isomers in each stoichiometry and the interaction energy of the gold core with the ligands. The dependence of basis set size and relativistic effects on energies are also compared. The following are some of the highlights. TPSS has shown accuracy, while mPWPW shows comparable speed and accuracy. For the relative energies of the clusters, the hybrid range-separated DFAs are the best option. CAM-B3LYP excels, whereas B3LYP performs poorly. Overall, LC-BLYP is a balanced performer considering both the geometry and relative stability of the structures, but it lacks diversity. The 3c-methods, although fast, are less impressive in relative stability.
Collapse
Affiliation(s)
- Maya Khatun
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Sayan Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Saikat Roy
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Subhasis Dey
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Anakuthil Anoop
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
5
|
Mehta N, Martin JML. Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12. J Phys Chem Lett 2022; 13:9332-9338. [PMID: 36178852 PMCID: PMC9575149 DOI: 10.1021/acs.jpclett.2c02620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 09/26/2022] [Indexed: 06/16/2023]
Abstract
Following earlier work [Mehta, N.; Martin, J. M. L. J. Chem. Theory Comput.2022, 10.1021/acs.jctc.2c00426] that showed how the slow basis set convergence of the double hybrid density functional theory can be obviated by the use of F12 explicit correlation in the GLPT2 step (second order Görling-Levy perturbation theory), we demonstrate here for the very large and chemically diverse GMTKN55 benchmark suite that the CPU time scaling of this step can be reduced (asymptotically linearized) using the localized pair natural orbital (PNO-L) approximation at negligible cost in accuracy.
Collapse
|
6
|
Zaboli M, Shariatinia Z. Highly efficient hole transport derivatives based on fluoranthene core for application in perovskite solar cells. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
7
|
Mehta N, Martin JML. Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database. J Chem Theory Comput 2022; 18:5978-5991. [PMID: 36099641 PMCID: PMC9558368 DOI: 10.1021/acs.jctc.2c00426] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Double-hybrid density functional theory (DHDFT) offers
a pathway
to accuracy approaching composite wavefunction approaches such as
G4 theory. However, the Görling–Levy second-order perturbation
theory (GLPT2) term causes them to partially inherit the slow ∝L–3 (with L the maximum
angular momentum) basis set convergence of correlated wavefunction
methods. This could potentially be remedied by introducing F12 explicit
correlation: we investigate the basis set convergence of both DHDFT
and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large
and chemically diverse general main-group thermochemistry, kinetics,
and noncovalent interactions (GMTKN55) benchmark suite. The B2GP-PLYP-D3(BJ)
and revDSD-PBEP86-D4 DHDFs are investigated as test cases, together
with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets
as large as cc-pVQZ-F12. We show that F12 greatly accelerates basis
set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12
basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPPD
in orbital-based approaches, and in fact comparable in quality to
cc-pV(5+d)Z. Somewhat surprisingly, aug-cc-pVDZ-F12 is not required
even for the anionic subsets. In conclusion, DHDF-F12/VDZ-F12 eliminates
concerns about basis set convergence in both the development and applications
of double-hybrid functionals. Mass storage and I/O bottlenecks for
larger systems can be circumvented by localized pair natural orbital
approximations, which also exhibit much gentler system size scaling.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
8
|
Informing geometric deep learning with electronic interactions to accelerate quantum chemistry. Proc Natl Acad Sci U S A 2022; 119:e2205221119. [PMID: 35901215 PMCID: PMC9351474 DOI: 10.1073/pnas.2205221119] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Predicting electronic energies, densities, and related chemical properties can facilitate the discovery of novel catalysts, medicines, and battery materials. However, existing machine learning techniques are challenged by the scarcity of training data when exploring unknown chemical spaces. We overcome this barrier by systematically incorporating knowledge of molecular electronic structure into deep learning. By developing a physics-inspired equivariant neural network, we introduce a method to learn molecular representations based on the electronic interactions among atomic orbitals. Our method, OrbNet-Equi, leverages efficient tight-binding simulations and learned mappings to recover high-fidelity physical quantities. OrbNet-Equi accurately models a wide spectrum of target properties while being several orders of magnitude faster than density functional theory. Despite only using training samples collected from readily available small-molecule libraries, OrbNet-Equi outperforms traditional semiempirical and machine learning-based methods on comprehensive downstream benchmarks that encompass diverse main-group chemical processes. Our method also describes interactions in challenging charge-transfer complexes and open-shell systems. We anticipate that the strategy presented here will help to expand opportunities for studies in chemistry and materials science, where the acquisition of experimental or reference training data is costly.
Collapse
|
9
|
Semidalas E, Martin JML. Automatic generation of complementary auxiliary basis sets for explicitly correlated methods. J Comput Chem 2022; 43:1690-1700. [PMID: 35852227 PMCID: PMC9544771 DOI: 10.1002/jcc.26970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Revised: 06/13/2022] [Accepted: 07/02/2022] [Indexed: 11/15/2022]
Abstract
Explicitly correlated calculations, aside from the orbital basis set, typically require three auxiliary basis sets: Coulomb‐exchange fitting (JK), resolution of the identity MP2 (RI‐MP2), and complementary auxiliary basis set (CABS). If unavailable for the orbital basis set and chemical elements of interest, the first two can be auto‐generated on the fly using existing algorithms, but not the third. In this paper, we present a quite simple algorithm named autoCABS; a Python implementation under a free software license is offered at Github. For the cc‐pVnZ‐F12 (n = D,T,Q,5), the W4‐08 thermochemical benchmark, and the HFREQ2014 set of harmonic frequencies, we demonstrate that autoCABS‐generated CABS basis sets are comparable in quality to purpose‐optimized OptRI basis sets from the literature, and that the quality difference becomes entirely negligible as n increases.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science Weizmann Institute of Science Reḥovot Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry and Materials Science Weizmann Institute of Science Reḥovot Israel
| |
Collapse
|
10
|
Ehlert S, Grimme S, Hansen A. Conformational Energy Benchmark for Longer n-Alkane Chains. J Phys Chem A 2022; 126:3521-3535. [PMID: 35616628 DOI: 10.1021/acs.jpca.2c02439] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the first benchmark set focusing on the conformational energies of highly flexible, long n-alkane chains, termed ACONFL. Unbranched alkanes are ubiquitous building blocks in nature, so the goal is to be able to calculate their properties most accurately to improve the modeling of, e.g., complex (biological) systems. Very accurate DLPNO-CCSD(T1)/CBS reference values are provided, which allow for a statistical meaningful evaluation of even the best available density functional methods. The performance of established and modern (dispersion corrected) density functionals is comprehensively assessed. The recently introduced r2SCAN-V functional shows excellent performance, similar to efficient composite DFT methods like B97-3c and r2SCAN-3c, which provide an even better cost-accuracy ratio, while almost reaching the accuracy of much more computationally demanding hybrid or double hybrid functionals with large QZ AO basis sets. In addition, we investigated the performance of common wave function methods, where MP2/CBS surprisingly performs worse compared to the simple D4 dispersion corrected Hartree-Fock. Furthermore, we investigate the performance of several semiempirical and force field methods, which are commonly used for the generation of conformational ensembles in multilevel workflows or in large scale molecular dynamics studies. Outstanding performance is obtained by the recently introduced general force field, GFN-FF, while other commonly applied methods like the universal force field yield large errors. We recommend the ACONFL as a helpful benchmark set for parametrization of new semiempirical or force field methods and machine learning potentials as well as a meaningful validation set for newly developed DFT or dispersion methods.
Collapse
Affiliation(s)
- Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany
| |
Collapse
|
11
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
12
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2208-2232. [PMID: 35313106 DOI: 10.1021/acs.jctc.1c01128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- MALTA Consolider Team, Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, E-33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
13
|
Beran GJO, Greenwell C, Rezac J. Spin-component-scaled and dispersion-corrected second-order Møller-Plesset perturbation theory: A path toward chemical accuracy. Phys Chem Chem Phys 2022; 24:3695-3712. [DOI: 10.1039/d1cp04922d] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Second-order Moller-Plesset perturbation theory (MP2) provides a valuable alternative to density functional theory for modeing problems in organic and biological chemistry. However, MP2 suffers from known limitations in the description...
Collapse
|
14
|
Rodrigues AC, de M. Camargo LT, Francisco Lopes Y, Sallum LO, Napolitano HB, Camargo AJ. Aqueous solvation study of melatonin using ab initio molecular dynamics. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117451] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
15
|
Ferrero R, Pantaleone S, Delle Piane M, Caldera F, Corno M, Trotta F, Brunella V. On the Interactions of Melatonin/β-Cyclodextrin Inclusion Complex: A Novel Approach Combining Efficient Semiempirical Extended Tight-Binding (xTB) Results with Ab Initio Methods. Molecules 2021; 26:5881. [PMID: 34641423 PMCID: PMC8512077 DOI: 10.3390/molecules26195881] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 09/10/2021] [Accepted: 09/21/2021] [Indexed: 11/17/2022] Open
Abstract
Melatonin (MT) is a molecule of paramount importance in all living organisms, due to its presence in many biological activities, such as circadian (sleep-wake cycle) and seasonal rhythms (reproduction, fattening, molting, etc.). Unfortunately, it suffers from poor solubility and, to be used as a drug, an appropriate transport vehicle has to be developed, in order to optimize its release in the human tissues. As a possible drug-delivery system, β-cyclodextrin (βCD) represents a promising scaffold which can encapsulate the melatonin, releasing when needed. In this work, we present a computational study supported by experimental IR spectra on inclusion MT/βCD complexes. The aim is to provide a robust, accurate and, at the same time, low-cost methodology to investigate these inclusion complexes both with static and dynamic simulations, in order to study the main actors that drive the interactions of melatonin with β-cyclodextrin and, therefore, to understand its release mechanism.
Collapse
Affiliation(s)
- Riccardo Ferrero
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, via P. Giuria 7, IT-10125 Torino, Italy; (R.F.); (S.P.); (F.C.); (F.T.)
| | - Stefano Pantaleone
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, via P. Giuria 7, IT-10125 Torino, Italy; (R.F.); (S.P.); (F.C.); (F.T.)
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, via Elce di Sotto 8, I-06123 Perugia, Italy
| | - Massimo Delle Piane
- Politecnico di Torino, Department of Applied Science and Technology (DISAT), Corso Duca degli Abruzzi, 24, 10129 Torino, Italy;
| | - Fabrizio Caldera
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, via P. Giuria 7, IT-10125 Torino, Italy; (R.F.); (S.P.); (F.C.); (F.T.)
| | - Marta Corno
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, via P. Giuria 7, IT-10125 Torino, Italy; (R.F.); (S.P.); (F.C.); (F.T.)
| | - Francesco Trotta
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, via P. Giuria 7, IT-10125 Torino, Italy; (R.F.); (S.P.); (F.C.); (F.T.)
| | - Valentina Brunella
- Dipartimento di Chimica and Nanostructured Interfaces and Surfaces (NIS) Centre, Università degli Studi di Torino, via P. Giuria 7, IT-10125 Torino, Italy; (R.F.); (S.P.); (F.C.); (F.T.)
| |
Collapse
|
16
|
Mehta N, Fellowes T, White JM, Goerigk L. CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions? J Chem Theory Comput 2021; 17:2783-2806. [PMID: 33881869 DOI: 10.1021/acs.jctc.1c00006] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| | - Thomas Fellowes
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia.,Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Victoria 3010, Australia
| | - Jonathan M White
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia.,Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
17
|
Gutten O, Jurečka P, Aliakbar Tehrani Z, Buděšínský M, Řezáč J, Rulíšek L. Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments. Phys Chem Chem Phys 2021; 23:7280-7294. [PMID: 33876088 DOI: 10.1039/d0cp05993e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.
Collapse
Affiliation(s)
- Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | | | | | | | | |
Collapse
|
18
|
Mehta N, Goerigk L. Assessing the Applicability of the Geometric Counterpoise Correction in B2PLYP/Double-ζ Calculations for Thermochemistry, Kinetics, and Noncovalent Interactions. Aust J Chem 2021. [DOI: 10.1071/ch21133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
We present a proof-of-concept study of the suitability of Kruse and Grimme’s geometric counterpoise correction (gCP) for basis set superposition errors (BSSEs) in double-hybrid density functional calculations with a double-ζ basis set. The gCP approach only requires geometrical information as an input and no orbital/density information is needed. Therefore, this correction is practically free of any additional cost. gCP is trained against the Boys and Bernardi counterpoise correction across a set of 528 noncovalently bound dimers. We investigate the suitability of the approach for the B2PLYP/def2-SVP level of theory, and reveal error compensation effects—missing London dispersion and the BSSE—associated with B2PLYP/def2-SVP calculations, and present B2PLYP-gCP-D3(BJ)/def2-SVP with the reparametrised DFT-D3(BJ) and gCP corrections as a more balanced alternative. Benchmarking results on the S66x8 benchmark set for noncovalent interactions and the GMTKN55 database for main-group thermochemistry, kinetics, and noncovalent interactions show a statistical improvement of the B2PLYP-gCP-D3(BJ) scheme over plain B2PLYP and B2PLYP-D3(BJ). B2PLYP-D3(BJ) shows significant overestimation of interaction energies, barrier heights with larger deviations from the reference values, and wrong relative stabilities in conformers, all of which can be associated with BSSE. We find that the gCP-corrected method represents a significant improvement over B2PLYP-D3(BJ), particularly for intramolecular noncovalent interactions. These findings encourage future developments of efficient double-hybrid DFT strategies that can be applied when double-hybrid calculations with large basis sets are not feasible due to system size.
Collapse
|
19
|
Williams ZM, Wiles TC, Manby FR. Accurate Hybrid Density Functionals with UW12 Correlation. J Chem Theory Comput 2020; 16:6176-6194. [PMID: 32820918 DOI: 10.1021/acs.jctc.0c00442] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In previous work, we suggested a single-parameter hybrid functional containing a novel correlation contribution based on the Unsöld approximation, UW12. This model resembles the explicitly correlated part of MP2-F12 theory and can be written as an explicit formula in terms of the single-particle reduced density matrix. Here, we further investigate hybrid functionals containing UW12 correlation and in particular look at functionals with a large fraction of exact exchange to reduce the self-interaction error. We suggest two new hybrid functionals B-LYP-osUW12 and fB-LYP-osUW12. On the test sets we use, our best hybrid functional overall (B-LYP-osUW12) is of similar accuracy to the best double hybrids considered while eliminating the need for virtual orbitals.
Collapse
Affiliation(s)
- Zack M Williams
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Timothy C Wiles
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| |
Collapse
|
20
|
Patra A, Jana S, Constantin LA, Samal P. Efficient yet accurate dispersion-corrected semilocal exchange–correlation functionals for non-covalent interactions. J Chem Phys 2020; 153:084117. [DOI: 10.1063/5.0011849] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Lucian A. Constantin
- Istituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, Italy
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
21
|
Förster A, Visscher L. Double hybrid DFT calculations with Slater type orbitals. J Comput Chem 2020; 41:1660-1684. [PMID: 32297682 PMCID: PMC7317772 DOI: 10.1002/jcc.26209] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Revised: 03/31/2020] [Accepted: 04/01/2020] [Indexed: 12/14/2022]
Abstract
On a comprehensive database with 1,644 datapoints, covering several aspects of main-group as well as of transition metal chemistry, we assess the performance of 60 density functional approximations (DFA), among them 36 double hybrids (DH). All calculations are performed using a Slater type orbital (STO) basis set of triple-ζ (TZ) quality and the highly efficient pair atomic resolution of the identity approach for the exchange- and Coulomb-term of the KS matrix (PARI-K and PARI-J, respectively) and for the evaluation of the MP2 energy correction (PARI-MP2). Employing the quadratic scaling SOS-AO-PARI-MP2 algorithm, DHs based on the spin-opposite-scaled (SOS) MP2 approximation are benchmarked against a database of large molecules. We evaluate the accuracy of STO/PARI calculations for B3LYP as well as for the DH B2GP-PLYP and show that the combined basis set and PARI-error is comparable to the one obtained using the well-known def2-TZVPP Gaussian-type basis set in conjunction with global density fitting. While quadruple-ζ (QZ) calculations are currently not feasible for PARI-MP2 due to numerical issues, we show that, on the TZ level, Jacob's ladder for classifying DFAs is reproduced. However, while the best DHs are more accurate than the best hybrids, the improvements are less pronounced than the ones commonly found on the QZ level. For conformers of organic molecules and noncovalent interactions where very high accuracy is required for qualitatively correct results, DHs provide only small improvements over hybrids, while they still excel in thermochemistry, kinetics, transition metal chemistry and the description of strained organic systems.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical ChemistryVrije UniversiteitAmsterdamThe Netherlands
| | - Lucas Visscher
- Theoretical ChemistryVrije UniversiteitAmsterdamThe Netherlands
| |
Collapse
|
22
|
Zhang D, Truhlar DG. Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory. J Chem Theory Comput 2020; 16:4416-4428. [PMID: 32525690 DOI: 10.1021/acs.jctc.0c00518] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe2+, Fe3+, and Co3+, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF2, FeCl2, and CoF2. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.
Collapse
Affiliation(s)
- Dayou Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
23
|
Chan B. Aqueous-Phase Conformations of Lactose, Maltose, and Sucrose and the Assessment of Low-Cost DFT Methods with the DSCONF Set of Conformers for the Three Disaccharides. J Phys Chem A 2020; 124:582-590. [PMID: 31927999 DOI: 10.1021/acs.jpca.9b10932] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
In the present study, we have examined a range of quantum chemistry methods for the calculation of conformers for lactose, maltose, and sucrose. We find that the DSD-PBE-P86/aug'-cc-pVTZ//B3-LYP-D3BJ/6-311+G(2d,p) protocol yields good relative energies in comparison with reference CCSD(T)/CBS//B3-LYP-D3BJ/maug-cc-pVTZ values. We have surveyed a total of ∼550 conformers for the three disaccharides with the chosen DSD-PBE-P86 method in conjunction with continuum aqueous solvation. In each case, the lowest free energy conformer is characterized by hydrogen bond(s) between the two rings. Another finding is that the major contributors to the overall variations in aqueous free energies are the electronic energies and the solvation energies. To facilitate investigations of larger systems, we have compiled the DSCONF set of conformers for the three disaccharides, and we have assessed lower cost methods with this set. We find MS1-D3/6-31+G(2d,p) to be cost-effective and accurate for both geometry optimization and the calculation of relative energies for disaccharides. In addition, we note that MS1-D3 has previously been found to yield good relative energies for the WATER27 set of water clusters. We thus deem this method to be appropriate for the study of saccharide conformations in both gas phase and aqueous solution.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering , Nagasaki University , Bunkyo 1-14 , Nagasaki 852-8521 , Japan
| |
Collapse
|
24
|
Matczak P, Mlostoń G, Hamera-Fałdyga R, Görls H, Weigand W. Structure of Diferrocenyl Thioketone: From Molecule to Crystal. MOLECULES (BASEL, SWITZERLAND) 2019; 24:molecules24213950. [PMID: 31683693 PMCID: PMC6864675 DOI: 10.3390/molecules24213950] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 10/28/2019] [Accepted: 10/29/2019] [Indexed: 01/07/2023]
Abstract
Ferrocenyl-functionalized thioketones have recently been recognized as useful building blocks for sulfur-containing compounds with potential applications in materials chemistry. This work is devoted to a single representative of such thioketones, namely diferrocenyl thioketone (Fc2CS), whose structure has been determined here for the first time. Both X-ray crystallography and a wide variety of quantum-chemical methods were used to explore the structure of Fc2CS. In addition to the X-ray structure determination, intermolecular interactions occurring in the crystal structure of Fc2CS were examined in detail by quantum-chemical methods. These methods were also an invaluable tool in studying the molecular structure of Fc2CS, from the gas phase to solutions and to its crystal. Intramolecular interactions governing the conformational behavior of an isolated Fc2CS molecule were deduced from quantum-chemical analyses carried out in orbital space and real space. Our experimental and theoretical results indicate that the main structural features of an isolated Fc2CS molecule in its lowest-energy geometry are retained both upon solvation and in the crystal. The tilt of ferrocenyl groups is only slightly affected by crystal packing forces that are dominated by dispersion. Nonetheless, a network of intermolecular interactions, such as H···H, C···H and S···H, was detected in the Fc2CS crystal but each of them is fairly weak.
Collapse
Affiliation(s)
- Piotr Matczak
- Department of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, 90236 Lodz, Poland.
| | - Grzegorz Mlostoń
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Róża Hamera-Fałdyga
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Helmar Görls
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| | - Wolfgang Weigand
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| |
Collapse
|
25
|
Hehre W, Klunzinger P, Deppmeier B, Driessen A, Uchida N, Hashimoto M, Fukushi E, Takata Y. Efficient Protocol for Accurately Calculating 13C Chemical Shifts of Conformationally Flexible Natural Products: Scope, Assessment, and Limitations. JOURNAL OF NATURAL PRODUCTS 2019; 82:2299-2306. [PMID: 31322872 DOI: 10.1021/acs.jnatprod.9b00603] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
An efficient protocol for calculating 13C NMR chemical shifts for natural products with multiple degrees of conformational freedom is described. This involves a multistep procedure starting from molecular mechanics and ending with a large basis set density functional model to obtain accurate Boltzmann conformer weights, followed by empirically corrected density functional NMR calculations for the individual conformers. The accuracy of the protocol (average rms <4 ppm) was determined by application to ∼925 diverse natural products, the structures of which have been confirmed either by X-ray crystallography or independent synthesis. The protocol was then applied to ∼ 2275 natural products, the structures of which were elucidated mainly by NMR and MS data. Five to ten percent of the latter compounds exhibited rms errors significantly in excess of 4 ppm, suggesting possible structural or signal assignment errors. Both data sets are available from an online browser ( NMR.wavefun.com ). The procedure can be and has been fully automated and is practical using present-generation personal computers, requiring a few hours or days depending on the size of the molecule and number of accessible conformers.
Collapse
Affiliation(s)
- Warren Hehre
- Wavefunction, Inc. , 18401 Von Karman Avenue, Suite 370 , Irvine , California 92612-1542 , United States
| | - Phillip Klunzinger
- Wavefunction, Inc. , 18401 Von Karman Avenue, Suite 370 , Irvine , California 92612-1542 , United States
| | - Bernard Deppmeier
- Wavefunction, Inc. , 18401 Von Karman Avenue, Suite 370 , Irvine , California 92612-1542 , United States
| | - Andy Driessen
- Wavefunction, Inc. , 18401 Von Karman Avenue, Suite 370 , Irvine , California 92612-1542 , United States
| | - Noritaka Uchida
- Wavefunction Japan , 3-5-2, Koji-Machi , Chiyoda-Ku , Tokyo 102-0083 , Japan
| | - Masaru Hashimoto
- Faculty of Agriculture and Life Science , Hirosaki University , 3-Bunkyo-cho , Hirosaki 036-8561 , Japan
| | - Eri Fukushi
- Graduate School of Agriculture , Hokkaido University , Sapporo 060-8589 , Japan
| | - Yusuke Takata
- Graduate School of Agriculture , Hokkaido University , Sapporo 060-8589 , Japan
| |
Collapse
|
26
|
Wappett DA, Goerigk L. Toward a Quantum-Chemical Benchmark Set for Enzymatically Catalyzed Reactions: Important Steps and Insights. J Phys Chem A 2019; 123:7057-7074. [DOI: 10.1021/acs.jpca.9b05088] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
27
|
Santra G, Sylvetsky N, Martin JML. Minimally Empirical Double-Hybrid Functionals Trained against the GMTKN55 Database: revDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4. J Phys Chem A 2019; 123:5129-5143. [PMID: 31136709 PMCID: PMC9479158 DOI: 10.1021/acs.jpca.9b03157] [Citation(s) in RCA: 223] [Impact Index Per Article: 44.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
![]()
We
present a family of minimally empirical double-hybrid DFT functionals
parametrized against the very large and diverse GMTKN55 benchmark.
The very recently proposed ωB97M(2) empirical double hybrid
(with 16 adjustable parameters) has the lowest WTMAD2 (weighted mean
absolute deviation over GMTKN55) ever reported at 2.19 kcal/mol. However,
refits of the DSD-BLYP and DSD-PBEP86 spin-component-scaled, dispersion-corrected
double hybrids can achieve WTMAD2 values as low as 2.33 with the very
recent D4 dispersion correction (2.42 kcal/mol with the D3(BJ) dispersion
term) using just a handful of adjustable parameters. If we use full
DFT correlation in the initial orbital evaluation, the xrevDSD-PBEP86-D4
functional reaches WTMAD2 = 2.23 kcal/mol, statistically indistinguishable
from ωB97M(2) but using just four nonarbitrary adjustable parameters
(and three semiarbitrary ones). The changes from the original DSD
parametrizations are primarily due to noncovalent interaction energies
for large systems, which were undersampled in the original parametrization
set. With the new parametrization, same-spin correlation can be eliminated
at minimal cost in performance, which permits revDOD-PBEP86-D4 and
revDOD-PBE-D4 functionals that scale as N4 or even N3 with the size of the system.
Dependence of WTMAD2 for DSD functionals on the percentage of HF exchange
is roughly quadratic; it is sufficiently weak that any reasonable
value in the 64% to 72% range can be chosen semiarbitrarily. DSD-SCAN
and DOD-SCAN double hybrids involving the SCAN nonempirical meta-GGA
as the semilocal component have also been considered and offer a good
alternative if one wishes to eliminate either the empirical dispersion
correction or the same-spin correlation component. noDispSD-SCAN66
achieves WTMAD2 = 3.0 kcal/mol, compared to 2.7 kcal/mol for DOD-SCAN66-D4.
However, the best performance without dispersion corrections (WTMAD2
= 2.8 kcal/mol) is reached by revωB97X-2, a slight reparametrization
of the Chai–Head-Gordon range-separated double hybrid. Finally,
in the context of double-hybrid functionals, the very recent D4 dispersion
correction is clearly superior over D3(BJ).
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
28
|
Bortoli M, Dalla Tiezza M, Muraro C, Pavan C, Ribaudo G, Rodighiero A, Tubaro C, Zagotto G, Orian L. Psychiatric Disorders and Oxidative Injury: Antioxidant Effects of Zolpidem Therapy disclosed In Silico. Comput Struct Biotechnol J 2019; 17:311-318. [PMID: 30867894 PMCID: PMC6396081 DOI: 10.1016/j.csbj.2019.02.004] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2018] [Revised: 02/01/2019] [Accepted: 02/02/2019] [Indexed: 01/09/2023] Open
Abstract
Zolpidem (N,N-Dimethyl-2-[6-methyl-2-(4-methylphenyl)imidazo[1,2-a]pyridin-3-yl]acetamide) is a well-known drug for the treatment of sleeping disorders. Recent literature reports on positive effects of zolpidem therapy on improving renal damage after cisplatin and on reducing akinesia without sleep induction. This has been ascribed to the antioxidant and neuroprotective capacity of this molecule, and tentatively explained according to a generic structural similarity between zolpidem and melatonin. In this work, we investigate in silico the antioxidant potential of zolpidem as scavenger of five ROSs, acting via hydrogen atom transfer (HAT) mechanism; computational methodologies based on density functional theory are employed. For completeness, the analysis is extended to six metabolites. Thermodynamic and kinetic results disclose that indeed zolpidem is an efficient radical scavenger, similarly to melatonin and Trolox, supporting the biomedical evidence that the antioxidant potential of zolpidem therapy may have a beneficial effect against oxidative injury, which is emerging as an important etiopathogenesis in numerous severe diseases, including psychiatric disorders.
Collapse
Key Words
- Antioxidant activity
- DFT calculations
- DFT, Density Functional Theory
- HAT, Hydrogen Atom Transfer (mechanism)
- M06-2X, Minnesota Hybrid functional with 54% Hartree-Fock exchange
- NBO, Natural Bond Orbitals
- NPA, Natural Population Analysis
- Oxidative stress
- PC, Product Complex
- Psychiatric disorders
- RAF, Radical Adduct Formation (mechanism)
- RC, Reactant Complex
- ROS, Reactive Oxygen Species
- Radical scavenging
- SMD, Solvation Model based on Density
- TS, Transition State
- Zolpidem
Collapse
Affiliation(s)
- Marco Bortoli
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| | - Marco Dalla Tiezza
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| | - Cecilia Muraro
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| | - Chiara Pavan
- Dipartimento di Medicina, Università degli Studi di Padova, Via Giustiniani 2, 35128 Padova, Italy
| | - Giovanni Ribaudo
- Dipartimento di Scienze del Farmaco, Università degli Studi di Padova, Via Marzolo 5, 35131 Padova, Italy
| | - Anna Rodighiero
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| | - Cristina Tubaro
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| | - Giuseppe Zagotto
- Dipartimento di Scienze del Farmaco, Università degli Studi di Padova, Via Marzolo 5, 35131 Padova, Italy
| | - Laura Orian
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, Via Marzolo 1, 35131 Padova, Italy
| |
Collapse
|
29
|
Cova TF, Milne BF, Pais AA. Host flexibility and space filling in supramolecular complexation of cyclodextrins: A free-energy-oriented approach. Carbohydr Polym 2019; 205:42-54. [DOI: 10.1016/j.carbpol.2018.10.009] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2018] [Revised: 10/03/2018] [Accepted: 10/03/2018] [Indexed: 12/20/2022]
|
30
|
PEPCONF, a diverse data set of peptide conformational energies. Sci Data 2019; 6:180310. [PMID: 30667382 PMCID: PMC6343515 DOI: 10.1038/sdata.2018.310] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Accepted: 11/30/2018] [Indexed: 12/13/2022] Open
Abstract
We present an extensive and diverse database of peptide conformational energies. Our database contains five different classes of model geometries: dipeptides, tripeptides, and disulfide-bridged, bioactive, and cyclic peptides. In total, the database consists of 3775 conformational energy data points and 4530 conformer geometries. All the reference energies have been calculated at the LC-ωPBE-XDM/aug-cc-pVTZ level of theory, which is shown to yield conformational energies with an accuracy in the order of tenths of a kcal/mol when compared to complete-basis-set coupled-cluster reference data. The peptide conformational data set (PEPCONF) is presented as a high-quality reference set for the development and benchmarking of molecular-mechanics and semi-empirical electronic structure methods, which are the most commonly used techniques in the modeling of medium to large proteins.
Collapse
|
31
|
Jin Y, Bartlett RJ. Perturbation Improved Natural Linear-Scaled Coupled-Cluster Method and Its Application to Conformational Analysis. J Phys Chem A 2019; 123:371-381. [PMID: 30585491 DOI: 10.1021/acs.jpca.8b07947] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The fragment-based coupled-cluster (CC) theory utilizing the transferable functional groups through natural localized molecular orbital (NLMO), that is, the natural linear-scaled coupled-cluster (NLSCC) has been further developed to take the extra-fragment interactions into account. The correction to the interaction energies sacrificed during the fragmentation process for the previous NLSCC method is computed by a computationally efficient perturbation theory that maintains the original scaling. The new linear-scaled coupled-cluster for the singles and doubles (CCSD) method is applied to the analysis of relative energies of delicate conformational problems of polypeptides. By adding a perturbation correction, results accurate to less than a kcal/mol are obtained for the alanine tetramer.
Collapse
Affiliation(s)
- Yifan Jin
- Quantum Theory Project and Departments of Chemistry and Physics , University of Florida , Gainesville , Florida 32611 , United States
| | - Rodney J Bartlett
- Quantum Theory Project and Departments of Chemistry and Physics , University of Florida , Gainesville , Florida 32611 , United States
| |
Collapse
|
32
|
Morgante P, Peverati R. Statistically representative databases for density functional theory via data science. Phys Chem Chem Phys 2019; 21:19092-19103. [DOI: 10.1039/c9cp03211h] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cluster analysis applied to quantum chemistry: a new broad database of chemical properties with a reasonable computational cost.
Collapse
|
33
|
Morgante P, Peverati R. ACCDB: A collection of chemistry databases for broad computational purposes. J Comput Chem 2018; 40:839-848. [DOI: 10.1002/jcc.25761] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 10/09/2018] [Accepted: 11/11/2018] [Indexed: 01/04/2023]
Affiliation(s)
- Pierpaolo Morgante
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| | - Roberto Peverati
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| |
Collapse
|
34
|
Sharapa DI, Genaev A, Cavallo L, Minenkov Y. A Robust and Cost‐Efficient Scheme for Accurate Conformational Energies of Organic Molecules. Chemphyschem 2018; 20:92-102. [DOI: 10.1002/cphc.201801063] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2018] [Indexed: 11/09/2022]
Affiliation(s)
- Dmitry I. Sharapa
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology (KIT) Hermann-von-Helmholtz Platz 1 Eggenstein-Leopoldshafen D-76344 Germany
| | - Alexander Genaev
- Vorozhtsov Novosibirsk Institute of Organic Chemistry Academician Lavrent'ev Ave., 9 Novosibirsk 630090 Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC)King Abdullah University of Science and Technology Thuwal- 23955-6900 Saudi Arabia
| | - Yury Minenkov
- Moscow Institute of Physics and Technology Institutskiy Pereulok 9, Dolgoprudny Moscow Region 141700 Russia
| |
Collapse
|
35
|
Najibi A, Goerigk L. The Nonlocal Kernel in van der Waals Density Functionals as an Additive Correction: An Extensive Analysis with Special Emphasis on the B97M-V and ωB97M-V Approaches. J Chem Theory Comput 2018; 14:5725-5738. [DOI: 10.1021/acs.jctc.8b00842] [Citation(s) in RCA: 108] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Asim Najibi
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
36
|
Mehta N, Casanova-Páez M, Goerigk L. Semi-empirical or non-empirical double-hybrid density functionals: which are more robust? Phys Chem Chem Phys 2018; 20:23175-23194. [PMID: 30062343 DOI: 10.1039/c8cp03852j] [Citation(s) in RCA: 89] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Parkville, Australia.
| | | | | |
Collapse
|
37
|
Xu J, Yi L, Mou Y, Cao J, Wang C. Effect of a molecule of imidazolium bromide ionic liquid on the structure and properties of cytosine by density functional theory. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.08.009] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
38
|
Nagy PR, Samu G, Kállay M. Optimization of the Linear-Scaling Local Natural Orbital CCSD(T) Method: Improved Algorithm and Benchmark Applications. J Chem Theory Comput 2018; 14:4193-4215. [DOI: 10.1021/acs.jctc.8b00442] [Citation(s) in RCA: 66] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Péter R. Nagy
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Gyula Samu
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
39
|
Mardirossian N, Head-Gordon M. Survival of the most transferable at the top of Jacob's ladder: Defining and testing the ωB97M(2) double hybrid density functional. J Chem Phys 2018; 148:241736. [PMID: 29960332 PMCID: PMC5991970 DOI: 10.1063/1.5025226] [Citation(s) in RCA: 117] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 05/09/2018] [Indexed: 01/19/2023] Open
Abstract
A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resulting functional, ωB97M(2), is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 data points clearly surpasses that of all of the tested density functionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
40
|
Řezáč J, Bím D, Gutten O, Rulíšek L. Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set. J Chem Theory Comput 2018; 14:1254-1266. [DOI: 10.1021/acs.jctc.7b01074] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Jan Řezáč
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Daniel Bím
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
41
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Atom-Centered Potentials with Dispersion-Corrected Minimal-Basis-Set Hartree–Fock: An Efficient and Accurate Computational Approach for Large Molecular Systems. J Chem Theory Comput 2018; 14:726-738. [DOI: 10.1021/acs.jctc.7b01158] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Viki Kumar Prasad
- Department
of Chemistry, University of British Columbia, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- Department
of Chemistry, University of British Columbia, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Gino A. DiLabio
- Department
of Chemistry, University of British Columbia, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
- Faculty
of Management, University of British Columbia, 1137 Alumni Avenue, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
42
|
Gould T. ‘Diet GMTKN55’ offers accelerated benchmarking through a representative subset approach. Phys Chem Chem Phys 2018; 20:27735-27739. [DOI: 10.1039/c8cp05554h] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The GMTKN55 benchmarking protocol allows comprehensive analysis and ranking of density functional approximations with diverse chemical behaviours. This work reports diet versions of GMTKN55 which reproduce key properties of the full protocol at substantially reduced numerical cost. ‘Diet GMTKN55’ can thus be used for benchmarking expensive methods, or in combination with solid state benchmarks.
Collapse
Affiliation(s)
- Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University
- Nathan
- Australia
| |
Collapse
|
43
|
Vamhindi BSDR, Karton A. Can DFT and ab initio methods adequately describe binding energies in strongly interacting C6X6⋯C2X π–π complexes? Chem Phys 2017. [DOI: 10.1016/j.chemphys.2017.05.020] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
44
|
Witte J, Neaton JB, Head-Gordon M. Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C. J Chem Phys 2017. [DOI: 10.1063/1.4986962] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
|
45
|
Mardirossian N, Head-Gordon M. Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1333644] [Citation(s) in RCA: 709] [Impact Index Per Article: 101.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
46
|
Chehelamirani M, Salahub DR. Effect of dispersion corrections on covalent and non-covalent interactions in DFTB calculations. Struct Chem 2017. [DOI: 10.1007/s11224-017-0976-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
47
|
Yoshizawa K, Semoto T, Hitaoka S, Higuchi C, Shiota Y, Tanaka H. Synergy of Electrostatic and van der Waals Interactions in the Adhesion of Epoxy Resin with Carbon-Fiber and Glass Surfaces. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2017. [DOI: 10.1246/bcsj.20160426] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Kazunari Yoshizawa
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Fukuoka 819-0395
| | - Takayuki Semoto
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Fukuoka 819-0395
| | - Seiji Hitaoka
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Fukuoka 819-0395
| | - Chisa Higuchi
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Fukuoka 819-0395
| | - Yoshihito Shiota
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Fukuoka 819-0395
| | - Hiromasa Tanaka
- Institute for Materials Chemistry and Engineering and IRCCS, Kyushu University, Fukuoka 819-0395
| |
Collapse
|
48
|
Witte J, Mardirossian N, Neaton JB, Head-Gordon M. Assessing DFT-D3 Damping Functions Across Widely Used Density Functionals: Can We Do Better? J Chem Theory Comput 2017; 13:2043-2052. [DOI: 10.1021/acs.jctc.7b00176] [Citation(s) in RCA: 56] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | | | - Jeffrey B. Neaton
- Kavli Energy
Nanosciences
Institute at Berkeley, Berkeley, California 94720, United States
| | | |
Collapse
|
49
|
Womack JC, Mardirossian N, Head-Gordon M, Skylaris CK. Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package. J Chem Phys 2017; 145:204114. [PMID: 27908114 DOI: 10.1063/1.4967960] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Finally, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.
Collapse
Affiliation(s)
- James C Womack
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
50
|
Bew SP, Stephenson GR, Rouden J, Godemert J, Seylani H, Martinez-Lozano LA. Gaining Insight Into Reactivity Differences Between Malonic Acid Half Thioesters (MAHT) and Malonic Acid Half Oxyesters (MAHO). Chemistry 2017; 23:4557-4569. [DOI: 10.1002/chem.201605148] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2016] [Indexed: 01/15/2023]
Affiliation(s)
- Sean P. Bew
- School of Chemistry; Norwich Research Park; University of East Anglia; NR4 7TJ UK
| | | | - Jacques Rouden
- Laboratoire de Chimie Moleculaire et thio-organique (LCMT); UMR CNRS 6507, Ensicaen; 6 Boulevard du Marechal Juin 14050 Caen France
| | - Jeremy Godemert
- School of Chemistry; Norwich Research Park; University of East Anglia; NR4 7TJ UK
| | - Haseena Seylani
- School of Chemistry; Norwich Research Park; University of East Anglia; NR4 7TJ UK
| | | |
Collapse
|