Idealer-Stoff-Netzwerkmodell + Phononen v4 — Phasen laden + Modell/Empirie-Vergleich

Was ist das Idealer-Stoff-Netzwerkmodell?

Klassische Lehrbücher beschreiben Gase durch freie Teilchen, Festkörper durch Gitterplätze und Flüssigkeiten meist als unbefriedigenden Sonderfall dazwischen. Das Idealer-Stoff-Netzwerkmodell wählt einen einheitlichen Ansatz: Es betrachtet nicht einzelne Teilchen, sondern die statistischen Eigenschaften der Nachbarschaftsbeziehungen — das Netzwerk der nächsten Nachbarn.

Kernidee

Der Ordnungsparameter ist die mittlere Koordinationszahl z = durchschnittliche Anzahl nächster Nachbarn pro Teilchen.

Modellhierarchie

Level 0: Einzelteilchen (kein φ)
Level 1: Paarpotential (LJ)
Level 2: Netzwerkstatistik ← hier
Level 3: Feldtheorie (DFT, MCT)

Ziel

Anschauliches, physikalisch konsistentes Modell, das Wärmekapazität, Diffusion, Viskosität und Phasenübergänge aus einer Zustandsvariablen erklärt.

Die drei idealen Phasen

☁ Ideales Gas

Koordinationszahl	z ≈ 0
Netzwerk	nicht perkolierend
Orientierungsordnung Q₆	≈ 0
Bindungslebensdauer τ_b	sehr kurz → 0
C_V	(3 + f/2) R
Diffusion	frei, D ~ T^{3/2}/p
Viskosität η	steigt mit T^{1/2}

Exakter Grenzfall: ε → 0 oder ρ → 0. Alle Netzwerkterme verschwinden.

Phononen: keine stabilen Moden (z < z_c) → C_V^akust = (3/2)·R (rein kinetisch). θ_F → ∞, alle TA unterdrückt.

☬ Ideale Flüssigkeit

Koordinationszahl	z_c < z < z_krist
Netzwerk	perkolierend, ungeordnet
Orientierungsordnung Q₆	≈ 0
Bindungslebensdauer τ_b	endlich (ps–ns)
C_V	(3+f/2)R + Netzwerkterm
Diffusion	Sprünge, D ~ a²Γ
Viskosität η	fällt mit T (Arrhenius)

Mean-Field-Näherung — kein exakter Grenzfall wie Gas oder Kristall. Fehler ~30 % typisch.

Phononen (Frenkel): LA-Zweig voll aktiv. TA-Zweige nur für ω > ω_F = 1/τ_b. θ_F = 7.64/τ_b [K]. C_V^akust = C_V^LA + C_V^TA(Frenkel).

⧧ Idealer Kristall

Koordinationszahl	z = z_krist (fest)
Netzwerk	periodisch, geordnet
Orientierungsordnung Q₆	> 0.45 (FCC/HCP)
Bindungslebensdauer τ_b	→ ∞
C_V	3R (Dulong-Petit, exakt)
Diffusion	D ≈ 0
Viskosität η	→ ∞ (starr)

Exakter Grenzfall: dz/dT = 0 → C_V = 3R folgt zwingend.

Phononen (Debye, voll): 3 Zweige (1 LA + 2 TA), keine Frenkel-Grenze. C_V^akust = 9R·(T/θ_D)³·∫ … → 3R bei T ≫ θ_D (Dulong-Petit). T³-Gesetz: C_V ∝ T³ bei T ≪ θ_D (Quantenunterdrückung).

Kernaussagen — was das Modell leistet

Aussage

Begründung im Modell

Dulong-Petit gilt exakt für den idealen Kristall.

z = const → dz/dT = 0 → C_V = 3R.

Wasser hat anomal hohe Wärmekapazität.

Großes |dz/dT| · ε → großer Netzwerkterm in C_V.

Gas-Viskosität steigt mit T, Flüssig-Viskosität sinkt.

Gas: kinetische Gastheorie. Flüssigkeit: Arrhenius-τ_b.

Tripelpunkt ist eindeutig (ein Punkt im p-T-Diagramm).

Gibbs-Phasenregel: F = 1 − 3 + 2 = 0 Freiheitsgrade.

Flüssig-Kristall-Übergang ist stets 1. Ordnung (latente Wärme).

Kubischer Term c · Q₆³ im Landau-Funktional erzwungen durch Symmetrie.

ΔH_verd ≈ z_liq · ε / 2 (Trouton-konsistent).

Vollständiges Aufbrechen des Netzwerks beim Verdampfen.

ΔH_schmelz ≈ 0.10 · z_krist · ε (klein).

Nur Orientierungsordnung Q₆ bricht zusammen; z bleibt fast gleich.

Brownschen Bewegung = Projektion des zeitabhängigen Netzwerkgraphen.

D ~ a² · Γ / z, wobei Γ = Rekonfigurationsrate des Netzwerks.

Vergleich aller 8 Referenzstoffe — Modellcharakter^{(KI-generiert, Halluzinationen und Fehler möglich)}

Stoff	Bindungstyp	ε [kJ/mol]	z_liq	\|dz/dT\| [K⁻¹]	ΔH_verd: mod/exp	C_p: mod/exp	Modell-Eignung
Argon	Van der Waals (LJ)	1.0	12	0.050	0.92 ✓	1.47	Sehr gut
Quecksilber	Metallisch	5.0	8	0.001	0.34	0.98 ✓	C_p gut, T_siede schlecht
Cyclohexan	Van der Waals (ring)	6.6	10	0.003	1.00 ✓	~0.44	ΔH_verd exzellent
Benzol	VdW + π-Stapelung	6.0	10	0.003	0.98 ✓	~0.47	ΔH_verd gut
Propan	Van der Waals (kette)	2.5	8	0.002	0.53	0.75	ε zu klein
Ethanol	H-Brücken + VdW	4.5	6	0.006	0.35	0.59	Vibr. fehlen
Wasser	H-Brückennetzwerk	18.0	3.6	0.007	0.80 ✓	1.40	Gut kalibriert
DMSO	Dipol-Dipol (stark)	7.5	6	0.004	0.43	—	Modellgrenze

Legende: mod/exp = Modellwert / Experimentalwert (Faktor). Wert 1.0 = perfekt. Grün ≤ 1.25, Gelb ≤ 1.8, Rot > 1.8 (oder < 0.55).

Brownsche Bewegung im Netzwerkmodell

Die Brownsche Bewegung ist die sichtbare Projektion der thermischen Netzwerkrekonfiguration in den Ortsraum: Jedes Brechen und Neubilden einer Bindung überträgt einen Impuls auf das betrachtete Teilchen. Über viele solche Ereignisse entsteht ein Zufallsspaziergang (Random Walk).

Netzwerk-Interpretation

Mittlere quadratische Verschiebung: <r²(t)> = 6·D·t (t >> τ_b, 3D) Selbstdiffusionskoeffizient (Sprungmodell): D = a² · Γ / z_liq Γ = (1/τ_b) · exp(−z·ε/(6RT)) Stokes-Einstein (Kontrollformel): D = k_B·T / (6π·η·r_H) r_H ≈ a/2

Der Rechner zeigt D_Modell (Sprungmodell) und D_SE (Stokes-Einstein aus berechneter η) zum Vergleich.

Empirische D-Werte bei 25 °C ^{(KI-generiert, Halluzinationen und Fehler möglich)}

Stoff	D [10⁻⁹ m²/s]	Quelle
Quecksilber	1.71	NMR
Argon (90 K)	2.30	MD-Sim.
Propan	3.65	MD-Sim.
Ethanol	1.00	NMR-PFG
Wasser	2.30	NMR-PFG
Cyclohexan	1.44	NMR-PFG
Benzol	2.20	NMR-PFG
DMSO	0.71	NMR-PFG

Holz et al., PCCP 2 (2000) 4740. doi:10.1039/b005319h

Kritischer Punkt — Berechnung und Bedeutung

Gas-Flüssigkeit-Kritikpunkt

Endpunkt der Koexistenzlinie Gas ↔ Flüssigkeit. Oberhalb von T_c ist keine Unterscheidung mehr möglich (superkritisches Fluid).

Bedingung: ρ_Gas(T_c, p_c) = ρ_Flüssig(T_c, p_c) = ρ_c

Guldberg-Regel (±10 %): T_c ≈ 1.6 · T_siede Kritische Dichte (MF-Näherung): ρ_c ≈ ρ_liq / 2.7 (≈ ρ_liq / 3 vdW) Kein zuverlässiger Modellwert für p_c.

Fest-Gas: kein kritischer Punkt

Die Sublimationskurve endet am Tripelpunkt — nicht an einem kritischen Punkt. Der Fest-Flüssig-Übergang ist stets 1. Ordnung (Q₆ springt diskontinuierlich). Eine kontinuierliche Verbindung Fest ↔ Gas existiert nicht.

Warum? Gas und Flüssigkeit haben dieselbe Symmetrie (Q₆ = 0), Kristall bricht die Translationssymmetrie (Q₆ > 0). Symmetriebrechung kann nicht kontinuierlich verschwinden.

Klarstellung: Was berechnet ΔH_verd im Modell?^{(KI-generiert, Halluzinationen und Fehler möglich)}
ΔH_verd = z_liq · ε / 2 → FLÜSSIG → GAS (Verdampfung bei T_siede)
ΔH_sub = ΔH_schmelz + ΔH_verd → FEST → GAS (Sublimation, unterhalb Tripelpunkt)
Beispiel Argon: ΔH_sub(Mod) = 1.18 + 6.0 = 7.18 kJ/mol vs. exp. 7.74 kJ/mol (Abweichung 7 %).

Stoff	T_c,exp [K]	T_c,Guldberg [K]	Faktor	p_c [bar]	ρ_c [kg/m³]	ρ_liq/ρ_c
Argon	150.87	139.7	0.93	49.0	536	2.6
Propan	369.9	369.7	1.00	42.5	220	2.2
Benzol	562.2	565.2	1.01	49.0	305	2.9
Cyclohexan	553.6	566.2	1.02	40.7	274	2.8
Wasser	647.1	597.0	0.92	220.6	322	3.1
Quecksilber	1750	1007.8	0.58	1670	5800	2.3

Guldberg-Regel versagt für Quecksilber (metallische Bindung, andere T_siede-Logik). Für organische und einfache Stoffe Fehler < 10 %. DMSO: T_c ≈ 729 K (Schätzung).

Systemat. Stärken und Grenzen des Modells

✓ Stärken

Dulong-Petit exakt reproduziert
Tripelpunkt zwingend eindeutig (Gibbs)
Phasenübergang Flüssig-Kristall stets 1. Ordnung
Arrhenius-Viskosität ohne VFT-Divergenz
ΔH_verd für VdW-Stoffe (Ar, Cyclohexan, Benzol) gut
Wasser-Anomalie qualitativ erklärt
Einheitlicher Rahmen: Gas, Flüssigkeit, Kristall

✗ Grenzen

Ideale Flüssigkeit: kein exakter Grenzfall
C_p für große Moleküle: Vibrationsmoden fehlen
ΔH_melt: Cyclohexan-Plastikkristall-Effekt nicht erfasst
T_siede: nur bei ε-Kalibrierung aus ΔH_vap exakt
DMSO/polare Stoffe: hohe Siedepunkt durch Dipol unkalibrierbar
π-Stapelung (Benzol): ΔH_melt unterschätzt
Keine Quanteneffekte (H₂, He, leichte Atome)

♫ Phonon-Kernaussagen (Ergänzung durch Debye + Frenkel)

Aussage

Begründung im Modell

Dulong-Petit ist der klassische Grenzfall des Debye-Modells.

C_V^Debye(T ≫ θ_D) → 3R; alle Moden klassisch besetzt.

Gas: C_V^akust = (3/2)·R (nicht 3R wie bisher).

Kein perkolierendes Netzwerk → keine Phononen; nur kinetische Energie.

Flüssigkeit hat weniger C_V als Kristall bei gleichem T.

Frenkel: TA-Moden (ω < 1/τ_b) diffusiv → kein Beitrag zu C_V.

Festkörper: C_V ∝ T³ bei T ≪ θ_D (Debye-T³-Gesetz).

Quantenunterdrückung: Moden mit ℏω ≫ k_BT sind nicht besetzt.

Phononen sind Bosonen: ihre Anzahl ist nicht erhalten (μ = 0).

Z_Phonon = ∏_k 1/(1 − e^{−ℏω_k/kT}) ≠ Z_mol (Maxwell-Boltzmann).

DMSO als Lehrbeispiel für Modellgrenzen: DMSO (Dimethylsulfoxid) hat einen ungewöhnlich hohen Siedepunkt (189 °C) trotz moderater Molekülmasse (78 g/mol). Ursache: sehr starkes Dipolmoment (3.96 D) führt zu geordneten Dipol-Dipol-Wechselwirkungen, die in einem einfachen isotropen ε nicht erfasst werden. Das Modell unterschätzt ΔH_verd um Faktor 2.4. Korrekte Behandlung erfordert richtungsabhängige Wechselwirkungen ε(r, Ω) — ein Niveau-3-Modell.

Benzol ΔH_melt als Lehrbeispiel: Benzols Schmelzenthalpie (9.87 kJ/mol) ist fast 4× größer als die von Cyclohexan (2.63 kJ/mol), obwohl beide ähnliche ΔH_verd haben. Ursache: Im Benzol-Kristall stabilisiert π-π-Stapelung die Kristallstruktur zusätzlich (gerichtete Wechselwirkung der aromatischen Ringe). Das Modell mit isotropem ε kennt nur δ = 0.10 und unterschätzt deshalb ΔH_melt(Benzol) deutlich (Faktor 1.6).

Vergleich mit etablierten physikalischen Modellen^{(KI-generiert, Halluzinationen und Fehler möglich)}

Die folgende Tabelle stellt rund 20 etablierte Theorien dem Idealer-Stoff-Netzwerkmodell gegenüber. Bewertet werden: anwendbare Phasen, mathematischer Aufwand, Schlüsselgröße, Stärken, Schwächen und die Beziehung zum Idealer-Stoff-Modell.

Bezug zum Idealer-Stoff-Modell: = Grenzfall Idealer-Stoff-Modell reproduziert dieses Modell exakt ≈ analog strukturell ähnlich → Grundlage dieses Modell fundiert Idealer-Stoff-Modell ← Vereinfachung Idealer-Stoff-Modell ist vereinfachte Version ⊕ ergänzend verschiedene Aspekte ★ Idealer-Stoff dieses Modell

Aufwand: ★ sehr gering ★★ gering ★★★ mittel ★★★★ hoch ★★★★★ sehr hoch

KI-generiert, Halluzinationen und Fehler möglich)
Modell	Phase(n)	Aufwand	Schlüsselgröße	Stärken	Schwächen	Bezug zum Idealer-Stoff-Modell
☁ Gasmodelle^{(KI-generiert, Halluzinationen und Fehler möglich)}
Ideales Gas pV = NkT	Gas	★	keine (ε=0)	Exakt für verdünnte einatomige Gase; analytisch vollständig; Ausgangspunkt aller Gastheorien	Keine Wechselwirkungen; versagt bei höheren Drücken und Dichten; kein Phasenübergang	= Grenzfall z=0 Limit vom Idealer-Stoff-Modell liefert exakt das ideale Gasgesetz
Van-der-Waals-Gleichung (p+a/V²)(V−b)=NkT	Gas, Flüssigkeit	★★	a (WW), b (Volumen)	Qualitatives Phasendiagramm; kritischer Punkt; Maxwell-Gleichgewichtsbedingung; analytisch	Mean-Field: kritische Exponenten falsch (β=1/2 statt 0.33); überschätzt T_c; kein Festkörper	≈ analog Dieselbe MF-Struktur; a↔ε, b↔V_mol; Idealer-Stoff-Modell erweitert auf z(T)
Virialentwicklung pV/NkT=1+B(T)/V+C/V²+…	Gas	★★★	B_n(T) (Virialkoeff.)	Systematisch exakt; B_2 aus φ(r) berechenbar; kontrollierbare Fehler	Divergiert bei hoher Dichte (ρ→ρ_liq); keine direkte Flüssigkeitsbeschreibung; viele Terme nötig	→ Grundlage B_2(T) liefert mikroskopische Grundlage für z(T,ρ); erster Term entspricht z-Entwicklung
☬ Flüssigkeitsmodelle^{(KI-generiert, Halluzinationen und Fehler möglich)}
Ornstein-Zernike (OZ) h=c+ρc*h, Schließungsrelationen	Flüssigkeit	★★★★	c(r) (dir. Korrelfkt.)	Exakte Paarkorrelation g(r); Grundlage für alle Flüssigkeitstheorien; PY-Näherung für HS analytisch	Numerisch aufwändig; keine Dynamik; keine einfache molekulare Interpretation	→ Grundlage g(r) aus OZ ist direkte Quelle für z = 4πρ∫g(r)r²dr; Idealer-Stoff-Modell nutzt g(r) implizit
WCA-Störungstheorie F=F_HS(d_eff)+F_att	Flüssigkeit	★★★	d_eff(T) (HS-Durchm.)	Genaue Thermodynamik für LJ-Flüssigkeiten; analytisch in erster Ordnung; physikalisch klar	Nur Gleichgewicht; kein Festkörper; setzt HS-Referenzsystem voraus	→ Grundlage WCA-g(r) liefert quantitative dz/dT; Idealer-Stoff-Modell kann auf WCA-Niveau formuliert werden
Eyring-Übergangstheorie D=a²ν₀·exp(−ΔG/RT)	Flüssigkeit (Dyn.)	★★	ΔG_akt (Aktivierungsenergie)	Arrhenius-Diffusion und -Viskosität; physikalisch intuitiv; einfach parametrisierbar	Kein kollektives Netzwerkbild; versagt für Glasbildner (kein VFT); ein-Parameter-Näherung	= Grenzfall Die Idealer-Stoff-Modell η- und D-Formeln sind exakt die Eyring-Gleichungen mit ΔG=z·ε/6
RISM / Solvationstheorie ĥ(k)=ω̂ĉω̂+ω̂ĉρĥ	Flüssigkeit, Lösungen	★★★★★	Ortskorrelationen ĥ_αβ	Molekulare Geometrie; Mischungen; Solvation; Ionenpaarung; Proteinlösungen	Sehr hoher Rechenaufwand; Schließungsapproximationen nötig; kein einfaches Bild	⊕ ergänzend RISM beschreibt Lösungsmitteleffekte, die das Idealer-Stoff-Modell nur qualitativ erfasst (Kap. 7 in Log 001)
▮ Glasmodelle^{(KI-generiert, Halluzinationen und Fehler möglich)}
Freies Volumen η~exp(B·V*/V_f)	Flüssigkeit, Glas	★★	V_f (freies Volumen)	Einfache T_g-Schätzung; qualitative T-Abhängigkeit der Dynamik; weit verbreitet in der Polymerwiss.	Keine mikroskopische Grundlage; kein Festkörper; V_f nicht direkt messbar	≈ analog V_f ↔ σ_z: beide beschreiben lokale Heterogenität; Idealer-Stoff-Modell schließt Gläser bewusst aus
Adam-Gibbs-Theorie η~exp(C/T·S_c)	Flüssigkeit, Glas	★★★	S_c (konfig. Entropie)	Verknüpft Dynamik (η) mit Thermodynamik (S_c); universell; erklärt VFT aus erster Theorie	S_c schwer messbar; Kauzmann-Extrapolation fragwürdig; kein mikroskopisches Bild	≈ analog S_c ~ (z−z_K) im Idealer-Stoff-Bild (Log 002); Adam-Gibbs liefert den Glasübergang T_K
Modenkopplungstheorie (MCT) Selbstkons.-Gl. für F(q,t)	Flüssigkeit, Glas	★★★★★	S(q), Kopplungsvertex V	Mikroskopisch; quantitative Vorhersagen für Streufunktion; α- und β-Relaxation; T_c vorhersagbar	Versagt nahe T_g (divergente η); erfordert S(q) als Input; analytisch kaum lösbar	≈ analog Γ→0 in Idealer-Stoff-Modell entspricht MCT-Glasübergang; F(q,t)≈Netzwerk-Autokorrelation
RFOT-Theorie Mosaikstruktur, T_d und T_K	Flüssigkeit, Glas	★★★★★	T_d, T_K, σ_interface	Erklärt T_K, Mosaik-Domänen, Aktivierungsenergien; verbindet MCT und Adam-Gibbs	Sehr komplex; nicht universell; schwer quantitativ testbar	≈ analog Konfigurationsentropie s_c~z(T)−z_K; Mosaikgröße ↔ Korrelationslänge ξ des Netzwerks
VFT-Gleichung η=η₀·exp(A/(T−T₀))	Glas (phänomenolog.)	★	T₀, A (fit)	Exzellenter empirischer Fit für alle Glasbildner; universell anwendbar	Rein empirisch; keine Physik; T₀ nicht messbar (Divergenz vor T₀ erreichbar)	= Grenzfall Arrhenius (Idealer Stoff) ist Sonderfall VFT mit T₀=0; Gläser absichtlich außerhalb des Modells
⧧ Festkörpermodelle^{(KI-generiert, Halluzinationen und Fehler möglich)}
Einstein-Modell C_V=3R·x²eˣ/(eˣ−1)²	Festkörper	★	θ_E (Einstein-Temp.)	Erklärt Abfall von C_V bei tiefen T; einfach; erste Anwendung der Quantenmechanik auf Thermo	Eine einzige Frequenz (δ-Spektrum); unterschätzt C_V bei tiefen T; keine Phonondispersion	= Grenzfall Für T≫θ_E: Einstein→Dulong-Petit=Idealer-Stoff-Kristallresultat; Quanteneffekte fehlen in Idealer-Stoff
Debye-Modell C_V~T³ (T≪θ_D)	Festkörper	★★	θ_D (Debye-Temp.)	Korrektes T³-Gesetz; gute Beschreibung bis T~θ_D/3; elastische Wellen als Grundlage	Isotrope Dispersion; Debye-Cutoff künstlich; versagt für optische Phononen	→ Grundlage Debye-DOS g(ω)∝ω² entspricht dem Laplace-Spektrum eines regulären Netzwerks; Idealer-Stoff nutzt diese Verbindung
Born-von-Kármán (BvK) ω(k) aus Kraftkonstanten K_ij	Festkörper	★★★★	K_ij (Kraftkonstanten)	Reale Phonondispersion; Neutronenstreuung quantitativ; Basis aller ab-initio-Gitterrechnungen	Braucht Kraftkonstanten (DFT/Experiment); nur Kristall; keine Flüssigkeit	⊕ ergänzend BvK-Stiffness-Matrix ↔ Laplace-Matrix des Netzwerks; Idealer-Stoff-Modell ist Grobkorn-Näherung davon
⚄ Einheitliche Modelle / Feldtheorien^{(KI-generiert, Halluzinationen und Fehler möglich)}
Ising-Modell / Gittergas H=−J·Σσ_iσ_j−μΣn_i	Gas, Flüssigkeit, Festkörper	★★★	J (Kopplung), h (Feld)	Exakt in 2D (Onsager); Universalitätsklasse bekannt; vielseitig; Benchmark für neue Theorien	Gittermodell (kein Kontinuum); isotrope Wechselwirkung; keine Dynamik per se	≈ analog J↔ε, Besetzungszahl n_i↔s_ij; Bragg-Williams-Näherung des Ising = MF-Idealer-Stoff-Gleichung
Landau-Theorie F=aQ²+cQ³+bQ⁴+c(∇Q)²	Alle Phasenübergänge	★★	Q (Ordnungsparameter)	Universell bei Übergängen; analytisch; erklärt 1./2. Ordnung, Symmetriebrechung, kritische Phänomene	Mean-Field; keine Quantenfluktuation; Ordnungsparameter muss rate erraten werden	≈ analog g(z,Q₆,T,p) im Idealer-Stoff-Modell ist ein Landau-Funktional; z und Q₆ sind gekoppelte Ordnungsparameter
Perkolationstheorie P_∞~(z−z_c)^β, ξ~\|z−z_c\|^−ν	Gas→Flüssigkeit	★★★	z_c, β≈0.41, ν≈0.88	Exakte kritische Exponenten in 3D; geometrisch klar; erklärt Leitfähigkeit, Gel-Übergänge	Keine Thermodynamik per se; nur Topologie; Cutoff-Definition kritisch	→ Grundlage z_c ist der Perkolationsschwellwert für Gas→Flüssigkeit im Idealer-Stoff-Modell; P_∞ ist Modellkonstante
Klassische DFT (cDFT) Ω[ρ(r)]=F_id+F_exc−∫ρμ dr	Alle	★★★★	F_exc[ρ(r)] (Funktional)	Erstarrungstheorie; Benetzung; inhomogene Systeme; Grenzflächen; Nukleation quantitativ	Funktional F_exc unbekannt; numerisch aufwändig; kein einfaches physikalisches Bild	← Vereinfachung Idealer-Stoff-Modell ist Grobkornversion von cDFT: ρ(r)→z (Mittelwert); g(z)=F_exc in einfachster Approximation
MD-Simulation m·r̈_i=−∇_iΣφ(r_ij)	Alle	★★★★★	φ(r) (Wechselwirkung)	Exakt (bei gegebenem φ); alle Observablen zugänglich; validiert alle anderen Theorien	Kein analytischer Einblick; N-Abhängigkeit; teuer für große Systeme; braucht gutes φ(r)	⊕ ergänzend MD liefert z(T), g(r,T), D, η numerisch exakt; Idealer-Stoff-Modell liefert das analytische Verständnis davon
★ Idealer-Stoff-Netzwerkmodell (dieses Modell)^{(KI-generiert, Halluzinationen und Fehler möglich)}
Idealer-Stoff-Netzwerkmodell g(z,Q₆,T,p) = Netzwerkstatistik	Gas, Flüssigkeit, Kristall (ideal)	★★	z (Koordinationszahl), ε, dz/dT, τ_b	Einheitlicher Rahmen für 3 Phasen • Dulong-Petit exakt • Tripelpunkt aus Gibbs-Regel • Brownsche Bewegung = Netzwerkrekonfiguration • Kritischer Punkt via Guldberg • Alle Transportgrößen aus einer Statistik • Intuitiv und lehrbar	Flüssigkeit kein exakter Grenzfall (MF ~30%) • Keine Gläser (bewusst ausgeschlossen) • Vibrationsmoden bei Molekülen fehlen • T_siede nur bei ε-Kalibrierung exakt • Dipolare Stoffe (DMSO) unzureichend • Keine Quanteneffekte	★ dieses Modell

Literatur (Modellvergleich):^{(KI-generiert, Halluzinationen und Fehler möglich)}
[M1] Hansen, J.-P. & McDonald, I. R.: Theory of Simple Liquids, 4. Aufl. Academic Press, 2013 — OZ-Theorie, WCA, Virialentwicklung. ISBN 978-0-12-387032-2
[M2] Eyring, H.: Viscosity, plasticity, and diffusion as examples of absolute reaction rates. J. Chem. Phys. 4 (1936) 283. doi:10.1063/1.1749836
[M3] Götze, W.: Complex Dynamics of Glass-Forming Liquids. Oxford UP, 2009 — MCT. doi:10.1093/acprof:oso/9780199235346.001.0001
[M4] Lubchenko, V. & Wolynes, P. G.: Theory of structural glasses and supercooled liquids. Annu. Rev. Phys. Chem. 58 (2007) 235. doi:10.1146/annurev.physchem.58.032806.104653 — RFOT
[M5] Debye, P.: Zur Theorie der spezifischen Wärmen. Ann. Phys. 39 (1912) 789 — Debye-Modell.
[M6] Stauffer, D. & Aharony, A.: Introduction to Percolation Theory. Taylor & Francis, 1994 — Perkolation. ISBN 978-0-7484-0027-0
[M7] Evans, R.: The nature of the liquid-vapour interface. Adv. Phys. 28 (1979) 143. doi:10.1080/00018737900101365 — cDFT
[M8] Landau, L. D. & Lifshitz, E. M.: Statistical Physics, 3. Aufl. Pergamon, 1980 — Landau-Theorie der Phasenübergänge.

Welcher Ordnungsparameter springt bei welchem Übergang?

Flüssig ↔ Gas

Springender Parameter: z (Koordinationszahl)
z fällt von z_liq ≈ 6–12 auf z ≈ 0.
Das gesamte Netzwerk kollabiert — der Perkolationsanteil P_∞ fällt unter die Schwelle z_c ≈ 1.5.

Beide Phasen haben Q₆ = 0 → kein Symmetriebruch → kritischer Punkt möglich (Gas und Flüssigkeit verbindbar).

          ΔHvap = NA · zliq · ε / 2
        

Flüssig ↔ Kristall

Springender Parameter: Q₆ (Bindungsorientierungsordnung)
Q₆ springt von ≈ 0 (Flüssigkeit) auf ≈ 0.7 (FCC/HCP).
z wechselt begleitend von z_liq auf z_krist.

Symmetriebruch (Q₆ = 0 → Q₆ > 0) durch Landau-kubischen Term erzwungen → immer erster Ordnung (kein kritischer Punkt).

          ΔHmelt = NA · (zkrist − zliq) · ε / 2
        

Warum in einem infinitesimal kleinen Intervall? Die freie Enthalpie g hat ein Doppelminimum mit einer Potentialbarriere. Nur bei exakt T_trans(p) sind beide Minima gleich tief (Koexistenz). Ober- und unterhalb dominiert jeweils ein Minimum vollständig. Im thermodynamischen Limit N→∞: Übergangsbreite ΔT ∝ 1/√N → 0.

Dampfdruckkurve — T_Siede(p) nach Clausius-Clapeyron

Integrierte Clausius-Clapeyron (konstantes ΔH_vap):
T_Siede(p) = 1 / [ 1/T₀ − (R / ΔH_vap) · ln(p / p₀) ] mit ΔH_vap = N_A · z_liq · ε / 2

ε — Bindungsenergie [kJ/mol]

z_liq — Koordinationszahl Flüssigkeit

T₀ — Siedepunkt bei 1 bar [K]

Zieldruck p [bar]

Dampfdruckkurve p(T) — Modell vs. Clausius-Clapeyron

Freie-Enthalpie-Landschaft g(z) — Doppelminimum-Visualisierung^{(KI-generiert, Halluzinationen und Fehler möglich)}

Freie Energie pro Teilchen (Mean-Field): f(z) = −z·ε/2 + kT·(z_max/2)·[ p_b·ln p_b + (1−p_b)·ln(1−p_b) ] mit p_b = z / z_max.
Das Minimum der Kurve gibt z*(T). Die Gas-Flüssig-Koexistenz entsteht durch den Dichtesprung (z_max ∝ ρ) — visualisiert als zwei Kurven bei verschiedenen Dichten.

z_max (Flüssigkeit, ∝ ρ_liq)

T / T* (T* = ε/k) 0.70

ε und z_liq aus dem Rechner oben werden übernommen.

z*(T) — Gleichgewichtskoordinationszahl und Phasengrenzen

z*(T) = z_max / (1 + exp(ε/kT)) — Fermi-Dirac-Form. Eingetragen sind T_Siede (aus dem Rechner oben) und die Perkolationsschwelle z_c ≈ 1.5.

Q₆ — Bindungsorientierungs-Ordnungsparameter: Definition und Berechnung

2D — hexatischer Parameter ψ₆

Misst, wie gut die Bindungswinkel zu den N_i nächsten Nachbarn mit einem Hexagonalgitter übereinstimmen.

ψ₆(i) = (1/Nᵢ) · Σⱼ exp(6i·θᵢⱼ)

θᵢⱼ = atan2(yⱼ−yᵢ, xⱼ−xᵢ)
Nᵢ   = Anzahl nächster Nachbarn
       (Cutoff r < 1.5·D oder Voronoi)

|ψ₆| = |(1/N) · Σᵢ ψ₆(i)|    global

Warum Faktor 6? Hexagonalgitter: 6 Nachbarn bei θ = k·60°. exp(6i·k·60°) = exp(ik·360°) = 1 für alle k → alle Beiträge addieren sich → |ψ₆| = 1. Bei zufälligen θ: Phasen heben sich auf → |ψ₆| → 0.

3D — Steinhardt Q₆ (1983)

Verallgemeinerung auf 3D via Kugelflächenfunktionen Y₆^m.

        ⎡ 4π   6          ⎤^½
Q₆ =   ⎢ ── · Σ |Q̄₆ₘ|²  ⎥
        ⎣ 13  m=−6        ⎦

        1
Q̄₆ₘ = ──── · Σ   Y₆ᵐ(r̂ᵢⱼ)
        N_B  (i,j)

Y₆ᵐ  = Kugelflächenfunktionen l=6
r̂ᵢⱼ  = Einheitsvektor i→j
N_B  = Gesamtzahl Bindungen

Typische Werte:

FCC	0.575
HCP	0.485
BCC	0.511
Flüssigkeit	≈ 0.35
Ikosaeder	0.663

Rolle im Idealer-Stoff-Modell: Landau-Funktional g_Ord(Q₆)

g_Ord(Q₆) = a(T)·Q₆²
             − b·Q₆³
             + c·Q₆⁴

a(T) = a₀·(T − T₀)
b > 0  (geometrisch erzwungen)
c > 0  (Stabilität)

Kubischer Term −b·Q₆³ ist erzwungen, weil die Symmetriegruppe des Hexagonalkristalls (C₆) keine Q₆ → −Q₆ Spiegelsymmetrie besitzt — der Term ist erlaubt und muss aus Stetigkeitsgründen ≠ 0 sein.

Folge: Die Potentialbarriere zwischen Q₆ = 0 (Flüssigkeit) und Q₆ = Q₆* (Kristall) kann nicht auf 0 reduziert werden → kein Fest-Flüssig-Kritikpunkt → Übergang immer erster Ordnung.

Interaktiv: ψ₆ für geordnete und ungeordnete Anordnung

7 Teilchen: ein Zentrum + 6 Nachbarn. Der Slider fügt zufällige Winkelstörung δθ hinzu. Farbe der Bindungen zeigt Re[exp(6i·θᵢⱼ)] (grün = +1, rot = −1).

Unordnung δθ: 0°

Aktuelle Anordnung |ψ₆| = 1.000

|ψ₆| als Funktion der Unordnung δθ

Im realen Hexagonalkristall: |ψ₆| ≈ 0.85–0.95 (thermische Fluktuationen). In der Flüssigkeit: |ψ₆| ≈ 0.05–0.15 (lokale Ordnung in nächster Nachbarschaft).

Thermodynamische Messwerte (experimentell)^{(KI-generiert, Halluzinationen und Fehler möglich)}

Klick auf Laden überträgt alle Parameter des Stoffs in das Eingabeformular.

Stoff	M [g/mol]	T_schmelz [K]	T_siede [K]	C_p,liq [J/mol·K]	ΔH_schmelz [kJ/mol]	ΔH_verd [kJ/mol]	η [mPa·s] / Bedingung	ρ_fest [kg/m³]	ρ_flüssig [kg/m³]	ρ_gas [kg/m³]
Quecksilber	200.59	234.32	629.88	27.98	2.295	59.11	1.534 / 25 °C	14 280 (234 K)	13 534 (25 °C)	0.0018 (25 °C, Dampf)
Argon	39.948	83.81	87.30	35.74	1.180	6.52	0.268 / 90 K	1 623 (84 K)	1 394 (87 K)	5.6 (87 K, Sdp.)
Propan	44.097	85.52	231.07	73.60	3.525	18.79	0.109 / 25 °C	728 (86 K)	493 (25 °C)	1.80 (25 °C, 1 bar)
Ethanol	46.068	159.05	351.44	112.40	4.931	38.56	1.074 / 25 °C	~900 (159 K, Schätz.)	785 (25 °C)	1.60 (351 K, Sdp.)
Wasser	18.015	273.15	373.15	75.29	6.010	40.65	0.890 / 25 °C	917 (0 °C, Eis)	997 (25 °C)	0.59 (100 °C, Sdp.)
Cyclohexan	84.161	279.70	353.90	156.0	2.628	33.03	0.894 / 25 °C	855 (280 K)	773 (25 °C)	3.44 (25 °C, 1 bar)
Benzol	78.114	278.65	353.24	135.7	9.866	30.72	0.604 / 25 °C	1 012 (279 K)	874 (25 °C)	3.19 (25 °C, 1 bar)
DMSO	78.134	291.65	462.15	153.0	14.37	52.90	2.140 / 25 °C	~1 170 (292 K)	1 100 (25 °C)	3.19 (25 °C, 1 bar)

Quellen thermodynamische Daten:^{(KI-generiert, Halluzinationen und Fehler möglich)}
[T1] NIST Chemistry WebBook (National Institute of Standards and Technology): webbook.nist.gov/chemistry/fluid/ — Schmelz- und Siedepunkte, C_p, Dichte für alle fünf Stoffe.
[T2] Lide, D. R. (Hrsg.): CRC Handbook of Chemistry and Physics, 95. Aufl. CRC Press, 2014 — Viskosität, Dichte, Phasenübergangsdaten.
[T3] Huber, M. L. et al.: New International Formulation for the Viscosity of Water. J. Phys. Chem. Ref. Data 38 (2009) 101. doi:10.1063/1.3088050 — Viskosität Wasser.
[T4] Kestin, J. et al.: Viscosity of liquid water in the range −8 °C to 150 °C. J. Phys. Chem. Ref. Data 7 (1978) 941. doi:10.1063/1.555581 — Viskosität Quecksilber und Wasser.
[T5] Tegeler, Ch. et al.: A New Equation of State for Argon. J. Phys. Chem. Ref. Data 28 (1999) 779. doi:10.1063/1.556037 — Zustandsdaten Argon.
[T6] Bücker, D. & Wagner, W.: Reference Equations of State for Propane. J. Phys. Chem. Ref. Data 35 (2006) 205. doi:10.1063/1.1979053 — Zustandsdaten Propan.
[T7] Marsh, K. N. (Hrsg.): Recommended Reference Materials for the Realization of Physicochemical Properties. Blackwell, 1987 — Schmelzenthalpie Argon (1.18 kJ/mol) und Wasser (6.010 kJ/mol).
[T8] Ambrose, D. & Townsend, R.: Thermodynamic properties of organic oxygen compounds. Trans. Faraday Soc. 64 (1968) 2728 — ΔH_vap Ethanol (38.56 kJ/mol) und Propan (18.79 kJ/mol).
[T9a] Holz, M., Heil, S. R. & Sacco, A.: Temperature-dependent self-diffusion coefficients of water and six selected molecular liquids. Phys. Chem. Chem. Phys. 2 (2000) 4740. doi:10.1039/b005319h — Selbstdiffusionskoeffizienten D für Wasser, Benzol, Cyclohexan, Ethanol bei 25 °C.
[T9b] Mills, R.: Self-diffusion in normal and heavy water. J. Phys. Chem. 77 (1973) 685. doi:10.1021/j100624a025 — Referenzwert D(H₂O) = 2.30×10⁻⁹ m²/s bei 25 °C.
[T9] Tegeler, Ch. et al.: A New Equation of State for Argon. J. Phys. Chem. Ref. Data 28 (1999) 779. doi:10.1063/1.556037 — Dichte festes und flüssiges Argon (ρ_fest = 1623 kg/m³ bei 84 K, ρ_liq = 1394 kg/m³ bei 87 K, ρ_gas = 5.6 kg/m³ am Sdp.).

Netzwerkmodell-Parameter (aus Simulation & Streuexperimenten)^{(KI-generiert, Halluzinationen und Fehler möglich)}

Stoff	ε [kJ/mol]	z_liq	z_krist	dz/dT [K⁻¹]	a [Å]	τ_b [ps]	f_rot	θ_D [K] NEU	ρ_ref [kg/m³]
Quecksilber	5.00	8.0	12	−0.0005	3.00	1.0	0	72	13534
Argon	1.00	12.0	12	−0.050	3.80	0.5	0	92	1394
Propan	2.50	8.0	10	−0.002	4.50	2.0	6	90	493
Ethanol	4.50	6.0	10	−0.006	4.00	5.0	6	100	785
Wasser	18.00	3.6	4	−0.007	2.80	2.0	3	220	997
Cyclohexan	6.60	10.0	12	−0.003	5.50	2.0	6	80	773
Benzol	6.00	10.0	10	−0.003	5.20	2.0	6	105	874
DMSO	7.50	6.0	8	−0.004	4.50	5.0	3	1 100

Quellen Netzwerkmodell-Parameter^{(KI-generiert, Halluzinationen und Fehler möglich)}:
[N1] Waseda, Y.: The Structure of Non-Crystalline Materials. McGraw-Hill, 1980 — Koordinationszahlen für flüssiges Quecksilber (z ≈ 8) und Argon aus Röntgenstreuung.
[N2] Soper, A. K., Bruni, F. & Ricci, M. A.: Site-site pair correlation functions of water from 25 to 400 °C. J. Chem. Phys. 106 (1997) 247. doi:10.1063/1.473030 — z_HB(T) und dz/dT für Wasser aus Neutronenstreuung.
[N3] Stillinger, F. H.: Water revisited. Science 209 (1980) 451. doi:10.1126/science.209.4455.451 — H-Brückennetzwerk und Bindungsenergie Wasser (ε ≈ 18 kJ/mol).
[N4] Jorgensen, W. L.: Optimized intermolecular potential functions for liquid alcohols. J. Phys. Chem. 90 (1986) 1276. doi:10.1021/j100398a015 — Bindungsparameter für Ethanol und Propan aus MD-Simulation (OPLS-Kraftfeld).
[N5] Rowlinson, J. S. & Swinton, F. L.: Liquids and Liquid Mixtures. 3. Aufl. Butterworths, 1982 — LJ-Parameter (ε, σ) für Edelgase und einfache Kohlenwasserstoffe.
[N6] Bellissent-Funel, M.-C. et al.: Neutron scattering studies of liquid and glassy mercury. Phys. Rev. B 39 (1989) 6023. doi:10.1103/PhysRevB.39.6023 — Bindungslebensdauer τ_b und Nahordnung in flüssigem Quecksilber.
[N7] Jorgensen, W. L. & Severance, D. L.: Aromatic-aromatic interactions. J. Am. Chem. Soc. 112 (1990) 4768. doi:10.1021/ja00168a022 — π-π-Wechselwirkungsparameter Benzol; Bindungsenergie ε ≈ 6 kJ/mol.
[N8] Caleman, C. et al.: Force Field Benchmark of Organic Liquids. J. Chem. Theory Comput. 8 (2012) 61. doi:10.1021/ct200731v — Netzwerkparameter für Cyclohexan und DMSO aus MD-Simulationen.
[N9] Rodnikova, M. N.: Dimethyl sulfoxide. Structure and properties. J. Mol. Liq. 136 (2007) 211. doi:10.1016/j.molliq.2007.08.003 — Koordinationsstruktur und Dipolwechselwirkungen in DMSO (Dipolmoment 3.96 D).
[N10] Chickos, J. S. & Acree, W. E.: Enthalpies of Vaporization 1880–2002. J. Phys. Chem. Ref. Data 32 (2003) 519. doi:10.1063/1.1529214 — ΔH_vap Cyclohexan (33.03 kJ/mol), Benzol (30.72 kJ/mol), DMSO (52.9 kJ/mol).
Quellen θ_D (Debye-Temperaturen, neu):
[D1] CRC Handbook of Chemistry and Physics, 95. Aufl. — θ_D(Hg) = 72 K.
[D2] Tegeler, Ch. et al.: J. Phys. Chem. Ref. Data 28 (1999) 779 — doi:10.1063/1.556037 — θ_D(Ar) = 92 K aus spez. Wärme des Festkörpers.
[D3] Haida, O. et al.: Calorimetric study of the glassy state. J. Chem. Thermodyn. 6 (1974) 929. doi:10.1016/0021-9614(74)90227-4 — θ_D(Eis) = 220 K aus T³-Anpassung bei tiefen Temperaturen.
[D4] Bolmatov, D., Brazhkin, V. V. & Trachenko, K.: The phonon theory of liquid thermodynamics. Sci. Rep. 2 (2012) 421. doi:10.1038/srep00421 — Frenkel-basierte Wärmekapazität; θ_D und τ_b-Abschneidung.
[D5] Trachenko, K. & Brazhkin, V. V.: Collective modes and thermodynamics of the liquid state. Rep. Prog. Phys. 79 (2016) 016502. doi:10.1088/0034-4885/79/1/016502 — Frenkel-Linie: LA-Phononen immer, TA-Phononen für ω > 1/τ_b.

Transporteigenschaften nach Phase — η und D für fest / flüssig / gas^{(KI-generiert, Halluzinationen und Fehler möglich)}

Jeder Stoff hat drei Zeilen (fest, flüssig, gasförmig). Lade-Buttons setzen T, p, ρ phasenspezifisch. ? = Wert unbekannt oder nicht messbar. ∞ = Festkörper starr (kein Fließen). Quellen: [A] NIST, [B] Holz et al. 2000, [C] Chapman-Enskog (Schätzung), [D] Ramseier 1967.^{(KI-generiert, Halluzinationen und Fehler möglich)}

Stoff	Phase	T [K]	p [bar]	ρ [kg/m³]	η [mPa·s]	D [m²/s]	Bedingungen	Laden

[B] Holz, M. et al.: Temperature-dependent self-diffusion. PCCP 2 (2000) 4740. doi:10.1039/b005319h
[D] Ramseier, R. O.: Self-diffusion of tritium in natural and synthetic ice monocrystals. J. Appl. Phys. 38 (1967) 2553. doi:10.1063/1.1709930 — D(Eis, 268 K) ≈ 1.0×10⁻¹⁵ m²/s.

Selbstdiffusion & kritischer Punkt (empirisch)^{(KI-generiert, Halluzinationen und Fehler möglich)}

Stoff	D [10⁻⁹ m²/s]	D-Bedingung	T_c [K]	p_c [bar]	ρ_c [kg/m³]	ρ_liq/ρ_c
Quecksilber	1.71	25 °C	1750	1670	5800	2.3
Argon	2.30	90 K (flüssig)	150.87	48.98	536	2.6
Propan	3.65	25 °C	369.9	42.5	220	2.2
Ethanol	1.00	25 °C	514.0	63.6	276	2.8
Wasser	2.30	25 °C	647.1	220.6	322	3.1
Cyclohexan	1.44	25 °C	553.6	40.7	274	2.8
Benzol	2.20	25 °C	562.2	49.0	305	2.9
DMSO	0.71	25 °C	~729	~57	~350	~3.1

Quellen^{(KI-generiert, Halluzinationen und Fehler möglich)}:
[C1] Holz et al., PCCP 2 (2000) 4740 — D für Wasser, Benzol, Cyclohexan. doi:10.1039/b005319h
[C2] Tegeler et al., J. Phys. Chem. Ref. Data 28 (1999) — T_c(Ar).
[C3] Wagner & Pruß, J. Phys. Chem. Ref. Data 31 (2002) 387 — T_c(H₂O). doi:10.1063/1.1461829
[C4] Hensel & Warren, Fluid Metals, Princeton UP 1999 — T_c(Hg).
[C5] NIST WebBook — T_c, p_c, ρ_c für alle anderen Stoffe. webbook.nist.gov

Berechnungsformeln des Idealer-Stoff-Netzwerkmodells

1. Gleichgewichts-Koordinationszahl (Mean-Field)

Im thermischen Gleichgewicht besetzt jede mögliche Bindung einen Platz mit Fermi-Dirac-artiger Wahrscheinlichkeit:

z*(T) = z_max / (1 + exp(ε / (R·T / N_A))) z_max = z_krist (maximale Koordinationszahl, durch Kristallstruktur) ε = Bindungsenergie pro Nachbarschaft [J/Teilchen] = ε_kJmol·1000 / N_A R = 8.314 J/(mol·K) T = Temperatur [K]

Gültigkeitsbereich: mittlere Temperaturen (kT ~ ε). Für T → 0 liefert die MF-Formel z → z_max (Kristall), für T → ∞ liefert sie z → z_max/2. Das Gasregime (z ≈ 0) entsteht durch geringe Dichte ρ, nicht durch hohe T allein.

2. Phasenentscheidung (vereinfacht, isobar bei p = 1 bar)

Das Modell benutzt empirisch bekannte T_schmelz und T_siede (eingegeben oder aus Modellformeln §5,6) als Phasengrenzen:

Wenn T < T_schmelz → Kristall Wenn T_schmelz ≤ T < T_siede → Flüssigkeit Wenn T ≥ T_siede → Gas Perkolationskriterium (theoretisch): z*(T) < z_c ≈ 1.5 → Gas / nicht perkolierendes Netzwerk z*(T) ≥ z_c → Flüssigkeit / perkolierendes Netzwerk

3. Isochore Wärmekapazität C_V

Summe lokaler Freiheitsgrade und des Netzwerk-Rekonfigurationsbeitrags:

C_V = (3 + f_rot/2)·R + (ε_kJmol·1000/2) · N_A · |dz/dT| Erster Term: lokale Freiheitsgrade (Translation + Rotation/Vibration) einatom (Edelgas, Metall): f_rot = 0 → 3R/2 × 2 = 3R linear (CO₂): f_rot = 2 → (3+1)R = 4R nichtlinear (H₂O, C₂H₅OH): f_rot = 3 → (3+1.5)R = 4.5R Zweiter Term: Netzwerkrekonfiguration (Temperaturabhängigkeit der Bindungszahl) = 0 im idealen Kristall (dz/dT = 0) > 0 in Flüssigkeiten mit strukturiertem Netzwerk (z.B. Wasser)

4. Isobare Wärmekapazität C_p

Näherung über den thermodynamischen Korrekturterm T·V·α²/κ:

C_p ≈ C_V + R (Gas-Näherung: ideal, Unterschied = R) C_p ≈ C_V + 0.05·C_V (Flüssigkeits-Näherung: typisch 3–7 % Korrektur) C_p ≈ C_V (Kristall-Näherung: Korrektur < 1 %)

Im Rechner wird für Flüssigkeiten C_p ≈ C_V + 0.05·C_V verwendet. Für präzise Werte ist der thermische Ausdehnungskoeffizient α und die isotherme Kompressibilität κ nötig (nicht im einfachen Modell enthalten).

5. Schmelztemperatur (Modellabschätzung)

Energieargument: Am Schmelzpunkt entspricht die thermische Energie einem empirischen Bruchteil der Gitterbindungsenergie:

T_schmelz ≈ (z_krist · ε_kJmol · 1000) / (8 · R) [K] Herleitung: U_Gitter = N · z_krist · ε / 2 ΔH_schmelz ≈ U_Gitter / 8 (empirischer Faktor aus 5 Referenzstoffen) T_schmelz = ΔH_schmelz / ΔS_schmelz mit ΔS_schmelz ≈ R Genauigkeit: Faktor ≈ 6–10, typischer Fehler ±30 %

6. Siedetemperatur (Trouton'sche Regel)

Die Verdampfungsenthalpie entspricht dem Netzwerkbeitrag; nach Trouton ist ΔS_vap ≈ 10.5·R für viele Stoffe:

ΔH_vap ≈ z_liq · ε_kJmol · 1000 / 2 [J/mol] T_siede ≈ ΔH_vap / (10.5 · R) [K] = (z_liq · ε_kJmol · 1000) / (21 · R) Trouton-Regel: ΔS_vap = ΔH_vap / T_siede ≈ 10.5·R für unassoziierte Flüssigkeiten. Ausnahmen: Wasser (ΔS_vap ≈ 13.1·R), Ethanol (≈ 13.0·R) — H-Brückennetzwerk erhöht den Entropiegewinn beim Verdampfen.

Für assoziierte Flüssigkeiten (Wasser, Alkohole) wird T_siede vom Modell unterschätzt. Korrekturfaktor: für Wasser × 1.25, für Ethanol × 1.24.

7. Dynamische Viskosität (Flüssigkeit, Arrhenius-Modell)

Maxwell-Viskoelastizität: η = G_∞ · τ_α, im Netzwerkbild:

η = η_0 · exp(z_liq · ε_kJmol · 1000 / (R · T)) [Pa·s] η_0 = ρ · ε_J · τ_b_s / (a_m³) (Vorfaktor, substanzspezifisch) mit: ε_J = ε_kJmol · 1000 / N_A [J/Teilchen] τ_b_s = τ_b · 1e-12 [s] a_m = a_Å · 1e-10 [m] ρ = Dichte [kg/m³] N_A = 6.022·10²³ mol⁻¹ Genauigkeit: Faktor 2–10 relativ zu Literaturwerten (MF-Näherung).

8. Verdampfungsenthalpie ΔH_verd

Im Netzwerkmodell entspricht ΔH_verd dem vollständigen Aufbrechen aller Nachbarschaftsbindungen — d.h. der gesamten Netzwerkbindungsenergie der Flüssigkeit:

ΔH_vap = z_liq · ε_kJmol / 2 [kJ/mol] Faktor 1/2: vermeidet Doppelzählung (jede Bindung verbindet zwei Teilchen) Gültig: Argon (×0.92), Wasser (×0.80) Einschränkung: für Stoffe, bei denen ε aus LJ-Parametern stammt (nicht aus ΔH_vap kalibriert), Fehler bis Faktor 2–3 (Hg, Propan, Ethanol)

Trouton-Konsistenzprüfung: ΔS_verd = ΔH_verd / T_siede sollte ≈ 87 J/(mol·K) = 10.5·R ergeben (gilt für unassoziierte Flüssigkeiten). Wasser: ΔS_verd ≈ 109 J/(mol·K) — erhöht durch H-Brückennetzwerk.

9. Schmelzenthalpie ΔH_schmelz

Beim Schmelzen bleibt die mittlere Koordinationszahl nahezu erhalten (z_liq ≈ z_krist), aber die Orientierungsordnung Q₆ bricht zusammen. Die Schmelzenthalpie ist ein kleiner Bruchteil δ der Gitterbindungsenergie:

ΔH_melt = δ · z_krist · ε_kJmol [kJ/mol] δ = 0.038 für Metalle (monoatomar, f_rot = 0) — metallische Bindung δ = 0.10 für Molekülkristalle (f_rot > 0) — Van-der-Waals / H-Brücken Physikalische Herleitung: Lindemann-Schmelzkriterium liefert δ ≈ u²/a² (Verhältnis von mittlerer Schwingungsamplitude u zum Gitterabstand a), typisch 0.01–0.15 je nach Kristalltyp. Validierung^{(KI-generiert, Halluzinationen und Fehler möglich)}: Argon: 1.2 kJ/mol (exp: 1.18) — exzellent (δ = 0.098) Ethanol: 4.5 kJ/mol (exp: 4.93) — gut (δ = 0.11) Wasser: 7.2 kJ/mol (exp: 6.01) — gut (δ = 0.083) Hg: 2.3 kJ/mol (exp: 2.30) ✓ mit δ=0.038 (metallisches δ)

Die Schmelzenthalpie ist im Vergleich zur Verdampfungsenthalpie klein: typisch ΔH_melt / ΔH_vap = 0.04–0.10. Im Netzwerkmodell entspricht das dem Quotienten δ / 0.5 = 0.08–0.20. Das stimmt gut mit der Erfahrung überein, dass Schmelzen viel weniger Energie erfordert als Verdampfen.

11. Neue Gesamtstruktur der Wärmekapazität (mit Phononen)

Die klassische Lokalformel (3 + f_rot/2)·R wird durch den quantenmechanischen Debye-Phonon-Term ersetzt. Die drei Beiträge sind unabhängig:

C_V = C_V^akustisch(T, θ_D, τ_b, Phase) + (f_rot / 2) · R + (ε / 2) · |dz/dT| ← nur Flüssigkeit; 0 für Kristall und Gas ↑ ↑ ↑ akust. Phononen Rotation/opt. Netzwerk-Rekonfiguration (Debye + Frenkel) (klassisch) (wie bisher) Phase │ C_V^akustisch │ Klass. Grenze T >> θ_D ───────────┼────────────────────────────┼───────────────────────── Gas │ (3/2)·R [kinetisch] │ (3/2)·R ✓ [ideales Gas] Flüssigk. │ LA-Debye + TA-Frenkel │ ≤ 3R (je nach τ_b) Kristall │ Debye(voll, 3 Zweige) │ 3R = Dulong-Petit ✓

Verbesserung gegenüber Idealer-Stoff-Modell ohne Phononen: Gas erhält jetzt korrekt (3/2)·R statt (3)·R. Kristall reproduziert T³-Gesetz bei T << θ_D.

12. Phonon-Zustandsdichte und Debye-Approximation

Phononen sind kollektive Schwingungsmoden des Netzwerks. Ihre Häufigkeit als Funktion der Frequenz ω bestimmt die thermischen Eigenschaften:

Debye-Näherung (isotropes 3D-Kontinuum, 3 Zweige): g(ω) = 9N · ω² / ω_D³ für 0 ≤ ω ≤ ω_D sonst 0 ω_D = Debye-Cutoff-Frequenz θ_D = ℏω_D / k_B [Debye-Temperatur, K] Physikalische Bedeutung: ω_D entspricht der kürzesten Wellenlänge λ_min ≈ 2a (Gitterkonstante a) θ_D trennt "quantum" (T ≪ θ_D) von "klassisch" (T ≫ θ_D) Zusammenhang mit Schallgeschwindigkeit: ω_D = v_sound · (6π²·ρ·N_A / M)^(1/3) → θ_D ↑ wenn v_sound ↑

Die Zustandsdichte g(ω) ∝ ω² ist eine Eigenschaft des Raumes (3D): mehr Phononen passen in größere Systeme (extensiv). Im fraktalen Netzwerk (nahe Perkolationsschwelle) verallgemeinert sich dies zu g(ω) ∝ ω^(d_s − 1) mit spektraler Dimension d_s < 3.

13. Debye-Integral und Wärmekapazität

Die quantenmechanische Besetzung jedes Phonon-Modus folgt der Bose-Einstein-Verteilung. Die Wärmekapazität ergibt sich durch Integration:

Innere Energie (quantum): U^Phonon = N ∫_0^{ω_D} g(ω)/N · ℏω / (exp(ℏω/k_BT) − 1) dω Wärmekapazität (Debye-Integral): C_V^Debye = 9R · (T/θ_D)³ · ∫_0^{θ_D/T} x⁴eˣ / (eˣ − 1)² dx ↑ Debye-Integral I(0, θ_D/T) [numerisch] x = ℏω / (k_B·T) (dimensionslose Energie) Grenzfälle: T ≫ θ_D (klassisch): C_V^Debye → 3R (Dulong-Petit, exakt) T ≪ θ_D (quantum): C_V^Debye → (12π⁴/5)R·(T/θ_D)³ (Debye-T³-Gesetz) T = θ_D: C_V^Debye ≈ 2.85R (≈ 95% des klass. Wertes)

14. Frenkel-Abschneidung in der Flüssigkeit

In Flüssigkeiten haben Bindungen die endliche Lebensdauer τ_b. Transversale Phononen (TA) mit ω < ω_F werden diffusiv — sie speichern keine potenzielle Energie mehr:

Frenkel-Frequenz: ω_F = 1 / τ_b [rad/s] Frenkel-Temperatur: θ_F = ℏ / (k_B · τ_b) = 7.639 K·ps / τ_b[ps] Phonon-Beiträge der Flüssigkeit (1 LA-Zweig + 2 TA-Zweige): C_V^LA = 3R·(T/θ_D)³ · I(0, θ_D/T) ← kein Frenkel C_V^TA = 6R·(T/θ_D)³ · I(θ_F/T, θ_D/T) ← Frenkel-begrenzt C_V^akustisch(Flüssigk.) = C_V^LA + C_V^TA mit dem Debye-Integral: I(x_low, x_high) = ∫_{x_low}^{x_high} x⁴eˣ/(eˣ−1)² dx Grenzverhalten: τ_b → ∞ (Kristall): θ_F → 0 → C_V^TA → C_V^voll = Debye-Kristall τ_b → 0 (Gas): θ_F → ∞ → C_V^TA = 0 nur LA bleibt

Die Frenkel-Abschneidung erklärt, warum Flüssigkeiten eine kleinere Wärmekapazität haben als Kristalle bei gleicher Temperatur: TA-Phononen unterhalb ω_F fehlen. Die Wärmekapazität der Flüssigkeit liegt daher zwischen der des Gases (nur C_V^LA, 1 Zweig) und des Kristalls (alle 3 Zweige).

15. Phonon-Zustandssumme

Die Phononen sind Bosonen: jeder Mode kann beliebig oft besetzt sein. Die kanonische Zustandssumme eines bosonischen Felds ist das Produkt über alle Moden:

Z_Phonon = ∏_k 1 / (1 − exp(−ℏω_k / k_BT)) Freie Energie: F_Phonon = −k_BT · Σ_k ln(1 − exp(−ℏω_k / k_BT)) = Nk_BT · ∫_0^{ω_D} g(ω) · ln(1 − exp(−ℏω/k_BT)) dω Vergleich mit molekularer Zustandssumme (Maxwell-Boltzmann): Z_mol ∝ exp(−E_i / k_BT) ← diskrete Energieniveaus, N fixe Teilchen Z_Phonon = ∏_k 1/(1−exp(…)) ← kontinuierliches Spektrum, N_mode = f(V) Schlüsselunterschied: Phononen sind nicht konserviert (kein chemisches Potential μ=0). Anzahl der thermisch angeregten Phononen: ⟨N_ph⟩ = ∫_0^{ω_D} g(ω) / (exp(ℏω/k_BT)−1) dω ~ T³ für T ≪ θ_D

Die Zustandssumme Z_Phonon ist multiplikativ (nicht additiv): jeder Mode trägt unabhängig bei. Das macht das Phononenfeld zu einem variierenden Feld — im Gegensatz zu den N fixierten molekularen Freiheitsgraden.

16. Literatur (Phonon-Erweiterungen)^{(KI-generiert, Halluzinationen und Fehler möglich)}

[P1] Debye, P.: Zur Theorie der spezifischen Wärmen. Ann. Phys. 39 (1912) 789. — Grundlage des Debye-Modells; g(ω) ∝ ω²; T³-Gesetz.
[P2] Frenkel, J.: Kinetic Theory of Liquids. Oxford UP, 1946. — Transversale Phononen in Flüssigkeiten; ω_min = 1/τ.
[P3] Bolmatov, D., Brazhkin, V. V. & Trachenko, K.: The phonon theory of liquid thermodynamics. Sci. Rep. 2 (2012) 421. doi:10.1038/srep00421
[P4] Trachenko, K. & Brazhkin, V. V.: Collective modes and thermodynamics of the liquid state. Rep. Prog. Phys. 79 (2016) 016502. doi:10.1088/0034-4885/79/1/016502
[P5] Alexander, S. & Orbach, R.: Density of states on fractals: fractons. J. Phys. Lett. 43 (1982) L625. doi:10.1051/jphyslet:019820043017062500 — Frakton-Konzept; spektrale Dimension d_s.

10. Literatur zu den Formeln^{(KI-generiert, Halluzinationen und Fehler möglich)}

[F1] Eyring, H.: Viscosity, plasticity, and diffusion as examples of absolute reaction rates. J. Chem. Phys. 4 (1936) 283. doi:10.1063/1.1749836 — Grundlage der Arrhenius-Viskosität und Sprungdiffusionstheorie.
[F2] Trouton, F.: On molecular latent heat. Philos. Mag. 18 (1884) 54 — Trouton'sche Regel ΔS_vap ≈ 10.5·R.
[F3] Angell, C. A.: Formation of glasses from liquids and biopolymers. Science 267 (1995) 1924. doi:10.1126/science.267.5206.1924 — Arrhenius vs. VFT-Viskosität, Glasbildnerklassifikation.
[F4] Chandler, D., Weeks, J. D. & Andersen, H. C.: Van der Waals picture of liquids, solids, and phase transformations. Science 220 (1983) 787. doi:10.1126/science.220.4599.787 — Perturbationstheorie für Flüssigkeiten; Grundlage für z(T)-Herleitung.
[F5] Lindemann, F. A.: Über die Berechnung molekularer Eigenfrequenzen. Physik. Z. 11 (1910) 609 — Lindemann-Kriterium für Schmelzpunkt; Grundlage für den Faktor δ in der ΔH_melt-Formel.
[F6] Trouton, F.: On molecular latent heat. Philos. Mag. 18 (1884) 54 — Trouton'sche Regel: ΔS_verd ≈ 87 J/(mol·K) für unassoziierte Flüssigkeiten; Grundlage der T_siede- und ΔH_verd-Formeln.
[F7] Chickos, J. S. & Acree, W. E.: Enthalpies of Vaporization of Organic and Organometallic Compounds 1880–2002. J. Phys. Chem. Ref. Data 32 (2003) 519. doi:10.1063/1.1529214 — Umfassende Sammlung von ΔH_vap-Werten; Validierung der Netzwerkmodell-Vorhersagen für organische Substanzen.