Langevin-féle leírás

Langevin-féle leírása (1907 vagy 1908)

Langevin levezetése "fizikaibb", Einsteinnél nem lehet tudni, hogy a részecskével konkrétan mi történik, a Langevin-féle leírás ebből a szempontból szemléletesebb.

Langevin: részecskék lökdösik a megfigyelt részecskét, a mozgás a hőmérsékleti fluktuációból jön, tehát egy véletlen erő hozza létre. $m\frac{d^2 x}{dt^2} = -6\pi \eta a\frac{dx}{dt} + F_{veletlen}$

Ahol a bal oldal első tagja a fékezőerő ($\eta$ a viszkozitás (angolul, magyarul)).

Célunk, hogy a fenti egyenletből kiindulva megkapjuk az Einstein-féle leírásnál eredményül kapott $<x^2> = 2Dt$ összefüggést.

$m \frac{d^2x}{dt^2} = -6\pi\eta a\frac{dx}{dt} + F_{veletlen}$

Mindkét oldalt szorozzuk x-szel :

$\frac{1}{2}m\frac{d^2x^2}{dt^2} - m(\frac{dx}{dt})^2 = -3\pi \eta a \frac{dx^2}{dt} + xF_{veletlen}$

$\frac{1}{2}m\frac{d^2<x^2>}{dt^2} - k_B T = -3\pi \eta a \frac{d<x^2>}{dt} + <xF_{veletlen}>$

Vezessük be a $\frac{d<x^2>}{dt} = z$ új változót. Ekkor a fenti egyenlet a következőképpen módosul:

$\frac{dz}{dt} - \frac{2k_B T}{m} = -\frac{6\pi \eta a}{m}z$

A fenti egyenlet megoldása a homogén megoldás és egy partikuláris megoldás összege:

$z_{homogen} = c_1 e^{-\frac{6\pi \eta a}{m}t}$

$z_{partikularis} = \frac{2k_B T}{6\pi \eta a}$ --> z-t konstansnak tételeztük fel

Ekkor z eredeti értékét visszaírva és a két megoldást összeadva a következő eredményt kapjuk:

$\frac{d<x^2>}{dt} = \frac{k_B T}{3\pi \eta a} + c_1 e^{-\frac{6\pi \eta a}{m}t}$

Ám $\frac{6\pi \eta a}{m}$ $\eta, a, m$ konkrét értékeire és $t = 10^{-2}-10^{-3}$-ra ez egy nagyon nagy szám, emiatt $e^{-...}$ nagyon kicsi.

Például: $\eta_{viz}=10^{-3} \frac{kg}{ms}; m = 10^3 \frac{kg}{m^3}10^{-18}m^3 = 10^{-15} kg; a = 10^{-6} m$

Ekkor: $\frac{6\pi\eta a}{m}=\frac{20\centerdot10^{-3}\frac{kg}{ms} 10^6 m}{10^{-15} kg} = 2\centerdot10^7 \frac{1}{s}$ $c_1 e^{-2\centerdot10^7\frac{t}{s}}$ bármely rezonábilis t-re nullát ad, tehát elhagyhatom. Ekkor:

$<x^2> = \frac{k_B T}{3\pi\eta a}t = 2Dt$

$D = \frac{k_B T}{6\pi\eta a} = \frac{1.4\centerdot10^{-23}\frac{J}{K} 300 K}{20\centerdot10^{-3}\frac{kg}{ms}10^{-6}m} = 2.1\centerdot10^{-13} \frac{m^2}{a}$

A fenti képlet alapján becslés (úgy általában): $D_{molekulak} \approx 2\centerdot10^{-10} \frac{m^2}{s}$

---