Dear Gabriel, For the stationary case, you can rely on the method that uses a "fictitious lead" with zero self energy. (No integration is needed) Please, check these two posts [1] [2] . I hope this helps,
[1] https://mail.python.org/archives/list/kwant-discuss@python.org/thread/2VB4IKZSCBZZNZJOWS3QFSP2AT6FCYVL/ [2] https://www.mail-archive.com/kwant-discuss@kwant-project.org/msg01356.html On Wed, Sep 13, 2023 at 8:38 PM Gabriel Garcia <gqgarci...@gmail.com> wrote: > Hi devs, > > I have been reproducing section 2.7, “Green’s Functions”. In this section, > the tutorial shows us how to calculate Green’s functions for a 1D > monoatomic network with two impurities. We learn to calculate all different > types of time-dependent and independent Green functions. However, I have a > question. When we consider a stationary system, the Green functions will be > calculated depending on the Energy of the system. In the impurities > example, the Green functions were calculated by varying the onsite energy > of the impurities. How can I calculate the Green functions in terms of the > system’s energy? Remembering that stationary Green functions are integrals > of the system’s energy. > > Best, > Gabriel Garcia > -- Abbout Adel