def print_charge():
print(
devsim.get_contact_charge(device=device,
contact="top",
equation="PotentialEquation"))
print(
devsim.get_contact_charge(device=device,
contact="bot",
equation="PotentialEquation"))
def simulate_charge(device, contact, equation, solver_params):
#charge_factor=1e7 #from F/cm^2 to fF/um^2
dv = 0.001
v1 = ds.get_parameter(device=device, name=GetContactBiasName(contact))
q1 = ds.get_contact_charge(device=device,
contact=contact,
equation="PotentialEquation")
v2 = v1 + dv
ds.set_parameter(name=GetContactBiasName(contact), value=v2)
ds.solve(**solver_params)
q2 = ds.get_contact_charge(device=device,
contact=contact,
equation="PotentialEquation")
return (v1, (charge_factor * (q2 - q1) / dv))
