VISUALIZATIONV ISUALIZATION C ORNER Editor: Jim X. Chen, [email protected] VISUALIZING BOSE–EINSTEIN CONDENSATES By Peter M. Ketcham and David L. Feder OSE–EINSTEIN CONDENSATES ARE THE COLDEST SUB- ity to do so, and the result is the for- mation of tiny quantum vortices, each B STANCES EVER MADE. ALONG WITH THIS DISTINCTION, with a quantized amount of superflow around it. BECS HAVE OTHER INTRIGUING PROPERTIES: FOR EXAMPLE, THEY These quantum whirlpools are pro- duced in BECs by rotating the mag- BEHAVE MUCH LIKE LIGHT DOES IN AN OPTICAL FIBER. IN FACT, netic trap in which the BEC is con- tained—much like trying to stir your many research groups have produced Quantum Vortices coffee by squeezing the cup and spin- “atom lasers” formed from traveling Laboratory-produced BECs aren’t ning yourself around. BECs range BECs. They also behave similarly to quite as simple as the ones Bose and from a millimeter to a centimeter in superfluids, much like liquid helium or Einstein originally predicted because size, so you’d think it should be easy to superconductors. By studying BECs, they are made from real alkali atoms, detect the vortices simply by shining we’ve gained important insights into not ideal gases. For instance, unlike an light on the cloud, taking a photo- much more complicated, but closely ideal gas, atoms can’t pass through one graph, and noting the tiny holes in the related, systems. another. They interact strongly, and at- picture where the atom density has BECs created from dilute alkali va- tract each other if they get too close to- gone to zero in the vortex cores. pors were first produced in the labora- gether. Thus, keeping the gas dilute Unfortunately, in spite of many ef- tory in 19951–4 using laser cooling and enough is very important; otherwise, forts, it wasn’t until several years after trapping techniques developed over you would end up with a lump of metal BECs were first formed that re- the past 20 years. Creating a BEC re- instead of a BEC. Conversely, the searchers saw quantum vortices in quires temperatures between one bil- atoms also repel each other if kept suf- gases made of a single species of atom.8 lionth and one millionth of a degree ficiently far apart, and it’s this repulsive (The first quantum vortices in BECs above absolute zero. interaction that gives rise to the most were produced in two-species gases, so Once the cooling technology became interesting properties of BECs. For that one species filled the core of the available, an explosion of research in more on the physical intricacies of vortex formed in the other species. the area followed; over 30 experimen- BECs, see the “Quantum Degeneracy This made the vortex large enough to tal groups around the world have made and BECs” sidebar. image, but the cloud as a whole had no BECs. However, BECs are not a new One of the most interesting phe- density depression in the vortex core.) idea. Satyendra Nath Bose and Albert nomena exhibited by superfluids is the The inability to see single-species gas Einstein predicted the Bose–Einstein formation of quantum vortices. When quantum vortices wasn’t properly un- condensation phenomenon in dilute an ordinary fluid is rotated—for exam- derstood until full three-dimensional ideal gases over 75 years ago.5–7 ple, when you stir milk into a cup of numerical simulations and visualiza- In this article, we present a method coffee—a large vortex, or depression, tions of rotating BECs began to for visualizing numerical simulations is formed in the center. But a super- emerge in the 1990s. These simula- of BECs, discuss the visualization fluid, by its very nature, wants to flow tions demonstrated that the rotation process, which includes selecting a without friction. This means that if frequencies needed to produce quan- color scheme and the application of a you tried to stir it with a spoon, it tum vortices were higher than previ- volume rendering model, and describe would attempt to flow around the ously expected. our visualization’s impact on associated spoon with no resistance. Of course, The visualizations also revealed that scientific research. there are some limitations on its abil- the vortices themselves were generally 86 Copublished by the IEEE CS and the AIP 1521-9615/03/$19.00 © 2003 US Gov’t Work Not Protected by US Copyright COMPUTING IN SCIENCE & ENGINEERING Quantum Degeneracy and BECs A Bose–Einstein condensate forms when a gas of particles obeying Bose–Ein- stein statistics, bosons, becomes quantum degenerate. In his 1924 doctoral the- sis, Louis de Broglie proposed that all particles possess a quantum wave nature quite curved. Due to the small size of in addition to a quantum particle nature. When viewed as a wave, a particle the vortices, this curvature meant that has an associated quantum wavelength. This quantum wavelength is inversely any photograph of the entire BEC proportional to the particle’s mass and temperature, so the colder the temper- wouldn’t detect the holes because the ature and the lighter the mass, the longer the quantum wavelength. Particles position of the holes changes down the with long quantum wavelengths more readily manifest their wavy character. line of sight, thus the cloud itself would Quantum degeneracy occurs when the quantum wavelength becomes obscure the vortices. comparable to the average spacing between particles. In practice, this means Vortices were finally imaged experi- either very high particle densities or very cold temperatures. Metals are good mentally using two techniques. At conductors because the electrons carrying the charge are quantum degener- École Normale Supérieure, Jean Dal- ate, which is possible at room temperatures due to the high electron densities. ibard’s group found that if they re- Two kinds of particles exist in the known universe: bosons and fermions, the leased the atoms from their trap and al- latter of which obey Fermi–Dirac statistics. (In certain artificial two-dimensional lowed the BEC to expand, the vortex systems, such as semiconductor heterojunctions, particles can actually have frac- cores were magnified sufficiently to be tional statistics, intermediate between bosons and fermions. These instances give detectable.9,10 This approach works rise to the fractional quantum Hall effect at low temperatures and high magnetic well when the number of atoms isn’t fields.) Bosons, which have integer values of something called spin, include pho- too large, on the order of a few hun- tons (the quantum particles of light). Fermions, which have half-integer spin, in- dred thousand. Wolfgang Ketterle’s clude the electrons, protons, and neutrons making up ordinary atoms. BEC, however, contained several mil- Atoms themselves are either bosons or fermions, depending on the total lion atoms; his solution was to pump a number of their fermion constituents: an even number of fermions gives a thin section of the cloud into another bosonic atom. For example, hydrogen is a boson, with one proton and one quantum state using a laser and then electron. In fact, most efforts to produce a BEC over the past two decades fo- image only that slice.11 cused on hydrogen, and in 1998, a hydrogen BEC was finally produced in the laboratory by Dan Kleppner and Thomas Greytak’s group at Harvard Univer- Numerical Simulation sity.1 Hydrogen was long considered a good candidate because of its light We can numerically model the forma- mass. tion of quantum vortices in BEC su- In principle, the gas needed to be only a few thousandths of a degree above perfluids by solving differential equa- absolute zero to become a BEC. But achieving these temperatures in hydrogen tions very similar to those used to proved to be very difficult, and efforts turned instead to alkali atoms where laser model ocean or air currents. The cooling could be used. Alkali atoms are much heavier than hydrogen, however, main differences are the total lack of so the temperatures needed to reach BEC are perhaps 100 times lower. friction and the inclusion of terms that take the quantum nature of the atoms Reference into account. 1. D.G. Fried et al., “Bose–Einstein Condensation of Atomic Hydrogen,” Physical Rev. Letters, All the atoms in a BEC are in a sin- vol. 81, no. 18, 1998, pp. 3811–3814. gle quantum state (which is why it’s called quantum degeneracy). This state is represented by a wave func- tion ψ, which is a complex number characterized by two values: a density flow: the superfluid velocity ν is pro- The data values represent the wave φ ρ and a phase φ, so that ψρ= ei . portional to the gradient of the phase function ψ on a three-dimensional grid Here, i =−1 is the i m a g i n a r y ν ∝ ∇φ. with 200 grid points in each dimension. unit, so that exp (i × 0) = 1 and exp (i Visualization provides a way to examine ×π) = –1. The density gives the Visualization this data set at a high conceptual level. cloud’s overall profile, expressing the The product of the numerical simula- average number of atoms in a given tion of BEC vortex formation we used Process time. The phase is a quantum mechan- is a data set containing 8 million com- The visualization process includes ical variable that defines the superfluid plex numbers expressed in polar form. three main activities: JANUARY/FEBRUARY 2003 87 V ISUALIZATION C ORNER of complex num- expressed linearly in the sense that it is bers expressed in not expressed in an angular or radial polar form. The fashion. This allows a mapping from HSV (hue, satura- the BEC wave function’s density to the tion, value) color value component of the HSV color model offers such model.