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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 为什么在dos计算时，会出现WARNING:stress and forces are not correct
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
101/1返回列表
查看: 5665  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gougou531

禁虫 (正式写手)
本帖内容被屏蔽

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-12-05 15:06:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

gougou531(金币+1): 哦。这样啊，如果用静态计算的poscar计算就应该不会有这个提示了吧，呵呵谢谢 2011-12-05 15:24:18
这个没有关系的
用优化后的CONTCAR算band或dos都会出现这个提示
赞一下 回复此楼
2楼2011-12-05 15:21:19
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)
本帖内容被屏蔽
3楼2011-12-05 15:29:16
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★
gougou531(金币+2): 谢谢，我还有个问题就是原来晶胞是a=b=c，优化完a变长了，b和c相等但是变短了，这是怎么回事？ 2011-12-05 19:28:35
franch(金币+2): 谢谢回帖交流 2011-12-05 19:54:05
1.“为什么计算的TOTAL-FORCE这一栏全是0，并且发现原子的坐标没有变化”   
这个应该是原子坐标全部由对称性决定的。然后vasp默认是保持对称性的（ISYM/=0），所以原子位置不发生变化。
2.我一般直接ISIF=3优化(PREC=high,EDIFF=10-5,EDIFFG=-0.01),然后CONTCAR拷成POSCAR重新优化一遍（一般一步就结束了）
3.体积能量拟合有一个murn.f的小程序，小木虫内可以搜一下
赞一下(1人) 回复此楼
4楼2011-12-05 17:37:43
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ffzx

金虫 (正式写手)

【答案】应助回帖

gougou531(金币+1): 对，用的是GGA势，这和势有什么关系吗？ 2011-12-07 20:28:15
引用回帖:
1楼: Originally posted by gougou531 at 2011-12-05 15:06:50:
如题：在计算dos时，为什么会出现 WARNING:stress and forces are not correct。这是那些原因造成的，计算的是多铁物质。还有就是在优化时，TOTAL-FORCE这项怎么全是0啊，原子位置没有变化啊，这是怎么回事啊？谢 ...

恩。应该没关系的。另外跟你的势有关，你是不是用的GGA势？
赞一下 回复此楼
5楼2011-12-07 20:24:08
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ffzx

金虫 (正式写手)

【答案】应助回帖

gougou531(金币+1): 恩，谢谢了，呵呵，计算这个东西说不清啊，真是。。。。。 2011-12-07 21:37:33
LDA就没这个warning。我测试过了。同样体系，就是势换了一下。所以不用太担心。
赞一下 回复此楼
6楼2011-12-07 20:36:06
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ffzx

金虫 (正式写手)

【答案】应助回帖

又找到一个类似的情况，也是GGA才有这个warning：

http://www.mdbbs.org/viewthread. ... table&tid=25765
赞一下 回复此楼
7楼2011-12-10 11:43:12
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)
引用回帖:
2楼: Originally posted by fzx2008 at 2011-12-05 15:21:19
这个没有关系的
用优化后的CONTCAR算band或dos都会出现这个提示

你好，想请教你个问题。你回复LZ的问题也是我刚遇到的，我想接着问一个，我需要在P4vasp里面导入文件vasprun.xml来画DOS图，但是正常的情况下这个文件应该有几十MB那么大，出现这种情况后变成了几十KB，也就是变得非常小，使得根本画不出DOS图，请问大神这是什么情况
赞一下 回复此楼
8楼2014-04-28 15:55:09
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
引用回帖:
5楼: Originally posted by ffzx at 2011-12-07 20:24:08
恩。应该没关系的。另外跟你的势有关，你是不是用的GGA势？...

请问在计算的时候还出现了WARNING: random wavefunctions but no delay for mixing, default for NELMDL
怎么改，有没有影响啊
赞一下 回复此楼
希望有一天成为大神
9楼2015-05-10 11:02:14
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

因为你用了ICHARG=11吧，这个是进行non self consistent的计算，所以肯定会出现这个提示的啊。
如果你怕的话，在INCAR可以注释掉ICHARG那一行，然后比较ICHARG=11的差距，你会发现其实差距很小的。但是运算的electronic steps大大减少了，因此是一种高效的计算方法。然后有关于NELMDL的问题，我们知道DFT计算式self consistent的计算，那么一开始的WAVE FUNCTION应该是随即猜测的，所以可能猜测不是很合理，但如果设置NELMDL的话，一开始进行一定步数的non self consistent计算，然后得到一个相对合理的WAVE FUNCTION之后，再用这个进行self consistent计算的初始WAVE FUNCTION，就会使运算的步数大大减少，效率也就提高了

推荐LZ多看手册，然后自己可以比较试验一下

ICHARG:
For ICHARG=11 and density functional theory calculations, all k-points become essentially independent, because the charge density and the potential are kept fixed. If necessary it is even possible, to split up the k-points and calculating the eigenvalues individually for each k-point, although with present computing platforms this is rarely required.

ref：
ICHARG: http://cms.mpi.univie.ac.at/vasp ... e_calculations.html
NELMDL: http://cms.mpi.univie.ac.at/vasp/guide/node104.html
赞一下 回复此楼
10楼2015-05-16 04:10:43
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gougou531 的主题更新
101/1返回列表
信息提示
关闭
请填处理意见
关闭