Hi Patrik,
> 1- The wave function approach does not give the the up/down > conductance similar to total. I see that you do the following: > psi_start = psi[0 : 2] > psi_end = psi[2 * 61: 2 * 61 + 2] Do you know if the first and last sites sites as they are ordered in the finalized system are indeed the sites that you want? In your script you never refer to the 'sites' attribute of the finalized system, so I would guess that these are not the sites that you want; sites number 0 and 61 are probably not even connected. > 2- the current-operator approach gives the error > 'kwant.graph.core.EdgeDoesNotExistError'. This means that you are trying to calculate the current between 2 sites that are not connected by a hopping. > J_spin = kwant.operator.Current(syst, sigma_z, > where=[(amorphous_lat(0), amorphous_lat(N-2))], sum=True) I see that you are trying to calculate the current between 'amorphous_lat(0)' and 'amorphous_lat(N-2)', and if I colour these sites when plotting the system I see indeed that they are at opposite ends of your double helix structure, and that they are not connected by a hopping. Happy Kwanting, Joe
signature.asc
Description: OpenPGP digital signature
