t_mn exp(i k_m).
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,ShivangOn Thu, Jun 7, 2018 at 8:29 PM Shivang Agarwal <shivang.agarwal@iitgn.ac.in> wrote: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!ShivangOn Thu, Jun 7, 2018 at 11:53 AM Abbout Adel <abbout.adel@gmail.com> wrote:I hope this helps.Now, since your aim is to get the wave function, the module kwant.wavefunction gives you the wavefunction as a complex number (module and phase). So, your claim that you are unable to get the phase is confusing!t=kwant.smatrix(sys,energy, out_leads=[1], in_leads=[0]).data #I suppose you have just two leads.You can find the details in this answer by Joseph [1].Gamma_L = -2*imag(Sigma_L)Sigma_L = lead_L.selfenergy(energy)lead_L = sys.leads[0]sys = sys.finalized()Dear Shivang,To get the matrix Gamma you can do:
To get directly the transmission matrix t you can do:
AdelOn Wed, Jun 6, 2018 at 11:53 PM, Shivang Agarwal <shivang.agarwal@iitgn.ac.in> wrote: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 conductor (system) and '†' is the dagger operator. The equation is from the following paper: https://journals.aps.org/prb/pdf/10.1103/PhysRevB. 73.075429 (1) Now as far as I know, Kwant allows us to calculate transmission as a number T(E). What I need for my code is 't' where Trace(t†t) = T(E). Could somebody let me know how can I get the desired quantity 't'?. But I don't know how I can get the coupling matrix ГL(E) between the left (or right) lead and the conductor(2) Also, we know that t = ГL(E)½ GC(E) ГR(E)½ .But I don't know how I can get the coupling matrix ГL(E) between the left (or right) lead and the conductor. Is it possible to get too?PS - My aim is to find the wavefunctions inside the nanoribbon (which Kwant can do very conveniently) and also their phases! I have found the wavefunctions but am unable to find their phases. If there's any other way to find it that would also be extremely helpful.Any help would be greatly appreciated.Thanks and Regards,Shivang Agarwal--Shivang AgarwalJunior UndergraduateDiscipline of Electrical EngineeringIIT GandhinagarContact: +91-9869321451
--Abbout Adel--Shivang AgarwalJunior UndergraduateDiscipline of Electrical EngineeringIIT GandhinagarContact: +91-9869321451--Shivang AgarwalJunior UndergraduateDiscipline of Electrical EngineeringIIT GandhinagarContact: +91-9869321451