LOBSTER投影质量怎么看?

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

引言

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

LOBSTER 能把 VASP 的平面波结果转成 COHP、ICOHP、Mulliken 电荷和局域投影 DOS,但可靠性取决于 WAVECAR/CHGCAR、结构端点和局域基组是否匹配。LOBSTER 结果首先要看投影质量,再讨论某条键的成键、反键和积分强度。charge spilling、基组选择、自旋设置和能量对齐都会影响 ICOHP 的可解释范围。

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

LOBSTER 到底在投影什么?

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

平面波结果怎样转成局域轨道?

VASP 的波函数以平面波形式表示,LOBSTER 需要把这些波函数投影到原子局域基组上,再计算轨道重叠、哈密顿布居和键合贡献。在这一投影流程中,自洽计算提供本征值和波函数,局域基组决定后续 COOP/COHP 与电荷分析的分辨率。

投影过程还要求计算端点和后处理端点完全一致。POSCAR 中原子顺序、POTCAR 中元素顺序、WAVECAR 和 CHGCAR 对应的结构都不能混用。若优化后又改了原子排序或移动了坐标,LOBSTER 输出的键对编号就可能对应错误对象。

因此,ICOHP 不是新的自洽能量,而是对已有 VASP 电子结构的后处理分解。它能回答某对原子之间哪些能区以成键或反键贡献为主,却不能独立替代总能、吸附能、NEB 势垒或结构稳定性。若 VASP 静态计算本身未收敛,或 WAVECAR/CHGCAR 与 POSCAR 不对应,LOBSTER 后处理也不会给出可信解释。

LOBSTER投影质量怎么看?

图1 LOBSTER 从平面波波函数投影到局域基组并计算 COOP/COHP/电荷的原理示意。该图用于说明投影质量是后处理的起点。 DOI: 10.1002/jcc.24300。

charge spilling 说明什么?

charge spilling 描述平面波波函数中有多少成分未被所选局域基组充分表示。数值越大,说明投影缺失越明显。不同元素、价层轨道、赝势和基组组合都会改变这个指标;对含 d/f 轨道、强杂化或高能未占据态的体系,默认基组并不总能覆盖需要分析的轨道空间。

不同基组下的投影密度与 PAW 密度差异,能直接暴露局域轨道覆盖是否不足。charge spilling 过大时,ICOHP 的绝对值和轨道分解都需要重新检查。后处理开始时先比较 lobsterout 中 total spilling 和 absolute spilling,再检查 DOSCAR.lobster 是否能重现实验关注能区或 VASP 原始 DOS 的主要峰形。

LOBSTER投影质量怎么看?

图2 不同局域基组下 PAW 密度与 LOBSTER 投影密度差异。该图说明基组不完整会影响投影质量。 DOI: 10.1002/jcc.24300。

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

基组选择怎样影响 ICOHP?

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

哪些轨道要放进去?

基组选择要围绕化学问题。分析过渡金属氧化物 M-O 键时,金属 3d/4d 和 O 2p 是核心;分析 CO 吸附时,C/O 的 2s、2p 与金属 d 态都要被包含;分析卤化物或重元素体系时,还要确认相应 p/d 轨道是否进入基组。

对单原子催化剂,金属中心、配位 N/C/O 原子和关键吸附物都应进入同一套投影讨论。若只分析 M-N 键而忽略 M-ads 键,容易把锚定稳定性和反应中间体成键混在一起;若只分析吸附物内部键,又看不到载体如何调控金属中心。

如果只用过小基组,某些反键态可能被投影得很弱,ICOHP 积分会偏离真实键合趋势。若基组过大,也可能引入线性相关或投影不稳定。合适基组应同时降低 spilling 并保持 PDOS 峰形合理,而不是单纯追求更多轨道。

自旋和磁矩要不要分开看?

磁性体系必须把自旋设置纳入 LOBSTER 分析。图3 的 α-Fe 对比显示,自旋极化会改变 COHP 分布和交换分裂,非磁性后处理不能解释磁性自洽结果。对 Fe、Co、Ni、Mn 氧化物、氮化物、单原子位点和 Li-S 催化剂,自旋分辨结果常常是键合解释的关键。

若两个磁构型总能接近,应分别做 LOBSTER 后处理,而不是只分析最低能态。这样可以判断键合差异是否伴随磁矩翻转出现。

若 VASP 计算中 ISPIN = 2,LOBSTER 输出也要按自旋通道读取 COHPCAR.lobster 和 ICOHPLIST.lobster。上自旋与下自旋的 ICOHP 差异可能对应不同轨道占据。这类差异应与原子磁矩、PDOS 和差分电荷密度一起解释,不能用两个自旋通道的合并值直接概括键合增强。

LOBSTER投影质量怎么看?

图3 α-Fe 中非磁性与自旋极化 COHP 对比。该图说明自旋设置会改变化学键分析结果。 DOI: 10.1002/jcc.24300。

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

ICOHP怎样和吸附能对照?

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

ICOHP 更负就一定吸附更强吗?

ICOHP 通常越负代表该原子对之间成键贡献越强,但吸附能包含多个键、结构变形、表面重排、分子内键变化和长程相互作用。若 CO 与金属之间 ICOHP 更负,同时 C-O 内部反键占据增强,反应活化可能提高,但吸附能未必同步变得更负。

以 *OOH、*COOH 或 Li₂Sₓ 为例,一个中间体常同时存在 M-O、M-C、Li-S、S-S 等多类相互作用。某一根键的 ICOHP 变强,可能伴随另一根键削弱。若反应关键在 O-O 或 S-S 断裂,就要把内部键的 COHP 变化放到同一张图里,而不是只看表面接触键。

表面模型中的 pCOHP 分析如图4 所示,能够解析低维结构中不同原子对的成键贡献。ICOHP 适合解释局域键合,不适合单独承担整体热力学排序。当吸附能和 ICOHP 趋势不一致时,应把构型畸变能、多个接触键、分子内部键和电荷转移分开列出。

LOBSTER投影质量怎么看?

图4 Ag(001) 支撑 Ru 二维层的表面模型与 pCOHP 分析。该图用于说明表面和低维模型也能进行 LOBSTER 键合分析。 DOI: 10.1002/jcc.24300。

CO 吸附能怎样做案例判断?

CO 吸附是 LOBSTER 很适合处理的例子,因为 5σ 给电子、2π* 反馈、C-O 键长和 COHP 反键占据之间可以互相对应。图5 中 CO 吸附 Co 表面的 lm-PDOS 把 σ 与 π 贡献拆开,适合观察不同吸附位点下分子轨道如何参与成键。

若某位点 CO-C/M ICOHP 更负,同时 CO 2π* 态占据增强、C-O 键长增加、C-O 伸缩频率降低,可说明该位点更有利于 CO 活化。这类判断必须由几何、PDOS、COHP 和能量共同支持。若只有 ICOHP 一个数值,不宜直接推出反应势垒降低。

LOBSTER投影质量怎么看?

图5 CO 吸附 Co 表面时 lm-PDOS 与 COHP 分析。该图说明分子轨道投影、COHP 和反应势垒可以联合解释。 DOI: 10.1021/acs.jpcc.4c00144。

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

VASP 到 LOBSTER 怎样执行?

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

静态计算要准备什么?

LOBSTER 通常需要高质量静态计算结果,包括 POSCAR、POTCAR、WAVECAR、CHGCAR 和 vasprun.xml。静态计算应保持充分收敛,ISYM 建议按体系检查;对表面吸附、缺陷和低维材料,常用关闭对称性来避免等价原子被不合适地约化。LORBIT 与投影 DOS 可作为对照,但不替代 LOBSTER 投影。

静态计算还应保持结构端点清楚。若吸附态优化后发生位点迁移,LOBSTER 分析应针对最终构型重新定义键对;若 NEB 过渡态用于键合解释,则应说明该结构来自路径最高点还是频率确认后的过渡态。不同端点的键对命名应保持一致,便于比较。

VASP 端还要确认 NBANDS 是否覆盖目标能区,尤其是需要分析未占据反键态时。若能量窗口过窄,COHP 或 PDOS 只能解释费米能级附近的一部分状态。参数示例可从 ENCUT 和 k 点收敛后的静态输入开始,并保持与吸附能或 NEB 端点相同的结构。

结果输出怎样检查?

后处理完成后,先读 lobsterout,确认 spilling、基组、投影轨道和总电荷;再读 ICOHPLIST.lobster,检查目标键的积分值;最后用 COHPCAR.lobster、DOSCAR.lobster 和 VASP 原始 DOS 对照峰位。若使用 pCOHP,还要明确分析的是哪对原子、哪一类键和哪一个自旋通道。

对吸附体系,还应把吸附前后的同一根键分别输出。若只看吸附后的 M-C 或 M-O 键,容易忽略分子内部 C-O、N-N 或 O-O 键是否被削弱。CO、NO₃⁻、OOH* 和 Li₂Sₓ 等物种都适合同时比较表面-吸附物键与分子内部键。

当不同键的 ICOHP 趋势相反时,不能只选择其中一条支持结论。更合理的矩阵应列出目标吸附键、关键内部键、吸附能、关键键长、Bader 电荷和反应能垒。若 M-C 键增强但 C-O 键也增强,CO 活化判断就要回到 π* 占据和频率变化再确认。

基组调整后 spilling 下降时,某些 ICOHP 绝对值可能整体改变。此时应比较相对趋势是否保持,而不是把两套基组下的积分值直接混排。对公开结果,basisFunctions、能量积分范围和是否采用自旋分辨输出都需要写明。

LOBSTER 的可靠使用方式,是先证明投影足够代表 VASP 波函数,再用 COHP/ICOHP 解释局域成键。当 charge spilling 偏高、DOS 峰形对不上或自旋态不一致时,应回到基组、自洽输入和结构端点重新检查。由此得到的 ICOHP 才能和吸附能、自由能或反应势垒形成有意义的对照。

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

本章要点总结

LOBSTER投影质量怎么看?
LOBSTER投影质量怎么看?

要点:LOBSTER 是对 VASP 波函数的局域基组投影,投影质量决定后续 COHP/ICOHP 的解释范围。

要点:charge spilling、基组轨道、自旋设置和能量窗口都应在讨论 ICOHP 前检查。

要点:ICOHP 解释局域键合,吸附能包含整体热力学、结构重排和多键贡献。

要点:可靠判断应把几何、PDOS、COHP/ICOHP、Bader 电荷和总能结果放在同一组对照中。

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

相关推荐