pyiron: An integrated development environment for computational materials science

Abstract To support and accelerate the development of simulation protocols in atomistic modelling, we introduce an integrated development environment (IDE) for computational materials science called pyiron ( http://pyiron.org ). The pyiron IDE combines a web based source code editor, a job management system for build automation, and a hierarchical data management solution. The core components of the pyiron IDE are pyiron objects based on an abstract class, which links application structures such as atomistic structures, projects, jobs, simulation protocols and computing resources with persistent storage and an interactive user environment. The simulation protocols within the pyiron IDE are constructed using the Python programming language. To highlight key concepts of this tool as well as to demonstrate its ability to simplify the implementation and testing of simulation protocols we discuss two applications. In these examples we show how pyiron supports the whole life cycle of a typical simulation, seamlessly combines ab initio with empirical potential calculations, and how complex feedback loops can be implemented. While originally developed with focus on ab initio thermodynamics simulations, the concepts and implementation of pyiron are general thus allowing to employ it for a wide range of simulation topics.

[1]  Blazej Grabowski,et al.  Improved method of calculating ab initio high-temperature thermodynamic properties with application to ZrC , 2015 .

[2]  Blazej Grabowski,et al.  Efficient approach to compute melting properties fully from ab initio with application to Cu , 2017 .

[3]  Ronald Morrison,et al.  Orthogonally persistent object systems , 1995, The VLDB Journal.

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

[5]  F. Murnaghan The Compressibility of Media under Extreme Pressures. , 1944, Proceedings of the National Academy of Sciences of the United States of America.

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

[7]  Blazej Grabowski,et al.  Temperature-driven phase transitions from first principles including all relevant excitations: The fcc-to-bcc transition in Ca , 2011 .

[8]  Neville Reid Moody,et al.  COMMENT: Trapping of hydrogen to lattice defects in nickel , 1995 .

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

[10]  Chandler A. Becker,et al.  Evaluating variability with atomistic simulations: the effect of potential and calculation methodology on the modeling of lattice and elastic constants , 2018 .

[11]  B. Grabowski,et al.  Advancing density functional theory to finite temperatures: methods and applications in steel design , 2012, Journal of physics. Condensed matter : an Institute of Physics journal.

[12]  Kamal Choudhary,et al.  High-throughput Identification and Characterization of Two-dimensional Materials using Density functional theory , 2017, Scientific Reports.

[13]  Matthew Horton,et al.  Atomate: A high-level interface to generate, execute, and analyze computational materials science workflows , 2017 .

[14]  Kamal Choudhary,et al.  Evaluation and comparison of classical interatomic potentials through a user-friendly interactive web-interface , 2017, Scientific Data.

[15]  B. Grabowski,et al.  Understanding Anharmonicity in fcc Materials: From its Origin to ab initio Strategies beyond the Quasiharmonic Approximation. , 2015, Physical review letters.

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

[17]  Roger C. Reed,et al.  Temperature dependence of the Gibbs energy of vacancy formation of fcc Ni , 2018, Physical Review B.

[18]  James P. Sethna,et al.  The potential of atomistic simulations and the knowledgebase of interatomic models , 2011 .

[19]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

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

[21]  Boris Kozinsky,et al.  AiiDA: Automated Interactive Infrastructure and Database for Computational Science , 2015, ArXiv.

[22]  Brian E. Granger,et al.  IPython: A System for Interactive Scientific Computing , 2007, Computing in Science & Engineering.

[23]  Henry Wu,et al.  The MAterials Simulation Toolkit (MAST) for atomistic modeling of defects and diffusion , 2017 .

[24]  Albert V. Davydov,et al.  MPInterfaces: A Materials Project based Python tool for high-throughput computational screening of interfacial systems , 2016, 1602.07784.

[25]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[26]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[27]  Vladan Stevanovic,et al.  A Computational Framework for Automation of Point Defect Calculations , 2016, 1611.00825.

[28]  S. Curtarolo,et al.  AFLOW: An automatic framework for high-throughput materials discovery , 2012, 1308.5715.

[29]  Anubhav Jain,et al.  A high-throughput infrastructure for density functional theory calculations , 2011 .

[30]  Blazej Grabowski,et al.  Ab initio up to the melting point: Anharmonicity and vacancies in aluminum , 2009 .

[31]  Pierre de Buyl,et al.  H5MD: A structured, efficient, and portable file format for molecular data , 2013, Comput. Phys. Commun..

[32]  G. Kresse,et al.  Ab initio molecular dynamics for liquid metals. , 1993 .

[33]  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.

[34]  Alexander S. Rose,et al.  NGLview–interactive molecular graphics for Jupyter notebooks , 2018, Bioinform..