Wie berechnet man das Phasenraumvolumen des Liouville-Theorems

Berechnen Sie das 6N-dimensionale Phasenraumvolumen (q,p)^N erhalten durch Hamiltonscher Fluss (Joseph Liouville 1838). Idealgas, Synchrotron, BEC, Plasma — INFN-LNF DAFNE, LENS Florenz, Cineca Bologna.

Berechnungsparameter
War dieser Rechner nützlich?

Der Liouville-Theorem-Rechner bestimmt das Volumen im 6N-dimensionalen Phasenraum (q₁,...,q_3N, p₁,...,p_3N) eines N-Teilchen-Systems, gemäß dem fundamentalen Theorem der statistischen und klassischen Mechanik formuliert von Joseph Liouville (École Polytechnique Paris, 1838) in "Sur la théorie générale de la variation des constantes arbitraires" Journal de Mathématiques Pures et Appliquées 3:342. Das Theorem besagt, dass die Wahrscheinlichkeitsdichte ρ(q,p,t) im Phasenraum entlang der dynamischen Trajektorien des Systems invariant ist: dρ/dt = ∂ρ/∂t + {ρ,H} = 0. Mathematisches Fundament der statistischen Thermodynamik.

1. Liouville-Gleichung und Formel

Klassische Liouville-Gleichung: ∂ρ/∂t = -{ρ, H}. Äquivalent: div(ρ·v_phase) = 0 (Hamiltonscher Fluss inkompressibel in 6N dim). Phasenvolumen-Erhaltung: V_phase = ∫∫ d³ᴺq d³ᴺp invariant. Quantisierung: Volumen in h³-Zellen pro Teilchen. N_Zellen = V_position · V_momentum / h^(3N). Für Idealgas: V_momentum ≈ (2π m k_B T)^(3/2). S = k_B ln(N_Zellen) (Boltzmann 1877).

2. Typische Systeme und Parameter

  • Idealgas Ar NTP: N=6,02·10²³, V=22,4 L, T=300 K, λ_th=16 pm.
  • Elektronen in Cu: N=8,5·10²² /cm³, T=300 K, λ_th=4 nm → Fermi-Dirac.
  • Sonnenkorona-Plasma: T=10⁶ K, n=10⁹ /cm³, klassisch.
  • Synchrotron-Elektronenstrahl: γ=10⁴, ε=1 nm·rad (DAFNE Frascati 510 MeV).
  • BEC Rb-87: T=10⁻⁷ K → BEC (Cornell-Wieman Nobel 2001, LENS Florenz 1991).
  • Neutronen Neutronensterne: n=10⁵⁷ /m³ entartet.
  • Compton-Lichtstrahlen: Brilliance Liouville-eingeschränkt.
  • Atomwolken Stern-Gerlach: Phasenraum-Tomographie Wigner-Funktion.

3. Grenzen und Verallgemeinerungen des Theorems

Klassisches Liouville gilt für konservative Hamiltonsche Systeme. Grenzen: (a) Dissipative Systeme: ρ kontrahiert → seltsame Attraktoren, Fraktale D<2N. (b) Quantensysteme: von-Neumann-Gleichung iℏ∂ρ̂/∂t = [Ĥ,ρ̂]. (c) Nichtlineare chaotische Systeme: Trajektorien divergieren exponentiell (Lyapunov), aber Volumen erhalten (KAM 1954-1963). (d) Offene Systeme: Lindblad master equation. (e) Relativistische Systeme: Synge 1957. (f) BEC: nichtlineare Gross-Pitaevskii-Gleichung.

4. Numerische Beispiele

Beispiel 1 — Idealgas Ar NTP. N=6,022×10²³, V=22,4 L, T=300 K. λ_th = 1,60×10⁻¹¹ m. V_phase pro Teilchen = 4,5×10⁻⁷⁴ (J·s)³. Planck-Zellen = 8,2×10⁶ pro Teilchen. S/N = 18,5 k_B = 154,2 J/(K·mol) (vs NIST 154,8).

Beispiel 2 — Rb-87 BEC bei 100 nK. λ_th = 1,15 μm. n·λ³ = 1,5 → BEC. Cornell-Wieman 1995 Nobel 2001.

Beispiel 3 — DAFNE 510 MeV Synchrotron. γ=1000. Emittanz ε = 1 nm·rad durch Liouville invariant. Kein Laser-Cooling für Elektronen.

5. Praktische Anwendungen des Theorems

  • Statistische Thermodynamik: Basis Boltzmann S=k_B ln Ω, Sackur-Tetrode 1911-1912 Idealgas-Entropie, kanonische Zustandssumme Z.
  • Beschleunigerphysik: invariante Strahl-Emittanz, Luminositätserhaltung INFN-LNF DAFNE 510 MeV + KLOE-2, Bruno-Touschek-Effekt 1956 (ADONE Frascati 1969).
  • Kühltechniken: stochastic cooling Van der Meer Nobel 1984 CERN UA1, Synchrotronstrahlungsdämpfung Leptonen.
  • Kalte Atome LENS Florenz 1991: Doppler/Sisyphus-Laserkühlung → μK, Verdunstungskühlung → BEC, MPQ Garching + Universität Innsbruck + ETH Zürich.
  • Plasmaeinschluss: ITER + DTT Frascati 2028, ENEA Tokamak, Wendelstein 7-X Greifswald, ASDEX Upgrade Garching.
  • Stellardynamik: INAF Bologna + Triest 1898 + Capodimonte 1812 + MPA Garching N-Körper mit Liouville → f(q,p,t).
  • Quantenchaos: Berry-Tabor, Bohigas-Giannoni-Schmit, GOE/GUE/GSE.
  • Molekulardynamik: GROMACS, LAMMPS, NAMD, AMBER auf HPC Cineca Bologna 1969 Leonardo 250 PFlops + INFN-CNAF + Jülich Supercomputing + LRZ Garching.

6. Vergleich mit anderen Erhaltungssätzen

Liouville 1838: ρ + V_phase Erhaltung Hamiltonsch. Boltzmann H-Theorem 1872: H abnehmend (Entropie wachsend) für Nicht-Gleichgewicht. Von Neumann 1927: iℏ∂ρ̂/∂t = [Ĥ,ρ̂] erhält Tr(ρ̂²). Lindblad 1976: master equation dissipative Systeme. Wigner-Funktion W(q,p): Phasenraum-Quantenmechanik-Darstellung. Husimi-Q-Funktion: geglättete Wigner. P-Glauber-Sudarshan: diagonale Darstellung. Jeans 1915: kollisionsloses Liouville Stellardynamik. Vlasov 1938: Liouville + Mean-Field-Coulomb Plasma. Boltzmann-Transportgleichung: Liouville + Kollisionsintegral.

7. Moderne Entwicklungen und Anwendungen

Die modernen Anwendungen umfassen: Symplektische Integratoren (Verlet, Leapfrog, Yoshida 1990): erhalten Liouville → MD-Langzeitstabilität. KAM-Theorie 1954-1963 (Kolmogorov-Arnold-Moser): Persistenz invarianter Tori, Sonnensystem-Stabilität (Laskar 1989 chaotisch). Ergodizität + Mischung: ergodische Hypothese + Sinai 1970. Quantenchaos: Bohigas-Giannoni-Schmit 1984, Berry-Tabor 1977. Holomorphe Ergodizität: Ratner 1991 (Fields-Medaille 2014 Mirzakhani). Phasenraum-Tomographie: LIF kalte Atome. Kinetische Plasmatheorie: Vlasov-Poisson, Vlasov-Maxwell ITER + DTT Frascati 2028. Software: GROMACS, LAMMPS, NAMD, AMBER, Quantum ESPRESSO (SISSA Triest 1978 Stefano Baroni), Cantera, Modelica. HPC: Cineca Bologna 1969 Leonardo 250 PFlops + INFN-CNAF + Jülich Supercomputing + LRZ Garching + DKRZ Hamburg. Zukünftige Experimente: EIC Brookhaven 2030, SKA 2030.

Aktualisiert im Juli 2026

Häufig gestellte Fragen

Was ist das Liouville-Theorem?

Das Liouville-Theorem (Joseph Liouville École Polytechnique Paris 1838, J Math 3:342) besagt, dass die Wahrscheinlichkeitsdichte ρ(q,p,t) im Phasenraum entlang der dynamischen Trajektorien eines Hamiltonschen Systems invariant ist: dρ/dt = ∂ρ/∂t + {ρ,H} = 0. Äquivalent ist das 6N-dimensionale Phasenvolumen durch den Hamiltonschen Fluss erhalten. Fundament der statistischen Mechanik.

Wie wird es mathematisch formuliert?

Klassische Liouville-Gleichung: ∂ρ/∂t = -{ρ,H} = -∑(∂ρ/∂q_i ∂H/∂p_i - ∂ρ/∂p_i ∂H/∂q_i), wobei {·,·} Poisson-Klammern, H Hamiltonian. Quantum (von Neumann 1927): iℏ ∂ρ̂/∂t = [Ĥ,ρ̂]. Äquivalent: div(ρ·v_phase) = 0 (Inkompressibilität des Flusses). Volumen V_phase = ∫∫ d³ᴺq d³ᴺp invariant.

Warum ist es für die Thermodynamik fundamental?

Liouville ist die Basis der statistischen Mechanik: ohne es können mikro-kanonisches/kanonisches/groß-kanonisches Ensemble nicht definiert werden. Phasenvolumen quantisiert in h³-Zellen (Sackur-Tetrode 1911-1912) verbindet klassische Dynamik mit Entropie S=k_B ln Ω (Boltzmann 1877). Ohne Liouville wäre die Thermodynamik von Idealgas, BEC, Fermi-Gas, Plasmen inkonsistent.

Wann gilt Liouville NICHT?

Gilt nicht für: dissipative Systeme (Reibung, Viskosität, Strahlung → ρ kontrahiert, fraktale seltsame Attraktoren); offene Systeme (Austausch mit Umgebung → Lindblad master equation); Quantenmischzustände (klassisch von Neumann). Gilt streng nur für konservative Hamiltonians. KAM-Theorie 1954-1963 erweitert auf kleine Störungen.

Was sagt Liouville für die Beschleunigerphysik?

Für geladene Teilchenstrahlen impliziert Liouville Erhaltung der Emittanz ε = π·σ_x·σ_x' = invariant unter Hamiltonian. Konsequenz: man kann einen Strahl mit konservierten Ladungen nicht kühlen (kein Proton-Laser-Cooling). Workarounds: stochastic cooling (Van der Meer Nobel 1984), Synchrotron-Strahlungsdämpfung (nur Leptonen e+e-), Cooling für verdunstende Ionen oder Myonen (MUCOOL Fermilab).

Wie wird das Phasenvolumen quantisiert?

Heisenberg 1927 + Sackur-Tetrode 1911-1912: das elementare Phasenraumvolumen ist h^3 pro Teilchen (h Plancksche Konstante). Zellenanzahl = V_position · V_momentum / h^(3N). Für Idealgas bei T: V_momentum ≈ (2π m k_B T)^(3/2). Boltzmann 1877 Entropie-Verbindung: S = k_B ln(N_Zellen/N!) (Gibbs-Paradoxon gelöst).

Welche praktischen Anwendungen gibt es?

Statistische Thermodynamik (Sackur-Tetrode), Beschleuniger-Emittanz (DAFNE INFN-LNF + ADONE 1969, Touschek-Effekt), stochastic cooling (Van der Meer Nobel 1984), kalte Atome LENS Florenz 1991 (BEC), ITER + DTT Frascati 2028, Stellardynamik INAF, MD GROMACS HPC Cineca Leonardo 250 PFlops + Jülich + LRZ.

Kommentare (0)

Noch keine Kommentare. Sei der Erste!