Gas di Fermi

In fisica, in particolare in meccanica statistica, un gas di Fermi è un gas di fermioni. La statistica di Fermi-Dirac permette di determinare la distribuzione dell'energia per un gas di fermioni all'equilibrio termico conoscendone la densità, la temperatura e il set di stati energetici possibili.

Con questo modello possono essere in prima approssimazione descritti i nucleoni all'interno del nucleo atomico o gli elettroni di conduzione in un metallo.

La statistica di Fermi-Dirac

Un gas di Fermi composto da particelle identiche segue la statistica di Fermi-Dirac, dalla quale si deduce che:

${\overline {n_{k}}}={\frac {1}{e^{(\varepsilon _{k}-\mu )/kT}+1}}$

che rappresentano i valori medi dei numeri di occupazione per un gas di Fermi. Per un gas di Fermi tutti i numeri di occupazione sono ${\overline {n_{k}}}\leq 1$. La normalizzazione impone:

$N=\sum _{k}{\frac {1}{e^{(\varepsilon _{k}-\mu )/kT}+1}}$

dove N è il numero totale di particelle del gas. Da questa possiamo determinare il potenziale chimico.

L'hamiltoniana di un gas di Fermi costituito da N fermioni di massa m racchiuso all'interno di una scatola cubica di lato L è:

$H_{0}=\sum _{i=1}^{N}{\frac {p_{i}^{2}}{2m}}$

dove l'energia di singola particella è:

$\varepsilon ={\frac {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}{2m}}$

espressi in termini di autovalori, ovvero i valori di energia accessibili al sistema: $\varepsilon _{k}={\frac {\hbar ^{2}k^{2}}{2m}}$. Tenendo conto della degenerazione di spin $g=2s+1$ dove s è lo spin del fermione, il numero di particelle nell'elemento di volume dello spazio delle fasi è:

$gd\tau =g{\frac {dp_{x}dp_{y}dp_{z}dV}{(2\pi \hbar )^{3}}}$

cioè si ha per la distribuzione di Fermi:

$dN={\frac {gd\tau }{e^{(\varepsilon -\mu )/kT}+1}}$

Più precisamente, integrando quest'ultima in $dV$ si ottiene la distribuzione dell'impulso:

$dN_{p}={\frac {gVp^{2}dp}{2\pi ^{2}\hbar ^{3}\left[e^{(\varepsilon -\mu )/kT}+1\right]}}$

e siccome $\varepsilon =p^{2}/2m$, si deduce subito la distribuzione di energia:

$dN_{\varepsilon }={\frac {gVm^{3/2}}{{\sqrt {2}}\pi ^{2}\hbar ^{3}}}{\frac {{\sqrt {\varepsilon }}\,d\varepsilon }{e^{(\varepsilon -\mu )/kT}+1}}$

Le ultime due espressioni rappresentano le distribuzioni di Maxwell nel caso di un gas di Fermi. Dalla seconda si ricava l'energia:

$\int _{0}^{\infty }\varepsilon \,dN_{\varepsilon }={\frac {gVm^{3/2}}{{\sqrt {2}}\pi ^{2}\hbar ^{3}}}\int _{0}^{\infty }{\frac {\varepsilon ^{3/2}\,d\varepsilon }{e^{(\varepsilon -\mu )/kT}+1}}$

e il numero totale di particelle:

$N={\frac {gVm^{3/2}}{{\sqrt {2}}\pi ^{2}\hbar ^{3}}}\int _{0}^{\infty }{\frac {{\sqrt {\varepsilon }}\,d\varepsilon }{e^{(\varepsilon -\mu )/kT}+1}}$

che fornisce anche il potenziale termodinamico:

$\Omega =-{\frac {VgT{\sqrt {(m)^{3}}}}{{\sqrt {2}}\pi ^{2}\hbar ^{3}}}\int _{0}^{\infty }{\frac {\varepsilon ^{3/2}\,d\varepsilon }{e^{(\varepsilon -\mu )/kT}+1}}$

Esso coincide con l'energia a meno di un fattore:

$\Omega =-PV={\frac {2}{3}}E$

che è una relazione del tutto generale e vale per tutti i sistemi, sia di Bose, sia di Fermi e di Boltzmann.

Gas di Fermi completamente degenere

Supponiamo di avere un gas di fermioni di spin $s=1/2$ (quindi $g=2s+1=2$), ad esempio elettroni, a temperatura assoluta $T$ = 0 K. Gli elettroni a tale temperatura cercano di porsi negli stati a energia minore in modo che il totale dell'energia sia il valore più basso possibile, partendo dallo stato a energia nulla, fino a un certo valore.

Il numero di stati quantistici di un elettrone in un volume V con impulso compreso in $(p,p+dp)$ è dato dalla prima delle distribuzioni di Maxwell:

$2{\frac {4\pi Vp^{2}\,dp}{(2\pi \hbar )^{3}}}=V{\frac {p^{2}\,dp}{\pi ^{2}\hbar ^{3}}}$

Gli elettroni occupano tutti gli stati con impulso uguale a zero (notare che $\varepsilon =p^{2}/2m$) fino al valore $p=p_{F}$ detto impulso di Fermi, che equivale nello spazio degli impulsi, al raggio di una sfera detta sfera di Fermi.

Il numero totale di elettroni in questi stati è dato da:

$N={\frac {V}{\pi ^{2}\hbar ^{3}}}\int _{0}^{p_{F}}p^{2}\,dp={\frac {Vp_{F}^{3}}{3\pi ^{2}\hbar ^{3}}}$

da cui possiamo ricavare l'impulso di Fermi:

$p_{F}=(3\pi ^{2})^{1/3}\hbar \left({\frac {N}{V}}\right)^{1/3}$

e l'energia di Fermi:

$\varepsilon _{F}={\frac {p_{F}^{2}}{2m}}=(3\pi ^{2})^{2/3}{\frac {\hbar ^{2}}{2m}}\left({\frac {N}{V}}\right)^{2/3}$

Questo si può vedere anche dai numeri di occupazione medi: infatti al limite $T\to 0$:

$\lim _{T\to 0}{\overline {n_{\mathbf {p} }}}=\lim _{T\to 0}{\frac {1}{e^{(\varepsilon -\mu )/kT}+1}}=\left\{{\begin{matrix}1&\varepsilon <\mu \\0&\varepsilon >\mu \end{matrix}}\right.$

cioè i numeri di occupazione medi diventano una funzione a gradino facendoci pensare al fatto che per $p<p_{F}$ oppure $\varepsilon <\varepsilon _{F}$ gli elettroni si dispongono a partire dal livello $\varepsilon =0$ fino ai livelli $p=p_{F}$ oppure $\varepsilon =\varepsilon _{F}$ con la condizione che in un livello ci sia al massimo una particella secondo il principio di esclusione di Pauli. Dopo questi valori, per $p>p_{F}$ non vi sono più elettroni da sistemare.

Da notare che:

$\varepsilon _{F}=\mu \$

L'energia totale del gas di Fermi completamente degenere si ottiene dall'integrazione:

$E={\frac {V}{2m\pi ^{2}\hbar ^{3}}}\int _{0}^{p_{F}}p^{4}\,dp={\frac {Vp_{F}^{5}}{10m\pi ^{2}\hbar ^{3}}}$

che, sostituendo l'espressione dell'impulso di Fermi e, nel successivo passaggio, quella dell'energia di Fermi, diventa:

$E={\frac {3(3\pi ^{2})^{2/3}\hbar ^{2}}{10m}}\left({\frac {N}{V}}\right)^{2/3}N={\frac {3}{5}}{N\varepsilon _{F}}$

Infine, usando la relazione generale del potenziale termodinamico, si ottiene:

$P={\frac {(3\pi ^{2})^{2/3}\hbar ^{2}}{5m}}\left({\frac {N}{V}}\right)^{5/3}$

ovvero: la pressione è proporzionale alla densità secondo la potenza 5/3.

Voci correlate

• Statistica di Fermi-Dirac
• Energia di Fermi
• Gas di Bose
• Gas di fotoni
• Gas ideale quantistico
• Materia degenere
• Liquido di Fermi