计算化学公社

标题: 请问LDA加U无法收敛该如何调整 [打印本页]

作者
Author: jzjy1994 时间: 2017-9-4 09:16
标题: 请问LDA加U无法收敛该如何调整
目前已经换成了模守恒赝势，空带也调整过了，然后能量变化值一直在波动。用超软也无法收敛，之前不加u的情况下很快就能收敛，请问一下各位老师有什么好的解决方法

作者
Author: 卡开发发 时间: 2017-9-4 12:17
能否将CASTEP的输出文件（*.castep）上传一下？

一般castep文件内记录着计算过程涉及到的内容，包括你采用的计算参数，结构，采用的赝势等，还有迭代以及构型优化收敛的情况。有这些内容方便论坛坛友来更快速、准确地诊断你的问题到底发生在何处。尤其是像现在讨论的SCF难收敛或者构型优化难收敛等问题，很难直接了当诊断为某个单一因素造成。大家可以根据你所提供的文件来说明如何分析你的输出文件，或者进行重复性的测试便于排查最为可能的因素，以及一些其他的经验。

就比如SCF难收敛我们可以先参考http://sobereva.com/61帖子尝试进行解决，其中的2下面的某些点对于平面波计算固体都具有移植性：
(1)对于平面波方法而言增加动能截断和密度的FFT网格等也是很常规的手段。
(13)的移植方式对于三维周期性的固体而言可以适当去调整晶格参数，一般+U可能会造成平衡晶格结构偏大等，还有结构不合理的情况，比如我见过有人磷酸铁系统切表面的时候把磷酸根的磷-氧键切断，导致SCF很难做。对于固体或表面的计算。
(14)类似的做法就是体系内的原子猜测一个更合理的磁矩，针对DFT+U而言这些对磁矩有贡献的原子更加局域，更合理的磁矩SCF收敛将会更快。
并不排除(17)也有可能会遇到，特定的泛函和赝势确实有可能会有潜在的问题未被发现。

在上述参数合理的情况下，3中提到3.1（也就是MS下的density mixing）和3.2（就是MS下的DIIS）是MS较为常用的做法，但调整起来确实需要一些实际经验。一般而言，DIIS都是开启的，调整3.1比较会有效，对于d电子在Fermi能级附近有贡献且局域的体系可以适当降低权重“density mixing-charge”会比较奏效，少数体系也可能需要调整“density mixing-spin”。

作者
Author: jzjy1994 时间: 2017-9-5 20:59

卡开发发发表于 2017-9-4 12:17
能否将CASTEP的输出文件（*.castep）上传一下？

一般castep文件内记录着计算过程涉及到的内容，包括你采 ...

谢谢老师您花费时间回复。任务失败之后我就没下载输出文件（=0=下次一定保存）我先按照您和sob老师帖子中的方式做一些调整，如果还是有问题的话，我将输出文件上传。（另外关于k点-，-我因为想后期做收敛性测试，所以单点能计算用的k点确实不多。）非常感谢啦

作者
Author: 卡开发发 时间: 2017-9-5 22:58

jzjy1994 发表于 2017-9-5 20:59
谢谢老师您花费时间回复。任务失败之后我就没下载输出文件（=0=下次一定保存）我先按照您和sob老师帖子中 ...

对，下次注意计算失败的结果不要马上删除，应该对出错的信息做一些分析，有时候很难说你会从那些东西里面发现什么别的。

欢迎光临计算化学公社 (http://ccc.keinsci.com/)