| 你好，能详细说一下这个问题怎么解决的吗 |
jitou11 发表于 2020-12-7 10:32 按我理解应该是选择一个合适的参考原子（比如化学环境不变），然后以这个原子的值对两个体系进行参考。 |
jitou11 发表于 2020-12-7 10:18 手册有提到。 |
jitou11 发表于 2020-12-3 10:55 按手册的意思应该是LVTOT或LVHAR会打出来，我没做过测试，你可以试一下，实在不行再看代码本身的情况。 |
jitou11 发表于 2020-12-3 10:06 不排除这样的可能，但VASP不是我常用的程序，要确定的话最好去官方论坛或者看源码确定。 |
| 很简单，测试一下HF分子，看看哪个算得对 |
lonelycvavc 发表于 2020-12-2 09:21 那看来不是巧合，具体绘制或处理的时候我不确定是怎样的。 |
卡开发发 发表于 2020-12-2 09:06 谢谢前辈指点！平面波方法计算静电势能级会发生平移，想要比较不同体系之间的静电势可以用average electrostatic potential先对齐能级位置再比较，这个我有点明白了。那您看我算出来LVHAR的图像结果和Multiwfn计算恰好相反是我算错了吗？感觉应该不是巧合，我也算了几个其他体系，结果都是一样的，图像形状几乎相同但是颜色相反。 |
lonelycvavc 发表于 2020-12-2 08:45 LVHAR就是计算静电势的（不包含Vxc），区别是平面波计算方法在处理静电作用的时候会扔掉G=0的部分，导致能级都会发生一个平移，所以算出来静电势的绝对值不具备参考意义（通常表面系统我们可以参考远离表面的区域作为能量0点来确定静电势的相对值）。当然，你上面提到的average electrostatic potential也可以用来对齐不同体系的能级位置进行比较。 |
卡开发发 发表于 2020-12-1 23:41 谢谢！我采用的是LVHAR，LVTOT做出的图像和Multiwfn相去甚远，但LVHAR的结果和Multiwfn的静电势几乎完全相同（除了颜色）。请问LVHAR是否可以用来计算周期性体系的静电势呢？以及为什么它好像和静电势数值正好相反呢？ |
lonelycvavc 发表于 2020-12-1 17:48 average electrostatic potential不是做这个用的，而是对标能级的时候才用得到。 我不确定你是否采用的是LVHAR而不是LVTOT，后者会把Vxc计入进来。 |
sobereva 发表于 2020-12-1 17:25 感谢社长。VASP手册中将这个式子称为“the average electrostatic potential at each ion”，全文如下：VASP also calculates the average electrostatic potential at each ion. This is done, by placing a test charge with the norm 1, at each ion and calculating然后给出了帖中的式子。 我也认为氧附近的静电势不可能为正值，但从图像来看，两个计算的结果非常相似，几乎就是颜色反了一下，那么VASP给出的结果我能不能添加负号直接使用呢？ |
|
静电势的定义就一种。Multiwfn用的是最标准的静电势的定义 静电势的定义不可能是你最后一张图的样子，那都不是函数了，而是个标量 氧附近的静电势明显不可能是正值，正值说明算错了。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-19 04:05 , Processed in 0.194135 second(s), 31 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites