Boltzmann Statistical Mechanics & Energy Distribution Simulation

Visualizing hard-sphere thermalization, collision-driven kinetic redistribution, and the natural emergence of exponential population decay $N(E) \propto e^{-E/k_B T}$.

Total Ensemble Particles N = 180 Conserved System Size
Mean Kinetic Energy ⟨E⟩ 1.00 E₀ Equipartition Baseline
Low-Energy Fraction (E < ⟨E⟩) ~63.2% Dominant Base Population
High-Energy Tail (E > 2⟨E⟩) ~13.5% Exponentially Suppressed
40 a.u.
1.0x

1. 2D Hard-Sphere Molecular Chamber Real-Time Kinetic Motion

Low Energy ($E < \langle E \rangle$) Moderate Energy ($\approx \langle E \rangle$) High Energy Tail ($E > 2\langle E \rangle$)

2. Live Energy Population Distribution Continuous Binned Counts vs. Theoretical Curve

Observed Bins $N(E)$ Boltzmann Curve: $N(E) \propto e^{-E/k_B T}$

Physical Principles of Ludwig Boltzmann’s Distribution (1877)

• Probability Density: P(E) = (1 / k_B T) · e^(-E / k_B T)
• Boltzmann Multiplicity: S = k_B · ln(Ω)
• Conservation Law: Total Energy E_total = ∑ (1/2 · m · v_i²) = Constant