计算化学公社
标题: 求助：关于VASPsol 使用细节问题 [打印本页]
作者Author: Tanmin 时间: 2021-5-25 23:07
标题: 求助：关于VASPsol 使用细节问题
本帖最后由 Tanmin 于 2021-5-25 23:07 编辑
a. 在使用VASPsol时，设置关键词LSOL = .TRUE. ，默认是相对介电常数为80的水作为隐式溶剂，此处若研究特定有机溶剂体系，是否需将关键词EB_k改为对应溶剂的介电常数值? 如果不知所研究体系的介电常数，能否采用该默认值，讨论吸附强弱问题？
b. 此外关键词TAU是关于表面张力，将TAU=0 是忽视cavitation energycontribution (暂时还不太清楚该含义），那是不是如果有特定溶剂，需要将TAU设为该溶剂的表面张力值？
c. VASPsol 手册中有这么一段内容：
Furthermore, if you shift the electrostatic potential to align thepotential in the electrolyte region to zero, or any other value, you also needto add a correction to the energy of your system. This correction is given by QV, where Q is the net charge of thesimulation cell and V the shift in reference potential, e.g.,V = FERMI_SHIFT.
其中Q是否无需设置？比如根据POTCAR中各元素价电子数，我算出体系总价电子20，将NELECT＝22 以考虑-2 价的阴离子在表面的吸附。是否此处自动Q即设为 2 ，而无需在INCAR 单独设置Q=2？
d. 最后针对带电离子的吸附问题，在VASP常规优化任务基础上，是否加以下关键词即可：
NELECT = ** # 说明带电情况
LSOL = .TRUE. #开启溶剂化
V = FERMI_SHIFT # 开启能量校正，此处见论坛有群友谈及需要上下表面对称，如果研究对象是石墨烯，还需要考察其他问题么？
另外一个小的想法，不知是否正确，请各位指点：
假如研究无机熔盐体系，以NaF-AlF3为例，在1000℃下， 假设体系中形成 AlF5[2-], 熔盐整体其实为中性，AlF5[2-] 之所以为-2价是因为从Na那里得了2个电子，如果AlF5周边有其他物质与它存在电荷转移，它就无需从Na那里得电子，那么我是不是可以将AlF5[2-]吸附问题简化为AlF5的吸附问题，进而比较AlF3, AlF4, AlF5 在某材料表面吸附差异？
作者Author: 卡开发发 时间: 2021-5-26 09:27
a、默认是水的介电常数，按道理应该修改为有机溶剂的介电常数。
b、TAU应该是非静电贡献部分的参数，这个可能需要看VASPsol原始文献的定义。
c、Q应该不用设置，印象中Q似乎不是个参数，而是需要手动根据给出来的FERMI_SHIFT和你体系的过剩电荷值去计算。FERMI_SHIFT是计算完成会输出的结果，而不是设定值。
d、石墨烯的话其实搞上下对称的表面似乎不大容易，上下都吸附好像也不太符合实际的情况。
我的看法是，如果AlF5[2-]的两个电子来自于Na，可以用Na作为抗衡离子；如果其他离子作为电子的给体也可以考虑用其他离子，这样避免体系带电荷，而不是简单只是放AlF5。
作者Author: qins 时间: 2021-10-20 19:19
👍
作者Author: 阿司匹林泡腾片 时间: 2022-4-29 18:56
本帖最后由 阿司匹林泡腾片 于 2023-7-19 15:01 编辑
LRHOB=TRUE打开会产生一个RHOION文件类似CHGCAR，可以用VESTA打开看
作者Author: qmlearner 时间: 2023-11-9 15:31
本帖最后由 qmlearner 于 2023-11-9 15:40 编辑
请问卡老师，楼主这个体系，用Na作为抗衡离子的话，实际操作中就是在模型中放一个AlF5，一个Na在材料表面，然后按常规优化吸附模型吗？算AlF5吸附能时要怎么做呢？是Eads=E(surf+AlF5+Na)-Esurf-E(AlF5+Na)吗？（E(AlF5+Na）通过直接删除表面单点计算获得吗？）
另外如果加了抗衡离子的话，可以不考虑溶剂了吗？
作者Author: 卡开发发 时间: 2023-11-21 01:16
如果只是得到吸附构型，这没问题，如果要得到单独的AlF5这种办法其实也行不通。
抗衡离子只是维持电中性，溶剂化是另一回事。
作者Author: qmlearner 时间: 2023-11-22 10:39
谢谢卡老师的回复。最近需要计算一个阳离子在金属氧化物某个晶面的吸附能，看了很多帖子，感觉不太好解决啊。尝试用簇模型，但是金属氧化物表面结构的簇模型实在是门大学问，一时半会没解决。请问算离子在表面吸附能的问题有没有什么好的处理方法呢？
作者Author: 卡开发发 时间: 2023-11-22 16:54
对vasp来说可能有一些修正方案，比方SCPC，我没有详细测试过，到底多可靠我也很难说，按说存在溶剂的情况问题会少很多。
作者Author: qmlearner 时间: 2023-11-23 10:10
谢谢老师，我去学习一下。另外，这类问题用CP2K是否可靠呢？我的体系是一个一价的金属离子在氧化物表面的吸附，CHARGE设置为1，真空中优化（加了溶剂收敛不了
），然后计算吸附能，这样的话投文章有没有问题呢？
作者Author: 卡开发发 时间: 2023-11-23 10:49
我不太确定。
作者Author: qmlearner 时间: 2023-11-23 15:40
谢谢卡老师啦。保险起见我还是用vasp在模型中加一个抗衡离子，然后加隐式溶剂计算了：）
作者Author: Gris 时间: 2024-12-4 12:47
您好，请问您最后怎么解决的问题？带电离子的吸附能计算也是用的vasp算带抗衡离子的体系吗？您的计算公式是Eads=E(surf+AlF5+Na)-Esurf-E(AlF5+Na)吗？
作者Author: lin2206 时间: 2025-7-16 16:09
您好，请问您最后会了吗？
| 欢迎光临 计算化学公社 (http://ccc.keinsci.com/) |
Powered by Discuz! X3.3 |