VASP如何计算应变效应？

引言

应变效应听起来很直观：把晶格拉长或压短，再算吸附能变化。但真正做起来，问题会马上出现：到底固定面内晶格还是让原子完全弛豫？吸附前后是否用同一个应变基准？应变改变的是电子结构，还是只是改变了键长？

应变计算最怕变量没锁住。拉伸和压缩不是简单改 POSCAR 的晶格常数，而是在定义一个新的受力表面状态。吸附能怎么变，要放在这个表面状态里理解。

实际出问题的地方往往很细：有人只改吸附态的晶格，却忘了对应应变下的干净 slab；有人把 z 方向真空也按比例缩放，导致真空层跟着变；还有人让 ISIF 把晶格释放掉，最后得到的已经不是指定应变。应变稿要写得扎实，必须把这些边界条件交代清楚。

应变到底改变了什么？

面内晶格先变，局域键长跟着变

对 slab 来说，常见做法是改变面内 a、b 晶格，保持真空方向不参与应变。拉伸会让表面原子间距变大，压缩会让表面原子靠近；吸附物和表面之间的最佳键长也会跟着重新平衡。

这意味着吸附能变化不只来自电子能级移动，也来自几何匹配。拉伸后吸附变强或变弱都可能发生，方向不能预设。

举个很常见的例子：*H 在金属顶位吸附时，拉伸可能让表面金属原子更容易向外松弛，M-H 键变短，ΔG_*H 接近 0；但对 *O 或 *OH 这类更强吸附物，过度拉伸也可能削弱多中心配位，使吸附反而变差。同样是拉伸，吸附物不同，响应就不一定同向。

应变会改变表面态和 d 态分布

拉伸通常会减小轨道重叠，使 d 带变窄，d 带中心可能上移；压缩则可能增强重叠，使 d 带展宽并下移。但这只是常见趋势，不是铁律。合金、氧化物和二维材料还会受到局域配位和电荷重排影响。

因此，应变效应最好同时报告几何和电子结构。只说“应变调节 d 带中心”不够，还要说明关键吸附键长和局域态怎样变。

氧化物表面还要多看一个因素：应变可能改变金属–氧八面体畸变、氧空位形成倾向和表面羟基稳定性。此时吸附能变化未必来自金属 d 态整体移动，而可能来自局部 M-O-M 键角、M-O 键长和表面 O 2p 态的共同变化。

图1 RuO₂中应变异质性与酸性 OER 性能关联示例。 DOI：10.1038/s41467-025-58570-3

VASP 模型里哪些量要固定？

吸附前后必须处在同一个应变状态

算某个应变下的吸附能时，干净 slab 和吸附 slab 要使用同一个面内晶格。不能用未应变 slab 的能量去减应变吸附态的能量，否则吸附能里会混入 slab 变形能。

E_ads(ε) = E_slab+X(ε) − E_slab(ε) − E_X

这里的 ε 表示同一个拉伸或压缩状态，slab+X 和 slab 的面内晶格必须一致。

如果想把变形能也讲清楚，可以另外保存 E_slab(ε) − E_slab(0)。这不是吸附能的一部分，但它能告诉你某个应变状态本身是否代价太高。一个吸附能很漂亮、但 slab 变形能很大的模型，未必有现实意义。

原子可以弛豫，但晶格边界要按问题固定

如果研究外加应变对表面反应的影响，面内晶格通常固定在指定 ε，原子坐标允许弛豫。这样模拟的是受基底或外力约束的表面。若让晶格完全释放，应变就会在弛豫中消失，最后不再是应变模型。

真空方向可以保持足够厚，不要跟着面内应变乱缩放。对于二维材料，也要注意泊松效应：如果你允许某个方向自由响应，就要在方法里写清楚，否则拉伸/压缩的物理含义会变模糊。

VASP 设置上，常见做法是先手动改 POSCAR 的面内晶格，再让离子坐标弛豫，而不是让计算自己优化晶格。若只想固定底层，可以继续用 Selective dynamics 固定深层原子，但这和固定面内晶格不是同一件事：前者限制原子位移，后者定义应变边界。

单轴应变和双轴应变不要混着说

双轴应变同时改变 a、b 两个面内方向，适合模拟均匀外延或整体拉伸；单轴应变只改一个方向，另一个方向可能固定，也可能按泊松响应调整。两种应变会给出不同的位点距离和轨道重叠，不能把结果混在一条趋势里。

如果表面本身各向异性，比如台阶面、二维条带、正交晶胞或畸变氧化物，单轴应变方向尤其重要。沿金属–金属键拉伸和垂直于金属–金属键拉伸，可能对应完全不同的吸附位点变化。应变方向本身就是模型条件，不能只写一个百分比。

图2 应变调控 RuO₂结构和表面性质的表征示例。 DOI：10.1038/s41467-025-58570-3

吸附能变化怎样拆开看？

先分清变形能和吸附相互作用

应变下的 slab 本身能量会升高或降低，吸附物上去以后也会诱导局域弛豫。若只看 E_ads(ε)，很难知道变化来自表面预应变，还是来自吸附后结构进一步释放。

可以比较同一应变下的干净 slab、吸附态和关键键长，再看差分电荷或 PDOS。若键长变化主导，电子态证据不一定显著；若 PDOS 峰位和占据明显改变，电子效应才更有说服力。

更细一点，还可以看吸附前后表面原子的回弹幅度。若干净 slab 在 +4% 应变下已经很紧张，吸附物上去后局部原子明显内缩，说明吸附过程部分释放了预应力；这类吸附能变化就不能简单写成“电子结构优化”。

拉伸和压缩最好成对出现

只算 +2% 拉伸，很难判断趋势是不是偶然。更稳的序列是 -4%、-2%、0%、+2%、+4%，至少要有压缩、原始和拉伸三点。若吸附能随应变近似线性，说明应变响应比较平滑；若出现拐点，往往对应构型切换或表面重构。

应变效应不是单点结论，必须看一条随 ε 变化的趋势线。这一点比堆很多中间体更重要。

趋势线也不要只画 E_ads。最好同时给出关键键长或 ΔG_*H 这类反应指标。比如 E_ads 随拉伸单调增强，但 ΔG_*H 从正值跨过 0 后继续变负，这说明最优应变可能在中间，而不是越拉越好。

图3 不同应变环境下电催化活性变化示例。 DOI：10.1038/s41467-025-58570-3

什么时候应变结果不该继续用？

表面已经重构时，不能再叫小应变

若压缩后表面原子错位、吸附物跑到新位点，或拉伸后键断裂，这个模型已经不再是同一表面的弹性应变响应。它可以作为重构结果讨论，但不适合和 0% 表面直接画成一条应变趋势线。

判断标准可以看表面原子位移、关键键断裂、吸附位点是否改变。如果 +4% 和 -4% 的低能构型已经不是同一类，最好把趋势分段讨论。

还有一种不明显的重构，是吸附物没有跑位，但表面第一层和第二层相对滑移了。图上看吸附位点还在，实际上局部配位已经变了。遇到这种情况，最好补一张俯视/侧视结构，或者报告关键 M-M、M-O、M-X 距离，别只给一条能量曲线。

过大应变会制造不现实的活性

很多二维材料或薄膜理论上可以承受较大应变，但电催化工作条件下未必能稳定存在。若 +8% 应变给出很漂亮的 ΔG_*H，也要问这种应变是否有实验实现方式，比如外延基底、缺陷诱导、晶格失配或弯曲结构。

比较稳的做法是把小应变当成可调范围，把大应变当成机制探索。应变越大，越要证明结构仍然稳定，而不是只展示吸附能变优。

写法上可以把 0–4% 作为主要讨论区间，把更大应变放到补充分析里说明趋势极限。这样既不浪费计算结果，也不会让主结论依赖一个很难实现的结构状态。

图4 应变相关反应中间体和理论分析示例。 DOI：10.1038/s41467-025-58570-3

本章要点总结

要点：应变计算要固定面内晶格边界，吸附前后必须使用同一个应变状态。

要点：吸附能变化可能来自几何键长，也可能来自电子结构，不能只讲 d 带中心。

要点：应变效应要看压缩、原始、拉伸的趋势线，不建议只算一个应变点。

要点：若表面发生重构或大应变不可实现，就不能把结果当成普通弹性应变响应。