dft计算材料相容性

我做过一个挺典型的异种金属界面筛选案例。客户是做铝钛复合结构扩散焊的工程师，焊前最怕的就是两片金属压到一起，要么焊不上、要么焊完界面脆裂。他手里列了五组候选：Al/Ti、Cu/Ni、Mg/Al、Fe/Cr、Mo/W，想让我在动手做实验之前，先用 dft计算材料相容性 把”热力学上到底稳不稳”这件事回答掉，省掉一轮盲目试错。

我的思路不是直接去算一个界面结合能就交差，而是搭一个三层筛选框架：晶格失配 → 界面形成能 → 界面结合能，逐层过滤。第一关是晶格失配。两个晶体要拼界面，第一步得分别优化各自的体相，拿到平衡晶格常数 a，然后算失配度 ε = |a₁ − a₂| / a₂。经验上如果 ε 超过 5%，硬拼会引入几百分之十的应变，算出来的界面能基本是假的。所以我先画了一张”失配度 vs 界面形成能”的扫描曲线，找到 |ε| < 3% 的兼容窗口，只有落进窗口的体系才进下一轮。Al/Ti 的失配约 2.8%、Mg/Al 约 2.1%，都在窗口内;Mo/W 虽然失配小但本质同源，意义不大。

第二关是界面模型构建。这一步最容易翻车。我没有用两个半无限大 slab 硬怼，而是各自切出低指数表面(Al(111)/Ti(0001) 这类密排面优先)，再拼成三明治结构：每侧至少 6 层原子，真空层 ≥ 15 Å 防止上下表面相互耦合，底部固定 2–3 层模拟半无限衬底，顶部弛豫。计算参数上，平面波截断能取 500 eV，k 点网格按体系大小取到 4×4×1 到 6×6×1.收敛判据是原子受力 < 0.02 eV/Å、能量变化 < 1×10⁻⁵ eV。磁性也不能忽略——Fe/Cr 是反铁磁耦合，初始磁矩设错会让总能差出几十 meV/atom，结果直接失真。

第三关才是结合能。公式用 ΔG = (E_interface − E_slab1 − E_slab2) / N，N 是界面原子数。我把五组结果排成柱状图：Al/Ti 和 Mg/Al 的 ΔG 都在 −0.4 eV/atom 附近，Cu/Ni 约 −0.3.Fe/Cr 和 Mo/W 接近 −0.15.我把 −0.2 eV/atom 当成”稳定界面”的门槛，最终推荐 Al/Ti、Mg/Al 优先推进实验验证。为了佐证，我还补了界面电荷密度差分和投影态密度(PDOS)：Al/Ti 界面能看到明显的电荷重新分布，说明是化学键合而非简单物理接触，这和它更负的结合能是自洽的。

我这个项目里踩过几个具体的坑，值得单独说。其一是赝势和强关联。金属界面一旦含 Ti、V 这类 3d 过渡金属，标准 PBE 会低估 d 电子的强关联，结合能可能偏 20% 以上。我一开始用纯 PBE 跑 Al/Ti，界面能偏低、和实验趋势对不上，加了 DFT+U(U≈4 eV 作用在 Ti 的 3d 上)之后才和文献吻合。其二是真空层，有次为了省算时把真空层压到 10 Å，上下表面静电镜像导致总能抖动，后来加到 18 Å 才平稳。其三是固定层太少，底部只固定 1 层时出现了虚假的表面重构，加到 3 层才稳住。

从我的工程经验看，dft计算材料相容性 这件事最大的误区，是以为”结合能为负就万事大吉”。DFT 只回答热力学稳不稳，而实际扩散焊还要看扩散系数(分子动力学去补)、脆性金属间化合物相的生成(CALPHAD 去补)。所以我交付报告时，明确给客户标注了边界：Al/Ti、Mg/Al 热力学兼容性好，但 Al/Ti 在 400℃ 以上会生成 TiAl₃ 脆性层，工艺温度窗口要卡住。更多类似项目的完整复盘，可以看 [dft计算材料相容性](https://www.keyanxueshu.com/category/dft/) 的案例库;如果要从头搭建这套筛选流程，[dft计算材料相容性](https://www.keyanxueshu.com/) 上也有对应的服务入口可以参考。