# Interact Exercise 6

## Imports

Put the standard imports for Matplotlib, Numpy and the IPython widgets in the following cell.


```python
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
```


```python
from IPython.display import Image
from IPython.html.widgets import interact, interactive, fixed
```

## Exploring the Fermi distribution

In quantum statistics, the [Fermi-Dirac](http://en.wikipedia.org/wiki/Fermi%E2%80%93Dirac_statistics) distribution is related to the probability that a particle will be in a quantum state with energy $\epsilon$. The equation for the distribution $F(\epsilon)$ is:


```python
Image('fermidist.png')
```

In this equation:

* $\epsilon$ is the single particle energy.
* $\mu$ is the chemical potential, which is related to the total number of particles.
* $k$ is the Boltzmann constant.
* $T$ is the temperature in Kelvin.

In the cell below, typeset this equation using LaTeX:

\begin{equation}
F(\epsilon) = {\displaystyle \frac{1}{e^{(\epsilon - \mu)/kT}+1}}
\end{equation}

Define a function `fermidist(energy, mu, kT)` that computes the distribution function for a given value of `energy`, chemical potential `mu` and temperature `kT`. Note here, `kT` is a single variable with units of energy. Make sure your function works with an array and don't use any `for` or `while` loops in your code.


```python
def fermidist(energy, mu, kT):
    """Compute the Fermi distribution at energy, mu and kT."""
    return (np.exp((energy-mu)/kT)+1)**-1
```


```python
assert np.allclose(fermidist(0.5, 1.0, 10.0), 0.51249739648421033)
assert np.allclose(fermidist(np.linspace(0.0,1.0,10), 1.0, 10.0),
    np.array([ 0.52497919,  0.5222076 ,  0.51943465,  0.5166605 ,  0.51388532,
               0.51110928,  0.50833256,  0.50555533,  0.50277775,  0.5       ]))
```

Write a function `plot_fermidist(mu, kT)` that plots the Fermi distribution $F(\epsilon)$ as a function of $\epsilon$ as a line plot for the parameters `mu` and `kT`.

* Use enegies over the range $[0,10.0]$ and a suitable number of points.
* Choose an appropriate x and y limit for your visualization.
* Label your x and y axis and the overall visualization.
* Customize your plot in 3 other ways to make it effective and beautiful.


```python
np.arange(0,10.01,0.01)
```




    array([  0.  ,   0.01,   0.02, ...,   9.98,   9.99,  10.  ])




```python
def plot_fermidist(mu, kT):
    energy=np.arange(0,10.01,0.01)
    plt.figure(figsize=(10,6))
    plt.plot(energy,fermidist(energy,mu,kT))
    plt.tick_params(axis='x', top='off')
    plt.tick_params(axis='y', right='off')
    plt.xlabel('Energy')
    plt.xlim(left=0, right=10)
    plt.ylim(bottom=0.0,top=1.0)
    plt.ylabel('Fermi Distribution')
```


```python
plot_fermidist(4.0, 1.0)
```


```python
assert True # leave this for grading the plot_fermidist function
```

Use `interact` with `plot_fermidist` to explore the distribution:

* For `mu` use a floating point slider over the range $[0.0,5.0]$.
* for `kT` use a floating point slider over the range $[0.1,10.0]$.


```python
interact(plot_fermidist, mu=(0.0,5.0,0.1), kT=(0.1,10.0,0.1));
```

Provide complete sentence answers to the following questions in the cell below:

* What happens when the temperature $kT$ is low?
* What happens when the temperature $kT$ is high?
* What is the effect of changing the chemical potential $\mu$?
* The number of particles in the system are related to the area under this curve. How does the chemical potential affect the number of particles.

Use LaTeX to typeset any mathematical symbols in your answer.

$\bullet$ When $kT$ is low, the Fermi distribution starts off very high, around $1.0$, but then sharply drops off to zero, or near zero, as energy increases.

$\bullet$ When $kT$ is high, the Fermi distribution becomes more linear in its shape. 

$\bullet$ Changing the chemical potential of the equation "shifts" the graph left or right as $\mu$ decreases or increases, respectively.

$\bullet$ As the chemical potential increases from $0$ to $5$, the area under the curve inceases as well. This means the particles in the system are being distributed more throughout the energies of the system. However, that is if we hold $kT$ at a constant value. If we increase both $kT$ and $\mu$ together, the area under the curve will change constantly, thus the particles are even more distributed throughout the energies.

dc51994fe7435e7532f5e1bdc3e0923bada3a5a0
ipynb
Jupyter Notebook
midterm/InteractEx06.ipynb
brettavedisian/phys202-2015-work
baf9f3ca507873547cf142b87d56bf879031121f
MIT
midterm/InteractEx06.ipynb
brettavedisian/phys202-2015-work
baf9f3ca507873547cf142b87d56bf879031121f
MIT
midterm/InteractEx06.ipynb
brettavedisian/phys202-2015-work
baf9f3ca507873547cf142b87d56bf879031121f
MIT
Qwen/Qwen-72B
1. YES
2. YES