Learning #24 – Termodynamika: entropia, drugie prawo, przejścia fazowe i mechanika statystyczna

Termodynamika jest teorią nietypową wśród teorii fizycznych: nie czyni żadnych założeń o mikroskopowej budowie materii. Jej prawa to makroskopowe stwierdzenia dotyczące mierzalnych wielkości — temperatury, ciśnienia, objętości, ciepła, pracy — które obowiązują uniwersalnie, niezależnie od tego, czy substancją roboczą jest gaz doskonały, gumowa taśma, czarna dziura czy komórka biologiczna. Ta ogólność jest zarówno jej siłą, jak i jej tajemnicą.

Mechanika statystyczna, rozwinięta przez Boltzmanna, Gibbsa i Maxwella w drugiej połowie XIX wieku, uzupełnia obraz mikroskopowy: makroskopowe wielkości termodynamiczne wyłaniają się z uśredniania po ogromnej liczbie mikroskopowych stopni swobody. Związek między tymi dwoma poziomami opisu — wyrażony najzwięźlej w napisie na nagrobku Boltzmanna S = k log W — jest jednym z najgłębszych spostrzeżeń w całej fizyce.

1. Prawa termodynamiki — zwięzłe sformułowanie

Cztery prawa termodynamiki to jedne z najbardziej skondensowanych sformułowań zasad fizycznych, jakie kiedykolwiek stworzono. Każde z nich koduje głębokie ograniczenie na to, jakie procesy są możliwe w przyrodzie.

Prawo zerowe (równowaga termiczna):
  Jeżeli A jest w równowadze z C, a B jest w równowadze z C, to A i B są w
  równowadze ze sobą nawzajem. To definiuje temperaturę jako właściwość —
  wielkość, która jest równa dla układów w kontakcie termicznym.

Pierwsza zasada (zachowanie energii):
  dU = δQ − δW
  U = energia wewnętrzna (funkcja stanu)
  δQ = ciepło dostarczone do układu (zależy od drogi procesu)
  δW = praca wykonana przez układ (zależy od drogi procesu)
  Dla procesu odwracalnego: dU = T dS − P dV + µ dN + …

Druga zasada (entropia):
  W dowolnym procesie samorzutnym całkowita entropia układu izolowanego nigdy nie maleje.
  dS_total ≥ 0, z równością dla procesów odwracalnych.
  Sformułowanie Clausiusa: ciepło nie może samorzutnie płynąć od zimnego do gorącego.
  Sformułowanie Kelvina: żaden cykl nie może w całości zamienić ciepła w pracę, korzystając z jednego rezerwuaru.

Trzecia zasada (zero bezwzględne):
  Gdy T → 0, entropia idealnego kryształu dąży do stałej (zwykle przyjmowanej jako zero).
  S(T→0) = 0 dla niezdegenerowanego stanu podstawowego.
  Konsekwencja: T = 0 K jest nieosiągalne w skończonej liczbie kroków.

Funkcje stanu a funkcje procesu:
  Funkcje stanu: T, P, V, U, H = U+PV, A = U−TS (Helmholtza), G = H−TS (Gibbsa), S
  Funkcje procesu: Q (ciepło), W (praca) — zależą od sposobu przeprowadzenia procesu.
  Tylko różnice funkcji stanu są dobrze zdefiniowane dla danej zmiany stanu.

2. Entropia — od Clausiusa do Boltzmanna

Entropię wprowadził Rudolf Clausius w 1865 roku jako funkcję stanu zdefiniowaną przez stosunek ciepła odwracalnego do temperatury: dS = δQ_rev / T. Clausius potrafił obliczać zmiany entropii dla maszyn parowych, nie wiedząc nic o atomach. Dwadzieścia lat później Ludwig Boltzmann nadał entropii mikroskopowe znaczenie, które przekształciło ją z narzędzia księgowego w okno na naturę nieuporządkowania.

Entropia Clausiusa (1865):
  dS = δQ_rev / T
  S jest funkcją stanu; mierzalna jest tylko ΔS = ∫ δQ_rev/T.
  Dla wszechświata: ΔS_universe ≥ 0 (druga zasada)

Entropia Boltzmanna (1877):
  S = k_B ln W
  W = liczba mikrostanów zgodnych z makrostanem (krotność)
  k_B = 1,380649 × 10^−23 J/K (stała Boltzmanna)
  Wyryta na nagrobku Boltzmanna w Wiedniu.

Entropia Gibbsa (dla rozkładu prawdopodobieństwa nad mikrostanami):
  S = −k_B Σ_i p_i ln p_i
  Sprowadza się do wyrażenia Boltzmanna, gdy wszystkie W mikrostanów są jednakowo prawdopodobne.
  Stosowana w mechanice statystycznej i teorii informacji (entropia Shannona: zastąp k_B przez 1/ln 2).

Entropia i informacja:
  Entropia Shannona H = −Σ_i p_i log&sub2; p_i  (w bitach)
  Dla rozkładu jednorodnego na W wynikach: H = log&sub2; W
  Entropia fizyczna S = k_B ln 2 · H  (współczynnik przeliczenia między bitami a dżulami na kelwin)
  Skasowanie jednego bitu informacji rozprasza co najmniej k_B T ln 2 energii (granica Landauera).

Entropia i nieodwracalność (przykład: swobodne rozprężanie):
  Gaz doskonały, N cząsteczek, swobodnie podwaja swoją objętość:
  W_final / W_initial = 2^N  (każda cząsteczka niezależnie ma do dyspozycji dwukrotnie większą objętość)
  ΔS = k_B ln(2^N) = N k_B ln 2 = nR ln 2
  Dla 1 mola: ΔS ≈ 5,76 J/K — makroskopowy wzrost wynikający z mikroskopowej kombinatoryki.

Związek między entropią a prawdopodobieństwem rozwiązuje pozorny paradoks nieodwracalności. Prawa mechaniki są odwracalne w czasie: każdą trajektorię mikrostanu można odtworzyć wstecz bez naruszenia praw Newtona. A jednak procesy makroskopowe są nieodwracalne. Odpowiedź Boltzmanna: stany początkowe układów makroskopowych są niezwykle nieprawdopodobne w porównaniu ze stanami końcowymi — nie dlatego, że trajektoria odwrócona jest zabroniona, lecz dlatego, że jest astronomicznie mało prawdopodobna. Kostka lodu topniejąca w ciepłej wodzie mogłaby w zasadzie ponownie zestalić się; prawdopodobieństwo tego jest po prostu tak znikomo małe, że obserwowanie takiego zdarzenia zajęłoby wielokrotność wieku wszechświata.

3. Cykl Carnota — maksymalna sprawność

Sadi Carnot pokazał w 1824 roku, że żaden silnik cieplny pracujący między dwoma rezerwuarami o temperaturach T_H i T_C nie może być wydajniejszy niż silnik odwracalny pracujący między tymi samymi rezerwuarami. To nie jest stwierdzenie o tarciu czy niedoskonałych materiałach — to fundamentalne ograniczenie narzucone przez drugą zasadę.

Cykl Carnota (gaz doskonały):
  Krok 1: Rozprężanie izotermiczne przy T_H, pochłania Q_H
    W_12 = nRT_H ln(V_2/V_1) = Q_H
  Krok 2: Rozprężanie adiabatyczne T_H → T_C
    TV^(γ−1) = const,  W_23 = c_v(T_H − T_C)
  Krok 3: Sprężanie izotermiczne przy T_C, oddaje Q_C
    W_34 = −nRT_C ln(V_3/V_4) = −Q_C
  Krok 4: Sprężanie adiabatyczne T_C → T_H
    W_41 = c_v(T_C − T_H)

Sprawność Carnota:
  η_C = W_net / Q_H = 1 − T_C/T_H   (temperatury w kelwinach)
  To maksymalna sprawność dowolnego silnika cieplnego między T_H i T_C.
  Wszystkie silniki odwracalne między tymi samymi rezerwuarami mają tę samą sprawność.

Dlaczego η < 1:
  Pewna ilość ciepła Q_C = Q_H(1 − η) zawsze musi zostać oddana do zimnego rezerwuaru.
  Wymaga tego druga zasada: entropia oddana do zimnego rezerwuaru
  musi być co najmniej tak duża, jak entropia pobrana z gorącego rezerwuaru.
  ΔS_univ = −Q_H/T_H + Q_C/T_C ≥ 0  ⇒  Q_C/T_C ≥ Q_H/T_H  ⇒  η ≤ 1 − T_C/T_H

Typowe sprawności rzeczywistych elektrowni:
  Turbina węglowa / gazowa:   η ~ 40–55%  (T_H ~ 600°C, T_C ~ 30°C ⇒ η_C ~ 66%)
  Cykl kombinowany (CCGT): η ~ 60–63%
  Silnik benzynowy:      η ~ 25–35%  (duża ilość ciepła traconego w spalinach)
  Silnik diesla:        η ~ 40–50%
  Termika słoneczna:        η ~ 20–35%

4. Rozkład Maxwella-Boltzmanna — statystyka cząsteczek gazu

W gazie znajdującym się w równowadze termicznej cząsteczki poruszają się z szerokim zakresem prędkości. James Clerk Maxwell (1860) i Ludwig Boltzmann (1872) wyprowadzili rozkład prawdopodobieństwa prędkości cząsteczek od podstaw — przełomowy wynik, który otworzył drzwi do mechaniki statystycznej.

Gęstość prawdopodobieństwa prędkości v (3D):
  f(v) = 4π (m/2πk_BT)^(3/2) · v² · exp(−mv² / 2k_BT)

Prędkości charakterystyczne:
  Prędkość najbardziej prawdopodobna:  v_p = √(2k_BT/m) = √(2RT/M)
  Prędkość średnia:           <v> = √(8k_BT/πm) = √(8RT/πM)     (  ≈ 1,596 v_p / √(2) )
  Prędkość średniokwadratowa:            v_rms = √(3k_BT/m) = √(3RT/M)            (używana w ½mv² = ½k_BT ×3 )
  Kolejność: v_p < <v> < v_rms

W 300 K, powietrze (∼ N&sub2;, M = 0,028 kg/mol):
  v_p ≈ 422 m/s,  <v> ≈ 476 m/s,  v_rms ≈ 517 m/s

Rozkład energii (rozkład energii Maxwella-Boltzmanna):
  g(ε) = 2π (1/πk_BT)^(3/2) · √ε · exp(−ε/k_BT)
  Średnia energia kinetyczna na cząsteczkę: <ε> = ¾k_BT

Zasada ekwipartycji:
  Każdy kwadratowy stopień swobody wnosi ½k_BT średniej energii.
  Gaz doskonały jednoatomowy: 3 stopnie swobody translacyjnej → U = ¾Nk_BT = ¾nRT
  Dwuatomowy (sztywny):    5 stopni swobody (3 translacyjne + 2 rotacyjne) → U = 5/2nRT   (c_v = 5/2R)
  Dwuatomowy (drgający): 7 stopni swobody przy wysokiej T → U = 7/2nRT  (drgania aktywują się przy wysokiej T)

Istnienie rozkładu Maxwella-Boltzmanna oznacza, że w każdej danej temperaturze niewielka frakcja cząsteczek zawsze ma dość energii, by opuścić powierzchnię cieczy (parowanie), pokonać barierę aktywacji (reakcje chemiczne) lub przeniknąć przez barierę jądrową (fuzja termojądrowa w gwiazdach). To właśnie te zdarzenia z wysokoenergetycznego „ogona" rozkładu napędzają znaczną część chemii i astrofizyki, dlatego szybkości reakcji są tak wrażliwe na temperaturę, a synteza termojądrowa wymaga temperatur rzędu dziesiątków milionów stopni.

5. Demon Maxwella — informacja i druga zasada

W 1867 roku James Clerk Maxwell zaproponował eksperyment myślowy, który zdawał się obalać drugą zasadę termodynamiki. Maleńka inteligentna istota — nazwana później demonem Maxwella — siedzi przy klapie między dwiema komorami gazu. Obserwując poszczególne cząsteczki i selektywnie otwierając klapę, mogłaby sortować szybkie cząsteczki do jednej komory, a wolne do drugiej, tworząc różnicę temperatur z początkowo jednorodnego gazu bez wykonywania jakiejkolwiek pracy. Rozwiązanie tego paradoksu zajęło niemal stulecie.

Pierwotny paradoks (1867):
  Demon obserwuje każdą cząsteczkę i odpowiednio otwiera/zamyka klapę.
  Wynik: szybkie cząsteczki gromadzą się po prawej stronie, wolne po lewej.
  Powstaje różnica temperatur ΔT → silnik cieplny może wykonać pracę.
  Brak widocznego nakładu pracy → naruszenie drugiej zasady?

Rozwiązanie poprzez teorię informacji (Szilard 1929, Bennett 1982, Landauer 1961):
  Demon musi przechowywać w pamięci informację o prędkości każdej cząsteczki.
  Aby działać cyklicznie, demon musi ostatecznie skasować swoją pamięć.
  Zasada Landauera: skasowanie jednego bitu informacji rozprasza co najmniej k_B T ln 2:
    ΔQ_erase ≥ k_BT ln 2 na bit
  Koszt entropowy kasowania dokładnie kompensuje zmniejszenie entropii gazu.

Granica Landauera (minimalna energia na operację logiczną przy 300 K):
  k_BT ln 2 ≈ 2,87 × 10^−21 J ≈ 0,018 eV ≈ 0,0179 meV
  Współczesne procesory rozpraszają ~1000× więcej na operację → wciąż daleko od granicy Landauera.
  Ostateczna granica obliczeniowa (fizyczna): ~10^21 operacji na dżul przy 300 K.

Kluczowa myśl:
  Informacja ma fizyczną realność. Pozyskiwanie informacji o układzie nie
  wymaga pracy; kasowanie (lub nadpisywanie) informacji — tak. To właśnie
  stąd bierze się termodynamiczna strzałka czasu w obliczeniach.

Współczesne realizacje:
  Eksperymentalne silniki Szilarda (lata 2010.): jednoelektronowe skrzynki lub
  cząstki koloidalne w pułapkach ze sprzężeniem zwrotnym wydobyły pracę bliską
  kT na bit informacji, potwierdzając eksperymentalnie model Szilarda-Landauera.

6. Przejścia fazowe — porządek, symetria i zjawiska krytyczne

Przejście fazowe to gwałtowna, makroskopowa zmiana właściwości materiału wywołana płynną zmianą parametru kontrolnego, takiego jak temperatura czy ciśnienie. Lód topnieje, woda wrze, magnes traci namagnesowanie w temperaturze Curie, nadprzewodnik przechodzi w nowy stan przy T_c. Model Isinga — najprostszy model ferromagnetyka — uchwytuje istotną fizykę przejść fazowych i zjawisk krytycznych w postaci, która jest ściśle rozwiązywalna w dwóch wymiarach.

Klasyfikacja przejść fazowych (Ehrenfest):
  Pierwszego rzędu: nieciągłość pierwszej pochodnej G (objętość, entropia)
    → ciepło utajone, współistnienie faz (np. topnienie, wrzenie)
  Drugiego rzędu (ciągłe): nieciągłość drugiej pochodnej G (pojemność cieplna)
    → brak ciepła utajonego; parametr porządku rośnie w sposób ciągły od zera
    → rozbieżna długość korelacji i podatność przy T_c

Model Isinga (sieć spinów σ_i = ±1):
  Hamiltonian: H = −J Σ_{<ij>} σ_i σ_j − h Σ_i σ_i
  J > 0: sprzężenie ferromagnetyczne (preferowane spiny zgodne)
  h: pole zewnętrzne
  Rywalizacja między J (uporządkowanie) a k_BT (nieporządek termiczny) napędza przejście.

Temperatura krytyczna:
  Przybliżenie pola średniego: k_B T_c = zJ (z = liczba koordynacyjna, np. 4 dla 2D sieci kwadratowej)
  Rozwiązanie dokładne 2D (Onsager 1944): k_B T_c = 2J / ln(1 + √2) ≈ 2,269 J
  Przybliżenie pola średniego zawyża T_c; w pobliżu T_c kluczową rolę odgrywają fluktuacje.

Parametr porządku i wykładniki krytyczne (2D klasa uniwersalności Isinga):
  Namagnesowanie: m ~ (T_c − T)^β      β = 1/8
  Podatność: χ ~ |T − T_c|^−γ    γ = 7/4
  Długość korelacji: ξ ~ |T − T_c|^−ν   ν = 1
  Ciepło właściwe: C ~ |T − T_c|^−α    α = 0 (rozbieżność logarytmiczna)

Uniwersalność:
  Wykładniki krytyczne zależą wyłącznie od wymiarowości i symetrii, nie od
  szczegółów mikroskopowych. 2D klasa uniwersalności Isinga opisuje wiele
  układów fizycznych w pobliżu ich punktów krytycznych: gaz sieciowy
  (przejście ciecz-para), stopy binarne, sieci polimerowe, a nawet modele
  dynamiki społecznej.

Symulacja Monte Carlo (algorytm Metropolisa):
  1. Wybierz losowy spin σ_i
  2. Oblicz ΔE = zmianę energii przy odwróceniu.
  3. Zaakceptuj odwrócenie: zawsze, jeśli ΔE < 0; z prawdopodobieństwem exp(−ΔE/k_BT), jeśli ΔE > 0.
  4. Powtórz. Generuje zespół mikrostanów o rozkładzie Boltzmanna.

Dokładne rozwiązanie dwuwymiarowego modelu Isinga, uzyskane przez Larsa Onsagera w 1944 roku, było jednym z wielkich osiągnięć matematycznych fizyki XX wieku. Potwierdziło, że wykładniki krytyczne są nietrywialnymi liczbami wymiernymi (a nie przewidywaniami pola średniego), i ustaliło, że w pobliżu punktu krytycznego układy zachowują się w sposób uniwersalny, charakteryzowany wyłącznie przez ich symetrię i wymiarowość. Ta uniwersalność leży u podstaw grupy renormalizacji, za którą Kenneth Wilson otrzymał Nagrodę Nobla w dziedzinie fizyki w 1982 roku i która tworzy ramy pojęciowe dla współczesnej fizyki cząstek elementarnych, fizyki materii skondensowanej i statystycznej teorii pola.

Termodynamika poza fizyką

Zasięg pojęciowy termodynamiki sięga daleko poza maszyny parowe i gazy. Energia swobodna G = H − TS wyznacza kierunek reakcji chemicznych i procesów biologicznych: reakcja przebiega samorzutnie, gdy ΔG < 0, co oznacza, że spadek entalpii i/lub wzrost entropii wystarcza, by „popchnąć" układ w dół w skali energii swobodnej. Hydroliza ATP w komórkach, zwijanie białek, transport błonowy i replikacja DNA — wszystko to podlega temu jednemu kryterium.

Teoria informacji (Shannon, 1948) jest formalnie identyczna z mechaniką statystyczną (Jaynes, 1957): wnioskowanie oparte na zasadzie maksymalnej entropii jest podstawą bayesowskiej teorii prawdopodobieństwa i współczesnego uczenia maszynowego. Termodynamika czarnych dziur (entropia Bekensteina-Hawkinga proporcjonalna do powierzchni horyzontu) sugeruje, że sama czasoprzestrzeń może być zjawiskiem termodynamicznym o charakterze emergentnym. W każdej dziedzinie — od zimnych atomów po kosmologię — druga zasada, mówiąca, że entropia nigdy nie maleje, pozostaje jednym z najbardziej uniwersalnych stwierdzeń w całej nauce.