优化收敛问题

卡开发发

lujinhui396 发表于 2014-11-30 20:44
之所以采用这个去掉二价钴的模型，是因为，有关于四氧化三铁的100 表面通过低能电子衍射，发现表面是有铁 ...

如果是考虑缺陷这个问题我就再找找原因了，你那个结构拿过来我算了一下电子步不怎么容易收敛，所以猜测结构可能不太好。我再找找原因吧。

lujinhui396

卡开发发 发表于 2014-11-30 21:11
如果是考虑缺陷这个问题我就再找找原因了，你那个结构拿过来我算了一下电子步不怎么容易收敛，所以猜测结 ...

谢谢

卡开发发

lujinhui396 发表于 2014-12-1 09:00
谢谢

电子依旧不收敛，试试看只有一面有缺陷的情况吧。

lujinhui396

卡开发发 发表于 2014-12-1 14:25
电子依旧不收敛，试试看只有一面有缺陷的情况吧。

首先，十分感谢您的热心帮助。电子不收敛，改变算法也是这样吗？另外还用问题向您请教
(1) 离子弛豫的时候，如果设的NELM=60，但是算完60个电子步之后，依然没达到电子的收敛标准，就自动进行下一步的离子弛豫，且最后离子弛豫，达到收敛标准，请问这样的结果，可信吗？
（2）请问您一般怎么优化结构不太合理的构型呢？比如表面，缺陷，吸附分子的模型？我一般都是分别优化。可是当我的力的收敛判据达到-0.02时，和初始的粗的力的收敛判据优化收敛得到的构型，相差很大。
（3）最后问您个私人问题，您是在天津吗？我在天津的学校上学，总是感觉有时文字不能说清我的意思，如果您是天津的，我想当面向您请教。

卡开发发

恩，切换了不少代数方法都不行，都很难在60步以内收敛，原则上就是结构有问题了，因为现在您的晶格处理的情况我也不是很了解，加上考虑了+U，导致不收敛的因素有点复杂。不过以前我自己算CASTEP的时候有时处理两面的缺陷离子步也经常不容易收敛，电子步没收敛很有可能是初始结构有问题，我做了些微调无济于事。我没做什么参数的检查，我默认您是做过测试的。
（1）最好检查一下末尾那几步的电子步是否也是收敛的，否则得到的电子态就是错误的，力得到的肯定也是错误的，这样结果就没有意义了；如果您的后面n多步电子步都是收敛的，只有前面少数几步和预设值稍有点差异，这种情况的结果应该是可以接受的；
（2）一般也是分别优化，适当降低一点精度，但是并不建议降低太多，印象中手册上提供了电子步最低精度不应小于1E-4，默认相应的EDIFFG=EDIFF*10，这意味着至少是E-3的级别，我觉得收敛精度也许还应该进行提升才有可靠的结果；
（3）不在，目前人在上海。不过可以站内私信或是其他的方式。不过说起来帮别人检查离子步(scf)不收敛有时候确实是个相当头疼的问题，有时候是电子步导致的就得反过来检查INCAR。

lujinhui396

卡开发发 发表于 2014-12-2 16:45
恩，切换了不少代数方法都不行，都很难在60步以内收敛，原则上就是结构有问题了，因为现在您的晶格处理的情 ...

参数测试，KPOINTS采用的是文献中默认的，因为我做K点测试时，一直也是电子步不收敛。

卡开发发

lujinhui396 发表于 2014-12-2 20:12
参数测试，KPOINTS采用的是文献中默认的，因为我做K点测试时，一直也是电子步不收敛。

不妨把文献发来，我们来找找原因为啥电子步和离子步这么难收敛。

lujinhui396

这是两个模型，一个是带有二价钴缺陷的，一个是没有的

lujinhui396

卡开发发 发表于 2014-12-2 22:05
不妨把文献发来，我们来找找原因为啥电子步和离子步这么难收敛。

电子步这么难收敛，是不是与ISMEAR有关呢？

卡开发发

lujinhui396 发表于 2014-12-3 11:36
电子步这么难收敛，是不是与ISMEAR有关呢？

一般而言应该影响没这么大，我测试过ismear=0还是不行的。

卡开发发

lujinhui396 发表于 2014-12-3 09:21
这是两个模型，一个是带有二价钴缺陷的，一个是没有的

恩，这两篇文章我先看看。这几天有点感冒，在休息。有时间的话我把你的INCAR参数测试一下，重新确定一下晶格常数等，表面我在按照文献做一下看看，有消息我再找你私聊，等最终有结果我们再把问题解决的推演过程放到帖子来。

lujinhui396

卡开发发 发表于 2014-12-4 10:26
恩，这两篇文章我先看看。这几天有点感冒，在休息。有时间的话我把你的INCAR参数测试一下，重新确定一下 ...

非常感谢~~~多喝水，多休息，早点好起来。。。。

viva

您好，能请教一下您的Co3O4晶胞优化的INCAR文件吗？
我也是初学VASP，但是自己的INCAR算出来没有离子步出来，并且想问下您算Co3O4需要注意一些什么，我也是看了您的帖子来的。谢谢啦

流年2022

请问这个电子步收敛是看的dE嘛？

卡开发发

流年2022 发表于 2021-5-26 10:33
请问这个电子步收敛是看的dE嘛？

dE和d eps两个参量。参考手册：

The relaxation of the electronic degrees of freedom will be stopped if the total (free) energy change and the band structure energy change (’change of eigenvalues’) between two steps are both smaller than EDIFF

以及

The columns have the following meaning: Column N is counter for the the electronic iteration step, E is the current free energy, dE the change of the free energy between two steps, and d eps the change of the band-structure energy.