sobereva 发表于 2025-11-19 00:30
我不用M$，只能说，如果你的结构已经做了能量极小化（而且没带限制条件，因而结构可以自发地充分弛豫），而 ...

谢谢sob老师，我的结构是有限制条件的，因为想保证嵌套结构的存在，设置了左右两个端头一个胞长度的笛卡尔坐标constrain。我还有计算CNT-CNT嵌套结构，是负值且趋势正确，用同样的方法这个结构就不行，本意是想把这两个结构进行对比的

我不用M$，只能说，如果你的结构已经做了能量极小化（而且没带限制条件，因而结构可以自发地充分弛豫），而且又没有空间原因导致结构无法充分弛豫开（比如位阻作用），则算出来的相互作用能为正，肯定数据处理错了或者计算有硬伤。