大家用VASPsol只做静态计算吗

dingniu2

请问，VASP做催化反应自由能图计算时考虑隐式水溶剂效应，会用到VASPsol编译的VASP，VASPsol编译的VASP可以做结构优化（获得后续计算的稳定结构）、静态计算（获得电子能量）和频率计算（获得自由能校正）。官网上的算例只说了用VASPsol做静态计算。请问大家做催化反应自由能图计算时，通常的做法会在哪个阶段开LSOL = .TRUE.考虑溶剂计算？

ultramanm87

只用过cp2k的、还有qe的environ插件，遇到稍微复杂点的geometry（比如slab上摆一个带电吸附质），结构优化收敛就挺困难。不知道VASPsol怎么样。

白皛白

加了隐式溶剂化效应以后电子步的耗时翻倍了，而且还更难收敛，实在算不动。
之前还发现一个问题，就是不考虑隐式溶剂化效应优化出来的结构，如果频率计算的时候考虑了溶剂化效应就有可能存在一个或多个小的虚频（不考虑溶剂化效应的频率计算就不会有虚频），我个人认为这也说明结构优化加不加溶剂化确实有影响吧。

wwjie

按照官方算例上来就行了

dingniu2

白皛白 发表于 2024-3-4 10:00
加了隐式溶剂化效应以后电子步的耗时翻倍了，而且还更难收敛，实在算不动。
之前还发现一个问题，就是 ...

请问，即使结构优化不用vaspsol，而只在静态计算获得电子能量时用vaspsol，那么频率计算获得自由能校正时也不需要结构优化为什么不加vaspsol？谢谢！

白皛白

dingniu2 发表于 2024-3-31 19:25
请问，即使结构优化不用vaspsol，而只在静态计算获得电子能量时用vaspsol，那么频率计算获得自由能校正时 ...

您好，我个人在频率计算时之前也试着用过vaspsol，但是在这之间发现了一些问题，所以之后频率计算的时候就不用vaspsol了。
1，发现如果结构优化没有使用vaspsol，而频率计算使用vaspsol的话，有可能会出现小的虚频（而不使用vaspsol的频率计算则不会）。另外在过渡态的计算中，如果过渡态搜索的时候没有使用vaspsol，而频率计算的时候使用vaspsol的话，有可能出现多个虚频（过渡态计算中出现这个问题的概率似乎很大，而不适用vaspsol的频率计算则只有一个虚频），并且这样的虚频除非结构计算或过渡态搜索也用vaspsol，否则是消除不掉的。我个人对这个问题出现的原因的理解是，可能在真空中计算稳定的结构在添加了溶剂化效应的计算中是不稳定的，这种虚频的出现可能会导致自由能矫正的异常，这是我主要考虑的一个问题。
2，使用或不使用vaspsol计算频率，自由能矫正的值大概差0.01~0.03ev左右（如果原本自由能矫正数值比较小比如只有0.3ev左右的话），感觉误差不大。
所以综合来说，我频率计算也选择不加隐式溶剂化了，不知道这个问题有没有标准答案。如果有的话也希望有老师能指点迷津。

dingniu2

白皛白 发表于 2024-4-1 10:49
您好，我个人在频率计算时之前也试着用过vaspsol，但是在这之间发现了一些问题，所以之后频率计算的时候 ...

谢谢回复！
1. 我觉得你的第一条说的是对的，既然结构优化决定不用VASPsol，那么频率计算也不应该用VASPsol，频率计算要与结构化学同一级别下计算，如果一个加了溶剂化儿另一个不加时结果可能是不合理的。而单点计算是可以与结构优化不同级别，因此可以加VASPsol。
2. 有些情况确实优化时加不了VASPsol，会有算不动或难收敛等问题，这个帖子的投票也看出来用VASPsol做结构优化的人不多。
3. 如果绝大多数体系用VASPsol做自由能校正只比不用差0.05以内就好了，说明不用VASPsol做自由能校正也可以了。但不知道你的0.01~0.03ev左右是特例还是一般情况。

cywu

老师们好，想请教一下，如果用vaspsol计算净电荷为0的体系，那么分析功函数，费米能级，dos时候，还是按照没有vaspsol的方法（例如，用LOCPOT文件算功函数）直接分析吗？还是要把vaspsol输出的“fermi-shift”加在费米能级上？不太清楚在净电荷为零的情况下，这个“fermi-shift”到底有什么用...请老师们指点迷津，多谢各位老师！

din5g

cywu 发表于 2024-6-6 15:45
老师们好，想请教一下，如果用vaspsol计算净电荷为0的体系，那么分析功函数，费米能级，dos时候，还是按照 ...

我也蹲下，我用vaspsol算真空能级，他的图像不是平整的，是齿形，所以我的真空能级用的就是不用vaspsol的，不知道对不对