本帖最后由 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 在某材料表面吸附差异?
|