Statistical mechanics · two-level equilibrium

Boltzmann Distribution Calculator

Translate an energy gap and a thermodynamic temperature into an equilibrium state population. This independent physics workspace keeps energy conversion, degeneracy, the partition function, entropy, heat capacity, and the temperature response on one auditable page.

eV · J · kJ/mol · cm⁻¹ · K Degenerate energy levels Live thermal response Two-state thermodynamics

Define the two-level system

The upper level is assumed to lie ΔE above the lower level. Temperature is converted to kelvins and every energy unit is converted to joules per particle before solving.

Celsius and Fahrenheit are converted to thermodynamic temperature.
The model requires a positive upper-minus-lower energy difference.
Number or statistical weight of microstates at the lower energy.
Number or statistical weight of microstates at the upper energy.
The calculator solves the positive temperature required to reach this target, when finite.

Thermal population ledger

Updates with every input
Thermal energy kBT
Boltzmann factor
Population ratio N₁/N₀
Upper population p₁ Normalized across both levels
Lower population p₀ p₀ + p₁ = 100%
Relative partition function Z Lower energy chosen as zero
Average excitation energy Per mole above E₀
Molar entropy S Two-level equilibrium entropy
Molar heat capacity C Two-level thermal response
Dimensionless gap ΔE / kBT
High-temperature p₁ limit g₁ / (g₀ + g₁)
Target-population temperature

State occupancy and thermal response

The left panel connects the entered energy separation to the normalized populations. The right curve uses the universal thermal ratio kBT/ΔE, so the current marker moves without hiding the role of degeneracy.

Lower state Upper state / current point High-temperature limit

How the same energy gap responds when temperature changes

Every row preserves the entered gap and degeneracies. Only absolute temperature changes, making the exponential response, state population, and heat-capacity peak directly comparable.

Scenario Temperature (K) kBT (eV) ΔE/kBT Boltzmann factor Upper population C (J/mol·K)

From energy units to normalized populations

  1. Convert the entered temperature to kelvins and the entered gap to joules per particle.
  2. Calculate thermal energy kBT and the dimensionless separation x = ΔE/kBT.
  3. Evaluate the raw Boltzmann factor B = exp(−x), which describes energy suppression before degeneracy.
  4. Multiply B by the degeneracy ratio g₁/g₀ to obtain N₁/N₀.
  5. Normalize both statistical weights through Z = g₀ + g₁B to obtain p₀ and p₁.
  6. Use those normalized populations to calculate average excitation energy, entropy, and heat capacity.
kB = 1.380 649 × 10⁻²³ J·K⁻¹ (exact)
B = exp[−(E₁ − E₀)/(kBT)] = exp(−ΔE/kBT)
N₁/N₀ = (g₁/g₀)B
Z = g₀ + g₁B;   p₀ = g₀/Z;   p₁ = g₁B/Z
U − E₀ = p₁ΔE
S/R = ln Z + (ΔE/kBT)p₁
C/R = (ΔE/kBT)²p₁(1 − p₁)

Current values inserted step by step

Energy-gap cross-check

Meaning and unit of every term

SymbolMeaningUnit
kBBoltzmann constantJ/K
TThermodynamic temperatureK
E₀, E₁Lower and upper state energiesJ per particle
ΔEPositive energy gap E₁ − E₀J, eV, kJ/mol, cm⁻¹, or K equivalent
g₀, g₁Lower and upper state degeneraciesdimensionless
BRaw Boltzmann factor exp(−ΔE/kBT)dimensionless
ZPartition function relative to E₀ = 0dimensionless
p₀, p₁Normalized state populationsfraction or percent
RMolar gas constant NAkBJ/(mol·K)

Degeneracy competes with the energy penalty

The exponential term always penalizes the higher energy at positive temperature. Degeneracy pushes in the opposite direction by providing more upper-level microstates. That is why the raw Boltzmann factor and the normalized upper-state population are not interchangeable.

  • When ΔE ≫ kBT, the upper level is strongly suppressed.
  • When ΔE ≈ kBT, both energy and degeneracy materially affect population.
  • When kBT ≫ ΔE, energy becomes unimportant and degeneracy sets the limiting shares.
  • A population ratio can exceed one when upper-level degeneracy is sufficiently large.

What this page deliberately does not assume

This is an equilibrium two-level canonical model. It does not calculate a continuous molecular speed distribution, non-equilibrium relaxation, transition rates, inverted populations, chemical reactions, quantum coherence, level broadening, interactions, or a multi-level molecular partition function.

Enter degeneracies and energy differences that belong to the same physical model. Spectroscopic term symbols, vibrational ladders, rotational levels, and electronic states may require many more levels and additional selection or symmetry rules.

Questions about factors, populations, and temperature

What does the Boltzmann factor mean?

The Boltzmann factor exp(-ΔE/kBT) measures the thermal suppression associated with an energy gap ΔE at absolute temperature T. A smaller gap or a higher temperature moves the factor closer to one.

Why does degeneracy change the state population?

Degeneracy counts how many microstates share the same energy. The population weight of a level is its degeneracy multiplied by its Boltzmann factor, so a highly degenerate upper level can hold more population than the energy gap alone suggests.

Can the upper-state population exceed 50%?

Yes. If the upper level has greater degeneracy, its high-temperature limit is g1/(g0+g1), which can exceed 50%. For equal degeneracies, the positive-temperature upper-state population approaches 50% from below.

Does this calculate a Maxwell-Boltzmann speed distribution?

No. This page models equilibrium occupation of two discrete energy levels. A Maxwell-Boltzmann speed distribution is a continuous distribution that also depends on particle mass.

Why must temperature be expressed on an absolute scale?

The exponential uses thermodynamic temperature. Celsius and Fahrenheit entries are converted to kelvins before any calculation, and temperatures at or below absolute zero are rejected.

What is the high-temperature population limit?

As temperature grows, the energy penalty becomes negligible and population is divided according to degeneracy. The upper-state limit is therefore g1/(g0+g1), not automatically 50%.

Constants and physical interpretation

The calculation uses the exact SI defining values kB = 1.380 649 × 10⁻²³ J/K, h = 6.626 070 15 × 10⁻³⁴ J·s, c = 299 792 458 m/s, e = 1.602 176 634 × 10⁻¹⁹ C, and NA = 6.022 140 76 × 10²³ mol⁻¹.