24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者半寸日光将赠送您 20 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

半寸日光

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 800.2
帖子: 259
在线: 40.3小时
虫号: 23577820
注册: 2020-09-21
专业: 机械工程

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

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

@smutao @oxox6085 发自小木虫Android客户端

回复此楼

» 猜你喜欢

1楼 2021-06-09 19:50:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

引用回帖:

8楼: Originally posted by valenhou001 at 2021-06-10 08:01:24
我现在对lammps里的这个Albes势格式不确定的地方就是这个n的值。从Albe势的公式里，看不出这个n对应的是哪个物理量。lammps程序里附带的GaN的Albe势里，n的值都是设置1的。我测试过，将n=0与n=1，其他的都不变的话 ...

补充h对应的是cos(\theta_0)，还需要注意h的正负号。Albe势，在很多文献里写的是[h+cos(\theta_{ijk})]^2。Tersoff势在很文献里写成的是[cos(\theta_{ijk}-h]^2或[h-cos(\theta_{ijk}]^2。 lammps主页上对Tersoff势表达式写的是[cos(\theta_{ijk}-h]^2，这里h=cos(\theta_0)。