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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 【求助】VASP计算ENCUT大了10.能量就为负了？请大家帮忙看看。。。
北京石油化工学院2026年研究生招生接收调剂公告
202/2返回列表上一页12
查看: 1694  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

ron_hai

铜虫 (初入文坛)
又碰到一个问题了。。总能应该是ENCUT达到一定程度就收敛的。。。
我这个ENCUT都从350一直加10加到800了。。总能还是没有收敛。。555.。
请问大家是不是INCAR还需修改？
赞一下 回复此楼
11楼2009-04-19 17:42:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fegg7502

荣誉版主 (知名作家)

无

优秀版主
★ ★
wuchenwf(金币+2):3Q！
随着encut的增大，能量结果应该是波动的，，但是波动振幅逐渐减小，这点我看回帖的人说的不怎么正确
赞一下(1人) 回复此楼
心诚意正，方可始终；身无苦灭，永登光明。
12楼2009-04-19 22:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)
★
wuchenwf(金币+1):3Q！
引用回帖:
Originally posted by xwyan at 2009-4-16 14:47:
ECUT增加时能量应该是减小的，总能量应该是负值，
在赝势中有每种赝势 的ECUT，INCAR中设置ECUT时，是要根据赝势中的ECUT数值来设的，至少大于赝势中的ECUT（有多个赝势时取最大的），一般是1.3倍（好像是）。

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^非也非也
赞一下(1人) 回复此楼
http://sim.phys.ruc.edu.cn
13楼2009-04-19 23:15:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)
你这个总能根本不对，你先算个负数出来再说。
引用回帖:
Originally posted by ron_hai at 2009-4-19 09:37:
感谢上面几位。。为什么最后一个计算有出入我已经弄清楚了。是我拷结果的时候太急了。。还没算完。。。囧。。。。

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

补上
ENCUT=550时结果

  free energy    TOTEN  =    ...

赞一下 回复此楼
http://sim.phys.ruc.edu.cn
14楼2009-04-19 23:17:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)
★
wuchenwf(金币+1):3Q！ 紫霞上的老朋友？
还弄清楚了呢。远着呢！

你的计算要么出错，要么结构偏离平衡位置太远了。。。
引用回帖:
Originally posted by ron_hai at 2009-4-19 09:37:
感谢上面几位。。为什么最后一个计算有出入我已经弄清楚了。是我拷结果的时候太急了。。还没算完。。。囧。。。。

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

补上
ENCUT=550时结果

  free energy    TOTEN  =    ...

赞一下(1人) 回复此楼
15楼2009-04-20 08:51:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwyan

新虫 (正式写手)
引用回帖:
Originally posted by wei.ji at 2009-4-19 23:15:

                        ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^非也非也

楼上您好
我以前在计算遇到的总能，都是负值，于是感觉一个稳定的基态的能量应该是负值。我看了一下VASP中OUTCAR中的能量组成，确实对各部分的意义不理解，请教一下各部分的意义？ 另外，在什么情况下，总能为负值？
Free energy of the ion-electron system (eV)
  ---------------------------------------------------
  alpha Z        PSCENC =         2.18821992
  Ewald energy   TEWEN  =      -135.42777035
  -1/2 Hartree   DENC   =        -0.51242094
  -V(xc)+E(xc)   XCENC  =        -9.08154305
  PAW double counting   =        34.39875164      -16.63333566
  entropy T*S    EENTRO =        -0.00000038
  eigenvalues    EBANDS =        18.95066519
  atomic energy  EATOM  =       103.30142565
  ---------------------------------------------------
  free energy    TOTEN  =        -2.81600799 eV

  energy without entropy =       -2.81600761  energy(sigma->0) =       -2.81600780


十分感谢回复
赞一下 回复此楼
事了拂衣去，千里不留行
16楼2009-04-20 09:15:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwyan

新虫 (正式写手)
版主说随着ECUT的增加，能量应该是波动的，能否解释一下为什么波动，而不是一直下降的去接近一个极限值。
谢谢回复

[ Last edited by xwyan on 2009-4-20 at 09:21 ]
赞一下 回复此楼
事了拂衣去，千里不留行
17楼2009-04-20 09:20:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ron_hai

铜虫 (初入文坛)
本人错的很离谱。。
TOTEN找错位置了。。
寒。。都是负的
ENCUT       TOTEN
350        -21.84353
360        -21.82009
370        -21.8091
380        -21.80035
390        -21.79253
400        -21.78402
410        -21.77879
420        -21.77879
430        -21.76246
440        -21.75823
450        -21.75428
460        -21.75265
470        -21.75003
480        -21.74772
490        -21.74656
500        -21.74622
510        -21.74615
520        -21.7462
530        -21.74636
540        -21.74659
550        -21.74698

ENCUT到500以后就已经收敛了。。
赞一下 回复此楼
18楼2009-04-20 11:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limaojlu

金虫 (正式写手)
★ ★
zxzj05(金币+2,VIP+0):3Q！ 期待更多讨论！ ^_^ 4-22 11:41
这是VASP的总能分项
Alpha Z    300.56    300.56    300.56
  Ewald    -3225.62  -3223.88  -3207.89     84.76     55.64     46.52
  Hartree    367.42    453.46    405.96     41.11      6.01     49.42
  E(xc)    -1360.87  -1360.94  -1360.95      0.64      0.30     -0.24
  Local    -1735.97  -1846.47  -1801.28   -114.72    -49.21   -104.11
  n-local   4325.27   4325.40   4325.69     -3.12      0.12      0.29
  augment  -1303.10  -1308.21  -1305.18      2.22      2.57     -1.87
  Kinetic   2632.27   2660.12   2643.05    10.84    -15.38      9.97
Hartree，E(xc)  是电子与电子相互作用，Ewald是离子离子的相互作用，其他项不太明确，计算总能是正负值我想是不是通过这几项来说明了。
赞一下(10人) 回复此楼
19楼2009-04-21 15:05:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

langes

铜虫 (小有名气)
★ ★
wuchenwf(金币+2,VIP+0):谢谢 4-23 22:44
vasp手册中都有说明的，对于ENCUT选项，要首先对计算体系进行参数优化，直到体系的能量变化小于20meV（好像是）就可以判断其收敛了。
如果体系的能量为正值，我认为很可能是构型出错了，这种构型自然界是不存在的，是我们自己拍脑袋想出来的，比如另外一个课题组每天就是构建团簇，只要能量为正，马上舍弃此结构。
赞一下(10人) 回复此楼
20楼2009-04-23 22:32:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ron_hai 的主题更新
202/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +11 Mr. Z 2026-04-05 12/600 2026-04-08 16:15 by luoyongfeng
[考研] 材料调剂 +19 一样YWY 2026-04-02 20/1000 2026-04-08 11:02 by 不吃魚的貓
[考研] 304求调剂（085602，过四级，一志愿985） +25 化工人999 2026-04-04 26/1300 2026-04-07 22:06 by hemengdong
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +9 what张 2026-04-01 11/550 2026-04-07 09:45 by momo皓
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 一志愿太原理工大学计算机技术专硕348，求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 385分 生物学（071000）求调剂 +12 qf626 2026-04-01 12/600 2026-04-06 17:36 by 52305043001
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 348求调剂 +3 车厘子zzz 2026-04-05 3/150 2026-04-05 20:30 by 啵啵啵0119
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 材料调剂 +10 懒羊羊轻置玉臀 2026-04-02 11/550 2026-04-04 21:56 by laoshidan
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 312求调剂 +4 赊月色 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[考研] 261求B区调剂 +5 明仔· 2026-04-01 7/350 2026-04-02 11:17 by 邹尉尉
信息提示
关闭
请填处理意见
关闭