W szklance wody jest ok. 1025 cząsteczek, a każda z nich porusza się z prędkością kilkuset metrów na sekundę. Nie da się śledzić ich wszystkich, ale nie trzeba: ciśnienie, temperatura i ciepło właściwe okazują się średnimi, a rozkład prędkości wynika z jednego prawa Boltzmanna. Ten dział prowadzi od cząsteczki odbijającej się od ściany do ruchów Browna, dzięki którym policzono atomy.
Czas
ok. 4 h teorii + ok. 5 h zadań
Zadania
25 zadań w treści + 8 feynmanowskich, z odpowiedziami
Źródło
Feynmana wykłady z fizyki, t. I, rozdz. 39–43
0 Prolog: pyłki, które tańczą
W 1827 r. botanik Robert Brown oglądał pod mikroskopem drobinki wydostające się z ziaren pyłku zawieszone w wodzie. Drobinki bez przerwy drgały i błądziły, choć woda była spokojna. Brown sprawdził, że tak samo zachowują się drobinki szkła i minerałów, więc nie chodziło o życie. Przez prawie osiemdziesiąt lat zjawisko pozostawało ciekawostką.
W 1905 r. Einstein pokazał, że to chaotyczne uderzenia cząsteczek wody, i przewidział, jak szybko drobinka powinna się oddalać od punktu startu. Wynik zależał od liczby cząsteczek w molu. W latach 1908–1909 Jean Perrin zmierzył te ruchy i wyznaczył liczbę Avogadra z dokładnością kilkunastu procent. Wielu fizyków i chemików wciąż wątpiło wtedy, czy atomy istnieją naprawdę. Po pomiarach Perrina spór praktycznie się skończył; w 1926 r. dostał za nie Nagrodę Nobla. Część 5 i przykład 8.6 powtarzają jego rachunek.
Feynman zaczyna wykład o materii od twierdzenia, że gdyby cała wiedza miała zginąć i można było przekazać tylko jedno zdanie, powinno ono mówić, że wszystko składa się z atomów w nieustannym ruchu. W tym dziale wyprowadzimy z tego zdania prawa gazu doskonałego, ciepło właściwe, spadek ciśnienia z wysokością i prawo dyfuzji.
Pytanie na start. Ile cząsteczek powietrza uderza w ciągu sekundy w jeden centymetr kwadratowy skóry? Oszacuj rząd wielkości, zanim przeczytasz część 1 i zadanie feynmanowskie 1.
1 Ciśnienie z uderzeń cząsteczek
Zamknijmy N cząsteczek o masie m w sześciennym pudle o krawędzi L. Gaz jest rozrzedzony: cząsteczki zajmują ułamek objętości, a między zderzeniami poruszają się swobodnie. Cząsteczka o składowej prędkości vx uderza w prawą ścianę, odbija się z −vx i przekazuje ścianie popęd 2mvx (dział 3 „Pęd, zderzenia i środek masy”). Do tej samej ściany wraca po czasie 2L/vx. Średnia siła od tej cząsteczki to 2mvx/(2L/vx)=mvx2/L. Zderzenia z innymi cząsteczkami nie zmieniają wyniku: zachowują pęd, a w średnim bilansie tylko przekazują go między cząsteczkami.FLP I-39-2
Wynik · ciśnienie gazu
p=AF=L3Nm⟨vx2⟩⇒pV=31Nm⟨v2⟩=32N⟨Ek⟩.
Skorzystaliśmy z izotropii: ⟨vx2⟩=⟨vy2⟩=⟨vz2⟩=31⟨v2⟩. Ciśnienie to dwie trzecie gęstości energii kinetycznej ruchu postępowego. Kształt naczynia nie ma znaczenia.
Dla powietrza w warunkach pokojowych: n=N/V≈2,5⋅1025m−3, średnia odległość między cząsteczkami n−1/3≈3,4 nm (ok. 10 średnic cząsteczki), a gęstość ρ≈1,2kg/m3. Z p=31ρ⟨v2⟩ dostajemy ⟨v2⟩≈500 m/s: szybciej niż dźwięk (przykład 8.1).
1.1 Sprężanie adiabatyczne
Wsuwajmy powoli tłok z prędkością u. Cząsteczka odbija się od ruchomej ściany z prędkością większą o 2u, więc gaz się ogrzewa: praca tłoka trafia do energii cząsteczek. Jeśli ciepło nie ucieka (proces adiabatyczny), to dU=−pdV. Dla gazu jednoatomowego U=N⟨Ek⟩=23pV, skąd 23(pdV+Vdp)=−pdV, czyli 25VdV+23pdp=0:
Wykładnik adiabaty oznaczamy κ, jak w polskich podręcznikach; Feynman pisze γ, ale tę literę zajęły już tłumienie i czynnik Lorentza. Przy sprężaniu adiabatycznym rośnie też temperatura: TVκ−1=const.
Symulacja 1 Gaz dwuwymiarowy: twarde krążki w pudle z tłokiem. Ciśnienie (siła na jednostkę długości ściany) liczymy z popędów przekazywanych ścianom i porównujemy z NkT/A. Po prawej rozkład prędkości. Kliknij „Wszystkie z tą samą prędkością” i uruchom animację: zderzenia szybko odtwarzają rozkład Maxwella. Przesuń tłok w trakcie animacji i obserwuj temperaturę.
Zadania · część 1
5 zadań · ok. 30 min
Ile razy na sekundę cząsteczka azotu (vx≈300 m/s) uderzałaby w prawą ścianę pudła o krawędzi 10 cm, gdyby nie zderzała się z innymi? Jaką średnią siłę wywiera na tę ścianę?
Odpowiedź
vx/2L=1500 razy na sekundę. Siła mvx2/L=4,65⋅10−26⋅9⋅104/0,1≈4⋅10−20 N.
Gęstość powietrza to 1,2 kg/m³, ciśnienie 101 kPa. Oblicz średnią kwadratową prędkość cząsteczek.
Odpowiedź
3p/ρ≈503 m/s.
Gdyby wszystkie cząsteczki gazu podwoiły prędkość, ile razy wzrosłoby ciśnienie? Dlaczego nie dwa razy, skoro każde uderzenie przekazuje dwa razy większy popęd?
Odpowiedź
Cztery razy: każde uderzenie przekazuje dwa razy większy popęd, a uderzeń jest dwa razy więcej na sekundę.
Hel sprężono adiabatycznie do jednej dziesiątej objętości, zaczynając od 300 K i 1 atm. Oblicz końcowe ciśnienie i temperaturę. Porównaj z powietrzem (κ=1,4).
Odpowiedź
Hel: p=105/3≈46 atm, T=300⋅102/3≈1390 K. Powietrze: p=101,4≈25 atm, T=300⋅100,4≈754 K.
Wyprowadź pVκ=const dla gazu o energii U=2fpV (gdzie f to liczba stopni swobody, część 2) i pokaż, że κ=1+2/f.
Odpowiedź
2f(pdV+Vdp)=−pdV⇒(2f+1)VdV+2fpdp=0⇒pV1+2/f=const. Dla f=3: 5/3; f=5: 7/5.
2 Temperatura i ekwipartycja energii
Co mają wspólnego dwa gazy w tej samej temperaturze? Feynman rozważa mieszaninę lekkich i ciężkich cząsteczek, które się zderzają. Zderzenia przekazują energię od szybszych do wolniejszych, aż ustali się stan, w którym średnie energie kinetyczne ruchu postępowego są równe: ⟨21m1v12⟩=⟨21m2v22⟩. Ta wspólna wielkość to miara temperatury.FLP I-39-4
Definicja · temperatura gazu i równanie stanu
⟨21mv2⟩=23kT,pV=NkT=nRT.
k=1,380649⋅10−23 J/K (stała Boltzmanna), NA=6,02214076⋅1023mol−1, R=NAk≈8,314 J/(mol·K); od 2019 r. k i NA mają wartości dokładne z definicji. n to liczba moli. Na jeden stopień swobody ruchu postępowego przypada 21kT: ⟨21mvx2⟩=21kT.
Warto zapamiętać: w temperaturze pokojowej kT≈0,026 eV ≈401 eV. Energie wiązań chemicznych są rzędu kilku eV, więc w zwykłych temperaturach tylko nieliczne zderzenia są dość silne, by je zerwać (część 3). Z pV=NkT wynika prawo Avogadra: w tych samych warunkach równe objętości różnych gazów zawierają tyle samo cząsteczek. Mol gazu w 0 °C i pod ciśnieniem 1 atm zajmuje 22,4 l.FLP I-39-5
2.1 Ekwipartycja i ciepło właściwe
Prawo, że na każdy ruch postępowy przypada 21kT, uogólnia się: w równowadze termicznej na każdy wyraz kwadratowy w energii (składową prędkości, składową prędkości kątowej, kwadrat wychylenia oscylatora) przypada średnio 21kT. To zasada ekwipartycji. Stąd ciepła molowe przy stałej objętości:FLP I-40-5
Stopnie swobody i ciepło molowe CV (według fizyki klasycznej)
układ
wyrazy kwadratowe na cząsteczkę
U na mol
CV
κ=Cp/CV
gaz jednoatomowy (He, Ar)
3 (ruch postępowy)
23RT
23R≈12,5 J/(mol·K)
5/3
gaz dwuatomowy (N₂, O₂)
3 + 2 (obroty)
25RT
25R≈20,8
7/5
dwuatomowy z drganiami
3 + 2 + 2
27RT
27R
9/7
ciało stałe (atomy na sprężynkach)
3 + 3
3RT
3R≈25
—
Przy stałym ciśnieniu gaz wykonuje jeszcze pracę, rozszerzając się: Cp=CV+R. Ostatni wiersz to prawo Dulonga–Petita (1819): ciepło molowe większości pierwiastków stałych to ok. 25 J/(mol·K). Dlatego na kilogram aluminium (27 g/mol) potrzeba ok. 7,7 razy więcej ciepła niż na kilogram ołowiu (207 g/mol): ciepło idzie „na atom”, a nie „na kilogram”.
Zadania · część 2
5 zadań · ok. 30 min
Oblicz średnią energię kinetyczną ruchu postępowego cząsteczki w 300 K w dżulach i elektronowoltach. Czy zależy od rodzaju gazu?
Odpowiedź
23kT≈6,2⋅10−21 J ≈0,039 eV. Nie zależy: wodór i ksenon mają tę samą średnią energię, ale bardzo różne prędkości.
Ile cząsteczek jest w 1 m³ powietrza w 20 °C pod ciśnieniem 1 atm? Ile moli?
Odpowiedź
N=pV/kT≈2,5⋅1025, czyli ok. 41,6 mola.
Oblicz vrms=3kT/m dla wodoru, azotu i ksenonu w 300 K.
Odpowiedź
H₂: 1930 m/s, N₂: 517 m/s, Xe: 239 m/s. Prędkość maleje jak 1/m.
Oszacuj ciepło właściwe (na kilogram) miedzi, aluminium i ołowiu z prawa Dulonga–Petita. Wartości zmierzone to 385, 897 i 129 J/(kg·K).
Odpowiedź
3R/M: Cu 393, Al 924, Pb 120 J/(kg·K). Zgodność kilku procent.
Ile energii wewnętrznej (ruchu cząsteczek) zawiera 1 m³ powietrza pod ciśnieniem 1 atm? Czy zależy od temperatury przy stałym ciśnieniu?
Odpowiedź
U=25pV≈253 kJ. Przy stałym ciśnieniu nie zależy od temperatury: w cieplejszym powietrzu cząsteczek jest odpowiednio mniej. Ogrzewając pokój, nie zwiększamy energii zawartego w nim powietrza; nadmiar wypływa szparami razem z częścią powietrza.
3 Atmosfera i prawo Boltzmanna
Dlaczego powietrze nie opada całe na ziemię i nie ucieka w kosmos? Wyobraźmy sobie słup gazu w stałej temperaturze. Warstwa o grubości dh musi być podtrzymywana przez różnicę ciśnień: dp=−nmgdh. Ponieważ p=nkT i T jest stałe, kTdn=−nmgdh, więc gęstość maleje wykładniczo.FLP I-40-1
Wynik · atmosfera izotermiczna
n(h)=n0e−mgh/kT=n0e−h/H,H=mgkT.
Wysokość skali H to wysokość, na której gęstość i ciśnienie maleją e razy. Dla powietrza w 288 K: H≈8,4 km; ciśnienie spada o połowę co ok. 5,8 km.
Licznik wykładnika, mgh, to energia potencjalna cząsteczki na wysokości h. Boltzmann uogólnił ten wynik: w równowadze termicznej gęstość cząstek w miejscu, gdzie mają energię potencjalną U, jest proporcjonalna do e−U/kT, niezależnie od rodzaju siły. Jeszcze ogólniej: prawdopodobieństwo stanu o energii E jest proporcjonalne do czynnika Boltzmanna.FLP I-40-2
Czynnik Boltzmanna rządzi wszystkim, co wymaga pokonania bariery energetycznej Ea: parowaniem cieczy (cząsteczka musi mieć energię wystarczającą do wyrwania się z powierzchni), emisją elektronów z gorącej katody i szybkością reakcji chemicznych. Szybkość reakcji jest proporcjonalna do e−Ea/kT (prawo Arrheniusa). Dla typowego Ea≈50 kJ/mol (≈0,5 eV) wzrost temperatury o 10 K w pobliżu 300 K przyspiesza reakcję prawie dwukrotnie. Dlatego lodówka przedłuża trwałość żywności, a gorączka przyspiesza metabolizm.FLP I-42-1, 42-4
To samo prawo opisuje drobinki zawiesiny w wodzie: zamiast m trzeba wziąć masę pomniejszoną o wypór. Ich „atmosfera” ma wysokość skali rzędu mikrometrów, a Perrin, licząc drobinki na różnych wysokościach pod mikroskopem, wyznaczył z H stałą k niezależnie od pomiaru ruchów Browna.
Symulacja 2 Równowaga sedymentacyjna. Cząstki w naczyniu z cieczą opadają pod wpływem grawitacji, a uderzenia cząsteczek cieczy przypadkowo je podbijają. Po prawej histogram gęstości na różnych wysokościach i krzywa e−h/H. Kliknij „Wymieszaj równomiernie” i uruchom animację, żeby zobaczyć dochodzenie do równowagi.
Zadania · część 3
5 zadań · ok. 30 min
Oblicz wysokość skali dla powietrza (M=29 g/mol) w 288 K i w 250 K.
Odpowiedź
H=RT/(Mg): 8,4 km i 7,3 km.
Jakie ciśnienie (w ułamku ciśnienia na poziomie morza) przewiduje model izotermiczny na Rysach (2499 m) i na Mount Evereście (8849 m), przy H=8,4 km?
Odpowiedź
Rysy: e−0,30≈0,74. Everest: e−1,05≈0,35 (zmierzone: ok. 0,33; różnica wynika z niższej temperatury na górze).
Atomy w krysztale mogą przeskakiwać do sąsiednich położeń, pokonując barierę 1 eV. Ile razy szybciej zachodzą przeskoki w 600 K niż w 300 K?
Odpowiedź
e−E/kT2/e−E/kT1=e(1eV)(1/kT1−1/kT2)=e38,7−19,3≈2,5⋅108. Dlatego wyżarzanie i dyfuzja w ciałach stałych tak silnie zależą od temperatury.
Jaki ułamek cząsteczek wody ma w 300 K energię większą niż 0,44 eV (ciepło parowania przypadające na cząsteczkę), jeśli przyjąć, że ten ułamek to w przybliżeniu e−E/kT? Dlaczego mimo tak małego ułamka woda paruje zauważalnie?
Odpowiedź
e−0,44/0,0259≈4⋅10−8. Ale cząsteczek przy powierzchni jest bardzo dużo i zderzają się biliony razy na sekundę. W równowadze strumień parowania równa się strumieniowi skraplania, 41npary⟨v⟩≈1026m−2s−1 (przy ciśnieniu pary nasyconej 3,5 kPa), czyli ok. 107 warstw cząsteczek na sekundę. Parowanie netto w suchym powietrzu to niewielka część tego strumienia, ale wciąż zauważalna.FLP I-42-1
Wirówka obraca się z przyspieszeniem dośrodkowym 105g. Pokaż, że cząstki w niej rozkładają się jak w atmosferze, i oblicz wysokość skali dla białka o masie 10−22 kg (z uwzględnieniem wyporu; przyjmij, że wypór zmniejsza masę o 25%) w 300 K.
Odpowiedź
W układzie wirówki siła bezwładności działa jak grawitacja (dział 7 „Układy odniesienia i szczególna teoria względności”). H=kT/(mefa)=4,1⋅10−21/(0,75⋅10−22⋅9,8⋅105)≈5,6⋅10−5 m, czyli ok. 56 μm: białko zbiera się w cienkiej warstwie przy dnie probówki. Dlatego ultrawirówki rozdzielają cząsteczki według masy (a do pomiarów równowagi używa się mniejszych prędkości obrotów).
4 Rozkład prędkości Maxwella
Prawo Boltzmanna mówi też, jak rozkładają się prędkości. Energia ruchu w kierunku x to 21mvx2, więc liczba cząsteczek o składowej prędkości vx jest proporcjonalna do e−mvx2/2kT. To rozkład normalny (lekcja 5.1 Prawdopodobieństwo i rozkłady) o średniej zero i wariancji σ2=kT/m: zgodnie z ekwipartycją m⟨vx2⟩=kT. Feynman wyprowadza to z atmosfery: cząsteczka, która wzlatuje z wysokości 0 na wysokość h, musi mieć 21mvx2>mgh, a liczba takich cząsteczek musi maleć jak e−mgh/kT.FLP I-40-4
Składowe są niezależne, więc rozkład wektora prędkości to iloczyn trzech gaussów: e−m(vx2+vy2+vz2)/2kT. Zależy tylko od wartości prędkości v. Wektory o wartości z przedziału (v,v+dv) wypełniają powłokę kulistą o objętości 4πv2dv; stąd czynnik v2 w rozkładzie wartości prędkości.
Prawo · rozkład Maxwella
f(v)=4πv2(2πkTm)3/2e−mv2/2kT,∫0∞fdv=1.
Trzy charakterystyczne prędkości: najbardziej prawdopodobna vp=2kT/m, średnia ⟨v⟩=8kT/πm≈1,13vp, średnia kwadratowa vrms=3kT/m≈1,22vp.
Najważniejsza własność to ogon. Ułamek cząsteczek szybszych niż 2vp to 4,6%, szybszych niż 3vp już tylko 0,04%, a dalej spada jak e−v2/vp2. Dlatego Ziemia utrzymuje azot i tlen, a traci wodór i hel: w górnej atmosferze (ok. 1000 K) prędkość ucieczki 11,2 km/s to 3,9 vp dla wodoru, ale 14,5 vp dla azotu (przykład 8.4). Księżyc nie utrzymał gęstej atmosfery (zadanie feynmanowskie 2).
Symulacja 3 Rozkład Maxwella dla kilku gazów. Wybierz gaz i temperaturę; przerywaną linią narysowano azot w 300 K dla porównania. Przesuwaj próg prędkości: zacieniowany ogon to cząsteczki szybsze niż próg, a pole zacieniowane podaje ich ułamek. Pionowe kreski to vp, ⟨v⟩ i vrms.
Zadania · część 4
5 zadań · ok. 35 min
Oblicz vp, ⟨v⟩ i vrms dla tlenu w 300 K.
Odpowiedź
395, 446 i 484 m/s.
W jakiej temperaturze cząsteczki helu mają tę samą vrms co cząsteczki azotu w 300 K?
Odpowiedź
T/m musi być takie samo: T=300⋅4/28≈43 K.
Sprawdź przez różniczkowanie, że maksimum f(v) leży w vp=2kT/m.
Odpowiedź
dvd(v2e−av2)=(2v−2av3)e−av2=0 dla v2=1/a=2kT/m.
Pokaż, że ⟨v2⟩ z rozkładu Maxwella wynosi 3kT/m, korzystając z tego, że każda składowa ma wariancję kT/m. Dlaczego ⟨v⟩2=⟨v2⟩?
Odpowiedź
⟨v2⟩=⟨vx2⟩+⟨vy2⟩+⟨vz2⟩=3kT/m. ⟨v2⟩−⟨v⟩2 to wariancja wartości prędkości, która jest dodatnia, więc ⟨v⟩<vrms (stosunek 8/3π≈0,92).
Jaki ułamek cząsteczek ma prędkość mniejszą od vp? Użyj symulacji 3 albo wzoru P(v>vpx)=erfcx+π2xe−x2.
Odpowiedź
Dla x=1: P(v>vp)≈0,157+0,415=0,572, więc wolniejszych jest ok. 43%. Rozkład jest asymetryczny: ma długi ogon w stronę dużych prędkości.
5 Droga swobodna, dyfuzja i ruchy Browna
Cząsteczki powietrza pędzą z prędkością 500 m/s, a mimo to zapach otwartej butelki perfum rozchodzi się po spokojnym pokoju bardzo powoli. Powód: cząsteczka nie leci prosto, tylko co chwila zderza się z innymi. Cząsteczka o średnicy d zakreśla „rurę” o przekroju σ=πd2; zderzy się, gdy w tej rurze znajdzie się środek innej cząsteczki. Średnia droga między zderzeniami to droga swobodna:FLP I-43-2
ℓ=2nσ1.
Czynnik 2 uwzględnia ruch pozostałych cząsteczek. Dla powietrza (d≈0,37 nm) ℓ≈66 nm, ok. 180 średnic cząsteczki, a zderzeń jest ok. 7⋅109 na sekundę.
5.1 Błądzenie losowe i dyfuzja
Ruch od zderzenia do zderzenia to błądzenie losowe (lekcja 5.1 Prawdopodobieństwo i rozkłady). Po N krokach długości ℓ w przypadkowych kierunkach średni kwadrat przesunięcia to Nℓ2, bo przesunięcia w kolejnych krokach są nieskorelowane i ich wariancje się dodają. Ponieważ N rośnie proporcjonalnie do czasu, przesunięcie rośnie jak t.FLP I-41-4, 43-5
Prawo · dyfuzja
⟨x2⟩=2Dt(w kaz˙dym kierunku),⟨r2⟩=6Dt(w trzech wymiarach).
D to współczynnik dyfuzji (m²/s). Dla gazu D∼31ℓ⟨v⟩, czyli ok. 10−5m2/s; dla małych cząsteczek w wodzie ok. 10−9m2/s. Dyfuzja jest szybka na krótkich odległościach i bardzo wolna na długich: przez 1 mm cukier w wodzie przechodzi w ok. 20 minut, przez 1 cm w ok. dobę.
5.2 Ruchy Browna i wzór Einsteina
Drobinka zawiesiny o promieniu a w cieczy o lepkości η jest hamowana siłą Stokesa 6πηav. Zgodnie z ekwipartycją ma średnią energię ruchu postępowego 23kT, jak każda cząsteczka, tylko porusza się wolniej, bo jest ciężka. Einstein pokazał, że współczynnik dyfuzji drobinki zależy od temperatury i od tego, jak łatwo ciecz pozwala jej się poruszać:FLP I-41-1, 43-5
Wynik · wzór Einsteina–Stokesa
D=μkT=6πηakT.
μ to ruchliwość: prędkość unoszenia na jednostkę siły. Ten sam ruch cieplny, który rozpycha drobinki, odpowiada za opór, który je hamuje; dlatego D i μ są związane przez kT. Mierząc D, η i a, wyznaczamy k, a z niego NA=R/k.
Kulka o średnicy 1 μm w wodzie w 20 °C ma D≈4,3⋅10−13m2/s i przesuwa się w ciągu minuty średnio o ok. 7 μm w każdym kierunku: łatwo to zobaczyć w zwykłym mikroskopie. W zwykłym mikroskopie prędkości drobinki zmierzyć się nie da: tor jest tak poszarpany, że im częściej zaznaczamy położenie, tym dłuższą drogę i większą „prędkość” dostajemy. Perrin to właśnie odkrył i dlatego mierzył przesunięcia, a nie prędkości.
Symulacja 4 Ruchy Browna kulek w wodzie, w ten sam sposób co Perrin (on zaznaczał położenie co 30 s): położenie zaznaczone co 15 s i połączone odcinkami. Po prawej średni kwadrat przesunięcia 300 kulek w funkcji czasu i prosta 4Dt (ruch na płaszczyźnie obrazu). Zmieniaj promień kulek i temperaturę wody: lepkość wody silnie maleje z temperaturą.
Zadania · część 5
5 zadań · ok. 35 min
Porównaj drogę swobodną w powietrzu (66 nm) ze średnią odległością między cząsteczkami (3,4 nm). Ile odległości międzycząsteczkowych mieści się w drodze swobodnej i dlaczego jest ich tak dużo?
Odpowiedź
Droga swobodna (66 nm) to ok. 20 odległości międzycząsteczkowych. Cząsteczka mija po drodze wiele innych, zanim w którąś trafi, bo są małe w porównaniu z odstępami.
Kulka ma D=4,3⋅10−13m2/s. Jaka jest średnia kwadratowa przesunięcia w jednym kierunku po 1 min, 4 min i 1 h?
Odpowiedź
2Dt: 7,2 μm, 14 μm, 56 μm. Czterokrotnie dłuższy czas daje dwukrotnie większe przesunięcie.
Jak zmieni się D kulki, gdy podgrzejemy wodę z 20 °C do 80 °C (lepkość maleje z 1,00 do 0,35 mPa·s)?
Odpowiedź
D∝T/η: (353/293)⋅(1,00/0,35)≈3,4 razy większe. Zmiana lepkości daje większość efektu.
Ile czasu potrzebuje cząsteczka tlenu (D≈2⋅10−9m2/s w wodzie), by dyfundować na odległość 10 μm (komórka) i 1 mm? Dlaczego duże zwierzęta potrzebują krwiobiegu?
Odpowiedź
t≈x2/2D: 10 μm w 0,025 s, 1 mm w 250 s, 1 cm w ok. 7 h. W skali komórki dyfuzja wystarcza, w skali organu już nie; serce i naczynia przenoszą tlen konwekcją.
Pokaż, że błądzenie losowe w trzech wymiarach o kroku ℓ i czasie kroku τ daje D=ℓ2/6τ. Oszacuj D dla powietrza z ℓ=66 nm i τ=ℓ/⟨v⟩.
Odpowiedź
Po N=t/τ krokach ⟨r2⟩=Nℓ2=6Dt, więc D=ℓ2/6τ=ℓ⟨v⟩/6≈66⋅10−9⋅470/6≈5⋅10−6m2/s. Dokładniejsza teoria (rozkład długości kroków, przypadkowe kierunki) daje 31ℓ⟨v⟩≈10−5; zmierzone wartości dla gazów w powietrzu to (1–2)⋅10−5.
6 Metoda i przykłady rozwiązane
Metoda dla teorii kinetycznej
Przelicz na cząsteczki:n=p/kT, m=M/NA. Wiele zadań staje się prostych, gdy znamy n i m.
Energia na stopień swobody:21kT na każdy wyraz kwadratowy; sprawdź, czy dany stopień swobody nie jest „zamrożony” (skok energii ≫ kT).
Prędkości: skaluj jak T/m; wybierz właściwą średnią (vp, ⟨v⟩, vrms).
Bariery i pola: czynnik Boltzmanna e−E/kT; porównaj E z kT≈0,026 eV w 300 K.
Transport: droga swobodna ℓ=1/(2nσ), dyfuzja ⟨x2⟩=2Dt; czas rośnie z kwadratem odległości.
Sprawdzenie: jednostki, rząd wielkości (prędkości setki m/s, n∼1025m−3) i znane przypadki graniczne.
Przykład 8.1 Ciśnienie azotu z ruchu cząsteczek
Sprawdź wzór p=31ρ⟨v2⟩ dla azotu w 300 K pod ciśnieniem 1 atm, licząc vrms z temperatury, a gęstość z masy molowej.
Wynik: wychodzi dokładnie 1 atm, bo oba wzory to ta sama tożsamość pV=NkT zapisana na dwa sposoby.
Sprawdzenie: prędkość dźwięku w azocie, κkT/m≈353 m/s, jest tego samego rzędu co vrms, ale mniejsza (stosunek κ/3≈0,68). Fala dźwiękowa to zaburzenie przenoszone przez cząsteczki, więc nie może biec dużo szybciej od nich ✓.
Przykład 8.2 Zapłon w silniku Diesla
Silnik Diesla spręża powietrze (300 K, 1 atm) 20-krotnie. Jaką temperaturę i ciśnienie ma powietrze na końcu suwu sprężania, jeśli sprężanie jest adiabatyczne?
TVκ−1=const, κ=1,4: T=300⋅200,4≈994 K, czyli ok. 720 °C.
pVκ=const: p=201,4≈66 atm.
Wynik: ok. 1000 K i 66 atm. Wtryśnięte paliwo zapala się samo, bez świecy (temperatura samozapłonu oleju napędowego to ok. 500–600 K).
Sprawdzenie: izotermicznie ciśnienie wzrosłoby tylko 20 razy; adiabatycznie więcej, bo gaz jest też cieplejszy: 66=20⋅(994/300) ✓. W prawdziwym silniku ciepło częściowo ucieka do ścianek, więc temperatura jest nieco niższa.
Przykład 8.3 Powietrze na Evereście
Ile tlenu dostaje himalaista na szczycie Mount Everestu (8849 m) w porównaniu z poziomem morza? Porównaj model izotermiczny z pomiarem (ok. 33,7 kPa).
Model izotermiczny, H=8,4 km: p/p0=e−8849/8430≈0,35, czyli ok. 35 kPa.
Skład powietrza jest stały (mieszanie), więc ciśnienie cząstkowe tlenu też spada do ok. 0,33–0,35 wartości z poziomu morza.
Wdychane powietrze ogrzewa się w płucach do ok. 310 K niezależnie od temperatury na zewnątrz, więc liczba cząsteczek w oddechu o tej samej objętości płuc jest proporcjonalna do samego ciśnienia: ok. 0,33. Po odjęciu ciśnienia pary wodnej w płucach (6,3 kPa) ciśnienie cząstkowe tlenu we wdychanym powietrzu to (33,7−6,3)/(101,3−6,3)≈0,29 wartości z poziomu morza.
Wynik: każdy oddech dostarcza ok. 1/3 tlenu z poziomu morza, a ciśnienie cząstkowe tlenu w płucach to niecałe 0,3 wartości z poziomu morza. Dlatego większość wspinaczy używa butli.
Sprawdzenie: model izotermiczny nieco zawyża ciśnienie na dużej wysokości, bo w zimniejszym powietrzu H jest mniejsze ✓.
Przykład 8.4 Dlaczego Ziemia traci wodór, a zatrzymuje azot
W egzosferze (powyżej ok. 500 km) temperatura sięga ok. 1000 K, a zderzenia są tak rzadkie, że szybka cząsteczka lecąca w górę może uciec. Porównaj ułamek cząsteczek H₂, He i N₂ szybszych od prędkości ucieczki (11,2 km/s).
vp=2kT/m w 1000 K: H₂ 2870 m/s, He 2040 m/s, N₂ 770 m/s; stosunek x=vucieczki/vp: 3,9; 5,5; 14,5.
Ułamek P(v>xvp)=erfcx+π2xe−x2: H₂ 1,1⋅10−6, He 5⋅10−13, N₂ 10−90.
Ułamki są uzupełniane, bo zderzenia ciągle odtwarzają ogon rozkładu. Cząsteczek w egzosferze jest bardzo dużo, a czas geologiczny długi: przy 10−6 wodór ucieka z górnej atmosfery w czasie znikomym w skali geologicznej, przy 10−90 azot nie ucieknie nigdy.
Wynik: wodór i (wolniej) hel uciekają z Ziemi, azot i tlen nie. Dlatego w atmosferze praktycznie nie ma wolnego wodoru, a hel trzeba wydobywać z gazu ziemnego.
Sprawdzenie: Jowisz, z prędkością ucieczki 60 km/s i zimną atmosferą, zatrzymał wodór i składa się głównie z niego ✓.
Przykład 8.5 Ile ciepła potrzeba na kilogram?
Kilogram aluminium i kilogram ołowiu ogrzewamy o 10 K. Ile ciepła potrzeba w każdym przypadku według prawa Dulonga–Petita? Ile potrzeba na kilogram powietrza przy stałej objętości?
Wynik: 9,2 kJ, 1,2 kJ i 7,2 kJ. Liczy się liczba cząstek i ich stopni swobody, a nie masa.
Sprawdzenie: zmierzone ciepła właściwe: Al 897 J/(kg·K), Pb 129, powietrze (cV) 718; zgodność kilku procent ✓.
Przykład 8.6 Liczba Avogadra z ruchów Browna
Kulki o promieniu 0,50 μm w wodzie w 20 °C (η=1,00 mPa·s) obserwowano przez 60 s. Średni kwadrat przesunięcia wzdłuż osi x wyniósł 51,5 μm². Wyznacz k i NA.
Wynik:k≈1,38⋅10−23 J/K, NA≈6,0⋅1023: z obserwacji drobinek pod mikroskopem policzyliśmy atomy.
Sprawdzenie: dane są wymyślone tak, by dać dokładny wynik; Perrin z prawdziwych pomiarów (na wielu kulkach i z niepewnością rozmiaru) dostawał wartości ok. 6–7⋅1023. Wynik jest odwrotnie proporcjonalny do ⟨x2⟩, więc 10% błędu pomiaru przesunięć kwadratowych daje 10% błędu NA ✓.
7 Błędne intuicje, których trzeba się oduczyć
„Temperatura to ilość ciepła w ciele.”
Temperatura mierzy średnią energię na stopień swobody, a ciepło to energia przekazana z powodu różnicy temperatur. Wanna letniej wody zawiera więcej energii wewnętrznej niż rozżarzona iskra.
„W danej temperaturze wszystkie cząsteczki mają tę samą prędkość.”
Prędkości mają rozkład Maxwella: w powietrzu są cząsteczki prawie nieruchome i takie, które pędzą 2 km/s. Temperatura wyznacza tylko jego szerokość.
„Cięższe gazy opadają i zalegają przy ziemi.”
W skali atmosfery mieszanie przeważa nad segregacją. CO₂ gromadzi się tylko lokalnie, przy braku ruchu powietrza (piwnice, jaskinie, zagłębienia przy wulkanach).
„Zapach rozchodzi się przez dyfuzję.”
Na odległość metrów dyfuzja jest zbyt wolna; zapach przenosi ruch powietrza.
„Ciśnienie gazu bierze się z odpychania cząsteczek.”
W gazie rozrzedzonym ciśnienie to pęd przenoszony na ściany przez uderzające cząsteczki.
„Ciężkie cząsteczki mają w tej samej temperaturze mniejszą energię.”
Mają tę samą średnią energię kinetyczną ruchu postępowego, 23kT, tylko mniejsze prędkości.
„W temperaturze zera bezwzględnego wszystko staje w miejscu.”
Według fizyki klasycznej tak, ale mechanika kwantowa zostawia drgania zerowe; hel pod normalnym ciśnieniem nie krzepnie nawet w pobliżu 0 K.
8 Zadania feynmanowskie
Te zadania wymagają więcej myślenia niż liczenia.
Zadania feynmanowskie
8 zadań · ok. 2 godziny
Deszcz cząsteczek. Liczba cząsteczek uderzających w jednostkę powierzchni ściany w jednostce czasu wynosi 41n⟨v⟩ (czynnik 41 pochodzi z uśrednienia po kierunkach). Oblicz ją dla powietrza w 20 °C. W komorze próżniowej przy 10−8 Pa czysta powierzchnia ma pozostać czysta: po jakim czasie pokryje się jedną warstwą cząsteczek (ok. 1015 na cm²), jeśli każda uderzająca cząsteczka przylgnie?
Rozwiązanie
n=2,5⋅1025m−3, ⟨v⟩≈463 m/s: strumień ≈2,9⋅1027m−2s−1=2,9⋅1023 na cm² na sekundę. Pod ciśnieniem atmosferycznym warstwa powstaje w kilka nanosekund.
Przy 10−8 Pa: n jest 1013 razy mniejsze, strumień ok. 2,9⋅1010 na cm² na sekundę, a warstwa powstaje w ok. 10 godzin. Dlatego badania powierzchni wymagają ultrawysokiej próżni.
Dlaczego Księżyc nie ma atmosfery. Prędkość ucieczki z Księżyca to 2,38 km/s, a temperatura powierzchni w dzień ok. 390 K. Oblicz ułamek cząsteczek N₂ i H₂O szybszych od prędkości ucieczki i oceń, czy Księżyc mógł utrzymać atmosferę przez 4,5 mld lat.
Rozwiązanie
vp w 390 K: N₂ 480 m/s (x≈5,0), H₂O 600 m/s (x≈4,0). Ułamki: ≈10−10 i ≈7⋅10−7.
Księżyc nie ma gęstej atmosfery, więc każda szybka cząsteczka od razu ucieka, a ogon rozkładu uzupełniają odbicia od nagrzanego gruntu: przy każdym cząsteczka „losuje” prędkość na nowo. Jeden skok trwa ok. 2vz/g≈400 s, czyli kilkadziesiąt tysięcy losowań na rok; dla N₂ daje to czas życia rzędu 105 lat. Atmosfera zniknęłaby w czasie dużo krótszym niż wiek Księżyca. Dodatkowo cząsteczki zdmuchuje wiatr słoneczny.
Balon na ogrzane powietrze. Jaką masę unosi 1 m³ powietrza ogrzanego do 100 °C w otoczeniu o temperaturze 15 °C? Ile metrów sześciennych potrzeba, by unieść 400 kg (kosz, palnik, dwie osoby i powłoka)? Porównaj z helem.
Rozwiązanie
Przy tym samym ciśnieniu gęstość jest odwrotnie proporcjonalna do T: Δρ=ρ0(1−288/373)≈1,225⋅0,228≈0,28kg/m3.
400/0,28≈1400 m³: balon o średnicy ok. 14 m. Hel unosi ok. 1,06 kg na m³, prawie cztery razy więcej, ale jest drogi i ucieka.
Zamrożone drgania. Cząsteczka N₂ drga z częstotliwością 7,1⋅1013 Hz. Według mechaniki kwantowej energia drgań zmienia się porcjami hν. Porównaj hν z kT w 300 K i wyjaśnij, dlaczego azot ma CV=25R, a nie 27R. W jakiej temperaturze drgania zaczną się liczyć?
Rozwiązanie
hν≈0,29 eV, ok. 11 razy więcej niż kT≈0,026 eV. Ułamek cząsteczek wzbudzonych to w przybliżeniu e−hν/kT≈e−11≈10−5: drgania prawie nie przyjmują energii.
Drgania zaczynają się liczyć, gdy kT∼hν, czyli przy ok. 3400 K (w praktyce wyraźnie od ok. 1000 K). Ta porażka ekwipartycji była jedną z pierwszych wskazówek, że fizyka klasyczna jest niepełna (Feynman, I-40-6).
Dźwięk i cząsteczki. Prędkość dźwięku w gazie to κp/ρ. Pokaż, że stosunek prędkości dźwięku do vrms to κ/3, i oblicz go dla powietrza i helu. Dlaczego głos po wdechu helu brzmi wyżej?
Rozwiązanie
p/ρ=kT/m=vrms2/3, więc cs/vrms=κ/3: powietrze 0,68, hel 0,75.
W helu prędkość dźwięku to ok. 1000 m/s, prawie trzy razy więcej niż w powietrzu. Częstotliwości rezonansowe krtani i jamy ustnej rosną proporcjonalnie, więc barwa głosu przesuwa się w górę (podstawowa częstotliwość drgań strun głosowych zmienia się mniej).
Granica czułości wagi skręceń. Lusterko zawieszone na cienkim włóknie kwarcowym o sztywności skrętnej κs=10−9 N·m/rad drga przypadkowo. Oblicz średnie kwadratowe wychylenie kątowe w 300 K. Dlaczego nie da się go zmniejszyć lepszym wyważeniem?
Rozwiązanie
Ekwipartycja: 21κs⟨θ2⟩=21kT, więc θrms=kT/κs≈2⋅10−6 rad.
To ruchy Browna całego lusterka: uderzenia cząsteczek powietrza i drgania cieplne włókna. Zmniejszyć je można tylko obniżając temperaturę albo zwiększając sztywność (a więc tracąc czułość). Kappler w 1931 r. wyznaczył z takich pomiarów NA z dokładnością ok. 1%.
Lodówka i Arrhenius. Procesy psucia żywności mają energię aktywacji rzędu 50 kJ/mol. Ile razy wolniej zachodzą w lodówce (4 °C) niż w pokoju (20 °C)? A w zamrażarce (−18 °C), gdyby szybkość wciąż zależała tylko od czynnika Boltzmanna?
Rozwiązanie
exp[REa(2771−2931)]≈3,3 razy wolniej.
Zamrażarka: exp[REa(2551−2931)]≈21 razy wolniej z samego czynnika Boltzmanna; w rzeczywistości dużo bardziej, bo zamarznięta woda przestaje być środowiskiem reakcji i bakterie się nie rozmnażają.
Perfumy w pokoju. Oszacuj, ile czasu potrzebowałby zapach, by przez samą dyfuzję (D≈2⋅10−5m2/s) dotrzeć na odległość 5 m. Jakiej prędkości przepływu powietrza potrzeba, by konwekcja była szybsza?
Rozwiązanie
t≈x2/2D=25/(4⋅10−5)≈6⋅105 s, ok. tygodnia.
Konwekcja wygrywa, gdy x/u<x2/2D, czyli u>2D/x≈10−5 m/s, czyli ok. 1 cm na 15 minut. Takie przepływy powietrza są w każdym pokoju: od różnic temperatur, oddechu, otwierania drzwi.
9 Sprawdź się
Odpowiedz w myślach, zanim rozwiniesz odpowiedź.
1 Skąd bierze się ciśnienie gazu?
Z pędu przekazywanego ścianom przez uderzające cząsteczki: pV=31Nm⟨v2⟩.
2 Co mają równe dwa gazy w tej samej temperaturze?
Średnią energię kinetyczną ruchu postępowego na cząsteczkę, 23kT; nie prędkości.
3 Ile wynosi kT w temperaturze pokojowej?
Ok. 0,026 eV ≈1/40 eV, czyli 4⋅10−21 J.
4 Dlaczego CV azotu to 25R?
Trzy stopnie swobody ruchu postępowego i dwa obrotowe; drgania są „zamrożone”, bo hν≫kT.
5 Jak gęstość atmosfery izotermicznej zależy od wysokości?
Wykładniczo: n0e−mgh/kT, z wysokością skali ok. 8 km dla powietrza.
6 Jak zmienia się vrms, gdy temperatura rośnie czterokrotnie?
Rośnie dwukrotnie, bo vrms=3kT/m.
7 Jak średnie przesunięcie w błądzeniu losowym zależy od czasu?
Jak t: ⟨x2⟩=2Dt.
8 Co zmierzył Perrin i co z tego wyznaczył?
Przesunięcia (i rozkład wysokości) drobinek zawiesiny; z D=kT/6πηa wyznaczył k, a z niego liczbę Avogadra.
10 Karta wzorów i co dalej
Gaz doskonały i ekwipartycja
pV=31Nm⟨v2⟩=NkT=nRT
⟨21mv2⟩=23kT; 21kT na wyraz kwadratowy
k=1,381⋅10−23 J/K, R=8,314 J/(mol·K); kT300≈0,026 eV
CV: 23R, 25R, ciało stałe 3R; Cp=CV+R
adiabata: pVκ=const, TVκ−1=const
Rozkłady i transport
Boltzmann: P∝e−E/kT; atmosfera n0e−h/H, H=kT/mg
Maxwell: f(v)=4πv2(2πkTm)3/2e−mv2/2kT
vp:⟨v⟩:vrms=2:8/π:3 (razy kT/m)
ℓ=1/(2nσ); ⟨x2⟩=2Dt
Einstein–Stokes: D=kT/(6πηa)
Do przeczytania u Feynmana
Obowiązkowo: rozdz. 39 (kinetyczna teoria gazów), 40 (zasady mechaniki statystycznej).
Dla chętnych: rozdz. 41 (ruchy Browna), 42 (zastosowania: parowanie, kinetyka chemiczna), 43 (dyfuzja), a także rozdz. 1 (atomy w ruchu).
Następny dział: Termodynamika i entropia. Wiemy już, czym jest ciepło na poziomie cząsteczek. Teraz zapytamy, ile pracy można z niego odzyskać. Odpowiedź Carnota nie zależy od żadnych szczegółów budowy materii, a pojęcie entropii połączy oba punkty widzenia: termodynamiczny i statystyczny.