大家好,我是小马老师。
本文介绍一个lammps合金势的拟合程序。
做过金属模拟的应该都知道Zhou合金势,专栏也对其使用方法进行过详细的介绍。
早期的Zhou势用Fortran编写,在使用之前需要进行对源代码进行编译,很多同学并不熟悉Fortran编译过程,在使用时有些困难。
最新版的lammps提供了一个python版本的Zhou势拟合程序,用起来就简单多了。
无需编译,也不需要编写输入EAM.input文件,只需要运行python命令即可生成所需要的合金势。
如运行以下命令:
python create_eam.py -n Ta Cu
就可以得到TaCu合金势:TaCu.eam.alloy
最新的Zhou势支持的元素为:
Cu, Ag, Au, Ni, Pd, Pt, Al, Pb, Fe, Mo, Ta, W, Mg, Co, Ti, Zr, Cr
这个源代码在lammps官方安装包即可找到。
公号:lammps加油站