[Kwant] Coupling Matrix between Lead and Conductor

2018-06-06 Thread Shivang Agarwal
Hello authors, I am trying to perform an eigenchannel analysis of a graphene nanoribbon. For that I will be using the formula : *T(E) = ГL(E)½ GC†(E) ГR(E) GC(E) ГL(E)½ * where *ГL(E)* is the coupling matrix between the left lead and the conductor, *GC(E)* is the greens function matrix of the cond

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-06 Thread Shivang Agarwal
Just a correction: I am interested in finding the transmission matrix 't' from the following: t†t = *ГL(E)½ GC†(E) ГR(E) GC(E) ГL(E)½ * Best Regards, Shivang On Thu, Jun 7, 2018 at 2:23 AM, Shivang Agarwal wrote: > Hello authors, > > I am trying to perform an eigenchannel analysis of a graphe

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-06 Thread Abbout Adel
Dear Shivang, To get the matrix Gamma you can do: sys= sys.finalized() lead_L = sys.leads[0] Sigma_L= lead_L.selfenergy(energy) Gamma_L = -2*imag(Sigma_L) You can find the details in this answer by Joseph [1]. To get directly the transmission matrix t you can do: t=kwa

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-07 Thread Xavier Waintal
Hi Shivang, > Le 6 juin 2018 à 22:53, Shivang Agarwal a écrit > : > > Hello authors, > > I am trying to perform an eigenchannel analysis of a graphene nanoribbon. For > that I will be using the formula : > T(E) = ГL(E)½ GC†(E) ГR(E) GC(E) ГL(E)½ > where ГL(E) is the coupling matrix between

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-07 Thread Shivang Agarwal
Thanks Xavier. I now have a clearer picture. Regards, Shivang On Thu, Jun 7, 2018 at 11:33 AM Xavier Waintal wrote: > Hi Shivang, > > > Le 6 juin 2018 à 22:53, Shivang Agarwal a > écrit : > > Hello authors, > > I am trying to perform an eigenchannel analysis of a graphene nanoribbon. > For tha

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-07 Thread Shivang Agarwal
Hi Abbout, Thanks for swift response. Indeed, the kwant.wavefunction module gives me a complex number. I had been working on probability (wavefunction squared) and had overlooked the phase part. A noob mistake. Appreciate your help! Shivang On Thu, Jun 7, 2018 at 11:53 AM Abbout Adel wrote: >

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-08 Thread Shivang Agarwal
Hi, I had one more query. Would it be possible to calculate and plot the magnitude and phase of the complex wave function of each transmission mode separately? Regards, Shivang On Thu, Jun 7, 2018 at 8:29 PM Shivang Agarwal wrote: > Hi Abbout, > > Thanks for swift response. > Indeed, the kwant

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-10 Thread Abbout Adel
Hi Shivang, The transmitted wave function transmitted in the mode m (and coming from the mode n) is: t_mn exp(i k_m). > Since you have the transmitted amplitudes and the wavenumbers k_m, you have everything for your request. I hope this helps. Adel On Fri, Jun 8, 2018 at 11:09 PM, Shivang

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-11 Thread Shivang Agarwal
Hi Abbout, I have a matrix giving me the complex wave function at each site. How do I know which mode (m/n) do the amplitudes correspond to? Also, I am not aware of the method to get wavenumbers k_m. Could you kindly help me out with that as well? Thanks for the help! On Sun, Jun 10, 2018 at 1:4

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-11 Thread Joseph Weston
Hi, > > I have a matrix giving me the complex wave function at each site. How > do I know which mode (m/n) do the amplitudes correspond to? Where did you get this matrix from? Is this the wavefunction of a mode in a lead, from 'kwant.modes', or of a scattering wavefunction, from 'kwant.wave_func

Re: [Kwant] Coupling Matrix between Lead and Conductor

2018-06-11 Thread Shivang Agarwal
Thanks a lot Joseph. That really helps a lot! I now get what I need to do. On Mon, Jun 11, 2018 at 10:14 PM Joseph Weston wrote: > Hi, > > > > > > I have a matrix giving me the complex wave function at each site. How > > do I know which mode (m/n) do the amplitudes correspond to? > > Where did y