In tight binding analysis, we derive a hamiltonian representing the system (this derivation is usually manual).
My question is that by using this tight binding hamiltonian, how can I plot a band structure working with any DFT code like VASP or Wien2k, etc.?