卡开发发 发表于 2024-3-10 05:56
这个问题可能vaspsol的开发组没提过，所以我去看了一眼vaspsol的源码，如果下面的做法解决不了问题建议发 ...

好的，谢谢卡老师详细的解答

Jack 发表于 2024-3-9 17:44
请问卡老师，怎么在VASPsol输出2维的介电常数分布，谢谢！

这个问题可能vaspsol的开发组没提过，所以我去看了一眼vaspsol的源码，如果下面的做法解决不了问题建议发邮件问一下vaspsol的开发组：
1、5.4.1之后的版本应该只依赖新的solvation.F，在solvation.F开头有#include "symbol.inc"的行。
2、早期版本可能不在solvation.F，是三个文件pot_k.F、pot_lpcm.F和solvation.F，然后其中我们看pot_lpcm.F，开头也有#include "symbol.inc"。
然后在#include "symbol.inc"前后增加#define debugsol，重新编译后再进行运行，这样应该会产生一个EPS的文件，其格式和CHGCAR应当相同。二维的话你可能需要把这个文件改改名字让其他可视化程序识别这个文件的类型，然后拖到可视化程序建立切片显示（例如vesta，和看电子密度差/差分电荷密度一样的操作），这里面可能存在的坑点是我不确定单位和系数上会需要额外做什么操作。

卡开发发 发表于 2021-3-25 09:35
也许可以尝试的方法是这样，有CHGCAR的格式意味着可以使用类似于Bader程序进行积分来看总的电子数，-ref建 ...

请问卡老师，怎么在VASPsol输出2维的介电常数分布，谢谢！

卡开发发 发表于 2021-12-19 01:47
换个版本试试，实在不行也只好手工把格式调整成和chgcar一样。

谢谢老师

飞翔的猪 发表于 2021-12-17 13:56
您好，请问INCAR中设置了LRHOB = .TRUE.  为什么输出文件RHOB无法用vesta打开显示出溶剂电荷密度

换个版本试试，实在不行也只好手工把格式调整成和chgcar一样。

您好，请问INCAR中设置了LRHOB = .TRUE.  为什么输出文件RHOB无法用vesta打开显示出溶剂电荷密度

calamity 发表于 2021-3-25 14:16
请问用全电子做-ref指的是加LAECHG=.TURE.算单点，然后把两个电荷密度加起来得到的CHGCAR_sum吗？

对。

卡开发发 发表于 2021-3-25 09:35
也许可以尝试的方法是这样，有CHGCAR的格式意味着可以使用类似于Bader程序进行积分来看总的电子数，-ref建 ...

请问用全电子做-ref指的是加LAECHG=.TURE.算单点，然后把两个电荷密度加起来得到的CHGCAR_sum吗？

也许可以尝试的方法是这样，有CHGCAR的格式意味着可以使用类似于Bader程序进行积分来看总的电子数，-ref建议上还是用全电子的（我不确定会不会卡在盆划分上，所以保险起见最好用全电子的）。虽然自己写个处理程序也成，不过没必要去折腾。如果有其他方法也欢迎讨论。