| 查看: 2563 | 回复: 0 | ||
[求助]
关于MEAM势能的选择问题
|
|
最近想做纯铁的模拟，用lammps自带的library.meam中的参数，发现结果存在严重错误。 后搜到一组势能参数来自于篇文章http://journals.aps.org/prb/abstract/10.1103/PhysRevB.89.094102（见附件） 做的是纯Fe3C的势能，相关参数文件如下： # Fe3C (Cementite) library MEAM file used in Phys. Rev. B 89, 094102 (2014) # Fe potential referenced from J. Phys.: Condens. Matter 24 225404 # C potential referenced from Phys. Rev. B 81, 104103 # Used together with "Fe3C.meam" for Fe-C alloy potential parameters # elt lat z ielement atwt # alpha b0 b1 b2 b3 alat esub asub # t0 t1 t2 t3 rozero ibar 'Fe' 'bcc' 8 26 55.847000 5.027000 3.800000 2.000000 0.900000 0.000000 2.851000 4.280000 0.585000 1.000000 -0.800000 12.300000 2.000000 1.000000 0.000000 'C' 'dia' 4. 6 12.0111 4.38 4.26 5.00 3.20 3.980 3.567 7.37 1.490 1.0 7.5 1.04 -1.01 1.0 0 来自：http://www.ctcms.nist.gov/potent ... _Liyanage_2014.meam # Fe-C alloy MEAM potential parameters from Phys. Rev. B 89, 094102 # Developed using April 2011 version of LAMMPS (lammps.sandia.gov) # Should work on recent release of LAMMPS # Contact information: lsiliyanage@gmail.com (Dr. Laalitha S. I. Liyanage) augt1=0 ialloy=1 nn2(1,1)=1 nn2(2,2)=1 nn2(1,2)=1 rc=4.500000 repuls(1,1)=0.300000 Cmin(1,1,1)=0.680000 Cmax(1,1,1)=1.900000 Cmin(2,2,2)=2.000000 Cmax(2,2,2)=2.800000 lattce(1,2)='b1' repuls(1,2)=0.125301 attrac(1,2)=0.125301 delta(1,2)=0.002348 alpha(1,2)=4.750000 re(1,2)=1.920000 rho0(1)=1.000000 rho0(2)=5.491125 Cmin(1,1,2)=0.059531 Cmax(1,1,2)=2.800000 Cmin(1,2,1)=2.000000 Cmax(1,2,1)=2.800000 Cmin(1,2,2)=2.000000 Cmax(1,2,2)=2.800000 Cmin(2,2,1)=0.507375 Cmax(2,2,1)=2.800000 来自：http://www.ctcms.nist.gov/potent ... _Liyanage_2014.meam 这两个势能参数文件是作者针对Fe3C的计算拟合结果。 如果用这两个势能参数文件来做纯Fe的模拟是否合适？ 如果是Fe-C合金呢？（成分在0~Fe3C之间） 希望大家多给意见，更希望能得到高手的指点! |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : PhysRevB.89.pdf
2014-08-28 11:53:56, 3.23 M
» 猜你喜欢
一志愿安大生物学07初试322、本科二本、调剂求助
已经有7人回复
-
0856专硕求调剂 希望是a区院校
已经有25人回复
-
085600材料与化工，求调剂
已经有12人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有19人回复
-
材料工程085601，270求调剂
已经有33人回复
-
复试调剂
已经有12人回复
-
新一代电子信息294求调剂 不挑学校
已经有8人回复
-
290调剂生物0860
已经有29人回复
-
279学硕食品专业求调剂院校
已经有15人回复
-
310求调剂
已经有14人回复
找到一些相关的精华帖子，希望有用哦~
LAMMPS和DL-POLY是否能同时调用两个力场
已经有4人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
10