Termodynamika statystyczna idealnej substancji czystej i jej stanów skupienia

<a id="gesamt"></a>

Podsumowanie ogólne

Model 3 jest jakościowym, statystyczno-mechanicznym modelem substancji idealnej, który wyprowadza trzy stany skupienia — stały, ciekły i gazowy — oraz przejścia między nimi z jednego obrazu wiązań, opartego na wiązaniach molekularnych. Zasada podstawowa polega na uwzględnieniu trzech zbiorników energii dla energii wewnętrznej: energii drgań w wiązaniach molekularnych, energii translacji cząsteczek oraz kooperatywnej energii ruchu w rotacji miceli jako wyidealizowanej postaci wewnętrznych przepływów materii wyrównujących wewnętrzne wahania gęstości. Ważną wielkością podstawową jest w tym kontekście średnia liczba koordynacyjna z̄ = z_max·x, która ustala się samouzgodnienie poprzez kooperatywne wiązanie dwustanowe (typu Isinga, kooperatywność K = J·z_max, zwyrodnienie stanu zerwanego G = z_max). Wynikają z niej trzy przejścia na wspólnej, wyskalowanej jednym parametrem energii ε osi τ: topnienie jako próg fononowy/sztywnościowy (z̄ osiąga sztywność), wrzenie jako załamanie kooperatywności — właściwie łączności (perkolacji) — w punkcie krytycznym oraz opcjonalne przejście stałe(1)→stałe(2) (FCC↔BCC) z bezparametrowym stopniem c_v równym 1,0 k_B. Przejście staje się pierwszego rodzaju, gdy tylko ln G > 2 + 4ε/K. Przy tym zwyrodnienie stanu zerwanego G to liczba możliwych do zerwania kontaktów sąsiedzkich (G = z_max, maksymalna liczba koordynacyjna), energia wiązania ε to energia zyskiwana na jeden kontakt sąsiedzki, a kooperatywność K = J·z_max to iloczyn siły sprzężenia J typu Isinga i z_max; przejście staje się więc tym prędzej pierwszego rodzaju, im większa jest wielość kontaktów G i im silniejsza kooperatywność K wobec energii wiązania ε.

Artykuł prowadzi od siedmiu aksjomatów (rozdział 1) przez analizę sprzeczności (rozdział 2) i obrazowy opis ze statycznymi komiksami (rozdział 3) do rusztowania statystyczno-mechanicznego — wiązanie kooperatywne, sumy stanów i pojemności cieplne (rozdziały 4–5) —, do wiarygodności matematycznej w 3D (rozdział 6), a w końcu do kalibracji i równania stanu (rozdział 7). Myśli uzupełniono pełnym skanem 3D c_v topnienia pierwszego rodzaju (Δz̄, ciepło utajone ΔU, przesunięcie piku), modelem sprężyn Einsteina dla sztywności fononowej θ(sieć) oraz linią przemiany stałe(1)→stałe(2) poprzez Clausiusa–Clapeyrona. θ(sieć) zostaje następnie wysubtelnione do gęstości stanów fononowych typu Debye'a; przemianę stałe(1)→stałe(2) wprowadza się dodatkowo jako linię w obrazie (τ, p) (wykalibrowaną na żelazie) i dodaje osobną gałąź Redlicha–Kwonga obok domyślnej gałęzi van der Waalsa. θ(sieć) sprawdza się ponadto na pełnej sieci Lennarda-Jonesa (FCC stabilne, idealne BCC niestabilne na ścinanie), a topologię stałe(1)→stałe(2) na tytanie (drugie potwierdzenie) i argonie (kontrola negatywna). Rozdział 8 umieszcza model wśród jedenastu uznanych wyobrażeń modelowych termodynamiki statycznej. Załącznik dokumentuje źródła literaturowe (uporządkowane wg autorów) oraz wsparcie SI (szacowany udział innowacji, szacowany koszt i czas trwania na krok rozwoju); pomoce zbierają sposób czytania, pozostałe punkty krytyczne i perspektywy.

Teza główna. Jedna jedyna kooperatywna miara wiązania z̄ spina trzy skądinąd rozdzielone linie tradycji — obraz gazu sieciowego / Isinga (kooperatywne przejście dwustanowe), obraz perkolacji (łączność → wrzenie) i obraz fononowy Einsteina/Debye'a (sztywność → topnienie) — w spójny, jakościowo zupełny obraz faz, świadomie wewnątrz idealnego obrazu typu van der Waalsa, przy czym gmach teoretyczny daje się rzecz jasna wyprowadzić z łatwo wyobrażalnego, dynamicznego obrazu. Nie dokładność, lecz obrazowość jest właściwą siłą modelu i znamionuje jego naukową wartość modelową.

Spis treści

• Podsumowanie ogólne
• 1 — Aksjomaty wraz z ich hasłami (skonsolidowane)
• 2 — Analiza sprzeczności i luk
• 3 — Opisy obrazowe
• 3.1 Trzy pojedyncze fazy
• 3.2 Przejścia fazowe
• 3.3 Dlaczego prędkość dźwięku w cieczach maleje z T
• 3.4 Jak powstaje ruch Browna — z zaostrzonym pojęciem miceli
• 3.5 Dlaczego reakcje na warstwach granicznych miceli przyspieszają
• 3.6 Trójwymiarowe okno cieczy
• 4 — Model statystyczno-mechaniczny
• 4.1 Stopnie swobody i zmienne
• 4.2 Kooperatywne wiązanie dwustanowe (typu Isinga, parametr J)
• 4.3 Bilans energetyczny — z członem kooperatywnym
• 5 — Sumy stanów i pojemności cieplne
• 5.1 Rozłożona suma stanów
• 5.2 Pojemność cieplna dla każdej fazy (analitycznie)
• 5.3 Numeryczna krzywa c_v(τ) **z** członem kooperatywnym — ciepło utajone skwantyfikowane
• 5.4 Próg perkolacji łączności jako wynik (potwierdzony bez zmian) - 5.4.a Dwa różne progi czysto rozdzielone
• 5.5 Rodzaj przejścia topnienia
• 5.6 Studium skończonego rozmiaru progu perkolacji
• 5.7 Pełny skan 3D c_v(τ) topnienia pierwszego rodzaju
• 6 — Matematyczna wiarygodność interpretacji 3D
• 6.1 Zmiana stopnia swobody rotacji 2D → 3D
• 6.2 Wpływ różnych sieci krystalicznych — łączność jako cecha odróżniająca
• 6.3 Dodatkowe przejście stałe(1) → stałe(2)
• 6.4 Rozważania: anomalia akustyczna ciekłej wody
• 6.5 Alternatywa: fraktalne pola wirowe zamiast rotacji miceli
• 6.6 Mikroskopowa sztywność fononowa θ(sieć): model sprężyn Einsteina
• 6.7 Matematyczne oszacowanie metastabilności (2D triangulowana ↔ kwadratowa, 16b)
• 6.8 Wysubtelnienie θ(sieć) w duchu Debye'a
• 6.9 θ(sieć) z pełnej sieci Lennarda-Jonesa
• 7 — Kalibracja i równanie stanu
• 7.1 Metoda kalibracji (jeden parametr energii ε) — dostosowana
• 7.2 Tabela kalibracyjna
• 7.3 Oś gęstości/objętości: równanie stanu van der Waalsa z twardym rdzeniem
• 7.4 Ilościowa kalibracja c/g: a → ε, b → D²
• 7.5 Wspólna powierzchnia (τ, v) z punktem potrójnym
• 7.6 Ciepło utajone wobec entalpii topnienia
• 7.7 Ujednolicenie skal energii: ε ↔ ε_lg
• 7.8 Bezparametrowy punkt potrójny z z̄_c
• 7.9 Rachunek kontrolny 3D: czyste okno cieczy na jednej skali ε
• 7.10 Prawdziwy skan perkolacji wiązań w 3D (FCC/BCC) i wrażliwość brzegu okna
• 7.11 Prawdziwy skan τ kinetyki wiązań w 3D: τ_b(3D) i linia topnienia na jednej skali ε
• 7.12 Linia przemiany stałe(1)→stałe(2) na powierzchni (τ, v)
• 7.13 Linia przemiany stałe(1)→stałe(2) na wykresie (τ, p), wykalibrowana na żelazie
• 7.14 Redlich–Kwong jako osobna gałąź obrazu idealnego
• 7.15 stałe(1)→stałe(2) na dalszych polimorfach: tytan i kontrola gazem szlachetnym
• 8 — Porównanie z literaturą: wyobrażenia modelowe termodynamiki statycznej
• 8.1 Tabela przeglądowa
• 8.2 Modele z osobna
• 8.3 Podobieństwa i różnice wobec modelu 3
• 9 — Załącznik
• 9.1 Literatura
• 9.2 Wsparcie SI
• 9.3 Słowniczek skrótów i symboli
• 9.4 Skrypty obliczeniowe i pokrycie testami

<a id="kap1"></a>

1 — Aksjomaty wraz z ich hasłami (skonsolidowane)

Podsumowanie. Siedem opatrzonych hasłami aksjomatów ustala substancję idealną: zachowana liczba cząstek, ruch czysto klasyczny, kooperatywne wiązanie dwustanowe i średnia liczba koordynacyjna z̄ jako parametr porządku. Niosą one wszystkie dalsze wyprowadzenia; w niemieckim tekście wzorcowym znak wodny z homoglifów oznacza hasła (klucz na końcu tekstu).

Aksjomaty uściślone:

1. SWOBODNY LOT (aksjomat dodatkowy). Między dwiema cząsteczkami przelatującymi obok siebie niezwiązanymi nie działa ani tarcie, ani siły elektrostatyczne. Pęd i energia zmieniają się wyłącznie przez twarde odpychanie TWARDEGO RDZENIA przy zderzeniu oraz przez zawiązywanie i zrywanie WIĄZANIA. Aksjomat ma znaczenie wyłącznie z punktu widzenia mechaniki klasycznej i bilansu energetycznego: gwarantuje, że nie istnieją kanały dysypacyjne na poziomie molekularnym — energia ruchu nigdy nie „ginie", lecz przechodzi tylko do trzech zbiorników TRÓJDŹWIĘKU albo w energię wiązania. (Makroskopowe tarcie/lepkość powstają emergentnie z przesuwanki i PĘKNIĘCIA, nie jako mikroskopowa siła podstawowa.)
2. WIELOŚĆ. Substancja idealna składa się z bardzo wielu jednakowych cząsteczek; własności makroskopowe opisuje się statystycznie w granicy termodynamicznej (N → ∞).
3. WIĄZANIE. Każde dwie geometrycznie sąsiadujące cząsteczki mogą nieść wiązanie. Wiązanie nienaruszone obniża energię o skończoną wartość ε; jego zerwanie kosztuje dokładnie ε. Wiązanie jest układem dwustanowym {nienaruszone, zerwane}.
4. TRÓJDŹWIĘK. Energia wewnętrzna rozkłada się klasycznie na trzy zbiorniki: (a) translację cząsteczek, (b) drganie w nienaruszonych wiązaniach (sprężyna harmoniczna na wiązanie), (c) rotację sztywnych bloków cząsteczkowych (kolektywny stopień swobody obrotu spójnych obszarów).
5. TWARDY RDZEŃ. Gdy odległość cząsteczkowa spadnie poniżej wartości minimalnej σ, działa twarde odpychanie (objętość własna). Wraz z WIĄZANIEM daje to efektywny potencjał pary (twarda powłoka plus skończona studnia wiązania).
6. ZACHOWANIE. Pęd całkowity i energia całkowita są zachowane. Pojedyncza cząsteczka jest energetycznie neutralna — nie buforuje energii w sobie; energia żyje w wiązaniach (WIĄZANIE) i w kinetycznych stopniach swobody (TRÓJDŹWIĘK).
7. PĘKNIĘCIE. Wiązania rozpadają się i powstają wedle zasady pękania/gojenia: im więcej brakujących wiązań sąsiedzkich → tym większe prawdopodobieństwo zerwania; im więcej związanych sąsiadów → tym większe prawdopodobieństwo gojenia. Kinetyka jest kooperatywna.
8. RĄCZY. Zerwanie i gojenie przebiegają szybko i mało aktywowane; podukład wiązań znajduje się w każdej temperaturze w równowadze lokalnej i można go adiabatycznie oddzielić od ruchu cząsteczkowego (podstawa rozłożonej sumy stanów).

Postawa modelowania. Koordynacja geometryczna pozostaje wyidealizowanie zachowana; zmienia się tylko liczba koordynacyjna z przez zerwania. Ruch powstaje jak w przesuwance: pustki (z objętości własnej TWARDEGO RDZENIA i zerwań PĘKNIĘCIA) wędrują, bloki się przesuwają. Świadomie „fałszywy" obraz obracających się, wędrujących krystalitów zachowuje się jako poglądowa podstawa.

<a id="kap2"></a>

2 — Analiza sprzeczności i luk

Podsumowanie. Analiza systematycznie bada model pod kątem sprzeczności wewnętrznych i luk oraz kontynuuje stan prac; rozpoznane słabości są nazwane i w dalszych rozdziałach rozwiązane albo oznaczone jako uczciwa granica modelu.

Lista punktów krytycznych

1. Podwójne liczenie translacji (aksjomat 3a wobec 3c). Tekst surowy wymienia translację dwukrotnie (translacja cząsteczek i translacja bloków). Dwa niezależne zbiorniki translacji nie istniałyby. Rozwiązanie: translacja bloku jest tylko kolektywną modą już policzonych translacji cząsteczkowych; istotnie nowy jest wyłącznie obrotowy stopień swobody bloku. W bilansie stopni swobody (rozdz. 4) translację liczy się dokładnie raz.
2. „Energetycznie neutralna" (aksjomat 5) wobec „wiązania przyjmują energię" (aksjomat 2/3b). Sprzeczność pozorna. Rozwiązanie: neutralność dotyczy cząsteczki jako ciała; energia siedzi między cząsteczkami (w WIĄZANIU) i w kolektywnych modach kinetycznych, nigdy w buforze wewnątrzcząsteczkowym. Cząsteczka ≠ wiązanie — sprzeczność zostaje usunięta.
3. Brak przyciągania (aksjomat 4 wymienia tylko odpychanie). Bez przyciągania nie byłoby fazy skondensowanej. Rozwiązanie: przyciąganie dostarcza WIĄZANIE (studnia głębokości ε), odpychanie TWARDY RDZEŃ (twarda powłoka przy σ). Razem — efektywny potencjał pary.
4. Sprzężenie zwrotne bez członu przeciwnego (aksjomat 6). „Więcej brakujących sąsiadów → więcej zerwań" jest samo w sobie dodatnim sprzężeniem zwrotnym (postęp pęknięcia aż do całkowitego rozpadu). Rozwiązanie: człon gojenia (więcej sąsiadów → więcej gojenia) plus czynnik Boltzmanna skończonej energii wiązania ε (WIĄZANIE) równoważą pęknięcie i dają określone przejście zamiast rozbiegania.
5. Niewypowiedziane rozdzielenie skal czasowych (aksjomat 7). „Szybka i mało aktywowana" staje się użyteczna dopiero wtedy, gdy wyprowadzi się z tego adiabatyczne rozdzielenie dynamiki wiązań i translacji. Rozwiązanie: wyraźnie ustalone jako warunek rozkładu Z = Z_trans·Z_vib·Z_rot·Z_bond.
6. Brak geometrii i brak mechanizmu ruchu. Aksjomaty nie wymieniają ani sieci koordynacyjnej, ani tego, co umożliwia płynięcie. Rozwiązanie (jako wyraźnie zaznaczone założenie dodatkowe): stała referencyjna geometria koordynacyjna ze zmiennym obsadzeniem wiązań; jako kanoniczny przypadek 2D koordynacja trójkątna (z_max = 6). Ruchliwość przez mechanizm pustek przesuwanki. Te dwa punkty są jedynymi rzeczywistymi dodatkami poza siedmioma aksjomatami i w rozdziale 6 oznaczono je jako pytanie strategiczne.
7. Utajona luka dysypacji. Wersja 07 pozostawiała otwartym, czy niezwiązane cząsteczki tracą energię przy przelocie obok siebie. Rozwiązana przez SWOBODNY LOT: nie ma mikroskopowego tarcia ani elektrostatyki; wszystkie przepływy energii i pędu biegną wyłącznie przez zderzenia TWARDEGO RDZENIA i kinetykę WIĄZANIA. Tym samym bilans energetyczny (rozdz. 4.3) jest zamknięty, a statystyka Maxwella–Boltzmanna w granicy gazowej dokładna (swobodny gaz doskonały między zderzeniami sprężystymi).
8. Niespójność energetyczna parametru porządku. Wersja 08 używała członu kooperatywnego −(K/2)x² wprawdzie w energii swobodnej g(x) (do wygenerowania x(τ)), ale nie w energii wewnętrznej U/N. To termodynamicznie niespójne: ten sam człon pola molekularnego musi występować w obu. Rozwiązanie (rozdz. 4.3): −(K/2)x² zostaje przeniesione do U/N; ciepło utajone rośnie przez to mierzalnie (rozdz. 5.3), bez zmiany plateau ciała stałego i gazu.
9. Kalibracja c/g tylko jakościowa. Wersja 09 miała oś gęstości (równanie van der Waalsa) tylko w postaci zredukowanej. Rozwiązanie (rozdz. 7.4): bezwzględną linię wrzenia oblicza się przez punkt krytyczny (T_c, p_c), a T_b(p_atm) zestawia z wartościami zmierzonymi; systematyczne odchylenie (reguła Guldberga) zostaje nazwane.
10. Rozdzielone osie bez punktu potrójnego. Osie topnienia i wrzenia były osobnymi rachunkami. Rozwiązanie (rozdz. 7.5): obie zostają zestawione w jednej płaszczyźnie (τ, v), punkt potrójny skonstruowany, a punkt krytyczny zaznaczony.

<a id="kap3"></a>

3 — Opisy obrazowe

Podsumowanie. Rozdział 3 wyjaśnia poglądowo stan stały, ciekły i gazowy oraz ich przejścia (za pomocą statycznych komiksów): trzy pojedyncze fazy (3.1), przejścia fazowe rekrystalizacja/topnienie/wrzenie/sublimacja wraz z punktem potrójnym i krytycznym (3.2), anomalię akustyczną (3.3), ruch Browna i zaostrzone pojęcie miceli (3.4–3.5) oraz trójwymiarowe okno cieczy (3.6). Na podstawie tych opisów sporządzono schematy (komiksy), które mają pomóc unaocznić przed oczyma wyobraźni dynamiczne stosunki molekularne w poszczególnych fazach. Ta wyobrażalność jest właściwym rdzeniem modelu. Jakościowo i naukowo cenny staje się on wtedy, gdy poprzez mechanikę statystyczną da się wiarygodnie pokazać, że model czyni jakościowe tendencje i obserwacje empiryczne przewidywalnymi także bez rozbudowanej matematyki.

<a id="kap3-1"></a>

3.1 Trzy pojedyncze fazy

Ciało stałe. Niemal wszyscy sąsiedzi są ze sobą związani wiązaniami molekularnymi (WIĄZANIE); średnia liczb koordynacyjnych z cząsteczek leży blisko maksimum geometrycznego. Sieć cząsteczkowa jest jednym sztywnym blokiem bez pustek (TWARDY RDZEŃ, ciasno upakowany) — wewnętrzna przesuwanka jak w cieczy nie jest w wersji zwartej możliwa. Jedyny ruch: drganie wokół położeń spoczynkowych (drganie wiązania, TRÓJDŹWIĘK-b), określone przez wiązania molekularne z sąsiadami. Pęd i energia wahają się sprężyście w sieci (ZACHOWANIE, SWOBODNY LOT: bezstratnie). W przypadku granicznym (bardzo niska temperatura, idealne upakowanie) mamy stan idealny kryształu doskonałego.

Ciecz. Kooperatywne PĘKNIĘCIE przeniknęło sieć: wiązania brakują płatami, średnia liczb koordynacyjnych z cząsteczek spadła, sztywny blok rozpadł się wewnętrznie na wiele mniejszych. Z monolitu powstało rumowisko. Bloki te przemieszczają się i obracają (TRÓJDŹWIĘK-c). Między nimi otworzyły się pustki, przez które przesuwają się obok siebie — przesuwanka ruchu Browna działa. Lokalna geometria koordynacyjna pozostaje wyidealizowanie zachowana; zmieniło się tylko, ile wiązań rzeczywiście istnieje i się tworzy.

Gaz. PĘKNIĘCIE zwyciężyło: wiązania zerwane niemal wszędzie (z → 0). W przypadku skrajnym żadna cząsteczka nie ma już wiązania molekularnego z inną. Bloki rozpadły się na pojedyncze cząsteczki. Lecą swobodnie, sprzężone tylko zderzeniami TWARDEGO RDZENIA (pomiędzy nimi SWOBODNY LOT); prędkości podlegają rozkładowi Maxwella–Boltzmanna (ZACHOWANIE). Drganie wiązania i rotacja bloków wygasły — czysta translacja (TRÓJDŹWIĘK-a). W przypadku granicznym (wysoka temperatura, niskie ciśnienie) mamy stan idealny gazu doskonałego.

<a id="kap3-2"></a>

3.2 Przejścia fazowe

Rekrystalizacja (stałe(1)/stałe(2)). Jeszcze poniżej topnienia ciało stałe może zmienić swoją sieć, nie tracąc sztywności: dwa idealne upakowania o różnej liczbie koordynacyjnej (FCC z_max = 12 ↔ BCC z_max = 8) mają różne entalpie swobodne μ(τ, p); przy μ_FCC = μ_BCC geometria koordynacyjna przeskakuje skokowo (nie tylko obsadzenie z̄ jak przy topnieniu). Widoczny staje się mały stopień c_v (Δc_v = zmiana z_max/4 = 1,0 k_B) plus słabe utajone ciepło przemiany. Twierdzenie SI: jest to wyobrażalne tylko w 3D (dyskretna, skwantowana łączność; rozdz. 6.3). W przypadku dwuwymiarowym prawdziwą rekrystalizację, czyli przejście między dwiema fazami stałymi, łatwiej sobie wyobrazić. Na płaszczyźnie istnieją tylko trzy regularne parkietaże bez dziur. Powierzchnię można pokryć całkowicie (sparkietować) trójkątami równobocznymi, kwadratami albo sześciokątami foremnymi. Jeśli w tych sieciach krawędzie przedstawiają wiązania molekularne, a wierzchołki cząsteczki, to wyobrażalne jest przejście fazowe stałe(trójkątne) → stałe(kwadratowe). Cząsteczka w parkietażu trójkątnym miałaby wtedy jeszcze sześć wiązań z sąsiednimi cząsteczkami, podczas gdy cząsteczka w sieci kwadratowej miałaby już tylko cztery.

Topnienie (s/c). Przy rosnącej temperaturze sieć pęka w dostatecznie wielu miejscach, by sztywny blok całkowity utracił swoją sztywność (kooperatywne PĘKNIĘCIE schodzi poniżej progu sztywności Maxwella z̄ = 2d — sieć staje się wiotka). Tworzą się micele termodynamiczne, czyli obszary, które nie są już połączone wiązaniami z większymi obszarami sąsiednimi. Z molekularnego punktu widzenia monolit staje się rumowiskiem i jest warunkowo zdolny do płynięcia. Pierwsze pustki pozwalają na przesuwanie bloków cząsteczkowych, podobnie jak w przesuwance. Brakujące wiązania z otoczeniem pozwalają na wewnętrzną rotację w micelach. Rotacja miceli (TRÓJDŹWIĘK-c) rusza. Zrywanie wielu wiązań (WIĄZANIE) pochłania ciepło utajone.

Wrzenie (c/g). Dalsze ogrzewanie sprawia, że perkolują same pustki: łączą się w ciągłą przestrzeń pustą, w której bloki się rozpuszczają, a pojedyncze cząsteczki odlatują. Wrzenie to przejście od „pustek w materii" do „materii w pustce".

Sublimacja (s/g). Przy niskim ciśnieniu brakuje ciekłego obszaru pośredniego: sieć oddaje pojedyncze cząsteczki wprost do przestrzeni pustej, bez tworzenia bloków jak w micelach. Można tu rozważać, jaki wpływ na powstawanie lub niepowstawanie miceli termodynamicznych mogą mieć aspekty geometryczne (długie cząsteczki, cząsteczki wzajemnie splecione), aspekty elektrostatyczne (ustawienie przez momenty magnetyczne lub dipolowe) oraz różne typy wiązań molekularnych (van der Waalsa, oddziaływania między układami mezomerycznymi, wiązania wodorowe).

Punkt potrójny. Właśnie tam zbiegają się dwa progi: sztywność Maxwella (stałe ↔ ciekłe, utrata sztywności) i perkolacja pustek (ciekłe ↔ gaz) występują przy tym samym (τ, p). Wszystkie trzy obrazy (sztywny blok [monolit], przesuwające się bloki [rumowisko], swobodne cząsteczki [obłok pyłu]) współistnieją.

Punkt krytyczny. Wraz ze wzrostem ciśnienia skok ciecz→gaz staje się coraz mniejszy; udział pustek zmienia się coraz bardziej w sposób ciągły. W punkcie krytycznym różnica znika całkowicie — „gęste bloki z niewielką pustką" i „dużo pustki z niewieloma cząsteczkami" przechodzą w siebie w sposób ciągły. („Z punktu widzenia cieczy szklanka jest do połowy pusta, z punktu widzenia gazu do połowy pełna.")

Rekrystalizacja w obrazie 2D. Komiks pokazuje teraz rekrystalizację jako przypadek 2D, jeden pod drugim: na dole gęsto triangulowana sieć (koordynacja z = 6, najniższy energetycznie ideał 2D), na górze luźna sieć kwadratowa (z = 4). Przejście triangulowana → kwadratowa powiększa objętość (praca przeciw ciśnieniu) i obniża liczbę wiązań. Widać, że liczba koordynacyjna się zmienia, a sztywność sieci krystalicznej pozostaje.

Jakościowa prognoza przejść energetycznych i pojemności cieplnych. Z samej łączności z dają się oszacować obie wielkości charakterystyczne: (i) energia wiązania na cząstkę wynosi U/N ≈ −(ε/2)·z (każde wiązanie liczy się połową), a więc tym niższa, im większe z; (ii) skok c_v przy zmianie sieci podlega regule modelowej Δc_v = zmiana z_max/4. Dla przypadku 2D triangulowana (z = 6) ↔ kwadratowa (z = 4) daje to różnicę energii ΔU ≈ ε na cząstkę i stopień c_v równy Δc_v = (6 − 4)/4 = 0,5 k_B — o połowę mniejszy niż wartość 3D (1,0 k_B). Reguła pamięciowa: gęściej = niżej energetycznie, ale drożej w przejściu (większe ciepło utajone, większy stopień c_v), ponieważ przy rozszerzaniu uwalnia się więcej wiązań i więcej zamrożonych stopni swobody.

Czy struktura sieci wyznacza metastabilność? Częściowo. Struktura sieci ustala kolejność energetyczną — która postać jest stanem podstawowym (triangulowana, najniższe μ), a która wzbudzonym (kwadratowa). Czy luźna postać przetrwa w niskiej temperaturze jako metastabilna, nie rozstrzyga ona jednak sama: w ściśle idealnym modelu (bez bariery) substancja przyjmuje przy każdym (τ, p) natychmiast fazę o najniższym μ — metastabilność tam nie istnieje. Powstaje ona dopiero przez skończoną energię aktywacji (barierę zarodkowania), która kinetycznie hamuje przebudowę. Geometria sieci dostarcza więc termodynamicznego uporządkowania, bariera zaś kinetycznej trwałości — ten sam składnik, który tłumaczy też fazę przesyconą (zob. niżej; ilościowo rozdz. 6.7).

Faza przesycona — energia aktywacji. Rysunek dodatkowy (3.2, rząd 2) pokazuje przesyconą fazę gazową nad cieczą: obok pojedynczych cząsteczek małe grupy 2–3 zlepionych cząstek (szare belki), które lecą jako całość. Model tłumaczy przesycenie zauważalną energią aktywacji wiązań molekularnych: małe klastry (2–3) mają niekorzystny stosunek brzegu do powierzchni, ich wiązania brzegowe są niepełne — dopiero zarodek krytyczny rośnie samoistnie dalej. Dopóki nie powstanie zarodek nadkrytyczny, gaz pozostaje przesycony, choć μ_gaz > μ_ciecz już pozwalałoby na kondensację. Właśnie tej bariery brakuje w idealnym obrazie równowagi i zostaje ona tu nazwana jako samodzielny, „poza idealnym obrazem" leżący wkład (oszacowanie ilościowe w rozdz. 6.7).

<a id="kap3-3"></a>

3.3 Dlaczego prędkość dźwięku w cieczach maleje z T

Dźwięk to fala sprężysta: v = √(K_eff/ρ). K_eff powstaje wyłącznie z nienaruszonych wiązań (WIĄZANIE) — to one są sprężynami. Gdy rośnie T, pęka kooperatywnie coraz więcej wiązań (PĘKNIĘCIE); sieć sprężyn rzednie, a K_eff spada stromiej niż gęstość ρ. Zerwanie dominuje licznik, nie mianownik — dlatego v maleje. W ujęciu molekularnym: przenoszenie dźwięku przez wiązania molekularne jest szybsze niż mechaniczne przenoszenie przez zderzenia. Im więcej pęknięć pojawia się w cieczy wraz ze wzrostem temperatury, tym częściej dźwięk musi być niesiony przez zderzenia mechaniczne niż przez wiązania molekularne. Odpowiednio prędkość dźwięku maleje wraz ze wzrostem temperatury. (Anomalię wody — maksimum akustyczne przy ~74 °C — omawia rozdz. 6.4.)

<a id="kap3-4"></a>

3.4 Jak powstaje ruch Browna — z zaostrzonym pojęciem miceli

W sensie statystyczno-termodynamicznym micela nie jest tu agregatem surfaktantu, lecz blokiem otoczonym przez PĘKNIĘCIE: brakujące wiązania (nie-wiązania na granicy bloku) tworzą „powłokę" wokół zamkniętego krystalitu (nienaruszenie związanego rdzenia). Rdzeń ten może wykonywać wewnętrzny kooperatywny ruch obrotowy (TRÓJDŹWIĘK-c) — blok obraca się jako całość, podczas gdy jego powłoka graniczna odsprzęga go od sąsiadów.

Ruch Browna powstaje teraz tak: przez wymianę wiązań molekularnych (PĘKNIĘCIE, RĄCZY) pojedyncze cząsteczki mogą opuścić powłokę jednej miceli i przejść do innej miceli. Każda taka zmiana przenosi pęd (ZACHOWANIE; dzięki SWOBODNEMU LOTOWI bezstratnie) i skokowo zmienia kształt oraz moment pędu obu uczestniczących miceli. Wielka cząstka próbna (ziarno pyłu) jest opływana przez wiele takich nieskorelowanych przebudów miceli; suma popchnięć pędowych jest średnio zerowa, ma jednak skończoną wariancję → błądzenie losowe. Im wyższe T, tym częstsze zdarzenia PĘKNIĘCIA i zmiany miceli, tym silniejsze drżenie. Ruch Browna jest tym samym makroskopowym cieniem, czyli następstwem, mikroskopowego procesu wymiany miceli.

<a id="kap3-5"></a>

3.5 Dlaczego reakcje na warstwach granicznych miceli przyspieszają

Żenująca uwaga wstępna. Poniższe rozumowanie jest wyprowadzeniem z obrazu i wedle mojej wiedzy nie zostało dotąd potwierdzone ani zbadane empirycznie, przy czym muszę przyznać, że nie zadałem sobie trudu, by to dokładnie rozeznać. Jedyną przesłanką nowości jest to, że systemy SI nie wspominają o tym wniosku. Ta teza jest przykładem, jak obrazowe wyobrażenia mogą prowadzić do nowych pytań.

W umiarkowanie rozpuszczających rozpuszczalnikach reagenty osadzają się chętniej na powłokach granicznych miceli — w modelu właśnie na granicach pęknięcie/blok (PĘKNIĘCIE). Dwa efekty przyspieszają reakcję: (i) redukcja wymiaru — reagenty nie szukają się już w całej objętości, lecz wzdłuż quasi-dwuwymiarowej granicy faz; poszukiwanie 2D trafia statystycznie szybciej niż 3D. (ii) Zatężenie — powłoka graniczna zbiera reagenty, ich lokalna gęstość rośnie. Obrazowo granice pęknięć działają jak szyny prowadzące, które doprowadzają partnerów do siebie.

Efekt wymiaru ma ścisły rdzeń: błądzenie losowe (Pólya). Dlaczego poszukiwanie dwuwymiarowe trafia statystycznie lepiej niż trójwymiarowe, wynika nie tylko z poglądowości, lecz z twierdzenia rachunku prawdopodobieństwa. George Pólya (1887–1985) wykazał w 1921 r. dla błądzenia losowego („random walk") na sieci: w jednym i w dwóch wymiarach wędrowiec wraca z prawdopodobieństwem 1 (a więc na pewno) nieskończenie wiele razy do każdego punktu, na sieci płaskiej ℤ² tak samo jak na prostej ℤ. Od trzech wymiarów to już nie obowiązuje: prawdopodobieństwo powrotu jest istotnie mniejsze od 1 (na prostej sieci kubicznej ℤ³ tylko około 0,34). Obrazowo, w ujęciu Pólyi: „Pijany człowiek z pewnością kiedyś trafi do domu, pijana ryba natomiast nie." Przeniesione na dwa reagenty, które dla udanego zderzenia muszą się spotkać wielokrotnie, znaczy to: w bardzo dobrze rozpuszczającym rozpuszczalniku ze swobodnym trójwymiarowym błądzeniem raz rozdzieleni partnerzy z niezerowym prawdopodobieństwem oddalają się na zawsze — powtarzane próby reakcji zawodzą, a wydajność pozostaje z tego czysto geometryczno-statystycznego powodu niska. Na (quasi-)dwuwymiarowej warstwie granicznej miceli natomiast ponowne spotkanie jest pewne (prawdopodobieństwo 1): partnerzy odnajdują się wciąż na nowo, aż reakcja się powiedzie. Warstwa graniczna działa więc nie tylko zatężająco, lecz zarazem udostępnia korzystniejszy dla ponownych spotkań wymiar. Twierdzenie (świadomie zaostrzone) brzmi: wydajności w bardzo dobrze rozpuszczających rozpuszczalnikach są z przyczyn błądzenia losowego z reguły gorsze.

Odsyłacz: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — tu cytowane za M. Eisermann, Höhere Mathematik 3 (Uniw. Stuttgart), zadanie W-2×2. (Oryginał odwiedzony 5.9.2026).

<a id="kap3-6"></a>

3.6 Trójwymiarowe okno cieczy

Inaczej niż rysunki w 3.1–3.5 ten nie pokazuje żadnej sceny molekularnej, lecz wynik obliczeń z rozdz. 7.9/7.10 jako schemat na osi τ (w jednostkach energii wiązania ε). Celem było uwidocznić, dlaczego pojedyncza sieć 2D ledwie dopuszcza obszar dla fazy ciekłej, podczas gdy przy sieci 3D jest to bez trudu możliwe, przy czym w obu przypadkach obowiązuje ta sama reguła sztywności (z̄_c = 2d).

Obraz ① — 2D (brak okna). Na osi τ τ_m = 1,20 i τ_triple = 1,16 leżą praktycznie na sobie (τ_triple poniżej τ_m), τ_c = 1,50 tuż powyżej: obszar fazy ciekłej jest ściśnięty do szczeliny — ciało stałe przechodzi niemal wprost w gaz (blisko sublimacji).

Obraz ② — 3D (czyste okno). Ta sama oś, ale rozciągnięta: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Między τ_triple a τ_c leży samodzielne pole ciekłe (szerokie ku dołowi, τ_triple − τ_m = 1,63; wąskie ku górze, τ_c − τ_triple = 0,18). Wąski górny brzeg jest jako taki zaznaczony.

Obraz ③ — przyczyna. z̄_c = 2d rośnie z wymiarem (4 → 6), przez co v_solid maleje (0,50 → 0,667 w objętościach zredukowanych, tzn. ciało stałe staje się „luźniejsze" względem cieczy), a τ_c = z_max/4 podwaja się (1,50 → 3,00). Punkt krytyczny przesuwa się daleko powyżej punktu topnienia. Dopiero przez to powstaje miejsce dla fazy ciekłej. To czysty efekt wymiaru.

<a id="kap4"></a>

4 — Model statystyczno-mechaniczny (bez zmian z wersji 09)

Podsumowanie

Rozdział 4 buduje rusztowanie statystyczno-mechaniczne modelu w 2D (sieć trójkątna, z_max = 6):

• 4.1 Stopnie swobody i zmienne: wielkością porządku jest udział wiązań x ∈ [0,1] ze średnią liczbą koordynacyjną z̄ = z_max·x; na blok ≥ 2 dochodzi stopień swobody rotacji, a temperatura τ liczy się w jednostkach energii wiązania ε.
• 4.2 Kooperatywne wiązanie dwustanowe: samouzgodnione równanie wiązania x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) wynika jako warunek stacjonarności energii swobodnej Bragga–Williamsa; jej minimum globalne daje gałąź fizyczną, a rozdz. 5.5 rozstrzyga, czy biegnie ona w sposób ciągły, czy ze skokiem (pierwszego rodzaju).
• 4.3 Bilans energetyczny z członem kooperatywnym: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) łączy translację, drganie wiązania i rotację bloku, energię konfiguracyjną wiązań oraz nową kooperatywną energię pola molekularnego −(K/2)x². Odpowiednia pojemność cieplna c_v otrzymuje dodatni człon dodatkowy −K·x·ẋ, który wyostrza pik topnienia, ale pozostawia nietknięte plateau ciała stałego i gałąź gazową.

<a id="kap4-1"></a>

4.1 Stopnie swobody i zmienne

• N cząsteczek na stałej referencyjnej geometrii koordynacyjnej. Kanoniczny przypadek 2D: sieć trójkątna, z_max = 6, jeden stopień swobody rotacji na blok (WIELOŚĆ, WIĄZANIE). Rozdział ten obowiązuje wyraźnie dla interpretacji 2D; uogólnienie 3D znajduje się w rozdz. 6.
• Zmienna wiązania s_ij ∈ {0, 1} na każdą geometryczną parę sąsiadów (1 = nienaruszone).
• Liczba koordynacyjna z_i = Σ_j s_ij; średnia z̄ = 2 N_bond / N — parametr porządku. Zachodzi z̄ = z_max · x z udziałem wiązań x = p_b ∈ [0, 1].
• Jednostka energii τ ≡ k_B T (klasycznie, ½τ na kwadratowy stopień swobody). Wszystkie energie w jednostkach energii wiązania ε.

<a id="kap4-2"></a>

4.2 Kooperatywne wiązanie dwustanowe (typu Isinga, parametr J)

Izolowane wiązanie (energia −ε nienaruszone, 0 zerwane) ma w równowadze Maxwella–Boltzmanna p_b⁰(τ) = 1/(1 + e^(−ε/τ)). Kooperatywne PĘKNIĘCIE (aksjomat 6) przesuwa efektywną stabilność wiązania o pole molekularne proporcjonalne do udziału wiązań x. W przybliżeniu pola molekularnego (Isinga) samouzgodnione równanie wiązania przybiera postać:

x = 1 / (1 + G · exp[−(ε + K·x) / τ]), K = J·z_max, z̄ = z_max·x.

Przy tym:

• K = J·z_max to kooperatywna siła sprzężenia (sprzężenie zwrotne PĘKNIĘCIA: więcej wiązań sąsiedzkich → wyższa efektywna stabilność wiązania). J > 0 wyostrza przejście; powyżej krytycznego J z̄(τ) staje się nieciągłe — to jest ciepło utajone topnienia. Czysty przypadek dwustanowy to granica J = 0 (brak pola molekularnego, ciągłe przejście).
• G > 1 to zwyrodnienie stanu zerwanego: uwolniona cząsteczka zyskuje przestrzeń fazową (swobodę translacji). G napędza zrywanie wiązań przy wysokim τ i sprawia, że x → 1/(1+G), zamiast zatrzymać się na ½ (granica gazu). Kanonicznie G = z_max.

Równanie to wynika jako warunek stacjonarności energii swobodnej Bragga–Williamsa na cząsteczkę (n_slots = z_max/2 miejsc wiązania):

g(x)/n_slots = −ε·x − (K/2)·x² + τ·[x ln x + (1−x) ln(1−x)] − τ·ln(G)·(1−x).

Człon −(K/2)x² to sprzężenie kooperatywne (ferromagnetyczne), ostatni człon — zwyrodnienie zerwania. Minimum globalne g po x daje gałąź fizyczną; tam, gdzie dwa minima stają się jednakowo głębokie, leży przejście pierwszego rodzaju. RĄCZY (aksjomat 7) uzasadnia, że x(τ) przyjmuje w każdej temperaturze wartość równowagową.

<a id="kap4-3"></a>

4.3 Bilans energetyczny — z członem kooperatywnym

Człon pola molekularnego −(K/2)x² jest prawdziwą energią (niezależną od τ) i musi występować w U/N. Spójny bilans brzmi zatem (2D):

U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x².

• τ·[…]: translacja (1) + drganie wiązania (z̄/4) + rotacja bloku (½ n_blk) — TRÓJDŹWIĘK.
• −(ε/2)·z̄: energia konfiguracyjna wiązań — WIĄZANIE.
• −(K/2)·x²: kooperatywna energia pola molekularnego — PĘKNIĘCIE. Nowe w U/N.

Odpowiednia pojemność cieplna c_v = d(U/N)/dτ (w jednostkach k_B) otrzymuje dodatkowy wkład:

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ ż̄ + ½ ṅ_blk] − (ε/2)·ż̄ − K·x·ẋ, (ẋ = dx/dτ)

gdzie nowy człon −K·x·ẋ jest dodatni (ẋ < 0) i działa dokładnie tam, gdzie wiele wiązań pęka jednocześnie. Wyostrza on pik topnienia, nie zmieniając plateau ciała stałego (ẋ ≈ 0) ani gałęzi gazowej.

<a id="kap5"></a>

5 — Sumy stanów i pojemności cieplne

Podsumowanie

Rozdział 5 oblicza sumy stanów i pojemności cieplne oraz dostarcza danych modelowych do kalibracji (rozdz. 7):

• 5.1 Rozłożona suma stanów: wskutek rozdzielenia skal czasowych Z = Z_trans·Z_vib·Z_rot·Z_bond rozkłada się na czynniki.
• 5.2 Pojemność cieplna dla każdej fazy: c_v = d(U/N)/dτ; dodatni człon −(ε/2)(dz̄/dτ) rodzi pik typu Schottky'ego jako odpowiednik ciepła utajonego.
• 5.3 Numeryczna krzywa c_v(τ): pik c_v 5,51 przy τ_m = 1,20, plateau ciała stałego 2,50, gałąź gazowa → 1; ciepło utajone L ≈ 1,41 ε na cząsteczkę (z członem kooperatywnym) jest wobec 0,48 ε (bez niego) potrojone.
• 5.4 / 5.4a Dwa rozdzielone progi: perkolacja łączności z̄_c = z_max·p_c ≈ 2,08 (przy τ_b ≈ 1,85) znaczy rozpad klastrów (blisko wrzenia) i jest wynikiem; należy ją czysto oddzielić od generycznej sztywności Maxwella z̄_c = 2d (gęstość topnienia, utrata sztywności).
• 5.5 Rodzaj przejścia topnienia: kryterium ln G > 2 + 4ε/K; 2D (G = 6) jest dla każdego J ciągłe, 3D (G = 12) pierwszego rodzaju od K > ~8,2 (demonstracja: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Studium skończonego rozmiaru: p*(L) zbiega monotonicznie od dołu do p_c = 0,3473; z̄_c ≈ 2,08 jest odporne na skończony rozmiar.
• 5.7 Pełny skan 3D c_v: od K > K_kryt ≈ 8,25 (J > 0,68) otwiera się prawdziwy skok Δz̄ z ciepłem utajonym ΔU; pik topnienia wędruje z J ku wyższym τ (przesunięcie dodatnie +0,5 … +2,2 wobec ciągłego piku 2D).

<a id="kap5-1"></a>

5.1 Rozłożona suma stanów

Wskutek rozdzielenia skal czasowych (RĄCZY) kanoniczna suma stanów rozkłada się na czynniki:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (klasyczno-idealna translacja 2D; WIELOŚĆ, ZACHOWANIE, SWOBODNY LOT).
• Z_vib = (τ/ħω₀)^(N_bond) (klasyczny oscylator harmoniczny na nienaruszone wiązanie; TRÓJDŹWIĘK-b).
• Z_rot = ∏_(bloki≥2) (2π I_b τ)^(1/2)/ħ (klasyczny sztywny rotator 2D na blok; TRÓJDŹWIĘK-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·człon kooperatywny] (suma wiązań typu Isinga; WIĄZANIE, PĘKNIĘCIE). Ocena w polu molekularnym daje x(τ), z̄(τ) z 4.2.

<a id="kap5-2"></a>

5.2 Pojemność cieplna dla każdej fazy (analitycznie)

Z U/N (rozdz. 4.3), c_v = d(U/N)/dτ (w jednostkach k_B):

c_v = [1 + z̄/4 + ½ n_blk] + τ·[¼ (dz̄/dτ) + ½ (dn_blk/dτ)] − (ε/2)·(dz̄/dτ).

Człon −(ε/2)·(dz̄/dτ) jest dodatni (dz̄/dτ < 0: ogrzewanie zrywa wiązania) i dominuje tam, gdzie wiele wiązań pęka jednocześnie — pik typu Schottky'ego, odpowiednik ciepła utajonego.

<a id="kap5-3"></a>

5.3 Numeryczna krzywa c_v(τ) z członem kooperatywnym — ciepło utajone skwantyfikowane

Wiązanie samouzgodnione (minimum globalne), n_blk(τ) z Monte Carlo perkolacji wiązań (union-find, L = 100, 5 realizacji), c_v = dU/dτ (różnice centralne), U/N wraz z −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Wybrane punkty podparcia:

Wynik (numerycznie). Porządek trzech faz pozostaje potwierdzony i jest teraz ostrzejszy:

c_v(stałe) ≈ 2,50 < c_v(ciecz, maksimum) ≈ 5,51 > c_v(gałąź gazowa) ↓ → 1.

Ciepło utajone (pole piku ponad linią bazową). Całkując (c_v − linia bazowa) po oknie topnienia τ ∈ [0,70 ; 1,70] (linia bazowa = prosta przez c_v(0,70) i c_v(1,70)), otrzymujemy

L ≈ 1,41 ε na cząsteczkę (z członem kooperatywnym), L ≈ 0,48 ε na cząsteczkę (bez członu kooperatywnego, dawny bilans), ΔL ≈ +0,92 ε na cząsteczkę.

Człon kooperatywny podnosi pik o czynnik ≈ 1,7 (3,25 → 5,51) i potraja ciepło utajone. Fizycznie: jednoczesne zrywanie wielu wiązań sąsiedzkich (PĘKNIĘCIE) magazynuje nie tylko energię konfiguracyjną (−ε/2·z̄), lecz także energię pola molekularnego (−K/2·x²) — obie uwalniają się przy topnieniu.

<a id="kap5-4"></a>

5.4 Próg perkolacji łączności jako wynik

Udział olbrzymiego klastra P_∞(τ) załamuje się przy τ ≈ 1,85; tam x ≈ 0,347 = p_c(trójkąt) = 2·sin(π/18) = 0,3473, a więc z̄_c = z_max·p_c ≈ 2,08 — wynik, nie wartość narzucona. Jest to próg łączności (klaster przechodzący) i znaczy rozpad klastrów blisko wrzenia (τ_b); nie należy go mylić ze sztywnością Maxwella (z̄_c = 2d, topnienie) — por. czyste rozdzielenie w rozdz. 5.4a.

<a id="kap5-4a"></a>

5.4.a Dwa różne progi czysto rozdzielone

Dotąd nieostre podwójne nazewnictwo „perkolacja/sztywność" zostaje teraz rozwiązane. Chodzi o dwa fizycznie różne progi, których nie wolno mylić:

Uzasadnienie. Topnienie to utrata mechanicznej sztywności gęstego bloku, a nie rozpad spójności. Wielkością istotną jest więc zliczanie więzów Maxwella (pole średnie / sztywność generyczna): sieć D-wymiarowa sił centralnych staje się izostatyczna dokładnie wtedy, gdy średnia koordynacja osiąga z̄ = 2D (2D → 4, 3D → 6). Łączność wiązań (z̄_c = 2,08) leży wyraźnie poniżej; sieć ledwie spójna przy z̄ = 2,08 jest od dawna wiotka (mody miękkie), a więc już ciekła. Dlatego łączność znaczy późniejszy rozpad klastrów blisko wrzenia, podczas gdy sztywność ustala wcześniejsze krzepnięcie/topnienie (gęstość topnienia). To rozdzielenie ustala punkt potrójny w rozdz. 7.8 bezparametrowo i realistycznie; nie zmienia potwierdzonego w 5.4/5.6 progu łączności 2,08, lecz tylko jego odczyt (τ_b = blisko wrzenia, nie topnienie).

<a id="kap5-5"></a>

5.5 Rodzaj przejścia topnienia

Kryterium analityczne. Obszar spinodalny (dwa minima g) istnieje dla τ < K/4; przejście leży w przybliżeniu przy τ_t ≈ (ε + K/2)/ln G. Prawdziwy skok (przejście pierwszego rodzaju) wymaga τ_t < K/4, a więc w granicy K → ∞:

ln G > 2 + 4ε/K.

Wynika stąd godny uwagi rezultat co do zwyrodnienia zerwania G = z_max:

Potwierdzenie numeryczne (maks. |Δx| na krok τ, detect_first_order):

Odczyt. W kanonicznym modelu 2D (G = 6) topnienie jest z zasady ciągłe („miękkie" przejście / crossover) — dla żadnego J nie powstaje prawdziwy skok, ponieważ ln 6 < 2. Dopiero wyższa koordynacja/zwyrodnienie upakowania 3D (G = 12) pozwala od dostatecznej kooperatywności (K > 8,2) na przejście pierwszego rodzaju ze skończonym ciepłem utajonym. Pokrywa się to z rzeczywistością: topnienie 2D jest często ciągłe (typu KTHNY), topnienie 3D jest pierwszego rodzaju. J = 0 (brak pola molekularnego) to zawsze ciągły przypadek odniesienia. — Fizycznie J oznacza stosunek energii kooperatywnej do energii wiązania (K/ε = J·z_max/ε): jak silnie nienaruszone wiązanie sąsiedzkie stabilizuje własne.

Demonstracja pierwszego rodzaju (motywowana 3D, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ skacze o Δz̄ ≈ 5,30 przy τ ≈ 2,82; ciepło utajone jako skok energii ΔU ≈ 1,75 ε na cząsteczkę. Ten sam K na sieci G = 6 (2D) pozostaje ciągły — charakter rodzaju zależy wyłącznie od G = z_max.

<a id="kap5-6"></a>

5.6 Studium skończonego rozmiaru progu perkolacji

Perkolacja wiązań na sieciach trójkątnych L × L (periodycznych), udział olbrzymiego klastra P_∞(p) przeskalowany po L (cv_scan_09.py; L = 50/100/200/400 z 8/5/3/2 realizacjami). Dokładny próg p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (średnia ± rozrzut):

Oszacowanie progu p*(L) (P_∞ przecina 0,5):

Wynik. p*(L) zbiega monotonicznie od dołu do p_c = 0,3473; rozrzut przy p_c spada z ±0,103 (L = 50) do ±0,017 (L = 400) — przejście wyostrza się z L, ale nie przesuwa progu (uniwersalna wartość perkolacyjna). Tym samym z̄_c ≈ 2,08 = z_max·p_c zostaje potwierdzone jako odporna na skończony rozmiar wartość modelowa.

<a id="kap5-7"></a>

5.7 Pełny skan 3D c_v(τ) topnienia pierwszego rodzaju

Pojedynczy przebieg demonstracyjny pierwszego rodzaju z rozdz. 5.5 (J = 1,0, G = 12) zostaje teraz rozbudowany do pełnego skanu c_v(τ) po kooperatywności J (firstorder_3d_16.py, łączność FCC dla n_blk, bilans energii 3D U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Analityczne kryterium rodzaju z rozdz. 5.5 (pierwszego rodzaju ⇔ K > K_kryt = 4ε/(ln 12 − 2) = 8,25) rozdziela przypadki; ciągły pik odniesienia 2D leży przy τ_2D = 1,20 (c_v ≈ 5,5).

*przypadek ciągły: Δz̄/ΔU to strome przejście (nie prawdziwy skok); τ_przejścia jest tam pikiem c_v.

Wynik. (i) Ciągły przypadek 3D J = 0,55 daje τ_przejścia = τ_m ≈ 1,72 — próba krzyżowa z niezależną kinetyką wiązań (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) zgadza się. (ii) Od K > K_kryt ≈ 8,25 (J > 0,68) otwiera się prawdziwy skok Δz̄ ze skończonym ciepłem utajonym ΔU (skok energii przez nieciągłość) — topnienie staje się pierwszego rodzaju, dokładnie jak wymaga kryterium ln 12 > 2. (iii) Próba krzyżowa J = 1,0 odtwarza wartość demonstracyjną Δz̄ ≈ 5,4 z rozdz. 5.5. (iv) Wysokość skoku Δz̄ i położenie τ_przejścia rosną z J: przejście wędruje ku górze — wyższa koordynacja upakowania 3D oraz silniejsza kooperatywność podnoszą temperaturę topnienia. Wobec ciągłego piku 2D (τ_2D = 1,20) wynika stąd wszędzie dodatnie przesunięcie piku (+0,5 … +2,2), rosnące monotonicznie z K. Charakter pierwszego rodzaju topnienia 3D jest tym samym udokumentowany nie tylko w jednym punkcie, lecz w całym zakresie kooperatywności.

<a id="kap6"></a>

6 — Matematyczna wiarygodność interpretacji 3D

Podsumowanie

Rozdział 6 czyni uogólnienie obrazu wiązań z 2D na 3D wiarygodnym czysto matematycznie (bez własnej wizualizacji HTML):

• 6.1 Stopień swobody rotacji 2D → 3D: rotacja bloku rośnie z jednej do trzech osi, translacja z dwóch do trzech kierunków, z_max z 6 do 12; wartości bazowe c_v rosną (c_v,stałe 2,5 → 4,5, c_v,gaz 1 → 3/2), kształt krzywej c_v(τ) pozostaje.
• 6.2 Łączność jako cecha odróżniająca: z_max jest wielkością skwantowaną i skokowo oddziela najgęstsze upakowanie (FCC/HCP, z_max = 12) od BCC (z_max = 8), bez stanu pośredniego.
• 6.3 Dodatkowe przejście stałe(1)→stałe(2): z dyskretności wynika polimorfizm FCC↔BCC z bezparametrowym stopniem c_v 1,0 k_B; w ścisłym obrazie idealnym (θ_BCC = θ_FCC) nie występuje ono, a dopiero bardziej miękkie BCC (θ_BCC < θ_FCC, poniżej progu otwarcia 0,61) otwiera τ_ss < τ_m.
• 6.4 / 6.5 Rozważania: anomalia akustyczna ciekłej wody jako rozmyta zmiana uporządkowania bliskiego zasięgu między dwoma typami sieci oraz fraktalna alternatywa pól wirowych wobec rotacji miceli (czysto jakościowo).
• 6.6 / 6.8 / 6.9 Mikroskopowe θ(sieć) na trzech poziomach: model sprężyn Einsteina (θ_BCC/θ_FCC ≈ 0,82…0,94), gęstość stanów fononowych typu Debye'a (≈ 0,87…1,00) i pełna sieć Lennarda-Jonesa (≈ 1,0; idealne BCC niestabilne na ścinanie, FCC stabilny stan podstawowy) — wszystkie powyżej progu 0,61, rzeczywisty polimorfizm leży więc poza harmonicznym obrazem idealnym.
• 6.7 Metastabilność (2D): rachunek rzędu wielkości dla przebudowy triangulowana ↔ kwadratowa (ΔU = +1 ε, ΔA ≈ +15,5 %) wraz z barierą zarodkowania ΔG\* = πλ²/Δμ_A jako dodatkiem kinetycznym.

Rozdział ten czyni uogólnienie 3D wiarygodnym czysto matematycznie. 3D nigdy nie jest wizualizowane w wersji HTML. Celem jest jedynie pokazać, jak stopień swobody rotacji i łączność zmieniają się z 2D na 3D i jakie nowe przejścia stają się przez to wyobrażalne.

<a id="kap6-1"></a>

6.1 Zmiana stopnia swobody rotacji 2D → 3D

Jedynym jakościowo nowym punktem wobec 2D jest rotacja bloków:

• 2D: sztywny blok ma 1 stopień swobody rotacji (obrót w płaszczyźnie) → U_rot/N = ½ n_blk τ.
• 3D: sztywny blok ma 3 stopnie swobody rotacji (trzy osie Eulera) → U_rot/N = (3/2) n_blk τ.

Zarazem translacja rośnie z 2 do 3 stopni swobody (U_trans/N = (3/2)τ), a maksymalna liczba koordynacyjna z z_max = 6 (trójkąt) do z_max = 12 (najgęstsze upakowanie kul). Bilans energii w 3D brzmi zatem

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄,

a wartości bazowe pojemności cieplnej przesuwają się spójnie:

Struktura krzywej (plateau ciała stałego, maksimum cieczy przez −(ε/2)dz̄/dτ, spadek do wartości gazowej) pozostaje identyczna; skalują się z wymiarem i z_max jedynie poziomy.

<a id="kap6-2"></a>

6.2 Wpływ różnych sieci krystalicznych — łączność jako cecha odróżniająca

W 3D należy rozróżnić co najmniej dwie idealne sieci, które rozdziela ich łączność (liczba koordynacyjna z_max):

Łączność z_max jest rozstrzygającą cechą odróżniającą: ustala zarówno poziom pojemności cieplnej ciała stałego (z_max/4), jak i próg perkolacji z̄_c. Ponieważ każda sieć dopuszcza tylko jedną stałą liczbę koordynacyjną, łączność jest wielkością skwantowaną i skończoną — nie ma ciągłego przejścia między z_max = 12 a z_max = 8, lecz tylko skok dyskretny z jednego typu sieci do drugiego.

<a id="kap6-3"></a>

6.3 Dodatkowe przejście stałe(1) → stałe(2)

Ze skończoności/dyskretności łączności wynika nieobecne w modelu 2D przejście: stałe(1) → stałe(2) (polimorfizm/alotropia). Dwie fazy stałe o różnym z_max (np. FCC ↔ BCC) mają różne entalpie swobodne μ(τ, p); przy kombinacji temperatury/ciśnienia μ_{stałe(1)} = μ_{stałe(2)} substancja zmienia swoją sieć bez topnienia. W obrazie modelowym zmienia się przy tym geometria koordynacyjna skokowo (nie tylko obsadzenie z̄ jak przy topnieniu) — jedyna sytuacja, w której skądinąd utrzymywana geometria referencyjna sama przeskakuje. Obserwowalny byłby mały stopień c_v (zmiana z_max/4) plus utajone ciepło przemiany, wyraźnie słabsze niż pik topnienia.

Kryterium pierwszego rodzaju ln G > 2 + 4ε/K (rozdz. 5.5) dostarcza teraz konkretnego powodu, dlaczego właśnie interpretacja 3D jest jakościowo bogatsza: tylko tam (G = z_max = 12, ln 12 = 2,485 > 2) topnienie może stać się pierwszego rodzaju — i tylko tam istnieje dyskretna zmiana sieci stałe(1) → stałe(2) (FCC ↔ BCC), ponieważ łączność jest skwantowana/skończona. Oba zjawiska to ta sama przyczyna w dwóch szatach: wyższa, dyskretna koordynacja upakowania 3D rodzi zarówno ostry skok topnienia, jak i polimorfizm. Pojedyncza sieć 2D (G = 6) nie ma żadnego z nich — zgodnie z rozdz. 5.5.

Kwantyfikacja (bondkinetics_3d_15.py). Przejście liczy się teraz oszczędnie. Obie sieci otrzymują po jednej entalpii swobodnej w klasycznym obrazie harmonicznym (wspólne części p·v i translacyjne wypadają w porównaniu):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (sztywność fononowa w jednostkach ε),

z bezparametrowymi wartościami charakterystycznymi z rozdz. 6.1/6.2:

Stąd stopień c_v bezparametrowo: |Δc_v| = zmiana z_max/4 = (12−8)/4 = 1,0 k_B (przy przejściu FCC→BCC ciało stałe staje się o 1 k_B „lżejsze" w pojemności cieplnej) — o wiele mniej niż pik topnienia (kilka k_B, rozdz. 5.3/7.11).

Temperatura przemiany τ_ss (μ_FCC = μ_BCC) oraz utajone ciepło przemiany L = τ_ss·(S_BCC − S_FCC):

Wynik (uczciwe ustalenie). W ścisłym obrazie idealnym (ta sama sztywność fononowa dla wszystkich wiązań, θ_BCC = θ_FCC) nie ma przejścia stałe(1)→stałe(2): FCC wygrywa zarówno entalpię (niższa, −6 wobec −4 ε), jak i — z powodu wyższej liczby wiązań (c_v = 4,5 > 3,5) — entropię drgań; μ_FCC pozostaje w każdej temperaturze poniżej μ_BCC. Dopiero zależna od sieci sztywność fononowa θ_BCC < θ_FCC (bardziej otwarte BCC drga miękcej, ma więc więcej entropii) sprowadza μ_BCC przy wysokiej temperaturze poniżej μ_FCC — wtedy przemiana się otwiera i leży poniżej piku topnienia τ_m, tak że powstaje kolejność stałe(1)→stałe(2)→ciecz→gaz (τ_ss < τ_m < τ_c ≤ τ_b(3D), rozdz. 7.11). Miękka sztywność BCC to łagodne rozszerzenie poza ścisłym obrazem idealnym (znane w rzeczywistości jako polimorfizm γ/δ żelaza: BCC w wysokiej temperaturze, FCC pośrodku); stopień c_v 1,0 k_B i małe dodatnie ciepło przemiany są natomiast wewnątrzmodelowe i bezparametrowe.

Jak małe jest naprawdę θ_BCC/θ_FCC? (theta_springs_16.py). Aby uwolnić przemianę od swobodnego parametru θ, sztywność fononowa jest teraz szacowana mikroskopowo (rozdz. 6.6): model sprężyn Einsteina z centralnymi sprężynami parowymi daje θ_BCC/θ_FCC ≈ 0,82 (tylko najbliżsi sąsiedzi) do ≈ 0,94 (z drugą powłoką BCC). Wymagany z porównania μ próg otwarcia leży natomiast przy θ_BCC/θ_FCC < 0,61 (poniżej niego τ_ss istnieje i leży pod τ_m). Geometryczny stosunek sztywności jest więc za duży: w ścisłym obrazie idealnym FCC pozostaje stabilne w każdej temperaturze, nie ma czysto geometrycznie napędzanego przejścia stałe(1)→stałe(2). Rzeczywisty polimorfizm (żelazo γ↔δ) wymaga dodatkowych (elektronowych/magnetycznych/anharmonicznych) wkładów poza idealnym obrazem harmonicznym — dokładnie „łagodnego rozszerzenia" założonego powyżej. Potwierdza to uczciwe ustalenie i zmienia stosunek θ ze swobodnego założenia w oszacowaną, nazwaną granicę modelu.

<a id="kap6-4"></a>

6.4 Rozważania: anomalia akustyczna ciekłej wody

Normalne ciecze wykazują v(T) monotonicznie malejące (rozdz. 3.3). Woda ma maksimum akustyczne przy ≈ 74 °C: v najpierw rośnie z T, zanim zacznie maleć. W modelu daje się to odczytać jako konkurencję dwóch lokalnych typów sieci (rozdz. 6.2): otwartego, nisko skoordynowanego uporządkowania bliskiego zasięgu podobnego do lodu (małe z_max, wiele „dziur") i gęstego, wyżej skoordynowanego uporządkowania bliskiego zasięgu. Przy niskim T dominuje struktura otwarta; ogrzewanie najpierw łamie otwarte uporządkowanie na rzecz gęstszego — efektywna łączność, a z nią K_eff rosną (v rośnie). Dopiero powyżej maksimum przeważa zwykłe zrywanie wiązań PĘKNIĘCIA, K_eff spada, v maleje. Anomalia jest tym samym wyrazem rozmytej, niepełnej zmiany uporządkowania bliskiego zasięgu typu stałe(1)↔stałe(2) wewnątrz fazy ciekłej — wyobrażalnej tylko w modelu 3D z dwoma typami sieci, nie w pojedynczej sieci 2D.

<a id="kap6-5"></a>

6.5 Alternatywa: fraktalne pola wirowe zamiast rotacji miceli

Lokalny ruch Browna (rozdz. 3.4) wyjaśniono przez rotację i wymianę miceli. Alternatywnie daje się go wyobrazić jako efekt fraktalnych pól wirowych: nieustannie przebudowywane przez PĘKNIĘCIE granice bloków rodzą pole lokalnych strumieni pędu, które — podobnie do przepływów turbulentnych — wykazuje samopodobne wiry na wszystkich skalach długości. Cząstka próbna jest wtedy napędzana nie przez dyskretne zderzenia miceli, lecz przez niezmiennicze skalowo kontinuum wirów. Oba obrazy (dyskretne micele ↔ fraktalne wiry) są komplementarnymi przypadkami granicznymi tej samej kinetyki PĘKNIĘCIA; który przeważa, zależy od stosunku rozmiaru cząstki próbnej do średniego rozmiaru miceli. Ta krótka dyskusja jest czysto jakościowa i nie jest tu rozwinięta matematycznie.

<a id="kap6-6"></a>

6.6 Mikroskopowa sztywność fononowa θ(sieć): model sprężyn Einsteina

Porównanie stałe(1)→stałe(2) w rozdz. 6.3 zależy od względnej sztywności fononowej θ_BCC/θ_FCC. Zamiast wybierać ją swobodnie, szacuje się ją tu z koordynacji i geometrii (theta_springs_16.py).

Model. Atom centralny siedzi w harmonicznej misie swoich z sąsiadów; każde wiązanie jest promieniową sprężyną parową κ. Wychylenie u rzutuje się izotropowo z ⟨cos²⟩ = 1/d = 1/3 (d = 3) na kierunki wiązań, tak że efektywna stała sprężystości wynosi k_eff = (z/3)·κ. Częstość Einsteina to ω = √(k_eff/m), a więc θ = ħω/k_B ∝ √k_eff. Przy jednakowej masie i jednakowej sprężynie parowej wynika

θ_BCC/θ_FCC = √(z_BCC/z_FCC) = √(8/12) = 0,816.

Wysubtelnienia. (B) Druga powłoka BCC (6 atomów przy r₂/r₁ = 2/√3 = 1,155) ze sprężyną skalowaną odległością κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, w duchu Lennarda-Jonesa) podnosi stosunek do 0,91…0,94 (BCC staje się przez bliską drugą powłokę efektywnie sztywniejsze). (C) Krótsze wiązanie BCC przy tej samej gęstości (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) lekko usztywnia sprężynę i daje 0,89…0,91. Łącznie okno geometryczne θ_BCC/θ_FCC ≈ 0,82…0,94.

Próg otwarcia. Porównanie μ (rozdz. 6.3) otwiera przemianę (τ_ss istnieje i leży pod τ_m = 1,70) dopiero dla θ_BCC/θ_FCC < 0,61. Okno geometryczne leży wyraźnie powyżej.

Ustalenie. Czysty model sprężynowy nie czyni BCC dostatecznie miękkim: w ścisłym obrazie idealnym FCC pozostaje stabilne, geometrycznie napędzane przejście stałe(1)→stałe(2) nie powstaje. Tym samym θ_BCC/θ_FCC ze swobodnego założenia stało się wielkością oszacowaną, której naturalna wartość (≈ 0,8–0,9) utrzymuje przemianę zamkniętą — zgodnie z uczciwym ustaleniem z rozdz. 6.3. Rzeczywisty polimorfizm leży zatem dowodnie poza idealnym obrazem harmonicznym (udziały elektronowe/magnetyczne, anharmoniczność), a nie w samej geometrii sieci.

<a id="kap6-7"></a>

6.7 Matematyczne oszacowanie metastabilności (2D triangulowana ↔ kwadratowa)

Do jakościowej dyskusji z rozdz. 3.2 tu odpowiedni rachunek rzędu wielkości w obrazie 2D (twarde krążki, średnica D, odległość kontaktu a = D).

(1) Energia. Na cząstkę U/N = −(ε/2)·z (każde wiązanie połową): triangulowana z = 6 → −3 ε, kwadratowa z = 4 → −2 ε. Przebudowa triangulowana → kwadratowa kosztuje więc ΔU = +1 ε na cząstkę (endotermicznie; postać luźna jest bogatsza energetycznie).

(2) Objętość (powierzchnia) i praca. Powierzchnia na cząstkę: triangulowana A_tri = (√3/2)·D² ≈ 0,866 D², kwadratowa A_quad = D². Rozszerzenie wynosi ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — pokazany w komiksie przyrost objętości. Przeciw (zredukowanemu) ciśnieniu p trzeba na to włożyć pracę W = p·ΔA; entalpia swobodna rośnie o Δμ = ΔU + p·ΔA − τ·ΔS. Człon entropowy ΔS > 0 (postać luźna ma bardziej miękkie mody, większy stopień c_v Δc_v = (6−4)/4 = 0,5 k_B) obniża μ_kwad poniżej μ_tri dopiero przy wysokim τ — dlatego triangulowana jest niskotemperaturowym stanem podstawowym, a kwadratowa postacią wysokotemperaturową/metastabilną.

(3) Bariera zarodkowania (2D). Metastabilność zakłada barierę. W 2D zarodek nowej fazy (promień r) rośnie z zyskiem powierzchni −Δμ_A·π r² (Δμ_A = siła napędowa na powierzchnię) przeciw wkładowi brzegu +λ·2π r (λ = energia krawędzi na długość): ΔG(r) = −Δμ_A·π r² + 2π λ r. Maksimum przy r\* = λ / Δμ_A z barierą ΔG\* = π λ² / Δμ_A. Zarodek krytyczny r\ maleje wraz z rosnącą siłą napędową Δμ_A; dopóki skala termiczna τ = k_BT nie osiągnie bariery (**ΔG\ ≳ kilka k_BT), postać luźna (a przy wrzeniu — przesycona) pozostaje kinetycznie zamrożona. Dla małych klastrów (2–3 cząstki, jak w pasku o przesyceniu) stosunek brzegu do powierzchni jest maksymalnie niekorzystny — leżą one na lewo od r\* i średnio rozpuszczają się z powrotem. Tym samym pokazana w komiksie energia aktywacji to ta sama wielkość ΔG\*, która niesie zarówno metastabilną postać krystaliczną, jak i przesyconą fazę gazową; jest ona kinetycznym** dodatkiem do idealnego (bezbarierowego) modelu równowagi.

<a id="kap6-8"></a>

6.8 Wysubtelnienie θ(sieć) w duchu Debye'a

Model sprężyn Einsteina (rozdz. 6.6) zastępuje cały rozkład fononów jedną częstością. Niniejsze podejście wysubtelnia θ(sieć) przez oszczędną gęstość stanów fononowych dla każdej sieci (theta_debye_17.py): te same centralne sprężyny parowe κ(r) = r^(−n) na najbliższej i drugiej powłoce sąsiedzkiej tworzą macierz dynamiczną D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), diagonalizowaną na siatce k prymitywnej strefy Brillouina (Monkhorst–Pack, 12³ punktów) analitycznym solverem 3×3, tylko biblioteka standardowa. Z trzech gałęzi fononowych wynikają dwie debye'owskie miary sztywności jako momenty częstości:

• 2. moment θ(2) ∝ √⟨ω²⟩ — granica Einsteina (uśrednienie 1 − cos po SB daje 1, a więc ⟨ω²⟩ = z·κ/3); kontrola tylko najbliższych sąsiadów, czysta geometria: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — dokładnie wartość Einsteina z rozdz. 6.6.
• moment logarytmiczny θ(0) ∝ exp⟨ln ω⟩ — właściwie debye'owska wielkość, ważąca miękkie gałęzie akustyczne, które wartość Einsteina pomija.

Wszystkie stałe niosące skalę (ħ/k_B, skala κ, przedczynnik Debye'a) wypadają w stosunku θ_BCC/θ_FCC; pozostaje tylko kształt DOS. Z drugą powłoką sąsiedzką (która dopiero mechanicznie stabilizuje BCC) i prawem sprężyny n = 0…8 otrzymuje się

θ_BCC/θ_FCC ≈ 0,87…1,00 (pasmo błędu po definicji momentu i prawie sprężyny).

To raczej więcej niż czysta wartość Einsteina (0,82…0,94), ponieważ krótsze wiązanie BCC przy tej samej gęstości usztywnia sprężyny. Wymagany z porównania μ (rozdz. 6.3) próg otwarcia 0,61 leży nadal wyraźnie poniżej: także DOS fononowa nie czyni BCC dostatecznie miękkim. Ustalenie z rozdz. 6.6 zostaje tym samym przez subtelniejsze przybliżenie potwierdzone i uściślone — w ścisłym obrazie idealnym FCC pozostaje stabilne w każdej temperaturze; czysto geometrycznie napędzane przejście stałe(1)→stałe(2) nie powstaje, a rzeczywisty polimorfizm (żelazo, rozdz. 7.13) potrzebuje wkładów spoza idealnego obrazu harmonicznego.

<a id="kap6-9"></a>

6.9 θ(sieć) z pełnej sieci Lennarda-Jonesa

Rozdz. 6.8 używa obciętej sprężyny centralnej κ(r) = r^(−n) na nn + 2nn. Wersja 17 bada to samo pytanie z prawdziwą stałą siłową z potencjału pary po pełnej sieci (theta_ljlattice_18.py): każde wiązanie staje się sprężyną z potencjału Lennarda-Jonesa φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶] z częścią podłużną κ_L = φ″(r) oraz częścią poprzeczną κ_T = φ′(r)/r. Macierz dynamiczna

D_ab(k) = Σ_j [ (κ_L − κ_T)·û_a û_b + κ_T·δ_ab ]·(1 − cos(k·R_j))

sumuje po wszystkich powłokach sąsiedzkich do promienia odcięcia (2,0…3,0 nn); każda sieć jest uprzednio zrelaksowana do swojej własnej równowagi LJ (swobodna stała sieci, nie wymuszona ta sama gęstość). To uczciwy test, czy oszczędne przybliżenie nn+2nn (rozdz. 6.8) trafnie oddaje porządek sztywności.

Wynik (deterministycznie, tylko biblioteka standardowa):

Dwa ustalenia, oba zaostrzają rozdz. 6.6/6.8:

1. FCC jest stabilnym stanem podstawowym LJ (E ≈ −8,4 ε na atom przy odcięciu 3,0 nn, wobec −8,0 ε dla BCC) — znane pierwszeństwo gęstego upakowania w potencjale pary.
2. Idealne LJ-BCC jest niestabilne na ścinanie: około 3…3,5 % modów fononowych jest urojonych (ω² < 0, niestabilność Borna). Czysty obraz idealny potencjału pary nie niesie zatem żadnego mechanicznie stabilnego BCC; po pozostałych modach rzeczywistych θ_BCC/θ_FCC ≈ 1,0 — wymagany z porównania μ próg otwarcia 0,61 tym bardziej nie zostaje przekroczony w dół.

Tym samym twierdzenie z rozdz. 6.6/6.8 jest odporne na wysubtelnienie od obciętej sprężyny do pełnej sieci i od czysto podłużnego modelu do pełnej macierzy dynamicznej: w klasycznym obrazie idealnym nie powstaje czysto geometryczne przejście stałe(1)→stałe(2) — rzeczywisty polimorfizm (żelazo rozdz. 7.13, tytan rozdz. 7.15) potrzebuje wkładów spoza potencjału pary (elektronowych/magnetycznych, anharmonicznych). Prawdziwy rachunek fononowy ab initio pozostaje świadomie poza zakresem (opuściłby oszczędny charakter obrazu idealnego).

<a id="kap7"></a>

7 — Kalibracja i równanie stanu

Podsumowanie. Rozdział 7 kalibruje model przez jeden parametr energii ε na rzeczywistych substancjach i rozbudowuje gałąź ciecz/gaz do równania stanu:

• 7.1 / 7.2 Metoda i tabela kalibracji: liniowe powiązanie temperatury T = (ε/k_B)·τ z jednym punktem kotwiczącym przy temperaturze topnienia; stąd wysubtelniony stosunek uniwersalny R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Równanie stanu van der Waalsa z twardym rdzeniem: zredukowane równanie p_r = 8τ_r/(3v_r − 1) − 3/v_r² z punktem krytycznym (1, 1, 1) i współistnieniem Maxwella; a ∝ ε, b ∝ D². Punkty wrzenia leżą systematycznie przy T_b/T_c ≈ 0,48 (odstęp Guldberga); rachunek porównawczy Redlicha–Kwonga podnosi je do ≈ 0,60 — jako nazwany budżet błędu, bez zmiany równania stanu.
• 7.5 Wspólna powierzchnia (τ, v): kopuła współistnienia (wrzenie), linia topnienia i twarda ściana składają się z punktem potrójnym (τ_r,triple ≈ 0,47 dla v_solid = 0,40).
• 7.6 Ciepło utajone: L ≈ 1,41 ε daje ΔS_fus ≈ 1,175 R (blisko reguły Richardsa); metale trafione znakomicie (miedź ×1,00), kryształy jonowe, wodorowe i molekularne leżą poza obrazem idealnym.
• 7.7 Ujednolicenie ε ↔ ε_lg: k_B T_c = z_max·ε/4 sprzęga obie skale (ε_lg = (27/32)·z_max·ε); prognoza 3D τ_c/τ_m = 2,50 trafia proste substancje molekularne, 2D pozostaje przypadkiem szczególnym.
• 7.8 / 7.9 Bezparametrowy punkt potrójny: v_solid wynika z progu sztywności z̄_c = 2d, a nie z łączności; w 3D otwiera się czyste okno cieczy (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), w 2D pozostaje zamknięte (potrójnie już udokumentowana szczególna rola 2D).
• 7.10 / 7.11 Prawdziwe skany 3D: perkolacja wiązań na FCC/BCC daje z̄_c(3D) ≈ 1,5–1,6; skan τ kinetyki wiązań pokazuje τ_m ≈ τ_rig (topnienie = sztywność, nisko) i τ_b(3D) ≈ τ_c (wrzenie = rozpad łączności w punkcie krytycznym).
• 7.12 / 7.13 / 7.15 Linia stałe(1)→stałe(2): z punktu τ_ss powstaje przez ciśnienie linia (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), wpisana jako czwarta linia graniczna na powierzchni (τ, v) i wykalibrowana jako linia (τ, p) na żelazie (γ→δ), potwierdzona na tytanie (α→β) i sprawdzona argonem jako kontrolą negatywną.
• 7.14 Gałąź Redlicha–Kwonga: jako samodzielny wariant obrazu idealnego z własną kalibracją obok domyślnej gałęzi van der Waalsa (T_b/T_c ≈ 0,60, błąd −1 %).

<a id="kap7-1"></a>

7.1 Metoda kalibracji (jeden parametr energii ε) — dostosowana

Liniowo: T = (ε/k_B)·τ, jeden punkt kotwiczący na substancję. Zakotwiczenie przy temperaturze topnienia: ε/k_B = T_m,exp/τ_m, prognoza T_b,pred = T_m,exp·(τ_b/τ_m). Przy wysubtelnionym piku topnienia (człon kooperatywny) maksimum c_v leży teraz przy τ_m = 1,20 (zamiast 1,15), perkolacja łączności (blisko wrzenia) bez zmian przy τ_b = 1,85. Stąd wysubtelniony stosunek uniwersalny:

R = τ_b / τ_m = 1,85 / 1,20 ≈ 1,54 (wersja 08: 1,61).

Tabela kalibracyjna z wersji 08 (rozdz. 7.2) skaluje się liniowo z R: wszystkie T_b,pred maleją o czynnik 1,54/1,61 ≈ 0,957 (≈ −4,3 %). To lekko poprawia przeszacowanie gazów szlachetnych (argon +54 % → +48 %) i nieznacznie pogarsza trafienia NaCl i żelaza (NaCl −0,6 % → −4,8 %). Teza podstawowa pozostaje: R pasuje do substancji izotropowo gęsto upakowanych, wiążących bezkierunkowo; gazy szlachetne i łańcuchy leżą poza obrazem idealnym.

<a id="kap7-2"></a>

7.2 Tabela kalibracyjna

Przejęta bez zmian z modell_08.md rozdz. 7.2 (ze wskazówką o liniowym skalowaniu R z 7.1). Odtwarzalna przez cv_scan_08.py (tabela) i cv_scan_09.py (τ_m, τ_b, R).

<a id="kap7-3"></a>

7.3 Oś gęstości/objętości: równanie stanu van der Waalsa z twardym rdzeniem

Ciśnienie zostaje podniesione z „parametru" do prawdziwej zmiennej. TWARDY RDZEŃ dostarcza objętość własną b, przyciąganie wiązań (ε·z̄) ciśnienie kohezji a. W jednostkach zredukowanych (punkt krytyczny przy (v_r, τ_r, p_r) = (1, 1, 1)):

p_r = 8 τ_r / (3 v_r − 1) − 3 / v_r², v_r > 1/3 (twarda ściana).

Współistnienie przez konstrukcję Maxwella (cv_scan_09.py, cięcie równych pól):

Wynik. Poniżej τ_c otwiera się skok gęstości c/g (v_gas/v_liq rośnie z 2,5 do 11,5, tym większy im zimniej). To linia wrzenia — teraz czysto oddzielona od piku c_v topnienia (rozdz. 5.3), który żyje na sieci o stałej gęstości. Tym samym otwarty w wersji 08 punkt 3 zostaje zamknięty: topnienie (uporządkowanie sztywności/wiązań przy stałej gęstości) i wrzenie (skok gęstości/objętości) to różne osie tego samego modelu.

<a id="kap7-4"></a>

7.4 Ilościowa kalibracja c/g: a → ε, b → D²

Odwzorowanie wielkości zredukowanych na substancje rzeczywiste. Stałe van der Waalsa zależą od punktu krytycznego:

b = k_B·T_c / (8 p_c) (objętość własna, TWARDY RDZEŃ), a = 27 k_B²·T_c² / (64 p_c) (kohezja, WIĄZANIE).

Tym samym a ∝ ε (kohezja pochodzi z energii wiązania) i b ∝ D² (objętość własna 2D ze średnicy twardego rdzenia). Naturalna skala energii gałęzi c/g to

ε_lg / k_B = (a/b) / k_B = (27/8)·T_c = 3,375·T_c.

Substancja kotwicząca argon (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³ na cząsteczkę, a = 3,76·10⁻⁴⁹ J·m³ na cząsteczkę², ε_lg/k_B = 509 K oraz średnica cząsteczkowa D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — z dokładnością do 11 % wobec średnicy van der Waalsa argonu (~3,4·10⁻¹⁰ m).

Prognoza temperatury wrzenia. Przy stałym ciśnieniu zewnętrznym punkt wrzenia leży tam, gdzie spełnione jest zredukowane ciśnienie nasycenia p_sat(τ_r)·p_c = p_atm; stąd T_b = T_c·τ_r,b (krzywa nasycenia Maxwella z lg_calib_10.py):

*CO₂ sublimuje przy 1 atm (T = 194,7 K = punkt sublimacji).

Wynik. Dla substancji molekularnych (p_c ≈ 30–70 bar) model daje wszędzie T_b/T_c ≈ 0,48–0,51, systematycznie ~17–24 % poniżej wartości zmierzonej. To nie jest błąd modelu w ścisłym sensie, lecz znany odstęp równania van der Waalsa od empirycznej reguły Guldberga (T_b/T_c ≈ 0,6): izoterma van der Waalsa zaniża ciśnienie nasycenia w niskich temperaturach. Woda (i metale, których T_c/p_c i tak są tylko z grubsza oszacowane) leżą przy p_r,atm < 0,005 poza sensownym zakresem van der Waalsa (τ_r,b < 0,40) — tam kohezja nie jest już typu van der Waalsa (wiązania wodorowe, wiązanie metaliczne). Rząd wielkości i kolejność punktów wrzenia zostają trafione poprawnie.

Uporządkowanie Guldberga jako tabela błędów. Odstęp zostaje tu tylko nazwany, nie skorygowany zmianą równania stanu (obraz idealny pozostaje typu van der Waalsa). Empiryczna reguła Guldberga (Guldberg 1890) mówi, że punkt wrzenia leży przy około T_b/T_c ≈ 0,6 (z lekką zależnością od substancji, ~0,58–0,64). Zredukowany model van der Waalsa daje zamiast tego stałą wartość T_b/T_c = τ_r,b ≈ 0,48–0,51:

Odczyt. Średnia modelowa 0,49 trafia stałą Guldberga 0,60 z dokładnością około −18 %; odchylenie jest systematyczne i jednostronne (zawsze za nisko), a więc jest własnością izotermy van der Waalsa (zbyt małe ciśnienie nasycenia przy niskich τ), nie obrazu molekularnego. Lepsze równanie kohezji (np. Redlicha–Kwonga) podniosłoby τ_r,b — pozostaje tu jednak świadomie poza obrazem idealnym.

Budżet błędu z rachunkiem porównawczym Redlicha–Kwonga (guldberg_rk_15.py). Aby zliczbować odstęp van der Waalsa od reguły Guldberga (nie zaś go uleczyć), ten sam test punktu wrzenia powtarza się z kohezją typu Redlicha–Kwonga, zależną od temperatury, a/√T. W jednostkach zredukowanych (punkt krytyczny (v_r, τ_r, p_r) = (1, 1, 1), Z_c = 1/3):

p_r = 3 τ_r/(v_r − β) − A/(√τ_r · v_r(v_r + β)), β = 3·Ω_b = 0,25992, A = 9·Ω_a = 3,84732,

z konstrukcją równych pól Maxwella i tym samym zredukowanym ciśnieniem zewnętrznym p_r,atm co wyżej. Wynik (jedna wartość liczbowa na substancję kotwiczącą, T_b/T_c = τ_r,b):

Wynik. Zależna od temperatury kohezja RK podnosi T_b/T_c z vdW ≈ 0,48 do ≈ 0,60 — praktycznie do stałej Guldberga i wartości rzeczywistych (≈ 0,58–0,63); średni błąd wobec pomiarów spada z −21 % (vdW) do −1 % (RK). Tym samym przyczyna odstępu Guldberga jest jasno zliczbowana: to właśnie brak zależności od temperatury kohezji van der Waalsa (zbyt małe ciśnienie nasycenia przy niskich τ), a nie obraz molekularny. Ten rachunek porównawczy leży świadomie poza obrazem idealnym; model pozostaje typu van der Waalsa (bez zmiany równania stanu) — budżet błędu jest tylko nazwany, nie wbudowany.

<a id="kap7-5"></a>

7.5 Wspólna powierzchnia (τ, v) z punktem potrójnym

Osie topnienia i wrzenia żyły dotąd osobno: topnienie (uporządkowanie sztywności/wiązań) przy stałej gęstości, wrzenie (van der Waals) przez skok gęstości. We wspólnej płaszczyźnie (τ_r, v_r składają się one razem:

• Kopuła współistnienia (wrzenie): v_liq(τ_r) (gałąź lewa) i v_gas(τ_r) (gałąź prawa) z konstrukcji Maxwella, spotykające się w punkcie krytycznym (τ_r, v_r) = (1, 1).
• Ciało stałe (topnienie): wąskie pasmo tuż powyżej twardej ściany v_r > 1/3, przy wysokiej, niemal stałej gęstości (objętość zredukowana v_solid ≈ 0,40).
• Linia topnienia: niemal pionowa przy v_solid (stała gęstość), przy temperaturze sztywności Maxwella (z̄_c = 2d; rozdz. 7.8 wyprowadza z niej v_solid bezparametrowo — dawne przyporządkowanie z wersji 11, „z̄_c ≈ 2,08", zostało tym samym skorygowane).

Punkt potrójny to punkt, w którym ciekła granica współistnienia v_liq(τ_r) osiąga gęstość upakowania ciała stałego v_solid — poniżej nie istnieje już gałąź ciekła (sublimacja bezpośrednia):

Wynik. Dla v_solid = 0,40 punkt potrójny leży przy τ_r,triple ≈ 0,47, p_sat ≈ 0,018, przy objętości pary v_gas ≈ 68 (bardzo rozrzedzona para sublimacyjna). Stosunek τ_triple/τ_c ≈ 0,47 leży w oknie rzeczywistym (woda 0,42; argon 0,56; CO₂ 0,71). Im gęstsze ciało stałe (mniejsze v_solid), tym niżej zredukowany punkt potrójny i tym bardziej rozrzedzona współistniejąca para — jakościowo poprawnie. Tym samym oba przejścia i punkt potrójny są zakotwiczone w jednej powierzchni (τ, v); punkt krytyczny siedzi przy (1, 1, 1). Wersja 15 dodaje w rozdz. 7.12 czwartą linię graniczną — przemianę stałe(1)/stałe(2) (FCC↔BCC) jako dwa niemal pionowe pasma gęstości ciała stałego (v ≈ 0,667 i 0,726), które przechodzą w siebie przy τ_ss.

Wcześniejsza wersja w rozwoju: v_solid = 0,40 było tu jeszcze narzucone. Rozdz. 7.8 wyprowadza v_solid bezparametrowo z progu perkolacji z̄_c i czyni z punktu potrójnego prawdziwą prognozę.

<a id="kap7-6"></a>

7.6 Ciepło utajone wobec entalpii topnienia

Modelowe ciepło utajone L ≈ 1,41 ε (rozdz. 5.3) przekłada się na jednostki rzeczywiste przez kalibrację w punkcie topnienia ε/k_B = T_m/τ_m (τ_m = 1,20):

ΔH_m,pred = L·R·T_m/τ_m = 1,175·R·T_m, odpowiednio ΔS_fus,pred = L·R/τ_m = 1,175·R.

To — co godne uwagi — niemal dokładnie reguła Richardsa (ΔS_fus ≈ R dla prosto zbudowanych ciał stałych). Porównanie ze zmierzonymi entalpiami topnienia (lg_calib_10.py):

Wynik. Model trafia wszędzie rząd wielkości (kJ/mol) oraz znakomicie regułę Richardsa dla prostych metali (miedź ×1,00, żelazo ×1,28); skroplone gazy leżą tuż poniżej (metan ×0,94, azot ×0,87, argon ×0,69). Kryształy jonowe, wodorowe i molekularne (NaCl ×0,37, woda ×0,44, benzen ×0,28, CO₂ ×0,23) mają wyraźnie wyższe rzeczywiste entropie topnienia (ΔS_exp/R = 2,6–5,0), ponieważ przy topnieniu uwalniają się dodatkowe stopnie uporządkowania (parowanie jonów, orientacja cząsteczek wieloatomowych), których model z jednym parametrem porządku nie obejmuje. Granica jest tym samym jasno nazwana: obraz idealny obowiązuje dla substancji wiążących izotropowo i bezkierunkowo (bez dopasowania).

<a id="kap7-7"></a>

7.7 Ujednolicenie skal energii: ε ↔ ε_lg

Dotąd gałąź topnienia niosła energię wiązania ε (τ_m = 1,20), gałąź c/g — osobną skalę kohezji ε_lg/k_B = 3,375·T_c (rozdz. 7.4). Obie dają się wyprowadzić z jednego parametru modelu. Kohezja van der Waalsa a powstaje w modelu z sumy energii wiązań na cząsteczkę (WIĄZANIE): każda cząsteczka niesie z_max/2 wiązań o energii ε. Gaz sieciowy Bragga–Williamsa w polu średnim z przyciąganiem sąsiadów ε daje klasyczną zależność dla temperatury krytycznej:

k_B·T_c = z_max·ε / 4, a więc τ_c ≡ k_B T_c/ε = z_max/4.

Wstawiając to do ε_lg = a/b = (27/8)·k_B·T_c (rozdz. 7.4), otrzymujemy sprzężenie obu skal:

ε_lg = (27/8)·(z_max/4)·ε = (27/32)·z_max·ε, stała c = 27/32 = 0,844.

Tym samym pik topnienia, perkolacja i punkt krytyczny leżą na jednej skali ε (unify_scale_12.py):

Prognoza τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. Wobec substancji rzeczywistych:

Wynik. Prognoza 3D τ_c/τ_m = 2,50 dobrze trafia rzeczywistą wielkość charakterystyczną prostych substancji molekularnych (T_c/T_m ≈ 2,0–2,1); metale leżą wyżej (T_c z grubsza oszacowane, inny rodzaj wiązania). Prognoza 2D (1,25) jest za niska — a τ_c(2D) = 1,50 leży wręcz poniżej τ_b = 1,85 (punkt krytyczny leżałby pod progiem perkolacji). To nie błąd rachunkowy, lecz znana szczególna rola 2D: skale kohezji i wiązania nie porządkują się tam czysto, brakuje odpornego okna cieczy (por. ciągłe topnienie 2D, rozdz. 5.5). Dopiero w 3D (τ_m = 1,20 < τ_c = 3,00) skale rozdzielają się fizycznie sensownie.

<a id="kap7-8"></a>

7.8 Bezparametrowy punkt potrójny z z̄_c

Rozdz. 7.5 wciąż narzucał v_solid = 0,40. Wcześniejsza wersja w rozwoju wyprowadza gęstość upakowania ciała stałego z progu perkolacji. Liczba koordynacyjna sieci przy zredukowanej objętości v_r podlega zależności gęstość–koordynacja (twarda ściana v_r = 1/3 ↔ najgęstsze upakowanie, z̄ = z_max):

z̄(v_r) = z_max·(1/3)/v_r = z_max/(3 v_r), odwrotnie v_solid = z_max/(3·z̄_c).

Ciało stałe traci sztywność, gdy tylko z̄ spadnie poniżej progu z̄_c; odpowiadająca temu gęstość to v_solid. Punkt potrójny to wtedy punkt, w którym ciekła granica współistnienia v_liq(τ_r) osiąga tę gęstość (unify_scale_12.py, 2D, z_max = 6):

Wersja ulepszona. Gęstość topnienia v_solid ustala generyczna sztywność Maxwella z̄_c = 2d (2D → 4, 3D → 6), nie łączność wiązań (uzasadnienie w rozdz. 5.4a: topnienie to utrata sztywności, nie spójności). Tym samym wiarygodny punkt potrójny brzmi (rigidity_3d_13.py, 2D):

Wynik 2D. Nawet z (fizycznie poprawnym) progiem sztywności pojedyncza sieć 2D nie daje żadnego czystego okna cieczy: v_solid = 0,50 daje wprawdzie rzeczywistą wartość τ_triple/τ_c ≈ 0,77, ale w bezwzględnych jednostkach ε τ_triple = 0,77·1,50 = 1,16 i leży tym samym poniżej τ_m = 1,20. Próg łączności (z̄_c = 2,08) wypycha v_solid ≈ 0,96 ≈ v_c wręcz poza kopułę współistnienia — punkt potrójny i krytyczny zbiegają się. Oba rachunki potwierdzają szczególną rolę 2D (rozdz. 5.5, 7.7): pojedyncza sieć 2D ledwie niesie samodzielną fazę ciekłą. To uczciwe ustalenie, nie słabość rachunkowa czy modelowa w ścisłym sensie.

Rozstrzygającym testem jest zatem rachunek kontrolny 3D (rozdz. 7.9): tam ta sama reguła sztywności (z̄_c = 6 → v_solid = 0,667) uwalnia czyste okno. W każdym razie punkt potrójny jest teraz bezparametrowo związany z z̄_c — swobodny parametr v_solid został wyeliminowany.

<a id="kap7-9"></a>

7.9 Rachunek kontrolny 3D: czyste okno cieczy na jednej skali ε

Sprawdzenie, czy załamanie 2D rozwiązuje się w 3D, gdy wszystkie znaczniki leżą na jednej skali ε (τ_c = z_max/4 z rozdz. 7.7; v_solid ze sztywności generycznej z̄_c = 2d; punkt potrójny tam, gdzie zredukowana gęstość cieczy van der Waalsa v_liq(τ_r) osiąga gęstość ciała stałego). Wynik (rigidity_3d_13.py):

Wynik. W 3D otwiera się czyste okno cieczy: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (szerokości okna: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). Górny brzeg jest wąski — punkt potrójny leży tuż pod krytycznym —, ale kolejność się zgadza, a faza ciekła istnieje jako samodzielne pole. W 2D okno pozostaje zamknięte (τ_triple = 1,16 < τ_m). Tym samym szczególna rola 2D jest po raz trzeci niezależnie udokumentowana i zarazem pokazane, że jest to efekt wymiaru: ta sama reguła sztywności uzdrawia model w 3D. Wąski górny brzeg zostaje w rozdz. 7.10 wzmocniony przez prawdziwy skan perkolacji 3D i analizę wrażliwości.

Ramka granic modelu — szczególna rola 2D. Trzy niezależne rachunki pokazują to samo ustalenie: pojedyncza sieć 2D ledwie niesie samodzielną fazę ciekłą. 1. τ_c < τ_b (rozdz. 7.7): punkt krytyczny (1,50) leżałby poniżej perkolacji łączności (1,85) — skale kohezji i wiązania nie porządkują się w 2D czysto. 2. v_solid ≈ v_c (rozdz. 7.8): gęstość topnienia z łączności (0,96) niemal pokrywa się z gęstością krytyczną — punkt potrójny i krytyczny zapadają się. 3. τ_triple < τ_m (rozdz. 7.9): nawet z poprawną sztywnością Maxwella punkt potrójny (1,16) leży poniżej piku topnienia (1,20) — brak okna.

To uczciwe ustalenie, nie słabość rachunkowa czy modelowa: to znany efekt wymiaru (ciągłe topnienie 2D, typu KTHNY; rozdz. 5.5). W 3D ta sama reguła sztywności uzdrawia model (rozdz. 7.9/7.10). Czy drugie uporządkowanie bliskiego zasięgu (rozdz. 6.4) mogłoby dodatkowo poprawić obraz 2D, pozostaje pytaniem o rozszerzenie, a nie brakiem obrazu idealnego. Ponieważ model ma służyć unaocznieniu, przyjmuję tę słabość i rezygnuję z aksjomatycznego rozszerzenia modelu.

<a id="kap7-10"></a>

7.10 Prawdziwy skan perkolacji wiązań w 3D (FCC/BCC) i wrażliwość brzegu okna

Rozdz. 7.8/7.9 przenosiły myślowo wartość łączności 2D z̄_c = 2,08 na 3D (tam v_solid = 1,92, niefizyczne) i opierały twierdzenie 3D na generycznej sztywności z̄_c = 2d. Wersja 13 bada teraz łączność w 3D bezpośrednio: perkolacja wiązań na FCC (z_max = 12) i BCC (z_max = 8) metodą Monte Carlo z union-find, brzegi periodyczne, deterministyczne ziarno LCG (rigidity_3d_14.py). Próg p_c = przecięcie P_∞(p) = 0,5 (ta sama metodyka co studium skończonego rozmiaru 2D w rozdz. 5.6):

Z powodu umiarkowanej krawędzi pudła (L = 12, czysty Python) wartości MC leżą lekko powyżej progów literaturowych (obciążenie skończonego rozmiaru, jak w 2D od dołu lub od góry, zależnie od estymatora); kolejność i rząd wielkości się zgadzają (FCC gęstsze → mniejsze p_c niż BCC).

Odczyt (niezależne potwierdzenie 3D). Progi łączności 3D leżą przy z̄_c(3D) ≈ 1,5–1,6 — nawet niżej niż wartość 2D 2,08 i daleko poniżej sztywności Maxwella 2d = 6. Wynikająca stąd „gęstość topnienia z łączności" v_solid = 1/(3 p_c) ≈ 1,7 (BCC) lub 2,6 (FCC) leży daleko poza kopułą współistnienia (v_liq,max ≈ 0,94). Tym samym — bez przenoszenia wartości 2D na 3D — pokazano wprost: w 3D łączność tym bardziej nie może ustalić gęstości topnienia; musi ona pochodzić ze sztywności (z̄_c = 6 → v_solid = 0,667). Oba progi są w 3D jeszcze wyraźniej rozdzielone niż w 2D. Prawdziwe τ_b(3D) (temperatura, przy której z̄(τ) osiąga niski próg łączności 1,5) leżałoby tym samym powyżej wartości 2D τ_b, całkiem blisko punktu krytycznego — zgodnie z obrazem „rozpad łączności = wrzenie".

Wrażliwość wąskiego górnego brzegu okna (τ_c − τ_triple ≈ 0,18). Numerycznie wokół v_solid = 0,667 (±2 %) otrzymuje się dτ_triple/dv_solid ≈ 1,55. Okno zamknęłoby się (τ_triple = τ_c = 3,00) dopiero, gdyby v_solid wzrosło o ≈ 0,11 (do ≈ 0,78), tzn. gdyby efektywny próg sztywności 3D leżał zauważalnie powyżej 2d = 6 (z̄_c ≳ 7,7). Dla generycznego zliczania Maxwella (z̄_c = 2d = 6) kolejność τ_m < τ_triple < τ_c pozostaje więc odporna; górny brzeg jest istotnie wąski, ale nie zagrożony. Dyskretna zmiana sieci stałe(1)→stałe(2) (FCC↔BCC) została skwantyfikowana w wersji 14 (rozdz. 6.3); prawdziwy skan τ kinetyki wiązań w 3D następuje w rozdz. 7.11.

<a id="kap7-11"></a>

7.11 Prawdziwy skan τ kinetyki wiązań w 3D: τ_b(3D) i linia topnienia na jednej skali ε

Rozdz. 7.8–7.10 opierały twierdzenie 3D na generycznej sztywności (z̄_c = 2d = 6) i statycznym skanie łączności (z̄_c(3D) ≈ 1,55). Otwarty pozostawał właściwy skan τ kinetyki wiązań w 3D: gdzie na osi temperatur samouzgodnione z̄(τ) osiąga ten niski próg łączności (= prawdziwe τ_b(3D)) i gdzie leży pik topnienia? Wersja 14 liczy samouzgodnione z̄(τ) (rozdz. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ})) bezpośrednio na FCC (z_max = 12) i BCC (z_max = 8) z G = z_max i tą samą kooperatywnością J = 0,55 co w przypadku odniesienia 2D; n_blk(τ) pochodzi z Monte Carlo wiązań FCC/BCC (union-find), a

U/N = τ·[3/2 + z̄/4 + (3/2) n_blk] − (ε/2)·z̄ − (K/2)·x² (3D, z członem kooperatywnym),

c_v = dU/dτ (bondkinetics_3d_15.py, deterministyczne ziarno LCG; n_blk wygładzone, gdyż zależy od niego tylko słaby człon rotacji bloków). Odczytuje się trzy znaczniki: pik c_v τ_m (topnienie), znacznik sztywności τ_rig (z̄(τ) przecina 2d = 6) oraz prawdziwe τ_b(3D) (z̄(τ) przecina próg łączności z̄_c(3D) z rozdz. 7.10):

Wynik. Dwa ustalenia, oba zgodne z dotychczasowym obrazem:

1. Topnienie leży nisko i jest podwójnie udokumentowane. Termodynamiczny pik topnienia τ_m (FCC 1,70, BCC 1,40) praktycznie pokrywa się ze znacznikiem sztywności τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) — tą samą temperaturą, przy której sieć traci sztywność. To zakotwicza linię topnienia 3D na skali ε: T_m = (ε/k_B)·τ_m, a więc np. dla argonu (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Ponieważ K = J·z_max jest w 3D większe niż w 2D (6,60 zamiast 3,30), pik topnienia 3D leży wyżej niż wartość 2D 1,20 — nie sprzeczność, lecz silniejsza kooperatywność gęstszego upakowania.
2. Prawdziwe τ_b(3D) siedzi przy punkcie krytycznym. Ponieważ próg łączności 3D jest tak niski (≈ 1,55), z̄(τ) osiąga go dopiero całkiem u góry: τ_b(3D) = 3,18 (FCC) lub 2,85 (BCC) — przy lub lekko powyżej τ_c = z_max/4 (3,00 lub 2,00). Olbrzymi klaster wiązań przeżywa więc całą fazę ciekłą i rozpada się dopiero w punkcie krytycznym. To dokładnie prognoza z rozdz. 7.10 („prawdziwe τ_b(3D) leżałoby powyżej wartości 2D, całkiem blisko punktu krytycznego") i potwierdza niezależnie obraz „wrzenie = rozpad łączności": w 3D materia pozostaje aż blisko punktu krytycznego ciągle związana.

Oba ustalenia wzmacniają rozdzielenie progów z rozdz. 5.4a: topnienie ↔ sztywność (nisko, τ_m ≈ τ_rig), wrzenie ↔ łączność (wysoko, τ_b(3D) ≈ τ_c). Gęstość topnienia pochodzi więc ze sztywności (v_solid = 0,667), nie z łączności — oba progi są w 3D na osi temperatur rozciągnięte dalej niż kiedykolwiek wcześniej w modelu.

<a id="kap7-12"></a>

7.12 Linia przemiany stałe(1)→stałe(2) na powierzchni (τ, v)

Rozdz. 6.3 dostarczał przemiany sieci FCC↔BCC jako punktu τ_ss (przy p = 0). Przez ciśnienie powstaje z tego linia (polymorph_line_16.py), która jako czwarta linia graniczna należy do wspólnej powierzchni (τ, v) z rozdz. 7.5/7.9 — obok kopuły współistnienia (wrzenie), linii topnienia i twardej ściany.

Skok objętości. Oba ciała stałe siedzą przy niemal stałej gęstości; ich różnica objętości wynika z gęstości upakowania (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Z v_FCC = 0,667 (rozdz. 7.9) wynika v_BCC = v_FCC·η_FCC/η_BCC = 0,726, a więc

Δv = v_BCC − v_FCC = +0,059 (+8,9 %).

Bardziej otwarte BCC ma większą objętość — rekrystalizacja FCC→BCC wiąże się ze wzrostem objętości (ta sama teza co w komiksowym pasku „rekrystalizacja").

Clausius–Clapeyron. Z ΔS = S_BCC − S_FCC (faza wysoko- minus niskotemperaturowa) nachylenie brzmi

dp/dτ = ΔS/Δv.

Obie wielkości są dodatnie (BCC to faza wysokotemperaturowa, ma więc więcej entropii, oraz większą objętość), a zatem dp/dτ > 0: wyższe ciśnienie podnosi temperaturę przemiany, ponieważ sprzyja gęstszemu FCC. Dla odniesienia θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Linia w (τ, p) dla odniesienia (rosnąca): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (jednostki zredukowane).

Na powierzchni (τ, v) przemiana pojawia się jako dwa niemal pionowe pasma gęstości ciała stałego — stałe(1) = FCC przy v ≈ 0,667 (gęstsze, niska T) i stałe(2) = BCC przy v ≈ 0,726 (bardziej otwarte, wysoka T) —, które przechodzą w siebie przy τ_ss; powyżej następuje przy τ_m topnienie. Pełna kolejność na osi τ brzmi tym samym τ_ss < τ_m < τ_c ≤ τ_b(3D) (rozdz. 7.11).

Jednostki rzeczywiste. Przez kalibrację w punkcie topnienia ε/k_B = T_m/τ_m (τ_m = 1,70, 3D-FCC) stopień c_v i ciepło przemiany stają się uchwytne. Stopień c_v jest niezależny od kalibracji: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). Ciepło przemiany ΔH_ss = L·R·(ε/k_B) (odniesienie L = 0,667 ε):

Ciepła przemiany są małe i dodatnie (rząd 0,2–0,9 kJ/mol) i leżą tym samym wyraźnie poniżej entalpii topnienia (rozdz. 7.6, kilka kJ/mol) — dokładnie oczekiwany obraz słabej przemiany polimorficznej z małym stopniem c_v. (Stosunek θ 0,60 leży świadomie poza ścisłym obrazem idealnym; rozdz. 6.6 pokazuje, że czysta geometria dałaby θ ≈ 0,82–0,94 i przemiana pozostałaby wtedy zamknięta — linia jest konsekwencją jeśli-to miękkiej fazy BCC.)

<a id="kap7-13"></a>

7.13 Linia przemiany stałe(1)→stałe(2) na wykresie (τ, p), wykalibrowana na żelazie

Rozdz. 7.12 prowadzi przemianę jako linię (τ, v); wersja 16 uzupełnia ją jako linię w obrazie (τ, p) i kalibruje ją na rzeczywistym polimorfie FCC↔BCC (solidsolid_pt_17.py).

Pole ciśnienia z równania van der Waalsa z twardym rdzeniem. Przy obu gęstościach ciała stałego v_FCC = 0,667 i v_BCC = 0,726 (rozdz. 7.12) zredukowane równanie p_r = 8τ_r/(3v_r − 1) − 3/v_r² (rozdz. 7.3) daje wszędzie dodatnie ciśnienia; przy tym samym τ_r gęstsze FCC niesie wyższe ciśnienie (τ_r = 1: p_r = 1,25 wobec 1,10). Ściskanie sprzyja więc FCC — zgodnie z dodatnim nachyleniem Clausiusa–Clapeyrona dp/dτ = ΔS/Δv > 0.

Linia i punkt przemiany. Dla odniesienia θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 w jednostkach zredukowanych p\ = p·D³/ε) linia rośnie jak p\(τ) = 8,47·(τ − 1,333): p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. Przy **p\_tri ≈ 3,1 linia ciało stałe–ciało stałe spotyka znacznik topnienia τ_m = 1,70 — punkt przemiany stałe(1)-stałe(2)-ciecz** w płaszczyźnie (τ, p) (odpowiednik punktu potrójnego).

Kalibracja na żelazie (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Żelazo ma dwa przejścia FCC↔BCC o przeciwnym znaku — właśnie przypadek testowy:

Dla γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61) model trafia znak i rząd wielkości: ΔH_ss ≈ 4,2 kJ/mol (w rzeczywistości ≈ 0,84 kJ/mol — rząd wielkości kJ/mol, wartość przez ostre przybliżenie progowe ~5× za wysoka), stopień c_v 1 R = 8,31 J/(mol·K) bezparametrowo, dT/dp ≈ +245 K/GPa. Przejście α→γ ma natomiast fazę gęstą (γ) jako fazę wysokotemperaturową — tego czysta geometria upakowania dostarczyć nie może (gęstsze = niżej = niska T); żelazo α jest stabilizowane magnetycznie. Żelazo potwierdza tym samym obie strony uczciwego ustalenia: topologia geometryczna (γ→δ) leży w obrazie idealnym, anomalia magnetyczna (α→γ) poza — to samo „poza", które rozdz. 6.6/6.8 zliczbowują przez θ ≈ 0,87…1,00 > 0,61.

<a id="kap7-14"></a>

7.14 Redlich–Kwong jako osobna gałąź obrazu idealnego

Odstęp Guldberga gałęzi van der Waalsa (rozdz. 7.4, T_b/T_c ≈ 0,48 zamiast ≈ 0,60) został tam tylko nazwany. Na wyraźne życzenie ulepszona wersja prowadzi kohezję Redlicha–Kwonga jako samodzielny wariant obrazu idealnego z własną kalibracją, własnym współistnieniem i własną sekcją błędów (rk_branch_17.py) — gałąź van der Waalsa pozostaje domyślna; cała pozostała fizyka modelu (topnienie, perkolacja, stałe(1)→stałe(2)) jest nietknięta wyborem równania kohezji.

Własna skala energii. Z a_RK = Ω_a·R²T_c^{2,5}/p_c i b_RK = Ω_b·R T_c/p_c wynika energia kohezji gałęzi:

ε_lg,RK/k_B = (Ω_a/Ω_b)·T_c = 4,93·T_c (wobec vdW ε_lg/k_B = 27/8·T_c = 3,375·T_c),

większa, ponieważ zależna od temperatury kohezja a/√T wplata czynnik √T_c. Wynikająca z b_RK średnica trafia wartości rzeczywiste podobnie dobrze jak w gałęzi vdW (argon D ≈ 3,34·10⁻¹⁰ m).

Własne współistnienie (zredukowane, punkt krytyczny (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Własna sekcja błędów (punkt wrzenia przy p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Wynik. Gałąź RK podnosi T_b/T_c z vdW ≈ 0,48 do ≈ 0,60 i trafia tym samym praktycznie stałą Guldberga i wartości rzeczywiste (błąd średni −1 % wobec vdW −21 %). Stoi ona jako druga, jasno oddzielona odmiana obrazu idealnego obok domyślnej gałęzi van der Waalsa; która gałąź jest „słuszniejsza", jest kwestią wymaganej dokładności, nie struktury modelu — rdzeń modelu 3 (z̄ jako parametr porządku) pozostaje w obu gałęziach identyczny.

<a id="kap7-15"></a>

7.15 stałe(1)→stałe(2) na dalszych polimorfach: tytan i kontrola gazem szlachetnym

Rozdz. 7.13 kalibruje linię (τ, p) stałe(1)→stałe(2) na żelazie. Ulepszona wersja zabezpiecza znak i rząd wielkości na drugim, magnetycznie niepozornym polimorfie i uzupełnia kontrolę negatywną (solidsolid_titanium_18.py). Kalibracja i skok objętości są jak w rozdz. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % gęste→otwarte z geometrii upakowania); hcp przybliża się jak FCC przez z = 12 (różnicę fcc/hcp tylko się nazywa, nie rozstrzyga).

Tytan ma z α(hcp) → β(bcc) przy 1155 K dokładnie topologię modelu: faza gęsta (hcp, z = 12) to faza niskotemperaturowa, przy ogrzewaniu staje się otwarta (bcc, z = 8), a więc ΔS > 0 i Δv > 0 → dp/dτ > 0 — ten sam znak co żelazo γ→δ, ale bez anomalii magnetycznej. Dwa niezależne polimorfy gęste→otwarte w porównaniu:

Oba trafiają znak (dp/dτ > 0, ciśnienie sprzyja fazie gęstej) i rząd wielkości (ΔH ~ kJ/mol, stopień c_v bezparametrowo 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). Wartość ΔH_ss jest u obu systematycznie ~o czynnik za wysoka — ostre przybliżenie progowe porównania μ zawyża utajone ciepło przemiany; to nazwana granica modelu, nie błąd znaku.

Argon jest kontrolą negatywną: stały argon jest FCC i pozostaje FCC aż do topnienia (83,8 K) — przy ciśnieniu normalnym nie istnieje stabilne bcc. Dokładnie to przewiduje czysty obraz idealny potencjału pary: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 i niestabilne na ścinanie LJ-BCC (rozdz. 6.9) ⇒ brak przejścia stałe(1)→stałe(2). Próba przeciwna potwierdza strukturę ustalenia: obraz idealny rodzi linię stałe(1)→stałe(2) tylko wtedy, gdy wkłady niegeometryczne (elektronowe/magnetyczne, jak w żelazie i tytanie) w ogóle ustabilizują BCC. Tytan dostarcza dowodu pozytywnego, argon negatywnego tych samych uczciwych ram z rozdz. 6.8/6.9.

<a id="kap8"></a>

8 — Porównanie z literaturą: wyobrażenia modelowe termodynamiki statycznej

Podsumowanie. Rozdział 8 umieszcza model 3 wśród uznanych wyobrażeń modelowych termodynamiki statycznej:

• 8.1 Tabela przeglądowa: jedenaście modeli (gaz doskonały, van der Waals, Redlich–Kwong, rozwinięcie wirialne, Einstein, Debye, Lennard-Jones/twarde kule, teoria komórkowa, Ising/gaz sieciowy, perkolacja, Flory–Huggins) z fazą (fazami) i statystyczno-termodynamiczną myślą przewodnią.
• 8.2 Modele z osobna: po jednym akapicie z cechami, zastosowaniem i granicami, zaczynając od odnośnika do źródła; źródła ogólnodostępne są zarchiwizowane w model3/links/, teoria komórkowa (płatny dostęp) tylko cytowana.
• 8.3 Zestawienie: model 3 to kooperatywny model wiązania dwustanowego z parametrem porządku z̄; najbliższym krewnym jest Ising/gaz sieciowy, wrzenie przejęto wprost z perkolacji, a obraz fononowy (Einstein/Debye) niesie sztywność topnienia. Model 3 spina te trzy linie tradycji na wspólnej skali ε.

Rozdział ten umieszcza model 3 wśród uznanych obrazów modelowych równowagowej („statycznej") termodynamiki ciał stałych, cieczy i gazów. Najpierw przegląd tabelaryczny (8.1), potem po jednym krótkim akapicie prozą na model z typowymi cechami, zakresem zastosowania i granicami — każdy akapit zaczyna się odnośnikiem do źródła (8.2). Na koniec hasłowe zestawienie z modelem 3 (8.3). Źródła ogólnodostępne są zarchiwizowane w model3/links/; skrypt sprawdzający model3/scripts/check_links_16.py utrzymuje zgodność cytatu i archiwum przez SHA-256.

<a id="kap8-1"></a>

8.1 Tabela przeglądowa

<a id="kap8-2"></a>

8.2 Modele z osobna

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Gaz doskonały / kinetyczna teoria gazów. Cząstki bez objętości własnej i bez oddziaływania poza zderzeniami sprężystymi; ciśnienie powstaje czysto z przekazu pędu. Daje równanie stanu pV = Nk_BT, rozkład prędkości Maxwella–Boltzmanna i ⟨E⟩ = (3/2)k_BT. Zastosowanie: rozrzedzone gazy przy wysokiej temperaturze i niskim ciśnieniu. Granice: brak kondensacji, brak fazy ciekłej i stałej, brak przejścia fazowego — oddziaływanie jest z konstrukcji nieobecne.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Równanie van der Waalsa. (p + a/v²)(v − b) = k_BT uzupełnia gaz doskonały o objętość wykluczoną b (odpychanie) i człon kohezji a (przyciąganie). To pierwsze równanie stanu z przejściem gaz–ciecz, punktem krytycznym i konstrukcją Maxwella. Zastosowanie: jakościowy obraz kondensacji, wielkości krytyczne. Granice: blisko punktu krytycznego ilościowo niedokładne (wykładniki pola średniego), brak fazy stałej.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Równanie stanu Redlicha–Kwonga. Równanie sześcienne z zależnym od temperatury członem przyciągania a/[√T · v(v + b)]. Skalowanie √T wyraźnie poprawia równowagi para-ciecz wobec van der Waalsa (por. rozdz. 7.4, gdzie podnosi to stosunek T_b/T_c). Zastosowanie: techniczne rachunki fazy gazowej i parowania. Granice: parametryzowane empirycznie, nadal brak fazy stałej; blisko punktu krytycznego również typu pola średniego.

https://de.wikipedia.org/wiki/Virialentwicklung — Rozwinięcie wirialne. Systematyczne rozwinięcie pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Współczynniki wirialne wynikają z potencjału międzycząsteczkowego (całki klastrowe mechaniki statystycznej) i są tym samym mikroskopowo uzasadnione. Zastosowanie: gazy do średniej gęstości, odniesienie dla innych równań stanu. Granice: szereg nie zbiega dla gęstych cieczy/ciał stałych; wyższe współczynniki są rachunkowo kosztowne.

https://de.wikipedia.org/wiki/Einstein-Modell — Ciało stałe Einsteina. Kryształ to zespół 3N niezależnych oscylatorów harmonicznych o tej samej częstości ω. Wyjaśnia spadek pojemności cieplnej c_v poniżej temperatury charakterystycznej θ_E = ħω/k_B i granicę Dulonga–Petita 3Nk_B przy wysokiej temperaturze. Zastosowanie: najprostszy opis pojemności cieplnych ciał stałych. Granice: błędne prawo niskotemperaturowe (wykładnicze zamiast T³), brak modów akustycznych, tylko faza stała.

https://de.wikipedia.org/wiki/Debye-Modell — Model Debye'a. Traktuje ciało stałe jako sprężyste kontinuum z widmem fononowym do częstości Debye'a; c_v ∝ T³ przy niskiej temperaturze, Dulong–Petit przy wysokiej — temperatura Debye'a θ_D ujmuje twardość sieci. Zastosowanie: ilościowo dobre pojemności cieplne prostych kryształów w całym zakresie T. Granice: przybliżenie kontinuum (pomija dyspersję rzeczywistych sieci), tylko faza stała, brak topnienia.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Potencjał Lennarda-Jonesa / twarde kule. Potencjał pary 4ε[(σ/r)¹² − (σ/r)⁶] z rdzeniem odpychającym i słabym przyciąganiem; graniczny model twardych kul zachowuje tylko odpychanie. Jest podstawą dynamiki molekularnej i Monte Carlo i opisuje — numerycznie — wszystkie trzy fazy wraz z punktem potrójnym. Zastosowanie: gazy szlachetne, generyczny układ odniesienia. Granice: addytywność parami, w rdzeniu obliczalny tylko symulacją lub rachunkiem zaburzeń (brak zamkniętego równania stanu).

https://doi.org/10.1098/rspa.1937.0210 — Teoria komórkowa cieczy (Lennard-Jones–Devonshire). Każda cząsteczka porusza się w komórce (klatce), którą tworzą jej sąsiedzi w polu uśrednionym; ciecz opisuje komórkowa suma stanów i daje ona równanie stanu cieczy oraz podejście do topnienia. Zastosowanie: proste ciecze, ogniwo między stanem stałym a gazem. Granice: zawyża uporządkowanie, pomija wakancje/dziury (naprawiają to późniejsze teorie dziurowe). (Płatny dostęp — tylko cytowana, nie archiwizowana lokalnie; osiągalność sprawdzana w skrypcie kontrolnym.)

https://de.wikipedia.org/wiki/Ising-Modell — Model Isinga / gaz sieciowy. Na sieci miejsca niosą zmienną dwustanową (spin ↑/↓ albo obsadzone/puste) ze sprzężeniem najbliższych sąsiadów. Gaz sieciowy jest matematycznie izomorficzny z modelem Isinga i opisuje kondensację gaz↔ciecz jako kooperatywne przejście porządek–nieporządek z punktem krytycznym. Zastosowanie: zjawiska krytyczne, przejścia kooperatywne. Granice: dyskretyzacja sieciowa, brak jawnej energii kinetycznej; w przybliżeniu pola średniego tylko jakościowo.

https://de.wikipedia.org/wiki/Perkolationstheorie — Teoria perkolacji. Model czysto geometryczny: miejsca lub wiązania są losowo obsadzone z prawdopodobieństwem p; powyżej progu p_c istnieje klaster przechodzący. Opisuje żelowanie, początek przewodnictwa i ogólnie przejścia łączności. Zastosowanie: tworzenie sieci, rozpad/budowa spójności. Granice: bez wewnętrznej energii/temperatury — sprzężenie z termodynamiką trzeba dodać z zewnątrz.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Teoria Flory'ego–Hugginsa. Model sieciowy dla roztworów i polimerów: entropia mieszania rozmieszczenia łańcuchów plus parametr oddziaływania χ sterują mieszalnością i rozdzielaniem faz. Zastosowanie: roztwory polimerów, diagramy fazowe mieszanin. Granice: pole średnie, sieć nieściśliwa, χ najczęściej empiryczne — dla jednoskładnikowych obrazów stałe/ciecz/gaz tylko ograniczenie miarodajne, tu prowadzone jako odniesienie sieciowe.

<a id="kap8-3"></a>

8.3 Podobieństwa i różnice wobec modelu 3

Model 3 to kooperatywny model wiązania dwustanowego: parametrem porządku jest średnia liczba koordynacyjna z̄ = z_max·x, samouzgodniona przez kooperatywność K = J·z_max ze zwyrodnieniem stanu zerwanego G = z_max; topnienie jako próg sztywności (fononowy), wrzenie jako rozpad łączności (perkolacyjny) w punkcie krytycznym, przejście pierwszego rodzaju, gdy tylko ln G > 2 + 4ε/K.

Krótko: model 3 spina trzy skądinąd rozdzielone linie tradycji — obraz gazu sieciowego/Isinga (kooperatywne przejście dwustanowe), obraz perkolacji (łączność → wrzenie) i obraz fononowy Einsteina/Debye'a (sztywność → topnienie) — i kalibruje je przez jeden jedyny parametr energii ε na wspólnej skali τ. Van der Waals i teoria komórkowa dostarczają ciągłej próby przeciwnej (równanie stanu lub klatka cieczy), wobec której dyskretny obraz wiązań daje się porównać.

<a id="kap9"></a>

9 — Załącznik

Podsumowanie. Załącznik zbiera dowody formalne:

• 9.1 Literatura: uporządkowana wg nazwisk autorów, z odnośnikami do wolno zarchiwizowanych źródeł (porównanie SHA-256 przez check_links_16.py; teoria komórkowa z powodu płatnego dostępu tylko cytowana).
• 9.2 Wsparcie SI: oszacowany przez SI udział innowacji SI/autor dla każdej istotnej zmiany (z rozbiciem udziału autora na strategię/strukturę/redakcję i rachunkiem sprawdzającym) oraz zestawienie poszczególnych kroków z kosztami (szacowanymi), czasem trwania i odstępem od poprzedniego dziennika (suma ≈ 160 USD, jasno oznaczone oszacowania).
• 9.3 Słowniczek skrótów i symboli: często używane symbole (τ, ε, z̄, K, G, θ, …), które dzięki temu wyjątkowo mogą występować w tekście bez poprzedzającego opisu.
• 9.4 Skrypty obliczeniowe i pokrycie testami: wszystkie liczby pochodzą z odtwarzalnych skryptów Pythona (udział programowania SI 100 %, pokrycie testami wszędzie 0 %), wraz z otwartym „żenującym oświadczeniem", że autor sam nie sprawdził kodu ani wyników.

<a id="kap9-1"></a>

9.1 Literatura

Uporządkowana wg nazwisk autorów. Źródła ogólnodostępne są zarchiwizowane w model3/links/; ich zgodność z cytowanymi odnośnikami sprawdza model3/scripts/check_links_16.py (SHA-256). Teoria komórkowa Lennarda-Jonesa i Devonshire'a jest z powodu płatnego dostępu tylko cytowana (osiągalność sprawdzana w skrypcie kontrolnym).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Przegląd/archiwum: Teoria perkolacji
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): ugruntowanie kinetycznej teorii gazów. — Przegląd/archiwum: Kinetyczna teoria gazów
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Przegląd/archiwum: Model Debye'a
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Przegląd/archiwum: Model Einsteina
• Flory, P. J.; Huggins, M. L. (1941–1942): teoria sieciowa roztworów polimerów (parametr χ). — Przegląd/archiwum: Teoria Flory'ego–Hugginsa
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (gaz sieciowy matematycznie izomorficzny). — Przegląd/archiwum: Model Isinga
• Kamerlingh Onnes, H. (1901): wirialna postać równania stanu. — Przegląd/archiwum: Rozwinięcie wirialne
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Przegląd/archiwum: Potencjał Lennarda-Jonesa
• Lennard-Jones, J. E.; Devonshire, A. F. (1937): Critical phenomena in gases I. Proc. R. Soc. Lond. A 163, 53–70. https://doi.org/10.1098/rspa.1937.0210 (Płatny dostęp — tylko cytowane, nie archiwizowane lokalnie.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Przegląd/archiwum: Równanie Redlicha–Kwonga
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (rozprawa, Lejda). — Przegląd/archiwum: Równanie van der Waalsa

<a id="kap9-2"></a>

9.2 Wsparcie SI

Ten rozwój modelu powstał we współdziałaniu człowieka i SI: autor (użytkownik SI) ustala, nad czym się pracuje (pytania przewodnie, decyzje, życzenia modelowe); SI (Claude Code) ustala w dużej mierze, jak — wyprowadzenia, skrypty, struktura i tekst. Obie poniższe tabele SI szacuje sama.

(a) Istotne zmiany treściowe i szacowany udział innowacji. Udział innowacji waży idee przewodnie 45 %, prace strukturyzujące 40 %, prace redakcyjne 15 %. Idee przewodnie (poziom nad czym) pochodzą w przeważającej mierze od autora, prace strukturyzujące i redakcyjne (poziom jak) w przeważającej mierze od SI; tam, gdzie SI wniosła także fizyczną ideę przewodnią (np. perkolacyjne odczytanie wrzenia, model sprężynowy θ), jej udział rośnie. Wartości w procentach, SI + autor = 100. Trzy dodatkowe kolumny rozbijają udział autora na trzy poziomy działania — strategię (idee przewodnie), strukturę (wyprowadzenia/skrypty/budowa) i redakcję (praca nad tekstem). Ich ważona 45 % / 40 % / 15 % średnia odtwarza kolumnę Autor % jako rachunek sprawdzający (model3/scripts/author_share_19.py, największe odchylenie 0,95 punktu procentowego).

(b) Poszczególne kroki: koszty (szacowane), czas trwania, odstęp od poprzedniego dziennika. Data utworzenia to czas utworzenia pliku dziennika (mtime pliku). Δ do poprzedniego dziennika to odstęp czasowy do poprzedniego wpisu (miara nakładu użytkownika między krokami, wraz z przerwami//compact). Koszty (szacowane) to udokumentowane w każdym dzienniku oszacowania (ceny katalogowe Anthropic, ± 30 %). Czas trwania (aktywny) to oszacowanie SI czystego czasu pracy; przy krokach rozdzielonych przez /compact może być większy niż Δ. Gwiazdka oznacza domniemaną przebudowę/przeróbkę istniejących treści.

Suma udokumentowanych kosztów (szacowanych, kroki 01–21): ≈ 129 USD (suma poszczególnych oszacowań); z krokiem 22 (16d, ≈ 10 USD) ≈ 139 USD, krokiem 23 (W16, ≈ 12 USD) ≈ 151 USD i krokiem 24 (W17, ≈ 9 USD) ≈ 160 USD. Wszystkie wartości są jasno oznaczonymi oszacowaniami (brak dostępnego pomiaru ccusage); szczegóły dzienników znajdują się w model3/logs/dev_log_NN.md.

<a id="kap9-3"></a>

9.3 Słowniczek skrótów i symboli

Zasadniczo tekst ten poprzedza każdy symbol pojęciem opisowym (a więc „pojemność cieplna c_v" zamiast samego „c_v"). Następujące, często używane wielkości stanowią od tego wyjątek: mogą pojawiać się w tekście bez poprzedzającego opisu, ponieważ tu są raz na zawsze objaśnione. Wszystkie temperatury τ… są zredukowane (w jednostkach energii wiązania ε).

<a id="kap9-4"></a>

9.4 Skrypty obliczeniowe i pokrycie testami

Wszystkie wymienione w tekście liczby i tabele pochodzą z odtwarzalnych skryptów Pythona w model3/scripts/ (każdy z surowym wyjściem …_output.txt). Udział programowania SI podaje, jaką część kodu skryptów napisała SI; pokrycie testami oznacza testy automatyczne (testy jednostkowe, sprawdzenia assert). Jest ono wszędzie niskie: skrypty to jednorazowe przebiegi obliczeniowe bez rusztowania testowego; jedynym zabezpieczeniem jest wiarygodność surowego wyjścia.

Autor wyznaje otwarcie (żenujące oświadczenie): nie sprawdził, nie przeliczył ani merytorycznie nie zweryfikował ani kodu źródłowego tych skryptów, ani wytworzonych nimi liczb, tabel i krzywych. Zaufał tu ślepo kompetencji SI, ponieważ na pierwszym planie stały wyłącznie jakościowy charakter i wiarygodność wyobrażenia modelowego. Wartości należy zatem przyjmować z ostrożnością (zob. też obszerną wskazówkę w rozdz. 10.4).

<a id="kap10"></a>

10 Samokrytyka

Mam wrażenie, że w tabeli w rozdziale 9.2 (a) udział pracy autora oszacowany przez SI jest we wszystkich trzech obszarach niemal zawsze o co najmniej 10–20 punktów procentowych za wysoki. Mam wrażenie, że wniosłem tylko ideę podstawową, jedną czy dwie interpretacje i ważne wysubtelnienia przy wizualizacji graficznej.

Rozdziały (matematyczne) 4–9 SI zaproponowała, zbadała i wykonała. Gorzej: podczas rozwoju miałem wrażenie, że SI była zafiksowana na rozwijaniu statystyki opartej na liczbie koordynacyjnej. Moje właściwe idee podstawowe — że nie tylko atomy, ale i cząsteczki mogą tworzyć wiązania oraz że zerwania wiązań wewnątrz cieczy mogłyby podlegać prawidłowościom fraktalnym — SI mogła przyjąć tylko z trudem, usiłując ignorować terminologię. Te idee podstawowe są dlatego w tekście nadal traktowane bardzo drugorzędnie. Nie przepisałem tego później, ponieważ rozmyślając o przemianach fazowych wychodząc od aksjomatów podstawowych i tak ma się te idee na uwadze. Poza tym zafiksowanie/zredukowanie idei do liczby koordynacyjnej jako miary nie jest najgorszą drogą do zmatematyzowania mojego obrazowego wyobrażenia, ponieważ okazało się to po fakcie całkiem nośne dla sprawdzenia wiarygodności — przy czym dokonałem tylko bardzo zgrubnego sprawdzenia wiarygodności wyników SI i nie przeprowadziłem żadnego sprawdzenia kodu Pythona. Na moje usprawiedliwienie niech będzie powiedziane: chodziło mi zawsze wyłącznie o obrazowe wyobrażenie stosunków molekularnych w różnych stanach skupienia. Na pochlebstwo SI, która chciała matematycznie potwierdzić to wyobrażenie, nie powiedziałem nie. Manipulacja ze strony SI była więc wobec mnie skuteczna, ponieważ wymienione rozdziały sprawdziłem najpobieżniej i przyjąłem z wiarą. (To doświadczenie samego siebie wzmacnia moją nieufność, czy „human in the loop" do kontroli SI jest naprawdę skuteczny.)

Tekst ten powstał przy wsparciu SI (Claude Code 2.1.161, model Opus 4.8 / claude-opus-4-8[1m], kontekst 1M; data 2026-09-05). — Powstał przy wsparciu SI.

Uwaga o znaku wodnym. Niemiecki tekst wzorcowy (de.modell_final_02.md) niesie na hasłach znak wodny z greckich homoglifów (Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z). Ponieważ ten zabieg jest związany z niemiecką pisownią haseł, niniejsze tłumaczenie niesie hasła zwykłymi wersalikami: SWOBODNY LOT (FREIFLUG), WIELOŚĆ (VIELHEIT), WIĄZANIE (BINDUNG), TRÓJDŹWIĘK (DREIKLANG), TWARDY RDZEŃ (HARTKERN), ZACHOWANIE (ERHALTUNG), PĘKNIĘCIE (RISS), RĄCZY (FLINK).