# Layer 0 | 250.0 | Si | 0 | 1e16 | n |
# Layer 1 | 250.0 | Si | 0 | 1e16 | p |
#
# To input this list in Gallium, we use lists as:
material = [[200.0, 'AlGaAs', 0.2, 0.2, 0, 'n', 'b'],
[10.0, 'GaAs', 0, 0.2, 2e18, 'n', 'w'],
[200.0, 'AlGaAs', 0.2, 0.2, 0, 'n', 'b']]
#----------------------------------------
import numpy as np
x_max = sum([layer[0] for layer in material])
def round2int(x):
return int(x + 0.5)
n_max = round2int(x_max / gridfactor)
#----------------------------------------
dop_profile = np.zeros(n_max)
#----------------------------------------
Quantum_Regions = False
Quantum_Regions_boundary = np.zeros((2, 2))
#----------------------------------------
surface = np.zeros(2)
#surface[0]=-0.6
#----------------------------------------
if __name__ == "__main__": #this code allows you to run the input file directly
input_obj = vars()
import aestimo_eh
aestimo_eh.run_aestimo(input_obj)
[2.0, 'AlGaInN', 0.48, 0.23, 1e17, 'p', 'b'],
[2.0, 'InGaN', 0.25, 0.0, 0.0, 'i', 'w'],
[3.0, 'AlGaInN', 0.48, 0.23, 0.0, 'i', 'b'],
[2.0, 'InGaN', 0.25, 0.0, 0.0, 'i', 'w'],
[2.0, 'AlGaInN', 0.48, 0.23, 1e17, 'n', 'b'],
[200.0, 'AlGaN', 0.1, 0.0, 1e17, 'n', 'b']]
#This is accourding to interpolated Vegard’s law for quaternary ABxCyD1-x-y=NGaxAlyIn1-x-y
#----------------------------------------
import numpy as np
x_max = sum([layer[0] for layer in material])
def round2int(x):
return int(x + 0.5)
n_max = round2int(x_max / gridfactor)
#----------------------------------------
dop_profile = np.zeros(n_max)
#----------------------------------------
Quantum_Regions = False
Quantum_Regions_boundary = np.zeros((2, 2))
#----------------------------------------
surface = np.zeros(2)
#----------------------------------------
if __name__ == "__main__": #this code allows you to run the input file directly
input_obj = vars()
import aestimo_eh as aestimo
aestimo.run_aestimo(input_obj)
