Maxwell-Boltzmann-verdeling

Bestand:Maxwell-Distr.png

De Maxwell-Bolzmann-verdeling voor stikstof (N₂) bij drie verschillende temperaturen.

De Maxwell-Boltzmann-verdeling of snelheidsverdelingswet van Maxwell-Boltzmann geeft de verdeling van de snelheden van gasmoleculen in een ideaal gas weer, wanneer de moleculen als puntvormig kunnen worden opgevat en zij volkomen elastisch botsen, zodat impuls en energie behouden blijven. Er vinden tevens geen simultane botsingen plaats van 3 of meer moleculen. De Maxwell-Boltzmann-verdeling vervult een centrale rol in toepassingen van de kinetische gastheorie.

De dichtheid $f(v)$ van de snelheidsverdeling van de deeltjes wordt gegeven door:

$f(v)=4\pi {\left(\frac{m}{2\pi kT}\right)}^{\frac{3}{2}}{v}^2{e}^{-\frac{m{v}^2}{2kT}}$

Daarin is

• $m$ de massa van een deeltje van het gas (in SI-eenheden in kg)
• $k$ de Boltzmannconstante (in SI-eenheden geldt $k\approx 1,38\times 10^{-23}$ J K⁻¹)
• $v$ de snelheid van een deeltje (in SI-eenheden in m s⁻¹)
• $T$ de absolute temperatuur van het gas (in SI-eenheden in K)

De verdeling is genoemd naar James Clerk Maxwell, die haar als eerste in 1866 afleidde, en Ludwig Boltzmann, die het bewijs verscherpt heeft. De verdeling is een bijzonder geval van de algemene Boltzmann-verdeling.

Afleiding

In de stationaire toestand zijn de gasdeeltjes gelijkmatig verdeeld over het volume. De energie van een deeltje is zijn kinetische energie en omdat de totale energie $E$ van het gas vastligt, is de snelheid van een deeltje begrensd. Alle mogelijke snelheden worden opgedeeld in een eindig aantal ($m$) klassen, waarbinnen de snelheid weinig varieert.

Voorwaarden

Elk van de $N$ deeltjes valt wat zijn snelheid betreft binnen een van de $m$ klassen. De aantallen in de klassen zijn $n_1,n_2,\dots ,n_m$. Er geldt dus:

$\sum n_i=N$

Ook moet het totaal van de energie van de deeltjes gelijk zijn aan de totale energie $E$ van het gas, dus:

$\sum n_i{E}_i=E$

Waarschijnlijkheidsdichtheden van realisaties

De verdeling van de deeltjes over de snelheidsklassen kan op meer manieren gerealiseerd worden. Zijn alle deeltjes in één klasse dan is er maar één manier, maar zijn ze op een na alle in één klasse dan zijn er al $N$ mogelijke realisaties. Algemeen is het aantal realisaties bij de verdeling van $n_1,n_2,\dots ,n_m$ deeltjes over de $m$ klassen:

$A=\frac{N!}{n_1!n_2!\dots n_m!}$

Hoe meer realisaties een verdeling heeft, hoe waarschijnlijker het is dat het gas zich in een realisatie van die verdeling bevindt, uitgaande van het belangrijkste postulaat van de statistische mechanica, namelijk dat alle microtoestanden a priori gelijke waarschijnlijkheden hebben. De meest waarschijnlijke verdeling is dus de verdeling met het grootste aantal realisaties, zij het dat aan de genoemde voorwaarden moet zijn voldaan.

Optimalisatie

Onder deze voorwaarden wordt de verdeling bepaald waarvoor het aantal realisaties $A$ maximaal is. Om gemakkelijker te rekenen neemt men, in plaats van het aantal realisaties $A$ zelf, de logaritme daarvan. Dit is toegestaan omdat de logaritme monotoon stijgend is. Met de multiplicatorenmethode van Lagrange wordt de vergelijking:

$\frac{\partial \log A}{\partial n_k}+\frac{\partial }{\partial n_k}(a\sum n_i+b\sum n_i{E}_i)=0$

waarin $a$ en $b$ de multiplicatoren zijn.

Uitwerken levert:

$\frac{\partial \log (n_k!)}{\partial n_k}+a+b{E}_k=0$

Met behulp van de formule van Stirling krijgt men de benadering:

$\log n!\approx n\log n-n$,

zodat voor de vergelijking resulteert:

$\log n_k+a+b{E}_k=0$,

met als oplossing:

$n_k=e^{-a-b{E}_k}=B{e}^{-b{E}_k}$

Aangezien de energie in een klasse alleen de kinetische energie van een deeltje in die klasse is, geldt:

$E_k=\frac{1}{2}m{v}_k^2$

De snelheidsverdeling heeft dus de dichtheid:

$n(\overrightarrow{v})=B{e}^{-\frac{1}{2}bm{v}^2}$

Deze is alleen afhankelijk van de grootte van de snelheid.

Toestandsdichtheid

Voor de toestandsdichtheid als functie van de snelheid geldt in drie dimensies:

$D_3(v)=4\pi v^2$

De uitdrukkingen voor de hypothetische toestandsdichtheid in 2 dimensies en 1 dimensie zijn (voor de volledigheid):

$D_2(v)=2\pi v$

$D_1(v)=1$

De kanssdichtheid in $n$ dimensies wordt gegeven door:

$f(v)=D_n(v)B{e}^{-\frac{1}{2}bm{v}^2}$

Kansdichtheid

Omdat de integraal van de kansdichtheid in drie dimensies gelijk moet zijn aan 1, volgt voor de constante $B:$

$1=\int f(v)\mathrm{d}v=4\pi B\underset{0}{\overset{\mathrm{\infty }}{\int }}{v}^2{e}^{-\frac{1}{2}bm{v}^2}\mathrm{d}v=4\pi B\frac{1}{(bm{)}^{\frac{3}{2}}}\underset{0}{\overset{\mathrm{\infty }}{\int }}{z}^2{e}^{-\frac{1}{2}{z}^2}\mathrm{d}z=4\pi B\frac{\sqrt{\frac{\pi }{2}}}{(bm{)}^{\frac{3}{2}}}$

De dichtheid wordt dus:

$f(v)=\sqrt{\frac{2}{\pi }}(bm{)}^{\frac{3}{2}}{v}^2{e}^{-\frac{1}{2}bm{v}^2}$

Voor een ideaal gas geldt voor de verwachte kinetische energie van een deeltje:

$⟨\frac{1}{2}m{v}^2⟩=\frac{3}{2}kT$,

dus

$\frac{3}{2}kT=\frac{1}{2}m\underset{0}{\overset{\mathrm{\infty }}{\int }}{v}^{2}f(v)\mathrm{d}v=\frac{1}{2}m\underset{0}{\overset{\mathrm{\infty }}{\int }}\sqrt{\frac{2}{\pi }}(bm{)}^{\frac{3}{2}}{v}^4{e}^{-\frac{1}{2}bm{v}^2}\mathrm{d}v=\frac{1}{2b}\sqrt{\frac{2}{\pi }}\underset{0}{\overset{\mathrm{\infty }}{\int }}{z}^4{e}^{-\frac{1}{2}{z}^2}\mathrm{d}z=\frac{3}{b}$

Dus is:

$b=\frac{2}{kT}$

en de kanssdichtheid in drie dimensies:

$f(v)=4\pi {\left(\frac{m}{2\pi kT}\right)}^{\frac{3}{2}}{v}^2{e}^{-\frac{m{v}^2}{2kT}}$

Compacte representatie

Bestand:Maxwell-Boltzmann distribution pdf.svg

Kansdichtheid van de Maxwell-Boltzmann-verdeling

Bestand:Maxwell-Boltzmann distribution cdf.svg

Verdelingsfunctie van de Maxwell-Boltzmann-verdeling

De kanssdichtheid van de Maxwell-Boltzmann-verdeling wordt in een compacte representatie gegeven door:

$f(v)=\sqrt{\frac{2}{\pi }}\frac{v^2{e}^{-v^2/\left(2a^2\right)}}{a^3}$

met:

$a=\sqrt{\frac{kT}{m}}$

De verdelingsfunctie is:

$F(v)=\mathrm{erf}\left(\frac{v}{\sqrt{2}a}\right)-\sqrt{\frac{2}{\pi }}\frac{v{e}^{-v^2/\left(2a^2\right)}}{a}$

Verschillende snelheidsgemiddelden

De meest waarschijnlijke (modale) snelheid ligt bij het maximum van $f(v)$ waarvoor geldt:

$\frac{\mathrm{d}}{\mathrm{d}v}f(v)=0$

en wordt gegeven door:

$v_p=a\sqrt{2}=\sqrt{\frac{2kT}{m}}$

RMS-snelheid

De kwadratisch gemiddelde (rms) snelheid is:

$v_{\text{rms}}=\sqrt{\overline{v^2}}=\sqrt{\underset{0}{\overset{\mathrm{\infty }}{\int }}{v}^{2}f(v)\mathrm{d}v}$

en wordt gegeven door:

$v_{\text{rms}}=a\sqrt{3}=v_p\sqrt{\frac{3}{2}}=\sqrt{\frac{3kT}{m}}$

Gemiddelde snelheid

De gemiddelde snelheid is:

$v_{\text{av}}=\bar{v}=\underset{0}{\overset{\mathrm{\infty }}{\int }}vf(v)\mathrm{d}v$

en wordt gegeven door:

$v_{\text{av}}=2a\sqrt{\frac{2}{\pi }}=2\frac{v_p}{\sqrt{\pi }}=\sqrt{\frac{8kT}{\pi m}}$

Verhoudingen

De drie waarden verhouden zich onderling bij benadering als:

$v_{rms}:v_{\text{av}}:v_p\approx 1,00:0,92:0,82$

De gemiddelde snelheid ligt ongeveer midden tussen de rms-snelheid en de meest waarschijnlijke snelheid in.

Externe links

• afleiding + LBRoulette
• afleiding
• verdeling