应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 怎么算不了单原子
181/2返回列表12下一页
查看: 1798  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chenxin1992

银虫 (正式写手)

[求助] 怎么算不了单原子

我在算石墨烯吸附原子的能量，需要计算一次被吸附原子的能量
但是算不出来啊，算几秒钟就提示
LAPACK: Routine ZPOTRF failed! 5



我的INCAR文件是这样的
SYSTEM = Au
NELMDL = 5
ISMEAR = 0
SIGMA  = 0.05
PREC   = ACCURATE

POSCAR是这样的
Au
1.0
        1.0000000000         0.0000000000         0.0000000000
        0.0000000000         1.0000000000         0.0000000000
        0.0000000000         0.0000000000         1.0000000000
   Au
    1
Cartesian
    0        0        0




KPOINT是这样的
1x1x1 au
0
Monkhorst
1 1 1
0 0 0


POTCAR用的是PAW的


我用的是VASP5.2



问题出在哪里啊，求解
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-06 13:05:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

huangll99

木虫 (职业作家)
引用回帖:
3楼: Originally posted by chenxin1992 at 2012-09-06 13:23:00
这个是单个原子啊，难道我需要扩大点？...

我觉得太小了，至少也要10A以上吧，不过不清楚这是不是报错的原因
赞一下(1人) 回复此楼
4楼2012-09-06 13:34:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
引用回帖:
6楼: Originally posted by chenxin1992 at 2012-09-06 17:20:05
是这个原因，但是我算出来Au的free energy才 -0.01571833 eV
氢原子的能量都是-13.6ev,怎么回事啊？...

LZ先花点时间搞清楚，H原子的什么能量是-13.6eV。并且也搞清楚，这个Free Energy在VASP中是如何定义的。
LZ的问题看就是什么理论基础（是基础中的基础，就好像读书你要先认字，做三角函数你要先知道1+1=2一样基础）都没有学习，直接拿过软件来凭自己的想象就用的～就不要让别人给自己的懒惰买单了吧，这样不怎么好呢（就好像自己不学英语，却要求导师把文献翻译成中文一样，不翻译就不知道自己查词典去看）。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
12楼2012-09-07 14:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
引用回帖:
15楼: Originally posted by chenxin1992 at 2012-09-07 17:37:32
难道不是电子的动能和势能么？...

任何东西的能量都可以写成动能和势能两个部分，这是物理学的基础啊……
动能这东西好办，一般在给定的坐标体系中是绝对量。那势能呢？从初中的学习就知道，势能没有什么绝对量，只有选定基准后的相对量。
在第一性原理中，势能又是怎样定义的，这个能量数值究竟是怎么给出的，LZ不觉得该看看书而不是在这里猜么？
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
17楼2012-09-07 20:23:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

huangll99

木虫 (职业作家)
★
sunyang1988: 金币+1, 谢谢交流 2012-09-06 14:43:45
超胞尺寸只有1A吗，为啥这么小啊？
赞一下 回复此楼
2楼2012-09-06 13:21:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)
引用回帖:
2楼: Originally posted by huangll99 at 2012-09-06 13:21:43
超胞尺寸只有1A吗，为啥这么小啊？

这个是单个原子啊，难道我需要扩大点？
赞一下 回复此楼
3楼2012-09-06 13:23:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

06022043

铁杆木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-09-06 20:09:37
Au 应该有晶体结构把？为什么不算呢？气体才在一个大超胞中算，超胞至少要 10 A。
赞一下 回复此楼
专心做学问
5楼2012-09-06 15:17:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)
引用回帖:
4楼: Originally posted by huangll99 at 2012-09-06 13:34:48
我觉得太小了，至少也要10A以上吧，不过不清楚这是不是报错的原因...

是这个原因，但是我算出来Au的free energy才 -0.01571833 eV
氢原子的能量都是-13.6ev,怎么回事啊？
赞一下 回复此楼
6楼2012-09-06 17:20:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)
引用回帖:
5楼: Originally posted by 06022043 at 2012-09-06 15:17:36
Au 应该有晶体结构把？为什么不算呢？气体才在一个大超胞中算，超胞至少要 10 A。

我设了15埃
回复此楼
7楼2012-09-06 17:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hcdd1986

铜虫 (小有名气)
我用的是MS软件，对你这个不太了解~ MS就是建模吸附然后开算，一下就出来了。
赞一下 回复此楼
8楼2012-09-06 17:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxin1992

银虫 (正式写手)
引用回帖:
8楼: Originally posted by hcdd1986 at 2012-09-06 17:44:08
我用的是MS软件，对你这个不太了解~ MS就是建模吸附然后开算，一下就出来了。

我用MS建模的，用VASP算
赞一下 回复此楼
9楼2012-09-06 17:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

06022043

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by chenxin1992 at 2012-09-06 17:20:05
是这个原因，但是我算出来Au的free energy才 -0.01571833 eV
氢原子的能量都是-13.6ev,怎么回事啊？...

单原子的能量都很底，因为它于其他原子没有相互作用，我认为你应该算H2的能量。然后除以2得H的能量。
赞一下 回复此楼
专心做学问
10楼2012-09-07 08:32:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenxin1992 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +26 医学老男孩 2026-08-13 56/2800 2026-08-15 13:47 by 医学老男孩
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:41 by KxMI1BYBxWX1
[基金申请] 时间戳又变了8-15 +5 archvillain 2026-08-15 9/450 2026-08-15 12:06 by archvillain
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 8/400 2026-08-15 11:53 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[基金申请] filecode +7 cratir 2026-08-14 11/550 2026-08-15 06:19 by 学员8dgXkO
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
关闭
请填处理意见
关闭