PASTA: Python Algorithms for Searching Transition stAtes

Chemical reactions are often associated with an energy barrier along the reaction pathway which hinders the spontaneity of the reaction. Changing the energy barrier along the reaction pathway allows one to modulate the performance of a reaction. We present a module, Python Algorithms for Searching Transition stAtes (PASTA), to calculate the energy barrier and locate the transition state of a reaction efficiently. The module is written in python and can perform nudged elastic band, climbing image nudged elastic band and automated nudged elastic band calculations. These methods require the knowledge of the potential energy surface (and its gradient along some direction). This module is written such that it works in conjunction with density functional theory (OFT) codes to obtain this information. Presently it is interfaced with three well known OFT packages: Vienna Ab initio Simulation Package (VASP), Quantum Espresso and Spanish Initiative for Electronic Simulations with Thousands of Atoms (SIESTA). This module is easily extendable and can be interfaced with other OFT, force-field or empirical potential based codes. The uniqueness of the module lies in its user-friendliness. For users with limited computing resources, this module will be an effective tool as it allows to perform the calculations image by image. On the other hand, users with plentiful computing resources (such as users in a high performance computing environment) can perform the calculations for large number of images simultaneously. This module gives users complete flexibility, thereby enabling them to perform calculations on large systems making the best use of the available resources. Program summary Program Title: PASTA Program files doi : http://dx.doi.org/10.17632/rv7fdm5gkf.1 Licensing provisions: BSD 3-clause Programming language: Python External routines/libraries: numpy, matplotlib Nature of problem: Most of the reactions have an energy barrier on their reaction pathway. This energy barrier affects the progress of the reaction. Solution method: We implement the NEB, CI-NEB and AutoNEB method to locate transition state and estimate the energy barrier. Our module works with density functional theory codes: VASP, SIESTA and Quantum Espresso presently. (C) 2018 Elsevier B.V. All rights reserved.

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

[2]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[3]  H. Jónsson,et al.  Nudged elastic band method for finding minimum energy paths of transitions , 1998 .

[4]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

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

[6]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[7]  Hafner,et al.  Ab initio molecular dynamics for liquid metals. , 1995, Physical review. B, Condensed matter.

[8]  Emily A. Carter,et al.  Adsorption and dissociation of CO on Fe(1 1 0) from first principles , 2004 .

[9]  Michael Walter,et al.  The atomic simulation environment-a Python library for working with atoms. , 2017, Journal of physics. Condensed matter : an Institute of Physics journal.

[10]  Satoshi Itoh,et al.  Calculations of Silicon Self-Interstitial Defects , 1999 .

[11]  Lindsey J. Munro,et al.  DEFECT MIGRATION IN CRYSTALLINE SILICON , 1999 .

[12]  G. Scuseria,et al.  Restoring the density-gradient expansion for exchange in solids and surfaces. , 2007, Physical review letters.

[13]  Rabe,et al.  Optimized pseudopotentials. , 1990, Physical review. B, Condensed matter.

[14]  Anton Kokalj,et al.  Computer graphics and graphical user interfaces as tools in simulations of matter at the atomic scale , 2003 .

[15]  G. Henkelman,et al.  Methods for Finding Saddle Points and Minimum Energy Paths , 2002 .

[16]  Mousseau,et al.  Dynamics of lennard-jones clusters: A characterization of the activation-relaxation technique , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  G. vanRossum.  Python reference manual , 1995 .

[18]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[19]  Stefano de Gironcoli,et al.  QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[21]  Richard J. Needs.  First-principles calculations of self-interstitial defect structures and diffusion paths in silicon , 1999 .

[22]  Soler,et al.  Self-consistent order-N density-functional calculations for very large systems. , 1996, Physical review. B, Condensed matter.

[23]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[24]  Egill Skúlason,et al.  A systematic DFT study of hydrogen diffusion on transition metal surfaces , 2012 .

[25]  Stefano de Gironcoli,et al.  Reproducibility in density functional theory calculations of solids , 2016, Science.

[26]  Graeme Henkelman,et al.  A generalized solid-state nudged elastic band method. , 2012, The Journal of chemical physics.

[27]  E. L. Kolsbjerg,et al.  An automated nudged elastic band method. , 2016, The Journal of chemical physics.

[28]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[29]  Hafner,et al.  Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. , 1994, Physical review. B, Condensed matter.

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

[31]  G. Henkelman,et al.  Optimization methods for finding minimum energy paths. , 2008, The Journal of chemical physics.