无定形材料的DFT计算模型：从非晶结构构建到性质预测的挑战

无定形材料（非晶碳、金属玻璃、非晶氧化物、非晶硅等）的DFT计算有一个根本性的困难：没有唯一的最优结构。晶体的结构优化有一个明确的终点——总能量最低的构型就是平衡结构。但非晶材料的势能面上有无数个局部极小值，每个都对应一个”合理的”非晶构型。这意味着单个DFT计算结果不代表”非晶材料的性质”，而只代表”某个特定构型下的性质”。

非晶结构的构建：熔体淬火法

构建非晶结构的通用方法是熔体淬火（melt-quench）。先用分子动力学在高温（远高于熔点）下让体系充分熔化并混合，然后快速冷却（淬火）到目标温度。淬火速率决定了非晶结构偏离平衡态的程度——冷却越快，非晶结构的无序度越高。

非晶碳的计算中，常用的是密度约3.0 g/cm³的sp²/sp³混合碳网络。团队用LAMMPS构建了一个512原子的碳盒子：先在5000 K下跑5 ps的NPT模拟（让原子充分随机化），然后以10 K/ps的速率冷却到300 K（淬火过程约470 ps），最后在300 K下跑50 ps的NVT弛豫。整个MD过程用的是AIREBO势函数——对碳-碳相互作用的描述在密度和成键方面比较可靠。

淬火得到的结构需要做结构分析确认”非晶性”——径向分布函数（RDF）应该没有长程有序峰，只有第一和第二近邻的宽峰。XRD图谱的衍射峰应该弥散而非尖锐。如果RDF出现了第三甚至第四近邻峰，说明淬火速率太慢，体系发生了部分结晶。

DFT精度优化：从MD构型到电子结构

MD构建的非晶结构只是”粗模型”——AIREBO势的精度不足以描述电子结构相关的性质（带隙、态密度、光吸收）。需要在DFT层面做进一步的结构优化和性质计算。

但DFT优化非晶结构面临计算量的问题。一个512原子的碳盒子，用PBE泛函，ENCUT 400 eV，Gamma点（非晶没有周期性，k点采样只需要Gamma点），结构优化可能需要几天时间——因为非晶结构中有很多被”冻住”的高应力键，优化过程中需要大量的离子步来释放这些应力。

团队实际操作中，先用DFT对MD构型做了短时间（约100步离子步）的结构松弛——不是完全优化到收敛，而是把MD阶段积累的明显不合理结构（比如过短的C-C键、过大的键角）修正到物理合理的范围内。然后在松弛后的构型上做静态自洽计算，得到电子密度和态密度。

统计平均：多个构型的意义

前面说过，非晶材料的单个构型不代表”材料性质”。要做有意义的预测，需要构建多个独立的非晶构型（至少5-10个），对每个构型做DFT计算，然后对物理量取统计平均。

团队在这个非晶碳项目中构建了8个独立的512原子构型（每个用不同的随机种子初始化MD）。8个构型的HOMO-LUMO能隙分布范围是0.8-2.1 eV，平均值1.4 eV，标准差0.4 eV。实验上非晶碳（ta-C）的光学带隙约1.2-2.5 eV——计算的平均值和范围都与实验兼容。

如果不做统计平均，只算一个构型，得到的带隙可能是0.8 eV也可能是2.1 eV——跟实验的对比就失去了统计意义。多个构型的平均值和标准差才是非晶材料性质的”正确报告方式”。

金属玻璃的额外挑战

金属玻璃（如Cu-Zr、Pd-Si体系）的非晶结构构建比非晶碳更复杂。金属键没有方向性，原子堆积的局部有序性（短程有序、中程有序）对性质的影响更大。需要更大的模型（1000-2000原子）、更精细的退火协议（包括中间温度的等温退火步骤），以及EAM或MEAM势函数来描述金属-金属相互作用。

DFT处理金属玻璃的计算量更大——需要金属PAW势、更高的k点密度（即使是非晶，k点采样也不能只用Gamma点——因为金属的电子态密度在费米面附近需要更密的采样），以及ISMEAR=1（Methfessel-Paxton展宽）来处理费米面的收敛问题。

无定形材料的DFT计算在VASP/第一性原理方向上属于计算量大、统计要求高的细分领域。它要求研究者不仅能用DFT，还要能合理构建非晶初始结构、做统计平均、正确报告不确定度。更多非晶材料计算的案例，可参阅科研学术网首页。