MS计算介电函数：一个TiO2光催化项目的能带与光学响应复盘

介电函数的计算，在DFT领域算不上最难的，但绝对是最容易被轻视的环节之一。很多项目把介电函数当成能带计算的附属产物随手提取，结果拿到的吸收边位置跟实验差出0.5 eV以上，最后只能怪”杂化泛函没开”或者”GW修正没做”——其实问题往往出在更基础的地方。

Globe background with binary code flowing from a computer

项目背景：为什么介电函数比能带更难算准

去年团队接了一个TiO2光催化剂的光学响应评估项目。客户需要知道锐钛矿TiO2在紫外-可见光范围内的吸收系数和折射率，用来评估光催化效率。听起来是标准操作——用Materials Studio的CASTEP模块算能带，顺带提取介电函数，画出吸收谱就行。

但拿到实验数据之后发现，对方的UV-Vis吸收谱上吸收边在3.42 eV附近，而团队用PBE泛函跑出来的带隙只有3.21 eV，对应的吸收边自然也偏移了约0.2 eV。这个偏差量不算大，但客户要求吸收系数的绝对值和实验的偏差控制在10%以内，这意味着不能只靠”修带隙”来凑数，必须从介电函数的计算细节上做功课。

截断能和k点：光学计算的特殊要求

DFT计算中截断能和k点密度的选择是老生常谈，但介电函数对这两个参数的敏感度比总能高得多。TiO2的介电函数虚部在带边附近变化极其陡峭——从零迅速上升到峰值，然后在5-8 eV区间维持较高值。如果k点不够密，带边的精细结构会被抹平，导致峰值位置漂移和高度变化。

团队用PBE泛函做了两组收敛测试：一组是截断能从400 eV到700 eV以50 eV步长递增，k点固定为6×6×3；另一组是截断能固定在550 eV，k点从3×3×2到9×9×5递增。结果显示，截断能550 eV以上时总能变化小于0.01 eV/atom，但介电函数虚部的峰值在550 eV和600 eV之间仍有约3%的差异。k点方面，6×6×3是一个临界点——低于这个密度时，介电函数高频区（>6 eV）的振荡特征明显变弱，跟实验对比的吻合度直线下降。

最终选定的参数是截断能600 eV、k点网格6×6×3、smearing宽度0.05 eV。这个组合不是计算量最小的，但在介电函数的带边和高频区都能给出稳定的收敛结果。

从带隙修正到吸收谱的完整链路

PBE泛函低估带隙是已知问题，TiO2的PBE带隙约3.2 eV，而实验值3.42 eV。团队没有直接上杂化泛函（HSE06对TiO2的计算量太大，客户预算和时间都不允许），而是采用了scissor算子修正——将导带整体上移0.22 eV，使计算带隙与实验一致。

这个操作在CASTEP中可以直接在光学性质设置里指定，简单粗暴但在带边位置对齐方面的效果是可靠的。需要注意的是，scissor算子只修正了带隙，对有效质量没有影响，所以介电函数的峰位和峰形是由原始DFT计算决定的，scissor只是把整个响应沿能量轴平移了0.22 eV。

经过scissor修正后，计算得到的吸收边出现在3.41 eV，与实验的3.42 eV几乎完全吻合。吸收系数在4.0-5.5 eV区间的绝对值偏差约为7%，满足客户要求的10%以内。折射率的计算结果在可见光范围（1.6-3.1 eV）内与文献椭偏测量值的偏差在5%左右，也是可接受的水平。

值得注意的两个技术细节

第一，CASTEP计算介电函数时默认只考虑偶极近似，忽略局域场效应（local field effects）。对于像TiO2这样的宽禁带绝缘体，局域场效应对介电函数的影响通常较小，但在带边附近的精细结构上可能贡献1-2%的修正。如果项目对精度要求更高，需要在光学性质设置中勾选”Include local field effects”，计算量会增加到原来的3-4倍。

第二，介电函数的实部和虚部必须满足Kramers-Kronig关系——它们不是独立的，CASTEP内部通过线性响应理论计算的是虚部，实部是通过K-K变换得到的。如果提取的实部和虚部看起来不协调，先检查能量范围是否足够宽（建议从0覆盖到至少15 eV），截断过低会导致K-K变换在低能区产生不准确的尾巴。

这个项目最终交付的计算报告包含了完整的能带结构、态密度、介电函数（实部+虚部）、吸收系数和折射率五个部分。客户拿着这些数据和他们的UV-Vis、椭偏结果做了对比，认为计算精度满足后续论文投稿的要求。这种光学性质的计算不算多复杂，但每个参数选择背后都有其物理依据——介电函数不是一个”算完能带顺便提一下”的东西，它值得被认真对待。

如需了解更多的科研计算案例与项目经验，可以参考 https://www.keyanxueshu.com/ 上的技术文章。