求助CP2K优化Fe2O3板模型吸附亚砷酸，优化过程中发生自发质子转移

2429783189

我在使用CP2K程序计算一个Fe2O3表面吸附亚砷酸的体系，在优化过程中，原本连在亚砷酸的O上的H自发转移到Fe2O3表面的O上，下面第一张图是初猜结构，第二张图是收敛结构 我尝试了三种初猜，亚砷酸均在优化过程中发生脱质子，甚至其中一种结构还出现优化过程中质子飘走（下面第一张图是初始，第二张图是收敛结构） ，想请问一下老师这种情况是否说明我的计算有问题呢？可以直接用当前情况下使用我优化出来的结构的E（Fe2O3-亚砷酸）-E（Fe2O3）-E（亚砷酸）来计算Fe2O3吸附亚砷酸的吸附能吗？我的三种初猜的输入文件inp，优化轨迹xyz，输出文件out均以上传至附件（三种初猜以As1，As2，As3文件命名区分）

ultramanm87

先不考虑你算的对不对，吸附过程中质子转移本来就很正常，在水溶液环境中这种表面吸附无机离子本来就是高度pH相关的，具体的可以了解下表面络合模型的有关理论。
回到你计算的设定，你得想好你模拟的到底是不是溶液中的吸附，如果是，得在一开始加溶剂模型优化或者做aimd是更加合适的。