Feature articles
Argonne Update
Stars surrender secrets to supercomputers

Electricity controls nanocrystal architecture

Stars surrender secrets
to supercomputers

A computer simulation of materials in a white dwarf star.

STAR FLASH – This simulation from the Center for Astrophysical Thermo-nuclear Flashes shows the density changes in the materials of a white dwarf star.

University of Chicago and Argonne scientists are preparing to run the most advanced supercomputer simulation of an exploding star ever attempted. Tomasz Plewa, senior research associate at the university's Center for Astrophysical Thermonuclear Flashes, expects the simulation to reveal the mechanics of exploding stars, called supernovae, in unprecedented detail.

The simulation is made possible by a special allocation of 2.7 million hours of supercomputing time to the Flash Center from the U.S. Department of Energy (DOE). The Flash Center typically uses less than 500,000 hours of supercomputer time annually.

"This is like a big grant. This is beyond imagination," said Plewa, who submitted the Flash Center proposal on behalf of a research team at the university and Argonne. The Argonne team includes Argonne's Mathematics and Computing Science Division (MCS) Director Rick Stevens, and MCS researchers Justin Binns, Anthony Chan, Randy Hudson, Rusty Lusk, Michael Papka, Katherine Riley and Andrew Siegel.

The Flash Center project was one of three selected to receive supercomputer time allocations under a new competitive program announced last July by Secretary of Energy Spencer Abraham.

The supercomputer time will help the Flash Center more accurately simulate the explosion of a white dwarf star, one that has burned most or all of its nuclear fuel. These supernovae shine so brightly that astronomers use them to measure distance in the universe. Nevertheless, many details about what happens during a supernova remain unknown.

Simulating a supernova is computationally intensive because it involves vast scales of time and space. White dwarf stars gravitationally accumulate material from a companion star for millions of years, but ignite in less than a second. Simulations must also account for physical processes that occur on a scale that ranges from a few hundredths of an inch to the entire surface of the star, which is comparable in size to Earth.

Similar computational problems vex the DOE's nuclear weapons Stockpile Stewardship and Management Program. In the wake of the 1996 Comprehensive Test Ban Treaty, the reliability of the nation's nuclear arsenal must now be tested via computer simulations rather than in the field.

"The question ultimately is how is the nuclear arsenal aging with time, and is your code predicting that aging process correctly," Plewa said.

Flash Center scientists verify the accuracy of their code by comparing the results of their simulations both to laboratory experiments and to telescopic observations. Spectral observations of supernovae, for example, provide a sort of bar code that reveals which chemical elements are produced in the explosions. Those observations currently conflict with simulations.

"You want to reconcile current simulations with observations regarding chemical composition and the production of elements," Plewa said.

Scientists also wish to see more clearly the sequence of events that occurs immediately before a star goes supernova. It appears that a supernova begins in the core of a white dwarf star and expands toward the surface like an inflating balloon.

According to one theory, the flame front initially expands at a relatively "slow" subsonic speed of 60 miles per second. Then, at some unknown point, the flame front detonates, and accelerates to supersonic speeds. In the ultra-dense material of a white dwarf, supersonic speeds exceed 3,100 miles per second.

Another possibility: The initial subsonic wave fizzles when it reaches the outer part of the star, leading to a collapse of the white dwarf, the mixing of unburned nuclear fuel and then detonation.

Flash Center scientists already are on the verge of recreating this moment in their simulations. The extra computer time should push them across the threshold.

The center will increase the resolution of its simulations to one kilometer (six-tenths of a mile) for a whole-star simulation. Previously, the center could achieve a resolution of five kilometers (3.1 miles) for a whole-star simulation, or 2.5 kilometers (1.5 miles) for a simulation encompassing only one-eighth of a star. – Steve Koppes

For more information, please contact Catherine Foster (630/252-5580 or media@anl.gov) at Argonne.

Go to next story: Electricity controls nanocrystal architecture