• 回答数

    3

  • 浏览数

    1961

  • 收藏数

    0

作者:怀念你的爱 发表于 2018-11-19 11:49:43
跳转到指定楼层
我是用2NN做Cu-Si的,也是已知Cu,Si,Cu——Si这些meam参数。
可是问题在于:
library.meam中Si和Cu给出参数与我要用的的各种参数有些不一致,想做一个修改和补充,我是再写一个Si.meam和Cu.meam还是直接在library里修改?还是直接写一个SiCu.meam
这样会存在两个问题,第一个是如何编译library.meam?求大神指教,这个可能对我来说有些难。
第二个是如果编写一个Si.meam和Cu.meam或者SiCu.meam,那么是用什么格式呢?例如,这是library.meam中关于Si的参数
'Si'  'dia' 4.  14   28.086
4.87  4.8 4.8 4.8    4.8 5.431 4.63  1.
1.0   3.30    5.105  -0.80   1.  1
而我要用的参数是
#elt=Si
#lat=dia
#z=4
#ielement=14
#atwt=28.086
#alpha=4.87
b0=3.55
b1=2.50
b2=0.0
b3=7.5
#alat=5.431
#esub=4.63
asub=0.58
#t0=1.0
t1=1.8
t2=5.25
t3=-2.61
#rozero=1
#ibar=0
这些改变直接写在.meam里LAMMPS能够识辨吗?这样该如何处理呢?
谢谢各位指点!感激不尽

分享:
回复

使用道具

该用户从未签到

新手上路

Rank: 1

积分
16
极客币
47
主题
8
帖子
21
注册时间
2018-11-5
在线时间
1 小时
性别
保密
发表于 2018-11-19 11:49:52 | 显示全部楼层
按照手册说法,你可以修改library.meam文件中关于Cu Si 元素的参数,
lammps在从头开始扫描library.meam 文件时,遇到指定的元素符号时,就不再往下扫描了。
我注意到你是做合金体系的,对于合金体系要另写一个CuSi.meam参数文件,
这个文件中有一个rho0参数,也就是ρ0,手册上说这是第二种元素相对于第一种元素的相对密度,
但是在meam势能最先提出的文献里,以及相关的参数拟合文献里,他们用ρ0这个符号表示原子的电子密度。
lammps里rho0到底指的是什么?怎么算?
不知道你有什么看法
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
14
极客币
38
主题
6
帖子
18
注册时间
2018-11-5
在线时间
1 小时
性别
保密
发表于 2018-11-19 11:50:10 | 显示全部楼层
应该CuSi.meam和library.meam都要写,可以参照library.meam自己单独写一个只有两个元素的文件放在运行的文件夹里,具体参数可以在lammps官网手册里面查meam查到,然后我是把写好的和in文件放在一个文件夹里运行,可以辨识的……但现在还不知道算出的结果对不对……
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
15
极客币
43
主题
7
帖子
20
注册时间
2018-11-5
在线时间
1 小时
性别
保密
 楼主| 发表于 2018-11-19 11:50:22 | 显示全部楼层
回复

使用道具 举报

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册 微信登录