计算化学公社

标题: 有关LBO和Mayer分析的问题 [打印本页]

作者
Author: dzdhp 时间: 2023-6-21 10:24
标题: 有关LBO和Mayer分析的问题
用ORCA跑了一个有机弱酸和碱之间的质子转移，方法用的是B973C，此过程氢从氧转移到氮上，我保留了每帧的波函数文件，按照sob老师教的脚本用Multiwfn计算了每帧的mayer键级和LBO，然后绘成图（图中红线为O-H，黑线为N-H），从图上看mayer键级还算比较合理，在200ps左右O-H断裂形成N-H，并且之后比较稳定。但是LBO就有点奇怪了，N-H持续波动且范围较大，sob老师原文里N-H的LBO值在0.86左右。请教一下各位老师导致LBO波动的可能原因是啥，是因为转移过去后共用电子都是N提供的吗，我该如何解释此波动，还是说这本来就是正常的。

作者
Author: sobereva 时间: 2023-6-21 14:11
可能是默认的积分格点不够多，导致结构变动时数值误差波动大所致。把Multiwfn的settings.ini里的iautointgrid设为0，然后尝试，如果波动还大，把radpot和sphpot都适当改大再试。如果改大后还是波动大，可能是LBO计算涉及的原子间重叠空间区域随结构变化变动得比较厉害，这时候可以在比如origin里对数据做个smoothing对相邻帧几帧数据取平均来得到平滑化的曲线。

作者
Author: dzdhp 时间: 2023-6-21 15:02

sobereva 发表于 2023-6-21 14:11
可能是默认的积分格点不够多，导致结构变动时数值误差波动大所致。把Multiwfn的settings.ini里的iautointgr ...

感谢sob老师百忙中抽空回复，我这就去试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)