24小时热门版块排行榜    

查看: 3504  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿东华大学化学070300,求调剂 +5 2117205181 2026-03-21 6/300 2026-03-22 14:18 by ColorlessPI
[考研] 环境学硕288求调剂 +4 皮皮皮123456 2026-03-22 4/200 2026-03-22 11:07 by 杨杨杨紫
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
请填处理意见