My system contains more than two radicals and I want to calculate all possible magnetic exchange coupling constants. I am using g09 for optimization. Orca 3.0.1 is intended to be used for single point and broken symmetry calculations since it has the feature of "Flipspin". Will this approach be giving desirable and satisfactory results?