forked from qpv-research-group/solcore5
-
Notifications
You must be signed in to change notification settings - Fork 0
/
coherency_example.py
79 lines (68 loc) · 2.14 KB
/
coherency_example.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
import numpy as np
import matplotlib.pyplot as plt
from solcore import material, si
from solcore.solar_cell import SolarCell, Layer, Junction
from solcore.solar_cell_solver import solar_cell_solver, prepare_solar_cell
from solcore.state import State
GaInP = material("GaInP")(In=0.5)
GaAs = material("GaAs")()
Ge = material("Ge")()
optical_struct = SolarCell(
[
Layer(material=GaInP, width=si("5000nm")),
Junction(
[
Layer(material=GaAs, width=si("200nm")),
Layer(material=GaAs, width=si("5um")),
],
kind="DA",
),
Layer(material=Ge, width=si("50um")),
]
)
wl = np.linspace(250, 1700, 300) * 1e-9
options = State()
options.wavelength = wl
options.optics_method = "TMM"
options.no_back_reflection = False
options.BL_correction = True
options.recalculate_absorption = True
options.positions = [1e-8, 1e-9, 1e-8, 1e-7]
options.theta = 0
c_list = [
["c", "c", "c", "c"],
["c", "c", "c", "i"],
["c", "i", "i", "c"],
["i", "i", "i", "i"],
]
titles = [
"All coherent",
"Bottom Ge layer explicity incoherent",
"Both layers of GaAs junction incoherent",
"All layers incoherent",
]
for i1, cl in enumerate(c_list):
plt.figure(i1)
options.coherency_list = cl
solar_cell_solver(optical_struct, "optics", options)
plt.plot(wl * 1e9, optical_struct[0].layer_absorption)
plt.plot(wl * 1e9, optical_struct[1].layer_absorption)
plt.plot(wl * 1e9, optical_struct[2].layer_absorption)
plt.plot(wl * 1e9, optical_struct.reflected, "--")
plt.plot(wl * 1e9, optical_struct.transmitted, "--")
plt.plot(
wl * 1e9,
optical_struct[0].layer_absorption
+ optical_struct[1].layer_absorption
+ optical_struct[2].layer_absorption
+ optical_struct.reflected
+ optical_struct.transmitted,
)
plt.legend(["GaInP", "GaAs", "Ge", "R", "T", "R+A+T"], loc="upper left")
plt.title(titles[i1])
plt.xlabel("Wavelength (nm)")
plt.show()
# plt.figure(1)
# plt.plot(np.linspace(0, 3000, 100),
# optical_struct[1].absorbed(np.linspace(0, 3000, 100)*1e-9)[:,0])
# plt.show()