← SciSim / Physics

Bose-Einstein Condensate & Superconductivity

Quantum Mechanics #64
Section 01
Interactive Simulation
Bose-Einstein Condensate & Superconductivity — SciSim
Ready
T / Tc
0.00
Condensate fraction
0.00
λ_th
0.00μm
SC State
Resistance
0.00Ω
Cooper pairs
0
Controls
Parameters
Temperature T500 nK
# Atoms200
Magnetic Field B0.0 T
SC Tc90 K
Speed1.0 ×
Display
Section 02
The Idea, Step by Step

Cool any gas and the atoms slow down — that part is familiar. But cool a gas of the right kind of atoms extraordinarily cold, and something with no everyday parallel happens: a huge fraction of the atoms suddenly stop being separate specks and start acting like one single thing, moving in perfect lockstep. Picture a stadium crowd: at room temperature everyone is milling about randomly; in a Bose-Einstein condensate it is as if every person in the stadium suddenly does the exact same motion at the exact same instant. That shared, single quantum state is the condensate.

Why cold makes atoms merge

Here is the key picture for a high-schooler: every atom is also a tiny wave, and the colder it gets, the wider that wave spreads. The width is the thermal de Broglie wavelength $\lambda_{\text{th}}$, and it grows as the temperature $T$ falls:

Atom-Wave Width
$$\lambda_{\text{th}} = \frac{h}{\sqrt{2\pi m k_B T}}$$

For a rubidium atom at about $170$ nK — that is $0.00000017$ degrees above absolute zero — this width works out to roughly $0.5\ \mu\text{m}$. The atoms in the trap sit only about $0.3\ \mu\text{m}$ apart, so each atom's wave is now wider than the gap to its neighbour. The waves overlap, the atoms can no longer be told apart, and they collapse into one shared state. Condensation happens precisely when the waves get this crowded, the condition $n\,\lambda_{\text{th}}^3 \approx 2.6$.

The precise version, and the leap to superconductors

Solving that crowding condition for temperature gives the critical temperature $T_c \propto n^{2/3}/m$ (full form in the next section), and below it the share of atoms in the condensate climbs as $N_0/N = 1-(T/T_c)^{3/2}$ — zero just above $T_c$, rising to everything at absolute zero. The deep idea then carries straight over to superconductivity: electrons are loners (fermions, no sharing allowed), but in a cold metal two of them link up through the vibrating lattice into a Cooper pair. A pair behaves like a boson, so the pairs can all condense into one state — and a current made of that single locked-together state flows with exactly zero resistance. The sliders let you drive both stories: Temperature and SC Tc set how far below the transition you are, and Magnetic field B tests how much field a superconductor will throw out before it gives up.

Try this in the sim above

In BEC Cooling mode, drag Temperature down past $T_c=170$ nK and watch the orange thermal cloud collapse into the glowing cyan condensate at the trap centre while the condensate fraction $N_0/N$ climbs toward 100%. Switch to Momentum Distribution and cool again: a sharp central spike grows out of the broad thermal hump — the famous bimodal "fingerprint" that won the 2001 Nobel Prize. Finally, in Meissner Effect mode set the temperature below the SC $T_c$, then raise $B$ and watch the field lines bend around the sample instead of through it — until $B$ gets large enough to destroy the superconducting state.

Section 03
Equations & Derivation

Bose-Einstein Condensation (BEC) and superconductivity (SC) are macroscopic quantum phenomena that emerge when integer-spin particles ('bosons') pile into a single quantum state below a critical temperature. In atoms, the bosons are dilute neutral atoms (like ⁸⁷Rb). In a superconductor, two electrons (fermions) pair up via phonon exchange to form a composite boson — the Cooper pair — which then condenses.

1. The de Broglie Thermal Wavelength

At temperature $T$, an atom of mass $m$ has thermal momentum $p \sim \sqrt{m k_B T}$, giving a wavelength:

Thermal de Broglie Wavelength
$$\lambda_{\text{th}} = \frac{h}{\sqrt{2\pi m k_B T}}$$

BEC sets in when $\lambda_{\text{th}}$ becomes comparable to the inter-particle spacing $n^{-1/3}$ — i.e., when the wavefunctions overlap and the particles can no longer be distinguished.

2. BEC Critical Temperature

For an ideal gas of $N$ bosons in volume $V$, condensation begins when:

$$n\,\lambda_{\text{th}}^3 \approx 2.612$$

Solving for $T$ gives the critical temperature:

BEC Critical Temperature
$$\boxed{\;T_c = \frac{2\pi\hbar^2}{m k_B}\left(\frac{n}{2.612}\right)^{2/3}\;}$$

For ⁸⁷Rb at $n \approx 2.8\times 10^{13}$ cm⁻³, $T_c \approx 170$ nK — a hundred-millionth of a degree above absolute zero. Below $T_c$, the condensate fraction is:

$$\frac{N_0}{N} = 1 - \left(\frac{T}{T_c}\right)^{3/2}$$

3. The Macroscopic Wavefunction

All condensed atoms share a single wavefunction $\psi(\vec r,t)$ — typically the ground state of the trap. This wavefunction has a definite phase and obeys the Gross-Pitaevskii equation (for a weakly-interacting gas):

Gross-Pitaevskii Equation
$$i\hbar\frac{\partial\psi}{\partial t} = \left(-\frac{\hbar^2}{2m}\nabla^2 + V(\vec r) + g|\psi|^2\right)\psi$$

where $g = 4\pi\hbar^2 a_s/m$ is the contact interaction strength ($a_s$ is the s-wave scattering length).

4. BCS Theory of Superconductivity

Below $T_c$, electrons near the Fermi surface attract each other via virtual phonon exchange. Bardeen, Cooper, and Schrieffer (1957) showed that any net attractive interaction, however weak, makes the Fermi sea unstable to pair formation:

BCS Energy Gap (T = 0)
$$\Delta(0) = 2\hbar\omega_D \exp\!\left(-\frac{1}{N(0)V_0}\right)$$

where $\omega_D$ is the Debye frequency and $N(0)V_0$ is the dimensionless coupling strength. The full temperature-dependent gap satisfies:

$$\Delta(T) \approx \Delta(0)\sqrt{1 - T/T_c} \quad (T \to T_c)$$

The BCS prediction $2\Delta(0)/(k_B T_c) = 3.528$ is a universal ratio for weak-coupling superconductors and agrees well with most conventional SCs.

5. Meissner Effect & London Penetration

A superconductor expels magnetic field from its interior — the Meissner effect. The field decays exponentially over the London penetration depth:

$$B(x) = B_0\,e^{-x/\lambda_L}, \qquad \lambda_L = \sqrt{\frac{m_e}{\mu_0 n_s e^2}}$$

For Type II SCs, above a lower critical field $B_{c1}$, the field penetrates as quantised vortices each carrying flux $\Phi_0 = h/(2e) = 2.07\times 10^{-15}$ Wb (the 2e is the Cooper pair charge — direct evidence of pairing).

Symbol Table

SymbolMeaningSI Unit / Value
$T_c$Critical temperatureK
$\lambda_{\text{th}}$Thermal de Broglie wavelengthm
$N_0/N$Condensate fraction
$\Delta$Superconducting energy gapJ (often meV)
$\Phi_0 = h/(2e)$Magnetic flux quantum$2.07\times 10^{-15}$ Wb
$\lambda_L$London penetration depthm (~50-500 nm)
$\xi$BCS coherence lengthm (~1-1000 nm)
$g = 4\pi\hbar^2 a_s/m$BEC interaction strengthJ·m³

Mapping to the simulation

The BEC mode shows $N$ atoms in a harmonic trap; their kinetic energy is sampled from a Maxwell-Boltzmann distribution at temperature $T$. As you lower $T$ below $T_c$ (drawn from the formula above), atoms condense into the trap ground state — they collapse to the centre and move slowly. Cooper-pair mode visualises the phonon-mediated attraction; Meissner mode shows field expulsion below $T_c$ with the London decay; Phase diagram tracks $B$-$T$ for Type II SCs with the upper/lower critical-field curves.

Reference: Pethick & Smith — Bose-Einstein Condensation in Dilute Gases, 2nd Ed., Cambridge, Ch. 2 & 6; Tinkham — Introduction to Superconductivity, 2nd Ed., Dover, Ch. 1-3; Annett — Superconductivity, Superfluids and Condensates, Oxford, 2004, Ch. 3-5; Bardeen, Cooper & Schrieffer — Phys. Rev. 108, 1175 (1957) — original BCS paper.
Section 04
Frequently Asked Questions
BEC mode shows atoms in a harmonic optical/magnetic trap. Above Tc, atoms move chaotically (orange = thermal). As you lower T below Tc, the condensate fraction (cyan, denser at trap centre) grows according to N₀/N = 1−(T/Tc)^(3/2). Momentum mode shows the famous bimodal distribution: thermal Gaussian wings around a sharp central peak — exactly the signature observed in the 1995 Cornell-Wieman-Ketterle experiments. Cooper-pair mode shows two electrons (red) loosely bound through phonon exchange (yellow waves) at separations of ~ξ (coherence length, often >>atomic spacing). Meissner mode shows magnetic field lines being expelled below Tc.
BEC: Used in atomic clocks (~10⁻¹⁹ precision), atom interferometry for gravimetry, sensitive force/rotation sensing, and as 'quantum simulators' for studying lattice models that classical computers can't solve. Superconductors: MRI machines (~1.5-3 T magnets are wound from NbTi superconductor), the LHC's bending magnets at CERN, maglev trains (e.g., Shanghai's Transrapid), Josephson voltage standards, SQUID magnetometers (sensitive to 10⁻¹⁵ T), and emerging quantum computers (transmon qubits are superconducting circuits).
Bosonic statistics — multiple particles in the same state are favoured by enhancement factor (1+nₖ). At low T, the lowest-energy state becomes massively occupied because adding a particle to it is statistically weighted against a thinly-occupied excited state. There's no force pushing atoms together; it's a pure statistical effect. This is unique to integer-spin particles. Half-integer-spin particles obey the Pauli exclusion principle — only one per state — so they form a Fermi sea, not a condensate.
They do, via direct Coulomb repulsion. But in a metal, an electron locally distorts the positively-charged ionic lattice, leaving an over-screened region of net positive charge. A second electron, moving past after a delay (set by the slow phonon timescale), is attracted to that lingering excess. The retarded interaction, mediated by phonons, can overcome the bare Coulomb repulsion in the right energy window (within ℏω_D of the Fermi surface). The net effect is the BCS attraction.
It's all about density and mass. BEC Tc ∝ n^(2/3)/m. Atomic gases have very low density (~10¹⁴ cm⁻³ — a few hundred thousand times less dense than air) and atomic masses are huge (~10⁵ × electron mass). Both factors push Tc down. Conduction electrons in a metal have density ~10²² cm⁻³ and mass mₑ — a billion times denser and much lighter than rubidium atoms, giving much higher Tc. The two phenomena look similar but operate at vastly different scales.
Discovered 1986 (Bednorz & Müller, YBa₂Cu₃O₇ at Tc ≈ 92 K — above liquid-nitrogen 77 K). The mechanism is still debated; it's clearly NOT phonon-mediated BCS. Most theories invoke spin fluctuations in 2D copper-oxide planes. The pairing is d-wave (Δ has nodes) rather than BCS s-wave, and the normal state is a strange metal that violates ordinary Fermi-liquid theory. Room-temperature superconductivity remains a holy grail; recent claims (2020-2024) under high pressure are still controversial.
Closely related but not identical. BEC is a phenomenon of non-interacting (or weakly interacting) bosons macroscopically occupying one state. Superfluidity (zero viscosity) requires an additional ingredient: an excitation spectrum with a linear (sound-like) dispersion at low momentum, so the Landau criterion forbids viscous dissipation below a critical velocity. Pure non-interacting BEC is technically NOT superfluid. Liquid ⁴He below 2.17 K is superfluid AND has a BEC-like macroscopic wavefunction; dilute ⁸⁷Rb gas is a textbook BEC and shows superfluid behaviour because of its weak repulsive interaction.
Resources: MIT OpenCourseWare 8.422 'Atomic and Optical Physics II'; HyperPhysics — 'Superconductivity'; Khan Academy — 'States of Matter'; Cornell, E. A. & Wieman, C. E. — 'Nobel Lecture: Bose-Einstein Condensation in a Dilute Gas', Rev. Mod. Phys. 74, 875 (2002); Ketterle, W. — 'Nobel Lecture: When atoms behave as waves', Rev. Mod. Phys. 74, 1131 (2002).
Section 05
Common Misconceptions
❌ Misconception: 'A BEC is just a really cold gas where everything stops moving.'
✅ Correction: Atoms in the condensate are NOT motionless. They occupy the ground state of the trap, which has zero-point kinetic energy ~ℏω/2 — a real, observable energy. Their velocity distribution is the Fourier transform of the trap ground-state wavefunction (a narrow Gaussian for a harmonic trap), much narrower than thermal but not zero. Time-of-flight imaging of BECs reveals this characteristic narrow distribution.
📖 Reference: Pethick & Smith — Bose-Einstein Condensation in Dilute Gases, 2nd Ed., §6.1
❌ Misconception: 'Cooper pairs are tightly-bound little molecules of two electrons.'
✅ Correction: Cooper pairs are NOT compact molecules. The coherence length ξ — the spatial scale over which the pair correlation extends — is typically 100-1000 nm in conventional SCs, far larger than the Cu-Cu spacing (~0.3 nm) or even typical sample dimensions of nanometre-scale wires. Many millions of pairs overlap in any region. The 'pair' is a many-body correlation, not a localised diatomic structure.
📖 Reference: Tinkham — Introduction to Superconductivity, 2nd Ed., §1.6 'Coherence Length'
❌ Misconception: 'Magnetic field can't get into a superconductor at all.'
✅ Correction: Wrong for two reasons. (1) The field penetrates a thin surface layer of thickness λ_L (typically 50-500 nm) — that's the London depth. Inside, B exponentially decays but isn't zero at the surface. (2) Type II superconductors (most useful ones) admit quantised flux vortices above a lower critical field B_c1 — the field gets in via filaments each carrying Φ₀ = h/(2e). Only the bulk between vortices remains field-free.
📖 Reference: Annett — Superconductivity, Superfluids and Condensates, §3.6 'Type-I and Type-II Superconductors'
❌ Misconception: 'BEC is a fourth state of matter, distinct from solid/liquid/gas.'
✅ Correction: BEC is a thermodynamic phase but classifying it 'beyond solid/liquid/gas' is a popular-science oversimplification. It's a quantum-degenerate gas — still a gas thermodynamically, just one where quantum statistics dominate. The 'fourth state' label confuses macroscopic quantum coherence with the broader question of what counts as a phase of matter. Plasmas, superfluids, and various exotic phases are equally valid candidates if we're labelling phases beyond the classical three.
📖 Reference: Anderson, M. H. et al. — Science 269, 198 (1995) — original BEC observation paper
❌ Misconception: 'High-Tc superconductors will give us room-temperature SC any day now.'
✅ Correction: Discovery of cuprates in 1986 sparked exactly this hope. Nearly 40 years later, the highest reproducibly-confirmed Tc at ambient pressure remains ~133 K (HgBa₂Ca₂Cu₃O₈₊δ). Recent claims of room-temperature SC (e.g., LK-99 in 2023, hydrides under high pressure) have either failed independent verification or required pressures (~1 million atmospheres) that make them practically useless. The mechanism of high-Tc SC is itself unsolved, so engineering toward higher Tc remains largely empirical.
📖 Reference: Physics Today, 'The state of high-Tc theory' — multiple reviews; Norman, M. R. — Rep. Prog. Phys. 79, 074502 (2016)
❌ Misconception: 'Ohm's law breaks down in a superconductor — that's why R = 0.'
✅ Correction: Ohm's law isn't 'broken'; it's that the conductivity becomes formally infinite, so V = IR is degenerate (V→0 for any finite I). The proper description uses the Meissner-London equation: J = −(n_s e²/m) A (linear in vector potential, not in field). Persistent currents survive years without decay — direct evidence of zero resistance. The textbook 'V = IR' formalism doesn't fail; it just becomes trivial.
📖 Reference: Tinkham — Introduction to Superconductivity, 2nd Ed., §3 'The Magnetic Properties of Superconductors'
Misconception research: Bouchée, A. & Marrongelle, K. — 'Difficulties students have with the concept of BEC', J. Mod. Phys. review articles; Kaur & Stadermann — Eur. J. Phys. 39, 035702 (2018) on superconductivity teaching; Norman, M. R. — Rep. Prog. Phys. 79, 074502 (2016).