24小时热门版块排行榜    

查看: 1919  |  回复: 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文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:21 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
请填处理意见