求助模拟气体分子穿透晶体薄膜过程中是否需要对晶体原子位置设置限制

cheng1

我想使用Gromacs模拟气体分子穿透团簇晶体薄膜，构建了如图所示体系，在NVT系综下进行了分子动力学计算（NPT系综会崩溃），其中团簇晶体中的16个重原子Au被固定。顶部和底部各使用惰性原子Ar（组名为WAL）并在模拟中冻结，以防止气体分子从盒子边界处跑出到真空区。
在模拟过程中存在如下问题：
1.气体无法穿透晶体，是不是应当取消团簇晶体中重原子的位置限制？（PS：一开始全部原子均未固定，但是在运行过程中晶体结构中的8个团簇分子排列一边发生错位一边向真空区一侧运动，导致真空区消失，因此给团簇中的Au原子施加了位置限制）

2.体系运行了100000000 步， 200ns，能否说明气体分子难以穿透晶体薄膜？

3.是不是我的晶体体系过厚，需要减少一层晶体原子吗？

以下是相关文件，恳请老师指教