Hi,

First, NumForce converts the MOs at the original geometry to C1 symmetry. Then, these MOs are used as an initial guess for the displaced geometries (for details, see the script "runsingle", which runs the actual displacement calculations). So, I don't see any obvious options for further optimization in this respect.

Best,

Antti