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

好的，谢谢老师的解答

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

vaspwiki:

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

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先同样泛函优化好，再贴到晶胞去。

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

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

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

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

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

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

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

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

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

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

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

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

你有什么想不开的不用 Gaussian ORCA 算偏要用 VASP 算？