1. DFT介电与光学性质计算工程导论

在材料研发中，很多关键性能并不是由“总能量最低”直接决定的，而是由材料与电磁场如何相互作用决定。介电函数描述材料在外加电场下如何极化、储能与耗散，光学性质则进一步表现为折射、吸收、反射、透射与色散。对于光伏吸收层、显示发光材料、透明导电氧化物、介电储能薄膜、红外窗口、光催化与光电探测材料而言，介电函数和光学响应往往比形成能、带隙数值更接近器件指标。本课程的目标，是把密度泛函理论（DFT）从“能带计算工具”推进到“介电与光学性质计算工程”，让学习者能够独立完成从结构准备、参数设置、介电函数计算、光学量提取到结果判读与报告输出的完整链路。

介电函数通常写作复函数 ε(ω)=ε₁(ω)+iε₂(ω)。其中 ε₁ 是实部，反映极化与储能能力，决定折射率与介电常数；ε₂ 是虚部，反映吸收与耗散，直接关联光吸收。二者通过 Kramers-Kronig 关系相互约束，因此不能只算一个分量。由介电函数可以导出折射率 n、消光系数 k、吸收系数 α、反射率 R、能量损失函数以及光电导率等。工程上最常看的指标包括：静态介电常数 ε₁(0)、电子部分介电常数、带隙与吸收边位置、可见光区平均吸收系数、折射率色散曲线、反射率平台以及等离子体频率附近能量损失峰。需要特别说明的是，DFT 计算得到的是电子贡献的介电函数，若材料存在显著晶格极化，还需引入声子贡献或采用 DFPT 方法计算离子部分介电常数。

DFT 在其中的角色可以概括为“从电子结构到光学响应的桥梁”。第一性原理计算先给出基态电子密度、能带与波函数，再通过独立粒子近似、随机相位近似（RPA）或 Bethe-Salpeter 方程（BSE）计算介电响应。对于常规半导体和绝缘体，独立粒子近似加剪刀修正已能给出可用的吸收边与折射率趋势；对于激子效应显著的二维材料、有机晶体或宽禁带氧化物，则需要考虑局域场效应与电子-空穴相互作用。工程实践中，计算精度不仅取决于交换关联泛函，还取决于 k 点密度、空带数量、展宽参数与剪刀修正量。一个常见误区是只关注带隙是否与实验一致，却忽略介电函数对 k 点收敛和空带数极其敏感。

一个完整的工程工作流通常包含五个阶段。第一，结构准备：确认晶体结构、空间群、超胞尺寸与真空层厚度，二维材料需保证层间真空大于 15 Å，避免周期性镜像相互作用。第二，基态自洽计算：选择交换关联泛函，设定截断能与 k 点网格，收敛阈值建议总能量 1E-6 eV、力 0.01 eV/Å。第三，能带与态密度分析：判断带隙类型与带边位置，确定是否需要剪刀修正。第四，介电函数计算：设置足够空带数（通常为占据态数的 2～3 倍）、展宽参数 η（常用 0.05～0.1 eV）与频率网格。第五，光学量提取与判读：由 ε₁、ε₂ 导出 n、k、α、R，并与目标器件指标对照。

以典型教学示例为例，假设计算一种立方钙钛矿型氧化物，采用 PBE 泛函，截断能 520 eV，k 点网格 8×8×8，空带数设为 48，展宽 0.08 eV。计算得到带隙约 1.8 eV，而实验值约 3.2 eV，此时需施加约 1.4 eV 的剪刀修正。修正后 ε₂ 的吸收边从 1.8 eV 蓝移至约 3.2 eV，可见光区吸收系数显著下降，折射率在低能区约为 2.1。若 k 点降至 4×4×4，ε₂ 峰形出现明显毛刺，吸收边位置偏移可达 0.2 eV 以上，说明 k 点收敛不可省略。这些数值属于教学示例假设，实际计算需根据具体材料与收敛测试确定。

本课程后续将围绕这条工作流展开：从结构建模与泛函选择，到介电函数计算设置、光学量提取、激子效应处理、缺陷与掺杂影响、界面与薄膜体系，再到结果可视化与报告输出。学习者最终应能独立完成一套可复现、可判读、可交付的 DFT 介电与光学性质计算方案，并把计算结果转化为材料研发中的有效决策依据。