24小时热门版块排行榜    

查看: 714  |  回复: 6
当前主题已经存档。

yangjun102

银虫 (小有名气)

[交流] 【求助】关于分子动力学能的求解

分子动力学中有许多能量的求解:
1.内聚能。2.结合能。3.空位能。4.双空位形成能。5.自间隙原子形成能。
6.点缺陷形成能。7.表面能。
我列出的可以有的重复了,这里提供一个大家讨论的平台吧。
大家可以解释一下自己熟悉的,最好能结合具体的程序来给出具体的求解的方法。
谢谢!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

lei0736

荣誉版主 (职业作家)

优秀版主

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zeoliters(金币+5,VIP+0):多谢回帖交流! 1-8 23:26
这个没有现成程序的 因为没有必要
首先要明白各个能量的定义 一般都可以根据定义直接算
比如内聚能 以FCC铜Pt为例 实验值一般为5.8电子伏
计算值也差别不大 不同的势算出来会小有差别
具体计算就是建一个FCC结构 用相应的EAM或其他势描述相互作用
计算0K下驰豫的总能除以原子个数就行了
空位能就是取走一个原子 驰豫 得到的总能减去N-1个完整FCC原子的总能
其他依次类推
2楼2010-01-06 21:50:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangjun102

银虫 (小有名气)

引用回帖:
Originally posted by lei0736 at 2010-1-6 21:50:
这个没有现成程序的 因为没有必要
首先要明白各个能量的定义 一般都可以根据定义直接算
比如内聚能 以FCC铜Pt为例 实验值一般为5.8电子伏
计算值也差别不大 不同的势算出来会小有差别
具体计算就是建一个FCC ...

谢谢回答,那么表面能的求解怎么办呢?
例如(100),(110),(111)面的求解。
3楼2010-01-07 17:40:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

judy1931

铜虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 1-7 20:02
这些能量应该都可以从MD里头计算出来吧,对其他的不是太清楚,但是内聚能是可以从Amophocell里计算的
4楼2010-01-07 18:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zeoliters(金币+3,VIP+0):多谢回帖交流! 1-8 23:26
从表面能的定义出发
建一个表面体系 算总能 减去相同数目原子的结合能
得到的能量差或除以表面积 或除以表面原子总数
分部得到以电子伏每平方厘米或电子伏每原子为单位的表面能值
需要注意的是 如果建了两个表面就要除以两个表面的面积或原子数
5楼2010-01-07 20:04:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangjun102

银虫 (小有名气)

引用回帖:
Originally posted by lei0736 at 2010-1-7 20:04:
从表面能的定义出发
建一个表面体系 算总能 减去相同数目原子的结合能
得到的能量差或除以表面积 或除以表面原子总数
分部得到以电子伏每平方厘米或电子伏每原子为单位的表面能值
需要注意的是 如果建了两个表 ...

再次感谢lei0736 的解疑!
6楼2010-01-07 20:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangjun102

银虫 (小有名气)

引用回帖:
Originally posted by lei0736 at 2010-1-7 20:04:
从表面能的定义出发
建一个表面体系 算总能 减去相同数目原子的结合能
得到的能量差或除以表面积 或除以表面原子总数
分部得到以电子伏每平方厘米或电子伏每原子为单位的表面能值
需要注意的是 如果建了两个表 ...

还是有点问题,用什么系综去算呢!
我先是采用npt再用nvt,和先用npt再用nve相差还是挺大的。还是其他的方法呢?
继续请教中....
7楼2010-01-10 18:41:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangjun102 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见