给一个材料加电场，最简单的想象是在两侧放两块带电板，像平行板电容器那样。VASP确实也是这么做的——通过EFIELD参数施加一个锯齿形势（sawtooth potential），在真空层中建立匀强电场。但锯齿形势的物理后果不那么简单：它会同时产生一个宏观极化、改变功函数、并可能让电子在真空层中隧穿。
这些后果不会出现在OUTCAR的任何一句报错里——VASP对EFIELD的容错性很高，算错了也不吭声。所以VASP计算电场的难点不在”怎么设置参数”，而在”怎么确认结果是物理的而不是数值的”。
INCAR中的电场设置只需要一行：
EFIELD = 0.1 # 单位：eV/Å（VASP默认单位）
但这个0.1 eV/Å到底是什么意思？
方向。 EFIELD的正方向是沿着晶格的第三个格矢方向（即POSCAR中第三个晶格矢量的方向）。这是一个很容易搞错的地方——如果你的slab模型把真空层设在了z方向（第三个格矢），EFIELD就是沿着z方向的。如果真空层在x方向（第一格矢），那EFIELD的正方向仍然是第三个格矢，和真空层垂直的就不是你设的电场了。检查方法：看OUTCAR里”external force on ions”的输出，如果所有原子受到的力都是沿着z方向的，那你的对准是没问题的。
强度。 0.1 eV/Å是一个很大的电场——换算成国际单位是约10 MV/cm。作为参照：Si-SiO₂界面在MOSFET工作电压下的电场约5-8 MV/cm（介电击穿的临界区），液体水在STM针尖下的电场约1 V/Å（即0.1 eV/Å），二维材料门电压调控的典型电场约0.01-0.1 eV/Å。如果你的EFIELD超过了材料的介电击穿强度，算出来的电子结构对应的是”材料在被击穿之前那一瞬间”的状态——结果本身可以发表，但必须在文章里说明这个前提。
锯齿形势的物理。 VASP在真空层中施加的势能是线性变化的——从一侧的V(z₀)=0到另一侧的V(z₁)=E·|z₁-z₀|。在slab两侧的表面上，电子感受到一个势能差ΔV = E × L_vac，其中L_vac是真空层的厚度。这意味着你用同样的EFIELD=0.1 eV/Å，如果真空层从15 Å变成了30 Å，slab两侧的势能差就从1.5 eV变成了3.0 eV——两倍的差异。所以报告中EFIELD参数时，必须同时报告真空层厚度，否则不可复现。
加电场后最直接的物理效应是能带的偏移和倾斜。
表面功函数的改变。 外加电场在slab的两侧产生不同的静电势。一侧表面静电势被抬高，另一侧被压低——功函数在两侧出现差异。在OUTCAR中搜索”vacuum level”可以读出两侧的真空能级，它们的差值应该约等于EFIELD × L_vac（允许有小于10%的偏差，来自表面偶极矩的屏蔽效应）。如果差值偏差超过20%，检查你的slab是否够厚——薄slab两侧的电荷重新分布可能部分屏蔽外电场，导致实际感受到的电场比外加值小。
能带的刚性偏移与倾斜。 外电场改变了静电势分布，导致能带在k空间和实空间都发生移动。对于二维材料（MoS₂、黑磷等），一个常见的做法是沿着面外方向（c轴）施加电场，观察带隙的变化。MoS₂单层在EFIELD=0.1 eV/Å下的带隙从约1.7 eV降到约1.5 eV——这个约0.2 eV的缩减来自上下两个硫原子层分别被不同的静电势”拉扯”，导致价带顶和导带底的相对位置改变。
电子局域函数（ELF）的变化。 强电场可以诱导电荷重新分布——在EFIELD>0.2 eV/Å时，某些化学键的ELF分布会出现明显的极化特征。这对于理解电场辅助催化的微观机制很重要，但这种级别的电场在实验上通常需要STM针尖或电化学双电层才能实现，自由空间中是达不到的。
铁电材料的自发极化是VASP计算电场的一个重要应用场景。极化P是一个多值量——它依赖于晶格中选择的参考原点——但极化的变化量ΔP（铁电相到顺电相的极化差）是物理上良定的。
计算ΔP的标准路线是Berry相位方法：
第一步，建两个态的结构。 铁电态（极性结构，原子偏离中心对称位置）和参考顺电态（中心对称结构）。中心对称结构不一定是实验上存在的相——它只是一个理论参考，确保参考极化是零（或整数倍的极化量子）。
第二步，逐一收敛。 两个态分别做高精度SCF计算：EDIFF ≤ 1E-7，PREC = Accurate。Berry相位计算对电荷密度的精度非常敏感——如果EDIFF不够紧，极化值的数值噪声可以达到0.1-0.5 μC/cm²，足以掩盖物理信号。
第三步，LCALCPOL = .TRUE.。 在INCAR中设置此标签，VASP会在SCF收敛后自动做Berry相位计算，输出三个方向的极化矢量。注意：LCALCPOL可以和EFIELD同时使用吗？可以——VASP在加电场的SCF收敛后也可以用Berry相位算极化，但此时极化包含外电场的贡献（诱导极化+自发极化），需要和EFIELD=0的结果对比才能分离出自发极化。
第四步，极化量子修正。 铁电材料的极化变化量可能出现”极化量子跳跃”——ΔP可能差了n×eR/Ω（其中R是格矢，Ω是原胞体积）。这是因为极化作为Berry相位的定义只模2π是确定的。如果两个态之间的极化差大于极化量子的一半，需要手动加减极化量子来修正。这个细节在第一次算铁电极化时往往会忽略——导致极化值和文献差了一个常数。
电催化反应发生在电极-电解液界面的双电层中，该区域的电场强度可以达到1 V/Å量级。VASP计算电场在电催化中的一个直接应用是：在slab上施加垂直于表面的电场，模拟不同电极电位下的反应中间体稳定性。
但这个方法有一个适用范围：它只模拟了静电效应——电场直接改变吸附态的能量。真正的电化学电位还包括电荷转移（电子从电极流入或流出）、溶剂化效应、以及阳离子在双电层中的具体分布。EFIELD参数只覆盖了静电这一项——对于d带中心对电场不敏感的反应（比如某些金属上的HER），纯静电近似是可接受的；对于d带中心随外加电位明显移动的反应（比如Pt上的ORR），必须配合其他方法（如DFT+溶剂化模型或显式水层）。
另外一种更精细的做法是用VASPsol（隐式溶剂模型）+ EFIELD的组合：VASPsol处理溶剂化和离子屏蔽的背景电荷，EFIELD施加特定的电极电位差。两者的结合比纯EFIELD更接近真实的电化学界面——但参数组合也更难调。没有独立验证手段的话，建议在同一种参数组合下只做趋势分析（”A催化剂的过电位比B催化剂低多少”），而不直接预测绝对电位值。
科研学术网（https://www.keyanxueshu.com）计算提供实战参考。
VASP计算电场是一个门槛不高但陷阱不少的功能。一个参数错误的EFIELD设置不会报错——它只会静悄悄地生成一套看起来合理、实际上偏离了物理的结果。唯一的防线是：做完计算后停下来，问自己这个电场在物理上意味着什么，然后从电子结构的变化里找到自洽的物理解释。
