The MOLDY short-range molecular dynamics package

Abstract We describe a parallelised version of the MOLDY molecular dynamics program. This Fortran code is aimed at systems which may be described by short-range potentials and specifically those which may be addressed with the embedded atom method. This includes a wide range of transition metals and alloys. MOLDY provides a range of options in terms of the molecular dynamics ensemble used and the boundary conditions which may be applied. A number of standard potentials are provided, and the modular structure of the code allows new potentials to be added easily. The code is parallelised using OpenMP and can therefore be run on shared memory systems, including modern multicore processors. Particular attention is paid to the updates required in the main force loop, where synchronisation is often required in OpenMP implementations of molecular dynamics. We examine the performance of the parallel code in detail and give some examples of applications to realistic problems, including the dynamic compression of copper and carbon migration in an iron–carbon alloy. Program summary Program title: MOLDY Catalogue identifier: AEJU_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEJU_v1_0.html Program obtainable from: CPC Program Library, Queenʼs University, Belfast, N. Ireland Licensing provisions: GNU General Public License version 2 No. of lines in distributed program, including test data, etc.: 382 881 No. of bytes in distributed program, including test data, etc.: 6 705 242 Distribution format: tar.gz Programming language: Fortran 95/OpenMP Computer: Any Operating system: Any Has the code been vectorised or parallelized?: Yes. OpenMP is required for parallel execution RAM: 100 MB or more Classification: 7.7 Nature of problem: Moldy addresses the problem of many atoms (of order 10 6 ) interacting via a classical interatomic potential on a timescale of microseconds. It is designed for problems where statistics must be gathered over a number of equivalent runs, such as measuring thermodynamic properities, diffusion, radiation damage, fracture, twinning deformation, nucleation and growth of phase transitions, sputtering etc. In the vast majority of materials, the interactions are non-pairwise, and the code must be able to deal with many-body forces. Solution method: Molecular dynamics involves integrating Newtonʼs equations of motion. MOLDY uses verlet (for good energy conservation) or predictor–corrector (for accurate trajectories) algorithms. It is parallelised using open MP. It also includes a static minimisation routine to find the lowest energy structure. Boundary conditions for surfaces, clusters, grain boundaries, thermostat (Nose), barostat (Parrinello–Rahman), and externally applied strain are provided. The initial configuration can be either a repeated unit cell or have all atoms given explictly. Initial velocities are generated internally, but it is also possible to specify the velocity of a particular atom. A wide range of interatomic force models are implemented, including embedded atom, Morse or Lennard-Jones. Thus the program is especially well suited to calculations of metals. Restrictions: The code is designed for short-ranged potentials, and there is no Ewald sum. Thus for long range interactions where all particles interact with all others, the order- N scaling will fail. Different interatomic potential forms require recompilation of the code. Additional comments: There is a set of associated open-source analysis software for postprocessing and visualisation. This includes local crystal structure recognition and identification of topological defects. Running time: A set of test modules for running time are provided. The code scales as order N . The parallelisation shows near-linear scaling with number of processors in a shared memory environment. A typical run of a few tens of nanometers for a few nanoseconds will run on a timescale of days on a multiprocessor desktop.

[1]  G. Ackland,et al.  An improved N-body semi-empirical model for body-centred cubic transition metals , 1987 .

[2]  Keith Refson,et al.  Moldy: a portable molecular dynamics simulation program for serial and parallel computers , 2000 .

[3]  Y. P. Varshni.  Temperature Dependence of the Elastic Constants , 1970 .

[4]  R. Mclellan,et al.  Carbon-vacancy interactions in B.C.C. iron , 1988 .

[5]  G. Ackland,et al.  Development of an interatomic potential for the simulation of phase transformations in zirconium , 2007 .

[6]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[7]  A. Serra,et al.  Study of copper precipitates in α‐iron by computer simulation I. Interatomic potentials and properties of Fe and Cu , 1995 .

[8]  F. Claisse,et al.  Improvement of Kryukov and Zhukhovitskii method for measurement of diffusion. Application to self-diffusion in alpha iron , 1968 .

[9]  G. Ackland,et al.  Rescaled potentials for transition metal solutes in α-iron , 2009 .

[10]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[11]  Molecular dynamics study of self-diffusion in bcc Fe , 2009, 0902.1338.

[12]  Two-band second moment model for transition metals and alloys , 2005, cond-mat/0505060.

[13]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[14]  V. Heine,et al.  Many-atom interactions in solids , 1991, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[15]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[16]  H. Schaefer,et al.  Vacancy formation in iron investigated by positron annihilation in thermal equilibrium , 1977 .

[17]  D. Bacon,et al.  Point-defect and stacking-fault properties in body-centred-cubic metals with n-body interatomic potentials , 1986 .

[18]  R. Wentzcovitch,et al.  Invariant molecular-dynamics approach to structural phase transitions. , 1991, Physical review. B, Condensed matter.

[19]  Peter M. Derlet,et al.  A ‘magnetic’ interatomic potential for molecular dynamics simulations , 2005 .

[20]  D. Bacon,et al.  A molecular dynamics study of displacement cascades in α-iron , 1993 .

[21]  D. McDowell,et al.  Atomistic simulations of homogeneous dislocation nucleation in single crystal copper , 2007 .

[22]  E. Hairer,et al.  Geometric numerical integration illustrated by the Störmer–Verlet method , 2003, Acta Numerica.

[23]  M. Baskes,et al.  Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals , 1984 .

[24]  G. Ackland,et al.  Metallic-covalent interatomic potential for carbon in iron , 2008 .

[25]  M. Finnis,et al.  A simple empirical N-body potential for transition metals , 1984 .

[26]  Graeme Ackland,et al.  Development of an interatomic potential for phosphorus impurities in α-iron , 2004 .

[27]  Graeme Ackland,et al.  Applications of local crystal structure measures in experiment and simulation , 2006 .

[28]  L. Schepper,et al.  Positron annihilation on pure and carbon-dopedα-iron in thermal equilibrium , 1983 .

[29]  G. Ackland,et al.  Validity of the second moment tight-binding model , 1988 .

[30]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[31]  S. Rice,et al.  ADVANCES IN CHEMICAL PHYSICS , 2002 .

[32]  S. Nosé,et al.  Constant pressure molecular dynamics for molecular systems , 1983 .

[33]  C. Domain,et al.  Ab initio calculations of defects in Fe and dilute Fe-Cu alloys , 2001 .

[34]  A. Serra,et al.  Study of copper precipitates in α-iron by computer simulation , 1997 .

[35]  Chu-Chun Fu,et al.  Stability and mobility of mono- and di-interstitials in alpha-Fe. , 2004, Physical review letters.

[36]  H. Ledbetter,et al.  Elastic constants of monocrystal iron from 3 to 500 K , 2006 .

[37]  C. L. Cleveland.  New equations of motion for molecular dynamics systems that change shape , 1988 .

[38]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[39]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[40]  A. Vehanen,et al.  Vacancies and carbon impurities inα- iron: Electron irradiation , 1982 .

[41]  G. Ackland,et al.  Simple N-body potentials for the noble metals and nickel , 1987 .