I haven’t studied your definition in detail, but a couple of things I notice -
You have the energy threshold set to 1e-40
I don’t think there is ever any need to use anywhere near this low a value - the radius of an electron divided by the radius of the earth is several orders of magnitude bigger than that. The errors from the floating point operations will be bigger than that, so you are essentially guaranteeing the solver will never consider the result converged.
The other thing is the ratio of the strengths of the goals. The ratio of the strongest to weakest goals is very high - and this will slow down convergence.