7.2 Geschwindigkeitsverteilung

PIC

7.2.1 Maxwell-Ansatz

Die Wahrscheinlichkeitsdichte f(v) = f(vx,vy,vz) gibt die Wahrscheinlichkeit an, ein Teilchen pro Einheitsvolumen in Geschwindigkeitsraum [vx,vx + dvx], [vy,vy + dvy], [vz,vz + dvz] zu finden. Die Wahrscheinlichkeit erhält man dann durch Multiplikation der Wahrscheinlichkeitsdichte mit dem Volumenelement dvxdvydvz:
f(v)dv = f(vx,vy,vz)dvxdvydvz

Schließlich läßt sich die Gesamtwahrscheinlichkeit durch Integration über alle möglichen Geschwindigkeiten finden; diese muß gerade den Wert 1 haben:

      integral  oo    integral  oo    integral  oo                         integral  oo 
1 =                  f(vx,vy,vz)dvxdvy dvz =     f(v)dv
   vx=- oo  vy=- oo  vz= - oo                     v= - oo

PIC
Die Maxwell-Annahme ist nun, daß f(vx), f(vy) und f(vz) nicht gekoppelt sind, da das Gas isotrop ist. Diese lassen sich also separat betrachten. Wir nehmen an, daß sich f(v) somit als Produkt dreier unabhängiger Funktionen bezüglich einer einzigen Geschwindigkeitskomponente schreiben läßt:

f(v) = f(vx,vy,vz) = f(vx).f(vy).f(vz)

                               V~ -----------
v.v = |v|2 = v2x + v2y +v2z,v = |v|= v2x +v2y + v2z

Dann wenden wir den Logarithmus auf obige Gleichung an:

lnf(v) = lnf(vx)+ lnf(vy)+ ln f(vz)

Wir leiten die Funktion nach vx ab und erhalten:

(      )
  @ln(v)-     = d-lnf-(vx) + dln(vy)-
   @vx   vy,vz     @vx     --@x--
                            =0

(        )              (    )
  @lnf(v)      = dln(v).  dv-
    @vx    vy,vz     dv     dvx vy,vz

Mit v =  V~ -----------
  v2x + v2y + v2z ergibt sich nun:

(    )
  -dv      = 1  V~ --2vx-----=  vx
  dvx  vy,vz   2  v2x + v2y + v2z  v

(       )              (    )
 @-lnf-(v)        d-ln(v)   dv-        dlnf(v)  vx
   @vx    vy,vz =  dv  .  dvx vy,vz =   dv   . v

Wir multiplizieren das ganze mit v1x durch:

dlnf(v)=  dlnf(vx)= d-ln-f(vy) = d-lnf-(vz) != const.  /\ = - 2g
  vdv      vxdvx      vydvy      vzdvz

|---------------|
|d-ln-f(vi)       |
|  vidvi  = -2g |
----------------

Durch Integration dieser Differentialgleichung erhalten wir:

|-------------------|
|f(vi) = A .exp(- gv2)
------------------x--

Des weiteren werten wir die Normierungsbedingung aus:

   oo  integral                     integral  oo                     o integral  o 
      f(v )dv = 1 = A .     exp(- gv2)dv = A .2.   exp (-gv2)dv
        i   i                      i   i                 i    i
vi=- oo                 vi= - oo                      0

     V~ ---
2A .  -p-= 1
      4g

|---------------------|
|       V~ -g     (    )|
f (vx) =   -.exp - gv2x |
----------p------------

Dies funktioniert auch analog für vx und vz.

7.2.2 Druck in der kinetischen Gastheorie

Wir nehmen an, daß die Komponenten vy und vz gleich 0 sind und daß sich außerdem die Geschwindigkeit vx nach einem Stoß in -vx ändert.

PIC

Der übertragene Impuls berechnet sich dann durch Dp = D(vx) = 2mvx. Wir berechnen außerdem die Kraft als zeitliche Impulsänderung:

D(mvx)-  2mvx-   mv2x-
  Dt   =   2va  =  a  = F
            x

Wir dividieren F durch die Fläche b . c, womit der Druck resultiert:

    -F--      -mv2x-
p = b.c <==> p = a .b.c

Der Druck, den ein Teilchen i erzeugt, lautet pi =   2
mvVx. Durch Summation über alle Teilchen N erhalten wir den Gesamtdruck:

    sum N     sum N mv2ix   m  sum N 2
p =   pi =    -V---= V-   vix
   i=1     i=1          i=1

Mit der mittleren quadratischen Geschwindigkeit

         N sum 
<v2x> =-1    v2ix
      N  i=1

läßt sich dies schreiben als:

         <  >             <  >
p = m-.N  v2ix  , pV = m .N . v2ix
    V

Es gilt, da es sich um ein isotropes Gas handelt und somit der Betrag der Geschwindigkeit unabhängig von der Richtung ist:

 2   2   2    2
v = vx + vy + vz

<  >  <  >   < >   1<  >
 v2x =  v2y  = v2z  = 3 v2

pV = 1m  .N  .<v2>
     3

Das ideale Gasgesetz lautet, wie wir wissen pV = nRT;

       1  N      <  >   1          < >
nRT =  -m ---.NA  v2 =  -m .n.NA  .v2
       3  NA            3

Damit gilt für die mittlere quadratische Geschwindigkeit, wobei M die Molmasse ist:

< 2>   -3RT--   3RT-
 v   = m .NA =   M

Allgemein gilt für den Mittelwert eine Funktion g(x):

         integral  oo 
<g(x)> =   g(x).f(x)dx
       -  oo

Wir wenden dies auf unser Problem an:

        oo                    V~ --
< 2>    integral  2                  g-   (   2)
 vx =    vxf(vx)dvx,f(vx) =  p exp -gvx
      - oo

       integral  oo   V~                     integral  oo   V~                 V~ -    V~ --
RT- =      g-.exp(- gv2).v2.dvx = 2.   gv2.exp (-gv2)=  2  g-1-. p-= -1-
 M         p         x  x            p x        x       p4g   g   2g
      - oo                         0

Damit erhalten wir die eindimensionale MAXWELL-BOLTZMANN-Geschwindigkeitsverteilung:

|------- V~ ----------(------)--- V~ ----------(------)-|
|         -M---        M-v2x      --m--        mv2x  |
|f(vx) =  2pRT  .exp  - 2RT   =   2pRT  .exp - 2RT   |
----------------------------------------------------

PIC

7.2.3 Dreidimensionale Geschwindigkeitsverteilung

Erinnern wir uns daran, daß man f(v) als Produkt aus f(vx), f(vy) und f(vz) schreiben kann:
                                     (  m  )32    (   m  (           ))
f(v) = f(vx) .f(vy).f (vz) = f (vx,vy,vz) = -----   .exp - ---- v2x + vy2+ v2z  =
      |-------3-----(------)-|        2pRT           2RT
    = ( -m---)2 .exp  - mv2-  |
      --2pRT-----------2RT----
(7.1)
Für uns ist die Wahrscheinlichkeitsdichte von Teilchen mit einem Betrag der Geschwindigkeit zwischen v und v + dv, die unabhängig von der Richtung ist, interessant:
  oo  integral           integral  oo   integral  oo   integral  oo 
   F (v)dv =          f(vx,vy,vz)dvx dvydvz = 1
v=0         - oo  - oo  - oo

Zur Berechnung des Integrals führen wir Kugelkoordinaten ein:

PIC

x = r.sinh cosf, y = r .sinh sinf, z = rcosf

 oo  integral           integral  oo   integral p  integral p
  F (v)dv =            f(vx,vy,vz)v2 dvsin hdhdf

0          v=0 h=0 f=-p

Die Lösung ist:

|----------(-----)3----(-----2)-|
|F(v) = 4pv2 -m--- 2 exp -mv--  |
-------------2pkT----------2kt---

7.2.4 Temperaturabhängigkeit der Geschwindigkeitskonstanten
(Arrhenius-Gesetz)

Die Geschwindigkeitskonstante einer chemischen Reaktion hängt nach dem Arrhenius-Gesetz exponentiell von der Temperatur ab:
|-----------(-----)-|
|              Ea-  |
|k = A(T )exp - RT  |
---------------------

Exponentialfunktion  /\ = Harmonisches Ensemble

kbimolekular  oc  ps2 |vrel|
 max          AB

                   (      )
       ( -m---)12       mv2x-
f (vx) =  2pkT   exp  - 2kT

7.2.5 Zusammenfassung der Ergebnisse: Maxwell-Boltzmann-Ge schwindigkeitsverteilung

7.2.6 Mittelwerte der dreidimensionalen Maxwell Boltzmann Geschwindigkeitsverteilung

Beispiele:
Es sei T = 273K:



Gas <v> [m]
 s




H2 1692
He 1204
N2 454
Hg 170


     [(       )   (      )   (       ) ]1
|vk|=   vAx -vBx 2 + vAy - vBy 2-  vAz - vBz 2 2

Es gilt kbimolekular  oc |vrel|. |vrel| kann durch folgendes zunächst furchterregendes Integral berechnet werden:

----------------------------------------------------------------------
|        oo  integral   oo  integral   oo  integral   oo  integral   oo  integral   integral  oo                                           |
||v  | =                   |v |f(vAvAvA)f(vBvBvB )dvAdvB dvAdvB dvAdvB |
| rel                      R   x  y z    x y z    x  x   y  y   z  z |
-------- oo -- oo -- oo -- oo -- oo -- oo ----------------------------------------

Die genaue Lösung findet sich in MOORE/HUMMEL: Physical Chemistry.

|---------------|
|      (     )12|
||vrel|=   8kT--  |
---------pmAB----

mAB ist hierbei die reduzierte Masse der beiden Teilchen:

                 |-------------------(-----)1-|
      -mAmB---   |2.Ordnung         2   3kT-- 2|
mAB = mA + mB ==>  kmax     = NAV ps AB  pmAB   |
                 ------------------------------

Wir haben jetzt in unserer Formel eine T-Abhängigkeit aber noch kein exp(-Ea )
  kT.