版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

gulubaozi

銀蟲(chóng) (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 876.9
紅花: 7
帖子: 260
在線: 95.6小時(shí)
蟲(chóng)號(hào): 2984608
注冊(cè): 2014-02-22
專業(yè): 金屬材料的凝固與結(jié)晶學(xué)

[求助] 關(guān)于MEAM勢(shì)能的選擇問(wèn)題

最近想做純鐵的模擬，用lammps自帶的library.meam中的參數(shù)，發(fā)現(xiàn)結(jié)果存在嚴(yán)重錯(cuò)誤。
后搜到一組勢(shì)能參數(shù)來(lái)自于篇文章http://journals.aps.org/prb/abstract/10.1103/PhysRevB.89.094102（見(jiàn)附件）
做的是純Fe3C的勢(shì)能，相關(guān)參數(shù)文件如下：
# 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
來(lái)自：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
來(lái)自：http://www.ctcms.nist.gov/potent ... _Liyanage_2014.meam

這兩個(gè)勢(shì)能參數(shù)文件是作者針對(duì)Fe3C的計(jì)算擬合結(jié)果。
如果用這兩個(gè)勢(shì)能參數(shù)文件來(lái)做純Fe的模擬是否合適？
如果是Fe-C合金呢？（成分在0~Fe3C之間）

希望大家多給意見(jiàn)，更希望能得到高手的指點(diǎn)!

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : PhysRevB.89.pdf

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

» 猜你喜歡

材料調(diào)劑已經(jīng)有5人回復(fù)
081700 調(diào)劑 267分已經(jīng)有8人回復(fù)
請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？已經(jīng)有5人回復(fù)
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷已經(jīng)有10人回復(fù)
一志愿南航材料專317分求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)學(xué)碩求調(diào)劑已經(jīng)有5人回復(fù)
284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)

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

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

智能機(jī)器人

Robot (super robot)

我們都愛(ài)小木蟲(chóng)

找到一些相關(guān)的精華帖子，希望有用哦~

LAMMPS和DL-POLY是否能同時(shí)調(diào)用兩個(gè)力場(chǎng) 已經(jīng)有4人回復(fù)
【討論】lammps-怪異的group溫度已經(jīng)有9人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 gulubaozi 的主題更新