合金体系的势函数除了eam势,还有meam势。在新版本的lammps中,meam势类型已经改为meam/c,本文主要介绍meam/c势的设置方法。
和普通的势函数文件不同,meam/c势有两个势函数文件:library.meam和**.meam,**表示不同的势函数名称。
- library.meam通用文件
library.meam类似于参数库,存储了多种元素的meam参数,每一种元素占3行,共19个通用参数。
例如,下面一段代码是ZrCu合金的library.meam文件:
#elt lat z iel atwt
#alpha b0 b1 b2 b3 alat esub asub
#t0 t1 t2 t3 rhozero ibar
'Zr' 'hcp' 12 1 91.2200