Hi, part of links are broken. I also want to compute separately conductance 
up-up, down-down, up-down and down-up.
I read the threads where up, but I still not sure if i implemented in code 
correctly with these informations.

In my code i use rashba interaction (in code):
----------------------------------------------------------------------------------
def rashba(xi, yi, xj, yj, alpha, t):
    #d_ij = sitej.pos - sitei.pos
    x_ij = xj - xi
    y_ij = yj - yi
    rashb = 1j * alpha * (sigma_x * y_ij - sigma_y * x_ij)
    return (rashb - t * sigma_0)

-----------------------------------------------------------------------------------
i use norbs=1, because is only electrons
And in conductances computations i've implemented (in code):
-----------------------------------------------------------------------------------

    smat = kwant.smatrix(syst, energy)

    t = np.matrix(smat.submatrix(1, 0)) #trasmission block from S-matrix

    tt_mat = np.matmul(t, t.getH()) #  t*t_dagger matrix

    tt = np.trace(np.matmul(t, t.getH())) #trace of t*t_dagger = total 
transmission counts (up/down)
    
    #getting upup Transmission, downdown and downup
    #up-up t*t_dagger block
    t11 = np.trace(np.matrix([[tt_mat[i,j] for j in range(int(len(tt_mat)/2))] 
for i in range(int(len(tt_mat)/2))])) 
    #down-down t*t_dagger block
    t22 = np.trace(np.matrix([[tt_mat[i + int(len(tt_mat)/2), j + 
int(len(tt_mat)/2)] for j in range(int(len(tt_mat)/2))] for i in 
range(int(len(tt_mat)/2))])) 
    #down-up t*t_dagger block
    t21 = np.trace(np.matrix([[tt_mat[i + int(len(tt_mat)/2), j] for j in 
range(int(len(tt_mat)/2))] for i in range(int(len(tt_mat)/2))]))

-----------------------------------------------------------------------------------

It is completely incorrect compute using these form? 
I've tryed by smat.transmission((1,0), (0,0)) structure but i got error, and 
decide to try these aproach above.
by information the i got using smat.transmission((1,0), (0,0)) structure:
-----------------------------------------------------------------------------------
return self.out_block_coords(lead_out), self.in_block_coords(lead_in)
  File 
"C:\Users\Diego\AppData\Roaming\Python\Python39\site-packages\kwant\solvers\common.py",
 line 854, in in_block_coords
    self.in_block_offsets[lead_ind][block_ind + 1])
IndexError: index 2 is out of bounds for axis 0 with size 2
-----------------------------------------------------------------------------------

Thanks in advance!

Reply via email to