光写maxstep根本没用，必须再加上notrust
另外，势能面平坦的大体系建议加上gdiis
明显http://sobereva.com/164文中的知识你还没充分掌握，再一个字一个字反复多读几遍此文

opt=maxstep=3甚至更低都是可以考虑的。

可以考虑opt=modredundant，手动加几个冗余内坐标，算几个冗余内坐标的数值Hessian。但这个我没操作过，只是手册上看到过。

如果计算量和内存承受得起，建议用opt=recalc=x，效果是每x步计算一次Hessian,一般3-5比较合适。
再降低就是opt=calcfc只算一次Hessian。不过我感觉300多原子的Hessian计算量会很大很大，慎用。