1
|
Riedmiller K, Reiser P, Bobkova E, Maltsev K, Gryn'ova G, Friederich P, Gräter F. Substituting density functional theory in reaction barrier calculations for hydrogen atom transfer in proteins. Chem Sci 2024; 15:2518-2527. [PMID: 38362411 PMCID: PMC10866341 DOI: 10.1039/d3sc03922f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 01/10/2024] [Indexed: 02/17/2024] Open
Abstract
Hydrogen atom transfer (HAT) reactions are important in many biological systems. As these reactions are hard to observe experimentally, it is of high interest to shed light on them using simulations. Here, we present a machine learning model based on graph neural networks for the prediction of energy barriers of HAT reactions in proteins. As input, the model uses exclusively non-optimized structures as obtained from classical simulations. It was trained on more than 17 000 energy barriers calculated using hybrid density functional theory. We built and evaluated the model in the context of HAT in collagen, but we show that the same workflow can easily be applied to HAT reactions in other biological or synthetic polymers. We obtain for relevant reactions (small reaction distances) a model with good predictive power (R2 ∼ 0.9 and mean absolute error of <3 kcal mol-1). As the inference speed is high, this model enables evaluations of dozens of chemical situations within seconds. When combined with molecular dynamics in a kinetic Monte-Carlo scheme, the model paves the way toward reactive simulations.
Collapse
Affiliation(s)
- Kai Riedmiller
- Heidelberg Institute for Theoretical Studies Heidelberg Germany
| | - Patrick Reiser
- Institute of Theoretical Informatics, Karlsruhe Institute of Technology Engler-Bunte-Ring 8 Karlsruhe 76131 Germany
- Institute of Nanotechnology, Karlsruhe Institute of Technology Hermann-von-Helmholtz-Platz 1: 76344 Eggenstein-Leopoldshafen Germany
| | | | - Kiril Maltsev
- Heidelberg Institute for Theoretical Studies Heidelberg Germany
| | - Ganna Gryn'ova
- Heidelberg Institute for Theoretical Studies Heidelberg Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University Heidelberg Germany
| | - Pascal Friederich
- Institute of Theoretical Informatics, Karlsruhe Institute of Technology Engler-Bunte-Ring 8 Karlsruhe 76131 Germany
- Institute of Nanotechnology, Karlsruhe Institute of Technology Hermann-von-Helmholtz-Platz 1: 76344 Eggenstein-Leopoldshafen Germany
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies Heidelberg Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University Heidelberg Germany
| |
Collapse
|
2
|
Pitman SJ, Evans AK, Ireland RT, Lempriere F, McKemmish LK. Benchmarking Basis Sets for Density Functional Theory Thermochemistry Calculations: Why Unpolarized Basis Sets and the Polarized 6-311G Family Should Be Avoided. J Phys Chem A 2023; 127:10295-10306. [PMID: 37982604 DOI: 10.1021/acs.jpca.3c05573] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
Basis sets are a crucial but often largely overlooked choice in setting up quantum chemistry calculations. The choice of the basis set can be critical in determining the accuracy and calculation time of your quantum chemistry calculations. Clear recommendations based on thorough benchmarking are essential but not readily available currently. This study investigates the relative quality of basis sets for general properties by benchmarking basis set performance for a diverse set of 139 reactions (from the diet-150-GMTKN55 data set). In our analysis, we find the distributions of errors are often significantly non-Gaussian, meaning that the joint consideration of median errors, mean absolute errors, and outlier statistics is helpful to provide a holistic understanding of basis set performance. Our direct comparison of performance between most modern basis sets provides quantitative evidence for basis set recommendations that broadly align with the established understanding of basis set experts and is evident in the design of modern basis sets. For example, while zeta is a good measure of quality, it is not the only determining factor for an accurate calculation with unpolarized double- and triple-ζ basis sets (like 6-31G and 6-311G) having very poor performance. Appropriate use of polarization functions (e.g., 6-31G*) is essential to obtain the accuracy offered by double- or triple-ζ basis sets. In our study, the best performances for double- and triple-ζ basis sets are 6-31++G** and pcseg-2, respectively. However, the performances of singly polarized double-ζ and doubly polarized triple-ζ basis sets are quite similar with one key exception: the polarized 6-311G basis set family has poor parametrization, which means its performance is more like a double-ζ than a triple-ζ basis set. All versions of the 6-311G basis set family should be avoided entirely for valence chemistry calculations moving forward.
Collapse
Affiliation(s)
- Samuel J Pitman
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Alicia K Evans
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Robbie T Ireland
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Felix Lempriere
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Laura K McKemmish
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
3
|
Pederson R, Burke K. The difference between molecules and materials: Reassessing the role of exact conditions in density functional theory. J Chem Phys 2023; 159:214113. [PMID: 38054515 DOI: 10.1063/5.0172058] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Accepted: 11/12/2023] [Indexed: 12/07/2023] Open
Abstract
Exact conditions have long been used to guide the construction of density functional approximations. However, hundreds of empirical-based approximations tailored for chemistry are in use, of which many neglect these conditions in their design. We analyze well-known conditions and revive several obscure ones. Two crucial distinctions are drawn: that between necessary and sufficient conditions and that between all electronic densities and the subset of realistic Coulombic ground states. Simple search algorithms find that many empirical approximations satisfy many exact conditions for realistic densities and non-empirical approximations satisfy even more conditions than those enforced in their construction. The role of exact conditions in developing approximations is revisited.
Collapse
Affiliation(s)
- Ryan Pederson
- Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
| | - Kieron Burke
- Department of Physics and Astronomy, University of California, Irvine, California 92697, USA
- Department of Chemistry, University of California, Irvine, California 92697, USA
| |
Collapse
|
4
|
Hernández‐Castillo D, Nau REP, Schmid M, Tschierlei S, Rau S, González L. Mehrere Triplett-Metall-zentrierte Jahn-Teller-Isomere bestimmen die temperaturabhängigen Lumineszenzlebensdauern in [Ru(bpy) 3] 2. ANGEWANDTE CHEMIE (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023; 135:e202308803. [PMID: 38529088 PMCID: PMC10962581 DOI: 10.1002/ange.202308803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/22/2023] [Indexed: 03/27/2024]
Abstract
AbstractEin genaues Verständnis der Faktoren, welche die Lumineszenzlebensdauer von Übergangsmetallverbindungen bestimmen, ist für Anwendungen in der Photokatalyse und der photodynamischen Therapie von entscheidender Bedeutung. Die im Falle von [Ru(bpy)3]2+ (bpy=2,2’‐Bipyridin) allgemein akzeptierte Theorie besagt, dass die Emissionslebensdauer durch Optimierung der Energiebarriere zwischen dem emittierenden Triplett‐Zustand des Metall‐Liganden‐Ladungstransfers (3MLCT) und dem thermisch aktivierten Triplett‐Zustand des Metall‐Zentrums (3MC), oder der Energielücke zwischen beiden Zuständen gesteuert werden kann. Hier zeigen wir, dass dies nicht allgemeingültig ist. Darüber hinaus demonstrieren wir, dass die Betrachtung eines einzelnen Relaxationspfades, der vom energetisch niedrigsten Minimum aus bestimmt wird, zu falschen Vorhersagen der temperaturabhängigen Emissionslebensdauer führt. Stattdessen erhalten wir eine ausgezeichnete Übereinstimmung mit den experimentellen temperaturabhängigen Lebensdauern, wenn ein erweitertes kinetisches Modell herangezogen wird, welches alle Pfade im Zusammenhang mit mehreren Jahn–Teller‐Isomeren und ihren effektiven Reaktionsbarrieren beinhaltet. Diese Konzepte sind für das Design weiterer lumineszierender Übergangsmetallkomplexe mit individuell angepassten Emissionslebensdauern auf der Grundlage theoretischer Vorhersagen unerlässlich.
Collapse
Affiliation(s)
- David Hernández‐Castillo
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Str. 171090ViennaAustria
- Doctoral School in Chemistry (DoSChem)University of ViennaWähringer Straße 421090ViennaAustria
| | - Roland E. P. Nau
- Institute of Inorganic Chemistry IUlm UniversityAlbert-Einstein-Allee 1189081UlmGermany
| | - Marie‐Ann Schmid
- Technische Universität BraunschweigDepartment of Energy Conversion, Institute of Physical and Theoretical ChemistryRebenring 3138106BraunschweigGermany
| | - Stefanie Tschierlei
- Technische Universität BraunschweigDepartment of Energy Conversion, Institute of Physical and Theoretical ChemistryRebenring 3138106BraunschweigGermany
| | - Sven Rau
- Institute of Inorganic Chemistry IUlm UniversityAlbert-Einstein-Allee 1189081UlmGermany
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Str. 171090ViennaAustria
- Vienna Research Platform Accelerating Photoreaction DiscoveryUniversity of ViennaWähringer Straße 171090ViennaAustria
| |
Collapse
|
5
|
Hernández‐Castillo D, Nau REP, Schmid M, Tschierlei S, Rau S, González L. Multiple Triplet Metal-Centered Jahn-Teller Isomers Determine Temperature-Dependent Luminescence Lifetimes in [Ru(bpy) 3 ] 2. Angew Chem Int Ed Engl 2023; 62:e202308803. [PMID: 37433755 PMCID: PMC10962642 DOI: 10.1002/anie.202308803] [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: 06/22/2023] [Revised: 07/11/2023] [Accepted: 07/11/2023] [Indexed: 07/13/2023]
Abstract
Understanding the factors that determine the luminescence lifetime of transition metal compounds is key for applications in photocatalysis and photodynamic therapy. Here we show that for[ Ru ( bpy ) 3 ] 2 + ${[{\rm{Ru}}({\rm{bpy}})_{\rm{3}} ]^{{\rm{2 + }}} }$ (bpy = 2,2'-bipyridine), the generally accepted idea that emission lifetimes can be controlled optimizing the energy barrier from the emissive triplet metal-to-ligand charge-transfer (3 MLCT) state to the thermally-activated triplet metal-centered (3 MC) state or the energy gap between both states is a misconception. Further, we demonstrate that considering a single relaxation pathway determined from the minimum that is lowest in energy leads to wrong temperature-dependent emission lifetimes predictions. Instead, we obtain excellent agreement with experimental temperature-dependent lifetimes when an extended kinetic model that includes all the pathways related to multiple Jahn-Teller isomers and their effective reaction barriers is employed. These concepts are essential to correctly design other luminescent transition metal complexes with tailored emission lifetimes based on theoretical predictions.
Collapse
Affiliation(s)
- David Hernández‐Castillo
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Str. 171090ViennaAustria
- Doctoral School in Chemistry (DoSChem)University of ViennaWähringer Straße 421090ViennaAustria
| | - Roland E. P. Nau
- Institute of Inorganic Chemistry IUlm UniversityAlbert-Einstein-Allee 1189081UlmGermany
| | - Marie‐Ann Schmid
- Technische Universität BraunschweigDepartment of Energy Conversion, Institute of Physical and Theoretical ChemistryRebenring 3138106BraunschweigGermany
| | - Stefanie Tschierlei
- Technische Universität BraunschweigDepartment of Energy Conversion, Institute of Physical and Theoretical ChemistryRebenring 3138106BraunschweigGermany
| | - Sven Rau
- Institute of Inorganic Chemistry IUlm UniversityAlbert-Einstein-Allee 1189081UlmGermany
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Str. 171090ViennaAustria
- Vienna Research Platform Accelerating Photoreaction DiscoveryUniversity of ViennaWähringer Straße 171090ViennaAustria
| |
Collapse
|
6
|
Radoń M. Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data? Phys Chem Chem Phys 2023; 25:30800-30820. [PMID: 37938035 DOI: 10.1039/d3cp03537a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.
Collapse
Affiliation(s)
- Mariusz Radoń
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Krakow, Krakow, Poland.
| |
Collapse
|
7
|
Marques EA, De Gendt S, Pourtois G, van Setten MJ. Benchmarking First-Principles Reaction Equilibrium Composition Prediction. Molecules 2023; 28:molecules28093649. [PMID: 37175062 PMCID: PMC10179931 DOI: 10.3390/molecules28093649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Revised: 04/19/2023] [Accepted: 04/20/2023] [Indexed: 05/15/2023] Open
Abstract
The availability of thermochemical properties allows for the prediction of the equilibrium compositions of chemical reactions. The accurate prediction of these can be crucial for the design of new chemical synthesis routes. However, for new processes, these data are generally not completely available. A solution is the use of thermochemistry calculated from first-principles methods such as Density Functional Theory (DFT). Before this can be used reliably, it needs to be systematically benchmarked. Although various studies have examined the accuracy of DFT from an energetic point of view, few studies have considered its accuracy in predicting the temperature-dependent equilibrium composition. In this work, we collected 117 molecules for which experimental thermochemical data were available. From these, we constructed 2648 reactions. These experimentally constructed reactions were then benchmarked against DFT for 6 exchange-correlation functionals and 3 quality of basis sets. We show that, in reactions that do not show temperature dependence in the equilibrium composition below 1000 K, over 90% are predicted correctly. Temperature-dependent equilibrium compositions typically demonstrate correct qualitative behavior. Lastly, we show that the errors are equally caused by errors in the vibrational spectrum and the DFT electronic ground state energy.
Collapse
Affiliation(s)
- Esteban A Marques
- Department of Chemistry, KU Leuven (University of Leuven), Celestijnenlaan 200 F, 3001 Heverlee, Belgium
- IMEC, Kapeldreef 75, 3001 Leuven, Belgium
| | - Stefan De Gendt
- Department of Chemistry, KU Leuven (University of Leuven), Celestijnenlaan 200 F, 3001 Heverlee, Belgium
- IMEC, Kapeldreef 75, 3001 Leuven, Belgium
| | | | - Michiel J van Setten
- IMEC, Kapeldreef 75, 3001 Leuven, Belgium
- ETSF European Theoretical Spectroscopy Facility, Institut de Physique, Université de Liège, Allée du 6 août 17, 4000 Liège, Belgium
| |
Collapse
|
8
|
Samsonova I, Tucker GB, Alaal N, Brorsen KR. Hydrogen-Atom Electronic Basis Sets for Multicomponent Quantum Chemistry. ACS OMEGA 2023; 8:5033-5041. [PMID: 36777583 PMCID: PMC9910068 DOI: 10.1021/acsomega.2c07782] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 01/13/2023] [Indexed: 06/18/2023]
Abstract
Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.
Collapse
|
9
|
Manian A, Hudson RJ, Ramkissoon P, Smith TA, Russo SP. Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants. J Chem Theory Comput 2023; 19:271-292. [PMID: 36490305 DOI: 10.1021/acs.jctc.2c00888] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S1 → S0 transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 106-1016 s-1 compared with an experimental rate of the order 107 s-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.
Collapse
Affiliation(s)
- Anjay Manian
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne3000, Australia
| | - Rohan J Hudson
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Melbourne, Parkville3010, Australia
| | - Pria Ramkissoon
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Melbourne, Parkville3010, Australia
| | - Trevor A Smith
- ARC Centre of Excellence in Exciton Science, School of Chemistry, University of Melbourne, Parkville3010, Australia
| | - Salvy P Russo
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne3000, Australia
| |
Collapse
|
10
|
Liu Y, Zhang C, Liu Z, Truhlar DG, Wang Y, He X. Supervised learning of a chemistry functional with damped dispersion. NATURE COMPUTATIONAL SCIENCE 2023; 3:48-58. [PMID: 38177952 PMCID: PMC10766516 DOI: 10.1038/s43588-022-00371-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2022] [Accepted: 11/11/2022] [Indexed: 01/06/2024]
Abstract
Kohn-Sham density functional theory is widely used in chemistry, but no functional can accurately predict the whole range of chemical properties, although recent progress by some doubly hybrid functionals comes close. Here, we optimized a singly hybrid functional called CF22D with higher across-the-board accuracy for chemistry than most of the existing non-doubly hybrid functionals by using a flexible functional form that combines a global hybrid meta-nonseparable gradient approximation that depends on density and occupied orbitals with a damped dispersion term that depends on geometry. We optimized this energy functional by using a large database and performance-triggered iterative supervised training. We combined several databases to create a very large, combined database whose use demonstrated the good performance of CF22D on barrier heights, isomerization energies, thermochemistry, noncovalent interactions, radical and nonradical chemistry, small and large systems, simple and complex systems and transition-metal chemistry.
Collapse
Affiliation(s)
- Yiwei Liu
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| | - Cheng Zhang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, China
| | - Zhonghua Liu
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN, USA.
| | - Ying Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha, China.
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China.
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai, China.
| |
Collapse
|
11
|
Benchmarking the Computed Proton Solvation Energy and Absolute Potential in Non-aqueous Solvents. Electrochim Acta 2023. [DOI: 10.1016/j.electacta.2022.141785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
|
12
|
Teale AM, Helgaker T, Savin A, Adamo C, Aradi B, Arbuznikov AV, Ayers PW, Baerends EJ, Barone V, Calaminici P, Cancès E, Carter EA, Chattaraj PK, Chermette H, Ciofini I, Crawford TD, De Proft F, Dobson JF, Draxl C, Frauenheim T, Fromager E, Fuentealba P, Gagliardi L, Galli G, Gao J, Geerlings P, Gidopoulos N, Gill PMW, Gori-Giorgi P, Görling A, Gould T, Grimme S, Gritsenko O, Jensen HJA, Johnson ER, Jones RO, Kaupp M, Köster AM, Kronik L, Krylov AI, Kvaal S, Laestadius A, Levy M, Lewin M, Liu S, Loos PF, Maitra NT, Neese F, Perdew JP, Pernal K, Pernot P, Piecuch P, Rebolini E, Reining L, Romaniello P, Ruzsinszky A, Salahub DR, Scheffler M, Schwerdtfeger P, Staroverov VN, Sun J, Tellgren E, Tozer DJ, Trickey SB, Ullrich CA, Vela A, Vignale G, Wesolowski TA, Xu X, Yang W. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. Phys Chem Chem Phys 2022; 24:28700-28781. [PMID: 36269074 PMCID: PMC9728646 DOI: 10.1039/d2cp02827a] [Citation(s) in RCA: 59] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 08/09/2022] [Indexed: 12/13/2022]
Abstract
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.
Collapse
Affiliation(s)
- Andrew M. Teale
- School of Chemistry, University of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, 4 Place Jussieu, CEDEX 05, 75252 Paris, France.
| | - Carlo Adamo
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany.
| | - Alexei V. Arbuznikov
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7Straße des 17. Juni 13510623Berlin
| | | | - Evert Jan Baerends
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy.
| | - Patrizia Calaminici
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Eric Cancès
- CERMICS, Ecole des Ponts and Inria Paris, 6 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France.
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering and the Andlinger Center for Energy and the Environment, Princeton UniversityPrincetonNJ 08544-5263USA
| | | | - Henry Chermette
- Institut Sciences Analytiques, Université Claude Bernard Lyon1, CNRS UMR 5280, 69622 Villeurbanne, France.
| | - Ilaria Ciofini
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - T. Daniel Crawford
- Department of Chemistry, Virginia TechBlacksburgVA 24061USA,Molecular Sciences Software InstituteBlacksburgVA 24060USA
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | | | - Claudia Draxl
- Institut für Physik and IRIS Adlershof, Humboldt-Universität zu Berlin, 12489 Berlin, Germany. .,Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany. .,Beijing Computational Science Research Center (CSRC), 100193 Beijing, China.,Shenzhen JL Computational Science and Applied Research Institute, 518110 Shenzhen, China
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France.
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile.
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, USA.
| | - Giulia Galli
- Pritzker School of Molecular Engineering and Department of Chemistry, The University of Chicago, Chicago, IL, USA.
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China. .,Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | - Nikitas Gidopoulos
- Department of Physics, Durham University, South Road, Durham DH1 3LE, UK.
| | - Peter M. W. Gill
- School of Chemistry, University of SydneyCamperdown NSW 2006Australia
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Andreas Görling
- Chair of Theoretical Chemistry, University of Erlangen-Nuremberg, Egerlandstrasse 3, 91058 Erlangen, Germany.
| | - Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Gold Coast, Qld 4222, Australia.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany.
| | - Oleg Gritsenko
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark.
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie UniversityHalifaxNova ScotiaB3H 4R2Canada
| | - Robert O. Jones
- Peter Grünberg Institut PGI-1, Forschungszentrum Jülich52425 JülichGermany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623, Berlin.
| | - Andreas M. Köster
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav)CDMX07360Mexico
| | - Leeor Kronik
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth, 76100, Israel.
| | - Anna I. Krylov
- Department of Chemistry, University of Southern CaliforniaLos AngelesCalifornia 90089USA
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andre Laestadius
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Mel Levy
- Department of Chemistry, Tulane University, New Orleans, Louisiana, 70118, USA.
| | - Mathieu Lewin
- CNRS & CEREMADE, Université Paris-Dauphine, PSL Research University, Place de Lattre de Tassigny, 75016 Paris, France.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA. .,Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Neepa T. Maitra
- Department of Physics, Rutgers University at Newark101 Warren StreetNewarkNJ 07102USA
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.
| | - John P. Perdew
- Departments of Physics and Chemistry, Temple UniversityPhiladelphiaPA 19122USA
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.
| | - Pascal Pernot
- Institut de Chimie Physique, UMR8000, CNRS and Université Paris-Saclay, Bât. 349, Campus d'Orsay, 91405 Orsay, France.
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA. .,Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| | - Elisa Rebolini
- Institut Laue Langevin, 71 avenue des Martyrs, 38000 Grenoble, France.
| | - Lucia Reining
- Laboratoire des Solides Irradiés, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, F-91120 Palaiseau, France. .,European Theoretical Spectroscopy Facility
| | - Pina Romaniello
- Laboratoire de Physique Théorique (UMR 5152), Université de Toulouse, CNRS, UPS, France.
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA.
| | - Dennis R. Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS – Centre for Molecular Simulation, IQST – Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary2500 University Drive NWCalgaryAlbertaT2N 1N4Canada
| | - Matthias Scheffler
- The NOMAD Laboratory at the FHI of the Max-Planck-Gesellschaft and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195, Germany.
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand.
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western OntarioLondonOntario N6A 5B7Canada
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA.
| | - Erik Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - David J. Tozer
- Department of Chemistry, Durham UniversitySouth RoadDurhamDH1 3LEUK
| | - Samuel B. Trickey
- Quantum Theory Project, Deptartment of Physics, University of FloridaGainesvilleFL 32611USA
| | - Carsten A. Ullrich
- Department of Physics and Astronomy, University of MissouriColumbiaMO 65211USA
| | - Alberto Vela
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Giovanni Vignale
- Department of Physics, University of Missouri, Columbia, MO 65203, USA.
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, Université de Genève30 Quai Ernest-Ansermet1211 GenèveSwitzerland
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China.
| | - Weitao Yang
- Department of Chemistry and Physics, Duke University, Durham, NC 27516, USA.
| |
Collapse
|
13
|
Bursch M, Mewes J, Hansen A, Grimme S. Best-Practice DFT Protocols for Basic Molecular Computational Chemistry. Angew Chem Int Ed Engl 2022; 61:e202205735. [PMID: 36103607 PMCID: PMC9826355 DOI: 10.1002/anie.202205735] [Citation(s) in RCA: 114] [Impact Index Per Article: 57.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2022] [Indexed: 01/11/2023]
Abstract
Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.
Collapse
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für KohlenforschungKaiser-Wilhelm-Platz 145470Mülheim an der RuhrGermany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Andreas Hansen
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Stefan Grimme
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| |
Collapse
|
14
|
Bursch M, Mewes J, Hansen A, Grimme S. Best‐Practice DFT Protocols for Basic Molecular Computational Chemistry**. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202205735] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| |
Collapse
|
15
|
Chan B, Karton A. Assessment of DLPNO-CCSD(T)-F12 and its use for the formulation of the low-cost and reliable L-W1X composite method. J Comput Chem 2022; 43:1394-1402. [PMID: 35709311 DOI: 10.1002/jcc.26892] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 03/28/2022] [Accepted: 04/28/2022] [Indexed: 02/03/2023]
Abstract
In the present study, we have investigated the performance of RIJCOSX DLPNO-CCSD(T)-F12 methods for a wide range of systems. Calculations with a high-accuracy option ["DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12"] extrapolated to the complete-basis-set limit using the maug-cc-pV[D+d,T+d]Z basis sets provides fairly good agreements with the canonical CCSD(T)/CBS reference for a diverse set of thermochemical and kinetic properties [with mean absolute deviations (MADs) of ~1-2 kJ mol-1 except for atomization energies]. On the other hand, the low-cost "RIJCOSX DLPNO-CCSD(T)-F12D" option leads to substantial deviations for certain properties, notably atomization energies (MADs of up to tens of kJ mol-1 ). With the high-accuracy CBS approach, we have formulated the L-W1X method, which further includes a low-cost core-valence plus scalar-relativistic term. It shows generally good accuracy. For improved accuracies in specific cases, we advise replacing maug-cc-pV(n+d)Z with jun-cc-pV(n+d)Z for the calculation of electron affinities, and using well-constructed isodesmic-type reactions to obtain atomization energies. For medium-sized systems, DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12 calculations are several times faster than the corresponding canonical computation; the use of the local approximations (RIJCOSX and DLPNO) leads to a better scaling than that for the canonical calculation (from ~6-7 down to ~2-4 for our test systems). Thus, the DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12 method, and the L-W1X protocol that based on it, represent a useful means for obtaining accurate thermochemical quantities for larger systems.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Nagasaki, Japan
| | - Amir Karton
- School of Molecular Sciences, University of Western Australia, Perth, WA, Australia
| |
Collapse
|
16
|
Theoretical study about the hydrogen abstraction reactions on methyl acetate on combustion conditions. J Mol Model 2022; 28:226. [DOI: 10.1007/s00894-022-05227-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Accepted: 07/08/2022] [Indexed: 11/26/2022]
|
17
|
Levina EO, Khrenova MG, Astakhov AA, Tsirelson VG. Keto-enol tautomerism from the electron delocalization perspective. J Comput Chem 2022; 43:1000-1010. [PMID: 35411548 DOI: 10.1002/jcc.26858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 03/16/2022] [Accepted: 03/18/2022] [Indexed: 11/10/2022]
Abstract
The equilibrium between keto and enol forms in acetylacetone and its derivatives is studied using electron delocalization indices and delocalization tensor density. We demonstrate how electron delocalization governs the equilibrium between keto and enol forms. The less stable enols have more distinct double and single bond character in the CCC fragment, while electron delocalization in this fragment is more pronounced in more stable enols. Looking for the origin of such behavior, we considered the one-electron potentials entering the Euler equation for the electron density. We found that electron delocalization is mainly governed by the static exchange potential, which depends on the three-dimensional atomic structure. It, however, does not distinguish differences in electron delocalization in more and less stable enols, the effect arising from the kinetic exchange contribution, which reflects spin-dependent effects in the electron motion. The local depletion of kinetic exchange in the conjugated fragment yields the enhanced electron delocalization along the CCC bonds in more stable enols. Thus, a combination of considered descriptors allowed us to reveal the influence of electron delocalization on the equilibrium between keto and enol forms and showed the significant features of this phenomenon.
Collapse
Affiliation(s)
- Elena O Levina
- Bach Institute of Biochemistry, Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences, Moscow, Russia.,Laboratory of Supercomputer Methods in Condensed Matter Physics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Maria G Khrenova
- Bach Institute of Biochemistry, Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences, Moscow, Russia.,Chemistry Department, Lomonosov Moscow State University, Moscow, Russia
| | - Andrey A Astakhov
- Flerov Laboratory of Nuclear Reactions, Joint Institute for Nuclear Research, Dubna, Russia
| | - Vladimir G Tsirelson
- Department of Quantum Chemistry, Mendeleev University of Chemical Technology of Russia, Moscow, Russia.,Research Laboratory of Multiscale Modelling of Polyfunctional Compounds, South Ural State University, Chelyabinsk, Russia
| |
Collapse
|
18
|
Gould T, Dale SG. Poisoning density functional theory with benchmark sets of difficult systems. Phys Chem Chem Phys 2022; 24:6398-6403. [PMID: 35244641 DOI: 10.1039/d2cp00268j] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Large benchmark sets like GMTKN55 [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184] let us analyse the performance of density functional theory over a diverse range of systems and bonding types. However, assessing over a large and diverse set can miss cases where approaches fail badly, and can give a misleading sense of security. To this end we introduce a series of 'poison' benchmark sets, P30-5, P30-10 and P30-20, comprising systems with up to 5, 10 and 20 atoms, respectively. These sets represent the most difficult-to-model systems in GMTKN55. We expect them to be useful in developing new approximations, identifying weak points in existing ones, and to aid in selecting appropriate DFAs for computational studies involving difficult physics, e.g. catalysis.
Collapse
Affiliation(s)
- Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111, Australia.
| | - Stephen G Dale
- Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111, Australia.
| |
Collapse
|
19
|
Tarrago M, Ye S, Neese F. Electronic structure analysis of electrochemical CO2 reduction by iron-porphyrins reveals basic requirements to design catalysts bearing non-innocent ligands. Chem Sci 2022; 13:10029-10047. [PMID: 36128248 PMCID: PMC9430493 DOI: 10.1039/d2sc01863b] [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: 03/31/2022] [Accepted: 06/22/2022] [Indexed: 11/21/2022] Open
Abstract
Electrocatalytic CO2 reduction is a possible solution to the increasing CO2 concentration in the earth atmosphere, because it enables storage of energy while using the harmful CO2 feedstock as starting...
Collapse
Affiliation(s)
- Maxime Tarrago
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 D-45470 Mülheim an der Ruhr Germany
| | - Shengfa Ye
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 D-45470 Mülheim an der Ruhr Germany
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics, Chinese Academy of Sciences Dalian 116023 China
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 D-45470 Mülheim an der Ruhr Germany
| |
Collapse
|
20
|
Affiliation(s)
- Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie Vladimir-Prelog-Weg 2 8093 Zürich Switzerland
| |
Collapse
|
21
|
Kirkpatrick J, McMorrow B, Turban DHP, Gaunt AL, Spencer JS, Matthews AGDG, Obika A, Thiry L, Fortunato M, Pfau D, Castellanos LR, Petersen S, Nelson AWR, Kohli P, Mori-Sánchez P, Hassabis D, Cohen AJ. Pushing the frontiers of density functionals by solving the fractional electron problem. Science 2021; 374:1385-1389. [PMID: 34882476 DOI: 10.1126/science.abj6511] [Citation(s) in RCA: 110] [Impact Index Per Article: 36.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
[Figure: see text].
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Louis Thiry
- Département d'informatique, ENS, CNRS, PSL University, Paris, France
| | | | - David Pfau
- DeepMind, 6 Pancras Square, London N1C 4AG, UK
| | | | | | | | | | | | | | - Aron J Cohen
- DeepMind, 6 Pancras Square, London N1C 4AG, UK.,Max Planck Institute for Solid State Research, 70569 Stuttgart, Germany
| |
Collapse
|
22
|
Yalcin-Ozkat G. Molecular Modeling Strategies of Cancer Multidrug Resistance. Drug Resist Updat 2021; 59:100789. [PMID: 34973929 DOI: 10.1016/j.drup.2021.100789] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Revised: 11/20/2021] [Accepted: 11/22/2021] [Indexed: 02/07/2023]
Abstract
Cancer remains a leading cause of morbidity and mortality worldwide. Hence, the increase in cancer cases observed in the elderly population, as well as in children and adolescents, makes human malignancies a prime target for anticancer drug development. Although highly effective chemotherapeutic agents are continuously developed and approved for clinical treatment, the major impediment towards curative cancer therapy remains multidrug resistance (MDR). In recent years, intensive studies have been carried out on the identification of new therapeutic molecules to reverse MDR efflux transporters of the ATP-binding cassette (ABC) superfamily. Although a great deal of progress has been made in the development of specific inhibitors for certain MDR efflux pumps in experimental studies, advanced computational studies can accelerate this drug development process. In the literature, there are many experimental studies on the impact of natural products and synthetic small molecules on the reversal of cancer MDR. Molecular modeling methods provide an opportunity to explain the activity of these molecules on the ABC-transporter family with non-covalent interactions as well as it is possible to carry out studies for the discovery of new anticancer drugs specific to MDR with these methods. The coordinate file of the 3-dimensional (3D) structure of the target protein is indispensable for molecular modeling studies. In some cases where a 3D structure cannot be obtained by experimental methods, the homology modeling method can be applied to obtain the file containing the target protein's information including atomic coordinates, secondary structure assignments, and atomic connectivity. Homology modeling studies are of great importance for efflux transporter proteins that still lack 3D structures due to crystallization problems with multiple hydrophobic transmembrane domains. Quantum mechanics, molecular docking and molecular dynamics simulation applications are the most frequently used molecular modeling methods in the literature to investigate non-covalent interactions between the drug-ABC transporter superfamily. The quantitative structure-activity relationship (QSAR) model provides a relationship between the chemical properties of a compound and its biological activity. Determining the pharmacophore region for a new drug molecule by superpositioning a series of molecules according to their physicochemical properties using QSAR models is another method in which molecular modeling is used in computational drug development studies with ABC transporter proteins. There are also in silico absorption, distribution, metabolism, excretion, and toxicity (ADME/Tox) studies conducted to make a prediction about the pharmacokinetic properties, and drug-likeness of new molecules. Drug repurposing studies, which have become a trending topic in recent years, involve identifying possible new targets for an already approved drug molecule. There are few studies in the literature in which drug repurposing performed by molecular modelling methods has been applied on ABC transporter proteins. The aim of the current paper is to create a complete review of drug development studies including aforementioned molecular modeling methods carried out between the years 2019-2021. Furthermore, an intensive investigation is also conducted on licensed applications and free web servers used in in silico studies. The current review is an up-to-date guide for researchers who plan to conduct computational studies with MDR transporter proteins.
Collapse
Affiliation(s)
- Gozde Yalcin-Ozkat
- Recep Tayyip Erdogan University, Faculty of Engineering and Architecture, Bioengineering Department, 53100, Rize, Turkey; Max Planck Institute for Dynamics of Complex Technical Systems, Molecular Simulations and Design Group, Sandtorstrasse 1, 39106, Magdeburg, Germany.
| |
Collapse
|
23
|
Wu X, Hou Q, Huang J, Chai J, Zhang F. Exploring the OH-initiated reactions of styrene in the atmosphere and the role of van der Waals complex. CHEMOSPHERE 2021; 282:131004. [PMID: 34082313 DOI: 10.1016/j.chemosphere.2021.131004] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 05/19/2021] [Accepted: 05/23/2021] [Indexed: 06/12/2023]
Abstract
Reacting with OH provides a major sink for styrene in the atmosphere, with three possible pathways including OH-addition, H-abstraction and addition-dissociation reactions. However, the total rate coefficients of styrene + OH were measured as 1.2-6.2 × 10-11 cm3 molecule-1 s-1 under atmospheric conditions, varying by a maximum factor of 5. On the other hand, only one theoretical work reported this rate coefficient as 19.1 × 10-11 cm3 molecule-1 s-1, which exhibits up to 16 times that measured in laboratory studies. In the present study, the reaction kinetics of styrene + OH was extensively studied with high-level quantum chemical methods combined with RRKM/master equation simulations. In particular, we carried out theoretical treatments for the formation of pre-reaction Van der Waals complexes of styrene + OH, and examined their influence on the reaction kinetics. The total rate coefficient for styrene + OH is calculated to be 1.7 × 10-11 cm3 molecule-1 s-1 at 300 K, 1 atm. The main products are addβ (88.2%), add5 (6.9%), addα (1.9%) and add3 (1.7%). Using our computed rate coefficient and the global atmospheric hydroxyl radical concentration (2 × 106 radicals per cm3), the lifetime of styrene in the atmosphere is estimated at 8.0 h. The degradation of styrene might be negligible for the formation of ozone in the atmosphere based upon the photochemical ozone creation potentials calculation. The computed product yields indicate that addβ via subsequent reactions could significantly produce formaldehyde and benzaldehyde that were observed in previous experimental studies on styrene oxidation, and contribute to the formation of secondary organic aerosols.
Collapse
Affiliation(s)
- Xiaoqing Wu
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, PR China; National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei, Anhui, 230029, PR China
| | - Qifeng Hou
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, PR China
| | - Jiabin Huang
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, PR China
| | - Jiajue Chai
- Institute at Brown for Environment and Society, And Department of Earth, Environmental and Planetary Sciences, Brown University, 182 Hope St., Providence, RI, 02912, USA
| | - Feng Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui, 230026, PR China.
| |
Collapse
|
24
|
Abstract
Chemical compound space (CCS), the set of all theoretically conceivable combinations of chemical elements and (meta-)stable geometries that make up matter, is colossal. The first-principles based virtual sampling of this space, for example, in search of novel molecules or materials which exhibit desirable properties, is therefore prohibitive for all but the smallest subsets and simplest properties. We review studies aimed at tackling this challenge using modern machine learning techniques based on (i) synthetic data, typically generated using quantum mechanics based methods, and (ii) model architectures inspired by quantum mechanics. Such Quantum mechanics based Machine Learning (QML) approaches combine the numerical efficiency of statistical surrogate models with an ab initio view on matter. They rigorously reflect the underlying physics in order to reach universality and transferability across CCS. While state-of-the-art approximations to quantum problems impose severe computational bottlenecks, recent QML based developments indicate the possibility of substantial acceleration without sacrificing the predictive power of quantum mechanics.
Collapse
Affiliation(s)
- Bing Huang
- Faculty
of Physics, University of Vienna, 1090 Vienna, Austria
| | - O. Anatole von Lilienfeld
- Faculty
of Physics, University of Vienna, 1090 Vienna, Austria
- Institute
of Physical Chemistry and National Center for Computational Design
and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, 4056 Basel, Switzerland
| |
Collapse
|
25
|
Ramos C, Muehlbrad J, Janesko BG. Density functionals with full nonlocal exchange, nonlocal rung-3.5 correlation, and D3 dispersion: Combined accuracy for general main-group thermochemistry, kinetics, and noncovalent interactions. J Comput Chem 2021; 42:1974-1981. [PMID: 34387364 DOI: 10.1002/jcc.26728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Accepted: 07/25/2021] [Indexed: 11/11/2022]
Abstract
We introduce the HF-R35-D3(BJ) functional combining full nonlocal exact (Hartree-Fock-like, HF) exchange, inexpensive rung-3.5 correlation constructed from nonlocal one-electron operators, and nonlocal D3 dispersion corrections. HF-R35-D3(BJ) is among the first full-exact-exchange functionals offering competitive accuracy for general main-group thermochemistry, kinetics, and noncovalent interactions. HF-R35-D3(BJ) gives weighted mean absolute deviation WTMAD-2 8.5 kcal/mol across the entire GMTKN55 dataset, outperforming most dispersion-corrected semilocal functionals and approaching the accuracy of dispersion-corrected global hybrids. This requires six fitted parameters, three each in the nonlocal correlation and dispersion corrections. Full nonlocal exchange appears to help give accurate binding energies and reasonable energy orderings for water hexamers. These results motivate continued exploration of inexpensive nonlocal correlation corrections to nonlocal exchange.
Collapse
Affiliation(s)
- Chloe Ramos
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas, USA
| | - Jeremiah Muehlbrad
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas, USA
| | - Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas, USA
| |
Collapse
|
26
|
Bakowies D, von Lilienfeld OA. Density Functional Geometries and Zero-Point Energies in Ab Initio Thermochemical Treatments of Compounds with First-Row Atoms (H, C, N, O, F). J Chem Theory Comput 2021; 17:4872-4890. [PMID: 34260240 PMCID: PMC8437339 DOI: 10.1021/acs.jctc.1c00474] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2021] [Indexed: 12/02/2022]
Abstract
Density functionals are often used in ab initio thermochemistry to provide optimized geometries for single-point evaluations at a high level and to supply estimates of anharmonic zero-point energies (ZPEs). Their use is motivated by relatively high accuracy at a modest computational expense, but a thorough assessment of geometry-related error seems to be lacking. We have benchmarked 53 density functionals, focusing on approximations of the first four rungs and on relatively small basis sets for computational efficiency. Optimized geometries of 279 neutral first-row molecules (H, C, N, O, F) are judged by energy penalties relative to the best available geometries, using the composite model ATOMIC/B5 as energy probe. Only hybrid functionals provide good accuracy with root-mean-square errors around 0.1 kcal/mol and maximum errors below 1.0 kcal/mol, but not all of them do. Conspicuously, first-generation hybrids with few or no empirical parameters tend to perform better than highly parameterized ones. A number of them show good accuracy already with small basis sets (6-31G(d), 6-311G(d)). As is standard practice, anharmonic ZPEs are estimated from scaled harmonic values. Statistics of the latter show less performance variation among functionals than observed for geometry-related error, but they also indicate that ZPE error will generally dominate. We have selected PBE0-D3/6-311G(d) for the next version of the ATOMIC protocol (ATOMIC-2) and studied it in more detail. Empirical expressions have been calibrated to estimate bias corrections and 95% uncertainty intervals for both geometry-related error and scaled ZPEs.
Collapse
Affiliation(s)
- Dirk Bakowies
- Institute
of Physical Chemistry, Department of Chemistry, University of Basel, Klingelbergstr. 80, CH 4056 Basel, Switzerland
| | - O. Anatole von Lilienfeld
- Faculty
of Physics, University of Vienna, Kolingasse 14-16, A 1090 Vienna, Austria
- Institute
of Physical Chemistry and National Center for Computational Design
and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstr. 80, CH 4056 Basel, Switzerland
| |
Collapse
|
27
|
Grimme S, Hansen A, Ehlert S, Mewes JM. r 2SCAN-3c: A "Swiss army knife" composite electronic-structure method. J Chem Phys 2021; 154:064103. [PMID: 33588555 DOI: 10.1063/5.0040021] [Citation(s) in RCA: 235] [Impact Index Per Article: 78.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The recently proposed r2SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r2SCAN-3c. To this end, the unaltered r2SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r2SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r2SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.
Collapse
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
28
|
Semidalas E, Martin JML. Canonical and DLPNO-Based Composite Wavefunction Methods Parametrized against Large and Chemically Diverse Training Sets. 2: Correlation-Consistent Basis Sets, Core-Valence Correlation, and F12 Alternatives. J Chem Theory Comput 2020; 16:7507-7524. [PMID: 33200931 PMCID: PMC7735707 DOI: 10.1021/acs.jctc.0c01106] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A hierarchy
of wavefunction composite methods (cWFT), based on
G4-type cWFT methods available for elements H through Rn, was recently
reported by the present authors [2020, 16, 4238]. We extend this hierarchy
by considering the inner-shell correlation energy in the second-order
Møller–Plesset correction and replacing the Weigend–Ahlrichs
def2-mZVPP(D) basis sets used with complete basis
set extrapolation from augmented correlation-consistent core–valence
triple-ζ, aug-cc-pwCVTZ(-PP), and quadruple-ζ, aug-cc-pwCVQZ(-PP),
basis sets, thus creating cc-G4-type methods. For the large and chemically
diverse GMTKN55 benchmark suite, they represent a substantial further
improvement and bring WTMAD2 (weighted mean absolute deviation) down
below 1 kcal/mol. Intriguingly, the lion’s share of the improvement
comes from better capture of valence correlation; the inclusion of
core–valence correlation is almost an order of magnitude less
important. These robust correlation-consistent cWFT methods approach
the CCSD(T) complete basis limit with just one or a few fitted parameters.
Particularly, the DLPNO variants such as cc-G4-T-DLPNO are applicable
to fairly large molecules at a modest computational cost, as is (for
a reduced range of elements) a different variant using MP2-F12/cc-pVTZ-F12
for the MP2 component.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Rehovot, Israel
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Rehovot, Israel
| |
Collapse
|
29
|
Benali A, Shin H, Heinonen O. Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set. J Chem Phys 2020; 153:194113. [DOI: 10.1063/5.0026275] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Anouar Benali
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Hyeondeok Shin
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Olle Heinonen
- Material Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| |
Collapse
|
30
|
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
|
31
|
Wu X, Huang C, Niu S, Zhang F. New theoretical insights into the reaction kinetics of toluene and hydroxyl radicals. Phys Chem Chem Phys 2020; 22:22279-22288. [DOI: 10.1039/d0cp02984j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This work provides theoretical insights into the kinetics of toluene + OH, focusing on the anharmonic effect and the accuracy of barrier heights.
Collapse
Affiliation(s)
- Xiaoqing Wu
- National Synchrotron Radiation Laboratory
- University of Science and Technology of China
- Hefei
- P. R. China
- Hefei National Laboratory for Physical Sciences at the Microscale
| | - Can Huang
- Chair of Technical Thermodynamics
- RWTH Aachen University
- 52062 Aachen
- Germany
| | - Shiyao Niu
- Science and Technology on Combustion and Explosion Laboratory
- Xi'an Modern Chemistry Research Institute
- Xi'an
- P. R. China
- School of Chemistry and Materials Science
| | - Feng Zhang
- Hefei National Laboratory for Physical Sciences at the Microscale
- University of Science and Technology of China
- Hefei
- P. R. China
| |
Collapse
|
32
|
Martin JML, Santra G. Empirical Double‐Hybrid Density Functional Theory: A ‘Third Way’ in Between WFT and DFT. Isr J Chem 2019. [DOI: 10.1002/ijch.201900114] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Jan M. L. Martin
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| | - Golokesh Santra
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| |
Collapse
|
33
|
Rakshit A, Bandyopadhyay P, Heindel JP, Xantheas SS. Atlas of putative minima and low-lying energy networks of water clusters n = 3-25. J Chem Phys 2019; 151:214307. [PMID: 31822087 DOI: 10.1063/1.5128378] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report a database consisting of the putative minima and ∼3.2 × 106 local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as low-lying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation (MP2)/aug-cc-pVTZ level of theory. In total, we have identified 3 138 303 networks corresponding to local minima of the clusters n = 3-25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1-3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but not limited to neural networks and machine and deep learning) to understand the properties of water, nature's most important substance.
Collapse
Affiliation(s)
- Avijit Rakshit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
34
|
Moulder CA, Kafle K, Cundari TR. Tungsten-Ligand Bond Strengths for 2p Elements Including σ- and π-Bond Strength Components, A Density Functional Theory and ab Initio Study. J Phys Chem A 2019; 123:7940-7949. [PMID: 31240921 DOI: 10.1021/acs.jpca.9b03272] [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/29/2022]
Abstract
Three WVI crystal structures with multifarious metal-ligand bond types are used to theoretically predict homolytic metal-element bond enthalpies with 11 popular DFT functionals, MP2 wave function methods, and four common valence basis set/pseudopotentials in order to evaluate the accuracy and precision of the resultant bond enthalpy data. To our knowledge, for the first time, estimates of component metal-ligand σ- and π-bond strengths are computed. The WE (E = C, N, O) bond enthalpies have the consistent trend σ > second π > first π. In contrast, the element-element BDE trend for the 2p homologues is second π > first π > σ for nitrogen and oxygen, and σ > first π > second π for carbon. These differences may underpin the differences in stability trends and thus reactivity behavior for metal-element multiple bonds as compared to the element-element multiple bonds, and metal-element triple bonds versus their corresponding double bonded counterparts. For example, Odom et al. show that MeI nucleophilically attacks at the imide (M═N) rather than the nitride (M ≡ N) ligand; the relative π-bond strengths derived herein provide a thermodynamic rationalization for this site preference. In this study, it is deduced from the calculated thermodynamics that the W-oxo ligand is more congruous with a triple bond than a double bond, consistent with the bonding model set forth in the seminal 1961 Ballhausen-Gray paper.
Collapse
Affiliation(s)
- Catherine A Moulder
- Department of Chemistry & Center for Advanced Scientific Computing and Modeling (CASCaM) , University of North Texas , 1155 Union Circle, #305070 , Denton , Texas 76203-5017 , United States
| | - Kristina Kafle
- Department of Chemistry & Center for Advanced Scientific Computing and Modeling (CASCaM) , University of North Texas , 1155 Union Circle, #305070 , Denton , Texas 76203-5017 , United States
| | - Thomas R Cundari
- Department of Chemistry & Center for Advanced Scientific Computing and Modeling (CASCaM) , University of North Texas , 1155 Union Circle, #305070 , Denton , Texas 76203-5017 , United States
| |
Collapse
|
35
|
Chan B, Karton A, Raghavachari K. G4(MP2)-XK: A Variant of the G4(MP2)-6X Composite Method with Expanded Applicability for Main-Group Elements up to Radon. J Chem Theory Comput 2019; 15:4478-4484. [DOI: 10.1021/acs.jctc.9b00449] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo-machi 1-14, Nagasaki 852-8521, Japan
| | - Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
36
|
Kayanuma M, Shoji M, Furuya K, Kamiya K, Aikawa Y, Umemura M, Shigeta Y. First-Principles Study of the Reaction Mechanism of CHO + H on Graphene Surface. J Phys Chem A 2019; 123:5633-5639. [PMID: 31244121 DOI: 10.1021/acs.jpca.9b02345] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Many organic molecules observed in the interstellar medium are considered to be formed on dust grains and populated into the gas phase. We analyzed the reaction of HCO + H on a graphene surface using ab initio molecular dynamics simulations as a case study of the formation and desorption of organic molecules on interstellar dust particles. During the reactions of chemisorbed CHO (chemisorbed at the C atom) with free H, CO was generated and efficiently desorbed from the surface. These results suggest that the reactions, of which the reactant forms a covalent bond with the surface while the product does not, cause efficient desorption of the product upon reaction. In such reactions a repulsive force between the product and the surface would be generated and accelerate translation of the product in a specific direction. In addition, it was also shown that the branching ratio of the reactions between radical species on the surface would be affected by the form of the adsorption on the surface, e.g., when a free H reacted with the CHO chemisorbed at the C atom, CH2O was not generated.
Collapse
Affiliation(s)
- Megumi Kayanuma
- Research Center for Computational Design of Advanced Functional Materials , National Institute of Advanced Industrial Science and Technology , Central 2, 1-1-1 Umezono , Tsukuba , Ibaraki 305-8568 , Japan.,Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| | - Mitsuo Shoji
- Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| | - Kenji Furuya
- Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| | - Katsumasa Kamiya
- Center for Basic Education and Integrated Learning , Kanagawa Institute of Technology , 1030 Shimoogino , Atsugi , Kanagawa 243-0292 , Japan
| | - Yuri Aikawa
- Department of Astronomy , The University of Tokyo , 7-3-1 Hongo , Bunkyo-ku , Tokyo 113-0033 , Japan
| | - Masayuki Umemura
- Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan.,Institute of Space and Astronautical Science , Japan Aerospace Exploration Agency , 3-1-1 Yoshinodai, Chuo-ku , Sagamihara , Kanagawa 252-0222 , Japan
| |
Collapse
|
37
|
Li F, Sa R, Song X. TCA self-assembled fluorescence probe for Cu (II) ion based on the unique configuration of extra nuclear electrons of metal ions: A TDDFT study. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.04.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
38
|
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: 217] [Impact Index Per Article: 43.4] [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
|
39
|
Dral PO, Wu X, Thiel W. Semiempirical Quantum-Chemical Methods with Orthogonalization and Dispersion Corrections. J Chem Theory Comput 2019; 15:1743-1760. [PMID: 30735388 PMCID: PMC6416713 DOI: 10.1021/acs.jctc.8b01265] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Indexed: 12/31/2022]
Abstract
We present two new semiempirical quantum-chemical methods with orthogonalization and dispersion corrections: ODM2 and ODM3 (ODM x). They employ the same electronic structure model as the OM2 and OM3 (OM x) methods, respectively. In addition, they include Grimme's dispersion correction D3 with Becke-Johnson damping and three-body corrections E ABC for Axilrod-Teller-Muto dispersion interactions as integral parts. Heats of formation are determined by adding explicitly computed zero-point vibrational energy and thermal corrections, in contrast to standard MNDO-type and OM x methods. We report ODM x parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine that are optimized with regard to a wide range of carefully chosen state-of-the-art reference data. Extensive benchmarks show that the ODM x methods generally perform better than the available MNDO-type and OM x methods for ground-state and excited-state properties, while they describe noncovalent interactions with similar accuracy as OM x methods with a posteriori dispersion corrections.
Collapse
Affiliation(s)
- Pavlo O. Dral
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Xin Wu
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
40
|
Chiniforoush S, Cramer CJ. Quantum Chemical Characterization of Factors Affecting the Neutral and Radical-Cation Newman-Kwart Reactions. J Org Chem 2019; 84:2148-2157. [PMID: 30672704 DOI: 10.1021/acs.joc.8b03132] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Details of the electronic and geometric structures of stationary points along the reaction coordinate of the Newman Kwart rearrangement, which transforms an O-arylthiocarbamate into an S-arylcarbamothioate, are examined using quantum-chemical methods for a large number of compounds considering both the thermal reactions of uncharged substrates as well as the corresponding reactions of radical-cation substrates generated under photoredox conditions. The uncharged mechanism, which has intrinsically high 298 K free energies of activation (in excess of 30 kcal mol-1), has the character of nucleophilic aromatic substitution and is thus accelerated by electron-withdrawing substituents on the aromatic ring. The radical-cationic mechanism, by contrast, has 298 K free energies of activation that are typically below 20 kcal mol-1 and is accelerated by electron donating substituents on the aromatic ring, which stabilize the hole character that is transferred to this fragment from the thiocarbamate fragment as the reaction proceeds. Opportunities to further accelerate the radical-cation reaction are revealed by computational assessment of alternative amino groups for the thiocarbamate functionality.
Collapse
Affiliation(s)
- Sina Chiniforoush
- Department of Chemistry, Minnesota Supercomputing Institute, and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
| | - Christopher J Cramer
- Department of Chemistry, Minnesota Supercomputing Institute, and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
| |
Collapse
|
41
|
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
|
42
|
Goerigk L, Mehta N. A Trip to the Density Functional Theory Zoo: Warnings and Recommendations for the User. Aust J Chem 2019. [DOI: 10.1071/ch19023] [Citation(s) in RCA: 79] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
This account is written for general users of density functional theory (DFT) methods as well as experimental researchers who are new to the field and would like to conduct such calculations. Its main emphasis lies on how to find a way through the confusing ‘zoo’ of DFT by addressing common misconceptions and highlighting those modern methods that should ideally be used in calculations of energetic properties and geometries. A particular focus is on highly popular methods and the important fact that popularity does not imply accuracy. In this context, we present a new analysis of the openly available data published in Swart and co-workers’ famous annual ‘DFT poll’ (http://www.marcelswart.eu/dft-poll/) to demonstrate the existing communication gap between the DFT user and developer communities. We show that despite considerable methodological advances in the field, the perception of some parts of the user community regarding their favourite approaches has changed little. It is hoped that this account makes a contribution towards changing this status and that users are inspired to adjust their current computational protocols to accommodate strategies that are based on proven robustness, accuracy, and efficiency rather than popularity.
Collapse
|
43
|
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
|
44
|
Simm GN, Reiher M. Error-Controlled Exploration of Chemical Reaction Networks with Gaussian Processes. J Chem Theory Comput 2018; 14:5238-5248. [PMID: 30179500 DOI: 10.1021/acs.jctc.8b00504] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For a theoretical understanding of the reactivity of complex chemical systems, relative energies of stationary points on potential energy hypersurfaces need to be calculated to high accuracy. Due to the large number of intermediates present in all but the simplest chemical processes, approximate quantum chemical methods are required that allow for fast evaluations of the relative energies but at the expense of accuracy. Despite the plethora of benchmark studies, the accuracy of a quantum chemical method is often difficult to assess. Moreover, a significant improvement of a method's accuracy (e.g., through reparameterization or systematic model extension) is rarely possible. Here, we present a new approach that allows for the systematic, problem-oriented, and rolling improvement of quantum chemical results through the application of Gaussian processes. Due to its Bayesian nature, reliable error estimates are provided for each prediction. A reference method of high accuracy can be employed if the uncertainty associated with a particular calculation is above a given threshold. The new data point is then added to a growing data set in order to continuously improve the model and, as a result, all subsequent predictions. Previous predictions are validated by the updated model to ensure that uncertainties remain within the given confidence bound, which we call backtracking. We demonstrate our approach with the example of a complex chemical reaction network.
Collapse
Affiliation(s)
- Gregor N Simm
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| |
Collapse
|
45
|
Chan B. Formulation of Small Test Sets Using Large Test Sets for Efficient Assessment of Quantum Chemistry Methods. J Chem Theory Comput 2018; 14:4254-4262. [DOI: 10.1021/acs.jctc.8b00514] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
46
|
Herr JE, Yao K, McIntyre R, Toth DW, Parkhill J. Metadynamics for training neural network model chemistries: A competitive assessment. J Chem Phys 2018; 148:241710. [PMID: 29960377 DOI: 10.1063/1.5020067] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Neural network model chemistries (NNMCs) promise to facilitate the accurate exploration of chemical space and simulation of large reactive systems. One important path to improving these models is to add layers of physical detail, especially long-range forces. At short range, however, these models are data driven and data limited. Little is systematically known about how data should be sampled, and "test data" chosen randomly from some sampling techniques can provide poor information about generality. If the sampling method is narrow, "test error" can appear encouragingly tiny while the model fails catastrophically elsewhere. In this manuscript, we competitively evaluate two common sampling methods: molecular dynamics (MD), normal-mode sampling, and one uncommon alternative, Metadynamics (MetaMD), for preparing training geometries. We show that MD is an inefficient sampling method in the sense that additional samples do not improve generality. We also show that MetaMD is easily implemented in any NNMC software package with cost that scales linearly with the number of atoms in a sample molecule. MetaMD is a black-box way to ensure samples always reach out to new regions of chemical space, while remaining relevant to chemistry near kbT. It is a cheap tool to address the issue of generalization.
Collapse
Affiliation(s)
- John E Herr
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - Ryker McIntyre
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - David W Toth
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - John Parkhill
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| |
Collapse
|
47
|
Wong ZC, Fan WY, Chwee TS. Computational modelling of singlet excitation energy transfer: a DFT/TD-DFT study of the ground and excited state properties of a syn bimane dimer system using non-empirically tuned range-separated functionals. NEW J CHEM 2018. [DOI: 10.1039/c8nj02920b] [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/21/2022]
Abstract
Range-tuned DFT/TD-DFT improves predicted properties for the sequence of events leading to excitation energy transfer (EET) in bimanes.
Collapse
Affiliation(s)
- Z. C. Wong
- Institute of High Performance Computing
- Agency for Science, Technology and Research (A*STAR)
- Singapore
- NUS Graduate School for Integrative Sciences and Engineering
- National University of Singapore
| | - W. Y. Fan
- Department of Chemistry
- National University of Singapore
- Singapore
| | - T. S. Chwee
- Institute of High Performance Computing
- Agency for Science, Technology and Research (A*STAR)
- Singapore
| |
Collapse
|
48
|
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
|
49
|
McAnanama-Brereton S, Waller MP. Rational Density Functional Selection Using Game Theory. J Chem Inf Model 2017; 58:61-67. [DOI: 10.1021/acs.jcim.7b00542] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Suzanne McAnanama-Brereton
- Theoretische
Organische Chemie, Organisch-Chemisches Institut der Universität Münster, Correnstraße, D-48149 Münster, Germany
| | - Mark P. Waller
- Department
of Physics and International Centre for Quantum and Molecular Structures, Shanghai University, Shanghai 200444, People’s Republic of China
| |
Collapse
|
50
|
Mata RA, Suhm MA. Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction? Angew Chem Int Ed Engl 2017; 56:11011-11018. [PMID: 28452424 PMCID: PMC5582598 DOI: 10.1002/anie.201611308] [Citation(s) in RCA: 100] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2016] [Revised: 12/20/2016] [Indexed: 11/15/2022]
Abstract
Theoreticians and experimentalists should work together more closely to establish reliable rankings and benchmarks for quantum chemical methods. Comparison to carefully designed experimental benchmark data should be a priority. Guidelines to improve the situation for experiments and calculations are proposed.
Collapse
Affiliation(s)
- Ricardo A. Mata
- Institut für Physikalische ChemieUniversität GöttingenTammannstrasse 637077GöttingenGermany
| | - Martin A. Suhm
- Institut für Physikalische ChemieUniversität GöttingenTammannstrasse 637077GöttingenGermany
| |
Collapse
|