Bose–Einstein Distribution Calculator

What this calculator tells you

This calculator focuses on one of the central ideas of quantum statistics: bosons do not have to avoid one another. Unlike fermions, they can pile into the same quantum state, and at low enough temperature that tendency becomes dramatic. The page gives you two related outputs. First, it computes the occupation number for a single-particle state of energy E. That output tells you the expected number of bosons in that state when the system is at temperature T and chemical potential μ. Second, it estimates the critical temperature T_c for Bose–Einstein condensation in an ideal gas with the supplied particle mass and number density. Read together, those values let you ask a physically meaningful question: is this state strongly occupied, and is the gas cold enough for a condensate to be expected in the ideal model?

The calculator is intentionally compact, but the physics behind the inputs is easy to lose if you treat the form like a black box. A Bose–Einstein occupation number is not a generic score. It is the expected population of a specific quantum state. Likewise, the critical temperature is not a universal constant for a given atom; it changes with mass and with how densely the particles are packed. That is why the explanation here leans on interpretation as much as arithmetic. If you understand what each input means, the output becomes much easier to trust.

What each input means in plain language

Energy E (J) is the single-particle energy of the state whose occupation you want to evaluate. In many textbook problems this is one energy level in a ladder of allowed states. In laboratory contexts, it may represent a mode, trap level, or chosen reference state. Temperature T (K) is the absolute temperature of the boson gas. Lower temperature generally increases low-energy occupation because thermal agitation becomes less effective at spreading particles into excited states.

Chemical potential μ (J) controls how the system distributes particles among available states. For an ideal Bose gas, μ approaches the lowest energy level from below as the gas nears condensation. That “from below” language matters: if you set μ too high for the state you are evaluating, the denominator in the Bose–Einstein formula stops being physical. Particle mass m (kg) and number density n (m⁻³) are used only for the critical-temperature estimate. Heavier particles tend to lower the critical temperature, while higher density tends to raise it because packing more bosons into the same volume makes macroscopic occupation easier.

The units also matter more than they may seem at first glance. Energy and chemical potential must both be entered in joules, temperature in kelvin, mass in kilograms, and density in particles per cubic meter. If your source data is in electronvolts, microkelvin, atomic mass units, or particles per cubic centimeter, convert before entering values. A correct formula fed the wrong units will still give the wrong physical answer.

The core Bose–Einstein formula

The occupation number for a bosonic state of energy E is computed from the standard Bose–Einstein expression below. The result is dimensionless because it represents an expected count in one state, not an energy or temperature.

Here k is Boltzmann’s constant. The structure of the formula already tells you how the answer behaves. If T decreases, then the thermal scale kT becomes smaller. If μ rises toward E from below, then the numerator of the exponent, E − μ, becomes smaller. Both changes reduce the exponential penalty and increase the occupation number. In the limit where μ gets very close to the state energy from below, occupation can become very large. That is the mathematical fingerprint of bosons bunching into low-energy states.

The critical temperature estimate comes from the ideal-gas Bose condensation condition. This calculator uses the common three-dimensional expression with the numerical factor ζ(3/2) ≈ 2.612.

This is why the second result depends on mass and density, but not on the chosen state energy. The occupation number is about one state. The critical temperature is about the whole ideal gas.

Why the calculator rejects some inputs

The most important physical restriction on the form is the relationship between μ and E. For positive temperature, you need μ < E for the chosen state. If you violate that condition, the denominator of the Bose–Einstein formula becomes zero or negative, which is why the calculator reports a nonphysical situation instead of a number. This is not a software quirk. It is the model telling you the specified state and chemical potential cannot be combined that way in equilibrium.

The form also expects T > 0, m > 0, and n > 0. A zero or negative absolute temperature is outside the intended ideal-gas interpretation here, and nonpositive mass or density makes the condensation formula meaningless. If you are working with photons, remember that the equilibrium chemical potential is typically zero. If you are working with ultracold atoms, μ may be negative relative to your energy zero or may approach the ground-state energy from below as the condensate threshold is approached.

A worked example with realistic scales

Suppose you want a quick ideal-gas estimate for an ultracold rubidium-87 cloud. Use an energy level E = 2.0 × 10⁻³⁰ J, a temperature T = 1.5 × 10⁻⁶ K, a chemical potential μ = 1.6 × 10⁻³⁰ J, a particle mass m = 1.443 × 10⁻²⁵ kg, and a number density n = 1.0 × 10²¹ m⁻³. The gap E − μ is positive but small, which already hints that occupation may be substantial. Meanwhile, the number density is high enough that the critical temperature will be in the microkelvin range rather than effectively zero.

When those values are inserted into the occupation formula, the exponent is small, so the denominator is much less than one and the occupation number becomes much larger than unity. Using the calculator, the occupation number comes out on the order of 10¹ to 10², meaning that this particular state is heavily occupied compared with a classical dilute-gas intuition. The critical temperature estimate lands near 1.8 × 10⁻⁶ K. Because the chosen temperature in this example is slightly below that threshold, the result message reports that a condensate is expected in the ideal model.

That worked example illustrates a useful habit: interpret both outputs together. A large occupation number alone does not mean the entire gas is condensed. It means the chosen state is strongly populated. The comparison between T and T_c is what tells you whether the gas as a whole sits below the ideal condensation threshold.

How to read the result without over-interpreting it

If the occupation number is very small, the chosen state is not strongly populated under the supplied conditions. That commonly happens when the state energy is much larger than the chemical potential, when the temperature is relatively high, or both. If the occupation number is moderate or large, bosonic bunching is becoming important for that state. If it grows extremely large, you are usually looking at a state very close to the chemical potential in a cold gas, which is exactly where Bose-enhanced occupation shows up most clearly.

The critical temperature result is a threshold estimate, not a guarantee about a real apparatus. In the ideal homogeneous-gas model, T < T_c means a condensate is expected. In real experiments, finite trap geometry, interactions, dimensionality, and nonequilibrium effects can shift the transition or blur it. So the right interpretation is: the calculator tells you what the ideal model predicts for the supplied density and mass, and it gives you a useful first check against experimental or textbook scales.

How this page’s formulas fit the broader mathematics

At a high level, the calculator still follows the same structure as many other scientific tools: inputs go in, a function acts on them, and a result comes out. The general notation below is still useful because it reminds you that even when the physics is specific, the software is applying a defined mapping from variables to outputs.

Quantum-statistical models also often combine contributions from many states. That broader idea is captured by the summation form below. While this particular calculator reports a single-state occupation number and an overall condensation threshold, the same summation mindset appears whenever you build total particle counts or thermodynamic averages from a whole spectrum of levels.

Assumptions and limits to keep in mind

This calculator uses the ideal Bose gas picture. That means interactions between particles are neglected, the system is treated as being in thermal equilibrium, and the critical-temperature estimate is the standard textbook form for a three-dimensional gas. Those assumptions are excellent for learning, for rough scaling checks, and for many comparison tasks. They are not the last word for real trapped atomic gases, strongly interacting systems, lower-dimensional condensates, or driven photonic platforms.

Another practical limit is numerical scale. Because kT is tiny in SI units at ultracold temperatures, the exponent in the occupation formula can become extremely large when E − μ is not small. In that regime the occupation number is effectively zero for practical purposes, and the script safely reports a tiny value rather than attempting an unstable calculation. On the opposite side, as μ approaches E from below, the occupation number rises sharply, which is physically meaningful but also a sign that you should pay close attention to how the energy reference was chosen.

As a final check, ask yourself two quick questions after every run. First, are the units coherent? Second, does the trend make sense? If you lower the temperature or raise the chemical potential toward the state energy from below, the occupation number should increase. If you increase density or use lighter bosons, the critical temperature should increase. When those trends line up with intuition, the result is usually on firm footing.

Optional mini-game: learn the intuition by tuning a condensate

If you want a faster, more tactile feel for the same ideas, the mini-game below turns the formula into a timing-and-tuning challenge. You steer the chemical potential line and cool the trap while energy packets stream toward the condensate gate. The best scores come from the same logic as the calculator: keep the gas cold, place μ just below the incoming state energy, and avoid crossing into the nonphysical region where μ ≥ E. It is optional, separate from the calculator result, and meant as a memory aid rather than a substitute for the math.