Mechanika statystyczna: jak termodynamika wyłania się z atomów

Gaz zawiera około 10²⁴ cząsteczek. Nie sposób śledzić każdej z osobna. A jednak potrafimy z całkowitą pewnością przewidzieć, jak gaz się zachowa — jakie będzie jego ciśnienie, temperatura, pojemność cieplna. Mechanika statystyczna wyjaśnia ten cud: prawa makroskopowe nieuchronnie wyłaniają się ze statystyki ogromnej liczby mikroskopowych oddziaływań.

Mikrostany i makrostany

Mikrostan określa dokładne położenie i pęd każdej cząstki układu. Jeśli mamy pojemnik z gazem zawierającym N = 10²⁴ cząsteczek, mikrostan jest punktem w 6N-wymiarowej przestrzeni fazowej — po trzy współrzędne położenia i trzy składowe pędu na każdą cząstkę. Makrostan natomiast opisuje wyłącznie te własności zbiorcze, które możemy zmierzyć: temperaturę, ciśnienie i objętość.

Dla dowolnego makrostanu zazwyczaj istnieje astronomicznie duża liczba zgodnych z nim mikrostanów. Kiedy rozpylasz perfumy w pokoju, wszystkie cząsteczki teoretycznie mogłyby samoistnie wrócić do flakonika — z zasady każdy mikrostan jest jednakowo prawdopodobny. Jednak liczba mikrostanów, w których cząsteczki są rozproszone po całym pomieszczeniu, przewyższa liczbę mikrostanów, w których cząsteczki są zamknięte w flakoniku, o czynnik rzędu e^(10²³). Druga zasada termodynamiki nie jest fundamentalnym prawem fizyki, lecz stwierdzeniem o przytłaczającym prawdopodobieństwie.

Wzór Boltzmanna na entropię

W 1877 roku Ludwig Boltzmann zapisał równanie, które dziś widnieje na jego nagrobku na wiedeńskim Cmentarzu Centralnym: S = k ln(W), gdzie W (od niemieckiego Wahrscheinlichkeit, czyli „prawdopodobieństwo") to liczba mikrostanów zgodnych z danym makrostanem, a k = 1,38 × 10⁻²³ J/K to stała Boltzmanna.

To jeden z najgłębszych wzorów w fizyce. Entropia nie jest nieporządkiem w jakimś mglistym, metaforycznym sensie — jest logarytmem liczby możliwości mikroskopowych. Logarytm pojawia się tu naturalnie, ponieważ zamienia zliczanie multiplikatywne w wielkość addytywną: jeśli połączymy dwa niezależne układy o W₁ i W₂ mikrostanach, połączony układ będzie miał W₁ × W₂ mikrostanów, więc łączna entropia S = k ln(W₁ × W₂) = k ln(W₁) + k ln(W₂) = S₁ + S₂, co czyni entropię wielkością ekstensywną.

Równanie Boltzmanna łączy świat mikroskopowy (rozmieszczenie cząstek) ze światem makroskopowym (mierzalną entropią termodynamiczną). To fundament, na którym zbudowana jest cała mechanika statystyczna.

Rozkład Maxwella-Boltzmanna

W gazie znajdującym się w równowadze termicznej cząstki nie poruszają się z jednakową prędkością. W 1860 roku James Clerk Maxwell wyprowadził rozkład prędkości cząstek z pierwszych zasad, a Boltzmann później go doprecyzował. Rozkład prędkości Maxwella-Boltzmanna ma postać: f(v) ∝ v² · exp(−mv²/2kT), gdzie m to masa cząstki, T to temperatura, a k to stała Boltzmanna.

Z tego rozkładu wyłaniają się trzy charakterystyczne prędkości. Prędkość najbardziej prawdopodobna — szczyt rozkładu — wynosi √(2kT/m). Prędkość średnia jest nieco wyższa: √(8kT/πm). Prędkość średnia kwadratowa (RMS), która wyznacza średnią energię kinetyczną, wynosi √(3kT/m). W tym ujęciu temperatura to po prostu średnia energia kinetyczna przypadająca na stopień swobody: każda cząstka niesie średnio (3/2)kT energii kinetycznej ruchu postępowego.

Kluczową rolę odgrywa eksponencjalny „ogon" rozkładu. Nawet w temperaturze pokojowej niewielka frakcja cząsteczek powietrza porusza się na tyle szybko, by uciec z ziemskiego pola grawitacyjnego — dlatego lżejsze cząsteczki (wodór, hel) w skali geologicznej stopniowo wyparowują z atmosfery.

Zobacz rozkład Maxwella-Boltzmanna w akcji w symulacji gazu doskonałego. Prześledź, jak prędkości poszczególnych cząstek składają się w charakterystyczną krzywą dzwonową i jak podniesienie temperatury przesuwa oraz poszerza ten rozkład.

Funkcja rozdziału

Centralnym obiektem matematycznym mechaniki statystycznej jest funkcja rozdziału (suma statystyczna): Z = Σᵢ exp(−Eᵢ/kT), gdzie sumowanie przebiega po wszystkich dostępnych poziomach energetycznych Eᵢ układu. Czynnik Boltzmanna exp(−Eᵢ/kT) daje względne prawdopodobieństwo znalezienia układu w stanie i — stany o energii znacznie większej niż kT są eksponencjalnie tłumione.

Ze znajomości Z, poprzez różniczkowanie, wynikają wszystkie wielkości termodynamiczne. Średnia energia wynosi ⟨E⟩ = −∂ln(Z)/∂β, gdzie β = 1/kT. Energia swobodna Helmholtza to F = −kT ln(Z). Entropia to S = −∂F/∂T. Ciśnienie to P = −∂F/∂V. Funkcja rozdziału zawiera pełną mechanikę statystyczną układu — mając Z, masz wszystko.

Dla gazu doskonałego Z rozkłada się na jednocząstkowe funkcje rozdziału, które można obliczyć analitycznie. Daje to dokładnie prawo gazu doskonałego PV = NkT jako naturalną konsekwencję mechaniki kwantowej i kombinatoryki.

Przejścia fazowe

Kiedy wodę ogrzewa się od 99°C do 101°C, przechodzi ona z cieczy w parę. Własności zbiorcze — gęstość, pojemność cieplna, lepkość — zmieniają się nieciągle, mimo że leżące u ich podstaw oddziaływania międzycząsteczkowe (wiązania wodorowe, siły van der Waalsa) zmieniają się z temperaturą w sposób gładki. To przejście fazowe pierwszego rodzaju, charakteryzujące się ciepłem utajonym: energię trzeba dostarczyć, by rozerwać wiązania międzycząsteczkowe, bez podnoszenia temperatury.

Przejścia fazowe drugiego rodzaju (ciągłe) są dziwniejsze. W temperaturze Curie (770°C dla żelaza) ferromagnetyk traci swoje trwałe namagnesowanie. Namagnesowanie — parametr porządku — maleje w sposób ciągły do zera w punkcie przejścia, ale podatność magnetyczna i pojemność cieplna rozbiegają się do nieskończoności. Długość korelacji (odległość, na jakiej spiny wpływają na siebie nawzajem) w punkcie krytycznym dąży do nieskończoności.

Model Isinga uchwytuje istotę ferromagnetyzmu za pomocą zwodniczo prostego hamiltonianu: H = −J Σ sᵢsⱼ, gdzie każdy spin sᵢ = ±1 oddziałuje ze swoimi najbliższymi sąsiadami. Mimo tej prostoty dwuwymiarowy model Isinga (rozwiązany dokładnie przez Larsa Onsagera w 1944 roku) wykazuje ostre przejście fazowe i ujawnia uniwersalną matematykę leżącą u podstaw wszystkich przejść drugiego rodzaju.

Dlaczego druga zasada ma charakter statystyczny

Wyprowadzenie przez Boltzmanna nieodwracalności z odwracalnych praw mechaniki wywołało zażarte ataki filozoficzne. Ernst Zermelo zwrócił uwagę, że twierdzenie Poincarégo o powrocie gwarantuje, iż każdy układ mechaniczny w końcu powróci dowolnie blisko swojego stanu początkowego — jak więc entropia może kiedykolwiek trwale rosnąć? Odpowiedź Boltzmanna miała charakter probabilistyczny: tak, entropia mogłaby teoretycznie samoistnie zmaleć, ale skala czasowa istotnych fluktuacji jest rzędu e^N kroków, gdzie N ~ 10²³. Dla mola gazu jest to niepojęcie dłuższe niż wiek wszechświata.

Rozwiązanie polega na tym, że druga zasada nie jest bezwzględnym zakazem, lecz przytłaczającą tendencją statystyczną. Stosunek liczby mikrostanów „nieuporządkowanych" do „uporządkowanych" jest tak ogromny, że dla każdego układu makroskopowego obserwowanego w dowolnej skali czasu ludzkiego wzrost entropii jest praktycznie pewny. Mechanika statystyczna zamienia termodynamikę ze zbioru reguł empirycznych w konsekwencję wyprowadzoną ze zwykłego zliczania.

Najczęściej zadawane pytania

Czym jest mechanika statystyczna?

Mechanika statystyczna to dział fizyki, który wyjaśnia makroskopowe własności termodynamiczne (temperaturę, ciśnienie, entropię) na podstawie mikroskopowego zachowania atomów i cząsteczek. Zamiast śledzić każdą cząstkę z osobna, wykorzystuje teorię prawdopodobieństwa i statystykę, aby wyprowadzić uśrednione własności makroskopowe z rozkładu stanów mikroskopowych.

Czym jest entropia w mechanice statystycznej?

W mechanice statystycznej entropię S definiuje równanie Boltzmanna S = k_B ln(Ω), gdzie k_B to stała Boltzmanna, a Ω to liczba mikrostanów zgodnych z danym makrostanem. Układ, który ma więcej sposobów rozmieszczenia swoich cząstek, ma wyższą entropię. Dzięki temu entropia zyskuje konkretne znaczenie fizyczne: mierzy mikroskopowy nieporządek albo liczbę równoważnych konfiguracji.

Czym jest rozkład Maxwella-Boltzmanna?

Rozkład Maxwella-Boltzmanna opisuje rozkład prawdopodobieństwa prędkości cząstek gazu doskonałego w stanie równowagi termicznej. Pokazuje, że większość cząstek porusza się z prędkościami umiarkowanymi, bliskimi prędkości najbardziej prawdopodobnej (√(2k_BT/m)), a cząstek o bardzo niskich lub bardzo wysokich prędkościach jest znacznie mniej. Temperatura określa szerokość rozkładu — jej wzrost przesuwa szczyt w stronę wyższych prędkości.

Jaka jest różnica między mikrostanem a makrostanem?

Mikrostan określa dokładne położenie i pęd każdej cząstki układu. Makrostan opisuje makroskopowe wielkości obserwowalne — temperaturę, ciśnienie, objętość — zgodne z wieloma różnymi mikrostanami. Mechanika statystyczna łączy je za pomocą fundamentalnego postulatu: wszystkie dostępne mikrostany o tej samej energii są jednakowo prawdopodobne (hipoteza ergodyczna).

Czym są statystyki Bosego-Einsteina i Fermiego-Diraca?

Statystyka Maxwella-Boltzmanna dotyczy rozróżnialnych cząstek klasycznych. Statystyka Bosego-Einsteina rządzi identycznymi bozonami (cząstkami o spinie całkowitym, jak fotony czy atomy helu-4), które mogą zajmować ten sam stan kwantowy — prowadząc do kondensacji Bosego-Einsteina. Statystyka Fermiego-Diraca rządzi identycznymi fermionami (o spinie połówkowym, jak elektrony), które podlegają zakazowi Pauliego — żadne dwa nie mogą dzielić tego samego stanu kwantowego — co wyjaśnia przewodnictwo metaliczne i stabilność białych karłów.

Czym jest funkcja rozdziału?

Funkcja rozdziału Z to centralna wielkość w mechanice statystycznej, definiowana jako suma (lub całka) po wszystkich mikrostanach z czynnika Boltzmanna e^(-E/k_BT). Zawiera pełną informację termodynamiczną: energię swobodną F = -k_BT ln(Z), średnią energię ⟨E⟩ = -∂ln(Z)/∂β, entropię S = -∂F/∂T. Obliczenie funkcji rozdziału to kluczowe wyzwanie w większości zagadnień mechaniki statystycznej.

Czym jest przejście fazowe w mechanice statystycznej?

Przejście fazowe to nagła jakościowa zmiana zachowania układu, gdy parametr kontrolny (temperatura, ciśnienie) przekracza wartość krytyczną. Przejścia pierwszego rodzaju (topnienie, wrzenie) wiążą się z ciepłem utajonym i nieciągłymi zmianami parametrów porządku. Przejścia drugiego rodzaju (ciągłe) wiążą się z rozbieganiem się długości korelacji i podatności, a zachowanie w pobliżu punktu krytycznego charakteryzują uniwersalne wykładniki krytyczne, niezależne od szczegółów mikroskopowych.

Czym jest model Isinga?

Model Isinga to uproszczony model ferromagnetyzmu ze spinami na sieci, które mogą być skierowane w górę (+1) albo w dół (-1). Sąsiednie spiny „wolą" ustawiać się zgodnie (sprzężenie ferromagnetyczne). W niskiej temperaturze spiny ustawiają się w stan namagnesowany; w wysokiej temperaturze fluktuacje termiczne wprowadzają nieporządek. Dwuwymiarowy model Isinga został dokładnie rozwiązany przez Larsa Onsagera w 1944 roku i pozostaje wzorcowym przykładem przejścia fazowego.

Czym jest ergodyczność?

Ergodyczność to własność polegająca na tym, że średnia czasowa wielkości fizycznej jest równa średniej po zespole statystycznym (średniej po wszystkich możliwych konfiguracjach układu). Jeśli układ jest ergodyczny, z czasem eksploruje wszystkie dostępne mikrostany. Ergodyczność jest zakładana w klasycznej mechanice statystycznej, aby powiązać zachowanie dynamiczne z opisami statystycznymi równowagi. Wiele rzeczywistych układów — szkła, szkła spinowe, układy biologiczne — jest nieergodycznych w praktycznych skalach czasowych.

Czym jest twierdzenie o ekwipartycji?

Twierdzenie o ekwipartycji mówi, że w równowadze termicznej każdy kwadratowy stopień swobody wnosi (1/2)k_BT do średniej energii, gdzie T to temperatura, a k_B to stała Boltzmanna. Jednoatomowy gaz doskonały ma 3 translacyjne stopnie swobody: średnia energia kinetyczna = (3/2)k_BT na cząstkę. Twierdzenie to wyjaśnia ciepło właściwe gazów, ciał stałych (prawo Dulonga-Petita) oraz katastrofę w nadfiolecie, której fizyka klasyczna nie potrafiła rozwiązać (skorygowaną przez mechanikę kwantową).