光学带隙计算看着只是换个泛函的事,真做起来才发现空带数、k 点密度和带隙修正哪一项没调好,吸收边都能跑偏。我做光伏和光催化方向这些年,被”算出来的吸收谱和实验对不上”这个问题追着问过无数次,后来才摸清门道。

电子带隙是能带的直接差,光学带隙还牵涉激子效应和跃迁选择定则,二者经常不是一个值。很多用户拿着 PBE 算的电子带隙就去解释吸收边,结果差出零点几电子伏,论文里怎么圆都别扭。我们在做光催化体系时,几乎每次都要把光学带隙单独拎出来算,而不是顺手用电子带隙替代。一个 光学带隙计算 没做透的案子,吸收谱的峰位和强度都会失真。
光学性质的本质是计算频率依赖的介电函数,实部给出折射率、虚部给出吸收。VASP 里要先算一波函数的导带部分(空带),再求光学矩阵元,对离子做求和。这一步对空带数量极其敏感——空带不够,高能区直接塌方。我们习惯把 NBANDS 开到价带数两到三倍再测收敛,而不是用默认值碰运气。
第一是空带数,前面说了,必须收敛测试。第二是 k 点密度,光学量对 k 点远比总能达到的更密,稀疏 k 点会让谱线一顿一顿。第三是带隙修正,PBE 低估的带隙会把整个吸收边压低,得上 HSE06 或最后做 scissor/GW 修正。我们内部对每种材料都有”光学参数模板”,新体系照着测一遍收敛就能直接用,省掉大量试错。
落地步骤:几何优化→自洽电子密度→读取波函数算介电函数(开够空带、加密 k 点)→后处理得到吸收系数与反射率。当你认真做 光学带隙计算 的收敛验证,会发现谱形从”乱跳”变成”服帖”只差那几次测试。最后我们还会把算出来的吸收边和实验紫外可见数据叠在一起比对,偏差超过经验范围就回头查参数。
最高频的是空带不足导致高能吸收被砍掉,看着像材料不吸光。其次是 k 点太疏,谱线振荡大得像噪声。还有人忘了做带隙修正,吸收边整体偏低零点三五电子伏,正好卡在光催化判据的临界上,结论直接翻车。这些都不是软件问题,是参数没测收敛,靠规范的测试流程就能根治。
回过头看,光学带隙计算最怕”想当然”。空带要测、k 点要密、带隙要修,这三件事做不到,谱再好看也是假象。我们交付的光学性质报告,必附收敛曲线和实验比对,客户拿去写论文心里有底。把这三条钉死,光学带隙这关就稳了。
cp2k第一性原理计算
CP2K模拟计算
CP2K模拟计算
castep计算吸附能:Materials Studio里的表面建模与能量收敛复盘
CP2K计算能带:大体系电子结构模拟的混合基组方案
能带理论计算:固体能带结构的DFT模拟方法与工程应用
CP2K计算能带:混合基组DFT方法在周期性体系中的实战应用
CP2K分子动力学模拟详解:大体系加速策略与GPW方法实战