3. DFT计算光学性质的理论框架

在DFT框架下计算光学性质，核心困难在于：基态DFT给出的Kohn-Sham能级与波函数本身并不是严格意义上的准粒子激发态。要把基态计算结果映射到“光吸收/发射”这一激发态物理过程，必须引入不同层级的近似方法。本节按“物理假设—适用条件—工程选择依据”三条线索，梳理独立粒子近似（IPA）、随机相位近似（RPA）、时间依赖DFT（TDDFT）与Bethe-Salpeter方程（BSE）四种主流框架，并给出可操作的选型判据。

一、四种方法的物理图像与层级关系

四种方法并非彼此独立，而是沿“电子-空穴相互作用是否被考虑、交换关联效应如何处理”两个维度逐级递进。独立粒子近似把光学跃迁视为单电子从占据态到空态的独立激发，完全忽略电子-空穴之间的库仑吸引与局域场效应；RPA在IPA基础上加入电子-空穴对的集体响应（即局域场效应），但仍不包含电子-空穴的束缚态；TDDFT通过含时交换关联核描述电子体系的集体激发，可处理激子效应但依赖交换关联泛函；BSE则显式求解电子-空穴对的双粒子格林函数方程，同时包含屏蔽库仑吸引项与交换排斥项，是当前描述激子最严格的方法。

二、独立粒子近似（IPA）：最快的基线

IPA直接使用DFT基态的Kohn-Sham能级与波函数，通过费米黄金规则计算跃迁矩阵元，虚部介电函数写作对占据态v与空态c求和的表达式。其最大优势是计算量小，只需一次基态自洽计算加一次非自洽能带计算即可得到完整介电谱。但IPA存在两个系统性偏差：一是DFT带隙通常比实验值小30%～50%，导致吸收边红移；二是完全忽略电子-空穴相互作用，对激子峰无法描述。

工程上使用IPA时，必须做带隙修正。常用做法是“剪刀修正”：将空带能量整体上移Δ，使带隙对齐实验值或GW值。例如对Si，PBE带隙约0.6 eV，实验值1.1 eV，剪刀量约0.5 eV。修正后吸收边位置可显著改善，但谱形权重不变。IPA适合快速筛选、金属光学常数计算以及激子效应可忽略的体系。

三、随机相位近似（RPA）：加入局域场

RPA在IPA基础上引入电子密度响应，考虑局域场效应。其核心是求解Dyson方程 χ = χ₀ + χ₀ v χ₀ + …，其中v为库仑核。RPA能正确描述等离激元集体激发，对金属和重掺杂半导体的低频介电行为明显优于IPA。在半导体中，RPA对E₂峰（如Si的4.3 eV峰）的位置和强度有改善，但仍不包含激子束缚效应。

RPA有两种常见实现：忽略局域场（RPA-noLF）与包含局域场（RPA-LF）。前者等价于IPA的独立粒子响应，后者才真正体现RPA价值。工程中若使用VASP的LOPTICS或Yambo，需明确开启局域场效应开关。RPA计算成本约为IPA的2～5倍，仍属可接受范围。

四、时间依赖DFT（TDDFT）：激子效应的实用选择

TDDFT通过含时Kohn-Sham方程描述电子体系对光场的响应，其核心是交换关联核f_xc。在绝热近似下，f_xc只依赖瞬时密度，忽略记忆效应。TDDFT的优点是能同时处理基态与激发态，且计算量介于RPA与BSE之间。对分子体系，TDDFT（尤其杂化泛函）能给出与实验吻合良好的激子峰位与振子强度。

但TDDFT存在两个关键限制：一是对固体，标准局域或半局域泛函的f_xc衰减过快，导致激子束缚能严重低估；二是对电荷转移激发，绝热近似下误差较大。工程上对固体激子计算，更推荐使用BSE而非TDDFT。TDDFT更适合有限体系（分子、团簇、量子点）的光学性质计算。

五、Bethe-Salpeter方程（BSE）：激子的严格描述

BSE是双粒子格林函数方程，描述电子-空穴对的传播。其本征值给出激子能级，本征矢给出激子波函数。BSE的核心是电子-空穴相互作用核：屏蔽库仑项W（吸引）与交换项v（排斥）。求解BSE需要先做GW计算得到准粒子能级，再构建电子-空穴对基组，最后对角化BSE矩阵。

BSE的计算流程通常为：DFT基态 → GW准粒子修正 → BSE激子求解。以VASP为例，典型参数为：GW截断能约200～400 eV，BSE的k点网格需足够密（如12×12×12），电子-空穴对数目可达数千。BSE能准确给出激子束缚能、激子峰位与振子强度，对有机半导体、二维材料、宽禁带氧化物尤为必要。

六、工程选择决策表

总结而言，IPA是快速基线，RPA加入局域场，TDDFT适合有限体系，BSE是激子计算的黄金标准。工程选择的核心判据是激子束缚能与室温kT的比较，以及材料维度与目标物理量。实际计算中，建议按“IPA → RPA → BSE”逐级验证，确保每一步的修正量合理，最终结果与实验或更高精度计算交叉验证。