第3章：计算软件与流程：VASP、Quantum ESPRESSO等软件介绍，缺陷计算的标准流程

做缺陷计算这么多年，我经常被问到同一个问题：「到底该用哪个软件？」说实话，没有标准答案。但我会告诉你我的选择逻辑。

目前主流的DFT软件，说白了就两大阵营：VASP和Quantum ESPRESSO（QE）。我两个都用过，各有各的脾气。下面我结合自己的项目经验，给你拆开揉碎了讲讲。

3.1 VASP：工业界的「老大哥」

VASP是我用得最多的软件。为什么？因为它稳定、高效、生态好。我在做氧化物缺陷计算时，几乎首选VASP。

但VASP有个「小毛病」——它是商业软件。我记得刚入行时，导师说「能用VASP就别折腾别的」。嗯，这话现在看依然有道理。

3.2 Quantum ESPRESSO：开源界的「战斗机」

如果你预算有限，或者想完全掌控计算流程，QE是不二之选。我曾在课题组里用QE搭建过一套自动化缺陷计算流程，效果相当不错。

QE的优势在于：

• 完全开源，代码透明
• 支持多种泛函（包括杂化泛函）
• 输入文件灵活，适合二次开发

但说实话，QE的学习曲线比VASP陡一些。我刚开始用QE时，光调k点收敛就折腾了两天。为什么会这样？因为QE对输入参数的敏感度更高，你得理解每个参数背后的物理意义。

3.3 缺陷计算的标准流程

不管用哪个软件，缺陷计算的流程是相通的。我把它总结为三步：结构优化 → 静态计算 → 能量提取。下面我详细说说每一步。

3.3.1 第一步：结构优化

结构优化，说白了就是让原子「找到最舒服的位置」。对于缺陷体系，这一步尤其关键。

我的做法是：

1. 先优化完美超胞（作为参考）
2. 引入缺陷原子，再优化
3. 检查残余力是否小于0.01 eV/Å

你想想看，如果结构都没优化好，后面算出来的能量就是「空中楼阁」。我在项目中遇到过一位同学，他算的缺陷形成能总是偏高，后来发现是结构优化时只用了Gamma点——这显然不够。

EDIFF = 1E-6        # 电子步收敛标准
EDIFFG = -0.01      # 离子步收敛标准（负值表示力）
ISIF = 3            # 优化晶胞参数和原子位置
IBRION = 2          # 共轭梯度算法
NSW = 100           # 最大迭代步数

3.3.2 第二步：静态计算

结构优化完成后，别急着提取能量。你需要做一次静态计算——也就是固定原子位置，只优化电子结构。

为什么要多此一举？因为结构优化时，电子步的收敛标准通常设得比较宽松（比如1E-5）。静态计算时，我会把EDIFF收紧到1E-7甚至1E-8，这样得到的能量才够精确。

我记得有一次，一个学生问我：「老师，我优化完直接提取能量，和再做静态计算差多少？」我让他对比了一下，结果差了0.05 eV——对于缺陷形成能来说，这个误差已经不可接受了。

3.3.3 第三步：能量提取

最后一步，就是从输出文件中提取总能量。这里有个坑：不同软件的输出格式不同，提取方式也不同。

提取能量时，一定要确认单位。VASP默认是eV，QE默认是Ry（1 Ry = 13.6 eV）。我曾经见过有人把Ry当成eV用，结果算出来的形成能差了十几倍——嗯，这种错误犯一次就够了。

3.4 核心流程示意图

下面这张图是我自己画的，把整个缺陷计算流程串起来了。你看一眼就能明白：

3.5 软件选择建议

最后，我根据自己的经验，给你一个简单的选择建议：

• 如果你是新手：从VASP开始。生态好，遇到问题容易找到答案。
• 如果你做高通量计算：考虑QE。开源、灵活，可以写脚本批量跑。
• 如果你做杂化泛函计算：VASP的HSE06实现更成熟，收敛更快。

当然，这只是我的个人看法。你完全可以根据自己的需求选择。记住：软件只是工具，理解物理才是根本。

嗯，这一章的内容就到这里。下一章我们会深入讨论超胞大小对缺陷形成能的影响——这可是个容易踩坑的地方。