24小时热门版块排行榜    

查看: 3460  |  回复: 1

hhy618

新虫 (初入文坛)

[求助] 关于用lammps计算金属表面能的问题

这是我从课件上找来计算Fe金属表面能的的程序

units   metal
boundary p p p
atom_style atomic

variable     a0 equal 2.866
lattice     bcc ${a0}  orient x 1 0 -1 orient y 0 1 0 orient z 1 0 1
region       box block 0 3  0 3  0 3
create_box     1 box
create_atoms   1 box

pair_style     meam
pair_coeff     * * library.meam Fe Fe.meam Fe

min_style     cg
minimize     1.0e-10 1.0e-8 10000 10000

variable     Etmp equal pe
variable     E0 equal ${Etmp}

change_box    all z delta 0 5

minimize    1.0e-10 1.0e-8 1000 1000

variable  area equal lx*ly
variable  Esurf equal (${Etmp}-${E0})/(2*${area})*16022

print   "@ Esurf: ${Esurf}"

但是在计算相比较100 111面的时候发现  110面的表面能并非是最小的 与理论值不符啊`` 麻烦哪位指点一下 看看程序是不是有错误
Ps 如果我还想用该程序运算Cu Ni Au的表面能的话 除了需要改动a0 晶体结构 还有调用的meam势文件之外 还要改动别的地方么?
拜托了。 。 。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

刘雪鹏

木虫 (小有名气)

请问a0就是晶格常数吗?
2楼2013-12-30 00:09:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hhy618 的主题更新
信息提示
请填处理意见