Q1：朱老师，反磁材料计算，mag越跑越大，这是怎么回事呢？A：
①mag 越跑越大，说明体系在自发找磁性解，不像真正反磁、非磁态；
②若是抗磁、非磁，应设 ISPIN=1，不要给 MAGMOM；
③若是反铁磁，总磁矩可接近0，但局域磁矩会一正一负；
④只看总 mag 不够，要查 OUTCAR 里各原子 magnetic moment；
⑤不该有磁性，删 WAVECAR 重跑，并检查奇数电子、缺陷、U值、过渡金属设置等参数。
Q2：朱老师，计算过渡态的初态、模态模型的原子数目是不是必须对应，比如初态有C3H3，末态也得有C3H3才可以计算出过渡态？
A：必须对应，否则不能正常做过渡态、NEB：
①初态和末态原子数必须相同，如初态 C3H3，末态也必须含 C3H3；
②原子顺序也要一致，不能中途增删或调换原子编号；
③反应中“脱附、迁移”的原子，也要保留，只是位置改变；
④若反应前后确实有分子进出，要把环境物种一起放进模型。
A：正常，吸附 OH 后总 mag 变化很常见：OH 会引起表面电荷转移，改变金属局域磁矩；
OH 的 O/H 轨道可能与表面 d 轨道杂化，导致自旋重排；
只看总 mag 不够，要看各原子局域磁矩和差分电荷;
若变化很大，需检查不同磁态初猜，确认不是收敛到亚稳态。
Q4：朱老师，我有一个Cu2O（111）面的slab，结构优化不收敛，出现了这种情况：ZBRENT: can’t locate minimum, use default step
ZBRENT 是离子步线搜索失败，常见原因是初始结构差、原子过近、表面重构太剧烈;
Cu2O(111) 表面较敏感，先检查终止面、真空层、固定底层、是否有悬挂键异常；
参数上建议改稳：IBRION=2，POTIM=0.1~0.2，ISIF=2，NELM=120；
若仍不收敛，先固定更多底层或分步优化，先粗收敛，再放开表面继续优化。
Q5：朱老师，想请教一下我想算图中吸附分子的偶极矩请问如何计算呢？是通过Bader电荷分析获得电荷量之后再乘以距离吗？
A：不建议只用 Bader 电荷 × 距离，那只是近似估算
若算吸附分子自身偶极矩，应把吸附后的分子单独取出，保持构型，放大真空盒静态计算；
Bader 电荷乘距离只能估算电荷转移偶极，不能代表完整偶极矩；
若算界面偶极，应用差分电荷，再沿表面法向积分，这比 Bader 更合理；
偶极矩是整体结果，通过加偶极修正，软件会算出来，outcar中搜索dipole，结果是exA为单位的数值。
