| 查看: 3862 | 回复: 14 | |||
| 【悬赏金币】回答本帖问题，作者半寸日光将赠送您 20 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求助：根据文献参数写适合lammps中tersoff势的Albe格式问题 已有1人参与
|
|||
求助：请问有人根据文献里面的参数写过适合lammps中的tersoff势？可以帮忙看下我根据文献参数所写关于三种元素的tersoff中Albe格式正确吗？(参考lammps自带的SiC的Albe)，文献里面的参数，自己所仿照写的W-C-Co体系的势函数如下，感谢!   @smutao @oxox6085 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
（081700）化学工程与技术-298分求调剂
已经有33人回复
-
0703化学求调剂
已经有7人回复
-
一志愿西南交大，材料专硕317求调剂
已经有3人回复
-
315求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
293求调剂
已经有3人回复
-
求b区学校调剂
已经有3人回复
-
310求调剂
已经有3人回复
-
302求调剂
已经有6人回复
-
求调剂
已经有7人回复
valenhou001
至尊木虫 (职业作家)
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
感谢参与，应助指数 +1
|
我现在对lammps里的这个Albes势格式不确定的地方就是这个n的值。从Albe势的公式里，看不出这个n对应的是哪个物理量。lammps程序里附带的GaN的Albe势里，n的值都是设置1的。我测试过，将n=0与n=1，其他的都不变的话，两个势给的结果基本一样。这个我在我的github做过注释。 从GaN势的出处文献来看，找不到n在公式里对应的变量。 从10.1103/PhysRevB.86.115410文献里（针对的BNC体系的Albes或Tersoff势）的公式来看，n是一个指数因子，并也出现在一个倒数的分母里，所以n=0似乎是不合理的。这也是让我困惑的地方。 你参考一下BNC体系的Tersoff势： https://www.ctcms.nist.gov/poten ... B-N-C/1/BNC.tersoff 这个势文件里，n并不等于0或1. 从这些文献和相应的势文件来看，这个n的取值应根据文献里所给的公式，看它们是怎么定的b_{ij}或B_{ij}的表达式（也就是repulsive potential前面的一个系数）。 原本的Tersoff势表达式中b_{ij}表达式比较复杂的。而Albe势对这个b_{ij}简化了，简化为(1+\chi_{ij})^{-1/2}，这样n就取为1了。 你可以参考B-N-C体系的Tersoff势文件（我提到的那篇PRB文献）来试试写W-C-Co三元体系的Albe势文件。 与三体相关的项只是b_ij、\chi_{ij}、g(\theta_{ijk})、cos(theta_{ijk})的项。那篇W-C-Co的势参考文件中势参数列表中h对应的就是cos(\theta_{ijk})的值。2\mu对应着lammps的Tersoff势参数格式中的\lambda_{3}。 |
8楼2021-06-10 07:01:24
valenhou001
至尊木虫 (职业作家)
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2021-06-09 23:35:23
valenhou001
至尊木虫 (职业作家)
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
3楼2021-06-09 23:39:33
valenhou001
至尊木虫 (职业作家)
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
4楼2021-06-09 23:42:38
10