A high-throughput framework for determining adsorption energies on solid surfaces

In this work, we present a high-throughput workflow for calculation of adsorption energies on solid surfaces using density functional theory. Using open-source computational tools from the Materials Project infrastructure, we automate the procedure of constructing symmetrically distinct adsorbate configurations for arbitrary slabs. These algorithms are further used to construct and run workflows in a standard, automated way such that user intervention in the simulation procedure is minimal. To validate our approach, we compare results from our workflow to previous experimental and theoretical benchmarks from the CE27 database of chemisorption energies on solid surfaces. These benchmarks also illustrate how the task of performing and managing over 200 individual density functional theory calculations may be reduced to a single submission procedure and subsequent analysis. By enabling more efficient high-throughput computations of adsorption energies, these tools will accelerate theory-guided discovery of advanced materials for applications in catalysis and surface science.Surface chemistry: an automatic sense of attractionAn automated procedure for determining the energy required for a molecule to adhere to a surface is developed by researchers in the United States. Joseph Montoya from the Lawrence Berkeley National Laboratory and Kristin Persson from the University of California, Berkeley, introduce an algorithm for finding the adsorption sites on an arbitrary surface. Knowing the amount of energy required for molecular adsorption is crucial for identifying the best materials for use in electronics and catalysis. Density functional theory can predict adsorption energies but usually requires human intuition to tune the calculations. With so many combinations of surface and adsorbate, an automated method is required. Montoya and Persson use open-source computational tools from the Materials Project to present a workflow for performing high-throughput density functional theory calculations for arbitrary slabs and adsorbed species.

[1]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[2]  J. Nørskov,et al.  Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals , 1999 .

[3]  M. Hove,et al.  Adsorption of CO on Pd(100) , 1980 .

[4]  Wei Chen,et al.  FireWorks: a dynamic workflow system designed for high‐throughput applications , 2015, Concurr. Comput. Pract. Exp..

[5]  Klaus Christmann,et al.  Interaction of hydrogen with solid surfaces , 1988 .

[6]  Kristin A. Persson,et al.  Commentary: The Materials Project: A materials genome approach to accelerating materials innovation , 2013 .

[7]  Andrew A. Peterson,et al.  Global Optimization of Adsorbate–Surface Structures While Preserving Molecular Identity , 2014, Topics in Catalysis.

[8]  I. Chorkendorff,et al.  Biomimetic Hydrogen Evolution: MoS2 Nanoparticles as Catalyst for Hydrogen Evolution , 2005 .

[9]  W. A. Brown,et al.  Femtomole Adsorption Calorimetry on Single-Crystal Surfaces. , 1998, Chemical reviews.

[10]  Anubhav Jain,et al.  Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis , 2012 .

[11]  Cormac Toher,et al.  Charting the complete elastic properties of inorganic crystalline compounds , 2015, Scientific Data.

[12]  Ib Chorkendorff,et al.  Discovery of a Ni-Ga catalyst for carbon dioxide reduction to methanol. , 2014, Nature chemistry.

[13]  Roman Garnett,et al.  BASC: Applying Bayesian Optimization to the Search for Global Minima on Potential Energy Surfaces , 2016, ICML.

[14]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[15]  G. Ertl,et al.  Adsorption of CO on Pd single crystal surfaces , 1974 .

[16]  M. Grunze,et al.  Interaction of nitrogen with iron surfaces: I. Fe(100) and Fe(111) , 1977 .

[17]  Thomas Bligaard,et al.  Density functionals for surface science: Exchange-correlation model development with Bayesian error estimation , 2012 .

[18]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[19]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[20]  J. Nørskov,et al.  Towards the computational design of solid catalysts. , 2009, Nature chemistry.

[21]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[22]  K. Jacobsen,et al.  Real-space grid implementation of the projector augmented wave method , 2004, cond-mat/0411218.

[23]  D. King,et al.  Adsorption energetics and bonding from femtomole calorimetry and from first principles theory , 2000 .

[24]  Zachary W. Ulissi,et al.  Automated Discovery and Construction of Surface Phase Diagrams Using Machine Learning. , 2016, The journal of physical chemistry letters.

[25]  D. King,et al.  Calorimetric investigation of NO and CO adsorption on Pd{100} and the influence of preadsorbed carbon , 1997 .

[26]  V. Fiorin,et al.  Microcalorimetry of O2 and NO on flat and stepped platinum surfaces , 2009 .

[27]  Kristin A. Persson,et al.  Surface energies of elemental crystals , 2016, Scientific Data.

[28]  F. Abild‐Pedersen,et al.  CO adsorption energies on metals with correction for high coordination adsorption sites – A density functional study , 2007 .

[29]  N. A. Romero,et al.  Electronic structure calculations with GPAW: a real-space implementation of the projector augmented-wave method , 2010, Journal of physics. Condensed matter : an Institute of Physics journal.

[30]  Philippe Sautet,et al.  Introducing structural sensitivity into adsorption-energy scaling relations by means of coordination numbers. , 2015, Nature chemistry.