Commit 727f7a99 authored by Isidora Araya's avatar Isidora Araya
Browse files

Update 13_semiconductors_solutions.md

parent 1ec2fa55
......@@ -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
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment