计算化学公社

标题: 请问有人用VASPsol的时候算过溶剂电荷密度嘛，怎么用RHOB计算引入的离子数 [打印本页]

作者
Author: calamity 时间: 2021-3-24 19:34
标题: 请问有人用VASPsol的时候算过溶剂电荷密度嘛，怎么用RHOB计算引入的离子数
VASPsol会自动引入补偿电荷使体系保持电中性，可以通过开启LRHOB = .TRUE. 和LRHOION = .TRUE.来输出溶剂电荷密度文件RHOB和溶剂离子浓度RHOION，输出文件和CHGCAR类似，但是怎么样从RHOB积分从而获得体系中总共引入了多少补偿电荷呢？就是这篇文章中的做法（DOI:10.1021/acs.jpclett.6b00358 J. Phys. Chem. Lett. 2016, 7, 1471−1477），作者利用修改体系电子数来获得恒定电极电势，并利用VASPsol保持体系电中性，但是在计算能量时会有阳离子电荷数N+和阴离子电荷数N-，这个怎么积分得到的？请问有人做过相似的研究吗？ (, 下载次数 Times of downloads: 75)

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

作者
Author: calamity 时间: 2021-3-25 14:16

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

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

作者
Author: 卡开发发 时间: 2021-3-25 15:37

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

对。

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

作者
Author: 卡开发发 时间: 2021-12-19 01:47

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

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

作者
Author: 飞翔的猪 时间: 2021-12-22 15:58

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

谢谢老师

作者
Author: Jack 时间: 2024-3-9 17:44

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

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

作者
Author: 卡开发发 时间: 2024-3-10 05:56

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，和看电子密度差/差分电荷密度一样的操作），这里面可能存在的坑点是我不确定单位和系数上会需要额外做什么操作。

作者
Author: Jack 时间: 2024-3-10 17:05

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

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

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