24小时热门版块排行榜    

查看: 1526  |  回复: 17

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



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

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

huangll99

木虫 (职业作家)

引用回帖:
3楼: Originally posted by chenxin1992 at 2012-09-06 13:23:00
这个是单个原子啊,难道我需要扩大点?...

我觉得太小了,至少也要10A以上吧,不过不清楚这是不是报错的原因
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一样基础)都没有学习,直接拿过软件来凭自己的想象就用的~就不要让别人给自己的懒惰买单了吧,这样不怎么好呢(就好像自己不学英语,却要求导师把文献翻译成中文一样,不翻译就不知道自己查词典去看)。
思想重于技巧,内涵重于表象
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不觉得该看看书而不是在这里猜么?
思想重于技巧,内涵重于表象
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +8 小楼。 2026-03-12 11/550 2026-03-16 09:46 by 无际的草原
[考研] 材料专硕306英一数二 +3 z1z2z3879 2026-03-16 3/150 2026-03-16 09:19 by Demonsssss
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见