How to get surface Green function of ideal semi-infinite lead in real space when I have the detail structure of the lead and tight-binding Hamiltonian . In effect, I have tried the transfer matrix method to calculate it. However, there are some problem. is there some other efficient method?