Dear Anant, See the previous thread in this mailing list: http://thread.gmane.org/gmane.comp.science.kwant.user/276 ; it was about exactly the same question. Best, Anton On Wed, Mar 18, 2015 at 10:54 AM, ANANT VIJAY VARMA <avterminator@gmail.com> wrote:

Hello everyone! I am calculating the longitudinal resistance (SdHo oscillation) using a square lattice.The order of leads is as follows:

#Left Lead sys.attach_lead(leadA)# lead no. 0 #right lead sys.attach_lead(leadA.reversed())# lead no. 1 # First upper lead sys.attach_lead(leadB)# lead no. 2 # First lower lead sys.attach_lead(leadB.reversed())# lead no. 3 # Second upper lead sys.attach_lead(leadC) # lead no. 4 #second lower lead sys.attach_lead(leadC.reversed()) # lead no. 5

And using the Folowing piece of code to calculate the Longitudinal resistance:

energy =2 def plot_conductance(sys, energy, Bs): data = []

for B in Bs: smatrix = kwant.smatrix(sys, energy, args=[B]) Trans = np.array([[smatrix.transmission(i, j) for j in xrange(6)] for i in xrange(6)]) Trans -= np.diag(Trans.sum(axis=0)) cm = Trans[:-1, :-1]

resistance = np.linalg.solve(cm, [-1,1,0, 0, 0])[3] data.append(resistance)

My question is that it is correct also what is the unit of current being used here, and units of rest of the quantities?Also How I am supposed to write code gate voltage. I am trying the same way accept now my arguement is not B but G (gate voltage). But its wrong I know.Please help me about these questions. Thank you in advance.