Dear Hu Yongliang, > Hello, I am interested in Green Function, so I read related Kwant codes > (kwant.solver.greens_function, make_linear_system, solve_linear_system, > etc.). I found a lot of optimizations in Kwant, such as LU Decomposition and > Matrix Partition, which make Kwant more efficient. To calculate Green > Function, we need to calculate the self-energy of leads. Though I know how > to get self-energy by using command "lead.selfenergy(energy)". I want to > know how this command works, the codes in "leads.py" are too diffcult to me.
The inverse of the left-hand side matrix produced by make_linear_system in Kwant is the Green's function, i.e. it already includes the relevant self-energy. However computing the complete Green's function is extremely inefficient. Can you please be more specific: do you perhaps need only some matrix elements of the Green's function? What observable is your final goal? Best, Anton > > > Sincerely, > > > Hu Yongliang > > > >
