Gas ideale quantistico

In meccanica statistica il gas ideale quantistico è un gas ideale tale che le particelle che lo compongono sono indistinguibili, e vanno trattate nell'ambito della meccanica quantistica.

In natura esistono due tipi fondamentali di particelle identiche: i fermioni e i bosoni: un insieme di fermioni è descritto da un autostato antisimmetrico rispetto allo scambio di due particelle, mentre per i bosoni l'autostato è simmetrico; riguardo alla statistica, per i fermioni si usa la statistica di Fermi-Dirac, mentre per i bosoni la statistica di Bose-Einstein. La natura di questa differenza risiede in particolare nella proprietà dello spin delle particelle: secondo il teorema spin-statistica, particelle con spin intero sono bosoni e particelle con spin semiintero sono fermioni.

Gas ideali

In meccanica statistica quantistica parlando di un gas ideale quantistico intendiamo un gas le cui interazioni tra particelle possono essere trascurate. L'hamiltoniana di un gas ideale non interagente composto da N particelle identiche in un volume L 3 = V {\displaystyle L^{3}=V} è:

H = i = 1 N p i 2 2 m {\displaystyle {\mathcal {H}}=\sum _{i=1}^{N}{\frac {p_{i}^{2}}{2m}}}

dove p i {\displaystyle p_{i}} è l'impulso della particella i-esima. La meccanica statistica vuole che a partire dall'hamiltoniano H {\displaystyle {\mathcal {H}}} possiamo ricavare gli autovalori ε p {\displaystyle \varepsilon _{\mathbf {p} }} e le autofunzioni ϕ {\displaystyle \phi } per ciascuna particella, risolvendo l'equazione agli autovalori:

H ϕ n = ε p ϕ n {\displaystyle {\mathcal {H}}\phi _{n}=\varepsilon _{\mathbf {p} }\phi _{n}}

dove identifichiamo gli autovalori con le energie ε p {\displaystyle \varepsilon _{\mathbf {p} }} dipendenti solo dall'impulso poiché l'energia è solo cinetica. Ogni autovalore è quindi:

ε p = p 2 2 m {\displaystyle \varepsilon _{\mathbf {p} }={\frac {p^{2}}{2m}}}

dove:

p = 2 π L n {\displaystyle \mathbf {p} ={\frac {2\pi \hbar }{L}}\mathbf {n} }

n {\displaystyle \mathbf {n} } è un numero intero. Questi livelli energetici per un sistema macroscopico sono continui.

Lo stato di un sistema ideale è completamente individuato dai suoi numeri di occupazione. Questi numeri significano che nello stato in questione vi sono n p {\displaystyle n_{\mathbf {p} }} particelle con impulso p {\displaystyle \mathbf {p} } o equivalentemente che ci sono n ε {\displaystyle n_{\varepsilon }} particelle di energia ε {\displaystyle \varepsilon } . Ovviamente l'energia totale del sistema di N particelle è:

E = p ε p n p {\displaystyle E=\sum _{\mathbf {p} }\varepsilon _{\mathbf {p} }\cdot n_{\mathbf {p} }}
N = p n p {\displaystyle N=\sum _{\mathbf {p} }n_{\mathbf {p} }}

Ebbene per un sistema di bosoni n p = 0 , 1 , 2 , . . . {\displaystyle n_{p}=0,1,2,...} infatti non vi è alcuna restrizione al numero di particelle che possono avere uno stato descritto da un particolare valore dell'impulso, mentre per i sistemi di fermioni n p = 0 , 1 {\displaystyle n_{p}=0,1} , infatti per tali particelle non vi possono essere più di una sola particella che occupa uno stesso stato con quel valore dell'impulso, questa differenza fondamentale è dovuta al principio di esclusione di Pauli e all'indistinguibilità delle particelle. Per un gas di Boltzmann possiamo assumere che n p = 0 , 1 , . . . {\displaystyle n_{p}=0,1,...} analogamente ai bosoni, le particelle che possono essere trattate con Boltzmann non subiscono alcuna restrizione sui numeri di occupazione.

Abbiamo detto che i livelli energetici di un gas sono quasi-continui. Supponiamo che il nostro gas abbia un'energia E assegnata. Lo spettro di energia forma un continuo di livelli e pensiamo di dividere questo continuo in celle i, ognuna delle quali ha un numero g i {\displaystyle g_{i}} di livelli energetici. Il numero di occupazione di ogni cella: n i {\displaystyle n_{i}} , è la somma di n p {\displaystyle n_{\mathbf {p} }} particelle di impulso p {\displaystyle \mathbf {p} } contenuto nella cella. Allora chiamiamo W { n i } {\displaystyle W\{n_{i}\}} il numero degli stati del sistema corrispondenti all'insieme dei numeri di occupazione { n i } {\displaystyle \{n_{i}\}} , così:

Γ ( E ) = n i W { n i } {\displaystyle \Gamma (E)=\sum _{n_{i}}W\{n_{i}\}}

cioè

W { n i } = Π j ω j   {\displaystyle W\{n_{i}\}=\Pi _{j}\omega _{j}\ }

dove ω j {\displaystyle \omega _{j}} è il numero dei modi in cui possiamo sistemare n i {\displaystyle n_{i}} particelle nella cella i-esima.

  • Per un gas di Bose ogni livello energetico può essere occupato da qualsiasi numero di particelle, ogni cella può contenere quanto vogliamo particelle quindi si hanno le combinazioni con ripetizione, data l'indistinguibilità delle particelle:
ω j = ( n j + g j 1 n j ) = ( n j + g j 1 ) ! n j ! ( g j 1 ) ! {\displaystyle \omega _{j}={n_{j}+g_{j}-1 \choose n_{j}}={\frac {(n_{j}+g_{j}-1)!}{n_{j}!(g_{j}-1)!}}}

così:

W { n j } = Π j ( n j + g j 1 ) ! n j ! ( g j 1 ) ! {\displaystyle W\{n_{j}\}=\Pi _{j}{\frac {(n_{j}+g_{j}-1)!}{n_{j}!(g_{j}-1)!}}}
  • Per un gas di Fermi il numero di particelle che possono occupare un livello è 0 (nessuna) oppure 1. Pertanto si hanno le combinazioni:
ω j = g j ! n j ! ( g j n j ) ! {\displaystyle \omega _{j}={\frac {g_{j}!}{n_{j}!(g_{j}-n_{j})!}}}
W { n j } = Π j g j ! n j ! ( g j n j ) ! {\displaystyle W\{n_{j}\}=\Pi _{j}{\frac {g_{j}!}{n_{j}!(g_{j}-n_{j})!}}}
  • Per un gas di Boltzmann invece scambiare impulsi tra le particelle porta ad un nuovo stato ma non cambia { n p } {\displaystyle \{n_{\mathbf {p} }\}} per cui dovrebbe essere:
N ! Π j g j n j n j ! {\displaystyle N!\Pi _{j}{\frac {g_{j}^{n_{j}}}{n_{j}!}}}

ma per il paradosso di Gibbs, dobbiamo tenere conto del corretto conteggio di Boltzmann e dividiamo questa quantità per N ! {\displaystyle N!} , per cui:

W { n j } = Π j g j n j n j ! {\displaystyle W\{n_{j}\}=\Pi _{j}{\frac {g_{j}^{n_{j}}}{n_{j}!}}}

Entropia

Sulla base del conteggio dei numeri di occupazione possiamo ricavare l'entropia:

S = k log Γ ( E )   {\displaystyle S=k\log \Gamma (E)\ }

se sommiamo W { n i } {\displaystyle W\{n_{i}\}} su tutti gli stati { n i } {\displaystyle \{n_{i}\}} . Per fare questo supponiamo che vada bene:

S = k log W { n i ¯ } {\displaystyle S=k\log W\{{\overline {n_{i}}}\}}

infatti Γ ( E ) {\displaystyle \Gamma (E)} è ben approssimato da W { n i ¯ } {\displaystyle W\{{\overline {n_{i}}}\}} , per dimostrarlo bisognerebbe massimizzare la funzione W { n i } {\displaystyle W\{n_{i}\}} tenendo conto dei vincoli E = ε i n i {\displaystyle E=\sum \varepsilon _{i}n_{i}} ed N = n i {\displaystyle N=\sum n_{i}} . Il calcolo mostra che per la statistica di Bose:

n i ¯ = g i z 1 e β ε i 1 {\displaystyle {\overline {n_{i}}}={\frac {g_{i}}{z^{-1}e^{\beta \varepsilon _{i}}-1}}}

per quella di Fermi:

n i ¯ = g i z 1 e β ε i + 1 {\displaystyle {\overline {n_{i}}}={\frac {g_{i}}{z^{-1}e^{\beta \varepsilon _{i}}+1}}}

per quella di Boltzmann:

n i ¯ = g i z e β ε i {\displaystyle {\overline {n_{i}}}=g_{i}ze^{-\beta \varepsilon _{i}}}

presi con i rispettivi segni, dove z = e β μ {\displaystyle z=e^{\beta \mu }} è la fugacità e β = 1 / k T {\displaystyle \beta =1/kT} è la beta termodinamica. Oppure, per un gas di Bose:

n p ¯ = 1 z 1 e β ε p 1 {\displaystyle {\overline {n_{\mathbf {p} }}}={\frac {1}{z^{-1}e^{\beta \varepsilon _{\mathbf {p} }}-1}}}

per un gas di Fermi:

n p ¯ = 1 z 1 e β ε p + 1 {\displaystyle {\overline {n_{\mathbf {p} }}}={\frac {1}{z^{-1}e^{\beta \varepsilon _{\mathbf {p} }}+1}}}

e per un gas di Boltzmann:

n p ¯ = z e β ε p {\displaystyle {\overline {n_{\mathbf {p} }}}=ze^{-\beta \varepsilon _{\mathbf {p} }}}

Sulla base di questi valori medi si può calcolare esplicitamente l'entropia:

  • Bose:
S = k i g i [ log z β ε i z 1 e β ε i 1 log ( 1 z e β ε i ) ] {\displaystyle S=k\sum _{i}g_{i}\left[-{\frac {\log z-\beta \varepsilon _{i}}{z^{-1}e^{\beta \varepsilon _{i}}-1}}-\log \left(1-ze^{-\beta \varepsilon _{i}}\right)\right]}
  • Fermi:
S = k i g i [ log z β ε i z 1 e β ε i + 1 + log ( 1 + z e β ε i ) ] {\displaystyle S=k\sum _{i}g_{i}\left[{\frac {\log z-\beta \varepsilon _{i}}{z^{-1}e^{\beta \varepsilon _{i}}+1}}+\log \left(1+ze^{-\beta \varepsilon _{i}}\right)\right]}
  • Boltzmann:
S = z k i g i e β ε i ( β ε i log z ) {\displaystyle S=zk\sum _{i}g_{i}e^{-\beta \varepsilon _{i}}\left(\beta \varepsilon _{i}-\log z\right)}

Funzione di partizione per i gas ideali quantistici

Per avere una trattazione specifica a partire dai gas ideali quantistici si veda gli insiemi statistici Insieme Microcanonico, canonico e gran canonico trattati secondo la meccanica statistica quantistica e i gas di Fermi e di Bose.

Per i gas ideali trattati con la meccanica statistica quantistica la funzione di partizione canonica è data da:

Z ( V , T ) = { n p } g { n p } e β E { n p } {\displaystyle Z(V,T)=\sum _{\{n_{\mathbf {p} }\}}g\{n_{\mathbf {p} }\}e^{-\beta E\{n_{\mathbf {p} }\}}}

dove sussistono sempre le condizioni:

E { n p } = p ε p n p {\displaystyle E\{n_{\mathbf {p} }\}=\sum _{\mathbf {p} }\varepsilon _{\mathbf {p} }n_{\mathbf {p} }}
p n p = N {\displaystyle \sum _{\mathbf {p} }n_{\mathbf {p} }=N}

Per un gas di Boltzmann:

log Z = N log [ V N ( m k T 2 π 2 ) 3 / 2 ] {\displaystyle \log Z=N\log \left[{\frac {V}{N}}\left({\frac {mkT}{2\pi \hbar ^{2}}}\right)^{3/2}\right]}

Per un gas di Fermi e di Bose invece utilizziamo la funzione di partizione gran canonica:

Z = N = 0 z N Z ( V , T ) = N = 0 { n p } z N e β p ε p n p {\displaystyle {\mathcal {Z}}=\sum _{N=0}^{\infty }z^{N}Z(V,T)=\sum _{N=0}^{\infty }\sum _{\{n_{\mathbf {p} }\}}z^{N}e^{-\beta \sum _{\mathbf {p} }\varepsilon _{\mathbf {p} }n_{\mathbf {p} }}}

ma la doppia sommatoria si può esprimere anche come:

Z = n 0 n 1 [ ( z e β ε 0 ) n 0 ( z e β ε 1 ) n 1 ] = Π p [ n ( z e β ε p ) n ] {\displaystyle {\mathcal {Z}}=\sum _{n_{0}}\sum _{n_{1}}\cdots \left[\left(ze^{-\beta \varepsilon _{0}}\right)^{n_{0}}\cdot \left(ze^{-\beta \varepsilon _{1}}\right)^{n_{1}}\cdots \right]=\Pi _{\mathbf {p} }\left[\sum _{n}\left(ze^{-\beta \varepsilon _{\mathbf {p} }}\right)^{n}\right]}

quindi per un gas di Bose:

Z ( V , T , z ) = Π p 1 1 z e β ε p {\displaystyle {\mathcal {Z}}(V,T,z)=\Pi _{\mathbf {p} }{\frac {1}{1-ze^{-\beta \varepsilon _{\mathbf {p} }}}}}

per un gas di Fermi:

Z ( V , T , z ) = Π p ( 1 + z e β ε p ) {\displaystyle {\mathcal {Z}}(V,T,z)=\Pi _{\mathbf {p} }(1+ze^{-\beta \varepsilon _{\mathbf {p} }})}

Equazioni di stato e numero di particelle

A partire dalla funzione di partizione gran canonica possiamo ricavare immediatamente le equazioni di stato dei gas ideali quantistici:

  • Bose:
P V = k T log Z ( V , T , z ) = p log ( 1 z e β ε p ) {\displaystyle PV=kT\log {\mathcal {Z}}(V,T,z)=-\sum _{\mathbf {p} }\log \left(1-ze^{-\beta \varepsilon _{\mathbf {p} }}\right)}
  • Fermi:
P V = k T log Z ( V , T , z ) = p log ( 1 + z e β ε p ) {\displaystyle PV=kT\log {\mathcal {Z}}(V,T,z)=\sum _{\mathbf {p} }\log \left(1+ze^{-\beta \varepsilon _{\mathbf {p} }}\right)}

Il numero di particelle si ricava subito:

  • Bose:
N = z z log Z ( V , T , z ) = p z e β ε p 1 z e β ε p {\displaystyle N=z{\frac {\partial }{\partial z}}\log {\mathcal {Z}}(V,T,z)=\sum _{\mathbf {p} }{\frac {ze^{-\beta \varepsilon _{\mathbf {p} }}}{1-ze^{-\beta \varepsilon _{\mathbf {p} }}}}}
  • Fermi:
N = z z log Z ( V , T , z ) = p z e β ε p 1 + z e β ε p {\displaystyle N=z{\frac {\partial }{\partial z}}\log {\mathcal {Z}}(V,T,z)=\sum _{\mathbf {p} }{\frac {ze^{-\beta \varepsilon _{\mathbf {p} }}}{1+ze^{-\beta \varepsilon _{\mathbf {p} }}}}}

Voci correlate

  Portale Quantistica: accedi alle voci di Wikipedia che trattano di quantistica