分子转动能垒计算

dixin

各位老师，导师让模拟一个体系的转动能垒（转子为中间的卟啉分子，外围有柔性链，整个体系有手性，转子部分无手性），模型分子见附件。
我采用两种方式进行限制性优化：
（1）固定两个中心金属原子和周围的12个氮原子。
（2）固定两个二面角。
步长均为5度。
优化后，计算单点能。
理论上来说，转动以后的能量总是应该低于转动之前的能量。但是两种固定方式，都出现转动以后能量反而降低的问题，就我查阅的文献来看，这种结果是不能用的。不知道该怎么解决了，所以在这里想向大家求助。麻烦各位大神，帮忙看一下。
ps:方法为b97d，基组为混合基组。优化使用的基组为6-31G（d）等等，单点能使用的基组为6-311G（d）等。转动是通过二面角（18-109-110-51）实现的，具体的固定原子和二面角见附件。