|
|
各位老师、同学好，我想请教一个关于 VASP 表面催化体系频率计算中 partial Hessian 区域选择的问题。
我的体系是 Pd/CeO2 催化甲烷 C-H 活化相关模型，包含表面氧空位和次表层氧空位。体系总原子数约 140-150 个。前期已经完成了 IS 和 FS 的结构优化，现在想做频率计算，主要目的是确认 IS/FS 是否为局域极小态，同时获得与吸附物和表面反应区域相关的振动信息。
由于体系中含 CeO2，计算中考虑了 Ce 4f 电子的强关联效应，采用 DFT+U 处理 Ce 4f。主要设置如下：
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = -1 3 -1 -1 -1
LDAUU = 0.0 5.0 0.0 0.0 0.0
LDAUJ = 0.0 0.0 0.0 0.0 0.0
LMAXMIX = 6
LASPH = .TRUE.
ISPIN = 2
元素顺序对应 O Ce PD C H，即只对 Ce 的 4f 轨道加 U，Ueff = 5.0 eV。结构优化阶段使用了相同的 DFT+U 设置，并根据体系情况设置过 MAGMOM；频率计算中继承结构优化 INCAR 中的 MAGMOM 设置，不重新自动生成。由于涉及表面/次表层氧空位，部分 Ce 可能存在 Ce3+/Ce4+ 局域化和磁矩分布差异。
频率计算设置大致如下：
IBRION = 5
NSW = 1
NFREE = 2
POTIM = 0.015
EDIFF = 1E-06
ISYM = 0
我使用 POSCAR 的 Selective dynamics 做 partial Hessian，只放开部分原子，其余 slab 原子固定。根据 VASP Wiki，IBRION=5 下 Selective dynamics 会只计算标记为 T 的 Hessian 分量；NFREE=2 时每个可动原子大约需要 6 个位移计算。因此可动原子数对计算量影响很大。
目前我的处理方式是按 z 方向分层，只放开最上方 5 层原子。用类似 VASPKIT 402 的分层方式判断，阈值 0.02 时整个 slab 大约可以识别出 8-12 层。放开 top 5 层后，可动原子数大致为：
- IS：约 34-67 个可动原子；
- FS：约 53-99 个可动原子。
这个区域通常包括 CH4 或 CH3+H、Pd、表面 O/Ce，以及最靠近吸附物的表面反应区域。我的反应构型中，不论氧空位位于表层还是次表层，C-H 活化时 C-H 基本都是对准表面 O，因此我初步认为主要振动自由度应集中在吸附物、Pd 和表面 O 附近。
但我担心的问题是：有些模型涉及次表层氧空位。如果只放开 top 5 层，实际可能主要包括吸附物、Pd、表面 O/Ce 等区域，而没有充分放开次表层氧空位附近的 O/Ce。如果再增加一层 O/Ce，可动原子数会显著增加，频率计算量也会大幅增加。
因此想请教大家：
1. 对这种 DFT+U 的 Pd/CeO2 表面甲烷活化体系，用 partial Hessian 只放开吸附物 + Pd + 表面反应区域 + 表面几层原子，是否可以作为判断 IS/FS 局域极小态的合理近似？
2. 如果存在次表层氧空位，是否有必要把氧空位附近的 Ce/O 也加入可动区域？还是说只要结构优化时这些原子已经充分弛豫，并且频率计算保持相同的 DFT+U 与 MAGMOM 设置，主要关注吸附物和表面反应 O/Pd 区域通常就可以接受？
3. 从经验上看，对于表面反应频率验证，是更推荐：
- 按 z 方向统一放开 top N 层；
- 还是手动选择“吸附物 + Pd + 成键表面 O/Ce + 氧空位近邻 Ce/O”作为可动区域？
我目前的想法是先用 top 5 层（只包含CH4和上面一层O和一层Ce批量计算，后续如果某些含次表层氧空位模型出现虚频，或者关键构型比较敏感，再扩大可动区域复算。想听听大家是否认为这个策略合理，或者是否有更规范的处理方式。
感谢各位指教！
|
|
发表于 Post on 3 day ago
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author