氧气在K掺杂石墨相氮化碳上找不到合适吸附位点

phchhh

各位老师好。我在使用Cp2k计算K掺杂石墨相氮化碳的二电子步ORR,泛函使用PBE，添加了DFT-D3(BJ)色散矫正,基组选择DZVP-MOLOPT-SR-GTH,以下是我遇到的问题。
在计算吸附氧气时，我首先先找到一个靠近K的边缘C作为吸附位点，氧气成功被活化了（如图1，坐标文件为O1.xyz），此时氧气键长为1.45Å,底端氧气和吸附位点C的距离为1.26Å。


                            （图一）
此时我用算好的结构添加H计算OOH的结构，发生了OOH的分解（如图2，坐标文件为OOH.xyz）

                            （图二）
因为OOH分解，曾设置过限制势来固定O-O键长，结果每一步SCF都能收敛，但是受力一直不达标，出现振荡。但结构是合适的，图三是我kill掉任务后的最后一帧几何优化结构，

                            （图三）
输入文件中关于限制势的部分是
&CONSTRAINT

   &COLLECTIVE

      COLVAR 1                     

      TARGET [ANGSTROM] 1.377

       INTERMOLECULAR

       &RESTRAINT

        K [kcalmol^-1*angstrom^-2] 500.0   

      &END RESTRAINT

   &END COLLECTIVE

&END CONSTRAINT

我怀疑这个位点不是合适的吸附位点。首先氧气吸附后被活化到1.45Å，属于过氧根的键长，且吸附位点与氧原子距离过近，形成了双键，我怎么调试初始结构最后OOH都分解。

各位老师，我主要有两个问题，第一是如果OOH的分解是否可以通过有效方式解决，使OOH成功吸附，或者通过修改限制势来解决？第二，如果最初的吸附位点并不是有效的，我该怎么找到合适的吸附位点呢？烦请各位老师解惑，感谢。