Commit c460cab8 authored by Joseph Weston's avatar Joseph Weston
Browse files

Merge branch 'fix/selfenergy-clarification' into 'master'

Add a clarification to the self-energy calculation for leads

See merge request kwant/kwant!353
parents 75ef39f9 88d3d09e
......@@ -2,6 +2,7 @@ Thomas Kloss <kloss@itp.uni-frankfurt.de>
<joseph@weston.cloud> <joseph.weston08@gmail.com>
<joseph@weston.cloud> <joseph.weston@cea.fr>
<joseph@weston.cloud> <me@josephweston.org>
<joseph@weston.cloud> <v-josewe@microsoft.com>
Tómas Örn Rosdahl <torosdahl@gmail.com>
<slavoutich@ostroukh.me> <ostroukh@lorentz.leidenuniv.nl>
<m.t.wimmer@tudelft.nl> <wimmer@lorentz.leidenuniv.nl>
......
......@@ -194,8 +194,11 @@ class StabilizedModes:
there is a propagating mode, `Sigma` will definitely be complex.)
"""
v = self.sqrt_hop
vecs = self.vecs[:, self.nmodes:]
vecslmbdainv = self.vecslmbdainv[:, self.nmodes:]
# Only include the *outgoing* modes (propagating + evanescent),
# meaning this is the *retarded* self-energy.
outgoing = slice(self.nmodes, None)
vecs = self.vecs[:, outgoing]
vecslmbdainv = self.vecslmbdainv[:, outgoing]
return dot(v, dot(vecs, la.solve(vecslmbdainv, v.T.conj())))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment