EON: software for long time simulations of atomic scale systems

被引:58
作者
Chill, Samuel T. [1 ,2 ]
Welborn, Matthew [1 ,2 ]
Terrell, Rye [1 ,2 ]
Zhang, Liang [1 ,2 ]
Berthet, Jean-Claude [3 ,4 ]
Pedersen, Andreas [3 ,4 ]
Jonsson, Hannes [3 ,4 ]
Henkelman, Graeme [1 ,2 ]
机构
[1] Univ Texas Austin, Dept Chem, Austin, TX 78712 USA
[2] Univ Texas Austin, Inst Computat Engn & Sci, Austin, TX 78712 USA
[3] Univ Iceland, Fac Phys Sci, IS-107 Reykjavik, Iceland
[4] Univ Iceland, Inst Sci, IS-107 Reykjavik, Iceland
基金
美国国家科学基金会;
关键词
rare event dynamics; transition state finding; saddle points; adaptive kinetic Monte Carlo; accelerated molecular dynamics; ACCELERATED MOLECULAR-DYNAMICS; LENNARD-JONES CLUSTERS; GLOBAL OPTIMIZATION; ALGORITHMS; SURFACES; DEFECTS; SILICON; MODEL;
D O I
10.1088/0965-0393/22/5/055002
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
The EON software is designed for simulations of the state-to-state evolution of atomic scale systems over timescales greatly exceeding that of direct classical dynamics. States are defined as collections of atomic configurations from which a minimization of the potential energy gives the same inherent structure. The time evolution is assumed to be governed by rare events, where transitions between states are uncorrelated and infrequent compared with the timescale of atomic vibrations. Several methods for calculating the state-to-state evolution have been implemented in EON, including parallel replica dynamics, hyperdynamics and adaptive kinetic Monte Carlo. Global optimization methods, including simulated annealing, basin hopping and minima hopping are also implemented. The software has a client/server architecture where the computationally intensive evaluations of the interatomic interactions are calculated on the client-side and the state-to-state evolution is managed by the server. The client supports optimization for different computer architectures to maximize computational efficiency. The server is written in Python so that developers have access to the high-level functionality without delving into the computationally intensive components. Communication between the server and clients is abstracted so that calculations can be deployed on a single machine, clusters using a queuing system, large parallel computers using a message passing interface, or within a distributed computing environment. A generic interface to the evaluation of the interatomic interactions is defined so that empirical potentials, such as in LAMMPS, and density functional theory as implemented in VASP and GPAW can be used interchangeably. Examples are given to demonstrate the range of systems that can be modeled, including surface diffusion and island ripening of adsorbed atoms on metal surfaces, molecular diffusion on the surface of ice and global structural optimization of nanoparticles.
引用
收藏
页数:16
相关论文
共 41 条
[1]  
Allen M. P., 1987, COMPUTER SIMULATION
[2]   BOINC: A system for public-resource computing and storage [J].
Anderson, DP .
FIFTH IEEE/ACM INTERNATIONAL WORKSHOP ON GRID COMPUTING, PROCEEDINGS, 2004, :4-10
[3]  
[Anonymous], 2003, J THEORET BIOL
[4]   Structural relaxation made simple [J].
Bitzek, Erik ;
Koskinen, Pekka ;
Gaehler, Franz ;
Moseler, Michael ;
Gumbsch, Peter .
PHYSICAL REVIEW LETTERS, 2006, 97 (17)
[5]   EMBEDDED-ATOM METHOD - DERIVATION AND APPLICATION TO IMPURITIES, SURFACES, AND OTHER DEFECTS IN METALS [J].
DAW, MS ;
BASKES, MI .
PHYSICAL REVIEW B, 1984, 29 (12) :6443-6453
[7]   A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives [J].
Henkelman, G ;
Jónsson, H .
JOURNAL OF CHEMICAL PHYSICS, 1999, 111 (15) :7010-7022
[8]   Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table [J].
Henkelman, G ;
Jónsson, H .
JOURNAL OF CHEMICAL PHYSICS, 2001, 115 (21) :9657-9666
[9]   METHODS OF CONJUGATE GRADIENTS FOR SOLVING LINEAR SYSTEMS [J].
HESTENES, MR ;
STIEFEL, E .
JOURNAL OF RESEARCH OF THE NATIONAL BUREAU OF STANDARDS, 1952, 49 (06) :409-436
[10]  
Horn R.A., 2012, Matrix Analysis