1. 项目概述：为什么截断能测试是VASP计算的第一道门槛

做第一性原理计算的朋友，尤其是刚接触VASP的，拿到一个体系准备开算时，第一个要确定的参数是什么？很多人会说是结构优化，但在我看来，在点击“运行”按钮之前，有一个更基础、更关键的参数必须首先敲定，那就是截断能（ENCUT）。你可以把它理解为计算精度的一把“标尺”，或者更形象地说，是构建我们理论模型时所用的“画笔”的粗细。画笔太粗（ENCUT太低），画出来的图像模糊不清，能量、力、电子结构等关键物理量都不准，后续所有计算都是建立在沙滩上的城堡；画笔太细（ENCUT太高），计算量会呈立方级增长，可能一个简单的单点能计算就要跑上好几天，严重拖慢科研进度。

所以，这个“VASP(1)_参数测试_截断能(ENCUT)测试”项目，本质上就是一次针对特定材料体系的“标尺校准”工作。它不是可有可无的预备步骤，而是确保我们所有后续计算结果可靠、可重复、可比较的基石。我见过太多因为ENCUT设置不当导致的计算结果异常：晶格常数偏差百分之几，形成能符号都反了，或者声子谱出现虚频却找不到结构上的原因，回头一查，问题往往就出在这个最基础的参数上。因此，无论你是计算半导体、金属、二维材料还是表面催化，在开展任何严肃的计算工作前，花上几个小时系统地做一次ENCUT测试，绝对是性价比最高的时间投资。

2. 核心原理：截断能到底是什么，以及它如何主宰计算成本与精度

要理解为什么必须测试ENCUT，我们得先搞懂它在VASP的平面波赝势方法中扮演的角色。VASP求解Kohn-Sham方程时，将电子波函数用平面波基组展开。平面波有个很好的性质：它们构成一组完备的正交基。但“完备”意味着需要无穷多个平面波，这显然无法计算。因此，我们必须引入一个截断，只使用动能低于某个特定值的平面波。这个特定的动能值，就是截断能ENCUT，其单位是电子伏特（eV）。

数学上，平面波的动能是 (ħ²|k+G|²)/(2m)，其中k是倒空间波矢，G是倒格矢。ENCUT定义的截断条件就是：只保留所有满足 (ħ²|k+G|²)/(2m) < ENCUT 的平面波基。所以，ENCUT直接决定了我们基组的规模，也就是计算精度上限。

这里有一个至关重要的衍生参数：平面波基组的数量。它大致与ENCUT的3/2次方成正比。这意味着，如果你把ENCUT提高一倍，平面波数量大约会增加2.8倍。而VASP中大多数计算模块（如电子自洽、力/应力计算）的计算量，又与基组数量的平方甚至三次方相关。因此，ENCUT增加10%，计算时间可能增加30%以上；增加50%，计算时间可能翻好几倍。这就是精度与计算成本之间最直接的权衡。

那么，是不是ENCUT越高越好呢？理论上，在完备基组极限下，结果会收敛到真实值。但实践中，我们使用的是赝势（Pseudopotential）。赝势在生成时，本身就在某个动能截断下进行了测试和优化，这个值通常记录在赝势文件（POTCAR）中，称为“赝势的推荐截断能（ENMAX）”。每个元素的POTCAR里都有这个值。VASP的默认行为是：取你体系中所有元素ENMAX的最大值，作为默认的ENCUT。这是一个安全的起点，但通常不是最优解。

为什么？因为对于很多体系，特别是包含“硬”元素（如O、F、第一行过渡金属）和“软”元素（如Na、K、Au）的化合物，最硬元素决定的ENMAX可能会远高于其他元素实际需要的精度。对所有元素都使用这个统一的、最高的ENCUT，会造成对较软元素的“过描述”，浪费大量计算资源。因此，我们测试的目的，就是找到在保证整体计算精度的前提下，那个计算成本最低的ENCUT值，也就是能量的收敛点。

3. 测试方案设计：如何科学地扫描与判断收敛

知道了为什么测，接下来就是怎么测。一个严谨的ENCUT测试方案，需要明确测试对象、扫描范围、评估标准和输出结果。

3.1 测试对象与初始结构准备

测试不能在空中楼阁中进行，必须针对你实际要研究的具体体系。通常，我们使用该体系经过初步松弛（可以用较低精度或默认ENCUT快速跑一下）后的稳定结构。用这个稳定结构做测试，结果才对你的后续计算有指导意义。如果你研究的是缺陷，那就用包含缺陷的超胞；如果是表面，就用表面slab模型。文件准备上，你需要标准的四个输入文件：INCAR, POSCAR, POTCAR, KPOINTS。其中POSCAR是你的测试结构，KPOINTS可以设置一个较密的网格（确保k点采样误差不干扰ENCUT测试结果），POTCAR用你计划一直使用的赝势。

关键的INCAR设置如下：

SYSTEM = ENCUT_Test # 任务名，便于识别 ISTART = 0; ICHARG = 2 # 从头开始计算 PREC = Accurate # 高精度模式，影响一些算法细节 ISMEAR = 0; SIGMA = 0.05 # 对于半导体/绝缘体，用Gaussian展宽，设置小展宽 EDIFF = 1E-6 # 电子步收敛标准，设得紧一些，让能量值更精确 NSW = 0 # 不做离子弛豫，只做单点能计算 IBRION = -1 # 同上，固定离子位置 LREAL = .FALSE. # 在倒空间求投影算符，精度更高（虽然慢点，但测试时精度优先） # 注意：这里不设置ENCUT，我们通过脚本或多次提交来改变它

注意：LREAL = .FALSE.对于测试很重要。在实空间投影(LREAL = .TRUE.或Auto)虽然快，但会引入额外的、与ENCUT相关的误差，干扰我们对纯平面波基组收敛性的判断。测试阶段务必关掉它。

3.2 ENCUT扫描范围的确定

这是测试的第一步智慧。盲目地从100 eV扫到1000 eV既无必要又浪费资源。

1. 查找ENMAX：首先用grep ENMAX POTCAR命令查看你所用赝势的推荐值。假设你的体系是Si和O，可能会看到Si: ENMAX = 245.0 eV，O: ENMAX = 400.0 eV。那么默认ENCUT就是400 eV。
2. 设定扫描区间：通常，从最大ENMAX的0.7倍到1.5倍开始扫描是合理的范围。以上述为例，就是从 280 eV (4000.7) 扫到 600 eV (4001.5)。如果你的计算资源紧张，或者体系较大，可以聚焦在0.8倍到1.2倍这个更窄的区间（320 eV 到 480 eV）。
3. 设置步长：初始扫描可以用较大的步长，比如20 eV或40 eV，以快速定位收敛区域。在能量变化剧烈的低ENCUT区间，或者接近收敛的高ENCUT区间，可以改用更小的步长（如10 eV）进行精细扫描。

3.3 核心评估标准：能量收敛判据

我们主要观察体系的总自由能（在OUTCAR中查找free energy TOTEN）。通常，我们关注绝对能量随ENCUT的变化。收敛的判断标准没有国际硬性规定，但社区内有一些经验准则：

• 宽松标准：总能量变化 < 1 meV/atom。适用于对能量精度要求不极高的初步筛选或大体系。
• 严格标准：总能量变化 < 0.1 meV/atom。适用于计算形成能、吸附能、相变势垒等对能量差极其敏感的性质。
• 我的常用实践：我会要求连续三个递增的ENCUT算出的总能量，其差值都小于我设定的阈值（例如0.5 meV/atom）。这比只看最后两点更稳健。

除了总能量，对于涉及力、应力、晶格优化的计算，还需要关注力（forces）和应力（stress）的收敛情况。有时能量收敛了，但力的收敛要慢一些。如果你后续要做几何优化，最好在ENCUT测试时就检查一下主要原子上的力分量随ENCUT的变化是否也趋于平稳。

3.4 自动化测试与数据处理

手动修改INCAR并提交几十个作业是不可接受的。必须借助脚本。一个简单的Bash脚本示例如下：

#!/bin/bash # 文件名：run_encut_test.sh # 定义扫描范围 for ENCUT in {300,320,340,360,380,400,420,440,460,480,500}; do mkdir ENCUT_${ENCUT} cd ENCUT_${ENCUT} # 拷贝输入文件 cp ../INCAR ../POSCAR ../POTCAR ../KPOINTS . # 创建新的INCAR，在原有基础上添加或替换ENCUT行 sed "/ENCUT/d" ../INCAR > temp_incar echo "ENCUT = $ENCUT" >> temp_incar mv temp_incar INCAR # 提交作业（根据你的作业管理系统修改qsub/sbatch命令） sbatch ../vasp_job.slurm # 或者 qsub ../vasp_job.pbs cd .. done

计算完成后，再用一个脚本从各个目录的OUTCAR中提取总能量和每个原子的能量。用Python的pandas或简单的gnuplot都能轻松绘图。关键是生成“能量 vs ENCUT”和“能量差（相对于最高ENCUT的能量）vs ENCUT”的曲线图。后者能更直观地看到收敛情况。

4. 实操过程详解：从文件准备到结果分析

让我们以一个具体的例子走一遍流程：计算体相硅（Si）的ENCUT收敛性。我们使用PAW-PBE赝势。

4.1 初始设置与参数确认

首先，准备一个松弛后的硅晶胞（POSCAR），晶格常数约为5.43 Å。使用grep ENMAX POTCAR发现Si的ENMAX是245.0 eV。因此，我们计划从170 eV（~0.7245）扫到370 eV（~1.5245），初始步长用20 eV。

KPOINTS设置一个较密的网格，例如9x9x9的Monkhorst-Pack网格，以确保k点误差远小于ENCUT变化引起的误差。

INCAR文件如下：

SYSTEM = Si_ENCUT_Test ISTART = 0; ICHARG = 2 PREC = Accurate ISMEAR = 0; SIGMA = 0.05 EDIFF = 1E-6 NSW = 0; IBRION = -1 LREAL = .FALSE. # ENCUT will be set by script

4.2 执行扫描与数据提取

使用上述脚本提交从170, 190, 210, ..., 370 eV共11个计算任务。所有任务完成后，运行一个数据提取脚本。这里给出一个简单的Python示例：

import os import matplotlib.pyplot as plt encut_list = [] energy_list = [] for encut in range(170, 371, 20): # 注意range的右边界是开区间 dir_name = f'ENCUT_{encut}' outcar_path = os.path.join(dir_name, 'OUTCAR') if os.path.exists(outcar_path): with open(outcar_path, 'r') as f: lines = f.readlines() for line in lines: if 'free energy TOTEN' in line: # 提取能量，单位通常是eV energy = float(line.split()[-2]) encut_list.append(encut) energy_list.append(energy) break # 按ENCUT排序 encut_list, energy_list = zip(*sorted(zip(encut_list, energy_list))) # 计算相对于最高ENCUT的能量差（单位：meV/atom） # 假设是单原子晶胞，实际应根据原子数归一化 energy_max = max(energy_list) energy_diff_meV = [(e - energy_max) * 1000 for e in energy_list] # 转换为meV print("ENCUT(eV), Total Energy(eV), Diff(meV)") for enc, en, diff in zip(encut_list, energy_list, energy_diff_meV): print(f"{enc}, {en:.6f}, {diff:.3f}")

假设我们得到的数据如下表所示：

4.3 结果分析与收敛点判定

将上表数据绘图。X轴是ENCUT，Y轴是能量差（meV）。你会看到一条曲线从高处迅速下降，然后在某个点之后变得非常平坦。

如何选择最终的ENCUT？

1. 观察收敛平台：从上表看，从250 eV开始，能量差已经小于0.01 meV，进入了平台区。
2. 权衡精度与效率：ENMAX是245 eV，我们的收敛点在250 eV。这意味着使用默认值（245 eV）已经非常接近收敛平台，但严格来说，250 eV更保险。考虑到计算成本，245 eV和250 eV的计算时间几乎无差别。这里，我会选择250 eV作为我后续所有计算的ENCUT。它比ENMAX略高一点，确保了充分的收敛余地，又不会带来不必要的计算负担。
3. 考虑安全边际：在确定收敛点后，我通常会再加10-20%的安全余量。对于收敛点250 eV，加10%就是275 eV。有时为了确保极端情况（如高压、受力很大）下的稳定性，我会选择这个带余量的值。但在这个Si的例子中，从250 eV到270 eV能量已无变化，所以250 eV本身已足够安全。

实操心得：不要只看最后一个点是否满足阈值。要观察曲线的变化趋势。如果曲线在某个值之后变得完全平坦（如上例中270 eV后），那这个值就是可靠的收敛点。如果曲线还在以非常缓慢的斜率下降，你可能需要测试到更高的ENCUT，或者评估这个缓慢下降对你的目标物理量（如能量差）影响有多大。

5. 高级技巧与疑难问题排查

掌握了基本流程后，一些进阶技巧和常见坑点能让你事半功倍。

5.1 多元素体系与PRECFOCK参数

对于包含多种元素的体系，尤其是那些ENMAX相差很大的（比如NaCl：Na的ENMAX可能~150 eV，Cl的~400 eV），你需要格外小心。统一使用400 eV会对Na的计算造成浪费。VASP提供了一个解决方案：不同元素使用不同的截断能，通过PRECFOCK关键字控制。但这属于更高级的优化，在初始测试阶段，我建议仍然使用统一的、基于最高ENMAX的扫描。确定一个保守的统一值后，如果计算量仍然巨大，再考虑使用PRECFOCK进行“双网格”技术优化，这需要额外的测试。

5.2 测试结果与后续计算参数的关联

你测试得到的“最优ENCUT”是在一组特定参数下得到的：特定的KPOINTS、PREC=Accurate、LREAL=.FALSE.等。当你改变这些参数时，收敛的ENCUT可能会变！

• KPOINTS：如果你的测试用了很密的k点，但后续实际计算用了较疏的k点，理论上所需的ENCUT可能会略有不同，但通常影响很小，可以忽略。为保险起见，测试时的k点密度应不低于你主要计算任务所用的密度。
• PREC：PREC标志控制着许多算法内部的截断因子。PREC=Normal时，实际使用的平面波截断是ENCUT * ENCUTFACTOR（默认~0.75）。所以，如果你测试时用PREC=Accurate（对应ENCUTFACTOR=1.0），但后续计算改用PREC=Normal，那么你实际使用的精度是低于测试精度的。强烈建议：测试用什么PREC，正式计算就用什么PREC。通常，对于发表级计算，直接使用PREC=Accurate是省心的选择。
• LREAL：测试时务必用.FALSE.。正式计算中，为了速度可以开启LREAL=Auto，但要知道这会引入微小误差。对于精度要求极高的计算（如弹性常数、声子），保持.FALSE.更稳妥。

5.3 常见问题与排查清单

1. 能量不收敛，一直随ENCUT增加而下降：

   • 可能原因1：扫描上限不够高。特别是对于含有O、F、N等“硬”元素或过渡金属的体系，收敛可能需要较高的ENCUT。尝试扫到1.8倍甚至2倍的最大ENMAX。
   • 可能原因2：赝势问题。某些赝势（特别是早期或非标准的）可能在推荐ENMAX附近收敛性不佳。尝试换一种赝势（如从USPP换到PAW，或换一个来源）。
   • 检查：确认LREAL=.FALSE.，并且PREC=Accurate。

2. 能量曲线出现非单调的“跳动”：

   • 可能原因：这通常不是ENCUT本身的问题，而是电子自洽过程在某个ENCUT下陷入了局部极小或收敛困难。可以检查该ENCUT任务目录下的OUTCAR，看电子自洽迭代是否正常收敛（没有达到EDIFF就跳出）。可以尝试对该点使用更小的EDIFF（如1E-7）或更换ALGO（如ALGO=All）重新计算。

3. 测试时计算很快，但正式计算极慢：

   • 检查差异：对比测试和正式计算的INCAR。正式计算是否开启了离子弛豫（NSW>0）、分子动力学、或使用了更密的k点？这些都会极大增加计算量。ENCUT测试只是确定了基组大小，其他参数会叠加影响时间。

4. 如何为超胞或表面模型选择ENCUT？

   • 原则：与原胞相同。截断能是倒空间的动能截断，与实空间晶胞大小无关。你为原胞测试得到的ENCUT，直接用于任何放大后的超胞或表面模型都是适用的。这是平面波基组的一个巨大优势。

6. 测试报告的记录与传承

最后，但同样重要的是，养成记录的好习惯。为每个重要的材料体系建立一个简单的文本记录，内容应包括：

• 测试日期和VASP版本。
• 使用的赝势类型和ENMAX值。
• 测试的结构（POSCAR简要描述）。
• KPOINTS设置。
• 扫描的ENCUT范围和步长。
• 最终选择的ENCUT值及理由（附上能量收敛曲线图）。
• 任何观察到的异常及处理方式。

这份记录不仅是你工作的备份，未来当你或你的合作者需要复现或扩展这项工作时，它能节省大量重新摸索的时间。ENCUT测试看似基础繁琐，但它是构建可靠计算工作的第一步，把这步走扎实了，后面的路才能走得稳、走得快。我个人的习惯是，每开始一个全新的材料体系，无论多简单，都会强制自己先完成这个“标尺校准”流程，这几乎避免了我所有因基础参数不当导致的返工。