声子和电子结构计算

项目过程：石墨烯电子结构与动力学稳定性验证

客户做二维材料器件，需要一份石墨烯的“基础物理档案”：能带结构要有狄拉克点，声子谱要没有虚频，确认材料在计算层面是动力学稳定的。同时他还关心光学支和声学支的数值，用来和拉曼实验对比。

电子结构用 VASP 做 PBE 能带计算，取 Γ–M–K–Γ 高对称路径，K 点网格 12×12×1，截断能 500 eV。声子谱用 density functional perturbation theory（DFPT）在 phonopy 框架下计算，超胞取 4×4×1，位移 0.01 Å。

为了和常见紧束缚模型对比，电子能带我同时画出了 π 和 π* 两条轨道贡献的能带。

案例分析：狄拉克点与六支声子全部正频

第一张案例图是石墨烯的能带结构。π 和 π* 两条带在 K 点相交，形成狄拉克点，交点恰好落在费米能级。这是石墨烯半金属特性的直接体现。M 点和 Γ 点附近的能量色散也基本对称，说明 K 点路径选择合理。

第二张案例图是声子色散。6 支声子分支全部位于正频区，没有任何虚频，说明石墨烯在这个晶格常数下是动力学稳定的。其中 LO/TO E₂g 支在 Γ 点附近频率约 1580 cm⁻¹，对应实验中著名的 G 峰；三支声学支在 Γ 点附近线性趋近于零，符合声学声子的长波极限行为。

总结教训：二维材料计算容易掉进的三个坑

第一，真空层不够会让层间相互作用引入虚频。我测试了 10 Å、15 Å、20 Å，15 Å 以上声子谱没有明显变化，最终用 15 Å。

第二，DFPT 的位移不能太大。0.01 Å 是标准选择，增大到 0.03 Å 会出现个别声学支轻微下凹，容易被误判为不稳定。

第三，PBE 会低估费米速度约 10%。如果客户要定量费米速度，需要用 HSE06 或 GW 修正。

项目经验：能带+声子是最基础也最有效的验证组合

对于任何新材料，先算能带看电子性质，再算声子看动力学稳定性，是性价比最高的两步。石墨烯这个项目虽然简单，但它验证了方法论的可靠性；后面如果换到 h-BN、MoS₂ 或黑磷，流程完全一样，只需要调整赝势和超胞大小。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：声子和电子结构计算。