24小时热门版块排行榜    

查看: 3506  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境学硕288求调剂 +6 皮皮皮123456 2026-03-22 6/300 2026-03-22 16:52 by i_cooler
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见