# 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)