DFT计算材料：从结构预测到性质筛选的全链条工作流

DFT计算材料的价值不在单次计算，而在筛选效率。以二维拓扑材料筛选为例：从一个包含上百种候选组分的材料库中，通过DFT全链条筛选——结构预测、稳定性验证、能带拓扑分析——最终锁定少数几种具有非平庸Z2拓扑不变量的候选材料。如果用实验手段做同等规模的筛选，至少需要两年;而计算筛选可以在数周内完成。这就是计算材料学的杠杆效应。

筛选管线的第一层：结构生成与预过滤

127个候选组分来源于两个渠道：一半来自Materials Project数据库的已知结构变体(通过pymatgen API批量下载)，另一半通过离子替换生成——例如将已知的Bi2Se3中的Bi替换为Sb，Se替换为Te，生成Sb2Te3、Bi2Te3等衍生物。

预过滤阶段的标准很简单但有效：

价电子数：排除价电子数<8的体系(不足以形成有意义的能带结构)

原子半径兼容性：替换原子的半径偏差<15%(否则结构不稳定概率高)

磁矩初判：非磁性体系优先(拓扑性质的分析更直接)

经过预过滤后通常淘汰约30-40%的候选组分。这一步看似粗暴，但实际上节省了大量计算资源——后续的DFT优化每步至少消耗几十核小时。

第二层：结构优化与稳定性验证

剩余候选结构进入VASP结构优化流程。参数配置：

ENCUT=520 eV(高于所有元素的ENMAX最大值约30%)

EDIFF=1E-6 eV(电子步收敛标准，高精度要求)

EDIFFG=-0.01 eV/Å(离子步收敛标准)

ISIF=3(同时优化晶胞和原子位置)

K-mesh：以k点间距<0.03 Å⁻¹为标准自动生成

结构优化完成后，进入热力学稳定性验证。方法是计算形成能并与凸包图对比：

ΔH_f = E(AB) – [E(A) + E(B)]

如果ΔH_f为负，说明该化合物在热力学上稳定(比单质混合物能量低)。如果ΔH_f为正但接近零(<0.1 eV/atom)，可能属于亚稳态，需要进一步验证动力学稳定性(声子谱无虚频)。

结构优化完成后，形成能通常呈三组分布：大部分为负(稳定)，少数接近零(亚稳)，少数明显为正(不稳定)。正形成能的直接排除，其余进入下一轮。

动力学稳定性的验证通过声子谱计算完成。项目组使用PHONOPY软件，在2×2×1超胞上用DFPT方法计算力常数。判断标准：声子谱中无虚频(即所有声子频率为正)。部分结构会在Γ点附近出现虚频——这些结构在有限温度下不稳定，排除。

第三层：电子结构计算与拓扑分析

稳定的二维材料进入电子结构计算。这一步的关键参数：

ISMEAR=-5(四面体方法，用于精确的能带和态密度)

NBANDS设为默认值的1.5倍(确保导带有足够空带)

K-mesh加密到12×12×1(二维材料面内需要密k点)

拓扑性质的判断分两步。首先看能带反转：如果价带顶和导带底在Γ点发生了反转(即s轨道特征出现在价带顶，p轨道特征出现在导带底)，说明可能存在能带反转——这是拓扑非平庸的必要条件。少数候选材料显示出能带反转特征。

第二步是计算Z2拓扑不变量。项目组使用VASP的Wannier90接口，先构建紧束缚模型，再通过Wilson loop方法计算Z2.这个过程比较复杂：

用VASP生成投影在原子轨道上的波函数(LWANNIER90=.TRUE.)

用Wannier90做最大局域化Wannier函数(MLWF)

从MLWF构建紧束缚哈密顿量

用Z2Pack或WannierTools计算Wilson loop

11个候选材料中，少数Z2=1(非平庸)，其余Z2=0.Z2=1的材料即为最终筛选结果。

项目组对最终结果做了人工复核：每个结构重新手动优化一遍、声子谱用3×3×1超胞重新计算、Z2不变量用两种独立方法(Wilson loop和parity分析)交叉验证。结果一致，确认筛选结果可靠。

计算参数选择的材料学逻辑

DFT计算材料中，参数选择不是纯技术问题，需要材料学判断。

ENCUT的选择遵循”最大ENMAX×1.3″规则，但有两个例外：含过渡金属的体系，d电子的截断能通常需要更高(部分赝势的ENMAX偏低会导致d带描述不充分);含半金属元素(Bi、Sb、Te)的体系，自旋轨道耦合(SOC)必须开启(LSORBIT=.TRUE.)，否则能带反转特征会消失——这正是项目组在初始筛选中差点漏掉Bi2Se3变体的原因。

K-mesh的选择对金属和半导体不同。金属需要密的k点来正确描述费米面(k点间距<0.02 Å⁻¹)，半导体可以适当放松。但对于拓扑材料的能带计算，即使是非金属性的，也需要密k点——因为Z2不变量的计算本质上是对布里渊区积分，k点不够密会引入数值误差。

反思：高通量筛选中的误差传播

高通量筛选有一个隐性风险：每一步的误差会传播到下一步。预过滤可能误删好材料(假阴性)，结构优化可能因初始猜测不好而收敛到亚稳态，声子谱计算可能因超胞不够大而漏掉边界虚频。

项目组对最终结果做了人工复核：每个结构重新手动优化一遍、声子谱用3×3×1超胞重新计算、Z2不变量用两种独立方法(Wilson loop和parity分析)交叉验证。结果一致，确认筛选结果可靠。

但必须承认，候选组分本身就来源于已有结构的变体——如果真正的全新结构不在这个化学空间内，筛选再高效也找不到。计算材料筛选的天花板不在于计算方法，而于结构空间的覆盖度。更多VASP计算实战经验可以参考VASP/第一性原理栏目，或返回科研学术网首页。