**kwargs_IV_Response_John)
##############################################
##### Introduce a factor to the IV curve #####
##############################################
#factor = np.arange(.7,1.5,.05) Does not work since it gets weird small values
factor = np.divide(np.arange(7, 15, .5), 10)
IVbyFactor = dict()
for i in factor:
IVbyFactor[i] = np.vstack(
[IV.binedIVData[0],
np.multiply(i, IV.binedIVData[1])])
title = newfig('Sum_Factor_1_1')
plot(IVbyFactor[1.], label='Factor 1.')
plot(IVbyFactor[1.], label='Factor 1.')
plt.plot(IVbyFactor[1.][0], IVbyFactor[1.][1] + IVbyFactor[1.][1], label='Sum')
pltsettings(save=directory + title,
fileformat='.pdf',
disp=True,
close=True,
xlabel=lbl['mV'],
ylabel=lbl['uA'],
xlim=[0, 2.7],
ylim=[0, 30],
title=None,
legendColumns=1,
skip_legend=False)
'c-6.00-45-iv.csv': ([1.1299, 5.1355]),
'c-5.60-45-iv.csv': ([0., 0.]),
'c-6.30-45-iv.csv': ([0., 0.]),
'c-5.80-45-iv.csv': ([0., 0.]),
'c-5.00-45-iv.csv': ([-1.2346, .45755]),
'c-6.10-45-iv.csv': ([-1.12578, 5.355]),
'c-5.20-45-iv.csv': ([0., 0.])
}
for i in filenames['univ']:
kwargs_IV_Response_rawData['fixedOffset'] = offset[i]
Unpumped = IV_Response(i, **kwargs_IV_Response_rawData)
#find the index of the corresponding pumped iv curve
index = np.where(filenamesstr['iv'][:, 1] == i[2:6])[0][0]
Pumped = IV_Response(filenames['iv'][index], **kwargs_IV_Response_rawData)
title = newfig('Unpumped_Pumped_at_%s_K' % (i[2:6]))
plot(Unpumped.binedIVData, label='Unpumped ' + i)
plot(Pumped.binedIVData, label='Pumped ' + filenames['iv'][index])
pltsettings(save=directory + title,
fileformat='.pdf',
disp=True,
close=False,
xlabel=lbl['mV'],
ylabel=lbl['uA'],
title=None,
legendColumns=1,
skip_legend=False)
pltsettings(save=directory + title + 'zoom',
fileformat='.pdf',
disp=True,
close=False,
'D1_17_n.csv': [-200, 200]
}
r300K = {
'Backup_Old_m.csv': 90.5,
'Backup_Old_n.csv': 90.5,
'D1_15_m.csv': 70.4,
'D1_15_n.csv': 70.4,
'D1_17_m.csv': 85,
'D1_17_n.csv': 85
}
for i in filenames:
kwargs_IV_Response_rawData['normalResistance300K'] = r300K[i]
IV = IV_Response(i, **kwargs_IV_Response_rawData)
title = newfig(i.replace('.csv', '_Manual_Offset_Correction'))
IV.plot_IV_with_Info([-3.85, 10, -3.85, 10], linespacing=1.7, fontsize=10)
plt.tight_layout()
pltsettings(save=directory + title,
fileformat='.pdf',
disp=True,
close=True,
xlabel=lbl['mV'],
ylabel=lbl['uA'],
xlim=[-4, 4],
ylim=ylims[i],
title=None,
legendColumns=1,
skip_legend=True)
#Old implementation before IV_Class was updated on 14.02.2020.
parser.add_argument('-f',
'--folder',
action='store',
default='Default_Folder',
help='The folder in which the result is stored in.')
args = parser.parse_args()
directory = 'Subgap_Leakage_Current_Unit_Test/' + args.folder + '/'
if not os.path.exists(directory):
os.makedirs(directory)
tC = 9.2 # K for Nb
rN = 13
bias = np.arange(0, 3e-3, 0.001e-3)
title = newfig('Temperature_Dependence')
temperatures = [3, 4, 5, 6, 7]
for te in temperatures:
delta = cooperPair_Binding_Energy_over_T(
t_over_Tc=te / tC, delta0=cooperPair_Binding_Energy_0K(tC),
tDebye=276)[0]
current = subgapLeakageCurrent(v=bias, t=te, d=delta, rN=rN)
plt.plot(bias * 1e3, current * 1e6, label='%.1f K' % te)
pltsettings(save=directory + title,
fileformat='.pdf',
disp=True,
close=False,
xlabel=lbl['mV'],
ylabel=lbl['uA'],
xlim=None,
ylim=None,
sys.path.insert(0, '../Helper')
from Cooper_Pair_Tunnelling_Current import cooperPairTunnellingCurrent,magneticFluxQuantum
from Fundamental_BCS_Equations import cooperPair_Binding_Energy_over_T,cooperPair_Binding_Energy_0K
from plotxy import newfig,pltsettings,lbl,plot
parser = ArgumentParser()
parser.add_argument('-f', '--folder', action='store',default = 'Default_Folder', help='The folder in which the result is stored in.')
args = parser.parse_args()
directory = 'Cooper_Pair_Tunnelling_Current_Unit_Test/'+args.folder+'/'
if not os.path.exists(directory):
os.makedirs(directory)
title = newfig('Flux_Dependency')
fluxes = np.arange(0,15,.001) * magneticFluxQuantum
te = 4 #K, The actual temperature
tC = 9.2 # K for Nb
delta = cooperPair_Binding_Energy_over_T(t_over_Tc =te/tC,delta0 = cooperPair_Binding_Energy_0K(tC), tDebye = 276)[0]
rN =10
current = cooperPairTunnellingCurrent(fluxes,delta,te,rN)
plt.plot(fluxes,current)
pltsettings(save=directory+title,fileformat='.pdf',disp = True,close=True, xlabel=lbl['Wb'],ylabel=lbl['A'],
xlim=None,ylim=None,title=None,legendColumns=1,skip_legend=True)
title = newfig('Normal_Resistance_Dependency')
fluxes = 0
te = 4 #K, The actual temperature
delta = cooperPair_Binding_Energy_over_T(t_over_Tc =te/tC,delta0 = cooperPair_Binding_Energy_0K(tC), tDebye = 276)[0] #9.2 K for Nb
rN =np.arange(5,40,.01)