24小时热门版块排行榜

返回列表

gulubaozi

银虫 (小有名气)

应助: 5 (幼儿园)
金币: 876.9
红花: 7
帖子: 260
在线: 95.6小时
虫号: 2984608
注册: 2014-02-22
专业: 金属材料的凝固与结晶学

[求助] 关于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之间）

希望大家多给意见，更希望能得到高手的指点!

回复此楼