求助关于搭载在碳毡上的材料，如何选取原子冻结范围？

2877321934

各位老师，我在进行Cu-HTTP的电化学问题计算，优化好的结构如下

已知在实验上，这个材料本身是负载碳毡电极表面的，先前了解到需要固定一部分原子以模拟这种负载，我起初冻结的部分如下，将图中选中的黄色部分冻结

然而，这样处理体系，当我尝试优化OH或者O吸附在中心Cu原子上时，会出现较为严重的扭曲，如下图，说明冻结的部分不合理

我先前也翻阅到一篇文章研究体系跟我类似，但是他的做法是将中心一块区域的原子固定，即将下图红框中部分冻结

现在有点拿不准到底应该如何冻结我这个体系了，能否直接将这整个优化好的材料都进行冻结呢？还是说应该效仿这篇文章一样冻结中心部分呢？