金属表面吸附能计算方法提问

GHR1111

本人一般计算吸附能的步骤是先优化切好面的金属基底（完全弛豫状态）,在单独优化需要吸附的小分子，然后再优化整个吸附构型，在优化整个吸附构型的时候，对于基底我又会固定下面几层，只允许表面2-3层弛豫，已知最后计算公式是E=E总-Emolecular-Ebase，那如果我构建吸附构型的时候选择了固定基底的下面几层，那我单独优化基底的时候是不是也应该对应固定相应的层数计算能量，最后计算的能量才能对应呢？以及我的步骤是否有问题呢？

阿夸夸

常规的计算流程，没发现明显问题