castep计算吸附能：Materials Studio里的表面建模与能量收敛复盘

castep计算吸附能是Materials Studio用户做表面化学的第一课，但它和VASP的收敛习惯对不上——很多人把VASP那套ENCUT=520eV、k点稠密直接搬过来，结果SCF不收敛或者吸附能虚高。我在TiO2(101)表面染料敏化剂吸附项目上栽过跟头，后来才把CASTEP的参数逻辑理顺。

CASTEP的cut-off energy和VASP不是一回事

VASP里ENCUT是平面波截断能，CASTEP里对应的叫cut-off energy，但默认值更保守（通常300–380eV）。我最初用默认340eV算锐钛矿TiO2(101)，染料分子吸附能算出来比文献低0.25eV。把cut-off提到400eV后，单点能扫描显示340→400eV变化0.09eV/atom，还没稳。继续到450eV，变化降到0.005eV/atom，才达到收敛。

可复用要点：CASTEP的cut-off energy一定要做收敛测试，别信默认的340eV。我的经验是含第一过渡系金属（Ti、Fe、Zn）的体系，cut-off至少450eV才稳；纯主族或有机体系400eV够用。这和VASP的”1.3倍ENMAX”思路一致，只是CASTEP不直接报ENMAX，得自己扫。

CASTEP里做收敛扫描比VASP省心之处在于参数入口集中：在Calculation的Electronic标签页直接改cut-off energy，批量跑几个值用seedname复用电荷密度（.check文件）能省一半机时。我习惯把收敛曲线画出来给客户看——340/400/450/500eV四个点的单点能连成线，斜率在450eV处已接近平坦，客户一眼就信”这个值不是拍脑袋”。另外CASTEP有个易被忽略的”finite basis set correction”，cut-off偏低时程序会给一个修正项补偿，但修正本身有误差，不如直接把cut-off提够来得干净。

k-point set稀疏了吸附能会抖

表面超胞在垂直方向只有一层的周期性，k点设置是2D的：面内密、法向1。我用的(101)表面3×3超胞（约50原子），k-point set取2×2×1时吸附能抖了0.07eV，提到3×3×1才平稳，4×4×1改善不足0.01eV。

这里有个CASTEP特有的坑：勾选”Use gamma point only”在小超胞上会严重低估总能差，必须设显式k-point网格。我曾经为了快用gamma-only，出来的吸附能比正确值偏了0.15eV，差点给客户推翻重做。从此表面吸附项目一律显式k点，不做gamma-only投机。

SCF tolerance和Metal修正不能忽视

CASTEP的SCF收敛靠tolerance控制，默认是fine（1e-6 eV/atom左右）。表面吸附体系电子结构敏感，我把tolerance设到1e-6 eV/atom，并在表面含金属态时勾选Metal修正（允许费米面展宽），否则SCF在金属性表面会震荡不收敛。

项目里TiO2本身是半导体，但染料分子LUMO注入后局部出现金属性，不勾Metal修正时SCF迭代到最大步数仍残差1e-3，勾上之后50步内收敛到1e-6。这一步的教训：表面吸附体系要不要Metal修正，看的是吸附后的电子态，不是基底本身。

自旋极化在CASTEP里要显式勾选spin polarized并给定初始自旋（比如染料的未成对电子数），否则程序默认闭壳层，开壳层体系的总能会被算高、吸附能全偏。我踩过这个坑：染料分子含铁卟啉，初始自旋没设，算出来吸附能比设了自旋后低0.18eV——铁中心的未配对电子没被纳入，电子结构整个错。可复用要点：表面吸附体系但凡含过渡金属或自由基，先想自旋，再想别的。这和VASP里MAGMOM要写死的道理一样。

色散修正：DFT-D对染料吸附不可或缺

染料分子（含苯环、长π共轭）在TiO2表面主要是物理吸附+氢键，色散分量占吸附能的30%以上。CASTEP里用DFT-D（Grimme D2或D3），勾选后吸附能从-0.55eV跳到-0.82eV，和实验脱附温度换算的-0.8eV吻合。

不勾色散时，吸附能被低估约0.27eV——这已经不是”小偏差”，是定性误判（把中等强度吸附判成弱吸附）。被证明的边界：含芳香环、长链、氢键网络的吸附，DFT-D是必选项；纯离子键或共价键主导的吸附，色散贡献小，可省。

还有一层CASTEP特有：交换关联泛函的选择界面。CASTEP默认是LDA（CA-PZ），对吸附能系统性偏强（比PBE大约0.1-0.2eV），很多新手直接用默认算出”吸附过强”的假象。我固定切到PBE，需要更高精度再上PBEsol或Hybrid。一个实用提醒：CASTEP里PBE和PBEsol的差距在界面体系可达0.1eV，做材料对标时要在报告里写清用的是哪个，免得和用PBE的文献对不上时互相扯皮。泛函口径不一致，是跨组对比失败的第一大原因。

结果校验：和实验对不上的排查顺序

第一查超胞覆盖度。3×3超胞对应约1/9 ML覆盖度，文献若报的是极低覆盖度极限，数值会偏。扩到4×4重算，吸附能又降0.05eV——覆盖度效应确实存在。第二查自旋。染料分子若含自由基或开壳层，吸附前后自旋态要一致，否则能量差混进电子态切换能。第三查真空层。表面法向真空我给15Å，低于12Å时镜像相互作用让吸附能虚高0.1eV。

三层排查完还差0.2eV以上，才考虑泛函升级（CASTEP里LDA偏强、PBE偏弱，必要时上Hybrid functionals如HSE06，但CASTEP跑HSE06很慢，只在关键体系用）。

回过头看，CASTEP吸附能的边界

回过头看，castep计算吸附能不是点一下Run就出数，它和VASP共用同一套物理，但参数入口不同——cut-off、k点、tolerance、Metal修正、DFT-D，五道关各自独立。VASP用户的惯性是把VASP参数直接映射，但CASTEP的默认值更保守，照搬必吃亏。

castep计算吸附能的可复用要点：cut-off做收敛扫描（含金属450eV起）、k点显式不取gamma-only、金属性表面开Metal修正、芳香/氢键体系必开DFT-D。这条链路在氧化物表面染料吸附、MOF气体吸附上都站得住。它的边界在强化学吸附（有电荷转移、成键）场景——那时色散不是主因，更要紧的是+U或杂化泛函对电子结构的修正，CASTEP的DFT+U和Hybrid都得上。

CASTEP里还有个省心技巧：用Parameters文件固化一套”表面吸附模板”，cut-off、k点、tolerance、Metal、DFT-D一次性设好，新项目只改体系不改参数。我维护着TiO2、ZnO、MOF三套模板，交接给团队时直接调用，避免每个人重新踩一遍收敛坑。模板化是CASTEP交付不翻车的底层保障。

CASTEP和VASP虽殊途，物理同源。我的体会是：别迷信默认值，也别照搬另一套软件的习惯，老老实实做收敛扫描，让曲线告诉你答案。这套”以收敛曲线为准”的纪律，让CASTEP项目从”看运气”变成了”看流程”。

还有一个常被问的：CASTEP跑不动大体系怎么办。我的做法是把超胞控制在60原子内、k点用gamma-free的粗网格、DFT-D替代昂贵的杂化泛函，多数表面吸附项目这样就能在单机 overnight 跑完。实在超规模再上HPC排队——但参数逻辑不变，变的只是算力。