Cómo calcular el volumen del espacio de fases del teorema de Liouville

Calcula el volumen 6N-dimensional del espacio de fases (q,p)^N preservado por flujo Hamiltoniano (Joseph Liouville 1838). Gas ideal, sincrotrón, BEC, plasma — INFN-LNF DAFNE, LENS Florencia, Cineca Bologna.

Parámetros de cálculo
¿Te ha sido útil esta calculadora?

La calculadora del teorema de Liouville determina el volumen en el espacio de fases 6N-dimensional (q₁,...,q_3N, p₁,...,p_3N) de un sistema de N partículas, según el teorema fundamental de la mecánica estadística y clásica formulado por Joseph Liouville (École Polytechnique París, 1838) en "Sur la théorie générale de la variation des constantes arbitraires" Journal de Mathématiques Pures et Appliquées 3:342. El teorema afirma que la densidad de probabilidad ρ(q,p,t) en el espacio de fases es invariante a lo largo de las trayectorias dinámicas del sistema: dρ/dt = ∂ρ/∂t + {ρ,H} = 0, donde {·,·} son corchetes de Poisson y H el Hamiltoniano. Fundamento matemático de la termodinámica estadística.

1. Ecuación de Liouville y fórmula

La ecuación de Liouville clásica: ∂ρ/∂t = -{ρ, H}. Equivalentemente: div(ρ·v_phase) = 0 (flujo Hamiltoniano incompresible en 6N dim). Conservación del volumen de fase: V_phase = ∫∫ d³ᴺq d³ᴺp invariante. Cuantización: volumen en celdas h³ por partícula (h Planck). N_celdas = V_position · V_momentum / h^(3N). Para gas ideal: V_momentum ≈ (2π m k_B T)^(3/2). Conexión entropía: S = k_B ln(N_celdas) (Boltzmann 1877).

2. Sistemas típicos y parámetros

  • Gas ideal Ar CNTP: N=6,02·10²³, V=22,4 L, T=300 K, λ_th=16 pm.
  • Electrones en Cu: N=8,5·10²² /cm³, T=300 K, λ_th=4 nm → Fermi-Dirac.
  • Plasma corona solar: T=10⁶ K, n=10⁹ /cm³, clásico.
  • Beam sincrotrón: γ=10⁴, ε=1 nm·rad (DAFNE Frascati 510 MeV).
  • BEC Rb-87: T=10⁻⁷ K → BEC (Cornell-Wieman Nobel 2001, LENS Florencia 1991).
  • Neutrones estrellas neutrón: n=10⁵⁷ /m³ degenerado.
  • Compton light beams: brilliance Liouville-constrained.
  • Atomic clouds Stern-Gerlach: tomografía espacio fases Wigner.

3. Límites y generalizaciones del teorema

Liouville clásico vale para sistemas Hamiltonianos conservativos. Límites: (a) Sistemas disipativos: ρ contracts → atractores extraños, fractales D<2N. (b) Sistemas cuánticos: ecuación von Neumann iℏ∂ρ̂/∂t = [Ĥ,ρ̂]. (c) Sistemas no-lineales caóticos: trayectorias divergen exponencialmente (Lyapunov), pero volumen conservado (KAM 1954-1963). (d) Sistemas abiertos: master equation Lindblad. (e) Sistemas relativistas: Synge 1957. (f) BEC: ecuación Gross-Pitaevskii no-lineal.

4. Ejemplos numéricos resueltos

Ejemplo 1 — Gas ideal Ar CNTP. N=6,022×10²³, V=22,4 L, T=300 K. λ_th = 1,60×10⁻¹¹ m. V_phase por partícula = 4,5×10⁻⁷⁴ (J·s)³. Celdas Planck = 8,2×10⁶ por partícula. S/N = 18,5 k_B = 154,2 J/(K·mol) (vs NIST 154,8).

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

Ejemplo 3 — Sincrotrón DAFNE 510 MeV. γ=1000. Emittancia ε = 1 nm·rad invariante por Liouville. No laser cooling para electrones.

5. Aplicaciones prácticas del teorema

  • Termodinámica estadística: base Boltzmann S=k_B ln Ω, Sackur-Tetrode 1911-1912 entropía gas ideal, función partición Z.
  • Accelerator physics: emittancia beam invariante, conservación luminosidad INFN-LNF DAFNE 510 MeV + KLOE-2, Bruno Touschek effect 1956 (ADONE Frascati 1969).
  • Cooling techniques: stochastic cooling Van der Meer Nobel 1984 CERN UA1, sincrotrón radiation damping leptones.
  • Cold atoms LENS Florencia 1991: laser cooling Doppler/Sisyphus → μK, evaporative cooling → BEC.
  • Plasma confinement: ITER + DTT Frascati 2028, ENEA Tokamak.
  • Stellar dynamics: INAF Bologna + Trieste 1898 + Capodimonte 1812 N-body con Liouville → f(q,p,t).
  • Quantum chaos: Berry-Tabor, Bohigas-Giannoni-Schmit, GOE/GUE/GSE.
  • Molecular dynamics: GROMACS, LAMMPS, NAMD, AMBER en HPC Cineca Bologna 1969 Leonardo 250 PFlops + INFN-CNAF + BSC Barcelona.

6. Comparación con otras conservaciones

Liouville 1838: ρ + V_phase conservativos. Boltzmann H-theorem 1872: H decrece (entropía crece) para no-equilibrio. Von Neumann 1927: iℏ∂ρ̂/∂t = [Ĥ,ρ̂] preserva Tr(ρ̂²). Lindblad 1976: master equation sistemas disipativos. Wigner function W(q,p): representación phase space cuántica. Husimi Q-function: Wigner suavizado. P-Glauber-Sudarshan: representación diagonal. Jeans 1915: collisionless Liouville stellar dynamics. Vlasov 1938: Liouville + mean-field Coulomb plasma. Boltzmann transport equation: Liouville + integral colisiones.

7. Desarrollos y aplicaciones modernas

Las aplicaciones modernas incluyen: Integradores simplécticos (Verlet, leapfrog, Yoshida 1990): preservan Liouville → MD estabilidad. KAM theory 1954-1963 (Kolmogorov-Arnold-Moser): persistencia tori invariantes, estabilidad sistema solar (Laskar 1989 caótico). Ergodicidad + mezcla: hipótesis ergódica + Sinai 1970. Quantum chaos: Bohigas-Giannoni-Schmit 1984, Berry-Tabor 1977. Ergodicidad holomórfica: Ratner 1991 (Fields medal 2014 Mirzakhani). Phase space tomography: LIF cold atoms. Kinetic theory plasmas: Vlasov-Poisson, Vlasov-Maxwell ITER + DTT Frascati 2028. Software: GROMACS, LAMMPS, NAMD, AMBER, Quantum ESPRESSO (SISSA Trieste 1978 Stefano Baroni), Cantera, Modelica. HPC: Cineca Bologna 1969 Leonardo 250 PFlops + INFN-CNAF + BSC Barcelona + LRZ Garching. Experimentos futuros: EIC Brookhaven 2030, SKA 2030.

Actualizado a julio de 2026

Preguntas frecuentes

¿Qué es el teorema de Liouville?

El teorema de Liouville (Joseph Liouville École Polytechnique París 1838, J Math 3:342) afirma que la densidad de probabilidad ρ(q,p,t) en el espacio de fases es invariante a lo largo de trayectorias dinámicas de un sistema Hamiltoniano: dρ/dt = ∂ρ/∂t + {ρ,H} = 0. Equivalentemente, el volumen de fase 6N-dimensional es preservado por flujo Hamiltoniano. Fundamento de la mecánica estadística.

¿Cómo se formula matemáticamente?

Ecuación de Liouville clásica: ∂ρ/∂t = -{ρ,H} = -∑(∂ρ/∂q_i ∂H/∂p_i - ∂ρ/∂p_i ∂H/∂q_i), donde {·,·} corchetes de Poisson, H Hamiltoniano. Cuántica (von Neumann 1927): iℏ ∂ρ̂/∂t = [Ĥ,ρ̂]. Equivalentemente: div(ρ·v_phase) = 0 (flujo incompresible). Volumen V_phase = ∫∫ d³ᴺq d³ᴺp invariante.

¿Por qué es fundamental para la termodinámica?

Liouville es la base de la mecánica estadística: sin él no se pueden definir ensembles micro-canónico/canónico/gran-canónico. Volumen de fase cuantizado en celdas h³ (Sackur-Tetrode 1911-1912) conecta dinámica clásica a entropía S=k_B ln Ω (Boltzmann 1877). Sin Liouville, termodinámica gas ideal, BEC, Fermi-gas, plasmas inconsistente.

¿Cuándo NO vale Liouville?

No vale para: sistemas disipativos (fricción, viscosidad, radiación → ρ contracts, atractores extraños fractales); sistemas abiertos (intercambio con entorno → Lindblad master equation); sistemas cuánticos mezclas (von Neumann clásica). Vale rigurosamente solo para Hamiltonianos conservativos. KAM theory 1954-1963 extiende a perturbaciones pequeñas.

¿Qué dice Liouville para accelerator physics?

Para beam de partículas cargadas, Liouville implica conservación de la emitancia ε = π·σ_x·σ_x' = invariante bajo Hamiltoniano. Consecuencia: no se puede enfriar un beam con cargas conservadas (no laser cooling para protones). Workarounds: stochastic cooling (Van der Meer Nobel 1984), sincrotrón radiation damping (sólo leptones e+e-), cooling iones evaporativos o muones (MUCOOL Fermilab).

¿Cómo se cuantiza el volumen de fase?

Heisenberg 1927 + Sackur-Tetrode 1911-1912: el volumen elemental en espacio de fases es h^3 por partícula (h constante de Planck). Número de celdas = V_position · V_momentum / h^(3N). Para gas ideal a temperatura T: V_momentum ≈ (2π m k_B T)^(3/2). Conexión entropía Boltzmann 1877: S = k_B ln(N_celdas/N!) (paradoja Gibbs resuelta).

¿Cuáles son las aplicaciones prácticas?

Termodinámica estadística (Sackur-Tetrode), accelerator emitancia (DAFNE INFN-LNF + ADONE 1969, Touschek), stochastic cooling (Van der Meer Nobel 1984), cold atoms LENS Florencia 1991 (BEC), confinement ITER + DTT Frascati 2028, stellar dynamics INAF, MD GROMACS HPC Cineca Leonardo 250 PFlops + BSC Barcelona.

Comentarios (0)

Sin comentarios aún. ¡Sé el primero!