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