CP2K如何消除碳纳米管内的水分子的质心漂移？

Alset

老师们好，

我在使用cp2k-aimd模拟碳纳米管里结冰时，固定了碳原子，并且设置COMVEL_TOL和ANGVEL_TOL为0，但仍发现内部水分子存在质心漂移的情况，而碳原子往相反的方向发生了运动。
具体情况如下图所示：


目前认为这些水分子是以冰的形式存在，因为水分子并没有像液态一样发生明显的旋转，这也可以从旋转自相关函数得到印证。但由于水发生漂移，我在计算水分子扩散系数时，得到了不合理的偏大结果。
请问目前设置有什么不合理的地方，以及该如何解决以上问题呢？附件是该模拟的输入、输出以及轨迹文件。