24小时热门版块排行榜    

查看: 3708  |  回复: 14
【悬赏金币】回答本帖问题,作者半寸日光将赠送您 20 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

半寸日光

新虫 (小有名气)

[求助] 求助:根据文献参数写适合lammps中tersoff势的Albe格式问题 已有1人参与

求助:请问有人根据文献里面的参数写过适合lammps中的tersoff势?可以帮忙看下我根据文献参数所写关于三种元素的tersoff中Albe格式正确吗?(参考lammps自带的SiC的Albe),文献里面的参数,自己所仿照写的W-C-Co体系的势函数如下,感谢!

求助:根据文献参数写适合lammps中tersoff势的Albe格式问题


求助:根据文献参数写适合lammps中tersoff势的Albe格式问题-1


求助:根据文献参数写适合lammps中tersoff势的Albe格式问题-2


@smutao @oxox6085 发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

引用回帖:
10楼: Originally posted by 半寸日光 at 2021-06-10 09:49:28
感谢!
...

再补充一下:写Albe势的格式或Tersoff1势格式时:
以X, Y,  Z三个不同的元素为例,当然了,也可以类推到它们俩俩为为相同的元素。

X  Y  Z 的其中所有与二体势相关的参数(n, beta, lambda2, B, labmda1, A)都与X-Y对的列表值一样。而与三体势相关的参数(m, gamma, lambda3,c, d, costheta0, R, D)的值需要要考虑到与元素C相关的设置,或者说这些三体势相关的参数的值是与X-Z对的参数列表值给出的,而二体势相关的参数有X-Y对的参数列表值给出。
A,B, labmda1, lambda2的值是基于参数r0, D0, S, beta的值按公式计算出来的:
lambda1=beta*sqrt(2.0*S)
lambda2=beta*sqrt(2.0/S)
A= D0/(S-1.0) * exp(lambda1*r0)
B= S*D0/(S-1.0)*exp(lambda2*r0)
11楼2021-06-10 08:55:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

valenhou001

至尊木虫 (职业作家)

我把这个W-C二元的势文件根据这篇文献整理了的。见https://github.com/houzf/ABOP_lammps

发自小木虫Android客户端
2楼2021-06-09 23:35:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

三元的比较复杂。还没有摸清楚转换关系。

发自小木虫Android客户端
3楼2021-06-09 23:39:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

你参考的这篇文献对某个参数的写错了,与原始文献有出了。

发自小木虫Android客户端
4楼2021-06-09 23:42:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见