vasp计算homo

项目过程：并五苯前线轨道的 VASP 计算

客户研究有机半导体载流子传输，需要并五苯（pentacene）的 HOMO-LUMO 能隙和前线轨道空间分布。他之前用 Gaussian 算气相分子，现在想知道晶体周期性边界下 VASP 的结果有什么差别。

我建了并五苯单分子层/薄膜模型，用 VASP 做周期性 DFT 计算。由于 PBE 会严重低估有机半导体带隙，这里采用 HSE06 杂化泛函。截断能 400 eV，K 点取 Γ 点（分子层间相互作用弱，Γ 足够），真空层 20 Å 避免层间相互作用。

计算分两步：先结构优化，再算 DOS 和能带，最后用 WAVECAR 画 HOMO 和 LUMO 的等值面。

案例分析：能隙 2.20 eV，HOMO 和 LUMO 都是离域 π 轨道

第一张案例图是并五苯的态密度（DOS）。价带顶对应 HOMO，导带底对应 LUMO，两者之间的能量差标注为 E_gap=2.20 eV。这个数值比气相光学吸收实验值略高，但在周期性薄膜模型里加上 HSE06 修正后，属于合理范围。

第二张案例图给出 HOMO 和 LUMO 的空间分布。HOMO 是占据的 π 轨道，五个苯环上都有电子云分布，相位交替呈现“红蓝红蓝”的节面结构；LUMO 是虚的 π* 反键轨道，节点位置和 HOMO 互补。这种离域性说明并五苯的载流子迁移率在面内方向会比较高，而层间方向由于 π-π 堆叠距离影响会低很多。

总结教训：周期性模型要区分“分子”和“晶体”

第一，VASP 算单分子必须加大真空层。20 Å 是最低要求，15 Å 时相邻单胞的分子波函数会有明显交叠，带隙会被压低。

第二，有机半导体一定要用杂化泛函或 GW。PBE 给出的带隙可能只有 1.0 eV 出头，完全偏离实验，不能直接用。

第三，DOS 和能带结果要对照。如果 DOS 在带边很尖锐，说明能带色散小，载流子有效质量大；并五苯 HOMO 带的 DOS 峰较窄，符合其低迁移率半导体的定位。

项目经验：前线轨道图比数值更有说服力

对于有机半导体客户，他们往往不是做理论的，一张 HOMO/LUMO 等值面图比 2.20 eV 这个数字更容易理解。以后做类似项目，我会把 DOS、能带、轨道图、部分态密度分解放在一页里，让客户一眼看到“电子在哪里、能隙在哪里、传输通道在哪里”。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：vasp计算homo。