计算化学公社

标题: 求助VASP计算Li+能量，能量很奇怪 [打印本页]

作者
Author: 炫柠檬 时间: 2025-4-25 23:04
标题: 求助VASP计算Li+能量，能量很奇怪
各位老师和同学好，
计算目的：用vasp求碱金属离子能量
设置：INCAR (, 下载次数 Times of downloads: 1)
KPOINTS (, 下载次数 Times of downloads: 0)
POTCAR：Li_sv

目前结果： (, 下载次数 Times of downloads: 1)
问题：锂离子能量很奇怪，和其余碱金属离子有很大差别，多次改动设置也没有解决
(, 下载次数 Times of downloads: 0)
非常感谢大家提供给宝贵意见

作者
Author: ultramanm87 时间: 2025-4-26 23:58
1. 比较绝对能量没有意义
2. 周期性体系计算静电荷不为0的能量不可靠，需要各种校正，不是偶极校正能解决的。
VASP有SCPC (https://github.com/aradi/SCPC-Method/tree/main)。
qe+environ可以在某些方向上实现非周期性，有PCC，MT等校正。
cp2k也可以实现非周期性。
如果以上校正方法都不可用，可以尝试用外推（https://pubs.rsc.org/en/content/articlelanding/2018/cp/c7cp08209f）。
论坛里有很多讨论，如：http://bbs.keinsci.com/thread-9422-1-1.html，http://bbs.keinsci.com/thread-28546-1-1.html
3. 计算离子需要考虑是否需要溶剂

作者
Author: 卡开发发 时间: 2025-4-27 04:11

ultramanm87 发表于 2025-4-26 23:58
1. 比较绝对能量没有意义
2. 周期性体系计算静电荷不为0的能量不可靠，需要各种校正，不是偶极校正能解决 ...

2、LDIPOL=.True.不完全等于只开偶极校正。对于立方盒子的带电原子分子而言，但需要严格监测能量E和盒子边长L收敛情况，其他情况不适用。

作者
Author: ultramanm87 时间: 2025-4-27 13:18

卡开发发发表于 2025-4-27 04:11
2、LDIPOL=.True.不完全等于只开偶极校正。对于立方盒子的带电原子分子而言，但需要严格监测能量E和盒子 ...

got it! 正好想在这里问下卡开发发大佬，在对带电体系观察E_tot vs 1/L 时，每次计算时要不要带上 IDIPOL、 LDIPOL 这些校正？我发现带上这些校正后，E_tot和1/L之间的线性关系会变差。而且我看文章里（https://pubs.rsc.org/en/content/articlelanding/2018/cp/c7cp08209f）没有说明在全周期时使用了偶极校正

作者
Author: 卡开发发 时间: 2025-4-27 14:34

ultramanm87 发表于 2025-4-27 13:18
got it! 正好想在这里问下卡开发发大佬，在对带电体系观察E_tot vs 1/L 时，每次计算时要不要带上 IDIPOL ...

开了LDIPOL和IDIPOL之后E_tot和1/L的关系会被破坏：

Charged systems:
Quadrupole corrections are only correct for cubic supercells (this means that the calculated 1/L3 corrections are wrong for charged supercells if the supercell is non-cubic). In addition, we have found empirically that for charged systems with excess electrons (NELECT>NELECTneutral) more reliable results can be obtained if the energy after correction of the linear error (1/L) is plotted against 1/L3 to extrapolate results manually for L→∞. This is due to the uncertainties in extracting the quadrupole moment of systems with excess electrons.
Potential corrections are only possible for orthorhombic cells (at least the direction in which the potential is corrected must be orthogonal to the other two directions).

接下来也得看具体情况：
（1）开放边界体系（IDIPOL=4），可以把模型放到一个立方晶胞，剩余相互作用的高阶项收敛随盒子尺寸应该会较快，谨慎测试即可。
（2）其他情况需要慎重，使用修正或许不是好的选择：

Specifically, for super cells using periodic boundary conditions, it is often not possible to determine the dipole at the defect site accurately, whereas for 0D systems (i.e. atoms and molecules) the dipole can be determined accurately.

作者
Author: vasp冲冲 时间: 2025-4-27 21:42

ultramanm87 发表于 2025-4-26 23:58
1. 比较绝对能量没有意义
2. 周期性体系计算静电荷不为0的能量不可靠，需要各种校正，不是偶极校正能解决 ...

这个好

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