I have been trying to reproduce the results of the following paper: " Negative differential resistance in nanoscale transport in the Coulomb blockade regime - Prakash Parida et al 2009 J. Phys.: Condens. Matter 21 095301 ". A code is being developed for this, and all the parameters are taken to be same as in the paper.
With all the considerations according to the paper, the final current is coming out to be incorrect; there is some leakage current at zero source-drain bias (i.e. at Vsd=0). With further analysis of the issue, it appears that transitions are occuring between 4e states and 3e states, which is incorrect.
I have not found any issues with implementation of the theory in the code. Any suggestions in this regard would be extremely helpful.