企业官网建设流程全解析

1. 项目概览与适用场景

做材料计算这几年，我遇到最多的问题不是“怎么跑程序”，而是“跑完DFT之后带隙和实验对不上”。尤其碰到半导体和绝缘体体系，用PBE算出来的带隙普遍偏低，有的甚至直接给成了金属，这种情况做能带分析、缺陷态评估、光学性质预测，基本等于白算。后来接触到GW+BSE这套组合，才算是把“基态算得差不多、激发态完全没谱”这个老大难问题真正打通了。

VASP是目前第一性原理计算里用得最顺手的平面波软件之一，GW+BSE是它内置的高精度激发态计算方法。简单说，GW负责修正准粒子能级，把DFT那套被低估的带隙拉回实验值附近；BSE负责在GW基础上处理电子-空穴相互作用，也就是激子效应，最后给出可以和实验吸收光谱、损失函数直接对照的数据。这两个方法叠加，是当前计算材料光谱性质比较主流且结果可靠的一整套方案。

这篇文章面向的是已经能跑通VASP常规DFT计算、但对GW+BSE流程还不够熟悉的从业者和研究生。我会从原理、参数、实操到报错排查完整走一遍，尽量把每个参数为什么这么设、哪个步骤容易翻车、哪个文件必须保留这类细节讲透。跟着这套流程走，你至少能少踩一半的坑。

2. 核心原理深挖：GW和BSE分别解决什么问题

2.1 GW近似：把DFT的带隙误差拉回来

要理解GW，先得理解DFT的尴尬。Kohn-Sham DFT的能带本征值严格来说不是准粒子激发能，它的带隙误差主要来源于交换关联泛函对自相互作用的错误处理，以及缺少频率依赖的屏蔽效应。PBE这类泛函把电子间的排斥作用平均化了，导致电子觉得“自己没有被自己排斥”，费米能级附近的占据态被抬得不够高、空态被放得过低，带隙自然就偏小。

GW近似做的事情，是把自能算符写成一个动态屏蔽的交换作用形式，也就是Green函数G和屏蔽库仑作用W的乘积，在里面显式包含频率依赖的屏蔽效应。准粒子方程里面用这个自能算符替换DFT的交换关联势，解出来的本征值才是真正的电子加空穴激发能。实际VASP计算里最常见的是G0W0（也叫一次性GW），也就是用DFT的波函数和本征值固定不动，用它们构造G和W，只算一遍自能修正。G0W0对大多数半导体和绝缘体已经足够准，带隙修正效果非常好，而且比自洽GW便宜很多，是我个人默认选择。

不过G0W0依赖初始的DFT波函数质量，如果体系有强关联电子，比如含有过渡金属d电子或稀土f电子的材料，那么PBE初始波函数本身就不太对，一次性GW修正常常不收敛或者偏离实验很大。这种体系一般需要先做DFT+U或者改用杂化泛函提供更好的初始波函数，实在不行再考虑evGW或QSGW这类自洽方案。VASP里ALGO=GW0对应的是部分自洽GW0，它对本征值做了自洽但固定了W；ALGO=GW则是W也一起自洽，计算量会成倍上涨，普通体系性价比不高。

2.2 BSE：处理电子-空穴相互作用

GW算完，准粒子能级准了，但还缺一个关键物理：光吸收过程中电子被激发到导带之后，它和价带留下的空穴之间是有库仑吸引的，这个被束缚的电子-空穴对就是激子。GW描述的是单粒子激发，BSE才是处理双粒子（电子+空穴）激发的方程，两者物理图像完全不同。

BSE的核心思路是把电子-空穴格林函数的运动方程写成类似两点散射问题的形式，这里面包含三部分：电子和空穴各自传播的准粒子能级项、电子-空穴交换作用（对应亮激子/暗激子的选择规则），以及电子-空穴吸引项（主要由屏蔽库仑作用W构成）。VASP求解BSE的时候，会先用GW步骤得到W和准粒子波函数，然后把激子哈密顿量显式对角化，得到激子本征态和对应的光学跃迁强度，最终输出介电函数虚部，也就是吸收光谱。

为什么不能用DFT的波函数直接做BSE？因为DFT带隙本身就低估了，激子束缚能在很多有机半导体里能到0.5甚至1 eV以上，两个误差叠在一起，光谱整体位置会明显偏离实验。先做GW把带隙和波函数修正了，BSE里的单粒子项才是在一个可靠平台上，这样才能确保算出来的激子束缚能和吸收峰位置接近实验值。我算过几种有机-无机杂化钙钛矿，如果用PBE拿到的能级直接算光学性质，吸收峰普遍偏低0.5 eV以上，GW+BSE之后基本能和实验光谱对到误差0.1 eV以内。这套组合值钱就值钱在这里。

3. 环境准备与VASP编译常见问题

3.1 Ubuntu下VASP编译的前置依赖

VASP收费且License管理严格，这是绕不开的客观条件，我这里默认你已经拿到了合法授权，只聊编译和运行环境的坑。我自己的主力编译环境是Ubuntu 22.04 LTS + Intel oneAPI，这套组合比GCC+OpenBLAS更容易发挥VASP的峰值性能，尤其GW+BSE这种对BLAS库和FFT性能极其敏感的任务，快慢能差出30%以上。

编译VASP之前，有几个前置工具必须备齐：

• Intel oneAPI Base Kit和HPC Kit，里面包含ifort/ifx编译器、MKL数学库和MPI库
• 或者是GCC + OpenMPI + OpenBLAS + FFTW的组合，性能略低但在纯AMD平台上反而更稳定
• GNU Make和CMake取决于你用的VASP版本，老版本（5.4.4）一般用make，新版本（6.3+）也可以走CMake

实际编译时最容易翻车的地方是MKL库接口不匹配。VASP 5.4.4默认用ifort编译，但在新版oneAPI里ifort已经被ifx替代，很多旧makefile里的MKL链接写法在新版本编译器中会提示找不到mkl_blas95_ilp64这类库。我的解决方法是直接改用6.x系列的CMake构建流程，它会自动探测MKL接口，比手工改makefile省心得多。

3.2 新增配置与硬件估算

如果你用的是VASP 6.x，编译前需要确认把-Dtbdyn、-Dscalapack这些选项打开，GW和BSE本身不需要额外特殊的宏，但必须在编译时开启了MPI和SCALAPACK支持，否则并行跑大规模GW任务会非常吃力。我自己遇到过一个情况是集群管理员只编译了串行版VASP，BSE计算时内存直接爆炸，因为BSE激子哈密顿矩阵规模是(空带数 x k点数)的平方级别，串行下根本装不下。

硬件方面，GW+BSE和普通DFT对资源的需求差别很大。普通DFT结构优化，32核加64GB内存就能跑得像模像样。到了GW步骤，需要存储大量空带波函数和频率依赖的屏蔽函数，通常建议至少128GB内存起步；BSE步骤如果体系稍大（比如超胞包含50个原子以上、k点6x6x6左右），256GB内存是及格线，对存储IO的吞吐要求也非常高。我见过不少人在工作站上跑GW+BSE，最后不是让CPU等计算，而是让计算等IO，所以强烈建议用NVMe SSD或者并行文件系统来存放WAVEDER和WAVECAR这类大文件。

4. 完整实操流程：从结构优化到BSE光谱

4.1 第零步：结构优化和基态DFT计算

GW+BSE是一个计算链，不是单步任务。我在正式跑之前，一定会先把前两步基础工作做扎实，否则后面全白费。

第一步是结构优化。对于晶体体系，用PBE泛函做原子位置和晶胞参数的弛豫基本够用，但需要注意加适当的k点密度和截断能。我的经验是体系的ENMAX来自POTCAR文件，但结构优化时ENCUT要取ENMAX的1.2到1.3倍，这样做是为了避免由于截断能不足导致的基态波函数人为畸变。晶胞优化要用ISIF=3，原子位置优化用ISIF=2，在GW之前晶胞参数最好用实验值或者更高精度的泛函校正过，因为GW步骤对晶格常数的敏感性远大于普通DFT。

第二步是基态静态计算，这一步的输出是GW的输入基础。有几个关键点：

• 使用与后续GW一致的ENCUT，不要中途改动
• 采用较密的k点网格，GW和BSE的k点必须比一般的静态计算更密，通常建议至少保证单胞4x4x4以上
• 必须打开LWAVE=.TRUE.，保存WAVECAR文件
• 建议写入WAVEDER文件，这个文件包含了占据态和空态之间的动量矩阵元，GW步骤计算介电函数时离不开它

如果你在旧版本VASP里跑过，可能会注意到WAVEDER通常是和BSE强绑定，实际上GW也需要它，它描述的是能带间光学跃迁矩阵元在倒空间的变化行为。这个文件生成一次要花不少时间，一定记得用LOPTICS=.TRUE.提前生成并保留好。

4.2 第一步：GW0计算的核心参数

基态算完之后，进入GW0步骤。INCAR里最核心的设置是下面这些：

• ALGO=GW0：表示采用GW0近似，即自能计算中更新准粒子本征值但不自洽更新屏蔽作用W，比全自洽GW便宜且数值稳定
• NELM=200：GW自洽循环的外层电子迭代步数上限，GW0一般100到200步足够
• NBANDS：GW对空带数量极其敏感，需要从默认值往上加，通常至少取“价带数+导带数”的2到3倍，具体需要在后续收敛性测试中确认
• NOMEGA=50：频率网格点数，默认值往往不够，尤其计算能量损失谱时建议增加到80左右

很多人第一次跑GW0会踩同一个坑：直接在静态计算的INCAR文件上改ALGO=GW0，但忘了检查NBANDS。DFT静态计算默认的空带数极少，对于半导体大约只有十几个导带，GW需要几十甚至上百个空带才能把屏蔽作用W描述到收敛。NBANDS不足时，GW自能高频部分严重缺失，带隙结果会偏低且不随空带数稳定。判断标准很简单：把NBANDS翻倍，如果带隙变化超过0.1 eV，说明还没收敛，继续加，直到变化量小于0.05 eV为止。

GW0跑完后，OUTCAR里会输出每个k点、每个能带的准粒子修正量，同时会在vasprun.xml里写入修正后的能带。我一般习惯把GW0和BSE分两步跑，而不是用ALGO=BSE一步到底，这样可以在中间检查GW的准粒子带隙是否合理，如果这一步已经严重偏离实验预期，再往下跑BSE就没意义了，先得回去检查初始波函数。

4.3 第二步：BSE计算核心参数

GW0收敛后，进入真正吃硬件的BSE步骤。保持已有的WAVECAR、WAVEDER文件不变，INCAR改成BSE模式：

• ALGO=BSE：让VASP进入BSE求解流程
• NBANDS：这是BSE激发态空带数，BSE需要比GW更多的空带，因为它要建激子哈密顿量，空带太少会把激子束缚在高能态，压制吸收峰
• NBSE：指定BSE实际使用的能带数量，通常和NBANDS一致，但可以从NBANDS里截取
• LOPTICS=.TRUE.：开启光学性质计算，在BSE里必须配合WAVEDER
• NEDOS=2000：介电函数能量网格点数，光谱想要平滑就设大点
• CSHIFT=0.1：展宽系数，单位是eV，用于给光谱加洛伦兹展宽，取值太小光谱会像梳子一样离散，太大则过度平滑模糊了峰位，我一般用0.1到0.2

BSE计算可调参数里很容易被忽略的是LSPECTRAL=.TRUE.。当你设置ALGO=BSE时，打开LSPECTRAL可以让VASP先进行一次无相互作用的虚频介电函数扫描，这样得到的无相互作用光谱可以当作参考线，便于分析激子束缚能（即相互作用光谱峰值和无相互作用光谱峰值的能量差）。这个量在讨论激子效应强弱时非常有用，强烈建议打开。

4.4 文件保留与中间产物管理

整套流程里会产生几个超大的文件：WAVECAR动辄几十到数百GB，WAVEDER也不小，中途如果磁盘紧张，很多人会手滑删掉WAVEDER。这个文件是BSE输入的必需品，删了就必须重跑GW0重新生成，代价极大。我的建议是给每个步骤建立独立目录，保留上一层的WAVCAR和WAVEDER作为下一层的输入，并做软链接而不是把文件复制来复制去。

还有个经验是GW0和BSE都建议在同一个k点网格下运行。很少有人会注意到VASP对BSE计算中有无空带截断与k点对称性的要求：BSE对k点的对称性处理会利用晶格点群，如果k点网格非均匀或者用了太粗糙的网格，BSE激子色散会被严重扭曲。从GW到BSE，k点网格必须保持一致，这一点比带宽收敛还重要。

5. 关键参数表与收敛性测试策略

5.1 一套可直接套用的INCAR模板

下面给出一套我在实际项目里验证过、可以直接修改使用的INCAR参数模板。这里以单胞半导体体系为例，超胞或分子体系需按需调整NBANDS和k点。

# 基态静态计算 SYSTEM = my_system ENCUT = 1.3 * ENMAX ISMEAR = 0 SIGMA = 0.05 LREAL = .FALSE. LWAVE = .TRUE. LCHARG = .TRUE. LOPTICS = .TRUE. NEDOS = 2000 # GW0步骤 ALGO = GW0 NELM = 200 NBANDS = 128 NOMEGA = 60 LSPECTRAL = .FALSE. # BSE步骤 ALGO = BSE NBANDS = 160 NBSE = 120 NEDOS = 2000 CSHIFT = 0.1 LOPTICS = .TRUE. LSPECTRAL = .TRUE.

NBANDS=128和NBSE=120这个设定是针对大概20个原子以下的半导体单胞，价带数量在10个左右。如果你的体系原子更多、电子更多，128这个数远远不够，需要按收敛测试结果来调整，别直接复制套用。

5.2 收敛性测试的正确顺序

GW+BSE的计算成本很高，收敛性测试不可能每个参数都做一遍完整的矩阵扫描，得有优先级。我自己的测试顺序是：

• 先固定k点和ENCUT，扫描NBANDS。这是最重要的一步，因为空带数对GW带隙和BSE光谱位置的影响最大
• NBANDS收敛后，再测k点密度是否需要增加。这一步成本也很高，但至少不用担心空带不够造成的假象
• 最后再测NOMEGA和CSHIFT这类相对次要的参数

实际操作中，还有一个经常被忽略的步骤：用较便宜的设置（比如较少的k点、较大的展宽CSHIFT）先跑一遍全流程，确认输入文件和流程没毛病，再上完整参数全面跑。我第一次跑BSE就吃了这个亏，直接上全参数，结果静算一步就卡了两天，发现INCAR里忘了关掉结构优化相关设置，白白浪费计算资源。

5.3 不同体系的计算成本分级

GW+BSE的计算成本与体系尺寸、维度、k点数目关系很大，这里给一个粗略的成本分级参考：

这个表只是经验值，实际耗时受硬件、版本和参数影响很大。但有一点是共通的：如果你发现BSE比GW0还便宜很多，多半是哪里出了问题——BSE要显式对角化一个巨大的激子哈密顿矩阵，计算量只增不减。

6. 常见报错与排查技巧实录

6.1 报错速查表

我在VASP的GW+BSE实战中收集了不少报错和问题，这里挑出频率最高、最影响后续运行的几个：

6.2 排查思路与心得

遇到问题时我的第一反应不是改参数，而是回到中间文件检视整个计算链。GW步骤看OUTCAR中的准粒子修正量是否合理，特别关注VBM和CBM处的自能修正大小。如果VBM修正量异常地大（超过2 eV），基本可以断定是初始基态波函数问题或者空带不足。BSE步骤出问题，则优先检查WAVEDER有没有正确加载、KPOINTS是否和GW一致。

还有一个细节，VASP不同版本之间GW和BSE的内部实现有差异。比如5.4.4版本里BSE的NBANDS不能超过NBANDSGW，6.x版本则引入NBSE来显式控制BSE的能带选取。升级VASP版本后如果沿用旧INCAR，轻则警告，重则算错结果。我的习惯是每换一个新版本，先用一个已知基准体系跑一遍全流程，对比输出光谱和已知文献数据，确认没问题再上实际项目。

6.3 关于数值噪声与物理判据

计算光谱很容易陷入“数值上跑出来了，物理上不知道对不对”的状态。我自己的判据有两条：一是GW带隙与实验值差距是否在0.1-0.2 eV以内，二是BSE吸收峰与实验光谱峰位差距是否在0.3 eV以内。如果都满足，说明参数收敛、物理合理；如果GW已经准了但BSE峰位还是偏，那大概率要审视激子束缚能讨论：束缚能本质是相互作用谱和无相互作用谱的峰位差，用LSPECTRAL=.TRUE.输出两条曲线，差值就是激子束缚能，和实验激子结合能对比一下，偏差太大就要警惕是不是屏蔽W算得不准、空带截断不够，或者k点太稀。

还有一个老生常谈但值得再提的点：不要把CSHIFT黑影当成物理峰。展宽过大的光谱会把两个本应分离的激子峰糊成一个峰，有时候还会把连续谱边缘抬起来被误认为新峰。对比不同CSHIFT取值下的峰位是否稳定，是一个低成本的验证办法。

7. 实操中的几个个人经验与扩展建议

GW+BSE这套流程，我在实战里反复调、反复踩坑，积累了几个百试不爽的小经验，这里一并分享。

第一个是关于WAVEDER和WAVECAR的保存策略。磁盘空间吃紧时，很多人会在GW0结束后删除WAVEDER（因为还没到BSE步骤），等BSE步骤报缺文件才追悔莫及。我现在的做法是跑完基态静态计算后立刻把WAVEDER做一次压缩归档，放在独立的存储目录，然后建立软链接给后续步骤用。这套习惯救了我很多次，尤其是中途有人误清临时目录的场景。

第二个是关于BSE内存用量的预估。BSE激子哈密顿矩阵的存储规模大约是(NSEED^2)级别，NSEED是BSE实际使用的价带空带组合数。粗略估算时可以这样算：假如你设置NBSE=120价带空带各选了一部分，最终的矩阵规模可能在几十万乘几十万，存储量轻轻松松超过几十GB。如果内存不足但不想砍k点，还可以考虑把部分空带态用能带投影方法降维，这个在VASP 6.3之后有更多的控制参数支持。

第三个是后续扩展。GW+BSE不只是能做吸收光谱，还能扩展到圆二色谱、反射率、能量损失谱、介电函数实部虚部、甚至非线性光学。VASP的BSE模块输出的是介电函数张量，只要把不同极化方向写清楚，就能分析光学各向异性。我最近在算一个低对称性有机晶体的偏振吸收谱，BSE出来的三组对角分量差异非常明显，这类信息对实验上的偏振光吸收测量有直接参考价值。

算完一套GW+BSE，我的原则是永远保留完整的输入输出档案和版本信息。这套计算非常贵，算完不把数据整理好、不把中间过程记录清楚，之后写论文或者复现结果时会非常痛苦。用表格记录每步的INCAR、版本、时间和关键结果，长期来看是性价比最高的一笔投资。