24小时热门版块排行榜    

查看: 1853  |  回复: 19
当前主题已经存档。

ron_hai

铜虫 (初入文坛)

[交流] 【求助】VASP计算ENCUT大了10.能量就为负了?请大家帮忙看看。。。

这是计算氧化铅的bulk时出现的情况

ENCUT=530时
  free energy    TOTEN  =       455.43723802 eV

  energy without entropy =      455.54301228  energy(sigma->0) =      455.49012515

ENCUT=540时
  free energy    TOTEN  =       459.65204558 eV

  energy without entropy =      459.75624941  energy(sigma->0) =      459.70414750

ENCUT=550时
  free energy    TOTEN  =       -24.32559093 eV

  energy without entropy =      -24.32529110  energy(sigma->0) =      -24.32544101

前面都很正常。。就是ENCUT=550时能量就突然为负值了。。。
请大家看看帮帮忙。。谢了。。。

[ Last edited by wuchenwf on 2009-6-23 at 20:11 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)


wuchenwf(金币+1):谢谢
没有遇到过这种问题,你是否设错什么了
2楼2009-04-15 19:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)


wuchenwf(金币+1):谢谢
总能量不就应该是负的吗
3楼2009-04-15 20:28:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ron_hai

铜虫 (初入文坛)

引用回帖:
Originally posted by jghe at 2009-4-15 20:28:
总能量不就应该是负的吗

但是之前从350~540的全部都是递增的正值啊。。。
很郁闷的说。。。

进行计算的其他参数均未变化就改了ENCUT而已。
4楼2009-04-16 09:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化学小工

木虫 (正式写手)

★ ★
wuchenwf(金币+2):谢谢
ENCUT的使用需要考虑potcar,你好好看看,然后再改,这个值不是随便设置的,不过大家谁来回答一下,为什么他的能量是正的?是不是输出的设置,显示成正数?
有问题可以联系我QQ676869820
5楼2009-04-16 09:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ron_hai

铜虫 (初入文坛)

引用回帖:
Originally posted by 化学小工 at 2009-4-16 09:50:
ENCUT的使用需要考虑potcar,你好好看看,然后再改,这个值不是随便设置的,不过大家谁来回答一下,为什么他的能量是正的?是不是输出的设置,显示成正数?

POTCAR不是两个给定的文件CAT成的吗?
而且我做这个也只是做初步分析,做总能的收敛曲线。
就是为了以后做优化用的。

为何为正偶也不清楚呀。初学者。。请大家多指教。。
6楼2009-04-16 10:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ron_hai

铜虫 (初入文坛)

多发了一个。。呵呵。。斑斑见谅哈。。
7楼2009-04-16 10:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwyan

新虫 (正式写手)

★ ★
wuchenwf(金币+2):谢谢
ECUT增加时能量应该是减小的,总能量应该是负值,
在赝势中有每种赝势 的ECUT,INCAR中设置ECUT时,是要根据赝势中的ECUT数值来设的,至少大于赝势中的ECUT(有多个赝势时取最大的),一般是1.3倍(好像是)。
事了拂衣去,千里不留行
8楼2009-04-16 14:47:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)


zdhlover(金币+1,VIP+0):谢谢 4-17 18:14
引用回帖:
Originally posted by xwyan at 2009-4-16 14:47:
ECUT增加时能量应该是减小的,总能量应该是负值,
在赝势中有每种赝势 的ECUT,INCAR中设置ECUT时,是要根据赝势中的ECUT数值来设的,至少大于赝势中的ECUT(有多个赝势时取最大的),一般是1.3倍(好像是)。

"总能量应该是负值"
不要说的这么绝对,你可以摆一个不稳定的构型。

[ Last edited by y1ding on 2009-4-17 at 16:27 ]
9楼2009-04-17 16:25:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ron_hai

铜虫 (初入文坛)

★ ★
fegg7502(金币+2,VIP+0):鼓励交流! 4-19 16:05
感谢上面几位。。为什么最后一个计算有出入我已经弄清楚了。是我拷结果的时候太急了。。还没算完。。。囧。。。。

我这组结果都是递增的正数哦。。

补上
ENCUT=550时结果

  free energy    TOTEN  =       461.09460683 eV

  energy without entropy =      461.19944857  energy(sigma->0) =      461.14702770

[ Last edited by ron_hai on 2009-4-19 at 09:39 ]
10楼2009-04-19 09:37:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ron_hai 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 jRl3mE6ddZGq 2026-08-03 10/500 2026-08-04 09:26 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见