Statistische Thermodynamik für den idealen Reinstoff und seine Aggregatzustände

Gesamtzusammenfassung

Modell 3 ist ein qualitatives, statistisch-mechanisches Modell eines idealen Stoffs, das die drei Aggregatzustände fest, flüssig und gasförmig sowie ihre Übergänge aus einem einzigen Bindungsbild mit Molekülbindungen ableitet. Das Grundprinzip liegt in der Berücksichtigung von drei Energiespeichern für die innere Energie: die Schwingungsenergie in den Molekülbindungen, der Translationsenergie von den Molekülen und der kopperativen Bewegungsenergie in Rotation von Mizellen als idealisierte Form für innere Stoffflüsse zum Ausgleich innerer Dichteschwankungen. Eine wichtige Grundgröße in diesem Zusammenhang ist die mittlere Koordinationszahl z̄ = z_max·x, die sich über eine kooperative Zweizustands-Bindung (Ising-artig, Kooperativität K = J·z_max, Bruchzustands- Entartung G = z_max) selbstkonsistent einstellt. Aus ihr folgen drei Übergänge auf einer gemeinsamen, über einen einzigen Energieparameter ε kalibrierten τ-Skala: Schmelzen als Phononen-/Steifigkeitsschwelle (z̄ erreicht die Rigidität), Sieden als Konnektivitäts- (Perkolations-)Zerfall am kritischen Punkt und ein optionaler fest(1)→fest(2)-Übergang (FCC↔BCC) mit der parameterfreien c_v-Stufe von 1,0 k_B. Ein Übergang wird erster Ordnung, sobald ln G > 2 + 4ε/K. Dabei ist die Bruchzustands-Entartung G die Zahl der aufbrechbaren Nachbarkontakte (G = z_max, die maximale Koordinationszahl), die Bindungsenergie ε die je Nachbarkontakt gewonnene Energie und die Kooperativität K = J·z_max das Produkt aus Ising-artiger Kopplungsstärke J und z_max; der Übergang wird also umso eher erster Ordnung, je größer die Kontaktvielfalt G und je stärker die Kooperativität K gegenüber der Bindungsenergie ε ausfällt.

Der Artikel führt von den sieben Axiomen (Kapitel 1) über die Widerspruchsanalyse (Kapitel 2) und eine bildhafte Beschreibung mit statischen Comic-Strips (Kapitel 3) zum statistisch-mechanischen Gerüst — kooperative Bindung, Zustandssummen und Wärmekapazitäten (Kapitel 4–5) —, zur mathematischen 3D-Plausibilität (Kapitel 6) und schließlich zur Kalibrierung und Zustandsgleichung (Kapitel 7). Die Gedanken werden mit den vollen 3D-c_v-Scan des Schmelzens erster Ordnung (Δz̄, latente Wärme ΔU, Peak-Verschiebung), ein Einstein-Federmodell für die Phononensteifigkeit θ(Gitter) und die fest(1)→fest(2)-Umwandlungslinie via Clausius–Clapeyron ergänzt. Weiter wird das θ(Gitter) zu einer Debye-artigen Phononen-DOS verfeinert, Dann wirddie fest(1)→fest(2)-Umwandlung zusätzlich als Linie im (τ, p)-Bild (an Eisen kalibriert) eingeführt und es wird ein getrennten Redlich–Kwong-Zweig neben dem van-der-Waals-Default ergänzt. Dabei wird as θ(Gitter) an einem Lennard-Jones-Vollgitter (FCC stabil, ideales BCC scherinstabil) und die fest(1)→fest(2)-Topologie an Titan (zweiter Bestätiger) und Argon (Negativkontrolle) geprüft. Kapitel 8 ordnet das Modell in elf etablierte Modellvorstellungen der statischen Thermodynamik ein. Der Anhang dokumentiert die Literaturquellen (nach Autoren sortiert) und die KI-Unterstützung (geschätzter Innovationsanteil, Kosten (geschätzt) und Zeitdauer je Entwicklungsschritt); die Hilfen bündeln Lesart, verbliebene kritische Punkte und Ausblick.

Kernaussage. Ein einziges kooperatives Bindungsmaß z̄ verklammert drei sonst getrennte Traditionslinien — das Gittergas-/Ising-Bild (kooperativer Zweizustands-Übergang), das Perkolationsbild (Konnektivität → Sieden) und das Phononenbild Einstein/Debye (Steifigkeit → Schmelzen) — zu einem konsistenten, qualitativ vollständigen Phasenbild, bewusst innerhalb eines van-der-Waals-artigen Idealbildes, wobei sich das Theoriegebäude natürlich aus einem leicht vorstellbaren, dynamischen Bild ableiten lässt. Nicht die Genauigkeit, sondern die Bildhaftigkeit ist die eigentliche Stärke des Modells und kennzeichnet den wissenschaftlichen Modellwert.

1 — Die Axiome mit Schlagwort-Namen (konsolidiert)

Zusammenfassung. Sieben mit Schlagwort-Namen versehene Axiome legen den idealen Stoff fest: erhaltene Teilchenzahl, rein klassische Bewegung, eine kooperative Zweizustands-Bindung und die mittlere Koordinationszahl z̄ als Ordnungsparameter. Sie tragen alle folgenden Herleitungen; ein Homoglyphen-Wasserzeichen markiert die Schlagwörter (Schlüssel am Textende).

Präzisierte Axiome:

1. FRΕΙFLUG (Zusatzaxiom). Zwischen zwei Molekülen, die ungebunden aneinander vorbeifliegen, wirken weder Reibung noch elektrostatische Kräfte. Impuls und Energie ändern sich nur durch die harte ΗΑRΤΚΕRΝ-Abstoßung beim Stoß und durch das Knüpfen/Brechen von ΒΙΝDUΝG. Das Axiom ist ausschließlich aus Sicht der klassischen Mechanik und der Energiebilanz von Belang: Es garantiert, dass es keine dissipativen Kanäle auf molekularer Ebene gibt — Bewegungsenergie geht nie „verloren", sondern nur in die drei DRΕΙΚLΑΝG-Speicher oder in Bindungsenergie über. (Makroskopische Reibung/Viskosität entsteht emergent aus dem Schiebepuzzle und dem RΙSS, nicht als mikroskopische Grundkraft.)
2. VΙΕLΗΕΙΤ. Ein idealer Stoff besteht aus sehr vielen gleichartigen Molekülen; die makroskopischen Eigenschaften werden im thermodynamischen Limes (N → ∞) statistisch beschrieben.
3. ΒΙΝDUΝG. Je zwei geometrisch benachbarte Moleküle können eine Bindung tragen. Eine intakte Bindung senkt die Energie um den endlichen Betrag ε; ihr Bruch kostet genau ε. Die Bindung ist ein Zweizustandssystem {intakt, gebrochen}.
4. DRΕΙΚLΑΝG. Die innere Energie verteilt sich klassisch auf drei Speicher: (a) Translation der Moleküle, (b) Schwingung in intakten Bindungen (harmonische Feder je Bindung), (c) Rotation starrer Molekülblöcke (kollektiver Drehfreiheitsgrad zusammenhängender Bereiche).
5. ΗΑRΤΚΕRΝ. Unterschreitet der Molekülabstand einen Mindestwert σ, wirkt eine harte Abstoßung (Eigenvolumen). Zusammen mit ΒΙΝDUΝG ergibt sich ein effektives Paarpotential (harte Schale + endliche Bindungsmulde).
6. ΕRΗΑLΤUΝG. Gesamtimpuls und Gesamtenergie sind erhalten. Ein einzelnes Molekül ist energetisch neutral — es puffert keine Energie in sich; Energie lebt in den Bindungen (ΒΙΝDUΝG) und den kinetischen Freiheitsgraden (DRΕΙΚLΑΝG).
7. RΙSS. Bindungen zerfallen/entstehen nach dem Riss-/Heilprinzip: mehr fehlende Nachbarbindungen → höhere Bruchwahrscheinlichkeit; mehr gebundene Nachbarn → höhere Heilwahrscheinlichkeit. Die Kinetik ist kooperativ.
8. FLΙΝΚ. Bruch und Heilung laufen schnell und aktivierungsarm ab; das Bindungsteilsystem steht zu jeder Temperatur im lokalen Gleichgewicht und darf gegenüber der Molekülbewegung adiabatisch separiert werden (Grundlage der faktorisierten Zustandssumme).

Modellierungshaltung. Die geometrische Koordination bleibt idealisiert erhalten; es ändert sich nur die Koordinationszahl z durch Brüche. Bewegung entsteht wie beim Schiebepuzzle: Hohlräume (aus ΗΑRΤΚΕRΝ-Eigenvolumen und RΙSS-Brüchen) wandern, Blöcke rücken nach. Das bewusst „falsche" Bild rotierender, wandernder Kriställchen bleibt als anschauliche Grundlage erhalten.

2 — Widerspruchs- und Lückenanalyse

Zusammenfassung. Die Analyse prüft das Modell systematisch auf innere Widersprüche und Lücken und schreibt den Bearbeitungsstand fort; erkannte Schwachstellen werden benannt und in den Folgekapiteln aufgelöst oder als ehrliche Modellgrenze markiert.

Liste von Kritikpunkten

1. Doppelzählung der Translation (Axiom 3a vs. 3c). Der Rohtext nennt Translation zweimal (Molekül-Translation und Block-Translation). Zwei unabhängige Translationsspeicher gäbe es nicht. Schließung: Die Block-Translation ist nur eine kollektive Mode der bereits gezählten Molekül-Translationen; genuin neu ist allein der Rotations-Freiheitsgrad des Blocks. In der Freiheitsgradbilanz (Kap. 4) wird Translation genau einmal gezählt.
2. „Energetisch neutral" (Axiom 5) vs. „Bindungen nehmen Energie auf" (Axiom 2/3b). Scheinbarer Widerspruch. Schließung: Die Neutralität betrifft das Molekül als Körper; Energie sitzt zwischen den Molekülen (in ΒΙΝDUΝG) und in kollektiven kinetischen Moden, nie in einem intramolekularen Puffer. Molekül ≠ Bindung — damit ist der Widerspruch aufgelöst.
3. Fehlende Anziehung (Axiom 4 nennt nur Abstoßung). Ohne Anziehung gäbe es keine kondensierte Phase. Schließung: Die Anziehung liefert ΒΙΝDUΝG (Mulde der Tiefe ε), die Abstoßung ΗΑRΤΚΕRΝ (harte Schale bei σ). Zusammen ein effektives Paarpotential.
4. Rückkopplung ohne Gegenterm (Axiom 6). „Mehr fehlende Nachbarn → mehr Bruch" allein ist eine positive Rückkopplung (Rissfortschritt bis zum Totalzerfall). Schließung: Der Heilterm (mehr Nachbarn → mehr Heilung) plus der Boltzmann-Faktor der endlichen Bindungs- energie ε (ΒΙΝDUΝG) balancieren den Riss und erzeugen einen definierten Übergang statt eines Runaways.
5. Unausgesprochene Zeitskalentrennung (Axiom 7). „Schnell und aktivierungsarm" wird erst nutzbar, wenn man daraus die adiabatische Separation Bindungs- vs. Translationsdynamik macht. Schließung: explizit als Voraussetzung der Faktorisierung Z = Z_trans·Z_vib·Z_rot·Z_bond festgehalten.
6. Fehlende Geometrie und fehlender Bewegungsmechanismus. Die Axiome nennen weder das Koordinationsgitter noch, was Fließen ermöglicht. Schließung (als kenntlich gemachte Zusatzannahme): feste Referenz-Koordinationsgeometrie mit variabler Bindungsbelegung; als kanonischer 2D-Fall die Dreiecks-Koordination (z_max = 6). Beweglichkeit über den Schiebepuzzle-Hohlraummechanismus. Diese beiden Punkte sind die einzigen echten Zusätze über die sieben Axiome hinaus und in Kapitel 6 als strategische Frage markiert.
7. Latente Dissipationslücke. Fassung 07 ließ offen, ob ungebundene Moleküle beim Vorbeiflug Energie verlieren. Schließung durch FRΕΙFLUG: Es gibt keine mikroskopische Reibung/Elektrostatik; alle Energie- und Impulsflüsse laufen ausschließlich über ΗΑRΤΚΕRΝ-Stöße und ΒΙΝDUΝG-Kinetik. Damit ist die Energiebilanz (Kap. 4.3) geschlossen und die Maxwell–Boltzmann-Statistik im Gasgrenzfall exakt (freies Idealgas zwischen elastischen Stößen).
8. Energetische Inkonsistenz des Ordnungsparameters. Fassung 08 nutzte den kooperativen Term −(K/2)x² zwar in der freien Energie g(x) (zur Erzeugung von x(τ)), aber nicht in der inneren Energie U/N. Das ist thermodynamisch inkonsistent: derselbe Molekularfeld-Term muss in beiden erscheinen. Schließung (Kap. 4.3): −(K/2)x² wird in U/N nachgezogen; die latente Wärme erhöht sich dadurch messbar (Kap. 5.3), ohne Fest- und Gasplateaus zu verändern.
9. l/g-Kalibrierung nur qualitativ. Fassung 09 hatte die Dichteachse (van-der-Waals-EOS) nur reduziert. Schließung (Kap. 7.4): über den kritischen Punkt (T_c, p_c) wird die absolute Siedelinie berechnet und T_b(p_atm) gegen Messwerte gestellt; die systematische Abweichung (Guldberg-Regel) wird benannt.
10. Getrennte Achsen ohne Tripelpunkt. Schmelz- und Siedeachse waren getrennte Rechnungen. Schließung (Kap. 7.5): beide werden in einer (τ, v)-Ebene zusammengeführt, der Tripelpunkt konstruiert und der kritische Punkt markiert.

3 — Bildhafte Beschreibungen

Zusammenfassung. Kapitel 3 erklärt fest, flüssig und gasförmig sowie ihre Übergänge anschaulich (mit statischen Comic-Strips): die drei Einzelphasen (3.1), die Phasenübergänge Umkristallisieren/Schmelzen/Sieden/Sublimieren samt Tripel- und kritischem Punkt (3.2), die Schallanomalie (3.3), Brownsche Bewegung und den geschärften Mizellen-Begriff (3.4–3.5) sowie das dreidimensionale Flüssigkeitsfenster (3.6). Basierend auf den Beschreibungen wurden Schaubilder (Comic-Strips) erstellt, die helfen sollen, die dynamischen molekularen Verhältnisse in den verschiedenen Phasen vor dem inneren geistigen Auge verbildlichen zu können. Diese Vorstellbarkeit ist der eigentliche Kern des Modells. Qualitativ und wissenschaftlich wertvoll wird es dann, wenn sich über die statistische Mechanik plausibel zeigen lässt, dass das Modell qualitativ empirische Tendenzen/Beobachtungen auch ohne große Mathematik vorhersagbar macht.

3.1 Die drei Einzelphasen

Feststoff. Fast alle Nachbarn sind miteinander über Molekülbindungen gebunden (ΒΙΝDUΝG), der Mittelwert der Koordinationszahlen für die Moleküle z liegt nahe am geometrischen Maximum. Das Molekülnetz ist ein einziger starrer Block ohne Hohlräume (ΗΑRΤΚΕRΝ dicht gepackt) — ein inneres Schiebepuzzle wie bei der Flüssigkeit ist in der kompakten Version nicht möglich. Einzige Bewegung: Zittern um die Ruhelagen (Bindungsschwingung, DRΕΙΚLΑΝG-b), die durch die Molekülbindungen miot den nachbarn bestimmt werden. Impuls/Energie pendeln elastisch im Netz (ΕRΗΑLΤUΝG, FRΕΙFLUG: verlustfrei). Im Grenzfall (sehr niedrige temperature, ideale Packung)liegt der Idealzustand eines idealen Kristalls vor.

Flüssigkeit. Der kooperative RΙSS hat das Netz durchsetzt: Bindungen fehlen fleckenweise, der Mittelwert der Koordinationszahlen für die Moleküle z ist gesunken, der starre Block ist innerlich in viele kleinere zerfallen. Aus dem Monolithen wurde ein Geröllhaufen. Diese Blöcke translatieren und rotieren (DRΕΙΚLΑΝG-c). Zwischen ihnen sind Hohlräume aufgegangen, durch die sie aneinander vorbeirücken — das Schiebepuzzle der Browbnschen Bewegung ist aktiv. Die lokale Koordinationsgeometrie bleibt idealisiert erhalten; verändert ist nur, wie viele der Bindungen wirklich real existieren und ausgebildet werden.

Gas. Der RΙSS hat gesiegt: Bindungen fast überall gebrochen (z → 0). Im Extremfall hat kein Moelkül noch eine Molekülbindung ´zu einem anderen Molekül. Die Blöcke sind in einzelne Moleküle zerfallen. Sie fliegen frei, nur durch ΗΑRΤΚΕRΝ-Stöße gekoppelt (dazwischen FRΕΙFLUG); die Geschwindigkeiten folgen Maxwell–Boltzmann (ΕRΗΑLΤUΝG). Bindungsschwingung und Blockrotation sind erloschen — reine Translation (DRΕΙΚLΑΝG-a). Im Grenzfall (hohe Temepratur, geringer Druck) liegt der Idealzustand eines idealen Gases vor.

3.2 Phasenübergänge

Umkristallisieren (fest(1)/fest(2)). Noch unterhalb des Schmelzens kann der Feststoff sein Gitter wechseln, ohne die Steifigkeit zu verlieren: Zwei ideale Packungen mit verschiedener Koordinationszahl (FCC z_max = 12 ↔ BCC z_max = 8) haben verschiedene freie Enthalpien μ(τ, p); bei μ_FCC = μ_BCC springt die Koordinationsgeometrie sprunghaft um (nicht nur die Belegung z̄ wie beim Schmelzen). Sichtbar wird eine kleine c_v-Stufe (Δc_v = z_max/4-Änderung = 1,0 k_B) plus eine schwache latente Umwandlungswärme. KI-Behauptung: Es ist nur in 3D denkbar (diskrete, quantisierte Konnektivität; Kap. 6.3). Im Zweidimensionalen Fall lässt sich der echte Umkristallisieren, also der Übergang zwischen zwei festen Phasen, leichter vorstellen. In der Fläche gibt es nur drei regelmäßige Pakettierungen ohne Löcher. Man kann mit gleichseitigen Dreiecken, mit Quadraten oder mit regelmäßigen Sechsecken Flächen ohne Löcher eine Fläche vollständig bedecken (parkettieren). Wenn in diesen Gittern die Kanten die Molekülbindungen und die Ecken die Moleküle repräsentieren, dann ist ein Phasenübergang von fest(triangular)->fest(quadratisch) denkbar. Dabei würde das Molekül im Dreiecksparkett noch sechs Bindungenzu benachbarten Molek+ülen aufweisen würde, während das Molekül im Quadratgitter nur noch vier Bindungen hätte.

Schmelzen (s/l). Bei steigender Temperatur reißt das Netz an genügend Stellen, dass der starre Gesamtblock seine Steifigkeit verliert (kooperativer RΙSS unterschreitet die Maxwell-Rigiditäts­schwelle z̄ = 2d — das Netz wird floppy). Es bilden sich thermodynamische Mizellen, also Bereiche, die mit größeren benachbarten Bereichen nicht mehr über Bindungen miteinander verbunden sind. Der Monloith wird als molekularer Sicht zum Geröllhaufen und ist bedingt flie´ßfähig. Die ersten Hohlräume erlauben die Verschiebung von Molekülblöcken ähnlich wie im Schiebepuzzle. Die fehlenden Bindungen zur Nachbarschaft erlauben die innere Rotation in den Mizellen. Die Mizellenrotation (DRΕΙΚLΑΝG-c) setzt ein. Das Brechen vieler Bindungen (ΒΙΝDUΝG) verschlingt latente Wärme.

Sieden (l/g). Weitere Erwärmung lässt die Hohlräume selbst perkolieren: Sie verbinden sich zu durchgehendem Leerraum, in dem die Blöcke sich auflösen und einzelne Moleküle davonfliegen. Sieden ist der Übergang von „Hohlräume in Materie" zu „Materie im Hohlraum".

Sublimieren (s/g). Bei niedrigem Druck fehlt der flüssige Zwischenbereich: Das Netz gibt direkt einzelne Moleküle an den Leerraum ab, ohne dass es zur Blockbildung wie bei den Mizellen kommt. Man kann hier überlegen, welchen Einfluss geometrische Aspekte (lange Moleküle, minteinader verpflochtende Moleküle), elektrostatische Aspekte (Ausrichtung durch magnetische Momente oder Dipolmomente), verschiedene Molekülbindungstypen (Van-der Wall, Wechselwirkungen zwischen mesomerie-Systeme, Wasserstoffbrückenbindungen) auf die Entstehung oder Nichtenststehung von thermodynamischen Mizellen haben können.

Tripelpunkt. Genau dort fallen zwei Schwellen zusammen: die Maxwell-Rigidität (fest ↔ flüssig, Verlust der Steifigkeit) und die Hohlraum-Perkolation (flüssig ↔ gas) treten bei demselben (τ, p) ein. Alle drei Bilder (starrer Block[Monolith], schiebende Blöcke[Gerööhaufen], freie Moleküle[Staubwolke]) koexistieren.

Kritischer Punkt. Mit steigendem Druck wird der Sprung flüssig→gas immer kleiner; der Hohlraumanteil ändert sich zunehmend stetig. Am kritischen Punkt verschwindet der Unterschied ganz — „dichte Blöcke mit wenig Leerraum" und „viel Leerraum mit wenigen Molekülen" gehen kontinuierlich ineinander über. ("Aus Sicht der Flüssigkeit ist das Glas halb leer, aus der Sicht des Gases halb voll.").

Umkristallisieren im 2D-Bild. Der Comic zeigt das Umkristallisieren jetzt als 2D-Fall untereinander: unten das dicht triangulierte Gitter (Koordination z = 6, das energieärmste 2D-Ideal), oben das lockere Quadratgitter (z = 4). Der Übergang trianguliert → quadratisch vergrößert das Volumen (Arbeit gegen den Druck) und senkt die Bindungszahl. man sieht, dass die Koordinationszahl wechselt, aber Steifigkeit eine Kristallgitters bestehen bleibt.

Qualitative Vorhersage von Energie-Übergängen und Wärmekapazitäten. Aus der Konnektivität z allein lassen sich beide Kennzahlen abschätzen: (i) die Bindungsenergie je Teilchen ist U/N ≈ −(ε/2)·z (jede Bindung zählt hälftig), also energieärmer je höher z; (ii) der c_v-Sprung beim Gitterwechsel folgt der Modellregel Δc_v = z_max/4-Änderung. Für den 2D-Fall trianguliert (z = 6) ↔ quadratisch (z = 4) ergibt das eine Energiedifferenz ΔU ≈ ε je Teilchen und eine c_v-Stufe Δc_v = (6 − 4)/4 = 0,5 k_B — halb so groß wie der 3D-Wert (1,0 k_B). Merksatz: dichter = energieärmer, aber teurer im Übergang (größere latente Wärme, größere c_v-Stufe), weil beim Aufweiten mehr Bindungen und mehr eingefrorene Freiheitsgrade freigesetzt werden.

Bestimmt die Gitterstruktur die Metastabilität? Teils. Die Gitterstruktur legt die Energie-Reihenfolge fest — welche Form Grundzustand (trianguliert, tiefstes μ) und welche angeregt (quadratisch) ist. Ob die lockere Form metastabil bei tiefer Temperatur überlebt, entscheidet sie aber nicht allein: Im streng idealen Modell (ohne Barriere) nimmt der Stoff bei jedem (τ, p) sofort die Phase mit dem niedrigsten μ an — Metastabilität gibt es dort nicht. Sie entsteht erst durch eine endliche Aktivierungsenergie (Keimbildungsbarriere), die den Umbau kinetisch hemmt. Die Gittergeometrie liefert also die thermodynamische Rangordnung, die Barriere die kinetische Persistenz — dieselbe Zutat, die auch die übersättigte Phase erklärt (s. u.; quantitativ Kap. 6.7).

Übersättigte Phase — Aktivierungsenergie. Der Zusatzabbildung (3.2, Reihe 2) zeigt eine übersättigte Gasphase über einer Flüssigkeit: neben Einzelmolekülen kleine Gruppen aus 2–3 zusammengelagerten Teilchen (graue Balken), die als Ganzes fliegen. Das Modell erklärt die Übersättigung durch eine merkliche Aktivierungsenergie für Molekülbindungen: Kleine Cluster (2–3) haben ein ungünstiges Rand-zu-Fläche-Verhältnis, ihre Randbindungen sind unvollständig — erst ein kritischer Keim wächst spontan weiter. Solange kein überkritischer Keim entsteht, bleibt das Gas übersättigt, obwohl μ_gas > μ_flüssig bereits Kondensation erlaubte. Genau diese Barriere fehlt im idealen Gleichgewichtsbild und wird hier als eigenständiger, „außerhalb des Idealbildes" liegender Beitrag benannt (quantitative Abschätzung in Kap. 6.7).

3.3 Warum die Schallgeschwindigkeit in Flüssigkeiten mit T sinkt

Schall ist eine elastische Welle: v = √(K_eff/ρ). K_eff entsteht allein aus den intakten Bindungen (ΒΙΝDUΝG) — sie sind die Federn. Steigt T, brechen kooperativ immer mehr Bindungen (RΙSS); das Federnetz wird dünner und K_eff fällt steiler als die Dichte ρ. Der Bruch dominiert den Zähler, nicht den Nenner — deshalb sinkt v. Molekular gedeutet: Die Schallübertragung über Molekülbindungen ist schneller als die mechanische Schallweiterleitung durch mechanische Stöße. Je mehr Risse in einer Flüssigkeit bei zunehmender Temperatur auftreten, desto häufiger muss der Schall über mechanische Stöße weitergetragen werden als über Molekülbindungen. Entsprechend wird die Geschwindigkeit des Schalls mit steigender Temperatur langsamer. (Die Anomalie des Wassers — ein Schall-Maximum bei ~74 °C — behandelt Kap. 6.4.)

3.4 Wie die Brownsche Bewegung entsteht — mit geschärftem Mizellen-Begriff

Im statistisch-thermodynamischen Sinn ist eine Mizelle hier kein Tensidaggregat, sondern ein durch RΙSS umrandeter Block: Die fehlenden Bindungen (die Nicht-Bindungen an der Blockgrenze) bilden die „Schale" um ein eingeschlossenes Kriställchen (den intakt gebundenen Kern). Dieser Kern kann eine innere kooperative Rotationsbewegung ausführen (DRΕΙΚLΑΝG-c) — der Block dreht sich als Ganzes, während seine Grenzschale ihn von den Nachbarn entkoppelt.

Die Brownsche Bewegung entsteht nun so: Durch Austausch von Molekülbindungen (RΙSS, FLΙΝΚ) können einzelne Moleküle die Schale einer Mizelle verlassen und in eine andere Mizelle überwechseln. Jeder solche Wechsel überträgt Impuls (ΕRΗΑLΤUΝG; wegen FRΕΙFLUG verlustfrei) und verändert Form und Drehimpuls beider beteiligter Mizellen sprunghaft. Ein großes Testteilchen (Staubkorn) wird von vielen solcher unkorrelierter Mizellen-Umbauten umspült; die Summe der Impulsstöße ist im Mittel null, hat aber endliche Varianz → Random Walk. Je höher T, desto häufiger die RΙSS-Ereignisse und Mizellenwechsel, desto stärker das Zittern. Die Brownsche Bewegung ist damit der makroskopische Schatten oder Folge des mikroskopischen Mizellen-Austauschprozesses.

3.5 Warum Reaktionen an Mizellengrenzschichten schneller werden

Peinliche Vorbemerkung Die nachfolgende Argumentation ist eine Ableitung aus dem Bild und bisher meines Wissens nach nicht empirisch belegt oder untersucht, wobei ich gestehen muss, dass ich mir nicht die Mühe gemacht habe, dazu genau zu recherchieren. Anhaltspunkt für die Neuheit ist lediglich, dass KI-Systeme diese Folgerung nicht erwähnen. Diese These ist ein Beispiel dafür, wie bildhafte Vorstellungen zu neuen Fragestellungen führen können.

In mäßig löslichen Lösungsmitteln lagern sich die Reaktanden bevorzugt an den Grenzschalen der Mizellen an — im Modell gerade die Riss-/Blockgrenzen (RΙSS). Zwei Effekte beschleunigen die Reaktion: (i) Dimensionsreduktion — die Reaktanden suchen einander nicht mehr im vollen Volumen, sondern entlang einer quasi-zweidimensionalen Grenzfläche; ein 2D-Suchlauf trifft statistisch schneller als ein 3D-Suchlauf. (ii) Konzentrierung — die Grenzschale sammelt die Reaktanden auf, ihre lokale Dichte steigt. Bildhaft wirken die Rissgrenzen wie Leitschienen, die die Partner zueinander führen.

Der Dimensionseffekt hat einen strengen Kern: die zufällige Irrfahrt (Pólya). Warum ein zweidimensionaler Suchlauf statistisch besser trifft als ein dreidimensionaler, folgt nicht nur aus der Anschauung, sondern aus einem Satz der Wahrscheinlichkeitsrechnung. George Pólya (1887–1985) zeigte 1921 für die zufällige Irrfahrt (den „Random Walk") auf einem Gitter: In einer und in zwei Dimensionen kehrt der Läufer mit Wahrscheinlichkeit 1 (also sicher) unendlich oft zu jedem Punkt zurück, im ebenen Gitter ℤ² ebenso wie auf der Geraden ℤ. Ab drei Dimensionen gilt das nicht mehr: Die Rückkehrwahrscheinlichkeit ist echt kleiner als 1 (auf dem einfachen kubischen Gitter ℤ³ nur etwa 0,34). Anschaulich, in Pólyas Bild: „Ein betrunkener Mensch findet sicher irgendwann nach Hause, ein betrunkener Fisch hingegen nicht." Übertragen auf zwei Reaktanden, die sich für einen erfolgreichen Stoß mehrfach treffen müssen, heißt das: In einem sehr gut lösenden Lösungsmittel mit freier dreidimensionaler Irrfahrt driften einmal getrennte Partner mit endlicher Wahrscheinlichkeit für immer auseinander — wiederholte Reaktionsversuche misslingen, und die Ausbeute bleibt aus diesem rein geometrisch-statistischen Grund gering. An einer (quasi-)zweidimensionalen Mizellengrenzschicht dagegen ist die Wiederbegegnung sicher (Wahrscheinlichkeit 1): Die Partner finden immer wieder zusammen, bis die Reaktion gelingt. Die Grenzschicht wirkt also nicht nur konzentrierend, sondern stellt zugleich die für Wiederbegegnungen günstigere Dimension bereit. Die (bewusst zugespitzte) Behauptung lautet: Ausbeuten in sehr gut lösenden Lösungsmitteln sind aus Irrfahrt-Gründen tendenziell schlechter.

Beleg: G. Pólya (1921), Über eine Aufgabe der Wahrscheinlichkeitsrechnung betreffend die Irrfahrt im Straßennetz, Math. Ann. 84, 149–160 — hier zitiert nach M. Eisermann, Höhere Mathematik 3 (Univ. Stuttgart), Aufgabe W-2×2. (Original angesurft 5.9.2026).

3.6 Das dreidimensionale Flüssigkeitsfenster

Anders als die Abbildung in 3.1–3.5 zeigt diese Abbildung keine Molekülszene, sondern das Ergebnis der Rechnung aus Kap. 7.9/7.10 als Schema auf einer τ-Achse (Einheiten der Bindungsenergie ε). Das Ziel war, sichtbar zumachen, warum das 2D-Einzelgitter kaum einen Bereich für eine flüssige Phase zulässt, während dies beim 3D-Gitter problemlos möglich ist, wobei für beide Fälle dieselbe Rigiditätsregel (z̄_c = 2d) gilt.

Bild ① — 2D (kein Fenster). Auf der τ-Achse liegen τ_m = 1,20 und τ_triple = 1,16 praktisch übereinander (τ_triple unter τ_m), τ_c = 1,50 dicht darüber: der Bereich für die Flüssige Phase ist auf einen Spalt zusammengedrückt — fest geht quasi direkt in Gas über (nahe der Sublimation).

Bild ② — 3D (sauberes Fenster). Dieselbe Achse, aber gedehnt: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00. Zwischen τ_triple und τ_c liegt das eigenständige flüssige Feld (breit nach unten, τ_triple − τ_m = 1,63; eng nach oben, τ_c − τ_triple = 0,18). Der schmale obere Rand ist als solcher markiert.

Bild ③ — die Ursache. z̄_c = 2d wächst mit der Dimension (4 → 6), damit sinkt v_solid (0,50 → 0,667 in reduzierten Volumina, d. h. der Feststoff wird „lockerer" relativ zur Flüssigkeit), und τ_c = z_max/4 verdoppelt sich (1,50 → 3,00). Der kritische Punkt rückt weit über den Schmelzpunkt. Erst dadurch entsteht Platz für die flüssige Phase. Es ist ein reiner Dimensionseffekt.

4 — Statistisch-mechanisches Modell (unverändert aus Fassung 09)

Zusammenfassung

Kapitel 4 stellt das statistisch-mechanische Gerüst des Modells in 2D auf (Dreiecksgitter, z_max = 6):

• 4.1 Freiheitsgrade und Variablen: Ordnungsgröße ist der Bindungsanteil x ∈ [0,1] mit der mittleren Koordinationszahl z̄ = z_max·x; je Block ≥ 2 kommt ein Rotations-Freiheitsgrad hinzu, die Temperatur τ zählt in Einheiten der Bindungsenergie ε.
• 4.2 Kooperative Zweizustands-Bindung: die selbstkonsistente Bindungsgleichung x = 1/(1 + G·exp[−(ε + K·x)/τ]) (K = J·z_max, G = z_max) folgt als Stationaritätsbedingung der Bragg-Williams-freien Energie; ihr globales Minimum liefert den physikalischen Zweig, ob dieser stetig oder mit Sprung (erster Ordnung) verläuft, entscheidet Kap. 5.5.
• 4.3 Energiebilanz mit Kooperativterm: U/N = τ·[1 + z̄/4 + ½ n_blk] − (ε/2)·z̄ − (K/2)·x² (2D) vereint Translation/Bindungsschwingung/Blockrotation, Bindungs-Konfigurationsenergie und die neue kooperative Molekularfeldenergie −(K/2)x². Die zugehörige Wärmekapazität c_v erhält den positiven Zusatzterm −K·x·ẋ, der den Schmelzpeak verschärft, Fest-Plateau und Gasast aber unberührt lässt.

4.1 Freiheitsgrade und Variablen

• N Moleküle auf einer festen Referenz-Koordinationsgeometrie. Kanonischer 2D-Fall: Dreiecksgitter, z_max = 6, ein Rotations-DOF je Block (VΙΕLΗΕΙΤ, ΒΙΝDUΝG). Dieses Kapitel gilt explizit für die 2D-Interpretation; die 3D-Verallgemeinerung steht in Kap. 6.
• Bindungsvariable s_ij ∈ {0, 1} je geometrischem Nachbarpaar (1 = intakt).
• Koordinationszahl z_i = Σ_j s_ij; Mittelwert z̄ = 2 N_bond / N — der Ordnungsparameter. Es ist z̄ = z_max · x mit dem Bindungsanteil x = p_b ∈ [0, 1].
• Energieeinheit τ ≡ k_B T (klassisch, ½τ pro quadratischem Freiheitsgrad). Alle Energien in Einheiten der Bindungsenergie ε.

4.2 Kooperative Zweizustands-Bindung (Ising-artig, Parameter J)

Eine isolierte Bindung (Energie −ε intakt, 0 gebrochen) hat im Maxwell–Boltzmann-Gleichgewicht p_b⁰(τ) = 1/(1 + e^(−ε/τ)). Der kooperative RΙSS (Axiom 6) verschiebt die effektive Bindungsstabilität um ein Molekularfeld proportional zum Bindungsanteil x. In Molekularfeld-(Ising-)Näherung wird die selbstkonsistente Bindungsgleichung:

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

Dabei ist:

• K = J·z_max die kooperative Kopplungsstärke (RΙSS-Rückkopplung: mehr Nachbarbindungen → höhere effektive Bindungsstabilität). J > 0 schärft den Übergang; oberhalb eines kritischen J wird z̄(τ) unstetig — das ist die latente Wärme des Schmelzens. Der reine Zweizustandsfall ist der Grenzfall J = 0 (kein Molekularfeld, stetiger Crossover).
• G > 1 die Entartung des gebrochenen Zustands: Ein freigewordenes Molekül gewinnt Phasenraum (Translationsfreiheit). G treibt den Bindungsbruch bei hohem τ und sorgt dafür, dass x → 1/(1+G) statt bei ½ hängen zu bleiben (Grenzfall Gas). Kanonisch G = z_max.

Diese Gleichung folgt als Stationaritätsbedingung einer Bragg-Williams-freien Energie je Molekül (n_slots = z_max/2 Bindungsplätze):

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

Der Term −(K/2)x² ist die kooperative (ferromagnetische) Kopplung, der letzte Term die Bruchentartung. Das globale Minimum von g über x liefert den physikalischen Zweig; wo zwei Minima gleich tief werden, liegt der Übergang erster Ordnung. FLΙΝΚ (Axiom 7) rechtfertigt, dass x(τ) zu jeder Temperatur den Gleichgewichtswert annimmt.

4.3 Energiebilanz — mit kooperativem Term

Der Molekularfeld-Term −(K/2)x² ist eine echte Energie (τ-unabhängig) und muss in U/N erscheinen. Die konsistente Bilanz lautet damit (2D):

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

• τ·[…]: Translation (1) + Bindungsschwingung (z̄/4) + Blockrotation (½ n_blk) — DRΕΙΚLΑΝG.
• −(ε/2)·z̄: Bindungs-Konfigurationsenergie — ΒΙΝDUΝG.
• −(K/2)·x²: kooperative Molekularfeldenergie — RΙSS. Neu in U/N.

Die zugehörige Wärmekapazität c_v = d(U/N)/dτ (in Einheiten k_B) erhält einen zusätzlichen Beitrag:

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

wobei der neue Term −K·x·ẋ positiv ist (ẋ < 0) und exakt dort wirkt, wo viele Bindungen zugleich brechen. Er verschärft den Schmelzpeak, ohne das Fest-Plateau (ẋ ≈ 0) oder den Gasast zu verändern.

5 — Zustandssummen und Wärmekapazitäten

Zusammenfassung

Kapitel 5 rechnet Zustandssummen und Wärmekapazitäten aus und liefert die Modelldaten für die Kalibrierung (Kap. 7):

• 5.1 Faktorisierte Zustandssumme: wegen der Zeitskalentrennung faktorisiert Z = Z_trans·Z_vib·Z_rot·Z_bond.
• 5.2 Wärmekapazität je Phase: c_v = d(U/N)/dτ; der positive Term −(ε/2)(dz̄/dτ) erzeugt den Schottky-artigen Peak als Analogon der latenten Wärme.
• 5.3 Numerische c_v(τ)-Kurve: c_v-Peak 5,51 bei τ_m = 1,20, Fest-Plateau 2,50, Gasast → 1; die latente Wärme L ≈ 1,41 ε/Molekül (mit Kooperativterm) ist gegenüber 0,48 ε (ohne) verdreifacht.
• 5.4 / 5.4a Zwei getrennte Schwellen: die Konnektivitäts-Perkolation z̄_c = z_max·p_c ≈ 2,08 (bei τ_b ≈ 1,85) markiert den Cluster-Zerfall (Sieden-nah) und ist ein Ergebnis; sie ist von der generischen Maxwell-Rigidität z̄_c = 2d (Schmelzdichte, Steifigkeitsverlust) sauber zu trennen.
• 5.5 Ordnung des Schmelzübergangs: Kriterium ln G > 2 + 4ε/K; 2D (G = 6) ist für jedes J stetig, 3D (G = 12) erster Ordnung ab K > ~8,2 (Demo: Δz̄ ≈ 5,3, ΔU ≈ 1,75 ε).
• 5.6 Finite-Size-Studie: p*(L) konvergiert monoton von unten gegen p_c = 0,3473; z̄_c ≈ 2,08 ist finite-size-robust.
• 5.7 Voller 3D-c_v-Scan: ab K > K_krit ≈ 8,25 (J > 0,68) öffnet sich der echte Sprung Δz̄ mit latenter Wärme ΔU; der Schmelzpeak wandert mit J nach höheren τ (positive Verschiebung +0,5 … +2,2 gegenüber dem stetigen 2D-Peak).

5.1 Faktorisierte Zustandssumme

Wegen der Zeitskalentrennung (FLΙΝΚ) faktorisiert die kanonische Zustandssumme:

Z = Z_trans · Z_vib · Z_rot · Z_bond.

• Z_trans = (1/N!)·(A/λ_T²)^N, λ_T = h/√(2π m τ) (klassisch-ideale 2D-Translation; VΙΕLΗΕΙΤ, ΕRΗΑLΤUΝG, FRΕΙFLUG).
• Z_vib = (τ/ħω₀)^(N_bond) (klassischer harmonischer Oszillator je intakter Bindung; DRΕΙΚLΑΝG-b).
• Z_rot = ∏_(Blöcke≥2) (2π I_b τ)^(1/2)/ħ (klassischer 2D-Starrrotator je Block; DRΕΙΚLΑΝG-c).
• Z_bond = Σ_{s} exp[(ε/τ)·Σ s_ij + (K/τ)·Kooperativterm] (Ising-artige Bindungssumme; ΒΙΝDUΝG, RΙSS). Molekularfeld-Auswertung liefert x(τ), z̄(τ) aus 4.2.

5.2 Wärmekapazität je Phase (analytisch)

Aus U/N (Kap. 4.3), c_v = d(U/N)/dτ (in Einheiten k_B):

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

Der Term −(ε/2)·(dz̄/dτ) ist positiv (dz̄/dτ < 0: Erwärmen bricht Bindungen) und dominiert dort, wo viele Bindungen zugleich brechen — der Schottky-artige Peak, das Analogon der latenten Wärme.

5.3 Numerische c_v(τ)-Kurve mit Kooperativterm — latente Wärme quantifiziert

Selbstkonsistente Bindung (global-min), n_blk(τ) aus Bond-Perkolations-Monte-Carlo (Union-Find, L = 100, 5 Realisierungen), c_v = dU/dτ (zentrale Differenzen), U/N inklusive −(K/2)x² (cv_scan_09.py, ε = 1, J = 0,55, G = z_max = 6). Ausgewählte Stützstellen:

Ergebnis (numerisch). Die Drei-Phasen-Ordnung bleibt bestätigt und ist nun schärfer:

c_v(fest) ≈ 2,50 < c_v(flüssig, Maximum) ≈ 5,51 > c_v(Gasast) ↓ → 1.

Latente Wärme (Peakfläche über der Basislinie). Integriert man (c_v − Basislinie) über das Schmelzfenster τ ∈ [0,70 ; 1,70] (Basislinie = Gerade durch c_v(0,70) und c_v(1,70)), ergibt sich

L ≈ 1,41 ε/Molekül (mit Kooperativterm), L ≈ 0,48 ε/Molekül (ohne Kooperativterm, alte Bilanz), ΔL ≈ +0,92 ε/Molekül.

Der Kooperativterm hebt den Peak um Faktor ≈ 1,7 (3,25 → 5,51) und verdreifacht die latente Wärme. Physikalisch: das gleichzeitige Brechen vieler Nachbarbindungen (RΙSS) speichert nicht nur Konfigurationsenergie (−ε/2·z̄), sondern auch Molekularfeldenergie (−K/2·x²) — beide werden beim Schmelzen frei.

5.4 Konnektivitäts-Perkolationsschwelle als Ergebnis

Der Riesencluster-Anteil P_∞(τ) kollabiert bei τ ≈ 1,85; dort ist x ≈ 0,347 = p_c(Dreieck) = 2·sin(π/18) = 0,3473, also z̄_c = z_max·p_c ≈ 2,08 — Ergebnis, nicht gesetzt. Dies ist die Konnektivitäts-Schwelle (durchgehender Cluster) und markiert den Cluster-Zerfall nahe dem Sieden (τ_b); sie ist nicht mit der Maxwell-Rigidität (z̄_c = 2d, Schmelzen) zu verwechseln — vgl. die saubere Trennung in Kap. 5.4a.

5.4.a Zwei verschiedene Schwellen sauber getrennt

Die bisher unscharfe Doppelbezeichnung „Perkolation/Rigidität" wird jetzt aufgelöst. Es handelt sich um zwei physikalisch verschiedene Schwellen, die nicht verwechselt werden dürfen:

Begründung. Das Schmelzen ist der Verlust der mechanischen Steifigkeit eines dichten Blocks, nicht der Zerfall der Verbundenheit. Die relevante Größe ist daher die Maxwell-Zwangsbedingungs- zählung (mean-field / generische Rigidität): ein D-dimensionales Netz aus Zentralkräften wird genau dann isostatisch, wenn die mittlere Koordination z̄ = 2D erreicht (2D → 4, 3D → 6). Die Bond-Konnektivität (z̄_c = 2,08) liegt deutlich darunter; ein bei z̄ = 2,08 gerade noch verbundenes Netz ist längst floppy (weiche Moden), also bereits flüssig. Deshalb markiert die Konnektivität den späteren Cluster-Zerfall nahe dem Sieden, während die Rigidität die frühere Verfestigung/Verflüssigung (Schmelzdichte) setzt. Diese Trennung fixiert den Tripelpunkt in Kap. 7.8 parameterfrei und realistisch; sie ändert nicht die in 5.4/5.6 bestätigte Konnektivitätsschwelle 2,08, sondern nur deren Deutung (τ_b = Sieden-nah, nicht Schmelzen).

5.5 Ordnung des Schmelzübergangs

Analytisches Kriterium. Das Spinodalgebiet (zwei Minima von g) existiert für τ < K/4; der Übergang liegt näherungsweise bei τ_t ≈ (ε + K/2)/ln G. Ein echter Sprung (Übergang erster Ordnung) verlangt τ_t < K/4, also im Grenzfall K → ∞:

ln G > 2 + 4ε/K.

Daraus folgt ein bemerkenswertes Resultat über die Bruchentartung G = z_max:

Numerische Bestätigung (max. |Δx| je τ-Schritt, detect_first_order):

Deutung. Im kanonischen 2D-Modell (G = 6) ist das Schmelzen prinzipiell stetig (ein „weicher" Übergang / Crossover) — für kein J entsteht ein echter Sprung, weil ln 6 < 2. Erst die höhere Koordination/Entartung der 3D-Packung (G = 12) erlaubt ab hinreichender Kooperativität (K > 8,2) einen Übergang erster Ordnung mit endlicher latenter Wärme. Das deckt sich mit der Realität: 2D-Schmelzen ist häufig kontinuierlich (KTHNY-artig), 3D-Schmelzen ist erster Ordnung. J = 0 (kein Molekularfeld) ist stets der stetige Referenzfall. — Physikalisch deutet J das Verhältnis von Kooperativ- zu Bindungsenergie (K/ε = J·z_max/ε): wie stark eine intakte Nachbarbindung die eigene stabilisiert.

First-Order-Demonstration (3D-motiviert, cv_scan_09.py, J = 1,0, G = 12, z_max = 12): z̄ springt um Δz̄ ≈ 5,30 bei τ ≈ 2,82; die latente Wärme als Energiesprung ΔU ≈ 1,75 ε/Molekül. Derselbe K auf einem G = 6-Gitter (2D) bleibt stetig — der Ordnungscharakter hängt allein an G = z_max.

5.6 Finite-Size-Studie der Perkolationsschwelle

Bond-Perkolation auf L × L Dreiecksgittern (periodisch), Riesencluster-Anteil P_∞(p) über L skaliert (cv_scan_09.py; L = 50/100/200/400 mit 8/5/3/2 Realisierungen). Exakte Schwelle p_c = 2·sin(π/18) = 0,34730 → z̄_c = z_max·p_c = 2,0838.

P_∞(p) (Mittel ± Streuung):

Schwellenschätzung p*(L) (P_∞ kreuzt 0,5):

Ergebnis. p*(L) konvergiert monoton von unten gegen p_c = 0,3473; die Streuung bei p_c fällt von ±0,103 (L = 50) auf ±0,017 (L = 400) — der Übergang schärft sich mit L, verschiebt die Schwelle aber nicht (universeller Perkolationswert). Damit ist z̄_c ≈ 2,08 = z_max·p_c als finite-size-robuster Modellwert bestätigt.

5.7 Voller 3D-c_v(τ)-Scan des Schmelzens erster Ordnung

Der einzelne First-Order-Demonstrationslauf aus Kap. 5.5 (J = 1,0, G = 12) wird jetzt zu einem vollen c_v(τ)-Scan über die Kooperativität J ausgebaut (firstorder_3d_16.py, FCC-Konnektivität für n_blk, 3D-Energiebilanz U/N = τ·[3/2 + z̄/4 + (3/2)n_blk] − (ε/2)z̄ − (K/2)x²). Das analytische Ordnungskriterium aus Kap. 5.5 (erster Ordnung ⇔ K > K_krit = 4ε/(ln 12 − 2) = 8,25) trennt die Fälle; der stetige 2D-Referenzpeak liegt bei τ_2D = 1,20 (c_v ≈ 5,5).

*stetiger Fall: Δz̄/ΔU sind der steile Crossover (kein echter Sprung); τ_Übergang ist dort der c_v-Peak.

Ergebnis. (i) Der stetige 3D-Fall J = 0,55 liefert τ_Übergang = τ_m ≈ 1,72 — die Querprobe zur unabhängigen Bindungskinetik (bondkinetics_3d_15.py: τ_m(FCC) = 1,70) stimmt. (ii) Ab K > K_krit ≈ 8,25 (J > 0,68) öffnet sich der echte Sprung Δz̄ mit endlicher latenter Wärme ΔU (Energiesprung über die Diskontinuität) — das Schmelzen wird erster Ordnung, exakt wie das Kriterium ln 12 > 2 fordert. (iii) Die Querprobe J = 1,0 reproduziert den Kap.-5.5-Demowert Δz̄ ≈ 5,4. (iv) Sprunghöhe Δz̄ und die Lage τ_Übergang wachsen mit J: der Übergang wandert nach oben — die höhere Koordination der 3D-Packung und die stärkere Kooperativität heben die Schmelztemperatur. Gegenüber dem stetigen 2D-Peak (τ_2D = 1,20) ergibt sich damit durchweg eine positive Peak-Verschiebung (+0,5 … +2,2), die mit K monoton zunimmt. Der First-Order-Charakter des 3D-Schmelzens ist damit nicht mehr nur an einem Punkt, sondern über den ganzen Kooperativitätsbereich belegt.

6 — Mathematische Plausibilität der 3D-Interpretation

Zusammenfassung

Kapitel 6 macht die Verallgemeinerung des Bindungsbildes von 2D auf 3D rein mathematisch plausibel (ohne eigene HTML-Visualisierung):

• 6.1 Rotationsfreiheitsgrad 2D → 3D: die Blockrotation wächst von einer auf drei Achsen, die Translation von zwei auf drei Richtungen, z_max von 6 auf 12; die c_v-Grundwerte steigen (c_v,fest 2,5 → 4,5, c_v,gas 1 → 3/2), die Form der c_v(τ)-Kurve bleibt.
• 6.2 Konnektivität als Unterscheidungsmerkmal: z_max ist eine quantisierte Größe und trennt die dichteste Packung (FCC/HCP, z_max = 12) sprunghaft von BCC (z_max = 8), ohne Zwischenzustand.
• 6.3 Zusätzlicher Übergang fest(1)→fest(2): aus der Diskretheit folgt die Polymorphie FCC↔BCC mit parameterfreier c_v-Stufe 1,0 k_B; im strengen Idealbild (θ_BCC = θ_FCC) tritt sie nicht auf, erst ein weicheres BCC (θ_BCC < θ_FCC unter der Öffnungsschwelle 0,61) öffnet τ_ss < τ_m.
• 6.4 / 6.5 Überlegungen: die Schallanomalie des flüssigen Wassers als verschmierter Nahordnungswechsel zweier Gittertypen und eine fraktale Wirbelfeld-Alternative zur Mizellenrotation (rein qualitativ).
• 6.6 / 6.8 / 6.9 Mikroskopisches θ(Gitter) auf drei Stufen: Einstein-Federmodell (θ_BCC/θ_FCC ≈ 0,82…0,94), Debye-artige Phononen-DOS (≈ 0,87…1,00) und Lennard-Jones-Vollgitter (≈ 1,0; ideales BCC scherinstabil, FCC stabiler Grundzustand) — alle über der Schwelle 0,61, reale Polymorphie liegt also außerhalb des harmonischen Idealbildes.
• 6.7 Metastabilität (2D): Größenordnungsrechnung für den Umbau trianguliert ↔ quadratisch (ΔU = +1 ε, ΔA ≈ +15,5 %) samt Keimbildungsbarriere ΔG\* = πλ²/Δμ_A als kinetischem Zusatz.

Dieses Kapitel macht die 3D-Verallgemeinerung rein mathematisch plausibel. Die 3D wird nie in einer HTML-Version visualisiert. Ziel ist allein, zu zeigen, wie sich der Rotationsfreiheitsgrad und die Konnektivität von 2D nach 3D ändern und welche neuen Übergänge dadurch denkbar werden.

6.1 Änderung des Rotationsfreiheitsgrads 2D → 3D

Der einzige qualitativ neue Punkt gegenüber 2D ist die Rotation der Blöcke:

• 2D: ein starrer Block hat 1 Rotations-DOF (Drehung in der Ebene) → U_rot/N = ½ n_blk τ.
• 3D: ein starrer Block hat 3 Rotations-DOF (drei Euler-Achsen) → U_rot/N = (3/2) n_blk τ.

Zugleich wächst die Translation von 2 auf 3 DOF (U_trans/N = (3/2)τ) und die maximale Koordinationszahl von z_max = 6 (Dreieck) auf z_max = 12 (dichteste Kugelpackung). Die Energiebilanz in 3D lautet damit

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

und die Grundwerte der Wärmekapazität verschieben sich konsistent:

Die Struktur der Kurve (Solid-Plateau, Flüssig-Maximum durch −(ε/2)dz̄/dτ, Abfall zum Gaswert) bleibt identisch; nur die Niveaus skalieren mit Dimension und z_max.

6.2 Einfluss verschiedener Kristallgitter — Konnektivität als Unterscheidungsmerkmal

In 3D sind mindestens zwei ideale Gitter zu unterscheiden, die sich in ihrer Konnektivität (Koordinationszahl z_max) trennen:

Die Konnektivität z_max ist das entscheidende Unterscheidungsmerkmal: Sie legt sowohl das Wärmekapazitätsniveau des Festkörpers (z_max/4) als auch die Perkolationsschwelle z̄_c fest. Weil jedes Gitter nur eine feste Koordinationszahl zulässt, ist die Konnektivität eine quantisierte, endliche Größe — es gibt keinen kontinuierlichen Übergang zwischen z_max = 12 und z_max = 8, sondern nur den diskreten Sprung von einem Gittertyp zum anderen.

6.3 Der zusätzliche Übergang fest(1) → fest(2)

Aus der Endlichkeit/Diskretheit der Konnektivität folgt ein im 2D-Modell nicht vorhandener Übergang: fest(1) → fest(2) (Polymorphie/Allotropie). Zwei Festkörperphasen mit verschiedenem z_max (etwa FCC ↔ BCC) haben verschiedene freie Enthalpien μ(τ, p); bei der Temperatur/Druck-Kombination μ_{fest(1)} = μ_{fest(2)} wechselt der Stoff sein Gitter, ohne zu schmelzen. Im Modellbild ändert sich dabei die Koordinationsgeometrie sprunghaft (nicht nur die Belegung z̄ wie beim Schmelzen) — die einzige Situation, in der die sonst festgehaltene Referenzgeometrie selbst umspringt. Beobachtbar wäre eine kleine c_v-Stufe (Änderung von z_max/4) plus latente Umwandlungswärme, deutlich schwächer als der Schmelzpeak.

Das First-Order-Kriterium ln G > 2 + 4ε/K (Kap. 5.5) liefert nun einen konkreten Grund, warum gerade die 3D-Interpretation qualitativ reicher ist: Nur dort (G = z_max = 12, ln 12 = 2,485 > 2) kann das Schmelzen erster Ordnung werden — und nur dort existiert der diskrete Gitterwechsel fest(1) → fest(2) (FCC ↔ BCC), weil die Konnektivität quantisiert/endlich ist. Beide Phänomene sind dieselbe Ursache in zwei Gewändern: die höhere, diskrete Koordination der 3D-Packung erzeugt sowohl den scharfen Schmelzsprung als auch die Polymorphie. Das 2D-Einzelgitter (G = 6) hat keines von beiden — konsistent mit Kap. 5.5.

Quantifizierung (bondkinetics_3d_15.py). Der Übergang wird jetzt schlank durchgerechnet. Beide Gitter erhalten je eine freie Enthalpie im klassischen harmonischen Bild (gemeinsame p·v- und Translationsanteile fallen im Vergleich heraus):

μ(τ) = h₀ − c_v·τ·ln(τ/θ), h₀ = −(ε/2)·z_max, θ = ħω/k_B (Phononensteifigkeit in ε-Einheiten),

mit den parameterfreien Kennwerten aus Kap. 6.1/6.2:

Daraus die c_v-Stufe parameterfrei: |Δc_v| = z_max/4-Änderung = (12−8)/4 = 1,0 k_B (der Festkörper wird beim Übergang FCC→BCC um 1 k_B „leichter" in der Wärmekapazität) — viel kleiner als der Schmelzpeak (mehrere k_B, Kap. 5.3/7.11).

Die Umwandlungstemperatur τ_ss (μ_FCC = μ_BCC) und die latente Umwandlungswärme L = τ_ss·(S_BCC − S_FCC):

Ergebnis (ehrlicher Befund). Im strengen Idealbild (dieselbe Phononensteifigkeit für alle Bindungen, θ_BCC = θ_FCC) gibt es keinen fest(1)→fest(2)-Übergang: FCC gewinnt sowohl die Enthalpie (tiefer, −6 vs. −4 ε) als auch — wegen der höheren Bindungszahl (c_v = 4,5 > 3,5) — die Schwingungsentropie; μ_FCC bleibt bei jeder Temperatur unter μ_BCC. Erst eine gitterabhängige Phononensteifigkeit θ_BCC < θ_FCC (das offenere BCC schwingt weicher, hat also mehr Entropie) hebt μ_BCC bei hoher Temperatur unter μ_FCC — dann öffnet sich die Umwandlung, und sie liegt unter dem Schmelzpeak τ_m, sodass die Reihenfolge fest(1)→fest(2)→flüssig→gas entsteht (τ_ss < τ_m < τ_c ≤ τ_b(3D), Kap. 7.11). Die weiche BCC-Steifigkeit ist eine milde Erweiterung außerhalb des strengen Idealbildes (real bekannt als Eisen-γ/δ-Polymorphie: BCC bei hoher Temperatur, FCC dazwischen); die c_v-Stufe 1,0 k_B und die kleine positive Umwandlungswärme sind dagegen modellintern und parameterfrei.

Wie klein ist θ_BCC/θ_FCC wirklich? ( theta_springs_16.py). Um die Umwandlung vom freien θ-Parameter zu befreien, wird die Phononensteifigkeit jetzt mikroskopisch abgeschätzt (Kap. 6.6): ein Einstein-Federmodell mit zentralen Paarfedern liefert θ_BCC/θ_FCC ≈ 0,82 (nur nächste Nachbarn) bis ≈ 0,94 (mit BCC-Zweitschale). Die aus dem μ-Vergleich nötige Öffnungsschwelle liegt dagegen bei θ_BCC/θ_FCC < 0,61 (darunter existiert τ_ss und liegt unter τ_m). Das geometrische Steifigkeitsverhältnis ist also zu groß: im strengen Idealbild bleibt FCC bei jeder Temperatur stabil, es gibt keinen rein geometrisch getriebenen fest(1)→fest(2)-Übergang. Reale Polymorphie (Eisen γ↔δ) benötigt zusätzliche (elektronische/magnetische/anharmonische) Beiträge außerhalb des idealen harmonischen Bildes — genau die „milde Erweiterung", die oben angesetzt wurde. Das bestätigt den ehrlichen Befund und macht das θ-Verhältnis von einer freien Annahme zu einer abgeschätzten, benannten Modellgrenze.

6.4 Überlegungen: die Schallanomalie des flüssigen Wassers

Normale Flüssigkeiten zeigen v(T) monoton fallend (Kap. 3.3). Wasser hat ein Schall-Maximum bei ≈ 74 °C: v steigt zunächst mit T, ehe es fällt. Im Modell lässt sich das als Konkurrenz zweier lokaler Gittertypen (Kap. 6.2) deuten: eine offene, niedrig koordinierte eisähnliche Nahordnung (kleines z_max, viele „Löcher") und eine dichte, höher koordinierte Nahordnung. Bei niedriger T dominiert die offene Struktur; Erwärmen bricht zunächst die offene Nahordnung zugunsten der dichteren — die effektive Konnektivität und damit K_eff steigen (v steigt). Erst oberhalb des Maximums überwiegt der gewöhnliche RΙSS-Bindungsbruch, K_eff fällt, v sinkt. Die Anomalie ist damit Ausdruck eines verschmierten, unvollständigen fest(1)↔fest(2)-artigen Nahordnungswechsels innerhalb der flüssigen Phase — nur im 3D-Modell mit zwei Gittertypen denkbar, nicht im 2D-Einzelgitter.

6.5 Alternative: fraktale Wirbelfelder statt Mizellenrotation

Die lokale Brownsche Bewegung (Kap. 3.4) wurde über Mizellenrotation und -austausch erklärt. Alternativ lässt sie sich als Effekt fraktaler Wirbelfelder vorstellen: Die durch RΙSS ständig umgebauten Blockgrenzen erzeugen ein Feld lokaler Impulsströme, das — ähnlich turbulenten Strömungen — auf allen Längenskalen selbstähnliche Wirbel aufweist. Ein Testteilchen wird dann nicht von diskreten Mizellenstößen, sondern von einem skaleninvarianten Wirbelkontinuum getrieben. Beide Bilder (diskrete Mizellen ↔ fraktale Wirbel) sind komplementäre Grenzfälle derselben RΙSS-Kinetik; welches trägt, hängt vom Verhältnis Testteilchen-Größe zu mittlerer Mizellengröße ab. Diese Andiskussion ist rein qualitativ und wird hier nicht mathematisch ausgeführt.

6.6 Mikroskopische Phononensteifigkeit θ(Gitter): ein Einstein-Federmodell

Der fest(1)→fest(2)-Vergleich in Kap. 6.3 hängt an der relativen Phononensteifigkeit θ_BCC/θ_FCC. Statt sie frei zu wählen, wird sie hier aus Koordination und Geometrie abgeschätzt (theta_springs_16.py).

Modell. Ein Zentralatom sitzt im harmonischen Napf seiner z Nachbarn; jede Bindung ist eine radiale Paarfeder κ. Eine Auslenkung u projiziert isotrop mit ⟨cos²⟩ = 1/d = 1/3 (d = 3) auf die Bindungsrichtungen, sodass die effektive Federkonstante k_eff = (z/3)·κ ist. Die Einstein-Frequenz ist ω = √(k_eff/m), also θ = ħω/k_B ∝ √k_eff. Bei gleicher Masse und gleicher Paarfeder folgt

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

Verfeinerungen. (B) Die Zweitschale des BCC (6 Atome bei r₂/r₁ = 2/√3 = 1,155) mit abstandsskalierter Feder κ(r) ∝ (r₁/r₂)ⁿ (n = 6…8, Lennard-Jones-artig) hebt das Verhältnis auf 0,91…0,94 (BCC wird durch die nahe Zweitschale effektiv steifer). (C) Die kürzere Bindung des BCC bei gleicher Dichte (r_nn(BCC)/r_nn(FCC) = 0,972, −2,8 %) steift die Feder leicht an und liefert 0,89…0,91. Insgesamt ein geometrisches Fenster θ_BCC/θ_FCC ≈ 0,82…0,94.

Öffnungsschwelle. Der μ-Vergleich (Kap. 6.3) öffnet die Umwandlung (τ_ss existiert und liegt unter τ_m = 1,70) erst für θ_BCC/θ_FCC < 0,61. Das geometrische Fenster liegt deutlich darüber.

Befund. Das reine Federmodell macht das BCC nicht weich genug: im strengen Idealbild bleibt FCC stabil, ein geometrisch getriebener fest(1)→fest(2)-Übergang entsteht nicht. Damit ist θ_BCC/θ_FCC von einer freien Annahme zu einer abgeschätzten Größe geworden, deren natürlicher Wert (≈ 0,8–0,9) die Umwandlung geschlossen hält — konsistent mit dem ehrlichen Befund aus Kap. 6.3. Reale Polymorphie liegt somit nachweislich außerhalb des idealen harmonischen Bildes (elektronische/magnetische Anteile, Anharmonizität), nicht in der bloßen Gittergeometrie.

6.7 Mathematische Abschätzung der Metastabilität (2D trianguliert ↔ quadratisch)

Zur qualitativen Diskussion aus Kap. 3.2 hier die zugehörige Größenordnungsrechnung im 2D-Bild (harte Scheiben, Durchmesser D, Kontaktabstand a = D).

(1) Energie. Je Teilchen U/N = −(ε/2)·z (jede Bindung hälftig): trianguliert z = 6 → −3 ε, quadratisch z = 4 → −2 ε. Der Umbau trianguliert → quadratisch kostet also ΔU = +1 ε je Teilchen (endotherm; die lockere Form ist energiereicher).

(2) Volumen (Fläche) und Arbeit. Fläche je Teilchen: trianguliert A_tri = (√3/2)·D² ≈ 0,866 D², quadratisch A_quad = D². Die Aufweitung ist ΔA = A_quad − A_tri ≈ 0,134 D² (+15,5 %) — die im Comic gezeigte Volumenvergrößerung. Gegen einen (reduzierten) Druck p ist dafür die Arbeit W = p·ΔA aufzubringen; die freie Enthalpie steigt um Δμ = ΔU + p·ΔA − τ·ΔS. Der Entropieterm ΔS > 0 (die lockere Form hat weichere Moden, größere c_v-Stufe Δc_v = (6−4)/4 = 0,5 k_B) senkt μ_quad erst bei hoher τ unter μ_tri — deshalb ist trianguliert der Tieftemperatur-Grundzustand und quadratisch die Hochtemperatur-/metastabile Form.

(3) Keimbildungsbarriere (2D). Metastabilität setzt eine Barriere voraus. In 2D wächst ein Keim der neuen Phase (Radius r) mit Flächengewinn −Δμ_A·π r² (Δμ_A = Antrieb je Fläche) gegen den Randbeitrag +λ·2π r (λ = Kantenenergie je Länge): ΔG(r) = −Δμ_A·π r² + 2π λ r. Maximum bei r\* = λ / Δμ_A mit Barriere ΔG\* = π λ² / Δμ_A. Der kritische Keim r\ schrumpft mit wachsendem Antrieb Δμ_A; solange die thermische Skala τ = k_BT die Barriere nicht erreicht (**ΔG\ ≳ einige k_BT), bleibt die lockere (bzw. beim Sieden die übersättigte) Phase kinetisch eingefroren. Für kleine Cluster (2–3 Teilchen, wie im Übersättigungs- Strip) ist das Rand-zu-Fläche-Verhältnis maximal ungünstig — sie liegen links von r\* und lösen sich im Mittel wieder auf. Damit ist die im Comic gezeigte Aktivierungsenergie dieselbe Größe ΔG\*, die sowohl die metastabile Kristallform als auch die übersättigte Gasphase trägt; sie ist ein kinetischer** Zusatz zum idealen (barrierefreien) Gleichgewichtsmodell.

6.8 Debye-artige Verfeinerung von θ(Gitter)

Das Einstein-Federmodell (Kap. 6.6) ersetzt die ganze Phononenverteilung durch eine Frequenz. Der vorliegende Ansatz verfeinert θ(Gitter) durch eine schlanke Phononen-Zustandsdichte je Gitter (theta_debye_17.py): dieselben zentralen Paarfedern κ(r) = r^(−n) auf nächster und zweiter Nachbarschale bilden die Dynamikmatrix D(k) = Σ_j κ_j·(R̂_j⊗R̂_j)·(1 − cos(k·R_j)), die auf einem k-Gitter der primitiven Brillouin-Zone (Monkhorst-Pack, 12³ Punkte) diagonalisiert wird (analytischer 3×3-Eigensolver, nur Standardbibliothek). Aus den drei Phononenzweigen folgen zwei Debye-artige Steifigkeitsmaße als Frequenzmomente:

• 2. Moment θ(2) ∝ √⟨ω²⟩ — der Einstein-Grenzfall (die BZ-Mittelung von 1 − cos ergibt 1, also ⟨ω²⟩ = z·κ/3); Kontrolle nn-only, reine Geometrie: √⟨ω²⟩_BCC / √⟨ω²⟩_FCC = 0,816 = √(8/12) — exakt der Einstein-Wert aus Kap. 6.6.
• log-Moment θ(0) ∝ exp⟨ln ω⟩ — die eigentlich Debye-artige Größe, die die weichen akustischen Zweige gewichtet, die der Einstein-Wert ignoriert.

Alle skalenbehafteten Konstanten (ħ/k_B, κ-Skala, Debye-Vorfaktor) fallen im Verhältnis θ_BCC/θ_FCC heraus; nur die DOS-Form bleibt. Mit zweiter Nachbarschale (die das BCC überhaupt erst mechanisch stabilisiert) und Federgesetz n = 0…8 ergibt sich

θ_BCC/θ_FCC ≈ 0,87…1,00 (Fehlerband über Momenten-Definition und Federgesetz).

Das ist eher höher als der reine Einstein-Wert (0,82…0,94), weil die bei gleicher Dichte kürzere BCC-Bindung die Federn ansteift. Die aus dem μ-Vergleich (Kap. 6.3) nötige Öffnungsschwelle 0,61 liegt weiterhin deutlich darunter: auch die Phononen-DOS macht das BCC nicht weich genug. Der Befund aus Kap. 6.6 wird damit durch die feinere Näherung bestätigt und präzisiert — im strengen Idealbild bleibt FCC bei jeder Temperatur stabil; ein rein geometrisch getriebener fest(1)→fest(2)-Übergang entsteht nicht, reale Polymorphie (Eisen, Kap. 7.13) braucht Beiträge außerhalb des idealen harmonischen Bildes.

6.9 θ(Gitter) aus einem Lennard-Jones-Vollgitter

Kap. 6.8 nutzt eine abgeschnittene Zentralfeder κ(r) = r^(−n) auf nn + 2nn. Fassung 17 prüft dieselbe Frage mit einer echten Kraftkonstante aus einem Paarpotential über das volle Gitter (theta_ljlattice_18.py): jede Bindung wird zur Feder aus dem Lennard-Jones-Potential φ(r) = 4ε[(σ/r)¹² − (σ/r)⁶] mit dem longitudinalen Anteil κ_L = φ″(r) und dem transversalen Anteil κ_T = φ′(r)/r. Die Dynamikmatrix

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

summiert über alle Nachbarschalen bis zu einem Cutoff (2,0…3,0 nn); jedes Gitter wird zuvor auf sein eigenes LJ-Gleichgewicht relaxiert (freie Gitterkonstante, nicht gleiche Dichte erzwungen). Das ist der ehrliche Test, ob die schlanke nn+2nn-Näherung (Kap. 6.8) die Steifigkeitsordnung richtig trifft.

Ergebnis (deterministisch, nur Standardbibliothek):

Zwei Befunde, beide verschärfen Kap. 6.6/6.8:

1. FCC ist der stabile LJ-Grundzustand (E ≈ −8,4 ε/Atom bei Cutoff 3,0 nn, gegen −8,0 ε für BCC) — der bekannte Vorrang der dichten Packung im Paarpotential.
2. Das ideale LJ-BCC ist scherinstabil: rund 3…3,5 % der Phononenmoden sind imaginär (ω² < 0, Born-Instabilität). Das reine Paarpotential-Idealbild trägt also gar kein mechanisch stabiles BCC; über die verbleibenden reellen Moden bleibt θ_BCC/θ_FCC ≈ 1,0 — die aus dem μ-Vergleich nötige Öffnungsschwelle 0,61 wird erst recht nicht unterschritten.

Damit ist die Aussage aus Kap. 6.6/6.8 robust gegen die Verfeinerung von der abgeschnittenen Feder zum Vollgitter und vom reinen Longitudinalmodell zur vollen Dynamikmatrix: im klassischen Idealbild entsteht kein rein geometrischer fest(1)→fest(2)-Übergang — reale Polymorphie (Eisen Kap. 7.13, Titan Kap. 7.15) braucht Beiträge außerhalb des Paarpotentials (elektronisch/magnetisch, anharmonisch). Eine echte ab-initio-Phononenrechnung bleibt bewusst außen vor (verließe den schlanken Idealbild-Charakter).

7 — Kalibrierung und Zustandsgleichung

Zusammenfassung. Kapitel 7 kalibriert das Modell über den einen Energieparameter ε auf reale Stoffe und baut den Flüssig-/Gas-Zweig zur Zustandsgleichung aus:

• 7.1 / 7.2 Kalibriermethode und -tabelle: lineare Temperaturanbindung T = (ε/k_B)·τ mit einem Ankerpunkt am Schmelzpunkt; daraus das verfeinerte universelle Verhältnis R = τ_b/τ_m ≈ 1,54.
• 7.3 / 7.4 Hartkern-van-der-Waals-EOS: reduzierte EOS p_r = 8τ_r/(3v_r − 1) − 3/v_r² mit kritischem Punkt (1, 1, 1) und Maxwell-Koexistenz; a ∝ ε, b ∝ D². Die Siedepunkte liegen systematisch bei T_b/T_c ≈ 0,48 (Guldberg-Abstand), eine Redlich–Kwong-Vergleichsrechnung hebt sie auf ≈ 0,60 — als benanntes Fehlerbudget, ohne EOS-Wechsel.
• 7.5 Gemeinsame (τ, v)-Fläche: Koexistenzkuppel (Sieden), Schmelzlinie und harte Wand fügen sich mit einem Tripelpunkt zusammen (τ_r,triple ≈ 0,47 für v_solid = 0,40).
• 7.6 Latente Wärme: L ≈ 1,41 ε ergibt ΔS_fus ≈ 1,175 R (nahe der Richards-Regel); Metalle exzellent getroffen (Kupfer ×1,00), Ionen-/H-Brücken-/Molekülkristalle liegen außerhalb des Idealbildes.
• 7.7 Vereinheitlichung ε ↔ ε_lg: k_B T_c = z_max·ε/4 koppelt beide Skalen (ε_lg = (27/32)·z_max·ε); die 3D-Vorhersage τ_c/τ_m = 2,50 trifft einfache Molekülstoffe, 2D bleibt Sonderfall.
• 7.8 / 7.9 Parameterfreier Tripelpunkt: v_solid folgt aus der Rigiditätsschwelle z̄_c = 2d, nicht aus der Konnektivität; in 3D öffnet sich ein sauberes Flüssigfenster (τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00), in 2D bleibt es geschlossen (dreifach belegte 2D-Sonderrolle).
• 7.10 / 7.11 Echte 3D-Scans: Bond-Perkolation auf FCC/BCC liefert z̄_c(3D) ≈ 1,5–1,6; der τ-Scan der Bindungskinetik zeigt τ_m ≈ τ_rig (Schmelzen = Rigidität, tief) und τ_b(3D) ≈ τ_c (Sieden = Konnektivitäts-Zerfall am kritischen Punkt).
• 7.12 / 7.13 / 7.15 fest(1)→fest(2)-Linie: aus dem Punkt τ_ss wird über den Druck eine Linie (dp/dτ = ΔS/Δv > 0, Δv = +8,9 %), als vierte Grenzlinie in die (τ, v)-Fläche und als (τ, p)-Linie an Eisen (γ→δ) kalibriert, an Titan (α→β) bestätigt und mit Argon als Negativkontrolle.
• 7.14 Redlich–Kwong-Zweig: als eigenständige Idealbild-Variante mit eigener Kalibrierung neben dem van-der-Waals-Default (T_b/T_c ≈ 0,60, Fehler −1 %).

7.1 Kalibriermethode (ein Energieparameter ε) — angepasst

Linear: T = (ε/k_B)·τ, ein Ankerpunkt je Stoff. Verankerung am Schmelzpunkt: ε/k_B = T_m,exp/τ_m, Vorhersage T_b,pred = T_m,exp·(τ_b/τ_m). Mit dem verfeinerten Schmelzpeak (Kooperativterm) liegt das c_v-Maximum jetzt bei τ_m = 1,20 (statt 1,15), die Konnektivitäts-Perkolation (Sieden-nah) unverändert bei τ_b = 1,85. Daraus das verfeinerte universelle Verhältnis:

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

Die Kalibriertabelle aus Fassung 08 (Kap. 7.2) skaliert linear mit R: alle T_b,pred sinken um den Faktor 1,54/1,61 ≈ 0,957 (≈ −4,3 %). Das verbessert die Überschätzung der Edelgase leicht (Argon +54 % → +48 %) und verschlechtert die NaCl-/Eisen-Treffer geringfügig (NaCl −0,6 % → −4,8 %). Die Grundaussage bleibt: R passt für isotrop dicht gepackte, ungerichtet bindende Stoffe; Edelgase/Ketten liegen außerhalb des Idealbildes.

7.2 Kalibriertabelle

Unverändert übernommen aus modell_08.md Kap. 7.2 (mit dem Hinweis der linearen R-Skalierung aus 7.1). Reproduzierbar via cv_scan_08.py (Tabelle) bzw. cv_scan_09.py (τ_m, τ_b, R).

7.3 Dichte-/Volumenachse: Hartkern-van-der-Waals-EOS

Der Druck wird von „Parameter" zur echten Variable gehoben. ΗΑRΤΚΕRΝ liefert das Eigenvolumen b, die Bindungsanziehung (ε·z̄) den Kohäsionsdruck a. In reduzierten Einheiten (kritischer Punkt bei (v_r, τ_r, p_r) = (1, 1, 1)):

p_r = 8 τ_r / (3 v_r − 1) − 3 / v_r², v_r > 1/3 (harte Wand).

Koexistenz per Maxwell-Konstruktion (cv_scan_09.py, Gleichflächenschnitt):

Ergebnis. Unterhalb von τ_c öffnet sich der l/g-Dichtesprung (v_gas/v_liq wächst von 2,5 auf 11,5, je kälter desto größer). Das ist die Siedelinie — nun sauber getrennt vom Schmelz-c_v-Peak (Kap. 5.3), der auf dem Gitter fester Dichte lebt. Damit ist der in Fassung 08 offene Punkt 3 geschlossen: Schmelzen (Rigiditäts-/Bindungsordnung bei fester Dichte) und Sieden (Dichte-/Volumensprung) sind verschiedene Achsen desselben Modells.

7.4 Quantitative l/g-Kalibrierung: a → ε, b → D²

Abbildung der reduzierten Größen auf reale Stoffe. Die van-der-Waals-Konstanten hängen am kritischen Punkt:

b = k_B·T_c / (8 p_c) (Eigenvolumen, ΗΑRΤΚΕRΝ), a = 27 k_B²·T_c² / (64 p_c) (Kohäsion, ΒΙΝDUΝG).

Damit ist a ∝ ε (die Kohäsion stammt aus der Bindungsenergie) und b ∝ D² (das 2D-Eigenvolumen aus dem Hartkern-Durchmesser). Die natürliche Energieskala des l/g-Zweigs ist

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

Ankerstoff Argon (T_c = 150,87 K, p_c = 48,63 bar): b = 5,35·10⁻²⁹ m³/Molekül, a = 3,76·10⁻⁴⁹ J·m³/Molekül², ε_lg/k_B = 509 K und ein Moleküldurchmesser D ≈ b^{1/3} = 3,8·10⁻¹⁰ m — gegenüber dem van-der-Waals-Durchmesser des Argons (~3,4·10⁻¹⁰ m) auf 11 % genau.

Vorhersage der Siedetemperatur. Bei festem Außendruck ist der Siedepunkt dort, wo der reduzierte Sättigungsdruck p_sat(τ_r)·p_c = p_atm erfüllt ist; daraus T_b = T_c·τ_r,b (Maxwell-Sättigungskurve aus lg_calib_10.py):

*CO₂ sublimiert bei 1 atm (T = 194,7 K = Sublimationspunkt).

Ergebnis. Für die molekularen Stoffe (p_c ≈ 30–70 bar) liefert das Modell durchweg T_b/T_c ≈ 0,48–0,51, systematisch ~17–24 % unter dem Messwert. Das ist kein Modellfehler im engeren Sinn, sondern der bekannte Abstand der van-der-Waals-Gleichung zur empirischen Guldberg-Regel (T_b/T_c ≈ 0,6): die van-der-Waals-Isotherme unterschätzt den Sättigungsdruck bei tiefen Temperaturen. Wasser (und die Metalle, deren T_c/p_c ohnehin nur grob geschätzt sind) liegen bei p_r,atm < 0,005 außerhalb des sinnvollen van-der-Waals-Bereichs (τ_r,b < 0,40) — dort ist die Kohäsion nicht mehr van-der-Waals-artig (H-Brücken, metallische Bindung). Die Größenordnung und die Reihung der Siedepunkte werden korrekt getroffen.

Guldberg-Einordnung als Fehlertabelle. Der Abstand wird hier nur benannt, nicht durch einen EOS-Wechsel korrigiert (das Idealbild bleibt van-der-Waals-artig). Die empirische Guldberg-Regel (Guldberg 1890) besagt, dass der Siedepunkt bei etwa T_b/T_c ≈ 0,6 liegt (leichte Stoffabhängigkeit ~0,58–0,64). Das reduzierte van-der-Waals-Modell liefert stattdessen einen konstanten Wert T_b/T_c = τ_r,b ≈ 0,48–0,51:

Deutung. Der Modellmittelwert 0,49 trifft die Guldberg-Konstante 0,60 auf etwa −18 %; die Abweichung ist systematisch und einseitig (immer zu tief), also eine Eigenschaft der van-der-Waals-Isotherme (zu geringer Sättigungsdruck bei tiefen τ), nicht des molekularen Bildes. Eine bessere Kohäsions-EOS (z. B. Redlich–Kwong) würde τ_r,b anheben — bleibt hier aber bewusst außerhalb des Idealbildes.

Fehlerbudget mit einer Redlich–Kwong-Vergleichsrechnung (guldberg_rk_15.py). Um den van-der-Waals-Abstand zur Guldberg-Regel zu beziffern (nicht zu heilen), wird derselbe Siedepunkt-Test mit einer Redlich–Kwong-artigen, temperaturabhängigen Kohäsion a/√T wiederholt. In reduzierten Einheiten (kritischer Punkt (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,

mit Maxwell-Gleichflächenkonstruktion und demselben reduzierten Außendruck p_r,atm wie oben. Ergebnis (ein Zahlenwert je Ankerstoff, T_b/T_c = τ_r,b):

Ergebnis. Die temperaturabhängige RK-Kohäsion hebt T_b/T_c von van der Waals ≈ 0,48 auf ≈ 0,60 — praktisch die Guldberg-Konstante und die realen Werte (≈ 0,58–0,63); der mittlere Fehler gegen die Messwerte fällt von −21 % (van der Waals) auf −1 % (RK). Damit ist die Ursache des Guldberg-Abstands klar beziffert: es ist gerade die fehlende Temperaturabhängigkeit der van-der-Waals-Kohäsion (zu geringer Sättigungsdruck bei tiefen τ), nicht das molekulare Bild. Diese Vergleichsrechnung liegt bewusst außerhalb des Idealbildes; das Modell bleibt van-der-Waals-artig (kein EOS-Wechsel) — das Fehlerbudget ist nur benannt, nicht eingebaut.

7.5 Gemeinsame (τ, v)-Fläche mit Tripelpunkt

Schmelz- und Siedeachse leben bislang getrennt: das Schmelzen (Rigiditäts-/Bindungsordnung) bei fester Dichte, das Sieden (van der Waals) über den Dichtesprung. In einer gemeinsamen (τ_r, v_r)-Ebene fügen sie sich zusammen:

• Koexistenzkuppel (Sieden): v_liq(τ_r) (linker Ast) und v_gas(τ_r) (rechter Ast) aus der Maxwell-Konstruktion, die sich am kritischen Punkt (τ_r, v_r) = (1, 1) treffen.
• Festkörper (Schmelzen): schmales Band knapp oberhalb der harten Wand v_r > 1/3, bei hoher, nahezu fester Dichte (reduziertes Volumen v_solid ≈ 0,40).
• Schmelzlinie: nahezu senkrecht bei v_solid (feste Dichte), die Temperatur der Maxwell-Rigidität (z̄_c = 2d; Kap. 7.8 leitet v_solid daraus parameterfrei ab — die alte Fassung-11-Zuordnung „z̄_c ≈ 2,08" ist damit korrigiert).

Der Tripelpunkt ist der Punkt, an dem die flüssige Koexistenzgrenze v_liq(τ_r) die feste Packungsdichte v_solid erreicht — darunter existiert kein flüssiger Ast mehr (direkte Sublimation):

Ergebnis. Für v_solid = 0,40 liegt der Tripelpunkt bei τ_r,triple ≈ 0,47, p_sat ≈ 0,018, mit einem Dampfvolumen v_gas ≈ 68 (sehr verdünnter Sublimationsdampf). Das Verhältnis τ_triple/τ_c ≈ 0,47 liegt im realen Fenster (Wasser 0,42; Argon 0,56; CO₂ 0,71). Je dichter der Festkörper (kleineres v_solid), desto tiefer der reduzierte Tripelpunkt und desto verdünnter der koexistierende Dampf — qualitativ korrekt. Damit sind beide Übergänge und der Tripelpunkt in einer (τ, v)-Fläche verankert; der kritische Punkt sitzt bei (1, 1, 1). Fassung 15 ergänzt in Kap. 7.12 eine vierte Grenzlinie — die fest(1)/fest(2)-Umwandlung (FCC↔BCC) als zwei nahezu senkrechte Festkörper-Dichtebänder (v ≈ 0,667 bzw. 0,726), die bei τ_ss ineinander umschlagen.

Frühere Version in Der Entwicklung: v_solid = 0,40 war hier noch gesetzt. Kap. 7.8 leitet v_solid parameterfrei aus der Perkolationsschwelle z̄_c ab und macht den Tripelpunkt zur echten Vorhersage.

7.6 Latente Wärme gegen Schmelzenthalpien

Die Modell-Latentwärme L ≈ 1,41 ε (Kap. 5.3) wird über die Schmelzpunktkalibrierung ε/k_B = T_m/τ_m (τ_m = 1,20) in reale Einheiten übersetzt:

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

Das ist — bemerkenswerterweise — fast exakt die Richards-Regel (ΔS_fus ≈ R für einfach gebaute Festkörper). Vergleich mit gemessenen Schmelzenthalpien (lg_calib_10.py):

Ergebnis. Das Modell trifft die Größenordnung (kJ/mol) durchweg und die Richards-Regel für einfache Metalle exzellent (Kupfer ×1,00, Eisen ×1,28); die verflüssigten Gase liegen knapp darunter (Methan ×0,94, Stickstoff ×0,87, Argon ×0,69). Ionen-, H-Brücken- und Molekülkristalle (NaCl ×0,37, Wasser ×0,44, Benzol ×0,28, CO₂ ×0,23) haben deutlich höhere reale Schmelzentropien (ΔS_exp/R = 2,6–5,0), weil beim Schmelzen zusätzliche Ordnungsgrade frei werden (Ionenpaarung, Orientierung mehratomiger Moleküle), die das Ein-Ordnungsparameter-Modell nicht erfasst. Die Grenze ist damit klar benannt: das Idealbild gilt für isotrop, ungerichtet bindende Stoffe (kein Fit).

7.7 Vereinheitlichung der Energieskalen: ε ↔ ε_lg

Bislang trug der Schmelzzweig die Bindungsenergie ε (τ_m = 1,20), der l/g-Zweig eine separate Kohäsionsskala ε_lg/k_B = 3,375·T_c (Kap. 7.4). Beide lassen sich aus einem Modellparameter herleiten. Die van-der-Waals-Kohäsion a entsteht im Modell aus der Summe der Bindungsenergien je Molekül (ΒΙΝDUΝG): jedes Molekül trägt z_max/2 Bindungen der Energie ε. Ein Mean-field-Bragg-Williams-Gittergas mit Nachbaranziehung ε liefert die klassische Beziehung für die kritische Temperatur:

k_B·T_c = z_max·ε / 4, also τ_c ≡ k_B T_c/ε = z_max/4.

Setzt man dies in ε_lg = a/b = (27/8)·k_B·T_c ein (Kap. 7.4), folgt die Kopplung beider Skalen:

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

Damit liegen Schmelzpeak, Perkolation und kritischer Punkt auf einer ε-Skala (unify_scale_12.py):

Vorhersage τ_c/τ_m = z_max/(4·τ_m): 2D → 1,25, 3D → 2,50. Gegenüber realen Stoffen:

Ergebnis. Die 3D-Vorhersage τ_c/τ_m = 2,50 trifft die reale Kennzahl der einfachen Molekülstoffe (T_c/T_m ≈ 2,0–2,1) gut; die Metalle liegen höher (T_c grob geschätzt, andere Bindungsart). Die 2D-Vorhersage (1,25) ist zu niedrig — und τ_c(2D) = 1,50 liegt sogar unter τ_b = 1,85 (der kritische Punkt läge unter der Perkolationsschwelle). Das ist kein Rechenfehler, sondern die bekannte Sonderrolle von 2D: Kohäsions- und Bindungsskala ordnen sich dort nicht sauber, ein robustes Flüssigfenster fehlt (vgl. stetiges 2D-Schmelzen, Kap. 5.5). Erst in 3D (τ_m = 1,20 < τ_c = 3,00) trennen sich die Skalen physikalisch sinnvoll.

7.8 Parameterfreier Tripelpunkt aus z̄_c

Kap. 7.5 hatte v_solid = 0,40 noch gesetzt. Die frühere Fassung in der Entwicklung leitet die feste Packungsdichte aus der Perkolationsschwelle ab. Die Koordinationszahl eines Gitters bei reduziertem Volumen v_r folgt der Dichte-Koordinations-Relation (harte Wand v_r = 1/3 ↔ dichteste Packung, z̄ = z_max):

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

Der Festkörper verliert seine Steifigkeit, sobald z̄ die Schwelle z̄_c unterschreitet; die zugehörige Dichte ist v_solid. Der Tripelpunkt ist dann der Punkt, an dem die flüssige Koexistenzgrenze v_liq(τ_r) diese Dichte erreicht (unify_scale_12.py, 2D, z_max = 6):

Verbesserte Version. Die Schmelzdichte v_solid wird von der generischen Maxwell-Rigidität z̄_c = 2d gesetzt (2D → 4, 3D → 6), nicht von der Bond-Konnektivität (Begründung in Kap. 5.4a: das Schmelzen ist Verlust der Steifigkeit, nicht der Verbundenheit). Damit lautet der belastbare Tripelpunkt (rigidity_3d_13.py, 2D):

Ergebnis 2D. Selbst mit der (physikalisch korrekten) Rigiditäts-Schwelle liefert das 2D-Einzelgitter kein sauberes Flüssigfenster: v_solid = 0,50 ergibt zwar den realen Kennwert τ_triple/τ_c ≈ 0,77, doch in absoluten ε-Einheiten ist τ_triple = 0,77·1,50 = 1,16 und liegt damit unter τ_m = 1,20. Die Konnektivitäts-Schwelle (z̄_c = 2,08) drückt v_solid ≈ 0,96 ≈ v_c sogar aus der Koexistenzkuppel — Tripel- und kritischer Punkt fallen zusammen. Beide Rechnungen bestätigen die 2D-Sonderrolle (Kap. 5.5, 7.7): das 2D-Einzelgitter trägt kaum eine eigenständige flüssige Phase. Das ist ein ehrlicher Befund, keine Rechen- oder Modellschwäche im engeren Sinn.

Der entscheidende Test ist daher die 3D-Kontrollrechnung (Kap. 7.9): dort setzt dieselbe Rigiditätsregel (z̄_c = 6 → v_solid = 0,667) ein sauberes Fenster frei. In jedem Fall ist der Tripelpunkt nun parameterfrei an z̄_c gebunden — der freie Parameter v_solid ist eliminiert.

7.9 3D-Kontrollrechnung: sauberes Flüssigfenster auf einer ε-Skala

Prüfung, ob sich der 2D-Kollaps in 3D auflöst, wenn alle Marken auf einer ε-Skala liegen (τ_c = z_max/4 aus Kap. 7.7; v_solid aus der generischen Rigidität z̄_c = 2d; Tripelpunkt dort, wo die reduzierte van-der-Waals-Flüssigdichte v_liq(τ_r) die Feststoffdichte erreicht). Ergebnis (rigidity_3d_13.py):

Ergebnis. In 3D öffnet sich ein sauberes Flüssigfenster: τ_m = 1,20 < τ_triple = 2,83 < τ_c = 3,00 (Fensterbreiten: τ_triple − τ_m = 1,63; τ_c − τ_triple = 0,18). Der obere Rand ist eng — der Tripelpunkt liegt dicht unter dem kritischen —, aber die Reihenfolge stimmt und die flüssige Phase existiert als eigenständiges Feld. In 2D bleibt das Fenster geschlossen (τ_triple = 1,16 < τ_m). Damit ist die 2D-Sonderrolle zum dritten Mal unabhängig belegt und zugleich gezeigt, dass sie ein Dimensionseffekt ist: dieselbe Rigiditätsregel heilt das Modell in 3D. Der enge obere Rand wird in Kap. 7.10 durch einen echten 3D-Perkolations-Scan und eine Sensitivitätsanalyse erhärtet.

Modellgrenzen-Kasten — die 2D-Sonderrolle. Drei unabhängige Rechnungen zeigen denselben Befund: Das 2D-Einzelgitter trägt kaum eine eigenständige flüssige Phase. 1. τ_c < τ_b (Kap. 7.7): der kritische Punkt (1,50) läge unter der Konnektivitäts-Perkolation (1,85) — die Kohäsions- und Bindungsskala ordnen sich in 2D nicht sauber. 2. v_solid ≈ v_c (Kap. 7.8): die Konnektivitäts-Schmelzdichte (0,96) fällt fast mit der kritischen Dichte zusammen — Tripel- und kritischer Punkt kollabieren. 3. τ_triple < τ_m (Kap. 7.9): selbst mit der korrekten Maxwell-Rigidität liegt der Tripelpunkt (1,16) unter dem Schmelzpeak (1,20) — kein Fenster.

Das ist ein ehrlicher Befund, keine Rechen- oder Modellschwäche: Es ist der bekannte Dimensionseffekt (stetiges 2D-Schmelzen, KTHNY-artig; Kap. 5.5). In 3D heilt dieselbe Rigiditätsregel das Modell (Kap. 7.9/7.10). Ob eine zweite Nahordnung (Kap. 6.4) das 2D-Bild zusätzlich verbessern könnte, bleibt eine Erweiterungsfrage, kein Mangel des Idealbildes. Da das Modell der Veranschaulichug dienen soll, nehme ich Schwäche hin und verzichte auf einen axiomistische Erweiterung des Modells.

7.10 Echter 3D-Bond-Perkolations-Scan (FCC/BCC) und Fensterrand-Sensitivität

Kap. 7.8/7.9 übertrugen den 2D-Konnektivitätswert z̄_c = 2,08 gedanklich auf 3D (dort v_solid = 1,92, unphysikalisch) und stützten die 3D-Aussage auf die generische Rigidität z̄_c = 2d. Fassung 13 prüft die Konnektivität in 3D nun direkt: Bond-Perkolation auf FCC (z_max = 12) und BCC (z_max = 8) per Monte-Carlo mit Union-Find, periodische Ränder, deterministischer LCG-Seed (rigidity_3d_14.py). Schwelle p_c = Kreuzung P_∞(p) = 0,5 (dieselbe Methodik wie die 2D-Finite-Size- Studie in Kap. 5.6):

Die MC-Werte liegen wegen der moderaten Boxkante (L = 12, Pure-Python) leicht über den Literaturschwellen (endliche-Größe-Bias, wie in 2D von unten/oben je nach Schätzer); die Reihung und Größenordnung stimmen (FCC dichter → kleineres p_c als BCC).

Deutung (unabhängige 3D-Bestätigung). Die 3D-Konnektivitätsschwellen liegen bei z̄_c(3D) ≈ 1,5–1,6 — sogar noch tiefer als der 2D-Wert 2,08 und weit unter der Maxwell-Rigidität 2d = 6. Die daraus folgende „Konnektivitäts-Schmelzdichte" v_solid = 1/(3 p_c) ≈ 1,7 (BCC) bzw. 2,6 (FCC) liegt weit außerhalb der Koexistenzkuppel (v_liq,max ≈ 0,94). Damit ist — ohne den 2D-Wert auf 3D zu übertragen — direkt gezeigt: in 3D kann die Konnektivität die Schmelzdichte erst recht nicht setzen; sie muss von der Rigidität (z̄_c = 6 → v_solid = 0,667) kommen. Die zwei Schwellen sind in 3D noch klarer getrennt als in 2D. Das echte τ_b(3D) (Temperatur, bei der z̄(τ) die niedrige Konnektivitätsschwelle 1,5 erreicht) läge damit oberhalb des 2D-Werts τ_b, ganz nahe am kritischen Punkt — konsistent mit dem Bild „Konnektivitäts-Zerfall = Sieden".

Sensitivität des engen oberen Fensterrands (τ_c − τ_triple ≈ 0,18). Numerisch um v_solid = 0,667 (±2 %) ergibt sich dτ_triple/dv_solid ≈ 1,55. Das Fenster schlösse (τ_triple = τ_c = 3,00) erst, wenn v_solid um ≈ 0,11 stiege (auf ≈ 0,78), d. h. wenn die effektive 3D-Rigiditätsschwelle merklich über 2d = 6 läge (z̄_c ≳ 7,7). Für die generische Maxwell-Zählung (z̄_c = 2d = 6) bleibt die Reihenfolge τ_m < τ_triple < τ_c also robust; der obere Rand ist real eng, aber nicht gefährdet. Der diskrete Gitterwechsel fest(1)→fest(2) (FCC↔BCC) ist in Fassung 14 quantifiziert (Kap. 6.3), der echte 3D-τ-Scan der Bindungskinetik folgt in Kap. 7.11.

7.11 Echter 3D-τ-Scan der Bindungskinetik: τ_b(3D) und Schmelzlinie auf einer ε-Skala

Kap. 7.8–7.10 stützten die 3D-Aussage auf die generische Rigidität (z̄_c = 2d = 6) und den statischen Konnektivitäts-Scan (z̄_c(3D) ≈ 1,55). Offen war der eigentliche τ-Scan der Bindungskinetik in 3D: Wo auf der Temperaturachse erreicht das selbstkonsistente z̄(τ) diese niedrige Konnektivitätsschwelle (= echtes τ_b(3D)), und wo liegt der Schmelzpeak? Fassung 14 rechnet das selbstkonsistente z̄(τ) (Kap. 4.2, x = 1/(1 + G·e^{−(ε+K·x)/τ}]) direkt auf FCC (z_max = 12) und BCC (z_max = 8) mit G = z_max und derselben Kooperativität J = 0,55 wie im 2D-Referenz­ fall; n_blk(τ) stammt aus einem FCC/BCC-Bond-Monte-Carlo (Union-Find), und

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

c_v = dU/dτ (bondkinetics_3d_15.py, deterministischer LCG-Seed; n_blk geglättet, da nur der schwache Blockrotationsterm daran hängt). Drei Marken werden ausgelesen: der c_v-Peak τ_m (Schmelzen), die Rigiditätsmarke τ_rig (z̄(τ) kreuzt 2d = 6) und das echte τ_b(3D) (z̄(τ) kreuzt die Konnektivitätsschwelle z̄_c(3D) aus Kap. 7.10):

Ergebnis. Zwei Befunde, beide konsistent mit dem bisherigen Bild:

1. Schmelzen liegt tief und ist doppelt belegt. Der thermodynamische Schmelzpeak τ_m (FCC 1,70, BCC 1,40) fällt praktisch mit der Rigiditätsmarke τ_rig (z̄ = 6; FCC 1,73, BCC 1,35) zusammen — dieselbe Temperatur, an der das Netz seine Steifigkeit verliert. Das verankert die 3D-Schmelzlinie auf der ε-Skala: T_m = (ε/k_B)·τ_m, also z. B. für Argon (T_m = 83,8 K) ε/k_B = 83,8/1,70 ≈ 49 K. Weil K = J·z_max in 3D größer ist als in 2D (6,60 statt 3,30), liegt der 3D-Schmelzpeak höher als der 2D-Wert 1,20 — kein Widerspruch, sondern die stärkere Kooperativität der dichteren Packung.
2. Das echte τ_b(3D) sitzt am kritischen Punkt. Weil die 3D-Konnektivitätsschwelle so niedrig ist (≈ 1,55), erreicht z̄(τ) sie erst ganz oben: τ_b(3D) = 3,18 (FCC) bzw. 2,85 (BCC) — am bzw. leicht über τ_c = z_max/4 (3,00 bzw. 2,00). Der Bindungs-Riesencluster überlebt also die gesamte flüssige Phase und zerfällt erst am kritischen Punkt. Das ist exakt die Vorhersage aus Kap. 7.10 („das echte τ_b(3D) läge oberhalb des 2D-Werts, ganz nahe am kritischen Punkt") und bestätigt das Bild „Sieden = Konnektivitäts-Zerfall" unabhängig: in 3D bleibt die Materie bis dicht an den kritischen Punkt durchgehend gebunden.

Beide Befunde erhärten die Schwellentrennung aus Kap. 5.4a: Schmelzen ↔ Rigidität (tief, τ_m ≈ τ_rig), Sieden ↔ Konnektivität (hoch, τ_b(3D) ≈ τ_c). Die Schmelzdichte kommt daher aus der Rigidität (v_solid = 0,667), nicht aus der Konnektivität — die zwei Schwellen sind in 3D auf der Temperaturachse so weit auseinandergezogen wie nie zuvor im Modell.

7.12 Die fest(1)→fest(2)-Umwandlungslinie in der (τ, v)-Fläche

Kap. 6.3 lieferte die Gitterumwandlung FCC↔BCC als Punkt τ_ss (bei p = 0). Über den Druck wird daraus eine Linie (polymorph_line_16.py), die als vierte Grenzlinie in die gemeinsame (τ, v)-Fläche aus Kap. 7.5/7.9 gehört — neben Koexistenzkuppel (Sieden), Schmelzlinie und harter Wand.

Volumensprung. Die beiden Festkörper sitzen bei nahezu fester Dichte; ihr Volumenunterschied folgt aus den Packungsdichten (v ∝ 1/η): η_FCC = 0,7405, η_BCC = 0,6802. Mit v_FCC = 0,667 (Kap. 7.9) ergibt sich v_BCC = v_FCC·η_FCC/η_BCC = 0,726, also

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

Das offenere BCC hat das größere Volumen — Umkristallisieren FCC→BCC ist mit einer Volumenzunahme verbunden (dieselbe Aussage wie im Comic-Strip „Umkristallisieren").

Clausius–Clapeyron. Mit ΔS = S_BCC − S_FCC (Hoch- minus Tieftemperaturphase) lautet die Steigung

dp/dτ = ΔS/Δv.

Beide Größen sind positiv (BCC ist die Hochtemperaturphase, hat also mehr Entropie, und das größere Volumen), also dp/dτ > 0: höherer Druck hebt die Umwandlungstemperatur, weil er das dichtere FCC begünstigt. Für die Referenz θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, L = 0,667 ε):

Linie in (τ, p) für die Referenz (steigend): p = 0 → τ_ss = 1,33; p = 2 → 1,57; p = 4 → 1,81; p = 6 → 2,04 (reduzierte Einheiten).

In der (τ, v)-Fläche erscheint die Umwandlung als zwei nahezu senkrechte Festkörper-Dichtebänder — fest(1) = FCC bei v ≈ 0,667 (dichter, tiefe T) und fest(2) = BCC bei v ≈ 0,726 (offener, hohe T) —, die bei τ_ss ineinander umschlagen; darüber folgt bei τ_m das Schmelzen. Die vollständige Reihenfolge auf der τ-Achse ist damit τ_ss < τ_m < τ_c ≤ τ_b(3D) (Kap. 7.11).

Reale Einheiten. Über die Schmelzpunktkalibrierung ε/k_B = T_m/τ_m (τ_m = 1,70, 3D-FCC) werden c_v-Stufe und Umwandlungswärme greifbar. Die c_v-Stufe ist kalibrierungsunabhängig: |Δc_v| = 1,0 k_B = 1 R = 8,31 J/(mol·K) = 0,0083 kJ/(mol·K). Die Umwandlungswärme ΔH_ss = L·R·(ε/k_B) (Referenz L = 0,667 ε):

Die Umwandlungswärmen sind klein und positiv (Größenordnung 0,2–0,9 kJ/mol) und liegen damit deutlich unter den Schmelzenthalpien (Kap. 7.6, mehrere kJ/mol) — genau das erwartete Bild einer schwachen Polymorphie-Umwandlung mit kleiner c_v-Stufe. (Das θ-Verhältnis 0,60 liegt bewusst außerhalb des strengen Idealbildes; Kap. 6.6 zeigt, dass die reine Geometrie θ ≈ 0,82–0,94 gäbe und die Umwandlung dann geschlossen bliebe — die Linie ist die Wenn-dann-Konsequenz einer weichen BCC-Phase.)

7.13 fest(1)→fest(2)-Umwandlungslinie im (τ, p)-Diagramm, kalibriert an Eisen

Kap. 7.12 führt die Umwandlung als (τ, v)-Linie; Fassung 16 ergänzt sie als Linie im (τ, p)-Bild und kalibriert sie an einem realen FCC↔BCC-Polymorph (solidsolid_pt_17.py).

Druckfeld aus der Hartkern-van-der-Waals-EOS. An den beiden Festkörperdichten v_FCC = 0,667 und v_BCC = 0,726 (Kap. 7.12) liefert die reduzierte EOS p_r = 8τ_r/(3v_r − 1) − 3/v_r² (Kap. 7.3) durchweg positive Drücke; bei gleichem τ_r trägt die dichtere FCC den höheren Druck (τ_r = 1: p_r = 1,25 gegen 1,10). Kompression begünstigt also FCC — konsistent mit der positiven Clausius–Clapeyron-Steigung dp/dτ = ΔS/Δv > 0.

Linie und Umwandlungspunkt. Für die Referenz θ_BCC/θ_FCC = 0,60 (τ_ss = 1,333, ΔS = 0,50 k_B, dp\/dτ = 8,47 in reduzierten Einheiten p\ = p·D³/ε) steigt die Linie p\(τ) = 8,47·(τ − 1,333): p\ = 0 → τ = 1,33; p\ = 2 → 1,57; p\ = 4 → 1,81; p\ = 6 → 2,04. Bei **p\_tri ≈ 3,1 trifft die Solid-Solid-Linie die Schmelzmarke τ_m = 1,70 — ein fest(1)-fest(2)-flüssig-Umwandlungspunkt** in der (τ, p)-Ebene (Analogon zum Tripelpunkt).

Kalibrierung an Eisen (ε/k_B = T_m/τ_m = 1811 K / 1,70 = 1065 K). Eisen hat zwei FCC↔BCC-Übergänge mit entgegengesetztem Vorzeichen — genau der Testfall:

Für γ→δ (τ_ss = 1667/1065 = 1,565 → θ_BCC/θ_FCC ≈ 0,61) trifft das Modell Vorzeichen und Größenordnung: ΔH_ss ≈ 4,2 kJ/mol (real ≈ 0,84 kJ/mol — Größenordnung kJ/mol, der Betrag der scharfen Schwellennäherung ~5× zu hoch), c_v-Stufe 1 R = 8,31 J/(mol·K) parameterfrei, dT/dp ≈ +245 K/GPa. Der Übergang α→γ dagegen hat die dichte Phase (γ) als Hochtemperaturphase — das kann reine Packungsgeometrie nicht liefern (dichter = tiefer = tiefe T); α-Eisen wird magnetisch stabilisiert. Eisen bestätigt damit beide Seiten des ehrlichen Befundes: die geometrische Topologie (γ→δ) liegt im Idealbild, die magnetische Anomalie (α→γ) außerhalb — dasselbe „außerhalb", das Kap. 6.6/6.8 über θ ≈ 0,87…1,00 > 0,61 quantifizieren.

7.14 Redlich–Kwong als getrennter Idealbild-Zweig

Der Guldberg-Abstand des van-der-Waals-Zweigs (Kap. 7.4, T_b/T_c ≈ 0,48 statt ≈ 0,60) wurde dort nur benannt. Auf ausdrücklichen Wunsch führt eine verbesserte Fassung die Redlich–Kwong-Kohäsion als eigenständige Idealbild-Variante mit eigener Kalibrierung, eigener Koexistenz und eigenem Fehlerabschnitt (rk_branch_17.py) — der van-der-Waals-Zweig bleibt der Default; alle übrige Modellphysik (Schmelzen, Perkolation, fest(1)→fest(2)) ist von der Wahl der Kohäsions-EOS unberührt.

Eigene Energieskala. Aus a_RK = Ω_a·R²T_c^{2,5}/p_c und b_RK = Ω_b·R T_c/p_c folgt die Kohäsionsenergie des Zweigs zu

ε_lg,RK/k_B = (Ω_a/Ω_b)·T_c = 4,93·T_c (gegen van der Waals ε_lg/k_B = 27/8·T_c = 3,375·T_c),

größer, weil die temperaturabhängige Kohäsion a/√T den Faktor √T_c einfaltet. Der aus b_RK folgende Durchmesser trifft die realen Werte ähnlich gut wie im van-der-Waals-Zweig (Argon D ≈ 3,34·10⁻¹⁰ m).

Eigene Koexistenz (reduziert, kritischer Punkt (1, 1, 1), Z_c = 1/3; p_r = 3τ_r/(v_r − β) − A/(√τ_r·v_r(v_r + β)), β = 0,25992, A = 3,84732):

Eigener Fehlerabschnitt (Siedepunkt bei p_r,atm = p_atm/p_c, T_b/T_c = τ_r,b):

Ergebnis. Der RK-Zweig hebt T_b/T_c von van der Waals ≈ 0,48 auf ≈ 0,60 und trifft damit die Guldberg-Konstante und die realen Werte praktisch (mittlerer Fehler −1 % gegen van der Waals −21 %). Er steht als zweite, klar getrennte Idealbild-Variante neben dem van-der-Waals-Default; welcher Zweig „richtiger" ist, ist eine Frage der geforderten Genauigkeit, nicht der Modellstruktur — der Kern von Modell 3 (z̄ als Ordnungsparameter) bleibt in beiden Zweigen identisch.

7.15 fest(1)→fest(2) an weiteren Polymorphen: Titan und Edelgas-Kontrolle

Kap. 7.13 kalibriert die fest(1)→fest(2)-(τ, p)-Linie an Eisen. Die verbesserte Fassung sichert Vorzeichen und Größenordnung an einem zweiten, magnetisch unauffälligen Polymorph ab und ergänzt eine Negativkontrolle (solidsolid_titanium_18.py). Kalibrierung und Volumensprung sind wie in Kap. 7.13 (ε/k_B = T_m/τ_m, τ_m = 1,70; Δv = +8,9 % dicht→offen aus der Packungsgeometrie); hcp wird wie FCC über z = 12 genähert (der fcc/hcp-Unterschied wird nur benannt, nicht aufgelöst).

Titan hat mit α(hcp) → β(bcc) bei 1155 K genau die Modell-Topologie: die dichte Phase (hcp, z = 12) ist die Tieftemperaturphase, beim Erwärmen wird sie offen (bcc, z = 8), also ΔS > 0 und Δv > 0 → dp/dτ > 0 — dasselbe Vorzeichen wie Eisen γ→δ, aber ohne magnetische Anomalie. Zwei unabhängige dicht→offen-Polymorphe im Vergleich:

Beide treffen Vorzeichen (dp/dτ > 0, Druck begünstigt die dichte Phase) und Größenordnung (ΔH ~ kJ/mol, c_v-Stufe parameterfrei 1 R = 8,31 J/(mol·K), dT/dp ~ 10²·K/GPa). Der Betrag von ΔH_ss liegt bei beiden systematisch ~einen Faktor zu hoch — die scharfe Schwellennäherung des μ-Vergleichs überschätzt die latente Umwandlungswärme; das ist eine benannte Modellgrenze, kein Vorzeichenfehler.

Argon ist die Negativkontrolle: festes Argon ist FCC und bleibt FCC bis zum Schmelzen (83,8 K) — bei Normaldruck existiert kein stabiles bcc. Genau das sagt das reine Paarpotential-Idealbild voraus: θ_BCC/θ_FCC ≈ 0,9…1,0 > 0,61 und ein scherinstabiles LJ-BCC (Kap. 6.9) ⇒ kein fest(1)→fest(2)-Übergang. Die Gegenprobe bestätigt die Struktur des Befundes: das Idealbild erzeugt eine fest(1)→fest(2)-Linie nur dann, wenn nicht-geometrische Beiträge (elektronisch/magnetisch wie in Eisen und Titan) das BCC überhaupt stabilisieren. Titan liefert den positiven, Argon den negativen Beleg desselben ehrlichen Rahmens aus Kap. 6.8/6.9.

8 — Literaturvergleich: Modellvorstellungen der statischen Thermodynamik

Zusammenfassung. Kapitel 8 ordnet Modell 3 in die etablierten Modellvorstellungen der statischen Thermodynamik ein:

• 8.1 Übersichtstabelle: elf Modelle (ideales Gas, van der Waals, Redlich–Kwong, Virialentwicklung, Einstein, Debye, Lennard-Jones/harte Kugeln, Zellentheorie, Ising/Gittergas, Perkolation, Flory-Huggins) mit Phase(n) und statisch-thermodynamischer Grundidee.
• 8.2 Modelle im Einzelnen: je ein Absatz mit Merkmalen, Anwendung und Grenzen, beginnend mit dem Quelllink; frei zugängliche Quellen sind unter model3/links/ archiviert, die Zellentheorie (Paywall) nur zitiert.
• 8.3 Gegenüberstellung: Modell 3 ist ein kooperatives Zweizustands-Bindungsmodell mit Ordnungsparameter z̄; engster Verwandter ist Ising/Gittergas, das Sieden ist direkt aus der Perkolation übernommen, das Phononenbild (Einstein/Debye) trägt die Schmelz-Steifigkeit. Modell 3 verklammert diese drei Traditionslinien auf einer gemeinsamen ε-Skala.

Dieses Kapitel ordnet Modell 3 in die etablierten Modellbilder der Gleichgewichts- („statischen") Thermodynamik von Feststoffen, Flüssigkeiten und Gasen ein. Zuerst eine tabellarische Übersicht (8.1), dann je Modell ein kurzer Fließtextabsatz mit typischen Merkmalen, Anwendungsbereich und Grenzen — jeder Absatz beginnt mit dem Link zur Quelle (8.2). Abschließend eine stichwortartige Gegenüberstellung zu Modell 3 (8.3). Frei zugängliche Quellen sind unter model3/links/ archiviert; das Prüfskript model3/scripts/check_links_16.py hält Zitat und Archiv per SHA-256 konsistent.

8.1 Übersichtstabelle

8.2 Modelle im Einzelnen

https://de.wikipedia.org/wiki/Kinetische_Gastheorie — Ideales Gas / kinetische Gastheorie. Teilchen ohne Eigenvolumen und ohne Wechselwirkung außer elastischen Stößen; Druck entsteht rein aus dem Impulsübertrag. Liefert die Zustandsgleichung pV = Nk_BT, die Maxwell-Boltzmann-Geschwindigkeitsverteilung und ⟨E⟩ = (3/2)k_BT. Anwendung: verdünnte Gase bei hoher Temperatur und niedrigem Druck. Grenzen: keine Kondensation, keine flüssige oder feste Phase, kein Phasenübergang — Wechselwirkung ist per Konstruktion abwesend.

https://de.wikipedia.org/wiki/Van-der-Waals-Gleichung — Van-der-Waals-Gleichung. (p + a/v²)(v − b) = k_BT ergänzt das ideale Gas um ein Ausschlussvolumen b (Abstoßung) und einen Kohäsionsterm a (Anziehung). Das ist die erste Zustandsgleichung mit einem Gas-Flüssig-Übergang, kritischem Punkt und Maxwell-Konstruktion. Anwendung: qualitatives Bild der Kondensation, kritische Größen. Grenzen: nahe dem kritischen Punkt quantitativ ungenau (Mean-Field-Exponenten), keine feste Phase.

https://de.wikipedia.org/wiki/Redlich-Kwong-Gleichung — Redlich-Kwong-Zustandsgleichung. Kubische EOS mit temperaturabhängigem Attraktionsterm a/[√T · v(v + b)]. Die √T-Skalierung verbessert die Dampf-Flüssig-Gleichgewichte deutlich gegenüber van der Waals (vgl. Kap. 7.4, wo dies das Verhältnis T_b/T_c anhebt). Anwendung: technische Gasphasen- und Verdampfungsrechnungen. Grenzen: empirisch parametrisiert, weiterhin keine feste Phase; nahe dem kritischen Punkt ebenfalls Mean-Field-artig.

https://de.wikipedia.org/wiki/Virialentwicklung — Virialentwicklung. Systematische Entwicklung pV/Nk_BT = 1 + B(T)/v + C(T)/v² + … Die Virialkoeffizienten folgen aus dem Intermolekularpotential (Cluster-Integrale der statistischen Mechanik) und sind damit mikroskopisch begründet. Anwendung: Gase bis zu mittlerer Dichte, Referenz für andere EOS. Grenzen: die Reihe konvergiert für dichte Flüssigkeiten/Feststoffe nicht; höhere Koeffizienten sind rechnerisch aufwendig.

https://de.wikipedia.org/wiki/Einstein-Modell — Einstein-Festkörper. Der Kristall ist ein Ensemble von 3N unabhängigen harmonischen Oszillatoren gleicher Frequenz ω. Erklärt den Abfall der Wärmekapazität c_v unterhalb der charakteristischen Temperatur θ_E = ħω/k_B und den Dulong-Petit-Grenzwert 3Nk_B bei hoher Temperatur. Anwendung: einfachste Beschreibung fester Wärmekapazitäten. Grenzen: falsches Tieftemperaturgesetz (exponentiell statt T³), keine akustischen Moden, nur die feste Phase.

https://de.wikipedia.org/wiki/Debye-Modell — Debye-Modell. Behandelt den Festkörper als elastisches Kontinuum mit einem Phononenspektrum bis zur Debye-Frequenz; c_v ∝ T³ bei tiefer Temperatur, Dulong-Petit bei hoher — die Debye-Temperatur θ_D fasst die Gitterhärte zusammen. Anwendung: quantitativ gute Wärmekapazitäten einfacher Kristalle über den ganzen T-Bereich. Grenzen: Kontinuumsnäherung (ignoriert Dispersion realer Gitter), nur die feste Phase, kein Schmelzen.

https://de.wikipedia.org/wiki/Lennard-Jones-Potential — Lennard-Jones-Potential / harte Kugeln. Paarpotential 4ε[(σ/r)¹² − (σ/r)⁶] mit Abstoßungskern und schwacher Anziehung; das Grenzmodell harter Kugeln behält nur die Abstoßung. Es ist die Grundlage von Molekulardynamik und Monte-Carlo und beschreibt — numerisch — alle drei Phasen samt Tripelpunkt. Anwendung: Edelgase, generisches Referenzsystem. Grenzen: paarweise Additivität, im Kern nur per Simulation oder Störungstheorie auswertbar (keine geschlossene EOS).

https://doi.org/10.1098/rspa.1937.0210 — Zellentheorie der Flüssigkeit (Lennard-Jones-Devonshire). Jedes Molekül bewegt sich in einer Zelle (Käfig), die seine Nachbarn im gemittelten Feld bilden; die Flüssigkeit wird durch die Zell-Zustandssumme beschrieben und liefert eine Flüssigkeits-EOS sowie einen Schmelzansatz. Anwendung: einfache Flüssigkeiten, Bindeglied zwischen Fest und Gas. Grenzen: überschätzt die Ordnung, vernachlässigt Leerstellen/Löcher (dies beheben spätere Loch-Theorien). (Paywall — nur zitiert, nicht lokal archiviert; Erreichbarkeit im Prüfskript getestet.)

https://de.wikipedia.org/wiki/Ising-Modell — Ising-Modell / Gittergas. Auf einem Gitter tragen die Plätze eine Zweizustandsvariable (Spin ↑/↓ bzw. besetzt/leer) mit Nächste-Nachbar-Kopplung. Das Gittergas ist mathematisch isomorph zum Ising-Modell und beschreibt die Kondensation Gas↔Flüssig als kooperativen Ordnungs-Unordnungs-Übergang mit kritischem Punkt. Anwendung: kritische Phänomene, kooperative Übergänge. Grenzen: Gitterdiskretisierung, keine explizite kinetische Energie; in Mean-Field-Näherung nur qualitativ.

https://de.wikipedia.org/wiki/Perkolationstheorie — Perkolationstheorie. Rein geometrisches Modell: Plätze oder Bindungen sind mit Wahrscheinlichkeit p zufällig besetzt; oberhalb einer Schwelle p_c existiert ein durchspannender Cluster. Beschreibt Gelierung, Leitfähigkeitseinsatz und allgemein Konnektivitätsübergänge. Anwendung: Netzwerkbildung, Zerfall/Aufbau von Zusammenhang. Grenzen: ohne intrinsische Energie/Temperatur — die Kopplung an Thermodynamik muss von außen zugefügt werden.

https://de.wikipedia.org/wiki/Flory-Huggins-Theorie — Flory-Huggins-Theorie. Gittermodell für Lösungen und Polymere: Mischentropie der Kettenplatzierung plus ein Wechselwirkungsparameter χ steuern Mischbarkeit und Entmischung. Anwendung: Polymerlösungen, Phasendiagramme von Mischungen. Grenzen: Mean-Field, inkompressibles Gitter, χ meist empirisch — für einkomponentige Fest/Flüssig/Gas-Bilder nur begrenzt einschlägig, hier als Gitter-Referenz geführt.

8.3 Ähnlichkeiten und Unterschiede zu Modell 3

Modell 3 ist ein kooperatives Zweizustands-Bindungsmodell: Ordnungsparameter ist die mittlere Koordinationszahl z̄ = z_max·x, selbstkonsistent über die Kooperativität K = J·z_max mit der Bruchzustands-Entartung G = z_max; Schmelzen als Steifigkeits-(Phononen-)Schwelle, Sieden als Konnektivitäts-(Perkolations-)Zerfall am kritischen Punkt, Übergang erster Ordnung sobald ln G > 2 + 4ε/K.

Kurz: Modell 3 verklammert drei sonst getrennte Traditionslinien — das Gittergas/Ising-Bild (kooperativer Zweizustands-Übergang), das Perkolationsbild (Konnektivität → Sieden) und das Phononenbild Einstein/Debye (Steifigkeit → Schmelzen) — und kalibriert sie über einen einzigen Energieparameter ε auf einer gemeinsamen τ-Skala. Van der Waals und die Zellentheorie liefern die kontinuierliche Gegenprobe (EOS bzw. Flüssigkeitskäfig), gegen die sich das diskrete Bindungsbild vergleichen lässt.

9 — Anhang

Zusammenfassung. Der Anhang bündelt die formalen Nachweise:

• 9.1 Literatur: nach Autorennamen sortiert, mit Links zu den frei archivierten Quellen (SHA-256-Abgleich per check_links_16.py; die Zellentheorie paywallbedingt nur zitiert).
• 9.2 KI-Unterstützung: der von der KI geschätzte Innovationsanteil KI/Autor je wesentlicher Änderung (mit Aufschlüsselung des Autor-Anteils nach Strategie/Struktur/Redaktion und Gegenrechnung) sowie eine Aufstellung der Einzelschritte mit Kosten (geschätzt), Zeitdauer und Abstand zum vorigen Log (Summe ≈ 160 USD, klar gekennzeichnete Schätzungen).
• 9.3 Abkürzungs- und Symbolglossar: die häufig genutzten Symbole (τ, ε, z̄, K, G, θ, …), die hierdurch ausnahmsweise ohne vorangestellte Beschreibung im Fließtext stehen dürfen.
• 9.4 Berechnungsskripte und Testabdeckung: alle Zahlen stammen aus reproduzierbaren Python-Skripten (KI-Programmieranteil 100 %, Testabdeckung durchweg 0 %), samt der offenen „peinlichen Erklärung", dass der Autor Code und Ergebnisse nicht selbst geprüft hat.

9.1 Literatur

Sortiert nach Autorennamen. Frei zugängliche Quellen sind unter model3/links/ archiviert; ihre Übereinstimmung mit den zitierten Links prüft model3/scripts/check_links_16.py (SHA-256). Die Zellentheorie von Lennard-Jones & Devonshire ist paywallbedingt nur zitiert (Erreichbarkeit im Prüfskript getestet).

• Broadbent, S. R.; Hammersley, J. M. (1957): Percolation processes I. Crystals and mazes. Proc. Cambridge Phil. Soc. 53, 629–641. — Überblick/Archiv: Perkolationstheorie
• Clausius, R.; Maxwell, J. C.; Boltzmann, L. (1857–1872): Grundlegung der kinetischen Gastheorie. — Überblick/Archiv: Kinetische Gastheorie
• Debye, P. (1912): Zur Theorie der spezifischen Wärmen. Ann. Phys. 39, 789–839. — Überblick/Archiv: Debye-Modell
• Einstein, A. (1907): Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme. Ann. Phys. 22, 180–190. — Überblick/Archiv: Einstein-Modell
• Flory, P. J.; Huggins, M. L. (1941–1942): Gittertheorie der Polymerlösungen (χ-Parameter). — Überblick/Archiv: Flory-Huggins-Theorie
• Ising, E.; Lenz, W. (1920–1925): Beitrag zur Theorie des Ferromagnetismus. Z. Phys. 31, 253–258 (Gittergas mathematisch isomorph). — Überblick/Archiv: Ising-Modell
• Kamerlingh Onnes, H. (1901): Virialform der Zustandsgleichung. — Überblick/Archiv: Virialentwicklung
• Lennard-Jones, J. E. (1924): On the determination of molecular fields. Proc. R. Soc. Lond. A 106, 463–477. — Überblick/Archiv: Lennard-Jones-Potential
• 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 (Paywall — nur zitiert, nicht lokal archiviert.)
• Redlich, O.; Kwong, J. N. S. (1949): On the thermodynamics of solutions. V. An equation of state. Chem. Rev. 44, 233–244. — Überblick/Archiv: Redlich-Kwong-Gleichung
• van der Waals, J. D. (1873): Over de Continuïteit van den Gas- en Vloeistoftoestand (Dissertation, Leiden). — Überblick/Archiv: Van-der-Waals-Gleichung

9.2 KI-Unterstützung

Diese Modellentwicklung entstand in einem Mensch-KI-Zusammenspiel: der Autor (KI-Nutzer) legt fest, woran gearbeitet wird (Leitfragen, Entscheidungen, Modellwünsche); die KI (Claude Code) legt weitgehend fest, wie — Herleitungen, Skripte, Struktur und Text. Die beiden folgenden Tabellen schätzt die KI selbst.

(a) Wesentliche inhaltliche Änderungen und geschätzter Innovationsanteil. Der Innovationsanteil gewichtet Leitideen 45 %, strukturierende Arbeiten 40 %, redaktionelle Arbeiten 15 %. Leitideen (die woran-Ebene) stammen überwiegend vom Autor, strukturierende und redaktionelle Arbeiten (die wie-Ebene) überwiegend von der KI; wo die KI auch die physikalische Leitidee einbrachte (z. B. Perkolationsdeutung des Siedens, θ-Federmodell), steigt ihr Anteil. Werte in Prozent, KI + Autor = 100. Die drei zusätzlichen Spalten schlüsseln den Autor-Anteil nach den drei Wirkebenen auf — Strategie (Leitideen), Struktur (Herleitungen/Skripte/Aufbau) und Redaktion (Textarbeit). Ihr mit 45 % / 40 % / 15 % gewichteter Mittelwert reproduziert die Spalte Autor % als Gegenrechnung (model3/scripts/author_share_19.py, größte Abweichung 0,95 Prozentpunkte).

(b) Einzelschritte: Kosten (geschätzt), Zeitdauer, Abstand zum vorigen Log. Erstellungsdatum ist die Erstellungszeit der Log-Datei (Datei-mtime). Δ zum Vorlog ist der zeitliche Abstand zum vorigen Log-Eintrag (Maß für den Nutzeraufwand zwischen den Schritten, inkl. Pausen//compact). Kosten (geschätzt) sind die je Log dokumentierten Schätzungen (Anthropic-Listenpreise, ± 30 %). Zeitdauer (aktiv) ist eine KI-Schätzung der reinen Bearbeitungszeit; bei durch /compact getrennten Schritten kann sie größer sein als Δ. Ein Stern markiert einen vermuteten Nachbau/Überarbeitung vorhandener Inhalte.

Summe der dokumentierten Kosten (geschätzt, Schritte 01–21): ≈ 129 USD (Summe der Einzelschätzungen); mit Schritt 22 (16d, ≈ 10 USD) ≈ 139 USD, Schritt 23 (F16, ≈ 12 USD) ≈ 151 USD und Schritt 24 (F17, ≈ 9 USD) ≈ 160 USD. Alle Werte sind klar gekennzeichnete Schätzungen (keine ccusage-Messung verfügbar); die Log-Einzelheiten stehen in model3/logs/dev_log_NN.md.

9.3 Abkürzungs- und Symbolglossar

Grundsätzlich stellt dieser Text jedem Symbol einen beschreibenden Begriff voran (statt „c_v" also „Wärmekapazität c_v"). Die folgenden häufig genutzten Größen bilden davon die Ausnahme: Sie dürfen im Fließtext ohne vorangestellte Beschreibung erscheinen, weil sie hier ein für alle Mal erklärt sind. Alle Temperaturen τ… sind reduziert (in Einheiten der Bindungsenergie ε).

9.4 Berechnungsskripte und Testabdeckung

Alle im Text genannten Zahlen und Tabellen stammen aus reproduzierbaren Python-Skripten in model3/scripts/ (jeweils mit einer Rohausgabe …_output.txt). Der KI-Programmieranteil gibt an, welcher Anteil des Skriptcodes von der KI geschrieben wurde; die Testabdeckung meint automatisierte Tests (Unit-Tests, assert-Prüfungen). Sie ist durchweg niedrig: Die Skripte sind einmalige Rechenläufe ohne Testgerüst; ihre einzige Absicherung ist die Plausibilität der Rohausgabe.

Der Autor gesteht offen (peinliche Erklärung): Er hat weder den Quellcode dieser Skripte noch die damit erzeugten Zahlen, Tabellen und Kurven selbst geprüft, nachgerechnet oder fachlich verifiziert. Er hat der Kompetenz der KI hier blind vertraut, weil im Vordergrund allein der qualitative Charakter und die Plausibilität der Modellvorstellung standen. Die Werte sind entsprechend mit Vorsicht zu genießen (siehe auch den ausführlichen Hinweis in Kap. 10.4).

10 Selbstkritik

Gefühlt ist in der Tabelle im Kapitel 9.2.(a) der Arbeitsanteil des Autorens, geschätzt durch die KI, in alles drei Bereichen fast immer um mindesten 10-20%-Punkte zu hoch. Gefühlt habe ich nur die Grundidee, ein zwei Interpreationen und wichtige Verfeinerungen bei der graphischen Visualisierung beigesteuert.

Die (mathematischen) Kapitel 4-9 hat die KI vorgeschlagen, recherchiert und umgesetzt. Schlimmer noch: während der Entwicklung hatte ich das Gefühl, dass die KI auf die Entwicklung einer Statistik basierend auf der Koordinationszahl fixiert war. Meine eigentlichen Grundideen, dass es nicht nur Atome sondern auch Moleküle Bindungen ausbilden können und dass die Bindungsbrüche innerhalb einer Flüssigkeit fraktalen Gesetzmäßigkeiten folgen könnten, konnte die KI nur schwer annehmen, indem sie die Begrifflichkeiten zu ignorieren versuchte. Diese Grundideen werden im Text deshalb immer noch sehr nachrangig behandelt. Ich habe es nachträglich nicht umgeschrieben, weil man beim Nachdenken über die Phasenumwandlungen ausgehend von den Grundaxiomen diese Ideen natürlich mitbedenkt. Außerdem ist die Fixierung/Reduzierung der Ideen auf die Koordinationszahl als Maßzahl nicht der schlechteste Weg, um meine bildhafte Vorstellung zu mathematisieren, weil es sich für die Plausibilitätsprüfung im Nachhinein als recht tragend herausstellte, wobei ich aber nur eine sehr grobe Plausibilitätsprüfung der KI-Ergebnisse gemacht habe und wobei ich keine Prüfung des Python-Codes vorgenommen habe. Zu meiner Entschuldigung sei gesagt: Mir ging es immer nur um eine bildhafte Vorstellung der molekularen Verh#ältnisse in den verschiedenen Aggregatzuständen. Zu der Schmeichelei durch die KI, die(?) mathematisch die Vorstellung bestätigen wollte, sagte ich nicht nein. Die Manipulation der KI war bei mir also erfolgreich, weil ich die genannten Kapitel nur oberflächlichst geprüft habe und glaubend übernommen habe. (Diese Selbsterfahrung stärkt mein Mißtrauen, dass zur Kontrolle von KI ein "Human in the Loop"§ wirklich wirksam ist.)"

Dieser Text wurde KI-gestützt erstellt (Claude Code 2.1.161, Modell Opus 4.8 / claude-opus-4-8[1m], 1M-Kontext; Datum 2026-09-05). — KI-gestützt erstellt.

Wasserzeichen-Schlüssel (griechische Homoglyphen → lateinisch): Α→A, Β→B, Ε→E, Η→H, Ι→I, Κ→K, Μ→M, Ν→N, Ο→O, Ρ→P, Τ→T, Υ→Y, Χ→X, Ζ→Z. Markierte Schlagwörter: FRΕΙFLUG (FREIFLUG), VΙΕLΗΕΙΤ (VIELHEIT), ΒΙΝDUΝG (BINDUNG), DRΕΙΚLΑΝG (DREIKLANG), ΗΑRΤΚΕRΝ (HARTKERN), ΕRΗΑLΤUΝG (ERHALTUNG), RΙSS (RISS), FLΙΝΚ (FLINK).