24小時熱門版塊排行榜

返回列表

gulubaozi

銀蟲 (小有名氣)

應助: 5 (幼兒園)
金幣: 876.9
紅花: 7
帖子: 260
在線: 95.6小時
蟲號: 2984608
注冊: 2014-02-22
專業(yè): 金屬材料的凝固與結晶學

[求助] 關于MEAM勢能的選擇問題

最近想做純鐵的模擬，用lammps自帶的library.meam中的參數，發(fā)現結果存在嚴重錯誤。
后搜到一組勢能參數來自于篇文章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之間）

希望大家多給意見，更希望能得到高手的指點!

回復此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內容負責。
本內容由用戶自主發(fā)布，如果其內容涉及到知識產權問題，其責任在于用戶本人，如對版權有異議，請聯系郵箱：xiaomuchong@tal.com
附件 1 : PhysRevB.89.pdf

2014-08-28 11:53:56, 3.23 M

» 猜你喜歡

求調劑已經有3人回復
生物學學碩求調劑已經有5人回復
284求調劑已經有10人回復
一志愿山東大學藥學學碩求調劑已經有4人回復
07化學280分求調劑已經有4人回復
298-一志愿中國農業(yè)大學-求調劑已經有12人回復
求材料，環(huán)境專業(yè)調劑已經有3人回復
335求調劑已經有5人回復
求調劑已經有7人回復
一志愿吉大化學322求調劑已經有4人回復

1樓 2014-08-28 11:55:06

已閱關注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

LAMMPS和DL-POLY是否能同時調用兩個力場已經有4人回復
【討論】lammps-怪異的group溫度已經有9人回復

點擊這里搜索更多相關資源

科研從小木蟲開始，人人為我，我為人人

相關版塊跳轉我要訂閱樓主 gulubaozi 的主題更新