Is it possible to perform ROMP2 numerical optimizations in Gaussian09?
I've seen this error when I want to combine an optimization with an orbital deletion (NBODel) in Gaussian. It cannot handle an optimization with more than 50 degrees of freedom, so you need to freeze as much as you can, or pre-optimize your molecule with a more traditional method and then freeze whatever you can and do the ROMP2 optimization on the parameters of interest.