| 查看: 7644 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者lm20220613将赠送您 5 个金币 | ||
[求助]
二氧化铈功函数计算求助 已有1人参与
|
||
请教大家：我用VASP计算二氧化铈（200）晶面的功函数，表面静电势结果如下，波动峰的高度不一致，这是为啥呢？还有，我加了偶极矫正，但是真空能级的位置还是向下倾斜。我附上了我的INCAR和晶胞结构，烦请大家帮我看看，非常感谢大家 。   |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有9人回复
-
放榜前的不淡定
已经有14人回复
-
建议基金发布提前给出明确的时间点
已经有13人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
HER12025
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 134.8
- 散金: 35
- 帖子: 206
- 在线: 112.1小时
- 虫号: 36722173
- 注册: 2025-12-26
- 专业: 高分子材料结构与性能
2楼2026-01-06 16:38:11
【答案】应助回帖
|
Your slab model contains roughly 5 atomic planes. For cerium dioxide ($\text{CeO}_2$), this is quite thin. Without a thick enough "bulk core," the structural relaxations from both surfaces propagate through the entire slab. The interlayer distances and planar charge densities vary from layer to layer, preventing the electrostatic potential from stabilizing into a perfectly periodic oscillation. By default, VASP's dipole correction algorithm expects the slab to be continuous in the center of the cell and the vacuum to be at the boundaries ($Z=0$ and $Z=c$). Alternatively, if it evaluates the center of charge automatically without an explicit DIPOL tag, the split nature of your slab misleads the algorithm. The potential step/jump intended to cancel the dipole field is likely being placed in the wrong position potentially overlapping with the atoms at the top of your box ($Z > 24\text{ \AA}$), which causes the sharp drop-off and leaves an uncompensated field in the vacuum. Modify your POSCAR so that the entire slab is continuous and centered in the middle of the simulation box (around $Z = c/2$). The vacuum should be split evenly at the bottom ($Z = 0$) and the top ($Z = c$). This allows VASP to easily identify the vacuum region and properly apply the dipole step function. |
3楼2026-05-20 17:26:26
40