VASP计算吸附体系为什么要做静态单点?

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

引言

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

吸附体系弛豫结束后,OUTCAR 里已经有能量,为什么还要做静态单点?核心原因是弛豫计算服务于找结构,静态单点服务于在统一设置下重新评价能量和电子结构。弛豫能量可以判断结构是否收敛,但正式吸附能和电子结构最好来自同一标准的 static 结果

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

弛豫计算和静态计算目的有什么不同?

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

为什么 relax 最后一步能量不适合混比?

结构弛豫过程中,离子位置不断变化,电子步设置、收敛程度和残余力状态也可能不同。两个体系看起来都完成了 relax,但一个残余力很小,另一个还在慢慢摆动,直接拿最后一步能量比较会引入不必要误差。

静态单点的作用,是把最终 CONTCAR 固定下来,用统一的 ENCUTKPOINTSEDIFFISMEAR/SIGMA 和电子结构输出设置重新计算总能。relax 负责找到合理几何,static 负责给最终结构一个统一能量标尺

(1)  Eads = Eads/static – Eslab/static – Egas/static

1说明静态单点前必须先确认最终结构可靠,否则后面的能量和电子态都只是给错误构型做精修。

VASP计算吸附体系为什么要做静态单点?

1  电催化体系结构与活性位关系。DOI: 10.1039/d6sc02712a

VASP计算吸附体系为什么要做静态单点?

哪些输出更依赖 static?

功函数和表面电势尤其依赖 static 输出。LOCPOT 通常要在固定结构上计算,且最好配合偶极修正和足够真空层。若用 relax 过程中某一步的电势文件,平面平均曲线可能同时包含结构未收敛和电子未收敛的误差。

COHP 或投影态密度也类似。结构刚弛豫结束时,电子态可能已经可用,但为了不同体系可比,仍建议统一用 static 目录生成 DOSCARvasprun.xml 或后续投影文件。电子结构图最怕来源目录混乱

吸附能、差分电荷、BaderPDOS、功函数、平面平均电势、COHP 等结果都更适合基于 static 文件整理。因为这些后处理要么依赖高质量电荷密度,要么依赖更稳定的能带和态密度输出。

如果一边用 relax 的 CHGCAR,一边用 static 的 OUTCAR,再用另一个目录的 DOSCAR 画图,后期很难复查。能量、结构和电子态最好来自同一套目录逻辑

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

静态单点应该保持哪些设置一致?

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

INCAR 和 KPOINTS 怎么处理?

static 通常设置 NSW = 0,固定离子,只重新收敛电子结构。正式比较吸附能时,slabads 和 gas 的泛函、赝势、ENCUTKPOINTS值、色散修正和偶极修正策略要一致;否则差值里会混进参数差异。

金属体系常用合适的 ISMEAR 和 SIGMA,半导体或分子体系则要更谨慎。若前后使用不同 smearing,能量和 DOS 解释都可能偏。吸附能比较中,参照体系必须来自同一把计算尺

(2)  NSW = 0

2可以放在参数一致性部分:电子结构变化要从同一套 static 设置里读,不能把不同精度的输出混到一张图里。

VASP计算吸附体系为什么要做静态单点?

2  金属载体相互作用与电子结构关系。DOI: 10.1039/d6sc02712a

VASP计算吸附体系为什么要做静态单点?

气相分子也要 static 吗?

气相分子还要单独做频率或热修正,尤其是自由能路径。static 只给电子能,ZPE 和熵修正通常来自频率计算或文献数据。若把 static 总能直接当自由能使用,需要明确说明近似范围。

对于带电或偶极较强的吸附体系,static 还应检查偶极修正和真空层是否足够。若 slab 两侧不对称,LOCPOT 平台和总能都可能受电场影响。此时只做 static 还不够,必须保证模型边界条件本身合理。

如果体系含强关联金属,值也要在 slabads 和参照体系中保持一致。不要在 clean slab 用一套 U,在吸附态为了收敛又换另一套 U;这样吸附能差值会混入人为参数变化。

要。吸附能里的 gas 参照不能随手从其他项目拿一个 OUTCARH2O2COCO2H2等分子要单独建盒子,并使用与表面体系一致的泛函、ENCUT 和赝势。O2 还要注意自旋。

气相分子的 K 点通常可以是 Γ 点,但这不代表其它设置可以随意。若 gas 参照来自不同泛函或不同赝势,吸附能再精确也没有意义。

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

电子结构图为什么更依赖静态计算?

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

PDOS、Bader 和差分电荷分别需要什么?

PDOS 需要 LORBIT 和足够稳定的态密度输出;Bader 需要高质量 CHGCAR 或 AECCAR;差分电荷需要 adsslab 和 adsorbate 在统一网格上相减。static 阶段正好可以把这些输出设置统一起来。

这也是为什么很多工作会把 relax 和 static 分成两个目录。relax 目录保存结构演化和收敛过程,static 目录保存最终能量、电荷密度和电子结构。后处理文件越多,越需要目录分层

3可以说明结构调控和电子性质常常一起出现;如果 static 输出混乱,机制解释就很难复查。

VASP计算吸附体系为什么要做静态单点?

3  多维结构调控与催化性能关系。DOI: 10.1039/d6sc02712a

VASP计算吸附体系为什么要做静态单点?

什么时候可以先不做 static?

早期构型粗筛可以先看 relax 最终能量,快速淘汰明显不合理结构。比如几十个吸附初始构型同时弛豫时,可以先根据最终结构和大致能量筛掉跑偏构型,再对保留候选做 static

但进入正式表格、自由能路径或电子结构讨论前,建议统一补 static粗筛可以用 relax,正式报数应回到 static

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

目录和结果表怎样整理?

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

relax 和 static 两级目录怎么排?

推荐每个构型至少保留 relax 和 static 两级目录relax 里放 POSCARCONTCAROUTCAROSZICAR 和结构轨迹;static 里从最终 CONTCAR 复制 POSCAR,再保留最终 OUTCARvasprun.xmlCHGCARDOSCARLOCPOT 等后处理文件。

结果表也要只读 static 能量。表格里可以另列 relax 是否跑偏、最大残余力、是否补频率、static 是否收敛。这样后期发现某个点异常时,可以顺着目录回查,而不是在一堆 OUTCAR 里猜哪个才是最终数值。

4说明结构、能量和性能必须能回到同一套文件记录中复查。

VASP计算吸附体系为什么要做静态单点?

4  催化结构调控与反应性能关联。DOI: 10.1039/d6sc02712a

VASP计算吸附体系为什么要做静态单点?

静态结果异常时先查什么?

还有一种常见异常是吸附体系 static 后电子结构不收敛,但 relax 可以结束。这时不要急着相信 relax 最后一步能量;应检查是否需要更合适的 smearing、是否存在磁矩翻转、是否需要提高 NELM,或者是否初始波函数继承不合适。

目录命名也可以更细一点:relax_01 用于初始构型优化,static_E 用于最终能量,static_DOS 用于更密 点或更适合 DOS 的设置。只要清楚标注能量表读取哪个目录,就不会把 DOS 专用计算和能量比较混在一起。

若要做 Bader 或差分电荷,static_DENSITY 目录可以打开 LAECHG 或保存更密电荷网格。不同后处理可以有不同 static 目录,但能量比较必须有统一来源

如果 static 后发现电荷密度或 DOS 明显异常,也要回头看结构是否已经跑偏。静态单点不是替代结构检查,而是在结构可信之后给出统一电子解。static 的前提是 CONTCAR 本身已经合理

如果 static 能量和 relax 最后一步差异异常大,先查 POSCAR 是否来自最终 CONTCARMAGMOM 是否合理、ISMEAR/SIGMA 是否更换、KPOINTS 是否一致。若电荷密度不收敛,还要检查体系是否存在金属态、带隙、小磁矩翻转或带电 slab 问题。

不要看到 static 失败就直接改很多参数。更稳的处理是一次只改一个关键设置,并记录原因。这样最终结果既能收敛,也能解释为什么采用该设置。

正式交付结果前,可以做一个小检查表:static 是否从最终 CONTCAR 开始、是否 NSW = 0、是否电子收敛、是否保留 CHGCAR/LOCPOT/DOSCAR、能量表是否只读指定目录。这个检查表比事后追 OUTCAR 更省时间。

如果后续还要做频率计算,可以从 static 后的最终结构继续准备频率目录,但不要把 frequency 的能量再混入吸附能表。频率目录主要服务于 ZPE 和熵修正,能量零点仍应回到统一 static 表。staticfrequency 和后处理目录的职责要分开

5说明static 不是多跑一步,而是把最终能量和电子结构放到可复查的标准流程里。

VASP计算吸附体系为什么要做静态单点?

5  静态计算与结果复查流程示例。DOI: 10.1039/d6sc02712a

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

本章要点总结

VASP计算吸附体系为什么要做静态单点?
VASP计算吸附体系为什么要做静态单点?

要点1relax 用于找结构,static 用于统一最终能量和电子结构。

要点2:正式吸附能、自由能和电子结构图最好来自 static 目录。

要点3slabads 和 gas 的 static 设置必须保持同一标准。

要点4:粗筛可先看 relax,正式报数应回到 static

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 2小时前
下一篇 2小时前

相关推荐