The effective Hamiltonian is obtained by using a reverse-engineered inverse density functional theory (iDFT), or more precisely an inverse site occupation function theory (iSOFT), combined with an adaptation of the kinetic energy terms. The determination of the effective parameter utilizes a steepest-descent gradient search. The ikDFT is currently not part of the public interface.