Can you please provide me the script for the same curves you've shown in the fig.?
The example is as follows:
#---------------------------------------------------
bias_voltage = numpy.linspace(0.0,1.0,6)
current = []
for vb in bias_voltage:
transmission_spectrum = nlread('gnr_IV.nc', object_id="Transmission"+str(vb))[0]
current += [transmission_spectrum.current()]
print "-"*50
print "Bias Voltage [V]","\t","Current [nA]"
print "-"*50
for vb,i in zip(bias_voltage,current):
print vb,"\t",i.inUnitsOf(Ampere)*(10**9)
#---------------------------------------------------