Bose–Einstein Condensation Critical Temperature Calculator

Bose–Einstein Condensation in Harmonic Traps

Bose–Einstein condensation, usually shortened to BEC, occurs when a gas of bosonic particles is cooled so far that a macroscopic number of atoms collects in the lowest available quantum state. Instead of behaving like a classical cloud of independent particles, the gas begins to act as a coherent quantum system with a shared matter-wave character. In laboratory settings this transition is most often discussed for dilute atomic gases held in magnetic or optical traps. The temperature scale is extremely low, typically in the nanokelvin to microkelvin range, so even a rough estimate of the transition temperature is useful when planning an experiment or checking whether a measured cloud is close to degeneracy.

This calculator is designed for the standard ideal-gas estimate of the critical temperature in a three-dimensional harmonic trap. That is the model used in many first-pass calculations for ultracold bosonic atoms such as rubidium-87, sodium-23, or lithium-7. You enter the total atom number, the atomic mass, and the three trap frequencies. The calculator then converts those values into SI units, forms the geometric mean trap frequency, and evaluates the usual expression for the critical temperature T_c. The result is shown in both nanokelvin and microkelvin so it is easy to compare with common experimental temperature scales.

The page is intentionally practical. It does not try to replace a full many-body treatment, but it does give a reliable baseline estimate for weakly interacting trapped Bose gases. If you are learning the topic, the sections below explain what each input means, how the formula is built, and how to interpret the output. If you already know the physics, you can jump straight to the form and compute a value in a few seconds.

Introduction

In a harmonic trap, the external potential rises quadratically away from the center. That shape is a good approximation for many magnetic and optical confinement methods used in ultracold atom experiments. Because the trap is characterized by three oscillation frequencies, one along each spatial axis, the energy level spacing depends on how tightly the atoms are confined in the x, y, and z directions. Tighter confinement means larger frequencies, larger level spacing, and generally a higher critical temperature for a fixed atom number.

The transition to Bose–Einstein condensation is not defined by a sudden change in chemistry or structure. Instead, it is a statistical transition in the occupation of quantum states. Above the critical temperature, atoms are distributed over many excited states. As the gas is cooled, the excited states can no longer hold all particles while maintaining the Bose distribution, and the excess population accumulates in the ground state. That accumulation is what we call condensation. In the ideal trapped-gas model, the critical temperature depends mainly on two things: how many bosons are present and how dense the single-particle energy levels are, which is controlled by the trap frequencies.

Although the calculator mentions phase-space density in the explanatory text, the current script computes the critical temperature directly. That is the quantity most users need first. Once you know the estimated T_c, you can compare it with your measured cloud temperature. If your sample temperature is well above the result, the gas is still thermal. If it is near the result, you are approaching degeneracy. If it is below the result, a condensate fraction is expected in the ideal-gas picture.

How to Use

Using the calculator is straightforward, but the units matter. Enter the atom number N as a plain count of particles, not in scientific notation with hidden powers unless your browser input supports it directly. For example, one million atoms can be entered as 1000000 or 1e6. Enter the atomic mass in atomic mass units, often abbreviated amu or u. For rubidium-87, a common approximate value is 87. For sodium-23, use 23. The script converts that mass to kilograms internally.

For the three trap frequencies, enter values in hertz, not angular frequency in radians per second. This distinction is important. Experimental papers often quote trap frequencies as ordinary frequencies such as 50 Hz, 100 Hz, or 200 Hz, while theoretical formulas are often written using angular frequencies ω = 2πf. The calculator handles that conversion for you. If your trap is isotropic, simply enter the same value for all three directions. If it is anisotropic, enter each axis separately.

After you submit the form, the result area reports the estimated critical temperature. The first number is in nanokelvin, which is convenient for very cold trapped gases. The second number is the same temperature in microkelvin. These are just two ways of expressing the same result. A value of 280 nK, for example, is equal to 0.280 µK.

To get meaningful output, use physically sensible positive values. Negative frequencies or zero atom number do not correspond to a trapped Bose gas in this model. The script preserves the original calculator behavior, so it does not perform advanced validation beyond reading the numeric inputs. As a user, you should therefore treat the form as a scientific tool and supply realistic experimental parameters.

Formula

The ideal-gas critical temperature for Bose–Einstein condensation in a three-dimensional harmonic trap is obtained by summing or integrating the Bose occupation of excited states and asking when those states can no longer accommodate all particles. In the semiclassical limit, valid for sufficiently large atom numbers and level densities, the result is the familiar expression below.

Here k_B is Boltzmann's constant, ħ is the reduced Planck constant, N is the total number of bosons, and ζ(3) ≈ 1.20206 is the Riemann zeta function evaluated at 3. The quantity ω̄ is the geometric mean of the three angular trap frequencies:

Formula: ω̄ = ω_xω_yω_z^1/3

$\mathrm{\omega ̄}={{\omega }_x{\omega }_y{\omega }_z}^{\frac{1}{3}}$

Because users enter frequencies in hertz, the script first converts them to angular frequencies using ω = 2πf. It then computes the geometric mean and substitutes the result into the critical-temperature formula. The scaling is worth noticing. The critical temperature grows as N^1/3, so increasing the atom number by a factor of eight doubles T_c. It also grows linearly with the geometric mean trap frequency, so tighter confinement raises the transition temperature.

The atomic mass is included in the form because it is an important physical parameter in many BEC discussions and in related quantities such as the thermal de Broglie wavelength. In the current JavaScript implementation, however, the mass does not enter the final critical-temperature expression for the ideal harmonic-trap estimate. That is not a bug in the physics of the formula shown here; rather, it reflects the fact that for a harmonic trap written in terms of trap frequencies, the standard ideal-gas expression for T_c depends on N and the frequencies, not explicitly on mass. The mass field is still useful context for users and may support future extensions.

Worked Example

Suppose you are estimating the transition temperature for rubidium-87 in a nearly symmetric trap. Enter an atom number of 1000000, an atomic mass of 87 amu, and trap frequencies of 100 Hz, 100 Hz, and 100 Hz. The calculator converts each frequency to angular frequency by multiplying by 2π, so each axis corresponds to about 628.3 rad/s. Because the trap is symmetric, the geometric mean is the same value.

Substituting these numbers into the ideal-gas formula gives a critical temperature of roughly 450 nK, or about 0.450 µK. That means a thermal cloud at 1 µK would still be above the transition, while a cloud cooled to around 400 nK would be expected to begin condensing in the ideal trapped-gas picture. The exact onset observed in a real experiment can shift somewhat because of interactions, finite atom number, calibration uncertainty, and the details of how temperature is extracted from imaging data.

This example also shows how to think about scaling. If you keep the same trap but reduce the atom number from 10⁶ to 10⁵, the critical temperature falls by a factor of 10^1/3, which is about 2.15. If instead you keep the atom number fixed and double all three trap frequencies, the critical temperature doubles. Those trends are often more useful than the exact number because they tell you which experimental knobs most strongly affect the transition point.

Interpreting the Result

The output should be read as an estimate of the temperature at which the condensate begins to appear, not as a guarantee that every atom suddenly enters the ground state. Below T_c, only a fraction of the gas is condensed at first, and that fraction grows as the temperature is lowered further. In experiments, the onset is usually identified by the appearance of a narrow dense component in time-of-flight images or by fitting a bimodal density profile. The calculator therefore gives you the right temperature scale to watch for, but the actual condensate fraction depends on how far below the transition you are.

If your measured temperature is close to the calculated value, it is wise to treat the result as approximate rather than exact. Small uncertainties in trap frequencies can noticeably shift the estimate because the geometric mean enters linearly. Likewise, uncertainty in atom number matters through the cube-root dependence. A 30% error in atom number changes the predicted critical temperature by much less than 30%, but it still changes it enough to matter in precision work.

Limitations and Assumptions

This calculator uses the textbook ideal-gas expression for a harmonically trapped Bose gas. That means several real-world effects are intentionally neglected. The most important omission is interatomic interaction. Weak repulsive interactions, usually characterized by the s-wave scattering length, can shift the observed transition temperature and modify the density profile near the trap center. For many dilute alkali gases the shift is modest, but it is not always negligible if you need high accuracy.

Finite-size effects are also ignored. The derivation assumes a large enough atom number that the semiclassical density of states is a good approximation. For very small clouds, the discrete nature of the trap spectrum becomes more important and the simple formula becomes less precise. The model also assumes a stable three-dimensional harmonic trap. Strong anharmonicity, quasi-one-dimensional or quasi-two-dimensional confinement, lattice potentials, or rapidly time-varying traps require different treatments.

Another practical limitation is that the script reports only the critical temperature. It does not currently compute condensate fraction below T_c, interaction corrections, finite-size corrections, or phase-space density from an independently supplied temperature. The explanatory text mentions those ideas because they are part of the broader physical context, but the interactive output remains intentionally simple. That simplicity is useful for quick estimates, yet users should not mistake it for a complete thermodynamic model.

Finally, the calculator is only appropriate for bosons. Fermionic atoms do not undergo Bose–Einstein condensation as single particles because they obey the Pauli exclusion principle. Degenerate Fermi gases are instead described by the Fermi temperature and, in some cases, by paired superfluid states. If you are working with a Bose–Fermi mixture, this tool applies only to the bosonic component.

Further Context

The importance of BEC extends beyond a single formula. The first dilute-gas condensates created in the 1990s opened a route to studying quantum many-body physics with extraordinary control. Since then, trapped condensates have been used to investigate superfluid flow, quantized vortices, atom interferometry, optical lattices, synthetic gauge fields, and nonequilibrium dynamics. Even when a full experiment involves many corrections and calibrations, the ideal critical-temperature estimate remains one of the first numbers researchers compute because it sets the overall scale of the problem.

That is why a compact calculator like this is useful. It turns a theoretical expression into a quick planning tool. If you are comparing species, changing trap geometry, or estimating whether an evaporative cooling sequence is likely to reach degeneracy, the result gives immediate intuition. It is best used as a starting point: simple enough to be fast, but grounded in the standard physics of trapped bosonic gases.