手机版

VASP 计算 HOMO/LUMO

发布时间：2026-08-31 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

最高占据分子轨道（HOMO）与最低未占据分子轨道（LUMO）是描述分子电子结构的核心能级，其能量差即 HOMO-LUMO 带隙（Eg），直接决定分子的光学吸收、导电性、氧化还原电位与化学反应活性。对有机半导体材料，HOMO/LUMO 能级与电极功函数/电子亲和势的相对位置决定了电荷注入的势垒，是器件设计的关键参数。因此，准确计算 HOMO/LUMO 能级对材料筛选与器件优化具有直接价值。

周期性 DFT 计算（VASP）将分子置于周期性盒子中，通过电子结构计算获得能级与态密度，可同时考虑分子间相互作用与晶体环境的影响，是研究有机固体与分子晶体的标准方法。本项目以有机半导体分子为对象，利用 VASP 计算 HOMO/LUMO 能级、带隙与态密度，分析其电子结构特征。

二、计算方法与理论背景

2.1 HOMO/LUMO 与带隙

对分子体系，HOMO 对应最高占据能级（价带顶），LUMO 对应最低未占据能级（导带底）。在周期性 DFT 中，HOMO-LUMO 带隙可由总态密度中占据态与非占据态的能级差直接读取。需要指出的是，DFT（尤其 GGA 泛函）对带隙存在系统性低估，原因是 Kohn-Sham 本征值并非严格意义的准粒子能量；对有机分子体系，这一低估通常为 0.5–1.5 eV。

2.2 周期性盒子与真空层

对孤立分子，需将其置于大尺寸周期性盒子中，并在各方向添加足够真空层（通常 15–20 Å）以抑制分子与其镜像之间的虚假相互作用。此时分子能级表现为离散的尖锐峰，可直接读取 HOMO/LUMO。若研究分子晶体，则采用实验晶体结构并考虑分子间耦合。

2.3 能级对齐

HOMO/LUMO 能级需相对真空能级对齐以获得电离势（IP = −E_HOMO）与电子亲和势（EA = −E_LUMO）。VASP 中可通过静电势计算获得真空能级参考，将 Kohn-Sham 能级平移对齐。本案例中真空能级定为 0 eV，HOMO = −5.2 eV、LUMO = −3.4 eV，对应 IP ≈ 5.2 eV、EA ≈ 3.4 eV，带隙 1.8 eV。

三、计算设置与模型构建

软件：VASP；泛函：GGA-PBE，PAW 赝势。
体系：单个有机半导体分子置于 20 Å×20 Å×20 Å 盒子（真空层约 15 Å）。
截断能：400 eV；k 点：Gamma 点（孤立分子）。
收敛：能量 1×10⁻⁶ eV；结构优化至力 < 0.01 eV/Å。
输出：总态密度（DOS）、分子轨道能级、HOMO/LUMO 波函数。

四、结果分析与案例图解读

图 2 的左面板给出有机半导体分子的 HOMO/LUMO 能级图。真空能级定为 0 eV（灰色横线），LUMO 位于 −3.4 eV（蓝色能级条），HOMO 位于 −5.2 eV（红色能级条），两能级间以双向箭头标注带隙 Eg = 1.8 eV。能级图清晰显示：HOMO 深入占据区 1.8 eV 以下，LUMO 位于真空能级下方 3.4 eV 处。由此可提取 IP ≈ 5.2 eV、EA ≈ 3.4 eV，这些数值直接用于器件设计——例如与电极功函数匹配以降低空穴/电子注入势垒。若将 LUMO 与常见电极功函数（如 Al 4.3 eV、ITO 4.7 eV）比较，可评估电子注入的可行性。

图 2 的右面板给出分子在带隙附近的总态密度（DOS）曲线。紫色曲线呈现两组尖锐的峰：低能侧（−6 至 −4 eV）的占据态峰群以红色填充标示，对应 HOMO 及其下层占据轨道；高能侧（−4 至 −2 eV）的非占据态峰群以蓝色填充标示，对应 LUMO 及其上层空轨道。两峰群之间存在明显的能隙（DOS 为零的区域），宽度 1.8 eV，与左面板的带隙一致。DOS 的尖锐峰形说明分子间相互作用弱（近孤立分子），电子态定域化明显，这与有机分子晶体的弱范德华耦合特征相符。

两面板结合，左面板给出 HOMO/LUMO 的绝对能级（用于能级匹配设计），右面板给出态密度的整体分布（用于理解电子态结构与带隙来源），共同构成有机半导体分子电子结构的完整描述，为有机光电器件（OLED、OFET、OPV）的材料选型提供定量能级依据。