ASTRONOMER Edwin P. Hubble, for whom the Hubble Space Telescope is named, discovered in 1929 that the universe is expanding. Astronomical observations almost 70 years later revealed the startling news that not only is the universe expanding, but also the rate at which it is expanding is accelerating. Theories for the observed acceleration abound, and they differ in many ways. But all theories ascribe the acceleration to dark energy—matter of unknown composition that does not emit or reflect electromagnetic radiation and is therefore difficult to observe directly. Dark energy, which may constitute up to 70 percent of all matter and energy in the universe, is the greatest mystery in astronomy today.
News about the accelerating expansion, reported in 1998, was based on observations of exploding stars known as Type Ia supernovae (SNe Ia). (See the box below.) These thermonuclear explosions mark the death throes of white dwarf stars. SNe Ia are less common than Type II supernovae, which occur from gravitational collapse. Because most SNe Ia are observed to have quite similar characteristics, astronomers have used them as “standard candles,” which serve as “yardsticks” for calculating the distance to their host galaxies from Earth. Using the calculated distances and observed velocity (determined by the red shift of the light), cosmologists gauge the degree of acceleration of these galaxies and, by extension, the influence of dark energy.
The Department of Energy (DOE) Office of Science has identified “dark energy and the search for genesis” as one of its top priorities. Fusion energy, advanced computing, and nuclear matter at the extremes are among others. Solving the mystery of dark energy is a key component of Quantum Universe: The Revolution in 21st-Century Particle Physics, a recent report from the Office of Science and the National Science Foundation that identifies the most compelling questions facing contemporary particle physics research. DOE and the National Aeronautics and Space Administration (NASA) plan to launch the SuperNova Acceleration Probe (SNAP) before 2020 to more precisely measure the expansion of the universe and to investigate the nature of the dark energy accelerating this expansion. SNAP is one of several Earth- and space-based missions planned over the next 15 years to quantitatively observe supernovae, examine dark energy, and better understand the beginning and possible future of our universe.
The application of SNe Ia as standard candles for determining distances to galaxies has revolutionized cosmology, giving astronomers an essential tool to extend their knowledge of the universe’s expansion and acceleration. However, that knowledge depends on the assumption that all SNe Ia are truly “standard.” Does each explosion provide the same intrinsic brightness? With the planned space missions, the number of observed SNe Ia is expected to increase from a few thousand to hundreds of thousands in the near future. Accurately understanding the behavior of SNe Ia is essential to making the investments in science pay off.
Birth of a Candle
The study was a part of Livermore’s Computing Grand Challenge Program, which allocates millions of central-processing-unit (CPU) hours on Laboratory supercomputers to unclassified projects that support Livermore’s missions. The supernova simulations were performed by the Computational Astrophysical Consortium (CAC), with members from the University of California at Santa Cruz, Lawrence Berkeley National Laboratory, State University of New York at Stony Brook, and Lawrence Livermore. The team was allocated 4 million CPU hours on Atlas, one of Livermore’s workhorse supercomputers for unclassified research. Previously, the team had used supercomputers at Oak Ridge National Laboratory, the National Energy Research Scientific Computing Center at Lawrence Berkeley, and NASA Ames. Likening these supercomputer resources to a baseball lineup, Hoffman says, “Atlas is our ‘cleanup’ hitter, not because it is a particularly heavy hitter—they all are—but because it gets a hit every time.”
SNe Ia evolve from a complex of two stars, in which a dead stellar core called a white dwarf accretes material from a companion star. As the white dwarf collects more material, it becomes denser and hotter. Eventually, conditions become so extreme that carbon atoms fuse and ignite one or more miniscule nuclear flames, which erupt near the center of the dwarf. A runaway thermonuclear reaction ensues that disrupts the entire white dwarf, creating an extremely bright optical display that is visible billions of light years away. For a few weeks, SNe Ia are so bright they outshine their host galaxies. “Astrophysicist Subrahmanyan Chandrasekhar, who made enormous contributions to the study of stellar structure and dynamics, showed theoretically in the 1930s that a white dwarf could grow no larger than 1.38 solar masses before it would explode,” says Hoffman. “An exploding white dwarf burning a significant fraction of this mass provides the observed luminosity. A strong argument for using SNe Ia as standard candles is their initial identical configurations.”
A Type Ia supernova explosion generates a light curve, which is quite similar for nearly all SNe Ia. The brightness of a typical Type Ia supernova peaks about 20 days after the explosion, after which the light curve follows almost exactly the radioactive decay curve for nickel-56, an isotope that decays to iron. Although the light curves of all SNe Ia exhibit a similar pattern, they do differ in peak brightness and the brighter ones (with higher peaks) are wider.
Astronomers use a technique that reduces all SNe Ia light curves to a single width–luminosity relationship. In this way, all SNe Ia can be considered as a single parameter family. The width of a light curve indicates the supernova’s intrinsic brightness, according to the rule “broader is brighter.” Also similar are the emission spectra sampled at specific times after peak brightness, lending yet more credence for use of SNe Ia as standard candles. The spectra of SNe Ia indicate the presence of silica, sulfur, calcium, and nickel-56. “About one-third of all iron in the universe, and by inference the hemoglobin in our blood, comes from these exploding white dwarfs,” says Hoffman.
A major goal of the team’s simulations is to better understand the width–luminosity relationship. In particular, scientists want to rule out the possibility that distant SNe Ia, formed in the earliest epochs of the universe, might explode differently from SNe Ia formed more recently. For example, the star-forming environment may have been different in the early universe because far fewer heavy elements were present. “We want to know if variations in initial conditions influence supernova brightness as we go back in time,” says Hoffman. “If so, we would like to account for those variations to more accurately determine the distances, which would be required for the dark energy surveys. By accounting for variations, we can better use SNe Ia as standard candles.”
A Library for Observations
This same library of calculations and simulations also reveals for the first time why the amount of nickel-56 apparently controls the width of the light curve. Scientists have known for some time that the brightness of SNe Ia is determined by the amount of nickel-56 synthesized in the explosion. The decay of this radioactive isotope powers the luminous display. However, the problem of the light-curve width is challenging because the physics involves complex details of light transfer through the time-evolving, wavelength-dependent opacity of the supernova gas.
Because most SNe Ia have been observed from Earth, they have traditionally been seen through the electromagnetic spectrum of visual light. One of the team’s important discoveries is that the width–luminosity relationship is most apparent in the color blue. Scientists are investigating how the color of supernova light evolves over time. Simulation results show that as supernovae iron-group elements (iron, cobalt, and nickel) cool, they exhibit a strong fluorescent effect, whereby light absorbed at blue wavelengths is reemitted at red and infrared wavelengths. The degree of fluorescence thus largely determines the color of the emergent supernova’s spectrum.
Numerical models reveal that the fluorescence process is more efficient in lower temperature gas. Thus, as supernovae expand and cool, their color becomes increasingly red. In the dimmer supernovae, where the gas is initially cooler, this reddening effect evolves more quickly. When observed, the light curves of dim supernovae are found to drop rapidly. These findings are the basis of the width–luminosity relationship. “We are finding that not all SNe Ia are standard candles,” says Hoffman, “but most are within reasonable constraints.”
Many uncertainties remain regarding the evolution of a white dwarf to its full expansion and the observed variation in brightness. Still not clearly understood are the strong convection currents leading to the formation of hot spots that ignite the flame, the propagation of turbulent nuclear burning through the star, the explosion itself, the accompanying synthesis of new elements, and the light curve and spectra emitted from the explosion.
The Livermore simulations, however, did address several outstanding questions surrounding SNe Ia, each of which has eluded a solution for decades. These questions included whether the flame ignites at one or several points and whether it originates at or near the center of the dwarf. The simulations also showed how the flame propagates through the white dwarf, how it is affected by turbulence as it moves outward toward less dense regions, and how it may transition from subsonic to supersonic speeds.
Tools of the Trade
The code SNe is used to study the microphysics of nuclear flames and how these flames interact with turbulence. Another code, MAESTRO, incorporates the changes that occur as the dwarf begins to expand and release heat and as its hot core and less dense outer strata create buoyancy. Simulations from both codes concentrate on early events occurring at subsonic speeds.
When the thermonuclear explosion begins and materials start to move at supersonic speeds, the CASTRO hydrodynamic code takes over. Livermore has a long history of studying hydrodynamic behavior, which plays an important role in exploding nuclear weapons. In addition, the SEDONA code calculates the emergent spectra.
The team successfully combined these codes to model different phases of supernova evolution, ranging from the time when hot spots deep inside a dwarf’s core give birth to the tiny nuclear flame to the weeks following the explosion. Various computer runs included modeling a white dwarf both in a stationary position and rotating at different speeds, the dwarf’s initial flame igniting at different sites, flames burning at different speeds, and turbulence of varying amounts.
In two-dimensional simulations, a detonation with supersonic burning looks bright, and the spectra match most SNe Ia. If supersonic detonation and burning do not occur, the white dwarf explodes but is dimmer. Simulations showed that these spectra do not match those of most SNe Ia.
Scientists have debated whether an asymmetric explosion might modify both the luminosity and spectra of the material streaming out. The team’s simulations of an asymmetric explosion produced the startling result that the peak brightness can vary with the viewing angle by nearly 40 percent. Therefore, if the spectra indicate an asymmetric explosion, one can mathematically correct for the brightness level.
Simulation results show an initially thin nuclear flame—just 0.1 millimeter thick—moving through the dwarf’s core, being pushed around by convective currents, and traveling at less than the local speed of sound. At first the flame remains fairly coherent. As it moves into regions of cooler and less dense fuel, turbulent motion deforms the flame, making it appear thicker and more wrinkled. At a critical point, the flame passes into the so-called distributed burning regime, in which its thickness grows dramatically, its burning rate increases by a factor of five, and its speed transitions to supersonic. Hoffman notes that many aspects of the burning process are similar to those that drive the burning of fuel in an internal combustion engine. In fact, the simulation results have been aided by combustion modeling expertise of scientists at Sandia National Laboratories.
Uncovering the Truth
Many questions remain, perhaps the most significant being why some SNe Ia are good standard candles while others are not. Livermore and other institutions will be examining additional aspects of SNe Ia, such as metallicity, the ratio of carbon and oxygen at the white dwarf’s core, and binary properties that might be the cause of SNe Ia “outliers”—that is, those that do not conform to the width–luminosity relationship.
The findings from current and future observations, together with data from simulations, will also help astronomers and cosmologists determine if dark energy was stronger in the early universe and whether the expansion of the universe will continue to accelerate. “There are many theories and many conflicts but little guidance about the nature of dark energy,” says Hoffman. “We’re hoping our simulations will help uncover the truth.”
—Arnie Heller and Katie Walter
Key Words: CASTRO code, Computational Astrophysical Consortium (CAC), cosmology, dark energy, Grand Challenge computing, hydrodynamics, MAESTRO code, SEDONA code, SNe code, Type Ia supernovae (SNe Ia).
For further information contact Rob Hoffman (925) 424-6411 (firstname.lastname@example.org).
Lawrence Livermore National Laboratory
Privacy & Legal Notice | UCRL-TR-52000-08-9/10 | September 9, 2008