从结构到性质：VASP 进阶任务

体相计算让晶体在三个方向无限重复；真实表面却会终止。常用 slab 模型把若干原子层与真空放进同一个周期盒子，得到“薄片—真空—薄片”的重复体系。它并不是一块半无限晶体，层数、真空和侧向重复尺寸都会成为模型误差。建议先完成硅的自洽与优化练习，再进入本节。

建立一个可比较的模型家族

1. 确定材料和晶面。先用同一泛函优化体相，再切出所需 (hkl) 晶面。记录表面终止、表面化学计量、侧向超胞、层数、真空厚度。极性表面可能重构、改变终止或需要补偿，不能仅靠加大真空解决。
2. 将法向对准第三晶格矢量。对入门模型，让 c 与表面垂直，便于使用一维真空、k 网格和偶极校正。不要把倾斜晶胞里的 Cartesian z 和第三晶格方向不加区分地混用。
3. 固定与放松要服务于问题。例如固定底部若干层模拟体相支撑，放松顶部与吸附物。固定层数是需检验的模型选择，不能认定“底两层”对所有表面都成立；对称 slab 可采用中心层固定、两面相同终止。
4. 检查周期相邻原子。吸附分子可能从边界伸向邻胞；看单个盒子不重叠并不足够。将结构沿面内方向重复显示，检查最短距离、分子方向和覆盖度。

POSCAR 的 Selective dynamics 行启用逐原子自由度：F F F 固定，T T T 放松。通常的离子优化中，这些标志对应直接晶格坐标方向，即使位置用 Cartesian 写入也一样。斜晶胞里想“只固定垂直高度”，必须先弄清晶格方向。官方 POSCAR 定义。

中性、非磁 slab 的优化起点

下列模板不含实际结构、k 网格或 POTCAR；ENCUT 与中心坐标必须替换。若研究磁性过渡金属、DFT+U 或分子弱吸附，还须单独验证相应设置。ISMEAR=0 是初始占据方案，金属表面须测试展宽和能量差收敛。

# TEMPLATE: neutral, nonmagnetic slab; normal parallel to c
# Supply POSCAR, converged KPOINTS and licensed POTCAR.
# Replace ENCUT; add validated magnetic/U/vdW settings when needed.
SYSTEM = neutral-slab-relax
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-6
NELM = 120
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.05
IBRION = 2
NSW = 150
ISIF = 2
EDIFFG = -0.02
POTIM = 0.3
ISYM = 0
LDIPOL = .TRUE.
IDIPOL = 3
DIPOL = REPLACE_WITH_SLAB_CENTER_IN_FRACTIONAL_COORDINATES
LWAVE = .FALSE.
LCHARG = .FALSE.

ISIF=2 固定整个晶胞，避免把真空当成应被压缩的体积；原子仍按选择性约束移动。EDIFFG=-0.02 是 0.02 eV/Å 的示范性力阈值，最终比较很小的吸附能差时应进一步收紧。先用不带偶极势的自洽结果预收敛，再启用校正，常有助于困难体系；这不是跳过最终校正的理由。

对于有净偶极矩的中性非对称 slab，IDIPOL=3 指定沿第三晶格方向的能量校正，LDIPOL=.TRUE. 同时修正势与力；DIPOL 用分数坐标给出靠近 slab 质心的参考中心。不要机械填 0.5 0.5 0.5，应先把薄片放在对应位置。对称模型也应检查是否真的没有净偶极。参见 IDIPOL、LDIPOL 与 DIPOL。

当真空沿 c 方向且不同薄片已充分解耦时，常采用 Γ 中心 Nₓ × Nᵧ × 1 的 slab k 网格；面内 Nₓ、Nᵧ 要随超胞与倒格矢长度重新收敛。更大的面内超胞通常允许更稀的网格，但相同整数网格不代表相同采样密度。检查真空平台时，在最终静态任务用 LVHAR=.TRUE. 输出离子加 Hartree 电势，再做法向平面平均；该势不包含交换关联项。LVHAR。

吸附能不是只减两个数

本教程约定：相对所选参考态，负值表示电子能层面的放热吸附。如果文献把“结合能”定义成相反符号，必须先统一定义再比较。对于 H 吸附，可选 ½E(H₂) 作为每个 H 的参考；选孤立 H 原子会得到不同数值，不能混称同一吸附能。若盒子中吸附了 n 个 X，报告总吸附能还是每个 X 的平均值也须说明。

1. 分别优化清洁表面、多个吸附位点/取向、独立参考分子。气相分子使用足够大盒子和 Γ 点，并确认自旋态；O₂ 等开壳层分子不能按默认非磁态处理。
2. 以一致泛函、POTCAR 数据集、ENCUT、面内 k 点密度和能量提取定义做最终静态计算。slab 与 slab+X 保持相同侧向晶胞和固定层策略；气相参考应单独做盒子收敛。
3. 不要只测试总能；逐一增加层数、真空、侧向面积、k 网格，看吸附能与关键键长是否达到研究要求。真空可从约 15 Å 起步测试，不能把它当普适合格线。
4. 如比较实验温压、电极电势或吸附自由能，另加一致的零点能、热贡献、熵、溶剂与电化学参考。仅有 ΔE 不等于 ΔG，更不直接给出反应速率。

空位、间隙与替位缺陷会改变局域键合、电荷俘获和离子迁移。一个含空位的周期超胞代表有周期缺陷阵列的材料，并不自动等于无限稀释的孤立缺陷。首先用最简单的中性空位练习，再讨论带电缺陷。

从中性硅空位开始

1. 由充分优化的体相构造多个大小的超胞。为每个大小计算完整母体参考，不要把不同数值精度的原胞能量直接拼接。
2. 删除一个 Si 原子，更新元素计数与坐标行数；电子总数默认也随原子数改变。围绕空位给几组轻微扰动，避免高对称初态把结构锁在局部驻点。
3. 通常固定已优化母体的晶格、只放松原子，即 ISIF=2；这对应稀缺陷处于体相晶格环境。若研究有限浓度或缺陷引起的宏观体积变化，需要另行定义晶格放松方案。
4. 用一致静态设置取得 E_defect 与 E_host，同时记录局域几何、自旋、是否出现新的占据态，不能只保存一个总能量。

μ_Si 是每个 Si 原子返回参考库时的化学势；纯硅在零温近似中常用一致计算的体相每原子能量。这里的“+μ”来自从超胞拿走一个原子。更换成孤立原子能量会改变物理参考，不是一个数值修补办法。

带电形成能中每一项是什么意思？

这里 n_i>0 表示向母体加入第 i 种原子，移走则为负；q>0 表示从中性的该缺陷超胞移走电子。E_F 从母体价带顶起算，E_VBM 需处于一致电势基准。本式把所选方法的电势对齐与有限尺寸处理统一放入 E_corr；若另一种写法显式加入 qΔV，就必须检查是否与 E_corr 重复计数。

不同电荷态应分别优化，形成能随 E_F 的斜率是 q，最低包络的交点给出热力学电荷转变位置。光学跃迁是在不同结构约束下讨论的垂直过程，不能拿热力学交点直接当光谱峰。标准形成能框架与稀释极限问题可从 Freysoldt 等的缺陷综述继续学习。

# TEMPLATE: charged bulk defect at fixed host lattice
# Not a charged-slab template. Supply the actual supercell and references.
# Replace ENCUT and NELECT. Set a justified spin state / MAGMOM if required.
SYSTEM = bulk-defect-relax
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-7
NELM = 150
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.05
NELECT = REPLACE_WITH_NEUTRAL_DEFECT_ELECTRONS_MINUS_Q
IBRION = 2
NSW = 200
ISIF = 2
EDIFFG = -0.01
ISYM = 0
LWAVE = .FALSE.
LCHARG = .TRUE.
# No universal correction tag: specify and validate correction externally.

NELECT 应由当前缺陷结构和 POTCAR 的价电子数求出中性 N₀，再取 N₀−q；不能用未删原子的母体电子数来判断 q。中性计算通常删除 NELECT 行，让程序按元素计数确定。非中性周期计算引入补偿背景，程序能够运行不代表带电总能可直接比较。参见 NELECT 官方说明。

优化结束后将该电荷态的 CONTCAR 复制到单独静态目录，以 IBRION=-1; NSW=0 获取最终能量。如校正流程需要电势，可在缺陷和母体静态任务中一致设置 LVHAR=.TRUE. 并保存 LOCPOT；同时保证用于比较的网格、势定义和参考区域匹配。仅保存 CHGCAR 不能替代所选校正工具要求的电势数据。官方电势输出定义。

三个互相独立的误差来源

• 化学势范围：化合物的“富 A / 贫 A”条件受到母相与竞争相稳定性约束；不能同时把所有元素取成任意最高化学势。先列质量守恒和不析出竞争相的条件。
• 有限尺寸：增大超胞、检查缺陷电荷是否局域，并采用适用于晶胞形状和介电张量的校正。FNV 及其扩展各有假设，软件输出校正数字后仍需比较尺寸趋势。参见 各向异性与弛豫相关的电荷校正研究。
• 电子结构：PBE 的带隙与局域化误差会影响电荷态排序。换 HSE 或 U 值后要一致重算相关参考，并说明改变的是近似方法，不只是“把能带画得更好看”。

优化得到“力很小”的结构，只说明一阶导数接近零。它可能是谷底，也可能是鞍点。谐振声子研究平衡位置附近的二阶导数，把原子集体振动写成波矢 q 与模式支上的频率。它能讨论动力学稳定性、振动自由能和热容，但仅有谐声子不能给出真实热导率、强非谐相变或所有有限温度效应。

有限位移方法从“推一下后，各原子受力怎样变化”求力常数。位移太小会放大力噪声，太大则偏离谐区间；所以电子收敛、结构残余力、位移幅度和超胞大小都需检验。

路线 A：VASP 自己生成位移

IBRION=5 不借助对称性减少位移，IBRION=6 利用对称性；NFREE=2 采用正负中心差分，POTIM=0.015 在这里是 Å，绝不是 MD 的时间步。下面只求输入晶胞 Γ 点的振动。若输入是超胞，结果是该超胞 Γ 点，不是已经完成原胞整个布里渊区的色散。

# TEMPLATE: internal finite-displacement Gamma modes of the INPUT cell
# Start from a tightly optimized, nonmagnetic semiconductor/insulator.
# Supply POSCAR / KPOINTS / POTCAR and replace ENCUT.
SYSTEM = gamma-phonon-internal
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-8
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.01
IBRION = 6
NFREE = 2
POTIM = 0.015
NSW = 1
ISIF = 2
LWAVE = .FALSE.
LCHARG = .FALSE.

在 OUTCAR 检查本征模式及频率，保留原子质量和本征矢。不要对本该完整振动的晶体随意固定原子；局域吸附物振动的部分 Hessian 是另一种近似，且内置有限位移的选择性约束支持与方向定义和普通优化不同，当前官方文档说明仅 IBRION=5 支持该约束，作用在 Cartesian Hessian 分量。有限位移官方说明、NFREE。

路线 B：Phonopy 生成超胞，VASP 只算力

这条路线易于并行：每个位移结构是独立静态任务。Phonopy 把多组力重建为力常数，再作傅里叶变换得到任意 q 处的动力学矩阵。关键是保留位移文件与输出的一一对应，不能重新排序原子、重新标准化晶胞，或把有位移的结构再次优化掉。

# TEMPLATE: one externally displaced Phonopy supercell, SINGLE POINT ONLY
# Do not relax POSCAR-xxx. Use identical settings for every displacement.
SYSTEM = phonopy-force
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-8
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.01
IBRION = -1
NSW = 0
ISYM = 0
LWAVE = .FALSE.
LCHARG = .FALSE.

# TEMPLATE: Phonopy 4.5-style workflow. Check your installed --version/--help.
# Need optimized POSCAR-unitcell; this does not run VASP or supply POTCAR.
phonopy-init -d --dim 2 2 2 --pa auto -c POSCAR-unitcell
# Copy each generated POSCAR-xxx to its own disp-xxx/POSCAR.
# Run STATIC VASP forces with the supplied INCAR-force.template there.
# Replace this sample list with ALL files, in phonopy_disp.yaml order:
phonopy-init --sp -f disp-001/vasprun.xml disp-002/vasprun.xml
# --band auto requires seekpath in the SAME active Python environment.
# Check: python -c "import seekpath"
# If missing, follow the official install guide; e.g. conda install -c conda-forge seekpath
# When phonopy_params.yaml has been created successfully:
phonopy --band auto -p phonopy_params.yaml
# Example post-processing q mesh; independently converge it:
phonopy --mesh 20 20 20 -t phonopy_params.yaml
# Older releases: see the manual; setup may use phonopy -d / phonopy -f,
# and postprocessing may use phonopy-load. Do not mix CLI generations.

截至本教程核验的 Phonopy 4.5 文档，准备与收集阶段使用 phonopy-init，后处理使用 phonopy；旧版本常把准备阶段放在 phonopy，后处理可能用 phonopy-load。先记录安装版本和 --help，不要把两代命令拼接。命令迁移说明、VASP 接口、后处理命令。

--band auto 还要求同一 Python 环境中可用 seekpath；安装了 Phonopy 并不代表每种安装方式都已包含这个可选依赖。先用 python -c "import seekpath" 检查，缺失时按官方安装说明补齐，再生成自动路径；自行指定路径也要与实际原胞的倒易基矢一致。

极性绝缘体的长程偶极相互作用还会影响 Γ 附近 LO–TO 分裂。需要同一方法计算的 Born 有效电荷与电子介电张量，按 Phonopy 的非解析项修正要求提供 BORN 数据；不能从另一个泛函或无关文献随意借值。Phonopy 输入文件与 BORN。

看见“负频率”后怎么做？

1. 先区分绘图约定：负轴常表示 ω²<0 的虚频，并不是“振动向负方向”。
2. 检查结构残余力、应力和磁态，增加 ENCUT/k 点、收紧 EDIFF，测试位移幅度与超胞。Γ 附近很小的声学偏移可能含数值或转动不变性误差，但不能一律删掉。
3. 观察对应本征矢，把结构沿软模正负方向做小扰动后再优化，寻找可能的低能畸变；波矢非 Γ 时必须用容纳该畸变的超胞。
4. 若收敛后仍存在显著虚频，应报告“在当前近似与结构条件下不具谐动力学稳定性”。若声称温度稳定了结构，还需要有限温度/非谐证据。

二维材料的面外声学支在 Γ 附近尤其敏感，不能仅用一张粗网格声子图下结论。官方虚频排查与软模操作。

弹性常数描述晶体对微小均匀形变的抵抗；从数学上看，它是能量关于应变的二阶导数。它回答“沿某方向拉伸多难”，与原子局域振动、塑性变形、断裂强度并非同一物理量。计算前先明确对象是无外应力体相、指定压力下的晶体，还是含真空的二维薄层。

两种计算思路

一种做法对晶胞施加若干独立、正负小应变，固定变形后的晶格，计算应力或能量并拟合；另一种使用 VASP 内置有限差分响应。前者便于看到线性区间，后者便于获取完整张量。无论哪条路线，钳位离子与允许内部原子位置响应的松弛离子弹性常数都要分清：后者通常更接近低频弹性测量。

1. 用适合应力精度的 ENCUT 和 k 网格，充分优化参考结构，确认残余应力相对目标精度很小。能量已收敛不代表应力一定收敛。
2. 内置路线以 IBRION=6 配合 ISIF=3 计算有限差分与内部应变响应。这里是在已优化结构附近求导，不是再做一次普通 CG 结构优化。IBRION 官方说明。
3. 显式应变路线可选几组正负小应变测试，例如 0.3%、0.5%、1%，这些只是测试起点。求松弛离子响应时在每个变形晶胞中用 ISIF=2 优化内部坐标；不能把晶胞也放松回原样。
4. 比较不同应变范围与计算精度下拟合斜率是否稳定；过小应变受噪声影响，过大应变混入非线性。不要仅取两个点而不给不确定性。

# TEMPLATE: internally relaxed elastic tensor of a 3D BULK solid
# Input must already be tightly relaxed at the target pressure.
# Not suitable for directly reporting a vacuum-containing slab in GPa.
SYSTEM = bulk-elastic-finite-difference
GGA = PE
ENCUT = REPLACE_WITH_STRESS_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-8
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.05
IBRION = 6
ISIF = 3
NFREE = 2
POTIM = 0.015
NSW = 1
LWAVE = .FALSE.
LCHARG = .FALSE.

读 OUTCAR 时最常见的三个误会

• kBar 不是 GPa。TOTAL ELASTIC MODULI (kBar) 要乘 0.1 才得到 GPa。不要把输出数字直接加上 GPa 标签。
• 看表头的剪切顺序。VASP 常列 XX YY ZZ XY YZ ZX，与常见 Voigt xx yy zz yz xz xy 不同。整理矩阵时行列都需同步映射，再使用对应应变约定；不能只把最后三个标签改名。
• 应力符号有约定。VASP 的对角正应力对应压缩倾向。自己用“拉伸为正”的连续介质定义拟合时应先转换，否则可能得到整体负的 C。ISIF 的应力定义。

内置输出可能分列刚性离子贡献、内部弛豫贡献和总弹性模量；报告时要注明选了哪一张表。官方输出示例可见 弹性矩阵列序讨论。对于含真空二维材料，GPa 会随人为盒子高度改变；如报告面内二维刚度，应使用清楚的厚度约定或换算成 N/m，并交代有效高度。

机械稳定不等于热力学最稳定

这些条件只适用于相应对称性和无外应力条件。更一般地，应检查与采用的应变约定一致的弹性能二次型是否正定；有有限外压时需用相应修正后的稳定性条件。不能把立方公式套给所有晶系，也不能只看对角元为正。Mouhat 与 Coudert 的稳定性判据。

两种稳定构型之间的能量差告诉我们反应热，却不告诉我们需要跨越多高的势垒。NEB 在已知初末态之间放置若干中间构型，沿路径维持间距、在垂直方向寻找低能通道。每个“image”是一份完整原子结构，不是图片文件。它寻找给定初始路径附近的最小能量路径，因此初猜与可能的替代机制很重要。

目录与结构必须一一对应

neb-native/
  INCAR       # rename a completed INCAR.template
  KPOINTS     # same converged sampling for all images
  POTCAR      # supplied locally by a licensed user
  00/POSCAR   # relaxed initial state
  01/POSCAR   # intermediate image 1
  02/POSCAR
  03/POSCAR
  04/POSCAR
  05/POSCAR   # intermediate image 5
  06/POSCAR   # relaxed final state

IMAGES=5 表示五个中间构型，所以共有七个结构目录。端点应事先独立充分优化并保存一致精度的能量；VASP 在父目录启动。不同 image 必须具有相同晶格、原子数、元素顺序、原子身份顺序和约束。仅元素总数相同还不够：如果插值时把原子 A 与另一个等价原子 B 交换了编号，路径可能穿过极高能重叠结构。

1. 选择初末态，分别确认力收敛和局部稳定。对于扩散，先考虑周期最短路径；直接对 0.99 与 0.01 分数坐标插值可能绕盒子一圈。
2. 构造少量中间结构并逐帧看动画，排除原子碰撞、整体平移和不必要转动。复杂反应可用更好的几何插值或分段机制；增加 image 不是修复错误原子对应的方法。
3. 先预收敛普通路径，再提高精度或增加关键区间采样。提交 MPI 任务时按集群规则让 image 并行划分合理，总进程数通常需可按 image 数分组；不要照抄别人的总核数。
4. 逐 image 检查电子收敛和 NEB 投影力；达到 NSW 上限不等于路径已经收敛。可继续运行，但需保存上一轮结果并让每个中间 CONTCAR 成为各自的新 POSCAR。

目录约定和固定晶胞限制参见 原生 NEB与 IMAGES。SPRING=-5 的负号是启用 NEB 形式的开关约定，不表示物理上的负刚度；参见 SPRING。

路线 A：原生 VASP 的普通 NEB

# TEMPLATE: NATIVE VASP NEB, fixed common cell; five intermediate images
# Supply 00/POSCAR through 06/POSCAR, parent KPOINTS and licensed POTCAR.
# Endpoints must be separately relaxed and calculated with matching settings.
SYSTEM = native-neb
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-7
NELM = 150
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.05
IMAGES = 5
SPRING = -5
IBRION = 1
NFREE = 2
POTIM = 0.3
ISIF = 2
NSW = 250
EDIFFG = -0.03
ISYM = 0
LWAVE = .FALSE.
LCHARG = .FALSE.
# Native template: NO LCLIMB, NO IOPT, and NOT POTIM=0.

此处 IBRION=1 与较短历史 NFREE=2 是预收敛起点，POTIM=0.3 是有限离子步长尺度，需按路径稳定性调整。NEB 的投影约束可能使常规 IBRION=2 收敛困难，不能默认照搬一般结构优化设置。原生也可研究合适的 IBRION=3 配置，但它同样需要非零的有效步长。

路线 B：确认安装 VTST 后再做 CI-NEB

CI-NEB 让最高能 image 沿路径爬向鞍点，同时在其他方向下降。先预收敛普通 NEB，待路径形状与最高点位置合理后再开启 climbing image，通常比一开始就爬坡更稳。所用 VASP、VTST 版本与编译配置都必须记录；这里只给扩展片段，不假定你的计算中心已安装。

# TEMPLATE: VTST-ENABLED VASP ONLY; not a standalone INCAR
# Replace native ionic-optimizer lines with these (do not leave duplicates).
# Retain validated electronic settings and fixed-cell NEB directory structure.
# Confirm the executable was built with a compatible VTST code version.
IBRION = 3
POTIM = 0
IOPT = 7
ICHAIN = 0
IMAGES = 5
SPRING = -5
ISIF = 2
NSW = 250
EDIFFG = -0.03
LCLIMB = .FALSE.
# First preconverge the path. For the separate climbing-image continuation:
# change the existing LCLIMB line to .TRUE. after inspecting the path.

IOPT=7 指 VTST 的 FIRE 优化器。该片段先显式关闭 LCLIMB，确认路径后再在单独续算阶段开启。来源：VTST NEB 文档、优化器文档目录、开发者关于 POTIM=0 的说明。

能垒的验收

绘图横轴应是累计构型位移等定义清楚的路径坐标，不要把 image 编号称为时间。普通 NEB 最高离散点只是势垒近似；鞍点应进一步精化，并在相关自由度上验证一个连接初末态的虚模。发现两个虚模可能说明仍未到一阶鞍点；局部约束 Hessian 的结论也仅对检查过的子空间成立。有限温度自由能垒、隧穿和动力学前因子需要额外理论，不能由零温 NEB 曲线直接推成实验速率。

几何优化沿能量下降方向寻找驻点；第一性原理分子动力学把每个时刻的电子结构力代入核运动方程。它可以观察局域重排、键断裂和扩散，但一段短轨迹没有发生某事件，并不证明该事件不会发生。首先明确温度、体积、组成、初始相、模拟时长和要测量的物理量。

固定晶胞的 Langevin NVT 起步

下列例子以恰好两种 POTCAR 原子类型示范 300 K 固定体积轨迹。POTIM=1.0 在 MD 中单位是 fs；含 H 的高频体系、强碰撞或高温体系可能需要更小步长，必须实际检验。示例 2000 步只对应 2 ps，适合短程排错，不能视作通用的生产采样长度。

# TEMPLATE: fixed-cell Langevin NVT; EXAMPLE for exactly TWO species
# Supply relaxed POSCAR, tested KPOINTS/POTCAR and electronic settings.
# Replace ENCUT and adapt one LANGEVIN_GAMMA value per POTCAR species.
SYSTEM = nvt-langevin-test
GGA = PE
ENCUT = REPLACE_WITH_FORCE_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-6
NELM = 120
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.05
IBRION = 0
MDALGO = 3
ISIF = 2
ISYM = 0
POTIM = 1.0
NSW = 2000
TEBEG = 300
TEEND = 300
LANGEVIN_GAMMA = 10.0 10.0
LWAVE = .FALSE.
LCHARG = .FALSE.
# No SMASS: it does not set the Langevin thermostat strength.
# 2000 * 1 fs = 2 ps; this is a short test, not a production guarantee.

LANGEVIN_GAMMA 的顺序跟 POTCAR 类型顺序一致；相同元素若被分成不同类型，也要按类型分别列。SMASS 用于其他热浴或某些离子算法，不能拿它代替 Langevin 摩擦参数。MD 支持依赖 -Dtbdyn 编译选项，当前官方模板通常包含它，仍应检查自己的可执行文件。参见 MDALGO、摩擦系数与 SMASS。

热平衡与生产采样分开

1. 短测试：先跑少量步，检查是否原子重叠、温度异常、每一步电子循环耗尽 NELM。输入结构、力精度或时间步不合格时，增加总时长只会扩大问题。
2. 热平衡：在目标条件下观察温度、势能、体积（若可变）、配位与结构量是否脱离初态。温度围绕目标值波动是正常的，单步恰好等于 300 K 不是平衡标准。
3. 生产段：丢弃有理由选定的平衡段后再统计。分块平均、不同时间起点和独立初始速度用于估计相关性与误差，不能把每个相邻 MD 帧当独立样本。
4. 验证步长：用更小时间步和更严电子收敛比较关键量；必要时在合适的 NVE 条件下检查能量漂移。恒温器会交换能量，不能要求 NVT 物理总能像理想 NVE 那样严格不变。

记录 OUTCAR/OSZICAR、REPORT、XDATCAR 与终态 CONTCAR 的对应关系；重启时若希望连续轨迹，要保存速度及程序所需的重启状态，而不是只复制几何后重新随机分配速度。官方分子动力学流程、初始速度与温度。

扩散系数与“看起来在动”不同

计算均方位移前要跨周期边界解包轨迹，并按需要去除整体漂移；分子体系还应区分质心扩散与分子内部振动。只有进入近似线性的长时扩散区，且经历足够多迁移事件后，斜率才有统计意义。三维用 d=3，严格二维扩散用 d=2；不能给各向异性体系只报一个数字而不交代方向。

计算参数有两类：一类决定数值精度，例如 ENCUT、k 点与 EDIFF；另一类定义近似物理模型，例如交换关联泛函、U、色散和自旋轨道耦合。提高前一类并不能自动纠正后一类的系统偏差。做研究时应先写下“希望改进什么物理”，再决定方法。

色散：分子、层状材料与弱吸附

半局域泛函往往不能充分描述长程色散相互作用。层间距离、分子取向、弱吸附能等问题可比较一致选择的色散修正或非局域 vdW 泛函。D3(BJ) 是一种几何相关的色散修正，本教程以 PBE 配套 IVDW=12 示范；IVDW=11 则对应另一种阻尼形式，数值不可当成“精度等级”随意互换。

# TEMPLATE OVERLAY ONLY: PBE with D3(BJ) dispersion
# Add to ONE otherwise complete, validated PBE INCAR.
# Do not stack onto another dispersion-corrected functional uncritically.
GGA = PE
IVDW = 12
# Apply the same method to all comparable structures / references.

使用时要检查两件事：实际输出是否识别了预期色散方案，以及该泛函是否有支持的阻尼参数。不要为了贴近实验单独调一组经验系数却仍称为标准 PBE-D3(BJ)。换成别的泛函，或使用需要外部库的新方案，必须重新核对版本与参数支持。D3 官方方法与参数说明。

1. 对同一体系比较 PBE 与选定 vdW 方法的结构，例如层间距或分子—表面高度；随后看目标能量差的变化，而非只比较两次总能数值。
2. 若最终报告采用 vdW 方法，应在该方法下取得与问题对应的结构或明确声明“在 PBE 几何上的单点修正”。几何响应与能量修正不是同一件事。
3. 计算吸附或相变能时，各项必须使用同一方法；不能只给吸附态加色散而给气相/清洁表面用不同定义。
4. 不要把已经含非局域色散的泛函再无条件叠加 D3，否则可能重复计数；含 U、磁性或表面模型的适用性也要独立验证。

HSE06：轨道依赖的屏蔽混合泛函

HSE06 用一定比例的短程精确交换替换半局域交换。典型定义是 PBE 基础、AEXX=0.25 与 HFSCREEN=0.2 Å⁻¹；只写 LHFCALC=.TRUE. 并不等于已经选择 HSE06。混合泛函常用于检查半导体带隙、局域电子与某些缺陷，但“更贵”不等于对每一种能量、金属体系或实验条件都更准确。HFSCREEN。

# TEMPLATE: HSE06 static self-consistent calculation from scratch
# For a nonmagnetic semiconductor; regular k mesh REQUIRED.
# Supply POSCAR, KPOINTS and suitable licensed PAW-PBE POTCAR.
SYSTEM = hse06-static
GGA = PE
LHFCALC = .TRUE.
AEXX = 0.25
HFSCREEN = 0.2
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-7
NELM = 150
ALGO = Damped
TIME = 0.4
LREAL = .FALSE.
LASPH = .TRUE.
ISMEAR = 0
SIGMA = 0.05
ISTART = 0
ICHARG = 2
IBRION = -1
NSW = 0
LWAVE = .TRUE.
LCHARG = .TRUE.
# PBE WAVECAR restart is possible with compatible settings and ISTART=1;
# review ICHARG/restart settings then. NEVER set ICHARG=11 for this task.

模板显式启用 LASPH=.TRUE.，包含 PAW 球内密度梯度的非球形成分；这对含 d/f 电子、过渡金属氧化物等体系的混合泛函尤其重要。用于方法对照的 PBE 计算也应一致启用该设置，避免把设置差异混入泛函差异。LASPH 官方说明。

ALGO=Damped 与 TIME 控制电子轨道优化；TIME 不等于离子优化或 MD 的 POTIM。发生振荡时可测试更保守的 TIME。官方不建议为混合泛函沿用 ALGO=Fast；具体优化器和版本性能应以小任务测试。混合泛函的电子优化建议。

为什么不能用 PBE 的 CHGCAR + ICHARG=11 算 HSE 能带？

半局域哈密顿量可由密度确定，而混合泛函还需要规则 k 网格上占据轨道构造非局域交换。固定 PBE 电荷并只沿高对称线求解，会遗漏必要信息。因此 HSE 必须按轨道依赖方法自洽，重启应保留兼容的 WAVECAR；PBE 波函数可作初猜，却不是最终 HSE 波函数。

带结构可用“有权重的规则网格 + 零权重路径点”，或使用支持版本的 KPOINTS_OPT 让规则网格与路径分开定义。两者都不能把规则网格删除，也不能设置 ICHARG=11。路径点的收敛检查与内存成本也不同于普通 SCF，先阅读官方 HSE 能带流程，再使用电子结构章节的绘图原则。

对于有带隙体系的 HSE 能带，官方建议研究 HFRCUT=-1 的交换库仑截断以避免某些路径不连续。显式零权点的轨道不直接决定 SCF 总能判据，从已收敛 HSE 波函数重启时可能过早停止；应独立检查路径能带的迭代收敛，必要时验证增加 NELMIN 的效果。这个设置属于能带专门流程，不能用“总能已收敛”代替路径轨道检查。

介电函数描述电场扰动下材料的响应，通常是依赖频率和方向的复张量。频率很高时原子核来不及运动；更低频时晶格振动也能贡献极化。因此“介电常数”至少要说明电子贡献还是电子加离子贡献、静态还是动态、哪个晶向，以及是否包含局域场与激子效应。

独立粒子光学谱的实操顺序

1. 用明确的几何与泛函建立电子基态。光学积分用密集规则 k 网格，不是能带图用的高对称路径。
2. 增加空带并查看最高带覆盖的跃迁能量范围。NBANDS 包含占据和空态，程序默认的空态余量往往不是光学收敛值。逐级增加直到目标能区谱形稳定，不能只看总能是否变化。
3. 分别加密 k 点、增加频率采样、改变展宽并比较。光谱更平滑可以只是展宽更大，不代表微观计算更准确。
4. 在 OUTCAR 或支持的输出中提取实部与虚部张量，标明分量、能量单位和展宽。每个曲线都要能追溯到同一套输入。

# TEMPLATE: independent-particle PBE optics of a nonmagnetic 3D insulator
# Supply optimized POSCAR, DENSE regular KPOINTS, POTCAR; replace two fields.
# This is a self-consistent static calculation, not a band-line calculation.
SYSTEM = pbe-independent-particle-optics
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-8
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.01
IBRION = -1
NSW = 0
LOPTICS = .TRUE.
NBANDS = REPLACE_WITH_CONVERGED_NUMBER_OF_BANDS
NEDOS = 2000
CSHIFT = 0.10
LWAVE = .TRUE.
LCHARG = .TRUE.

LOPTICS 在电子基态后计算频率响应，标准独立粒子近似忽略相应局域场变化；不能把谱峰直接称为实验激子峰。它可与混合泛函结合，但仍需采用混合泛函自洽与适用优化器。电子带隙、光学吸收起始与含声子/激子的实验阈值可能不同。LOPTICS 的近似与适用条件。

CSHIFT=0.10 的单位为 eV，用于 Lorentz 展宽；同时使用高斯占据时也有 SIGMA 影响，所以模板让 SIGMA 明显小于 CSHIFT，仍须检查谱尾与频率范围。NEDOS=2000 只是输出频率采样示例，不是 k 点数。参见 CSHIFT、NEDOS与 NBANDS。

电子静态响应与 Born 电荷

# TEMPLATE: electronic (ion-clamped) static dielectric tensor / Born charges
# For a nonmetallic system with supported semilocal DFT, NOT HSE.
# Supply optimized POSCAR, converged KPOINTS/POTCAR, replace ENCUT.
SYSTEM = pbe-dielectric-dfpt
GGA = PE
ENCUT = REPLACE_WITH_CONVERGED_ENCUT
PREC = Accurate
EDIFF = 1E-8
ALGO = Normal
LREAL = .FALSE.
ISMEAR = 0
SIGMA = 0.01
IBRION = -1
NSW = 0
LEPSILON = .TRUE.
LRPA = .FALSE.
LWAVE = .FALSE.
LCHARG = .FALSE.
# Not the full ionic+electronic static permittivity of a relaxed crystal.
# Do not append LOPTICS or excessive NBANDS to this separate DFPT task.

Born 有效电荷是原子位移与宏观极化变化之间的响应张量，不等于 Bader 净电荷，也不必等于形式氧化态。LEPSILON 的 DFPT 路线不需要像 LOPTICS 那样大量增加空带，LRPA=.FALSE. 对应保留适当交换关联响应的设置；两条任务分开做更清楚。LEPSILON。

需要完整静态介电响应时，进一步结合适用的振动与离子响应计算，并把电子 ε∞ 和离子贡献分别报告。HSE 等轨道依赖方法可研究有限电场路线，但场强需要做线性区测试，绝缘性也必须保持。官方当前说明 LEPSILON 不适用于混合泛函；参考 LCALCEPS。

进阶计算最容易出现的失败不是程序崩溃，而是完成了一套与科学问题不匹配的计算。每个任务开始前写一句目标，例如“比较相同覆盖度下两个吸附位点的电子吸附能”，再据此决定需要多严格的能量差、力、尺寸和温度采样标准。

从症状回到最小可诊断问题

排错时保存原输入与日志，在新目录一次只改变有理由的一组设置。可以用较便宜设置定位错误，但最终物理结果必须恢复并验证所需精度。不要往 INCAR 中同时堆入十几个来源不明的标签，使问题无法归因。官方电子收敛排查。

先做小规模计时，再估算总成本

• 声子：总成本取决于不等价位移数乘以每个超胞静态受力成本；每个位移可独立提交，但要避免文件排序错误。
• NEB：中间结构数量、每幅最困难的电子循环与离子步数都会影响等待时间。并行 image 不会让总计算资源消耗消失。
• AIMD：用试跑测每步平均成本，再乘生产步数，并为热平衡、重复轨迹与失败重启预留资源。
• HSE/光学：更多轨道、k 点与交换计算可能显著增加内存；把核数加倍不保证耗时减半，甚至可能因通信更慢。

使用本机型、本程序版本和近似相同体系测试并行。NCORE、KPAR、image 分组与 GPU/OpenMP 限制互有关联，不能把某个 CPU 节点上的最佳参数当通用模板。先核验 dry-run 或 ALGO=None 支持，再检查实际 bands、k 点、内存和短程计时。官方并行优化流程。

人工智能能把哪里做得更快？

机器学习势可从 DFT 的能量、力和应力学习势能面，在经过验证的构型范围内扩大动力学时间与空间尺度；主动学习则优先补算模型不熟悉的局部环境。VASP 本身提供基于在线学习的 MLFF 路线，具体标签随版本与模式而定，入门应先学会评估数据质量。官方 MLFF 总览。

合理的验证集应覆盖独立温度、结构扰动、缺陷或反应环境，而不是把同一条轨迹相邻帧随机分到训练和测试。低平均力误差不保证迁移能垒、罕见事件和新相也正确；关键候选结构和鞍点仍应返回相同基准 DFT 复核。模型继承训练泛函的偏差，无法靠数据量自动变成实验真值。

大语言模型适合解释输出、帮助生成检查脚本与整理文档。它也可能编造不存在的参数或把 VTST 当原生功能，所以将自动建议转成计算前，应核对官方定义、单位、兼容性，并检查真正生效的输出。研究报告分别标明“AI 预测”“DFT 复核”“有限温度模拟”与“实验验证”。

一份合格的任务记录

1. 问题与模型：结构来源、晶胞、组成、电荷、自旋、表面终止、固定层和参考态。
2. 软件与输入：VASP/外部工具版本、构建或扩展信息、POTCAR 标识及校验信息、泛函、ENCUT、k 点、展宽、优化与响应参数。
3. 证据：电子/离子收敛、超胞与数值精度测试、所用输出文件、后处理脚本和单位转换。
4. 结论范围：哪些是确定的数值趋势，哪些是当前模型的预测，哪些仍需更大模型、更多采样或实验检验。

不要只保留截图和最终表格。保存原始输入、软件输出、结构与分析流程；公开可复现材料时不分发受许可限制的 POTCAR 文件。回到课程地图选择练习，或在参数与检查清单逐项核对。

先读文件说明与对应章节，再补齐模型、POTCAR 与前置计算。后缀为 .template 的文件包含待改参数，不能原样提交。

正文为独立整理的学习说明。标签含义和版本差异以与你的程序版本对应的官方文档为准。

1. 官方 POSCAR 定义
2. IDIPOL
3. LDIPOL
4. DIPOL
5. LVHAR
6. Freysoldt 等的缺陷综述
7. NELECT 官方说明
8. 各向异性与弛豫相关的电荷校正研究
9. 官方带电体系讨论
10. 有限位移官方说明
11. NFREE
12. 命令迁移说明
13. VASP 接口
14. 后处理命令
15. Phonopy 输入文件与 BORN
16. 官方虚频排查与软模操作
17. IBRION 官方说明
18. ISIF 的应力定义
19. 弹性矩阵列序讨论
20. Mouhat 与 Coudert 的稳定性判据
21. 原生 NEB
22. IMAGES
23. SPRING
24. 官方排错说明
25. VTST NEB 文档
26. 优化器文档目录
27. 开发者关于 POTIM=0 的说明
28. MDALGO
29. 摩擦系数
30. SMASS
31. 官方分子动力学流程
32. 初始速度与温度
33. D3 官方方法与参数说明
34. HFSCREEN
35. 混合泛函的电子优化建议
36. 官方 HSE 能带流程
37. LOPTICS 的近似与适用条件
38. CSHIFT
39. NEDOS
40. NBANDS
41. LEPSILON
42. LCALCEPS
43. 官方电子收敛排查
44. 官方并行优化流程
45. 官方 MLFF 总览
46. Phonopy 安装与 seekpath 依赖
47. LASPH：PAW 非球形成分