| 查看: 3490 | 回复: 2 | ||
[求助]
EAM势能模拟的文件参数代表什么？ 已有1人参与
|
|
大家好！我在做Cu-Zr合金的模拟，使用lammps软件，EAM势能。在网站http://www.ctcms.nist.gov/potentials/Cu-Zr.html 找到相关的势参数文件CuZr.eam.fs，根据lammps手册介绍，是DYNAMO setfl格式。CuZr.eam.fs文件前6行我都看懂了。后面有很多数据，共5列，请问分别代表什么含义？ Sourse: M.I. Mendelev, unpublished Contact information: mendelev@ameslab.gov Thursday, May 29, 2008 The potential was taken from v5_14_CuZr_B2 (in C:\SIMULATION.MD\Cu-Zr\v5_14) 2 Cu Zr 10000 3.00000000000000E-0002 10000 7.60000000000000E-0004 7.60000000000000E+0000 29 6.35460000000000E+0001 3.63908700000000E+0000 fcc 0 -1.39576844400000E-0001 -1.97391466400000E-0001 -2.41754186100000E-0001 -2.79153688900000E-0001 -3.12103312200000E-0001 -3.41892048800000E-0001 -3.69285619100000E-0001 -3.94782932600000E-0001 -4.18730533200000E-0001 -4.41380736900000E-0001 -4.62924022300000E-0001 -4.83508372100000E-0001 -5.03251469300000E-0001 -5.22248730900000E-0001 -5.40578793900000E-0001 -5.58307377600000E-0001 -5.75490072400000E-0001 -5.92174399100000E-0001 -6.08401359600000E-0001 -6.24206624400000E-0001 -6.39621454800000E-0001 -6.54673430700000E-0001 -6.69387030300000E-0001 -6.83784097400000E-0001 -6.97884222000000E-0001 -7.11705053200000E-0001 -7.25262558200000E-0001 -7.38571238100000E-0001 -7.51644310400000E-0001 -7.64493861800000E-0001 -7.77130980100000E-0001 -7.89565865400000E-0001 -8.01807926700000E-0001 -8.13865865400000E-0001 -8.25747747400000E-0001 -8.37461066400000E-0001 -8.49012799300000E-0001 -8.60409454100000E-0001 -8.71657113900000E-0001 -8.82761473800000E-0001 -8.93727875400000E-0001 -9.04561336100000E-0001 -9.15266576600000E-0001 -9.25848044600000E-0001 -9.36309936500000E-0001 -9.46656216700000E-0001 -9.56890635400000E-0001 -9.67016744300000E-0001 -9.77037910800000E-0001 -9.86957331800000E-0001 -9.96778044600000E-0001 -1.00650293870000E+0000 -1.01613476520000E+0000 -1.02567614610000E+0000 -1.03512958230000E+0000 -1.04449746170000E+0000 -1.05378206620000E+0000 -1.06298557780000E+0000 -1.07211008490000E+0000 -1.08115758780000E+0000 -1.09013000380000E+0000 -1.09902917190000E+0000 -1.10785685720000E+0000 -1.11661475520000E+0000 |
回复此楼» 猜你喜欢
参与限项
已经有4人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps的minimize会丢失原子
已经有6人回复
- 有没有用 lj势函数 模拟lammps的朋友？
已经有6人回复
- 请问gromacs导入EAM势能函数的问题和lammps导入martini粗粒化力场的问题
已经有4人回复
- ms建模指定力场导入lammps后是否需要引入势函数
已经有3人回复
- 如何利用GULP拟合分子间势能曲线得出力场参数
已经有10人回复
- 刚开始学习LAMMPS，求助关于lattice的设置
已经有7人回复
- lammps里面的pair_coeff设置问题
已经有3人回复
- 关于势函数的求助
已经有23人回复
- 【求助】分子的LJ势能参数的拟合
已经有8人回复
- 【求助】如何根据势能曲线拟合力场参数
已经有19人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fyc801: 金币+50, ★★★很有帮助, 不好意思，上次忘了给金币。现在才想起来 2014-09-22 23:11:18
感谢参与，应助指数 +1
fyc801: 金币+50, ★★★很有帮助, 不好意思，上次忘了给金币。现在才想起来 2014-09-22 23:11:18
|
首先，你要知道EAM势包含两部分，一个是原子嵌入势，跟电子云密度ρ相关。另一个是 两体势，可以理解为范德华力，跟原子间距r相关。 因此，你可以看到第5行，第一个10000表示将电子云密度函数数值等分为10000份 第二个10000表示将距离r等分为10000份。这样就将连续函数离散化了（连续函数 即是你所用的势函数的解析表达式）。 下面的数值，即是F(ρ),ρ,U(r)的函数值，F表示原子嵌入势，ρ表示电子云密度 U表示两体势。F,ρ的数值个数是由第一个10000决定，U的个数是由第二个10000决定 排布方式为F,ρ,U 如果有两个元素以上，则排布方式是F1，ρ1，F2，ρ2，U1，U2.以此类推 |
2楼2014-09-03 10:33:06
3楼2016-06-13 16:54:01