大家好,这一阵一直忙乎表面吸附分子的第一原理计算。因为表面一层就有100个原子,要做DFT计算就只建了一层的表面,结果还是非常的慢。总算完成了一个吸附计算,但发现表面原子的起伏非常大。因为担心与真实的表面有较大出入,我把这层表面做了固定,就是下面这样的设定:
%block Geometry Constraints
position from 1 to 100
%endblock Geometry Contraints
但算了几天,我发现表面的原子位置并没有完全固定只是坐标变化小了一些,这是怎么回事?根据以前用其他程序的经验。我固定原子后再进行计算时该原子坐标会保持不变。这与我想像的不一样,想请问各位专家给予解答,不胜感激!