计算化学公社

标题: 求助：VASP杂化泛函优化结构电子步不收敛 [打印本页]

作者
Author: wangy0822 时间: 2023-2-22 20:36
标题: 求助：VASP杂化泛函优化结构电子步不收敛
大家好，我用B3LYP杂化泛函优化并计算频率，结构如下图所示，这个结构原本是单胞里有4个分子，然后我删除了3个，留了一个，目前优化非常难收敛，能量的收敛标准已经降到10E-4了，下面是我做优化的INCAR，电子步到了设置的200步上限都没有收敛，难道还要增加电子步上限（我觉着没太大意义了吧），而且优化之后还要计算频率，是不是会更难算。是不是这个泛函优化很困难，或者说大家有没有好的方法可以帮助收敛。谢谢大家

作者
Author: 北大-陶豫 时间: 2023-2-23 10:44
你有什么想不开的不用 Gaussian ORCA 算偏要用 VASP 算？

作者
Author: wangy0822 时间: 2023-2-23 11:10

北大-陶豫发表于 2023-2-23 10:44
你有什么想不开的不用 Gaussian ORCA 算偏要用 VASP 算？

谢谢您的解答，这是审稿人要求用VASP计算的，目的是和Gaussian做一下对比，但是优化确实很难，我想请教一下能不能行，如果不行我就给审稿人说明情况

作者
Author: jiangning198511 时间: 2023-2-23 11:20

wangy0822 发表于 2023-2-23 11:10
谢谢您的解答，这是审稿人要求用VASP计算的，目的是和Gaussian做一下对比，但是优化确实很难，我想请教一 ...

审稿人的建议完全不合理，即使要对比，也应该用类似的软件如ORCA

作者
Author: wangy0822 时间: 2023-2-23 14:25

jiangning198511 发表于 2023-2-23 11:20
审稿人的建议完全不合理，即使要对比，也应该用类似的软件如ORCA

谢谢您的解答，我们的文章主要内容是证明周期性边界条件在计算分子晶体太赫兹光谱中的重要性。构建了不同大小的团簇模型（Gaussian/B3LYP）以及周期性结构模型（VASP/opt-B88，opt-B86，opt-PBE），结果分别与实验谱对比。然后审稿人就提出了下面的一个问题。

comment：
虽然我同意周期性模型优于团簇模型，但我也认为，随着团簇规模的增加，越来越多的非局部相互作用被添加到模型中，其结果应该在某种程度上收敛于周期性研究。
因此，我建议使用两种模型的相同函数（VASP中的B3LYP或高斯中的opt-B88等，VASP中也有DFT-D3）重复/增加团簇研究的某些部分，以进行直接比较。

作者
Author: 北大-陶豫 时间: 2023-2-23 15:30

wangy0822 发表于 2023-2-23 14:25
谢谢您的解答，我们的文章主要内容是证明周期性边界条件在计算分子晶体太赫兹光谱中的重要性。构建了不同 ...

实在要比较的话，咱可以用个纯泛函吧，这样计算量小点。比如都用PBE（高斯里写作PBEPBE）

作者
Author: kubingxia23 时间: 2023-2-23 15:53
ALGO = Damped 或者 ALGO = All 试一下呢

vaspwiki:
The blocked-Davidson algorithm (ALGO=Normal) is generally rather slow, and in many cases the Pulay mixer will be unable to determine the proper ground-state. We hence recommend to select the blocked-Davidson algorithm only in combination with straight mixing or a Kerker like mixing

作者
Author: wangy0822 时间: 2023-2-23 15:54

北大-陶豫发表于 2023-2-23 15:30
实在要比较的话，咱可以用个纯泛函吧，这样计算量小点。比如都用PBE（高斯里写作PBEPBE）

感谢您的回复，如您所说，我在Gaussian和VASP中都用了PBE泛函做了计算，对比结果发现频率非常接近。但是审稿人建议用B3LYP，我就试了一下，发现太难收敛甚至根本收敛不了，直接劝退。。。

作者
Author: 卡开发发 时间: 2023-2-23 16:18
本帖最后由卡开发发于 2023-2-23 16:22 编辑

wangy0822 发表于 2023-2-23 15:54
感谢您的回复，如您所说，我在Gaussian和VASP中都用了PBE泛函做了计算，对比结果发现频率非常接近。但是 ...

1、7楼给的提示是可以的，但我建议可以先不要尝试ALGO=A或者D，大部分情况应该还是ALGO=N对角快一些，提示中建议不要用Pulay混合（DIIS）而换成Kerker混合，我个人经验基本上也是这样，PBC计算杂化泛函Kerker混合可能要稳定一些，你可以试试IMIX=1配合一个合适的AMIX。
2、此前可以纯泛函测试下盒子尺寸是否足够，盒子不必和晶胞的尺寸一样大，合理范围内取小一些节省计算量。动能截断最好也做好充分的测试。
3、结构可以用Gaussian先同样泛函优化好，再贴到晶胞去。

作者
Author: wangy0822 时间: 2023-2-23 16:30

kubingxia23 发表于 2023-2-23 15:53
ALGO = Damped 或者 ALGO = All 试一下呢

vaspwiki:

感谢您的回复，我尝试过了，仍然很难收敛

作者
Author: wangy0822 时间: 2023-2-23 16:31

卡开发发发表于 2023-2-23 16:18
1、7楼给的提示是可以的，但我建议可以先不要尝试ALGO=A或者D，大部分情况应该还是ALGO=N对角快一些，提 ...

好的，谢谢老师的解答

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