| 查看: 2597 | 回复: 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
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
LAMMPS和DL-POLY是否能同时调用两个力场
已经有4人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30