1
|
Eikey EA, Maldonado AM, Griego CD, Von Rudorff GF, Keith JA. Evaluating quantum alchemy of atoms with thermodynamic cycles: Beyond ground electronic states. J Chem Phys 2022; 156:064106. [DOI: 10.1063/5.0079483] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Emily A. Eikey
- Chemistry, University of Pittsburgh, United States of America
| | - Alex M. Maldonado
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, United States of America
| | | | | | - John A. Keith
- Dept. of Chemical & Petroleum Engineering, University of Pittsburgh, United States of America
| |
Collapse
|
2
|
Gómez T, Fuentealba P, Robles-Navarro A, Cárdenas C. Links among the Fukui potential, the alchemical hardness and the local hardness of an atom in a molecule. J Comput Chem 2021; 42:1681-1688. [PMID: 34121207 DOI: 10.1002/jcc.26705] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2021] [Revised: 05/24/2021] [Accepted: 05/31/2021] [Indexed: 02/01/2023]
Abstract
This paper presents a brief summary of the difficulty that resides in the definition of the elusive concept of local chemical hardness. We argue that a definition of local hardness should be useful to a reactivity principle and not just as a mere definition. We then continue with a formal discussion about the benefits and difficulties of using the Fukui potential, which is interpreted as an alchemical derivative (alchemical hardness), as descriptor of local hardness of molecules. Computational evidence shows that the alchemical hardness is at least as good a descriptor as the combination of other two well-stabilized descriptors of local hardness, such as the Fukui function and grand canonical local hardness. Although our results are auspicious for the alchemical hardness as descriptor of local hardness, we finish by calling the attention of the community on the importance of discussing the raison d'être of a local hardness function and its main characteristics. We suggest that an axiomatic construction of local hardness could be they way of constructing a local hardness which is both useful and free of arbitrariness.
Collapse
Affiliation(s)
- Tatiana Gómez
- Theoretical and Computational Chemistry Center, Institute of Applied Chemical Sciences, Faculty of Engineering, Universidad Autonoma de Chile, Santiago, Chile
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Santiago, Chile.,Centro para el Desarrollo de la Nanociencia y la Nanotecnología (CEDENNA), Santiago, Chile
| | | | - Carlos Cárdenas
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Santiago, Chile.,Centro para el Desarrollo de la Nanociencia y la Nanotecnología (CEDENNA), Santiago, Chile
| |
Collapse
|
3
|
von Rudorff GF, von Lilienfeld OA. Simplifying inverse materials design problems for fixed lattices with alchemical chirality. SCIENCE ADVANCES 2021; 7:eabf1173. [PMID: 34138735 PMCID: PMC8133750 DOI: 10.1126/sciadv.abf1173] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 03/25/2021] [Indexed: 05/03/2023]
Abstract
Brute-force compute campaigns relying on demanding ab initio calculations routinely search for previously unknown materials in chemical compound space (CCS), the vast set of all conceivable stable combinations of elements and structural configurations. Here, we demonstrate that four-dimensional chirality arising from antisymmetry of alchemical perturbations dissects CCS and defines approximate ranks, which reduce its formal dimensionality and break down its combinatorial scaling. The resulting "alchemical" enantiomers have the same electronic energy up to the third order, independent of respective covalent bond topology, imposing relevant constraints on chemical bonding. Alchemical chirality deepens our understanding of CCS and enables the establishment of trends without empiricism for any materials with fixed lattices. We demonstrate the efficacy for three cases: (i) new rules for electronic energy contributions to chemical bonding; (ii) analysis of the electron density of BN-doped benzene; and (iii) ranking over 2000 and 4 million BN-doped naphthalene and picene derivatives, respectively.
Collapse
Affiliation(s)
- Guido Falk von Rudorff
- University of Vienna, Faculty of Physics, Kolingasse 14-16, 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
| | - O Anatole von Lilienfeld
- University of Vienna, Faculty of Physics, Kolingasse 14-16, 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
|
4
|
Griego CD, Zhao L, Saravanan K, Keith JA. Machine learning corrected alchemical perturbation density functional theory for catalysis applications. AIChE J 2020. [DOI: 10.1002/aic.17041] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Charles D. Griego
- Department of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh Pittsburgh Pennsylvania USA
| | - Lingyan Zhao
- Department of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh Pittsburgh Pennsylvania USA
| | - Karthikeyan Saravanan
- Department of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh Pittsburgh Pennsylvania USA
| | - John A. Keith
- Department of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh Pittsburgh Pennsylvania USA
| |
Collapse
|
5
|
Muñoz M, Robles-Navarro A, Fuentealba P, Cárdenas C. Predicting Deprotonation Sites Using Alchemical Derivatives. J Phys Chem A 2020; 124:3754-3760. [PMID: 32286831 DOI: 10.1021/acs.jpca.9b09472] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
An alchemical transformation is any process, physical or fictitious, that connects two points in the chemical space. A particularly important transformation is the vanishing of a proton, whose energy can be linked to the proton dissociation enthalpy of acids. In this work we assess the reliability of alchemical derivatives in predicting the proton dissociation enthalpy of a diverse series of mono- and polyprotic molecules. Alchemical derivatives perform remarkably well in ranking the proton affinity of all molecules. Additionally, alchemical derivatives could be use also as a predictive tool because their predictions correlate quite well with calculations based on energy differences and experimental values. Although second-order alchemical derivatives underestimate the dissociation enthalpy, the deviation seems to be almost constant. This makes alchemical derivatives extremely accurate to evaluate the difference in proton affinity between two acid sites of polyprotic molecule. Finally, we show that the reason for the underestimation of the dissociation enthalpy is most likely the contribution of higher-order derivatives.
Collapse
Affiliation(s)
- Macarena Muñoz
- Facultad de Ingenierı́a y Ciencias, Universidad Adolfo Ibañez, Diagonal Las Torres 2640, Santiago 7941169, Chile
| | - Andrés Robles-Navarro
- Departamento de Fı́sica, Facultad de Ciencias, Universidad de Chile, Las Palmeras 3425, Santiago Casilla 653, Chile.,Centro para el Desarrollo de la Nanociencia y la Nanotecnologı́a (CEDENNA), Avda. Ecuador 3493, Santiago 9170124, Chile
| | - Patricio Fuentealba
- Departamento de Fı́sica, Facultad de Ciencias, Universidad de Chile, Las Palmeras 3425, Santiago Casilla 653, Chile.,Centro para el Desarrollo de la Nanociencia y la Nanotecnologı́a (CEDENNA), Avda. Ecuador 3493, Santiago 9170124, Chile
| | - Carlos Cárdenas
- Departamento de Fı́sica, Facultad de Ciencias, Universidad de Chile, Las Palmeras 3425, Santiago Casilla 653, Chile.,Centro para el Desarrollo de la Nanociencia y la Nanotecnologı́a (CEDENNA), Avda. Ecuador 3493, Santiago 9170124, Chile
| |
Collapse
|
6
|
Schütt KT, Gastegger M, Tkatchenko A, Müller KR, Maurer RJ. Unifying machine learning and quantum chemistry with a deep neural network for molecular wavefunctions. Nat Commun 2019; 10:5024. [PMID: 31729373 PMCID: PMC6858523 DOI: 10.1038/s41467-019-12875-2] [Citation(s) in RCA: 189] [Impact Index Per Article: 37.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Accepted: 09/25/2019] [Indexed: 12/03/2022] Open
Abstract
Machine learning advances chemistry and materials science by enabling large-scale exploration of chemical space based on quantum chemical calculations. While these models supply fast and accurate predictions of atomistic chemical properties, they do not explicitly capture the electronic degrees of freedom of a molecule, which limits their applicability for reactive chemistry and chemical analysis. Here we present a deep learning framework for the prediction of the quantum mechanical wavefunction in a local basis of atomic orbitals from which all other ground-state properties can be derived. This approach retains full access to the electronic structure via the wavefunction at force-field-like efficiency and captures quantum mechanics in an analytically differentiable representation. On several examples, we demonstrate that this opens promising avenues to perform inverse design of molecular structures for targeting electronic property optimisation and a clear path towards increased synergy of machine learning and quantum chemistry.
Collapse
Affiliation(s)
- K T Schütt
- Machine Learning Group, Technische Universität Berlin, 10587, Berlin, Germany
| | - M Gastegger
- Machine Learning Group, Technische Universität Berlin, 10587, Berlin, Germany
| | - A Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511, Luxembourg, Luxembourg.
| | - K-R Müller
- Machine Learning Group, Technische Universität Berlin, 10587, Berlin, Germany.
- Department of Brain and Cognitive Engineering, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea.
- Max-Planck-Institut für Informatik, Saarbrücken, Germany.
| | - R J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Road, CV4 7AL, Coventry, UK.
| |
Collapse
|
7
|
von Rudorff GF, von Lilienfeld OA. Atoms in Molecules from Alchemical Perturbation Density Functional Theory. J Phys Chem B 2019; 123:10073-10082. [PMID: 31647233 DOI: 10.1021/acs.jpcb.9b07799] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Based on thermodynamic integration, we introduce atoms in molecules (AIM) using the orbital-free framework of alchemical perturbation density functional theory (APDFT). Within APDFT, atomic energies and electron densities in molecules are arbitrary because any reference system and integration path can be selected as long as it meets the boundary conditions. We choose the uniform electron gas (jellium) as a reference and linearly scale up all nuclear charges, situated at any query molecule's atomic coordinates. Within the approximations made when calculating one-particle electron densities, this universal choice affords unambiguous and exact definitions of energies and electron densities of AIMs. Numerical results are presented for neutral small molecules (CO, N2, BF, CO2), various small molecules with different electronic hybridization states of carbon (CH4, C2H6, C2H4, C2H2, HCN), and all of the possible BN-doped mutants connecting benzene to borazine (C2nB3-nN3-nH6, 0 ≤ n ≤ 3). Our results, as well as comparison to atomic energy estimates resulting from either DFT trained neural network models or atomic basis set overlap within CCSD, suggest that APDFT based AIMs enable meaningful, interesting, and counterintuitive interpretations of chemical bonding and molecular electron densities.
Collapse
Affiliation(s)
- Guido Falk von Rudorff
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry , University of Basel , Klingelbergstrasse 80 , CH-4056 Basel , Switzerland
| |
Collapse
|
8
|
Freeze JG, Kelly HR, Batista VS. Search for Catalysts by Inverse Design: Artificial Intelligence, Mountain Climbers, and Alchemists. Chem Rev 2019; 119:6595-6612. [PMID: 31059236 DOI: 10.1021/acs.chemrev.8b00759] [Citation(s) in RCA: 99] [Impact Index Per Article: 19.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In silico catalyst design is a grand challenge of chemistry. Traditional computational approaches have been limited by the need to compute properties for an intractably large number of possible catalysts. Recently, inverse design methods have emerged, starting from a desired property and optimizing a corresponding chemical structure. Techniques used for exploring chemical space include gradient-based optimization, alchemical transformations, and machine learning. Though the application of these methods to catalysis is in its early stages, further development will allow for robust computational catalyst design. This review provides an overview of the evolution of inverse design approaches and their relevance to catalysis. The strengths and limitations of existing techniques are highlighted, and suggestions for future research are provided.
Collapse
Affiliation(s)
- Jessica G Freeze
- Department of Chemistry , Yale University , New Haven , Connecticut 06520 , United States.,Energy Sciences Institute , Yale University , West Haven , Connecticut 06516 , United States
| | - H Ray Kelly
- Department of Chemistry , Yale University , New Haven , Connecticut 06520 , United States.,Energy Sciences Institute , Yale University , West Haven , Connecticut 06516 , United States
| | - Victor S Batista
- Energy Sciences Institute , Yale University , West Haven , Connecticut 06516 , United States.,Department of Chemistry , Yale University , P.O. Box 208107 , New Haven , Connecticut 06520 , United States
| |
Collapse
|
9
|
Abstract
In silico design of new molecules and materials with desirable quantum properties by high-throughput screening is a major challenge due to the high dimensionality of chemical space. To facilitate its navigation, we present a unification of coordinate and composition space in terms of alchemical normal modes (ANMs) which result from second order perturbation theory. ANMs assume a predominantly smooth nature of chemical space and form a basis in which new compounds can be expanded and identified. We showcase the use of ANMs for the energetics of the isoelectronic series of diatomics with 14 electrons, BN doped benzene derivatives (C6-2 x(BN) xH6 with x = 0,1,2,3), predictions for over 1.8 million BN doped coronene derivatives, and genetic energy optimizations in the entire BN-doped coronene space. Using Ge lattice scans as reference, the applicability of ANMs across the periodic table is demonstrated for III-V and IV-IV semiconductors Si, Sn, SiGe, SnGe, SiSn, as well as AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs, and InSb. Analysis of our results indicates simple qualitative structure property rules for estimating energetic rankings among isomers. Useful quantitative estimates can also be obtained when few atoms are changed to neighboring or lower lying elements in the periodic table. The quality of the predictions often increases with the symmetry of system chosen as reference due to cancellation of odd order terms. Rooted in perturbation theory, the ANM approach promises to generally enable unbiased compound exploration campaigns at reduced computational cost.
Collapse
Affiliation(s)
- Stijn Fias
- General Chemistry (ALGC) , Vrije Universiteit Brussel (Free University Brussels - VUB) , Pleinlaan 2 , 1050 Brussel , Belgium
- Department of Chemistry & Chemical Biology , McMaster University , Hamilton , ON , Canada L8S 4L8
| | - K Y Samuel Chang
- 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
| | - O Anatole von Lilienfeld
- 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
|
10
|
Griego CD, Saravanan K, Keith JA. Benchmarking Computational Alchemy for Carbide, Nitride, and Oxide Catalysts. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800142] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Charles D. Griego
- Department of Chemical and Petroleum Engineering University of Pittsburgh 3700 O'Hara Street Pittsburgh PA 15261 USA
| | - Karthikeyan Saravanan
- Department of Chemical and Petroleum Engineering University of Pittsburgh 3700 O'Hara Street Pittsburgh PA 15261 USA
| | - John A. Keith
- Department of Chemical and Petroleum Engineering University of Pittsburgh 3700 O'Hara Street Pittsburgh PA 15261 USA
| |
Collapse
|
11
|
Al-Hamdani YS, Michaelides A, von Lilienfeld OA. Exploring dissociative water adsorption on isoelectronically BN doped graphene using alchemical derivatives. J Chem Phys 2018; 147:164113. [PMID: 29096500 DOI: 10.1063/1.4986314] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The design and production of novel 2-dimensional materials have seen great progress in the last decade, prompting further exploration of the chemistry of such materials. Doping and hydrogenating graphene are an experimentally realised method of changing its surface chemistry, but there is still a great deal to be understood on how doping impacts on the adsorption of molecules. Developing this understanding is key to unlocking the potential applications of these materials. High throughput screening methods can provide particularly effective ways to explore vast chemical compositions of materials. Here, alchemical derivatives are used as a method to screen the dissociative adsorption energy of water molecules on various BN doped topologies of hydrogenated graphene. The predictions from alchemical derivatives are assessed by comparison to density functional theory. This screening method is found to predict dissociative adsorption energies that span a range of more than 2 eV, with a mean absolute error <0.1 eV. In addition, we show that the quality of such predictions can be readily assessed by examination of the Kohn-Sham highest occupied molecular orbital in the initial states. In this way, the root mean square error in the dissociative adsorption energies of water is reduced by almost an order of magnitude (down to ∼0.02 eV) after filtering out poor predictions. The findings point the way towards a reliable use of first order alchemical derivatives for efficient screening procedures.
Collapse
Affiliation(s)
- Yasmine S Al-Hamdani
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
12
|
Abstract
The chemical space contains all possible compounds that can be imagined. Its size easily equals the number of fundamental particles in the observable universe. Rational design of compounds aims to find those sectors of the chemical space where compounds optimize a set of desired properties. Then, rational design demands tools to efficiently navigate the chemical space. Ab initio alchemical derivatives offer the possibility to navigate, without empiricism, the energy landscape through alchemical transformations. An alchemical transformation is any process, physical or fictitious, that connects to points in the chemical space. In this work, those transformations are constructed as a perturbative expansion of the energy with respect to perturbations in the stoichiometry. The response functions of that expansion are what is called alchemical derivatives. In this work we assess how effective alchemical derivatives are in predicting energy changes associated to changes in the composition. We do this by including in the expansion, for the first time, electrostatic, polarization and electron-transfer effects. The system we chose is one that challenges alchemical derivatives because none of these effects dominates its behavior. The transmutations studied here correspond to substitutional doping of Al13 with up to four atoms of Si, Al13-nSin. Two types of transformations are considered, those in which the number of electrons remains constant and those in which the number of electrons also changes. It is found that contrary to what has been reported before, polarization cannot be neglected. If polarization is not included, alchemical derivatives fail to predict the change of energy and the relative energy between isomers. For isoelectronic substitution of four or more atoms, the perturbative approach collapses because the strength of the perturbation becomes too strong to guarantee convergence of the series. It is shown, however, that if only one atom is mutated at a time, alchemical derivatives rank pretty well the isomers of Al13-nSin according to their energy. In the case of non-isoelectronic transformations, it is observed that the series rapidly diverges with increasing number of electrons. In this situation, it becomes more important to keep the degree of transmutation of the parent system small.
Collapse
Affiliation(s)
- Macarena Muñoz
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, 653-Santiago, Chile.
| | | |
Collapse
|
13
|
Balawender R, Lesiuk M, De Proft F, Geerlings P. Exploring Chemical Space with Alchemical Derivatives: BN-Simultaneous Substitution Patterns in C 60. J Chem Theory Comput 2018; 14:1154-1168. [PMID: 29300479 DOI: 10.1021/acs.jctc.7b01114] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
With the idea of using alchemical derivatives to explore in an efficient, computer- and cost-effective way Chemical Space was launched several years ago. In the context of Conceptual DFT response functions, these energies vs nuclear charge derivatives permit the estimatation of the energy of transmutants of a given starting or reference molecule showing different nuclear compositions. After an explorative study on small and planar molecules ( Balawender et al. J. Chem. Theory Comput. 2013 , 9 , 5327 ) by the present authors of this paper, the present study fully exploits the computational advantages of the alchemical derivatives in larger three-dimensional systems. Starting from a single reference calculation on C60, the complete BN substitution pattern, from single substituted C58BN via the belt (C20(BN)20 and the ball C12(BN)24 structures to the fully substituted (BN)30, is explored. Successive and simultaneous substitution strategies are followed and compared, indicating that both techniques yield identical results up to 13 substitutions but that for higher substitutions the simultaneous approach needs to be taken. Due to the cost-efficiency of the algorithm this path can indeed be followed as opposed to earlier work in the literature where for each step a full SCF calculation was at stake leading to prohibitively large computational demands for adopting the simultaneous approach. Previously formulated rules governing the substitution pattern by Kar and co-workers are scrutinized in this context and reformulated giving chemical insight in the gradual substitution process and the relative energies of the isomers. In its present form the method offers an interesting venue to study BN substitution patterns in higher fullerenes and graphene and in general paves the way for more efficient exploration of the Chemical Space.
Collapse
Affiliation(s)
- Robert Balawender
- Institute of Physical Chemistry, Polish Academy of Sciences , Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Michał Lesiuk
- Faculty of Chemistry, University of Warsaw , Pasteura 1, PL-02-093 Warsaw, Poland
| | - Frank De Proft
- Algemene Chemie, Vrije Universiteit Brussel, Faculteit Wetenschappen , Pleinlaan 2, 1050 Brussels, Belgium
| | - Paul Geerlings
- Algemene Chemie, Vrije Universiteit Brussel, Faculteit Wetenschappen , Pleinlaan 2, 1050 Brussels, Belgium
| |
Collapse
|
14
|
Saravanan K, Kitchin JR, von Lilienfeld OA, Keith JA. Alchemical Predictions for Computational Catalysis: Potential and Limitations. J Phys Chem Lett 2017; 8:5002-5007. [PMID: 28938798 DOI: 10.1021/acs.jpclett.7b01974] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Kohn-Sham density functional theory (DFT) is the workhorse method for calculating adsorbate binding energies relevant for catalysis. Unfortunately, this method is too computationally expensive to methodically and broadly search through catalyst candidate space. Here, we assess the promise of computational alchemy, a perturbation theory approach that allows for predictions of binding energies thousands of times faster than DFT. We first benchmark the binding energy predictions of oxygen reduction reaction intermediates on alloys of Pt, Pd, and Ni using alchemy against predictions from DFT. Far faster alchemical estimates yield binding energies within 0.1 eV of DFT values in many cases. We also identify distinct cases where alchemy performs significantly worse, indicating areas where modeling improvements are needed. Our results suggest that computational alchemy is a very promising tool that warrants further consideration for high-throughput screening of heterogeneous catalysts.
Collapse
Affiliation(s)
- Karthikeyan Saravanan
- Department of Chemical and Petroleum Engineering, Swanson School of Engineering, University of Pittsburgh , Pittsburgh, Pennsylvania 15261, United States
| | - John R Kitchin
- Department of Chemical Engineering, Carnegie Mellon University , Pittsburgh, Pennsylvania 15213, United States
| | - O Anatole von Lilienfeld
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials, Department of Chemistry, University of Basel , 4001 Basel, Switzerland
| | - John A Keith
- Department of Chemical and Petroleum Engineering, Swanson School of Engineering, University of Pittsburgh , Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
15
|
López-Rosa S, Molina-Espíritu M, Esquivel RO, Soriano-Correa C, Dehesa JS. Study of the Chemical Space of Selected Bacteriostatic Sulfonamides from an Information Theory Point of View. Chemphyschem 2016; 17:4003-4010. [DOI: 10.1002/cphc.201600790] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2016] [Indexed: 11/08/2022]
Affiliation(s)
- Sheila López-Rosa
- Departamento de Fsica Aplicada II; Universidad de Sevilla; 41012 Sevilla Spain
- Instituto Carlos I de Fsica Teórica y Computacional; Universidad de Granada; 18071 Granada Spain
| | | | - Rodolfo O. Esquivel
- Instituto Carlos I de Fsica Teórica y Computacional; Universidad de Granada; 18071 Granada Spain
- Departamento de Química; Universidad Autónoma Metropolitana; 09340- Mexico City México
| | - Catalina Soriano-Correa
- Química Computacional, FES-Zaragoza; Universidad Nacional Autónoma de México; 09230- Iztapalapa, Mexico City México
| | - Jésus S. Dehesa
- Instituto Carlos I de Fsica Teórica y Computacional; Universidad de Granada; 18071 Granada Spain
- Departamento de Fsica Atómica, Molecular y Nuclear; Universidad de Granada; 18071 Granada Spain
| |
Collapse
|
16
|
Solovyeva A, von Lilienfeld OA. Alchemical screening of ionic crystals. Phys Chem Chem Phys 2016; 18:31078-31091. [DOI: 10.1039/c6cp04258a] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We introduce alchemical perturbations as a rapid and accurate tool to estimate fundamental structural and energetic properties in pure and mixed ionic crystals.
Collapse
Affiliation(s)
- Alisa Solovyeva
- 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
| | - O. Anatole von Lilienfeld
- 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
|