We can set a device configuration for the parameter "configuration". But what will we get with MolecularEnergySpectrum for a two-probe system? The MPSH? We diagonalize the the Hamiltonian of the central region? Since the Hamiltonian is no Hermitian for a two probe system, what will be the wave functions? The real part?