From 727f7a99a870ef467bc1679754c7f249625f278b Mon Sep 17 00:00:00 2001
From: Isidora Melania Araya Day <isidora.araya@ug.uchile.cl>
Date: Mon, 30 Mar 2020 02:00:22 +0000
Subject: [PATCH] Update 13_semiconductors_solutions.md
---
src/13_semiconductors_solutions.md | 14 +++++++++-----
1 file changed, 9 insertions(+), 5 deletions(-)
diff --git a/src/13_semiconductors_solutions.md b/src/13_semiconductors_solutions.md
index 6ee43514..000e6716 100644
--- a/src/13_semiconductors_solutions.md
+++ b/src/13_semiconductors_solutions.md
@@ -59,21 +59,25 @@ $$ n_h = \int_{\varepsilon-2t}^{\varepsilon+2t} (1-f(\varepsilon)) g_h(\varepsil
## Exercise 3:
```python
-def dispersion(N=100, kmax=pi/2, EG, tcb, tvb):
-
+def dispersion(EG, tcb, tvb, N=100, kmax=np.pi/2):
+ a = 1
kx = np.linspace(-kmax, kmax, N)
Ecb = EG - 2*tcb*(np.cos(kx*a)-1)
Evb = 2*tvb*(np.cos(kx*a)-1)
# Plot dispersion
plt.figure(figsize=(6,5))
- plt.plot(kx, Ecb, label="Conduction B.")
- plt.plot(kx, Evb, label="Valence B.")
+ cb, = plt.plot(kx, Ecb, label="Conduction B.")
+ vb, = plt.plot(kx, Evb, label="Valence B.")
plt.xlabel('$k_x$', fontsize=20)
plt.ylabel('$E$', fontsize=20)
+ plt.title('E(k) for tcb:'+str(tcb)+' tvb:'+str(tvb))
+ plt.legend(handles=[cb, vb])
+
+ plt.show()
-dispersion(10, 2, 2)
+dispersion(10, 2, 8)
```
### Subquestion 1
--
GitLab