Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aewy_v1_0.tar.gz(55081 Kbytes)
Manuscript Title: MEAMfit: A reference-free modified embedded atom method (RF-MEAM) energy and force-fitting code
Authors: Andrew Ian Duff, M. W. Finnis, Philippe Maugis, Barend. J. Thijsse, Marcel H. F. Sluiter
Program title: MEAMfit
Catalogue identifier: AEWY_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 196(2015)439
Programming language: Fortran.
Computer: Linux based workstations.
Operating system: Linux.
RAM: 120 Megabytes
Keywords: Interatomic potential, Fitting, Lattice-vibrations.
PACS: 07.05.Tp.
Classification: 16.1.

External routines: TOMS611 Unconstrained Minimization [1] included in the MEAMfit code.

Nature of problem:
Fitting embedded atom method (EAM) and reference-free modified embedded atom method (RF-MEAM) potentials [2-3] to energies and forces produced by VASP [4-7].

Solution method:
A computer program is presented which uses a conjugate-gradient minimizer paired with a genetic algorithm to fit EAM and RF-MEAM potentials to energies and/or atomic forces read directly from VASP output files. Potentials produced by the code are directly usable with the LAMMPS [8] or Camelion [9] molecular-dynamics packages.

Additional comments:
User manual provided.

Running time:
The run-time depends on the required level of accuracy of the final potential. For an EAM potential fit to 670 energies, a few hours on a single core is usually sufficient to produce a potential with R = 12%-13% (see Equation. 9 in main-text for definition). To ensure a maximally optimized potential however (R = 12%), a run-time of 24 hours is recommended. To optimize a RF- MEAM potential, a further 24 hours should be allowed. One will already find an improvement over the EAM using just a single core, however to ensure a maximally optimized potential, one should run several instances of MEAMfit in parallel.

[1] J. E. Dennis, D. Gay and R. E. Welsch, ACM Trans. on Math. Soft., 7 (1981) 348-368.
[2] M. I. Baskes, Materials Science and Engineering A, 261 (1999), 165.
[3] M. Timonova and B. J. Thijsse, Modelling Simul. Mater. Sci. Eng., 19 (2011) 015003
[4] G. Kresse, J. Hafner, Phys. Rev. B 47 (1993) 558.
[5] G. Kresse, J. Hafner, Phys. Rev. B 49 (1994) 14251.
[6] G. Kresse, J. Furthmüller, Comput. Mat. Sci. 6 (1996) 15.
[7] G. Kresse, J. Furthmüller, Phys. Rev. B 54 (1996) 11169.
[8] S. Plimpton, J. Comp. Phys. 117 (1995) 1-19.
[9] https://dl.dropboxusercontent.com/u/54720544/Molecular% 20Dynamics%20pack/Aug%202011%20v%2015Feb2013%20%2B% 20Next%20Release.zip