24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3512  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂 +6 asdasdassda 2026-04-05 7/350 2026-04-12 00:26 by 丛麟ll
[考研] 一志愿211,0703化学305分求调剂 +26 严西西戏 2026-04-06 33/1650 2026-04-11 23:01 by 314126402
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 考研求调剂 +3 ban班小七 2026-04-11 3/150 2026-04-11 20:48 by may_新宇
[考研] 求调剂 +6 电气300求调剂不 2026-04-08 6/300 2026-04-11 20:14 by 逆水乘风
[考研] 308求调剂 +3 VvvvL 2026-04-10 3/150 2026-04-11 15:32 by xiangxu
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 080500求调剂 +17 黄宇博 2026-04-06 17/850 2026-04-11 08:36 by zhq0425
[考研] 调剂 化学 307 +21 73372112 2026-04-09 23/1150 2026-04-10 23:53 by wj165256
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 求调剂288 +6 ioodiiij 2026-04-10 8/400 2026-04-10 21:07 by zhouxiaoyu
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 调剂求助(生物与医药) +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 302分 085601求调剂推荐 +11 zyx上岸! 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
信息提示
请填处理意见