DFT计算材料相容性：界面的电子结构不匹配，宏观相容性无从谈起

DFT计算材料相容性的切入点和传统热力学相容性评估有本质区别——热力学看的是混合焓的正负，DFT看的是界面处原子尺度的电子结构是否匹配。材料相容性在工程语境里是一个被滥用的词。在工程语境里，”这两种材料相容吗”的意思是”它们在一起会不会分层、会不会反应、高温下会不会互相扩散”。这些问题追到底，答案在界面上的头几个原子层里。

一个做高温涂层的项目把问题带过来的方式是直接的——陶瓷涂层（Al₂O₃）镀在镍基高温合金基体上，服役温度900°C，100小时后涂层局部剥落。SEM截面显示剥落发生在界面处，不是涂层内部。失效模式指向界面结合本身。

界面建模的几何切入点

Al₂O₃(0001)/Ni(111)是最常见的取向关系——Al₂O₃的氧密排面和Ni的密排面形成半共格界面。晶格失配度约12%，需要一个4×4 Al₂O₃搭配5×5 Ni的超级胞来将失配应力控制在~3%以内。

界面建模这一步就在做判断。一种思路是删掉失配——把Ni层在面内拉伸到与Al₂O₃完美共格（coherent interface approximation）。这种做法物理上不对——12%的应变对Ni来说已经过了弹性极限，引入的应变能会系统性地高估界面结合能——但计算上简单。另一种思路是做半共格（semi-coherent），在Ni层中引入失配位错，这需要在面内方向做更大的超级胞（>400原子），计算量膨胀。

这个项目里选了折衷：保留4×4 Al₂O₃不变，Ni层取5×5拉伸3%——Ni的弹性模量约200GPa，3%应变对应应力约6GPa，虽然已经在塑性区边缘，但界面处的Ni原子可能通过局部弛豫释放部分应力。弛豫后界面处Ni原子偏离理想Ni(111)约0.15Å，残余应力降到1.2GPa。做这个折衷是因为：不做界面弛豫的界面结合能是纯粹的力学储存能加上化学键合能——两个量搅在一起没法分。做了弛豫之后分开看。

结合能只是第一步

Al₂O₃(0001)/Ni(111)界面的理想结合功（W_ad）在这个项目里算出来是3.82 J/m²。这个数字放在界面结合能的数据库里属于”中强”结合——高于石墨烯/Ni（约1.5 J/m²），低于TiN/MgO（约5.5 J/m²）。

但结合能是一个标量——只告诉你”结合有多强”，不告诉你”结合是什么类型的”。同是3.82 J/m²，如果结合来自范德华力、或者来自Al₂O₃表面氧和Ni的共价键、或者来自界面处的金属键——这三种键合类型的热稳定性和高温响应完全不同。

电荷密度差和静电势对齐

我们接着算了界面电荷密度差——Δρ = ρ_slab − ρ_Al2O3 − ρ_Ni。在Al₂O₃侧的界面氧和Ni侧的界面镍之间，出现了一个清晰的电子积累区（Δρ > 0），峰值位置偏Ni侧约0.6Å。这表明电荷从Al₂O₃的O-2p向Ni-3d转移。

同时算的静电势对齐也印证了这一点。Al₂O₃的功函数约5.8eV（绝缘体），Ni约5.2eV（金属）。接触后，Ni的费米能级等于Al₂O₃的价带顶加界面偶极修正——电子从Ni流向Al₂O₃侧，在界面处形成偶极层。这个偶极层的存在意味着界面键合有离子性成分——不是纯金属键。

结合起来：Al₂O₃/Ni界面结合以离子-共价混合键为主，带少量金属性。这个类型的键合在高温下（>800°C）面临两个退化路径：氧空位生成（Al₂O₃侧氧离子的热激活迁移）和Ni原子扩散进入Al₂O₃晶格。后者在我们算的界面态密度里有证据——Ni-3d在Al₂O₃侧的带隙中引入了间隙态，深度约1.2eV。

这个间隙态是Ni扩散的前驱信号——电子层面先出现”混合”，原子层面的扩散才跟得上。900°C下100小时的剥落，原因不在宏观热胀系数不匹配（我们也算了，膨胀系数差约3×10⁻⁶/K，不足以造成剥落），而在界面化学——Ni扩散进入Al₂O₃后形成的NiAl₂O₄尖晶石相体积膨胀了约3%，是热应力失效的源头，不是涂层设计问题。

DFT计算材料相容性的核心判据不是结合能，是界面化学——电子从哪里转移到哪里、形成了什么类型的键、在高温下这些键怎么退化。