优化一个70多原子的分子（只有CHON,Mg），用DFT-B3LYP方法，先用3-21g优化至收敛后，再用6-31g(d)优化，很快出现了能量震荡（见图1,2）不收敛，手动终止后选取最后一步的结构，用opt=(maxstep=1,calcfc,NoTrust)接着算，出现了力常数减小时能量在上升，和力常数上升时能量减小的情况（见图3,4），请问我接下来应当怎样做才能使结构收敛？刚注册账号没多久，还没有太多金币悬赏，先谢谢大家了。 返回小木虫查看更多