24小时热门版块排行榜    

查看: 1958  |  回复: 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 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见