'Vpds': 0.0724, 'Vpdp': 0.0772, 'Vdds': 0.1849, 'Vddp': 0.0429, 'Vddd': - 0.0333, }, 'MoMo': { 'Vsss': - 0.0575, 'Vsps': 0.0057, 'Vpps': 0.0296, 'Vppp': 0.0946, 'Vsds': - 0.1082, 'Vpds': 0.0212, 'Vpdp': - 0.0448, 'Vdds': - 0.0216, 'Vddp': - 0.0285, 'Vddd': 0.0432, } } system.with_eigenvectors = True system.num_of_bands = 8 for i in xrange(len(system.atoms)): system.atoms[i].orbitals = ['s', 'px', 'py', 'pz', 'dxy', 'dyz', 'dxz', 'dx2-y2', 'dz2'] system.just_do_main_magic() plt = Plotter(system.name) plt.plot_energy_bands_from_file()