做第一性原理计算的人都知道，赝势的质量基本决定了整个计算链路的可信度。镍（Ni）这个元素看起来简单，但在PAW_PBE赝势里藏着不少细节：3d电子局域性、磁矩、半芯态选择，哪一步没处理干净，后面晶格常数、声子频率、缺陷形成能都会加倍还回来。这篇文章就用ATOMPAW完整走一遍Ni的PAW_PBE赝势生成流程，从思路拆解、参数取舍、实际操作到常见问题排查，把我踩过的坑和验证过的方法一起讲清楚。
1. 为什么选择ATOMPAW生成Ni的PAW_PBE赝势
1.1 PAW和PBE到底是什么关系
PAW（Projector Augmented Wave）是目前处理电子与离子相互作用的主流方法之一，它通过全电子波函数与赝波函数之间的线性变换，把芯区快速振荡的电子行为转换成平滑、易处理的赝势形式。PBE（Perdew-Burke-Ernzerhof）则是广义梯度近似（GGA）里最常用的交换关联泛函，对绝大多数金属、半导体、氧化物体系都有不错的描述。
把这俩放到一起，就是“用PBE泛函计算全电子参考态，再用PAW变换生成一套可供第一性原理软件直接调用的赝势”。ATOMPAW干的活就是这件事：输入一个原子的基本信息，输出一套完整的PAW数据集，常见格式可以对接VASP等平面波软件。换句话说，你不需要手动去拟合任何一个数值，只需把Ni的电子组态、径向网格、截断半径这些参数告诉它，剩下的匹配工作在程序内部完成。
1.2 Ni的电子结构为什么这么“挑”
Ni的原子序数是28，基态电子排布是[Ar]3d^8 4s^2。问题在于3d轨道径向分布非常局域，和4s轨道在空间上有明显重叠。这种结构对赝势非常敏感：
- 如果直接把3d当芯态处理，磁性完全丢失，Ni的许多性质都没法算。
- 如果把3d当价态处理，又必须保证径向网格能够分辨3d轨道的快速振荡，否则会产生非物理的“鬼态”或数值不稳定。
- Ni是典型的铁磁金属，单个原胞里的本征磁矩大约在0.6 μB左右，这意味着自洽计算必须能稳定打破自旋简并。如果赝势的参考态选择不对，后续算出来磁矩可能直接归零，或者收敛到错误的亚稳态。
所以生成Ni的PAW_PBE赝势，核心就是做两个平衡：一个是保留3d电子细节，另一个是让赝波函数足够光滑、便于平面波展开；另一个是覆盖尽可能完整的物理环境，另一个是控制截断能不要过高。
1.3 自己生成赝势比直接下载标准库好在哪
很多人会直接下载现成的POTCAR，这当然可用，但标准赝势是面向通用场景优化的，它不会为某个具体体系做额外调整。自己用ATOMPAW生成，好处非常实际：
- 可以按需引入半芯态，例如把3p或3s2p一并划入价态，专门应对高压、磁性、X射线近边吸收这类对芯态敏感的计算。
- 可以调整截断半径rc，在不严重牺牲精度的前提下降低平面波截断能，这对大规模分子动力学或表面计算非常友好。
- 可以同时生成多个配置，对比同一元素不同组态之间的差异，搞清楚体系结果对赝势的敏感程度。
- 对刚接触第一性原理计算的人而言，自己从头生成一套赝势并做校验，是理解PAW方法最直观的一条路。
2. ATOMPAW生成Ni赝势的核心思路与参数拆解
2.1 ATOMPAW内部到底做了哪些事
ATOMPAW的工作流程可以拆成四步。第一步，对孤立原子做一次全电子计算，得到参考态的径向波函数和能量本征值。第二步，构造一个光滑的赝波函数，要求它在截断半径rc之外与全电子波函数几乎完全一致，在rc之内则平滑过渡。第三步，根据全电子波函数和赝波函数之间的差异构造PAW投影子和补偿电荷密度。第四步，把整套数据写入目标格式文件，比如VASP的POTCAR格式。
理解这个流程有个很关键的点：PAW不是简单的“赝势”，它保留了一部分全电子信息，所以对电子密度、磁矩、原子力这些量的描述比传统模守恒赝势和超软赝势更细腻。你在输入文件里看到的每一个数字，最终都会影响这套全电子-赝波函数变换的匹配质量。
2.2 价电子组态和半芯态怎么选
生成Ni赝势时，第一个要拍板的问题就是“哪些电子算价态，哪些算芯态”。
| 配置方案 | 价电子范围 | 优点 | 代价 |
|---|---|---|---|
| 标准配置 | 3d^8 4s^2 | 截断能低，计算快，适合常规晶格热力学 | 对强压、磁性反转等场景可能不够 |
| 带3p半芯态 | 3p^6 3d^8 4s^2 | 更准确地处理p-d杂化，适合磁性、高压、光谱 | 需要更高ENCUT，计算成本增加 |
| 全半芯态 | 3s^2 3p^6 3d^8 4s^2 | 最保守，适合极端压缩或芯层激发 | PAW数据集更大，速度最慢 |
我实际生成Ni赝势时，普通金属态计算用标准配置就够了，但一旦涉及高压相变、表面磁性或者Ni基合金中的电荷转移，我会至少把3p纳入价态。原因很直接：3p与3d之间的能级间隔不到100 eV，在压缩环境中3p轨道会和价带发生明显杂化，把它们放在芯态会让不可转移误差变大。
需要注意，原子基态是3d^8 4s^2，但在金属Ni中，由于s-d重排，实际d占据数会升高。因此有人会选用3d^9 4s^1作为参考组态。这没有绝对的对错，更稳妥的做法是两种组态各生成一遍，用目标体系的实验结果做交叉验证。如果算出的磁矩、晶格常数、弹性常数都接近实验值，那这组态就是合适的。
2.3 径向网格、截断半径和能量窗口，三个最容易翻车的参数
径向网格的作用是把波函数从原子核附近到无穷远的空间离散化。网格太疏，对数导数和匹配误差会变大；网格太密，计算量上升但精度提升有限。ATOMPAW里一般使用对数网格或指数网格，只需保证覆盖范围足够大，通常径向上限需要超过价态波函数尾部衰减到可忽略的位置。
截断半径rc是PAW里面最需要斟酌的参数。它不能小于价态波函数内部节点的最外侧位置，否则匹配会产生严重的数值误差；也不能太大，否则赝波函数会偏离真实化学环境，尤其在短程化键合重的体系里会失真。对Ni来说，我常用的rc大致在1.5到2.0 Bohr之间。注意指定的是原子单位，不是埃，很多新手第一次做就在这里看漏导致后面全部跑偏。
能量窗口则决定了程序在多大能量范围内构造投影子。毕竟实际固体中电子分布在一段能带范围内，不是孤立原子那几条能级。窗口太窄，转移性差；窗口太宽，匹配困难。ATOMPAW会在输出日志中给出对数导数的偏差值，这是判断能量窗口是否合理最直接的依据。
3. 实操记录：从零生成Ni的PAW_PBE赝势
3.1 编译ATOMPAW
第一步是准备环境。我习惯在Linux服务器上操作，需要Fortran编译器和基本的BLAS/LAPACK数学库。新版ATOMPAW的部分版本编译时需要libxc支持，如果你输入文件里明确选PBE，而程序又配置了libxc，那么libxc版本也必须匹配，否则会报函数签名错误。
编译流程大概是这样：
tar -xzf atompaw-x.y.z.tar.gz cd atompaw-x.y.z ./configure FC=gfortran make如果configure找不到BLAS库，可以显式指定：
./configure FC=gfortran LDFLAGS="-L/usr/local/lib -lblas -llapack"编译完成后运行：
./atompaw --version确认能正常输出版本号。这一步卡住的人不少，最常见的原因就是Fortran编译器版本太老，不支持代码里用到的某些新语法，升级到gfortran 10以上基本能解决。
3.2 编写Ni的atom.info输入文件
ATOMPAW通过一个输入文件来读入原子参数，我常用的Ni输入骨架包含下面这些关键信息。具体字段顺序在不同版本里可能有调整，建议先对照自带examples目录里的模板。
# ATOMPAW input for Ni, PAW-PBE Ni # 元素符号 28 # 原子序数 Z 3d 8 4s 2 # 价电子组态 PBE # 交换关联泛函 1.70 # 核心半径 rc（原子单位） 0.04 # 径向网格步长 12.0 # 径向网格最大半径这不是一份严格意义上的完整输入文件，而是把最重要的几个参数列出来。实际使用时需要把每组轨道对应一个投影通道，比如s通道、p通道、d通道各自设置截断半径，有的版本还要指定每个通道的参考能量个数。若想复现标准库类似的Ni_pv赝势，就把价电子组态改成3p^6 3d^8 4s^2，同时给p通道单独选一个合适的rc，通常比d通道略大一点。
3.3 运行生成命令
把输入文件保存为atom.info之后执行：
./atompaw --vasp < atom.info有些版本支持的选项不太一样，实在不确定就先用./atompaw -h看帮助。运行结束后，工作目录会出现对应的POTCAR，或者带元素标识的VASP格式数据文件。打开POTCAR看头部，应该能看到类似PAW_PBE和元素名Ni的标记，中间还有Z值、电子质量等元信息。如果文件头部信息不对，说明参数解析或格式选项没匹配上，回头检查版本。
顺便提一句，ATOMPAW产生的结果是否适合特定软件，取决于输出格式是否符合该软件版本的要求。比如VASP不同小版本对POTCAR头的解析可能有细微差异，生成之前先确认你的计算软件用到的是VASP5格式还是更老的格式。
3.4 用Ni的POTCAR做一次完整自洽校验
赝势生成的完成，不代表可以扔进计算目录直接跑。我每次生成完Ni的赝势，都要用面心立方Ni原胞做一轮自洽计算，看四件事：总能随截断能是否收敛、磁矩是否接近0.6 μB、晶格常数是否接近3.52 Å、应力张量是否干净。
一个很简单的收敛性扫描脚本片段：
#!/bin/bash for encut in 350 400 450 500 550 600 do sed -i "s/ENCUT.*/ENCUT = ${encut}/" INCAR mpirun -np 8 vasp_std > log_${encut} echo "ENCUT=${encut} E=$(grep 'E0' OUTCAR | tail -1)" >> encut_scan.dat done扫描完之后画总能-ENCUT曲线，通常到400~500 eV区间差值应小于1 meV/原子。如果到600 eV还在振荡，第一件事不是加截断能，而是回头检查rc是不是取得太小，或者径向网格密度不够。
磁矩检查也很关键。Ni的fcc结构用PBE计算，磁矩大约在0.6 μB附近波动都算正常，但如果算出来是0.0或者0.1这样的值，就要怀疑参考态、自旋初始设置或者半芯态选择出了问题。注意VASP等软件里需要设置ISPIN=2并给一个合理的MAGMOM初始值，否则Ni这种弱磁金属很容易收敛到无磁态。
4. 常见问题与排查技巧实录
4.1 编译时出现未定义引用
最常见的是链接BLAS/LAPACK时少库，报错关键词是dgemm_、dgesv_这类。解决办法是确认系统装了对应的数学库，并在configure时传入库路径。如果不想费劲，也可以在Makefile或编译命令里直接加-lblas -llapack。
还有一次我遇到configure检查通过但make中途报错，后来定位到是libxc版本太新，和程序默认接口约定不兼容。这种情况我会把ATOMPAW换到官方补丁较新的版本，或改用老版本libxc，尽量避免混合版本踩坑。
4.2 输入文件里的组态导致自洽失败
写3d10 4s0或者漏掉某个壳层，很可能在生成过程中出现电子密度不收敛或者能量发散。Ni的孤立原子基态确实存在多种组态竞争，程序日志里会提示自洽迭代是否成功。
我遇到过一种情况：用3d^8 4s^2可以顺利生成，但用3d^9 4s^1时前几步电子密度振荡很剧烈。这时不要急着加迭代步数，先减小径向网格步长，或把初始占据做一个平均化处理，比如设置分数占据让自洽更稳。这也反过来说明，不同参考组态对数值稳定性有明显影响，做高精度计算前，值得花时间把组态敏感性试出来。
4.3 生成的赝势算出的磁矩和实验对不上
如果磁矩明显偏低，最先怀疑半芯态。只把3d^8 4s^2作为价态，p-d杂化被截断，磁性往往被轻微低估。改用3p^6 3d^8 4s^2后，磁矩会回到更合理的区间。另一个容易被忽略的点是rc设得太大，导致d轨道在芯区被压得太平滑，局域磁矩形成不了。
排查方法很简单：同时用标准库的Ni赝势和你生成的赝势跑同一套POSCAR和INCAR，对比总能、磁矩和原子受力。如果差异集中在磁矩，多半就是半芯态或自旋初始值的问题；如果连总能都差异巨大，那要重新审视输入参数和格式兼容性。
4.4 输出文件格式与计算软件版本不兼容
用VASP时，首先要保证POTCAR文件里就一个元素的PAW数据，不要把一个以上的元素粘贴串联弄混。其次是考虑VASP版本对POTCAR头部的解析方式，有时更新软件后旧赝势没有损坏但头部识别异常，看起来直接报“reading POTCAR”相关错误。
我的习惯是在正式跑大任务之前，先用同版本VASP跑一个Ni原子的自洽计算，或者跑一个单胞计算，确认log里读取的原子类型、Z值、磁矩都符合预期。这一步花的时间不超过半小时，但能省下后面几十小时踩格式坑的精力。
5. 从Ni这份示例里延伸出的几点经验
现在数据库确实很发达，生成赝势也用得不算频繁，但自己做一遍之后，对整个PAW框架的理解会深很多。尤其是Ni这种过渡金属，3d电子、磁性、半芯态三个词几乎就是前期最容易反复试错的地方。成年人的学习方式有时候就是这么绕，亲手踩一遍坑比读十遍理论更管用。
如果打算把这一套方法复用到相邻元素，比如Co、Fe、Cu，强烈建议把每个元素的atom.info、生成日志、校验结果按版本管理起来。同一批参数在不同版本的ATOMPAW和不同数学库环境下都可能产生微小差异，有了记录才能追溯结果变化的原因。
还有个小技巧：生成Ni的POTCAR后，我会顺手用一个带Ni原子和真空层的表面模型跑一跑原子力，验证赝势在非均匀环境下的表现。这一步比单纯跑体材料更能暴露短程匹配问题。你在正式计算前如果也能多做这样一个低成本冒烟测试，很多问题都能在项目早期暴露出来。
