def equation1(vb, gammaL, gammaR, deltaE1, eta, sigma):
T = 300
c = 0
vg = 0
gammaC = gammaL * gammaR
gammaW = gammaL + gammaR
return models.tunnelmodel_singleLevel(vb, gammaC, gammaW, deltaE1, eta,
sigma, c, vg, T)
def fxn(vb, gammaL, gammaR, deltaE, eta, sigma):
c = 0
vg = 0
T = Temp[col]
n = 1000
gammaC = gammaL * gammaR
gammaW = gammaL + gammaR
I = models.tunnelmodel_singleLevel(vb, gammaC, gammaW, deltaE, eta,
sigma, c, vg, T)
return n * I
def memoryMolecule_gauss(vb, pos, width, p0, gammaC1_H, gammaW1_H,
deltaE1_H, eta1_H, sigma1_H, gammaC2_H, gammaW2_H,
deltaE2_H, eta2_H, sigma2_H, gammaL_L, gammaR_L,
deltaE_L, eta_L, width_L):
c = 0
vg = 0
T = 300
args_L = (gammaL_L, gammaR_L, deltaE_L, eta_L, width_L, c, vg, T)
P_H = models.sigmoid(vb, pos, width) + (vb - 1) * p0
P_L = 1 - P_H
I_H = models.tunnelmodel_2level(vb, c, vg, T, gammaC1_H, gammaW1_H,
deltaE1_H, eta1_H, sigma1_H, gammaC2_H,
gammaW2_H, deltaE2_H, eta2_H, sigma2_H)
I_L = models.tunnelmodel_singleLevel(vb, *args_L)
return P_H * I_H + P_L * I_L