24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3514  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 一志愿211生物学280分 求调剂 +7 李rien 2026-04-05 7/350 2026-04-11 11:16 by 逆水乘风
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 一志愿 江南大学 085602 化工专硕 338分求调剂 +16 路痴小琪 2026-04-05 16/800 2026-04-10 08:08 by kangsm
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 生物学363调剂求助 +7 fanzhang6666 2026-04-06 9/450 2026-04-07 17:37 by lijunpoly
信息提示
请填处理意见