if $guessed_energy >= 0.0 it is calculated internally from the hydrogenic formula
The guessed_energy for valence orbitals better be good, since the program gives some dumb answers if the hydrogenic default is used.$iof governs some internal step, the mixing weight between previous and next iteration