Boltzmann-fordeling

Boltzmann var I slutten av 20-årene da han fant sin statistiske fordeling.

Boltzmann-fordelingen er en statistisk sannsynlighetsfordeling for å beskrive energien til identiske partikler i et system som er i termisk likevekt. Den ble funnet av den østerrikske fysikeren Ludwig Boltzmann i siste halvdel av 1800-tallet og danner grunnlaget for den videre utviklingen av statistisk mekanikk.

I sin opprinnelige form angir fordelingen sannsynligheten p  for å finne en partikkel med energi E  i et systemet av frie partikler med likevektstemperatur T,

$p=\frac{1}{Z}{e}^{-E/k_{B}T}$

Her bestemmer Z  normeringen av fordelingen og k_B  er en konstant som I ettertid har fått navnet Boltzmanns konstant. For en ideell gass der hver partikkel har energi Mal:Nowrap går Boltzmann-fordelingen over til å gi Maxwell-fordelingen for hastighetene i gassen. Det er derfor vanlig å si at den beskrives ved Maxwell-Boltzmann statistikk.

Fordelingsfunksjonen viste seg å ha større gyldighet enn for energien til enkelte partikler i en gass. Spesielt arbeidene til den amerikanske fysiker Willard Gibbs viste at samme sannsynlighetsfordeling også gjelder for den totale energien til et system av partikler med gjensidige vekselvirkninger. Det beskrives da i det som kalles det kanoniske ensemblet av identiske systemer.

Den statistiske mekanikken til Boltzmann spilte en avgjørende rolle da Max Planck kunne forklare egenskapene ved varmestråling i 1900. Dette markerte begynnelsen på kvantemekanikken. Det var Planck som ga navnet til Boltzmanns konstant. Fem år senere gjorde Albert Einstein lignende, men enda dypere betraktninger rundt systemet av stråling og kom frem til at den besto av fotoner.

Bakgrunn og utledning

Etter at James Clerk Maxwell i 1860 hadde publisert sin utledning av sannsynlighetsfordelingen for hastighetene i en gass, begynte Ludwig Boltzmann en lang periode med intense studier av den generelle gyldigheten av denne Maxwell-fordelingen og hvordan den kan opptre for systemer som opprinnelig ikke befinner seg i termisk likevekt. Med årene ble disse spørsmålene avklart, og i dag eksisterer det mange forskjellige måter å utlede denne mer generelle Boltzmann-fordelingen.

En statistisk beskrivelse av energien E  til én partikkel som befinner seg i en bestemt bevegelse eller mikrotilstand sammen med et større antall av like partikler, er gitt ved å finne en sannsynlighetsfordelinge p(E). Den relative sannsynligheten p(E₁)/p(E₂) for å finne den med to forskjellige energier, kan ikke være avhengig av hvordan man setter nullpunktet for disse energiene. Derfor må man ha at dette forholdet bare avhenger av differensen mellom de to energiene,

$\frac{p(E_1)}{p(E_2)}=f(E_1-E_2)$

Ved å betrakte en tredje energi E₃ kombinert ved observasjonen at p(E₁)/p(E₃) = p(E₁)/p(E₂)×p(E₂)/p(E₃), følger det dermed at

$f(E_1-E_3)=f(E_1-E_2)f(E_2-E_3)$

Den eneste funksjonen som tilfredsstiller ligningen, er eksponensialfunksjonen. Dermed må fordelingsfunksjonen ha formen

$p(E)=\frac{1}{Z}{e}^{-\beta E}$

hvor β  er en konstant. Sammenligning med Maxwell-fordelingen viser at Mal:Nowrap der k_B  er Boltzmanns konstant og T  er likevektstemperaturen. Den ukjente faktoren Z  kan bestemmes ved normalisering av sannsynlighetsfordelingen.

Boltzmanns entropi

Forsiden av Boltzmanns forelesninger om gassteori, utgitt 1896.

I sitt arbeid med å forstå Maxwells hastighetsfordeling og hvordan den var kompatibel med termodynamikkens andre hovedsetning, kunne Boltzmann forbinde et systems entropi S  med antall mikrotilstander W  som det har tilgang til. Denne sammenhengen ble av Max Planck skrevet som

$S=k_B\log W$

hvor man benytter den naturlige logaritmen. Det var ved denne formelen han definerte Boltzmanns konstant. Da antall mikrotilstander prinsipielt kan beregnes også når systemet er ute av likevekt, kan formelen også forklare hvordan entropien hele tiden vokser ved at systemet søker mot en likevektstilstand inneholdende et maksimalt antall mikrotilstander.

For et system av partikler med total energi E_T  vil sannnsynligheten p(E) for at en partikkel har energien E  i én bestemt mikrotistand, være proporsjonal med antall tilstander W(E_T - E) som forblir tilgjengelig for de resterende partiklene. Ved å benytte at E << E_T, kan man knytte dette direkte til systemets entropi,

$\log W(E_T-E)=\log W(E_T)-\frac{E}{k_{B}T}$

Dette følger fra den vanlige, termodynamiske deriverte $\partial S/\partial U=1/T.$ Da den søkte sannsynligheten $p(E)\propto W(E_T-E),$ har man igjen resultatet

$p(E)=\frac{1}{Z}{e}^{-E/k_{B}T}$

Argumentet i eksponensialfunksjonen blir ved denne utledningen automatisk bestemt uten å måtte sammenligne med kinetisk teori for en ideell gass.

Generell gyldighet

Begge utledningene av Boltzmanns sannsynlighetsfordeling er basert på det mentale bildet at den gjelder for én partikkel i termisk likevekt med andre partikler i en gass. Men denne partikkelen kunne selv ha en indre struktur som for eksempel et molekyl som består av flere atomer. Disse atomene kunne igjen vibrere om en likevektstilling som en harmonisk oscillator. Sannsynlighetsfordelingen er derfor gyldig for mer generelle, systemer enn en strukturløs partikkel. Generelt gir Boltzmann-fordelingen sannsynligheten for å finne et vilkårlig system i en bestemt mikrotilstand når det er i termisk likevekt med et større varmereservoir som vedlikeholder det med en konstant temperatur. Det gjelder derfor ikke bare for en ideell gass, men også for en gass hvor partiklene vekselvirker med hverandre. Det faller da sammen med sannsynlighetsfordelingen i det kanoniske ensemblet som ble utarbeidet av Josiah Willard Gibbs omtrent på samme tid som Boltzmann.

Ising-spinn

Kanskje det enkleste system man kan beskrive ved Boltzmann-fordelingen, er et enkelt spinn som er lokalisert for eksempel i et fast stoff og kan innta to retninger Mal:Nowrap. Det kan derfor kun befinne seg i to mikrotilstander og kalles ofte i denne sammenhengen for et «Ising-spinn». Når det er i termisk likevekt, vil det ved observasjon peke like mange ganger i den ene som i den andre retningen. Men hvis det har et magnetisk moment Mal:Nowrap og er omgitt av et magnetfelt B, vil det ha en energi Mal:Nowrap som kan anta de to verdien ±μB. Sannsynligheten for å befinne seg i en bestemt tilstand er da

$p_{\sigma }=\frac{1}{Z}{e}^{-\beta E_{\sigma }}$

hvor $\beta =1/k_{B}T.$ Normeringskonstanten Z  er bestemt ved at summen over de to sannsynlighetene er én, det vil si

$Z=e^{\mu B/k_{B}T}+e^{-\mu B/k_{B}T}=2\cosh \frac{\mu B}{k_{B}T}$

Generelt involverer den en sum over alle mikrotilstander og kalles for systemets partisjonsfunksjon. Betegnelsen Z  stammer fra det tyske ordet Zustandssumme for tilstandssum.

Den midlere energien til spinnet er på tilsvarende måte

$⟨E⟩=\frac{1}{Z}\sum _{\sigma =\pm 1}{E}_{\sigma }{e}^{-\beta E_{\sigma }}=-\frac{\partial }{\partial \beta }\ln Z$

Dette er en generell relasjon som følger fra definisjonen av partisjonsfunksjonen. Da man her kan skrive middelverdien $⟨E⟩=\mu ⟨m⟩,$ blir den midlere magnetiseringen av spinnet $⟨m⟩=\mu \tanh (\mu B/k_{B}T).$ Når det ytre feltet er så sterkt at $\mu B\gg k_{B}T,$ vil det derfor overvinne de termiske påvirkningene og spinnet peker alltid i samme retning langs feltet. I det motsatte tilfellet er disse så store at den midlere magnetiseringen forblir lik null.

Harmonisk oscillator

Energien til en partikkel som beveger seg som en harmonisk oscillator i én dimensjon med vinkelfrekvens ω  og utslag q, er

$E=\frac{p^2}{2m}+\frac{1}{2}m{\omega }^2{q}^2$

der p  er dens impuls. Hver mikrotilstand er nå gitt ved å angi dens posisjon og impuls. Men disse to variable q  og p  tar kontinuerlige verdier, og det er ikke opplagt hvordan man skal kunne telle opp antall slike tilstander i dette tilfellet. Dette problemet ble løst ved etableringen av kvantemekanikken og den resulterende Heisenbergs uskarphetsrelasjon. Den sier at disse to variable ikke kan betemmes samtidig med en mindre nøyaktighet enn hva som er gitt ved Plancks konstant h. Konkret betyr det at et infinitesemalt areal dq dp  i det kombinerte impuls-posisjonsrommet inneholder dq dp/h  mikrotilstander. En mer stringent, matematisk utledning gir samme resultat som også dimensjonsmessig er korrekt.

Summasjonen over alle mikrotilstander i partisjonsfunksjonen for oscillatoren kan dermed gjennomføres ved en integrasjon over alle mulige verdier for q  og p . Det gir

$\begin{array}{cc}Z& =\frac{1}{h}\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dq\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}dp{e}^{-(p^2+m^2{\omega }^2{q}^2)/2m{k}_{B}T}\\ & =\frac{k_{B}T}{\hslash \omega }\end{array}$

ved å benytte Gauss-integralet og den reduserte Planck-konstanten Mal:Nowrap. Den indre energien $U$ til oscillatoren er gitt ved middelverdien $⟨E⟩=-(\partial /\partial \beta )\ln Z$ og blir dermed $U=k_{B}T.$ Det er i overenstemmelse med ekvipartisjonsprinsippet som sier at hver kvadratisk variabel i energien til et system bidrar med $k_{B}T/2$ til den indre energien.

Kvantisert oscillator

Den samme Boltzmann-fordelingen kan også benyttes for en kvantisert harmonisk oscillator. Da kan den bare ha energiene $E_n=\hslash \omega (n+1/2)$ hvor kvantetallet n = 0, 1, 2 og så videre. Det betegner nå de kvantemekaniske mikrotistandene slik at den tilsvarende partisjonsfunksjonen blir

$Z=\sum _{n=0}^{\mathrm{\infty }}{e}^{-\beta \hslash \omega (n+1/2)}=\frac{e^{-\beta \hslash \omega /2}}{1-e^{-\beta \hslash \omega }}$

Ved så høye temperaturer at $k_{B}T\gg \hslash \omega ,$ går dette uttrykket over til å gi det klassiske resultatet $Z=k_{B}T/\hslash \omega .$

Oscilllatorens indre energien følger nå som

$U=-\frac{\partial }{\partial \beta }\ln Z=\frac{1}{2}\hslash \omega +\frac{\hslash \omega }{e^{\hslash \omega /k_{B}T}-1}$

Det første bidraget $\hslash \omega /2$ er dens nullpunktsenergi som er den indre energien i grensen T → 0. Ved høyere temperaturer bidar det andre leddet som Einstein benyttet i 1907 til å forklare hvordan varmekapasiteten til et fast stoff varierer med temperaturen. I sin modell tenkte han seg at det består av partikler i faste posisjoner som kan oscillere med samme frekvens om sine likevektsposisjoner.

Degenerasjon og tilstandstetthet

En partikkel har vanligvis mange forskjellige bevegelsestilstander som har samme energi. For eksempel vil en partikkel med kvantisert spinn s, innta i alt 2s + 1 forskjellige retninger i rommet. Hvis dens mulige energier E  er uavhengige av dette spinnet, vil derfor hvert energinivå inneholde et slikt antall mikrotilstander. Det sies derfor å ha en «degenerasjon» g(E) = 2s + 1. Også andre forhold som har å gjøre med symmetrier i systemet, kan bidra til denne.

På lignende vis vil en partikkel i klassisk mekanikk kunne ha en lignende degenerasjon. For eksempel varierer den kinetiske energien med kvadratet av hastigheten til partikkelen. Da den er uavhengig av retningen til hastigheten, vil den i utgangspunktet synes å ha et uendelig stort antall mikrotilstander. Men dette problemet blir løst av kvantemekanikken. På samme måte som for den harmoniske oscillatoren, er det ikke mulig å angi presise verdier for posisjonen r = Mal:Nowrap eller impulsen p = Mal:Nowrap til en partikkel. De kan kun angis innen et infinitesemalt volumelement i det kombinerte rom-impulsrommet som da vil inneholde i alt d³r d³p /h³ mulige kvantetillstander. På dette viset kan man beregne en «tilstandstetthet» g(E) slik at antall mikrotilstander mellom E  og E + dE  kan skrives som g(E)dE.

Antall partikler som finnes med energier mellom E  og E + dE  vil nå kunne skrives på samme vis som n(E)dE. Det må være lik med sannsynligheten p(E) for å finne partikkelen i én mikrotilstand med energi E  multiplisert med antall mikrotilstander med denne energien og det totale antall partikler N  i systemet. På denne måten har man

$\frac{n(E)}{N}=\frac{g(E)}{Z}{e}^{-\beta E}$

hvor nå partisjonsfunksjonen er gitt ved integralet

$Z=\underset{0}{\overset{\mathrm{\infty }}{\int }}dEg(E)e^{-\beta E}$

På samme vis beregnes middelverdien av energien til én partikkel fra

$⟨E⟩=\frac{1}{Z}\underset{0}{\overset{\mathrm{\infty }}{\int }}dEg(E)E{e}^{-\beta E}$

slik at den midlere totalenergien $U=N⟨E⟩$ er systemets indre energi.

Ideell gass

Hver partikkel i en ideell gass har energien E = p²/2m  når den ikke har andre frihetsgrader enn de som skyldes dens egen bevegelse. Da denne er uavhengig av partikkelens posisjon i volumet V, blir dens tilstandstetthet

$g(E)dE=\frac{V}{h^3}{d}^{3}p=\frac{V}{h^3}4\pi p^{2}dp$

Sannsynlighetsfordelingen for impulsene i denne gassen er dermed

$n(E)dE=\frac{N}{Z}\frac{V}{h^3}4\pi p^2{e}^{-p^2/2m{k}_{B}T}dp$

som tilsvarer Maxwells hastighetsfordeling da p = mv.

Indre energi

De termodynamiske egenskapene er inneholdt i partisjonsfunkjonen

$Z=\frac{V}{h^3}4\pi \underset{0}{\overset{\mathrm{\infty }}{\int }}dp{p}^2{e}^{-p^2/2m{k}_{B}T}$

Integralet kan gjøres på enkelt vis ved å definere $a=1/2m{k}_{B}T.$ Da er

$Z=\frac{V}{h^3}4\pi (-\frac{\partial }{\partial a})\underset{0}{\overset{\mathrm{\infty }}{\int }}dp{e}^{-a{p}^2}$

Dette er nå et Gauss-integral med verdien (1/2)√(π /a). Ved en enkel derivasjon finner man så resultatet $Z=V/{\mathrm{\Lambda }}^3$ hvor

$\mathrm{\Lambda }=\sqrt{\frac{h^2}{2\pi m{k}_{B}T}}$

kalles den «termiske bølgelengden» til partikkelen. Middelverdien av dens energi blir

$⟨E⟩=-\frac{\partial }{\partial \beta }\ln Z=\frac{3}{2}{k}_{B}T,$

i overensstemmelse med kinetisk teori. Har partiklene indre frihetsgrader som når de er sammensatt av flere atom, vil disse også bidra til deres indre energi.

Fri energi

I motsetning til kinetisk teori tillater statistisk mekanikk også en beregning av de termodynamiske egenskapene til systemet. De følger fra partisjonsfunksjonen når den skrives på formen

$Z=e^{-F/k_{B}T}$

hvor funksjonen F  er Helmholtz fri energi for systemet. Det følger fra den indre energien

$⟨E⟩=-\frac{\partial }{\partial \beta }\ln Z=F-T\frac{\partial F}{\partial T}$

hvor $S=-\partial F/\partial T$ er den termodynamiske entropien til partikkelen. Men en eksplisitt utregning viser at den ikke blir ekstensiv. Det betyr her at hvis man for eksempel dobler størrelsen både på volumet V  og antall partikler N, så skal også entropien til gassen fordobles. Problemet har sin forklaring i at man ennå ikke har tatt hensyn til at gassen består av like partikler.

Ekstensiv entropi

Partisjonsfunksjonen for hver partikkel i gassen er gitt ved $Z$ slik at den ville vært $Z^N$ hvis de kunne lokaliseres. Men partiklene er like og kan byttes om på $N!=N(N-1)(N-2)\cdots 1$ forskjellige måter som ikke kan adskilles Gassens partisjonsfunksjon er derfor

$Z_N=\frac{Z^N}{N!}=e^{-F_N/k_{B}T}$

hvor nå $F_N$ er den totale, frie energien for alle $N$ partikler. Denne justeringen er den samme som løser Gibbs paradoks for en blanding av to forskjellige gasser.

Da partikkeltallet $N\gg 1,$ kan man benytte Stirlings formel $\ln N!=N\ln N-N$ slik at den fri energien til gassen blir

$\begin{array}{cc}{F}_N& =-k_{B}T\ln \frac{Z^N}{N!}=-N{k}_{B}T(\ln \frac{Z}{N}+1)\\ & =N{k}_{B}T(\ln \rho {\mathrm{\Lambda }}^3-1)\end{array}$

hvor $\rho =N/V$ er tetthet av partikler i gassen. Dette kan skrives som $F_N=U_N-T{S}_N$ hvor $U_N=3N{k}_{B}T/2$ er den totale, indre energien till gassen. Det betyr at den har entropien

$S_N=N{k}_B[\frac{5}{2}-\ln \rho {\mathrm{\Lambda }}^3]$

som nå er proporsjonal med antall partikler den inneholder.

Sammensatte partikler

En ideell gass består av partikler som ikke vekselvirker med hverandre. Vanligvis er de molekyl som er sammensatt av flere atom. Energien til et slikt molekyl kan beregnes ved bruk av kvantemekanikk som gir diskrete verdier E₀ ≤ E₁ ≤ E₂ og så videre. Den laveste energien E₀ opptrer når molekylet er i «grunntilstanden», mens E₁ er energien til den første eksiterte tilstanden over denne.

Boltzmann-fordelingen sier at sannnsynlighten for å finne et molekyl i en tilstand med energi E_n, er gitt som Mal:Nowrap når gassen er i termisk likevekt med temperaturT. Hvis denne senkes, vil derfor flere og flere partikler befinne seg i de laveste energitilstandene. Til slutt når T → 0, vil de alle ende opp i grunntilstanden. Dette er et typisk, kvantemekanisk fenomen som eksperimentelt kan påvises ved måling av varmekapasiteten til gassen.

Den indre energien til hver partikkel kan nå beregnes som

$U=⟨E⟩=\frac{1}{Z}\sum _{n=0}^{\mathrm{\infty }}{E}_n{e}^{-\beta E_n}$

hvis man ser bort fra en mulig degenerasjon i hvert energinivå. Når temperaturen går mot null, er kun de første leddene i de to summene av betydning. Da har man

$\begin{array}{cc}U& =\frac{E_0{e}^{-E_0/k_{B}T}+E_1{e}^{-E_1/k_{B}T}+\cdots }{e^{-E_0/k_{B}T}+e^{-E_1/k_{B}T}+\cdots }\\ & =E_0+E_c{e}^{-E_c/k_{B}T}+\cdots \end{array}$

hvor energiforskjellen E_c = E₁ - E₀. Det betyr at under en «karakteristisk temperatur» Mal:Nowrap vil det første leddet dominere i uttrykket for den indre energien. Da dette er uavhengig av temperaturen, vil det tilsvarende bidraget til gassens varmekapasitet bli null. Det kommer tydelig frem fra

$c_V={\left(\frac{\partial U}{\partial T}\right)}_V=k_B{\left(\frac{T_c}{T}\right)}^2{e}^{-T_c/T}$

som blir forsvinnende lite i denne grensen.

Energinivåene til partikler i et makroskopisk volum kan beregnes på tilsvarende måte som for en partikkel i en kubisk boks ved å løse Schrödinger-ligningen. De resulterende verdier ligger da meget nær hverandre slik at den karakteristiske temperaturen for bevegelsesenergien i praksis kan settes lik null. Dermed kan summasjonen over disse diskrete verdiene erstattes med en integrasjon og dermed gi det klassiske resultatet Mal:Nowrap. Dette er derimot ikke tilfellet for bidragene til den indre energien fra rotasjon og vibrasjon av molekylet, og de kvantemekaniske verdiene for energinivåene vil gjøre en forskjell ved tilstrekkelig lave temperaturer.

Diatomisk molekyl

De enkleste molekylene består kun av to atomer og sies derfor å være diatomiske. Typiske eksempel er nitrogenmolekylet N₂ og karbonmonoksid CO eller "kullos". Avstanden mellom de to atomene er gitt ved en balansering av forskjellige krefter som virker på dem. Men molekylet kan også befinne seg i eksisterte tilstander som tilsvarer svingninger rundt denne likevektsavstanden. De kan beskrives som en harmonisk oscillator som har kvantiserte energier. Energiforskjellen mellom grunntilstanden og første, eksiterte tilstand er Mal:Nowrap der ω  er oscillatorens vinkelfrekvens. Den resulterende, karakteristiske temperatur T_c  er for de aller fleste molekyl mye høyere enn typisk romtemperatur. Derfor vil ikke denne indre frihetsgraden bidra til den spesifikke varmekapasiteten ved vanlige temperaturer.

Klassisk rotasjon

Viktigere er bidraget fra rotasjon av molekylet. Denne kan foregå langs to akser som står vinkelrett på hverandre og på aksen mellom de to atomene. Grunnen er at rotasjoner om denne har et veldig lite treghetsmoment slik at den resulterende rotasjonsenergeni blir neglisjerbar. En klassisk beskrivelse av molekylet sier da at det har to kvadratiske frihetsgrader som derfor bidrar tilsammen med k_B  til den spesifikke varmekapasiteten.

Dette er i overensstemmelse med hva som følger av en direkte beregning av partisjonsfunksjonen for denne bevegelsen. Hvis orienteringen til molekylet er beskrevet ved kulekoordinater (θ, φ), er dets energi gitt som

$E=\frac{1}{2m{r}^2}(p_{\theta }^2+\frac{p_{\varphi }^2}{{\sin }^2\theta })$

hvor m  er den reduserte massen til de to atomene og r  er deres likevektavstand. Dermed er Mal:Nowrap molekylets treghetsmoment. Partisjonsfunksjonen er nå gitt ved integralet

$Z=\frac{1}{h^2}\int e^{-\beta E}d\varphi d{p}_{\varphi }d\theta d{p}_{\theta }=\frac{8{\pi }^{2}I}{\beta h^2}$

Bidraget til den indre energien fra dette ene molekylet er derfor

$U=-\frac{\partial }{\partial \beta }\ln Z=k_{B}T$

Det gir den spesifikke varmen c_V = k_B  som forventet i denne klassiske grensen.

Kvantemekanisk rotasjon

Energien til en roterende partikkel er i kvantemekanikken gitt ved egenverdiene til kvantisert dreieimpuls. For et lineært molekyl er disse gitt som

$E_{\ell }=\frac{{\hslash }^2}{2I}\ell (\ell +1)$

hvor I  er dets treghetsmoment og kvantetallet ℓ = 0, 1, 2 og så videre. Hvert av disse energinivåene har en degenerasjon Mal:Nowrap som tilsvarer de forskjellige retningene den kvantiserte dreieimpulsen kan ha i rommet.

Den kvantemekaniske partisjonspartisjonen til et molekyl er dermed

$Z=\sum _{\ell =0}^{\mathrm{\infty }}(2\ell +1)e^{-\ell (\ell +1)T_c/2T}$

når man definerer en karakteristisk temperatur ved k_BT_c = ħ²/I. Denne uendelige summen kan ikke gjøres eksakt. Men ved høye temperaturer forandrer hvert ledd seg lite når ℓ øker med 1. Man kan da med god tilnærmelse erstatte summen med et integral. Det kan lett beregnes ved å innføre den nye variable Mal:Nowrap slik at Mal:Nowrap Dermed blir

$\begin{array}{cc}Z& =\underset{0}{\overset{\mathrm{\infty }}{\int }}d\ell (2\ell +1)e^{-\ell (\ell +1)T_c/2T}\\ & =\frac{2T}{T_c},(T\gg T_c)\end{array}$

i overensstemmelse med det klassiske resultatet. Derimot for temperaturer vel under den karakterristiske temperaturen finner man nå et helt annet resultat. De to første leddene dominerer da summen og man har

$Z=1+3e^{-T_c/T},(T\ll T_c)$

Bidraget til den indre energien blir

$U=-\frac{\partial }{\partial \beta }\ln Z=3k_B{T}_c{e}^{-T_c/T}$

som igjen gir et sterkt redusert bidrag til den spesifikke varmekapasiteten sammenlignet med den klassiske verdien.

Den karakteristiske temperaturen som skyldes rotasjon, er for de fleste diatomiske molekyl så lav at den klassiske verdien for deres spesifikke varmekapasiteter som følger fra ekvipartisjonsprinsippet, er gyldig i et vidt område over og under vanlige romtemperaturer. Men når ett av atomene er hydrogen som i HCl, er dette så lett at det resulterende treghegtsmomentet I  til molekylet blir tilsvarende mindre. Dermed blir den karakteristiske temperaturen høyere og kvanteeffektene mer synlige. Dette gjelder spesielt for det letteste molekylet, H₂ hvor også andre kvanteeffekter spiller en viktig rolle.

Hydrogenmolekylet

Et "heteronukleært" molekyl som CO må roteres 360° om en av sine to symmetriakser for å komme tilbake til seg selv. Derimot behøves det bare en rotasjon på 180° for at det skal skje med et "homonukleært" molekyl som N₂. Dette vil halvere den klassiske partisjonsfunksjonen for et slikt molekyl.

Kvantemekanisk representerer en rotasjon på 180° av et homonukleært molekyl et ombytte av to identiske atomkjerner. Hvis de har heiltallig spinn, er de bosoner, og den komplette bølgefunksjonen til molekylet må være symmetrisk under denne operasjonen. Alternativt er de fermioner med halvtallig spinn og bølgefunksjonen må være antisymmetrisk. Dette er tilfelle for hydrogenmolekylet H₂ hvor atomkjernene er protoner som har spinn Mal:Nowrap. Da dette molekylet også har en høy, karakteristisk temperatur, vil denne kvanteeffekten komme til syne i den spesifikke varmen til hydrogengass.

Det totale spinnet til to partikler, hver med spinn Mal:Nowrap, er enten en triplett med Mal:Nowrap elller en singlett med Mal:Nowrap Da bølgefunksjonene til de tre triplett-tilstandene er symmetriske, må bølgefunksjonen til den sammenkoblede rotasjonsbevegelsen være antisymmetrisk for at den totale bølgefunksjonen skal være antisymmetrisk. Det betyr at de har rotasjonskvantetall Mal:Nowrap, 3, 5 og så videre. I motsatt fall for singletter hvis bølgefunksjon er antisymmetrisk, må man på samme vis ha Mal:Nowrap, 2, 4 og så videre. Den totale partisjonsfunksjonen for hydrogenmolekylet må derfor splittes opp i to ulike summer tilsvarende bidragene fra singletter og tripletter. Disse to delene omtales som «parahydrogen» og «ortohydrogen». Den laveste rotasjonsenergien Mal:Nowrap kan bare parahydrogen ha. Det betyr at ved temperatur Mal:Nowrap vil hele hydrogengassen bestå av slike singlettmolekyler når den befinner seg i termisk likevekt.

Denne forklaringen av den spesifikke varmen til hydrogen ble først gitt av Werner Heisenberg i 1927. Den var basert på Paulis eksklusjonsprinsipp som han også benyttet til å gi en tilsvarende forklaring av energinivåene til heliumatomet. Han viste hvordan de splittes opp i kvantetilstander som enten er singletter eller tripletter som i dette tilfellet dannes av de to elektronene i atomet. I 1929 ble de to formene til hydrogenmolekylet eksperimentelt påvist. Dette ble gitt som en del av begrunnelsen da Heisenberg ble tildelt nobelprisen i fysikk i 1932.

Reelle gasser

I motsetning til den ideelle gassen vil partiklene i en reell gass vekselvirke seg i mellom. Inneholder den N  partikler og vekselvirkningen antas å være uavhengig av deres hastigheter, kan den beskrives ved et potensial $U({𝐫}_1,{𝐫}_2,\dots ,{𝐫}_N).$ Energien til hele gassen er dermed

$E(𝐫,𝐩)=\sum _{a=1}^N\frac{{𝐩}_a^2}{2m}+U({𝐫}_1,{𝐫}_2,\dots {𝐫}_N)$

hvor en klassisk mikrotilstand til gassen er gitt ved å spesifisere alle posisjonene $𝐫=({𝐫}_1,{𝐫}_2,\dots ,{𝐫}_N)$ i et volum V  sammen med de tilsvarende impulsene $𝐩=({𝐩}_1,{𝐩}_2,\dots ,{𝐩}_N).$ Sannsynligheten for å finne gassen i denne tilstanden er nå gitt ved Bolltzmann-fordelingen

$p(𝐫,𝐩)=\frac{1}{Z_N}{e}^{-\beta E(𝐫,𝐩)}$

som er den samme som i det kanoniske ensemblet. Her er nå Z_N  partisjonsfunksjonen til hele gassen og er gitt ved integralet

$\begin{array}{cc}{Z}_N& =\frac{1}{N!}\frac{1}{h^{3N}}\int d^{3N}r{d}^{3N}p{e}^{-E(𝐫,𝐩)/k_{B}T}\\ & =e^{-F_N/k_{B}T}\end{array}$

hvor det er tatt hensyn til at partiklene er like. Funksjonen $F_N=F_N(T,V)$ er da Helmholtz fri energi til denne reelle gassen.

Konfigurasjonsintegral

Da den kinetiske energien til partiklene er kvadratisk i deres impulser, kan den delen av integralet i partisjonsfunksjonen gjøres som for den ideelle gassen. Dermed blir

$Z_N=\frac{1}{N!}\frac{Q_N}{{\mathrm{\Lambda }}^{3N}}$

hvor $\mathrm{\Lambda }$ er den termiske bølgelengden og

$Q_N(T,V)=\int d^{3N}r{e}^{-\beta U({𝐫}_1,{𝐫}_2,\dots ,{𝐫}_N)}$

kalles «konfigurasjonsintegralet» til gassen. Det inneholder effekten av alle vekselvirkningene mellom partiklene. For en ideell gass er det ganske enkelt $V^N.$ Da trykket i gassen er gitt ved den termodynamiske deriverte $P=-(\partial F/\partial V{)}_T,$ vil det nå kunne beregnes fra

$P=k_{B}T\frac{\partial }{\partial V}\ln Q_N(T,V)$

På denne måten kan tilstandsligningen for den reelle gassen finnes.

Eksterne lenker

• L. Susskind, Statistical Mechanics Lectures, YouTube from Stanford University (2011)