大家用VASP做电催化的恒电势计算都用哪些程序

dingniu2

通过查阅发现VASP做电催化的恒电势计算有用JDFTx（https://jdftx.org/），也有用CP-VASP这种VASP补丁重新编译VASP的（https://sites.utexas.edu/yuanyue-liu/codes/cp-vasp/），还有用新的vaspsol++做恒电势（https://github.com/VASPsol/VASPsol），还有ASE接口的FCP-vasp-ase（https://github.com/hellozhaoming/FCP-vasp-ase）。请问大家两个问题：1. 目前流行的恒电势计算只有这4种方式吗？还有其它的吗？2.这4种做恒电势有什么区别？大家都选的是哪一种？

nullzy

现在vaspsol++也能做恒电势了，CP-VASP用过还可以，功能非常全面，JDFTx很早之前测试过，奈何速度实在是太慢了。如果非要用VASP的话清华大学肖海课题组开发的FCP算法也可以算恒电势，不用VASP的话可以用GPAW软件的SJM模块算恒电势，国产软件PWmat内置了高国平老师的巨正则固定电势方法也能算恒电势。

nullzy

还有Nørskov组的电荷外推法跟晶胞外推法，前者必须得加溶剂算过渡态，后者计算量太大很少有人用了。还有只需要装vaspsol就可以用的double-reference method，用的人也越来越多了。

dingniu2

nullzy 发表于 2024-12-11 14:36
还有Nørskov组的电荷外推法跟晶胞外推法，前者必须得加溶剂算过渡态，后者计算量太大很少有人用了。还 ...

感谢回答，请问：这些做恒电势的方法中是不是CP-VASP是目前最流行的，认可度比较高的？谢谢！

nullzy

dingniu2 发表于 2024-12-11 18:05
感谢回答，请问：这些做恒电势的方法中是不是CP-VASP是目前最流行的，认可度比较高的？谢谢！

CP-VASP很有前景，因为也是最近才开放大家使用的，所以目前用的人还不是很多，但我不确定能不能跑恒电势CI-NEB，毕竟他们算的能垒都是AIMD增强采样得到的。要论用户数量的话，电荷外推法用的非常多，毕竟是老方法了，而且不需要额外装插件。

dingniu2

nullzy 发表于 2024-12-12 01:13
CP-VASP很有前景，因为也是最近才开放大家使用的，所以目前用的人还不是很多，但我不确定能不能跑恒电势C ...

谢谢回答，请问电荷外推法是不是一定要加显式溶剂算过渡态？如果是这样，算的会不会很慢，且不容易收敛？

nullzy

dingniu2 发表于 2024-12-12 09:17
谢谢回答，请问电荷外推法是不是一定要加显式溶剂算过渡态？如果是这样，算的会不会很慢，且不容易收敛？

是的，必须要加，实测确实很慢，你可以多看看肖建平老师的文章，他们以前经常用这种方法，但最近他们也开发了基于电场调控的恒电势方法，也是要给VASP打补丁，不过貌似代码没有开源。

dingniu2

nullzy 发表于 2024-12-12 15:37
是的，必须要加，实测确实很慢，你可以多看看肖建平老师的文章，他们以前经常用这种方法，但最近他们也开 ...

谢谢回答，谷歌学术查了下用CP-VASP计算的文献，发现只有5篇左右，是因为这个代码刚开源吗？还是用的人不多？谢谢！

Trinity

nullzy 发表于 2024-12-11 14:36
还有Nørskov组的电荷外推法跟晶胞外推法，前者必须得加溶剂算过渡态，后者计算量太大很少有人用了。还 ...

您好，请问这里提到的double-reference method是Phys. Rev. B 73, 165402这个工作吗？或者双参考指的是隐式溶剂区域的泊松方程要求静电势为0和VASP要求盒子整体静电势为0两个静电势参考点之间的差值？谢谢！

机智的咸鱼

dingniu2 发表于 2024-12-13 10:49
谢谢回答，谷歌学术查了下用CP-VASP计算的文献，发现只有5篇左右，是因为这个代码刚开源吗？还是用的人不 ...

恒电势下跑AIMD，收敛困难计算速度也很慢，不是所有课题组的计算资源都能受得了，这其中还要算上学生基本概念不明确导致乱算浪费机时的成本。有隐式模型的时候涉及电荷密度精度，ENCUT原则上600 eV起步，现在大家都用这么三四百跑我也是没什么话说了。如果你只是要算个恒电势加台阶图，JDFTx足矣，而且里面先进的隐式溶剂模型更多。VASPSOL在++之前的版本，说实话，那个隐式溶剂模型完全就是个残废。