and follow the same process for the other cases.
Of course an easier option is to use one of the EAM or MEAM models that comes with GULP in the libraries meam, meam_2nn, meidavenport, suttonchen etc.
It looks like these files have the same format as the ones used in Lammps. Indeed, if these 3 files are put together in a single file, then this would correspond to the Lamps potential file.
Should I use the following?
manubody numeric
Al AlAl.gulp.pair
eam_fn numeric
Al Al.gulp.eamfn
eam_density numeric
Al Al.gulp.denfn
Apologies but my memories were a bit hazy since it’s been a while since I’ve looked at that part of the code. It appears that at present only the numerical format is supported for the functional. While there are options to input spline points for the density and potential they use GULP’s input format and so the files would have to be converted.
The focus of GULP has tended to be on analytic forms (it makes fitting easier) and so it would be better to work with one of these if possible as in the library files. Without looking at the specific model you are trying to use, it’s generally the case that MEAM-2NN improves on EAM for a lot of things.
I will use the meam_2nn library. Is it possible to calculate phonon dispersions with this MEAM potential? I get the following error:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! ERROR : Finite difference phonons only allowed for gamma point
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Unfortunately it’s not currently possible to do phonon dispersion with MEAM. This is because the analytic second derivatives are not yet available (it gets complicated due to the use of a manybody screening term within a manybody formalism) and so only the numerical gamma point force constants are available. Sorry. I’ll add it to the list of things to try to implement, but it’s not a trivial thing.