I am trying to optimize some molecules (~50 atoms including boron and Fe) with ROHF. In Gaussian16, I always get error /convergence failure while in Q-Chem, I can optimize the molecule, but always get very large imaginary frequencies.

For my initial guess, I am using my HF optimized molecule, still no luck. Any suggestion will be very helpful!