计算化学公社

标题: vasp采用PBE0算的能隙是负值，请问如何解决？ [打印本页]

作者
Author: 西装暴徒 时间: 2024-6-23 11:27
标题: vasp采用PBE0算的能隙是负值，请问如何解决？
本帖最后由西装暴徒于 2024-6-23 12:14 编辑

各位老师好，我在PBE下优化，用PBE0计算能带，但是能带是负值，请问是为什么？
KPOINTS是用vaspkit251生成的

体系是氧化银

作者
Author: 乐平 时间: 2024-6-23 21:56
本帖最后由乐平于 2024-6-24 11:00 编辑

(1) 结构优化完毕后，先用 PBE 泛函自洽计算得到 WAVECAR, CHGCAR

(2) 假设你的体系是氧化银块体，用 VASPKIT —> 303 得到 KPATH.in

(3) 用 VASPKIT ——> 251 得到杂化泛函的 KPOINTS， VASPKIT ——> 101 设置 ST 自洽计算 INCAR，先使用 PBE 泛函，按杂化泛函的 KPOINTS，读取第（1）步的 WAVECAR, CHGCAR （即 ISTART = 1, ICHARG = 1）得到 PBE 泛函的能带。并新建一个目录，将此能带结果移动到新建的目录里。

(4) 在第（3）步的 INCAR 里设置 PBE0 杂化泛函的参数，其余不变，读取第（3）步里得到的 WAVECAR, CHGCAR 信息（即 ISTART = 1, ICHARG = 1），重新计算 PBE0 杂化泛函的能带。

作者
Author: 西装暴徒 时间: 2024-6-24 10:49

乐平发表于 2024-6-23 21:56
(1) 结构优化完毕后，先用 PBE 泛函自洽计算得到 WAVECAR, CHGCAR

(2) 假设你的体系是氧化银块体，用 VA ...

谢谢老师，我可能就是vaspkit设置的302导致能带是负值，303算了之后是正值。老师，这两个会产生很大差异吗？

作者
Author: 乐平 时间: 2024-6-24 10:53

西装暴徒发表于 2024-6-24 10:49
谢谢老师，我可能就是vaspkit设置的302导致能带是负值，303算了之后是正值。老师，这两个会产生很大差异 ...

VASPKIT 302 是针对二维材料的，你没说你的体系氧化银是二维还是块体

作者
Author: 西装暴徒 时间: 2024-6-24 17:04

乐平发表于 2024-6-24 10:53
VASPKIT 302 是针对二维材料的，你没说你的体系氧化银是二维还是块体

老师，我该如何看是否是二维或者是块体？老师这个二维类似是晶体的切面吗？如果是的话我的氧化银体系就是块体。

作者
Author: 西装暴徒 时间: 2024-6-24 17:48

乐平发表于 2024-6-24 10:53
VASPKIT 302 是针对二维材料的，你没说你的体系氧化银是二维还是块体

老师，我看了您B站的CASTEP教学视频，故再请教您一个CASTEP的态密度投影问题，我如何知道投影的轨道是哪个原子的S或者P轨道？这些投影的轨道是不是这些原子的价层电子轨道？

作者
Author: 乐平 时间: 2024-6-24 17:58

西装暴徒发表于 2024-6-24 17:04
老师，我该如何看是否是二维或者是块体？老师这个二维类似是晶体的切面吗？如果是的话我的氧化银体系就是 ...

三维是 x, y, z 三个方向都按周期性重复并且铺满了重复单元的原子。

二维就是比三维少一个维度。虽然计算的时候仍然是 x, y, z 三个方向周期性重复，但是某一个维度上并不是铺满了原子，而是有较大的真空层间隔。比如你说的切了表面。

作者
Author: 乐平 时间: 2024-6-24 18:01

西装暴徒发表于 2024-6-24 17:48
老师，我看了您B站的CASTEP教学视频，故再请教您一个CASTEP的态密度投影问题，我如何知道投影的轨道是哪 ...

先选中某个原子，重新 analysis 画 PDOS 就行了

作者
Author: 西装暴徒 时间: 2024-6-24 21:51

乐平发表于 2024-6-24 18:01
先选中某个原子，重新 analysis 画 PDOS 就行了

好的，谢谢老师。

作者
Author: 西装暴徒 时间: 2024-6-24 21:52

乐平发表于 2024-6-24 17:58
三维是 x, y, z 三个方向都按周期性重复并且铺满了重复单元的原子。

二维就是比三维少一个维度。虽然 ...

谢谢老师，那就是块体。

作者
Author: Autumn899 时间: 2025-3-15 21:35
您好，请问用PBE0计算完能带应该怎么处理数据呢？我处理完的能带在费米能级以下出现很多刺，能带不是很平滑，是怎么回事啊？

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