Home/Lezioni/Lezione 7

Struttura matematica: le relazioni fondamentali

Relazioni fondamentali e potenziali termodinamici.

Sistemi apertiRelazioni fondamentaliRappresentazione dell'energiaRappresentazione dell'entropiaRelazioni di MaxwellRelazione di Gibbs-DuhemTeorema di EuleroTerzo principioStabilità termodinamicaConvessità

Nella lezione precedente abbiamo introdotto l'entropia e il secondo principio e ne abbiamo dedotto la relazione fondamentale dU=T dS−P dVdU=T\,dS-P\,dV per un sistema chiuso. Vediamo ora che cosa ci insegna questa relazione sulle variabili che descrivono un sistema termodinamico e sulle loro interdipendenze. Considereremo un sistema semplice: una sostanza pura (un solo componente chimico) e omogenea (le stesse proprietà macroscopiche locali in ogni punto), il cui unico lavoro è quello delle forze di pressione. Nel corso della lezione mostreremo anche come questa descrizione si generalizzi alle miscele, costituite da più specie chimiche.

1. Comprendere la relazione fondamentale

Prima di interpretare questa relazione, estendiamola ai sistemi aperti per affrontare direttamente il caso in cui la quantità di materia può variare.

1.1. Estensione ai sistemi aperti

Nella lezione 6 abbiamo considerato un percorso localmente reversibile tra due stati di equilibrio vicini di un sistema chiuso per ottenere δQrev=T dS\delta Q_{\mathrm{rev}}=T\,dS e δWrev=−P dV\delta W_{\mathrm{rev}}=-P\,dV, e poi

dU=T dS−P dV,dU=T\,dS-P\,dV,

tramite il primo principio. In un sistema aperto, il numero di particelle NN può variare. Il differenziale dell'energia contiene in questo caso un termine aggiuntivo:

Definizione 1 (Relazione fondamentale e potenziale chimico)
Per un sistema aperto costituito da una sola specie chimica e soggetto soltanto al lavoro delle forze di pressione, la relazione fondamentale tra stati di equilibrio vicini è:

dU=TdS−PdV+μdN\boxed{dU=T dS-P dV+\mu dN}
(1)

dove il coefficiente μ\mu è detto potenziale chimico.

Si noti che qui trattiamo NN come una variabile continua, sebbene, a rigore, sia un numero intero. Questa approssimazione è adatta ai sistemi macroscopici, per i quali N≫1N\gg1.

Nonostante il nome, il coefficiente μ\mu interviene anche in assenza di reazioni chimiche. Per capire che cosa misuri, fissiamo SS e VV: rimane dU=μ dNdU=\mu\,dN. Il potenziale chimico misura quindi la variazione dell'energia interna per particella aggiunta, a entropia e volume costanti. Ha le dimensioni di un'energia.

Possiamo svilupparne un'intuizione per analogia con la temperatura: come il calore passa spontaneamente da un mezzo caldo a uno freddo, così le particelle di una specie tendono a passare, a temperatura comune, dal mezzo in cui il loro potenziale chimico è maggiore a quello in cui è minore, se la parete ne consente il passaggio. L'uguaglianza dei potenziali chimici è dunque una condizione di equilibrio rispetto agli scambi di materia, come l'uguaglianza delle temperature è una condizione di equilibrio termico. Lo dimostreremo esplicitamente nella sezione 6.

1.2. Variabili fondamentali ed equazioni di stato

Abbiamo ormai molte variabili termodinamiche: UU, SS, VV, NN, TT, PP e μ\mu. Questa proliferazione è spesso fonte di confusione per gli studenti. Sorge spontanea una domanda: possiamo scegliere queste variabili in modo indipendente? In caso contrario, quali dobbiamo assegnare per determinare le altre? Abbiamo già visto che alcune sono legate da un'equazione di stato, come nel caso del gas perfetto. Ma quante equazioni di stato possiede un sistema semplice e quante variabili realmente indipendenti rimangono?

La relazione fondamentale permette di rispondere a questa domanda cruciale. Ricordiamo infatti la definizione del differenziale di una funzione f(x,y,z)f(x,y,z):

df=(∂f∂x)y,zdx+(∂f∂y)x,zdy+(∂f∂z)x,ydz.df=\drp{f}{x}{y,z}dx+\drp{f}{y}{x,z}dy+\drp{f}{z}{x,y}dz.

Dal confronto con la relazione dU=T dS−P dV+μ dNdU=T\,dS-P\,dV+\mu\,dN seguono due conseguenze.

Proprietà 1 (Variabili naturali ed equazioni di stato)
Anzitutto, l'energia interna si esprime in funzione di sole tre variabili indipendenti: SS, VV e NN. Sono dette variabili fondamentali in questa rappresentazione, chiamata rappresentazione energetica. Si dice anche che sono le variabili naturali di UU, e si scrive U=U(S,V,N)U=U(S,V,N). Inoltre, identificando i coefficienti di dSdS, dVdV e dNdN, otteniamo tre relazioni:

T(S,V,N)=(∂U∂S)V,N,P(S,V,N)=−(∂U∂V)S,N,μ(S,V,N)=(∂U∂N)S,V\boxed{ T(S,V,N)=\drp{U}{S}{V,N}, P(S,V,N)=-\drp{U}{V}{S,N}, \mu(S,V,N)=\drp{U}{N}{S,V} }
(2)

Queste tre relazioni sono le equazioni di stato: danno TT, PP e μ\mu in funzione di SS, VV e NN.

Occorre comprendere bene il ruolo delle diverse variabili in dU=T dS−P dV+μ dNdU=T\,dS-P\,dV+\mu\,dN. La presenza di TT, PP e μ\mu non significa che debbano essere aggiunte a SS, VV e NN come tre ulteriori variabili indipendenti. Sono i coefficienti del differenziale e sono a loro volta funzioni di SS, VV e NN: T=T(S,V,N)T=T(S,V,N), e analogamente per PP e μ\mu. Il contenuto non banale della relazione fondamentale è che queste tre funzioni sono le derivate di una stessa funzione U(S,V,N)U(S,V,N). Le loro variazioni sono quindi già incluse attraverso quelle di SS, VV e NN; non aggiungono tre variabili indipendenti a questa descrizione.

Quanto precede mostra che i sette parametri introdotti non sono tutti indipendenti. Per un dato sistema semplice, tre parametri estensivi, SS, VV e NN, bastano a determinare lo stato di equilibrio: la funzione U(S,V,N)U(S,V,N) dà l'energia interna, e le sue tre derivate danno temperatura, pressione e potenziale chimico.

Bisogna dunque distinguere l'assegnazione dei parametri S,V,NS, V, N, che fissa lo stato del sistema, dall'assegnazione della funzione U(S,V,N)U(S,V,N), che descrive il sistema stesso. Questa funzione è chiamata equazione fondamentale del sistema nella rappresentazione energetica. La relazione differenziale (1) è soddisfatta da tutti i sistemi semplici, ma la funzione UU cambia da un sistema all'altro: non ha la stessa espressione per un gas perfetto, un gas reale, un liquido o un solido.

Osservazione 1 (Generalizzazione a rr specie chimiche)
Per una miscela di rr specie chimiche, occorre assegnare il numero di particelle NiN_i di ciascuna. La relazione fondamentale diventa dU=T dS−P dV+∑i=1rμi dNidU=T\,dS-P\,dV+\sum_{i=1}^{r}\mu_i\,dN_i

Ogni specie possiede così il proprio potenziale chimico. La funzione fondamentale si scrive allora U(S,V,N1,...,Nr)U(S,V,N_1,...,N_r), con

μi=(∂U∂Ni)S,V,Nj≠i.\mu_i=\drp{U}{N_i}{S,V,N_{j\ne i}}.

La derivata è calcolata mantenendo costanti l'entropia, il volume e i numeri di particelle delle altre specie.

1.3. Esempio del gas perfetto monoatomico

Prendiamo l'esempio del gas perfetto monoatomico, per il quale conosciamo già l'espressione

U=32NkBT.U=\frac{3}{2}Nk_BT.

Questa espressione dà UU in funzione di TT e di NN. Per scriverla nella forma U(S,V,N)U(S,V,N), occorre ancora esprimere TT in funzione di SS, VV e NN. Finché non conosciamo questa dipendenza, non possiamo calcolare direttamente le ultime due equazioni di stato mediante le derivate di UU riportate sopra. In particolare, sarebbe un errore concludere che la pressione è nulla perché VV non compare esplicitamente nell'espressione di UU: a SS e NN fissati, TT dipende dal volume. Analogamente, per calcolare μ\mu, bisogna tenere conto della dipendenza di TT da NN a SS e VV fissati. Le equazioni di stato danno in realtà

P=−32NkB(∂T∂V)S,N,μ=32kBT+32NkB(∂T∂N)S,V.P=-\frac{3}{2}Nk_B\drp{T}{V}{S,N}, \qquad \mu=\frac{3}{2}k_BT+\frac{3}{2}Nk_B\drp{T}{N}{S,V}.

Calcoleremo queste derivate parziali nella sezione 5: la prima permette di ritrovare l'equazione di stato PV=NkBTPV=Nk_BT, mentre la seconda fornisce il potenziale chimico del gas perfetto.

1.4. Le relazioni di Maxwell

Le relazioni di Maxwell sono identità termodinamiche, valide qualunque sia la sostanza considerata. Derivano dall'uguaglianza delle derivate parziali miste: se le derivate seconde di UU sono continue, il teorema di Schwarz permette di scambiare l'ordine delle derivazioni.

Applichiamolo alle variabili SS e VV, a NN fissato:

∂2U∂V∂S=∂2U∂S∂V.\frac{\partial^2U}{\partial V\partial S} =\frac{\partial^2U}{\partial S\partial V}.

Il membro sinistro indica che si deriva prima UU rispetto a SS, a VV e NN costanti, poi la funzione ottenuta rispetto a VV, a SS e NN costanti. Ora, secondo le equazioni di stato (2), (∂U∂S)V,N=T\drp{U}{S}{V,N}=T. Questo membro vale quindi

∂2U∂V∂S=(∂∂V)S,N[(∂U∂S)V,N]=(∂T∂V)S,N.\frac{\partial^2U}{\partial V\partial S} =\left(\frac{\partial}{\partial V}\right)_{S,N} \left[\drp{U}{S}{V,N}\right] =\drp{T}{V}{S,N}.

A destra si deriva prima UU rispetto a VV, a SS e NN costanti, poi la funzione ottenuta rispetto a SS, a VV e NN costanti. Poiché (∂U∂V)S,N=−P\drp{U}{V}{S,N}=-P, si ha

∂2U∂S∂V=(∂∂S)V,N[(∂U∂V)S,N]=−(∂P∂S)V,N.\frac{\partial^2U}{\partial S\partial V} =\left(\frac{\partial}{\partial S}\right)_{V,N} \left[\drp{U}{V}{S,N}\right] =-\drp{P}{S}{V,N}.

L'uguaglianza dei due membri dà così la prima relazione di Maxwell:

(∂T∂V)S,N=−(∂P∂S)V,N.\drp{T}{V}{S,N}=-\drp{P}{S}{V,N}.

Procedendo allo stesso modo con le coppie (S,N)(S,N) e (V,N)(V,N), otteniamo altre due relazioni. Le tre relazioni di Maxwell nella rappresentazione energetica si scrivono quindi

(∂T∂V)S,N=−(∂P∂S)V,N,(∂T∂N)S,V=(∂μ∂S)V,N,(∂P∂N)S,V=−(∂μ∂V)S,N\boxed{\begin{aligned} \drp{T}{V}{S,N} &=-\drp{P}{S}{V,N},\\ \drp{T}{N}{S,V} &=\drp{\mu}{S}{V,N},\\ \drp{P}{N}{S,V} &=-\drp{\mu}{V}{S,N} \end{aligned}}

Queste relazioni discendono direttamente dalla relazione fondamentale quando UU è almeno due volte continuamente differenziabile. Costituiscono quindi una previsione forte della teoria termodinamica, indipendente dal modello: qualsiasi sostanza descritta da una funzione fondamentale di questo tipo deve soddisfarle (sotto queste ipotesi di regolarità di UU).

Viceversa, se le funzioni TT, PP e μ\mu soddisfano le relazioni di Maxwell (e sono almeno C1C^1), esiste localmente una funzione UU tale che dU=T dS−P dV+μ dNdU=T\,dS-P\,dV+\mu\,dN. Si può quindi risalire dalle equazioni di stato alla funzione fondamentale per integrazione, a meno di una costante additiva.

Questo procedimento inverso è importante nella pratica: quando si studia una sostanza in laboratorio, non se ne conosce fin dall'inizio la funzione fondamentale U(S,V,N)U(S,V,N). Misurare l'entropia stessa non è immediato; torneremo sull'argomento nella prossima lezione. Gli esperimenti forniscono piuttosto relazioni tra grandezze misurabili, come pressione, volume e temperatura. Lo abbiamo visto nella lezione storica: l'equazione del gas perfetto PV=NkBTPV=Nk_BT fu stabilita dapprima sperimentalmente. Occorre poi integrare queste relazioni per risalire all'equazione fondamentale della sostanza studiata. Effettueremo questa ricostruzione per il gas perfetto nella sezione 5, ottenendo così la sua equazione fondamentale U(S,V,N)U(S,V,N).

2. La rappresentazione entropica

Nella rappresentazione energetica, UU è la variabile dipendente: la esprimiamo come funzione di SS, VV e NN. Possiamo anche scegliere l'entropia come variabile dipendente, purché sia possibile invertire U(S,V,N)U(S,V,N) per ottenere S(U,V,N)S(U,V,N) a VV e NN fissati. Nel seguito ci collocheremo in regioni in cui questa inversione è possibile, senza ricordarlo ogni volta.

Isolando dSdS nella relazione fondamentale, per T≠0T\ne0, si ottiene

dS=1TdU+PTdV−μTdN\boxed{dS=\frac{1}{T}dU+\frac{P}{T}dV-\frac{\mu}{T}dN}
(3)

Questa è la relazione fondamentale nella rappresentazione entropica. Con lo stesso ragionamento usato prima, le variabili naturali di SS sono dunque UU, VV e NN:

S=S(U,V,N)\boxed{ S=S(U,V,N) }

L'identificazione dei coefficienti del differenziale fornisce ancora tre equazioni di stato

1T=(∂S∂U)V,NPT=(∂S∂V)U,N−μT=(∂S∂N)U,V\boxed{\frac{1}{T}=\drp{S}{U}{V,N}} \qquad \boxed{\frac{P}{T}=\drp{S}{V}{U,N}} \qquad \boxed{-\frac{\mu}{T}=\drp{S}{N}{U,V}}
(4)

Attenzione: non si tratta di tre equazioni di stato indipendenti da quelle ottenute nella rappresentazione energetica; esprimono la stessa informazione, ma con una diversa scelta di variabili. La funzione S(U,V,N)S(U,V,N) si chiama equazione fondamentale nella rappresentazione entropica. Contiene la stessa informazione di U(S,V,N)U(S,V,N).

L'uguaglianza delle derivate miste di SS fornisce le relazioni di Maxwell associate a questa rappresentazione:

(∂1/T∂V)U,N=(∂P/T∂U)V,N,(∂1/T∂N)U,V=−(∂μ/T∂U)V,N,(∂P/T∂N)U,V=−(∂μ/T∂V)U,N\boxed{\begin{aligned} \drp{1/T}{V}{U,N} &=\drp{P/T}{U}{V,N},\\ \drp{1/T}{N}{U,V} &=-\drp{\mu/T}{U}{V,N},\\ \drp{P/T}{N}{U,V} &=-\drp{\mu/T}{V}{U,N} \end{aligned}}

Queste relazioni non sono indipendenti da quelle ottenute nella rappresentazione energetica; sono semplicemente quelle adatte a questa scelta di variabili.

Da ricordare
Ciascuna delle funzioni U(S,V,N)U(S,V,N) e S(U,V,N)S(U,V,N) costituisce un'equazione fondamentale: conoscere l'una o l'altra basta a determinare tutta la termodinamica di equilibrio del sistema, in particolare le sue variabili intensive, le sue equazioni di stato e le sue relazioni di Maxwell. Geometricamente, queste due funzioni descrivono la stessa ipersuperficie tridimensionale nello spazio delle variabili estensive (U,S,V,N)(U,S,V,N). Questa superficie di equilibrio rappresenta lo spazio degli stati di equilibrio E\mathcal E, introdotto nella lezione 3 e ripreso nella lezione 4. In quest'ultima si utilizzava, per esempio, il grafico di U(T,V,N)U(T,V,N) nello spazio (U,T,V,N)(U,T,V,N): gli stati rappresentati sono gli stessi, cambia soltanto la scelta delle variabili.

3. Il ruolo dell'estensività

Del gas perfetto conosciamo ancora soltanto due equazioni di stato esplicite. Non bastano per ricostruire interamente l'equazione fondamentale per integrazione, come prospettato nella sezione 1.4. Possiamo però sfruttare una proprietà aggiuntiva: l'estensività. Essa permette di applicare il teorema di Eulero a UU o a SS e di dedurne una nuova relazione non banale, la relazione di Gibbs-Duhem.

La studiamo qui in forma generale, prima di applicarla al caso del gas perfetto.

3.1. Estensività e teorema di Eulero

Definizione 2 (Funzione omogenea)
In matematica, una funzione f(x1,...,xn)f(x_1,...,x_n) si dice omogenea di grado kk se soddisfa f(λx1,...,λxn)=λkf(x1,...,xn)f(\lambda x_1,...,\lambda x_n)=\lambda^k f(x_1,...,x_n)

per ogni λ>0\lambda>0.

Se ff è differenziabile, il teorema di Eulero afferma che

∑i=1nxi∂f∂xi=kf\boxed{\sum_{i=1}^{n}x_i\frac{\partial f}{\partial x_i}=k f}

(5)

Dimostrazione.
Basta derivare la relazione di omogeneità rispetto a λ\lambda, mantenendo fissi gli xix_i: ∑i=1nxi∂f∂xi(λx1,...,λxn)=kλk−1f(x1,...,xn).\sum_{i=1}^{n}x_i\frac{\partial f}{\partial x_i} (\lambda x_1,...,\lambda x_n) =k\lambda^{k-1}f(x_1,...,x_n).

Si pone poi λ=1\lambda=1.

Abbiamo discusso l'estensività di UU e di SS nelle lezioni 4 e 6. Questa proprietà vale per i sistemi macroscopici nei quali si possono trascurare gli effetti di superficie e le interazioni a lungo raggio. Supporremo che sia così in questo caso. Allora:

U(λS,λV,λN)=λU(S,V,N),S(λU,λV,λN)=λS(U,V,N),\begin{aligned} U(\lambda S,\lambda V,\lambda N)&=\lambda U(S,V,N),\\ S(\lambda U,\lambda V,\lambda N)&=\lambda S(U,V,N), \end{aligned}

cioè UU e SS sono funzioni omogenee di grado uno. Il teorema di Eulero applicato a U(S,V,N)U(S,V,N) con k=1k=1 dà quindi:

U=S(∂U∂S)V,N+V(∂U∂V)S,N+N(∂U∂N)S,V,U=S\drp{U}{S}{V,N} +V\drp{U}{V}{S,N} +N\drp{U}{N}{S,V},

e, sostituendo le equazioni di stato (2), otteniamo la relazione di Eulero:

U=TS−PV+μN\boxed{U=TS-PV+\mu N}
(6)

Osservazione 2
La relazione fondamentale differenziale non presuppone l'estensività. Questa forma integrata, invece, ne dipende. Un calcolo analogo nella rappresentazione entropica fornisce la stessa equazione (esercizio: verificarlo).

3.2. La relazione di Gibbs-Duhem

La relazione di Eulero conduce a un'importante identità termodinamica per i sistemi estensivi. Differenziandola, troviamo

dU=T dS+S dT−P dV−V dP+μ dN+N dμ.dU=T\,dS+S\,dT-P\,dV-V\,dP+\mu\,dN+N\,d\mu.

Raggruppiamo i termini a sinistra come segue:

dU−T dS+P dV−μ dN=S dT−V dP+N dμ.dU-T\,dS+P\,dV-\mu\,dN=S\,dT-V\,dP+N\,d\mu.

Riconosciamo allora che il membro sinistro è nullo in virtù della relazione fondamentale. Rimane quindi la relazione

SdT−VdP+Ndμ=0\boxed{S dT-V dP+N d\mu=0}
(7)

chiamata relazione di Gibbs-Duhem. Mostra che le tre variabili intensive TT, PP e μ\mu non possono variare indipendentemente in una fase omogenea costituita da una sola specie.

Osservazione 3 (Nella rappresentazione entropica)
Nella rappresentazione entropica, la relazione di Eulero si scrive S=UT+PTV−μTN.S=\frac{U}{T}+\frac{P}{T}V-\frac{\mu}{T}N.

Differenziandola e usando la relazione fondamentale (3), otteniamo analogamente

Ud(1T)+Vd(PT)−Nd(μT)=0\boxed{U d\left(\frac{1}{T}\right) +V d\left(\frac{P}{T}\right) -N d\left(\frac{\mu}{T}\right)=0}
(8)

Osservazione 4 (Generalizzazione a più specie)
Lo stesso calcolo, a partire dalla relazione fondamentale di una miscela di rr specie chimiche, dà U=TS−PV+∑i=1rμiNi,U=TS-PV+\sum_{i=1}^{r}\mu_iN_i,

e la relazione di Gibbs-Duhem diventa

S dT−V dP+∑i=1rNi dμi=0\boxed{ S\,dT-V\,dP+\sum_{i=1}^{r}N_i\,d\mu_i=0 }

4. I gradi di libertà intensivi

Abbiamo descritto un sistema semplice mediante tre variabili estensive, SS, VV e NN. Quando il sistema è estensivo, queste tre variabili si riducono a due variabili intensive per descriverne lo stato intensivo. Infatti, l'estensività dà immediatamente

U(S,V,N)=Nu(s,v),U(S,V,N)=N u(s,v),

dove l'energia per particella u=U/Nu=U/N è dunque una funzione di sole due variabili intensive, s=S/Ns=S/N e v=V/Nv=V/N.

Definizione 3 (Gradi di libertà intensivi)
Si chiama numero di gradi di libertà di un sistema termodinamico il numero di variabili intensive indipendenti che si possono scegliere per determinarne lo stato intensivo di equilibrio. È pari a due per una sostanza pura in un'unica fase.

Per descrivere lo stato completo del sistema, occorre ancora specificarne la dimensione, per esempio assegnando NN o un'altra grandezza estensiva adatta. Non bisogna dunque confondere i due gradi di libertà intensivi con le tre variabili estensive necessarie per descrivere lo stato completo.

Generalizzeremo questo conteggio ai sistemi con più componenti e più fasi nella lezione 10, dove stabiliremo la regola delle fasi di Gibbs.

5. L'equazione fondamentale del gas perfetto

5.1. L'equazione di Sackur-Tetrode

Le equazioni del gas perfetto monoatomico sono PV=NkBTPV=Nk_BT e U=32NkBTU=\frac{3}{2}Nk_BT, dove kBk_B è la costante di Boltzmann, legata alla costante dei gas perfetti RR da kB=R/NAk_B=R/N_A, con NAN_A numero di Avogadro. Queste due equazioni e la relazione di Gibbs-Duhem permettono di trovare l'equazione fondamentale del gas perfetto.

È più semplice cercare SS che UU. Partiamo quindi dalla prima equazione di stato nella rappresentazione entropica:

(∂S∂U)V,N=1T=3NkB2U.\drp{S}{U}{V,N} =\frac{1}{T} =\frac{3Nk_B}{2U}.

Integrando rispetto a UU a VV e NN fissati, con un'energia di riferimento arbitraria U0>0U_0>0, troviamo

S(U,V,N)=32NkBln⁡(UU0)+F(V,N),S(U,V,N)=\frac{3}{2}Nk_B\ln\left(\frac{U}{U_0}\right)+F(V,N),
(9)

dove FF è una funzione ancora sconosciuta. L'argomento del logaritmo è così adimensionale; alla fine U0U_0 sarà raggruppata con le altre costanti. La seconda equazione di stato dà poi

(∂S∂V)U,N=PT=NkBV.\drp{S}{V}{U,N} =\frac{P}{T} =\frac{Nk_B}{V}.

Sostituendo (9) e integrando rispetto a VV, otteniamo

F(V,N)=NkBln⁡(VV0)+Z(N),F(V,N)=Nk_B\ln\left(\frac{V}{V_0}\right)+Z(N),

dove anche V0>0V_0>0 è un volume di riferimento arbitrario, e quindi

S(U,V,N)=32NkBln⁡(UU0)+NkBln⁡(VV0)+Z(N).S(U,V,N) =\frac{3}{2}Nk_B\ln\left(\frac{U}{U_0}\right) +Nk_B\ln\left(\frac{V}{V_0}\right) +Z(N).

Resta da determinare Z(N)Z(N) con l'aiuto della relazione di Gibbs-Duhem. Il calcolo è svolto nella dimostrazione qui sotto. Si ottiene, a meno di una costante AA,

S=NkB[ln⁡(AVN(UN)3/2)+52]\boxed{ S =Nk_B\left[ \ln\left( A\frac{V}{N}\left(\frac{U}{N}\right)^{3/2} \right) +\frac{5}{2} \right] }
(10)

dove AA è una costante dimensionale, indipendente da UU, VV e NN. A meno del valore di AA, questa è la formula di Sackur-Tetrode, ricavata indipendentemente da Otto Sackur e Hugo Tetrode nel 19121.

Note 1 : La termodinamica determina la dipendenza di SS da UU, VV e NN, ma non può fissare AA. È proprio questo il contributo di Sackur e Tetrode: dividendo lo spazio delle fasi in celle di volume h3h^3 per particella, ottennero, per un gas perfetto monoatomico di particelle indistinguibili di massa mm, senza degenerazione interna, A=(4πm3h2)3/2.A=\left(\frac{4\pi m}{3h^2}\right)^{3/2}. La presenza della costante di Planck hh mostra che questo valore appartiene alla fisica quantistica e non alla sola termodinamica.
Dimostrazione.
La terza equazione di stato dà

μT=−(∂S∂N)U,V=−32kBln⁡(UU0)−kBln⁡(VV0)−Z′(N).\frac{\mu}{T} =-\drp{S}{N}{U,V} =-\frac{3}{2}k_B\ln\left(\frac{U}{U_0}\right) -k_B\ln\left(\frac{V}{V_0}\right)-Z'(N).
(11)

D'altra parte, la relazione di Gibbs-Duhem nella rappresentazione entropica, equazione (8), si scrive

d(μT)=UN d(1T)+VN d(PT)=UN d(3NkB2U)+VN d(NkBV)=3kB2(dNN−dUU)+kB(dNN−dVV)=5kB2dNN−3kB2dUU−kBdVV.\begin{aligned} d\left(\frac{\mu}{T}\right) &=\frac{U}{N}\,d\left(\frac{1}{T}\right) +\frac{V}{N}\,d\left(\frac{P}{T}\right)\\ &=\frac{U}{N}\,d\left(\frac{3Nk_B}{2U}\right) +\frac{V}{N}\,d\left(\frac{Nk_B}{V}\right)\\ &=\frac{3k_B}{2}\left(\frac{dN}{N}-\frac{dU}{U}\right) +k_B\left(\frac{dN}{N}-\frac{dV}{V}\right)\\ &=\frac{5k_B}{2}\frac{dN}{N} -\frac{3k_B}{2}\frac{dU}{U}-k_B\frac{dV}{V}. \end{aligned}

Questa espressione si integra direttamente, dando

μT=−3kB2ln⁡(UU0)−kBln⁡(VV0)+5kB2ln⁡N−α,\frac{\mu}{T} =-\frac{3k_B}{2}\ln\left(\frac{U}{U_0}\right) -k_B\ln\left(\frac{V}{V_0}\right) +\frac{5k_B}{2}\ln N-\alpha,

dove α\alpha è una costante ancora indeterminata. Confrontando con l'equazione (11), otteniamo

Z′(N)=−5kB2ln⁡N+α.Z'(N)=-\frac{5k_B}{2}\ln N+\alpha.

Integriamo ancora una volta:

Z(N)=−5kB2Nln⁡N+5kB2N+αN+β,Z(N)=-\frac{5k_B}{2}N\ln N+\frac{5k_B}{2}N+\alpha N+\beta,

dove β\beta è una seconda costante di integrazione. Sostituendo questa espressione in quella dell'entropia e raggruppando i logaritmi, troviamo

S(U,V,N)=NkBln⁡[(UNU0)3/2VNV0]+52NkB+αN+β.S(U,V,N) =Nk_B\ln\left[ \left(\frac{U}{NU_0}\right)^{3/2}\frac{V}{NV_0} \right] +\frac{5}{2}Nk_B+\alpha N+\beta.

Tutti i termini tranne β\beta sono estensivi. Per avere S(λU,λV,λN)=λS(U,V,N)S(\lambda U,\lambda V,\lambda N)=\lambda S(U,V,N), deve quindi essere β=0\beta=0. Raggruppando U0U_0, V0V_0 e α\alpha nella costante A=eα/kB/(U03/2V0)A=e^{\alpha/k_B}/(U_0^{3/2}V_0), otteniamo

S(U,V,N)=NkB[ln⁡(AVN(UN)3/2)+52],S(U,V,N) =Nk_B\left[ \ln\left( A\frac{V}{N}\left(\frac{U}{N}\right)^{3/2} \right) +\frac{5}{2} \right],
(12)

dove AA possiede le dimensioni necessarie.

Osservazione 5 (Nella rappresentazione energetica)
Invertendo la relazione (10), si ottiene l'equazione fondamentale del gas perfetto nella rappresentazione energetica: U(S,V,N)=N(NAV)2/3exp⁡(2S3NkB−53)\boxed{ U(S,V,N)=N\left(\frac{N}{AV}\right)^{2/3} \exp\left(\frac{2S}{3Nk_B}-\frac{5}{3}\right) }

5.2. Ritrovare le equazioni di stato

Verifichiamo che l'equazione fondamentale S=S(U,V,N)S=S(U,V,N) data in (10) permetta di ritrovare le due equazioni di stato iniziali. Derivando rispetto a UU e a VV, troviamo

1T=(∂S∂U)V,N=3NkB2U,\frac{1}{T} = \drp{S}{U}{V,N}=\frac{3Nk_B}{2U},

e

PT=(∂S∂V)U,N=NkBV.\frac{P}{T} = \drp{S}{V}{U,N}=\frac{Nk_B}{V}.

Ritroviamo dunque U=3NkBT/2U=3Nk_BT/2 e PV=NkBTPV=Nk_BT. Resta soltanto da calcolare il potenziale chimico μ=−T(∂S∂N)U,V\mu=-T\drp{S}{N}{U,V}. Il calcolo della derivata dà (verificarlo):

μ=−kBTln⁡[AVN(UN)3/2]\boxed{\mu=-k_BT\ln\left[A\frac{V}{N}\left(\frac{U}{N}\right)^{3/2}\right]}

(13)

Osservazione 6
Si possono allora verificare esplicitamente le relazioni di Maxwell per il gas perfetto; si veda l'esercizio «Verificare le relazioni di Maxwell per il gas perfetto».

6. Equilibrio e stazionarietà dell'entropia

Nella lezione 6, nel paragrafo «L'entropia e l'equilibrio», abbiamo ammesso che l'entropia di un insieme isolato è stazionaria all'equilibrio: la sua variazione al primo ordine è nulla quando si fanno variare leggermente i parametri liberi del sistema rispettando i vincoli imposti. Riprendiamo questo calcolo per due sottosistemi che possono anche scambiare particelle. Si suppone che ciascuno sia in equilibrio interno e si trascurano i contributi della loro interfaccia. L'insieme è isolato: la sua energia UU, il suo volume VV e il suo numero di particelle NN sono fissati. I vincoli si scrivono quindi

U1+U2=U,V1+V2=V,N1+N2=N.U_1+U_2=U,\qquad V_1+V_2=V,\qquad N_1+N_2=N.

I parametri di ciascun sottosistema possono tuttavia variare, purché le loro variazioni si compensino:

dU2=−dU1,dV2=−dV1,dN2=−dN1.dU_2=-dU_1,\qquad dV_2=-dV_1,\qquad dN_2=-dN_1.

La parete può imporre ulteriori vincoli: se è fissa, dV1=0dV_1=0; se è impermeabile, dN1=0dN_1=0. Occorre quindi precisare gli scambi consentiti per sapere quali parametri restano liberi.

Esprimendo i parametri del secondo sottosistema mediante i vincoli, l'entropia totale diventa una funzione di U1U_1, V1V_1 e N1N_1:

Stot=S1(U1,V1,N1)+S2(U−U1,V−V1,N−N1).S_{\mathrm{tot}} =S_1(U_1,V_1,N_1)+S_2(U-U_1,V-V_1,N-N_1).

La relazione fondamentale dà

dStot=(1T1−1T2)dU1+(P1T1−P2T2)dV1−(μ1T1−μ2T2)dN1.\begin{aligned} dS_{\mathrm{tot}} ={}&\left(\frac{1}{T_1}-\frac{1}{T_2}\right)dU_1\\ &+\left(\frac{P_1}{T_1}-\frac{P_2}{T_2}\right)dV_1\\ &-\left(\frac{\mu_1}{T_1}-\frac{\mu_2}{T_2}\right)dN_1. \end{aligned}

Se i tre parametri liberi U1U_1, V1V_1 e N1N_1 possono variare indipendentemente in entrambi i sensi nell'intorno dell'equilibrio, la stazionarietà impone che il coefficiente di ciascuna delle variazioni dU1dU_1, dV1dV_1 e dN1dN_1 sia nullo. Si ottiene successivamente:

T1=T2P1=P2μ1=μ2\boxed{T_1=T_2} \qquad \boxed{P_1=P_2} \qquad \boxed{\mu_1=\mu_2}

L'uguaglianza dei potenziali chimici completa così gli equilibri termico e meccanico già incontrati. Queste condizioni valgono solo per gli scambi consentiti: una parete impermeabile, per esempio, impone dN1=0dN_1=0 e quindi non richiede l'ultima condizione all'equilibrio.

Per una parete permeabile ma fissa, a temperatura comune T1=T2=TT_1=T_2=T, si ha semplicemente

dStot=μ2−μ1T dN1,dS_{\mathrm{tot}}=\frac{\mu_2-\mu_1}{T}\,dN_1,

il che mostra che, se μ2>μ1\mu_2>\mu_1, un trasferimento di particelle da 2 a 1 (dN1>0dN_1>0) aumenta l'entropia e quindi avviene spontaneamente. Come anticipato, a temperatura uniforme la materia tende a passare dalle regioni a potenziale chimico maggiore a quelle a potenziale chimico minore, fino all'uguaglianza dei potenziali chimici all'equilibrio.

7. Riferimenti

  1. H. B. Callen, Thermodynamics and an Introduction to Thermostatistics, 2nd ed., Wiley (1985)