Time-of-flight and correlation imaging
Turn an absorption image into a calibrated observable, separating ballistic momentum mapping, thermal widths, pair correlations, and density-dependent loss.
On this page
The question for today
What you’ll understand
- State when expansion images represent an initial momentum distribution.
- Extract a thermal temperature from variances with explicit expansion assumptions.
- Distinguish calibrated pair correlations and three-body loss from raw atom-count fluctuations.
A cold-atom image is the endpoint of a preparation, release, expansion, illumination, and detection sequence. Inferring temperature or correlations requires a model for that entire sequence. Free expansion can magnify a momentum distribution, but interactions can transform it during flight. Absorption measures a column density with finite resolution rather than every atom's three-dimensional position. A persuasive claim therefore connects a fitted parameter to a calibrated observable and checks which alternative physical effects could produce the same image.
Ballistic expansion maps velocities only under conditions
After forces and collisions become negligible, each particle moves from its initial position with its initial velocity. At long flight times the displacement from velocity can exceed the initial cloud extent. Only in that far-field limit does the expanded density directly reproduce the initial momentum distribution, with a Jacobian needed to conserve particle number. A Fourier distribution defined in wavenumber has a different Jacobian from one defined in momentum; the normalization must be stated.
For a noninteracting classical thermal sample with no initial position–velocity covariance, the spatial variance is its initial variance plus flight time squared times velocity variance. This permits temperature extraction without assuming that the initial cloud was a point. Collective flow introduces covariance or an additional velocity component. A dense condensate can convert interaction energy into anisotropic expansion, so its final width is not automatically a thermal velocity width. Changing scattering length before release can help but also changes the preparation protocol.
A camera measures a calibrated column density
In weak-probe absorption imaging, the logarithm of reference intensity divided by transmitted intensity is optical depth. It equals an absorption cross section times density integrated along the imaging axis under ideal conditions. Probe detuning, polarization, transition strength, saturation, and optical pumping affect that cross section. Stray light and camera offsets can corrupt the logarithm, particularly in optically thick regions. Independent atom-number calibration is necessary before using a density-sensitive interaction or loss law.
The imaging point-spread function broadens narrow structures and averages short-distance fluctuations. A Gaussian thermal fit is appropriate only in its stated classical regime; quantum-degenerate Bose and Fermi distributions have different shapes. Fitting the thermal wings can reduce condensate contamination, but limited signal and interactions still need assessment. A bimodal distribution or interference peaks provide evidence for coherence or condensation within a model. Superfluidity additionally concerns response and phase rigidity, so one visually narrow peak does not establish every many-body property.
Correlations and loss test different parts of the state
Repeated images allow density–density correlations after removing the smooth mean profile and accounting for self-counting and detector noise. Bosonic bunching and fermionic antibunching reflect exchange statistics, while additional correlations can reveal lattice order or pairs. Column integration and finite resolution reduce observable contrast. Shot-to-shot number variation and common technical noise can generate correlations unrelated to the desired quantum effect. Report normalization, spatial binning, noise subtraction, and an appropriate uncorrelated reference alongside any correlation peak.
Lifetime measurements supply another diagnostic because three-body recombination weights the cube of local density. The atom-loss coefficient must include a declared counting convention, including how many atoms leave per event and any local correlation factors. As density changes, loss is generally nonexponential, even at fixed volume. Trapped profiles, heating, and changing correlations can further modify the curve. A fitted loss coefficient and a fitted temperature are complementary evidence, and neither replaces calibration of the evolving density.
Work it through
Infer temperature from a width difference
A noninteracting classical thermal cloud has mass m = 1.38 × 10⁻²⁵ kg, initial rms width 10 μm, and rms width 30 μm after 10 ms. Take kB = 1.38 × 10⁻²³ J/K and zero initial position–velocity covariance. Find its directional temperature.
The variance increase is . Subtract variances rather than widths.
Dividing by flight time squared gives .
The temperature is . Flow, interactions, or imaging broadening would require a modified inference.
The directional thermal temperature is 80 nK under the stated ballistic, covariance, and calibration assumptions.
The living notebook
Explore the model
Interactive lab · probability
Count variance needs its own normalization
Reset and set BOTH Success probability p to 0.5 and Trials n to 10; predict mean 5 and variance 2.5. Then Reset and set BOTH p to 0.5 and n to 20; predict mean 10 and variance 5. Relative variance therefore falls from 0.10 to 0.05 even though absolute variance grows. This exact binomial display illustrates counting moments only. It does not generate camera frames, quantum bunching, pair correlations, random samples, or a thermometer.
- Mean · np
- 5
- Variance · np(1−p)
- 2.5
- Total probability
- 1
These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forA fluctuation amplitude has meaning only after its mean, normalization, and detection model are specified.
Pause & reason
Check your understanding
Your turn
Try it yourself
An ideal uniform fixed-volume sample loses atoms only through a constant three-body coefficient L₃ = 2.0 × 10⁻⁴¹ m⁶/s, with density initially 1.0 × 10²⁰ m⁻³. Neglect heating and changes of correlations. Find the time until density halves. Explain why an exponential with the initial per-atom rate is inappropriate.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
At fixed volume, dn/dt = −L₃n³. Integrate n⁻³ dn before applying the halved-density condition.
Compare with the solution
Integration gives . Initially . For n/n₀ = 1/2 the square root equals two, so . The instantaneous per-atom rate decreases as n², so keeping its initial value in an exponential overestimates continuing loss. This fixed-volume illustration omits trapped-profile evolution and recombination heating.
Keep the idea
Take it with you
- Time-of-flight becomes momentum imaging only after ballistic and far-field conditions are justified.
- Temperature, atom number, and correlation contrast depend on optical calibration and resolution.
- Density-dependent loss requires an evolving-density model and an explicit atom-counting convention.
Read further
Original sources and deeper explanations to follow at your own pace.
- Bloch, Dalibard and Zwerger, II.3 and III.1–III.2: ballistic mapping, self-correlations, and resolution limits
- Dalfovo et al., IV.4: interaction-driven condensate expansion
- Grimm, Weidemüller and Ovchinnikov, III.2.2–III.3: diagnostics and density-dependent loss
- Reinaudi, Lahaye, Wang and Guéry-Odelin, Eqs. 1–3 and low-intensity limit: absorption calibration
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.