Scientists release most accurate simulation of the universe to date

September 30, 2011
Bolshoi High Res

High-resolution supercomputer simulation of the large-scale structure of the universe (credit: Stefan Gottlober, AIP)

The Bolshoi supercomputer simulation, announced Thursday, is the most accurate and detailed large cosmological simulation run to date, giving physicists and astronomers a powerful new tool for understanding such cosmic mysteries as galaxy formation, dark matter, and dark energy.

The simulation traces the evolution of the large-scale structure of the universe, including the evolution and distribution of the dark matter halos, in which galaxies coalesced and grew. Initial studies show good agreement between the simulation’s predictions and astronomers’ observations.

“What’s exciting is that we now have this highly accurate simulation that will provide the basis for lots of important new studies in the months and years to come,” said Joel Primack, distinguished professor of physics at the University of California, Santa Cruz.

Primack and Anatoly Klypin, professor of astronomy at New Mexico State University, lead the team that produced the Bolshoi simulation. Klypin wrote the computer code for the simulation, which was run on the Pleiades supercomputer at NASA Ames Research Center.

Primack, who directs the University of California High-Performance Astrocomputing Center (UC-HIPACC), said the initial release of data from the Bolshoi simulation began in early September. “We’ve released a lot of the data so that other astrophysicists can start to use it,” he said. “So far it’s less than one percent of the actual output, because the total output is so huge, but there will be additional releases in the future.”

Based on latest cosmic microwave background studies

The previous benchmark for large-scale cosmological simulations, known as the Millennium Run, has been the basis for some 400 papers since 2005. But the fundamental parameters used as the input for the Millennium Run are now known to be inaccurate.

The Bolshoi simulation is based on the latest cosmic microwave background (WMAP5) parameters, which are consistent with the later WMAP7 results. “The WMAP1 cosmological parameters on which the Millennium simulation is based are now known to be wrong,” Primack said. “Moreover, advances in supercomputer technology allow us to do a much better simulation with higher resolution by almost an order of magnitude. So I expect the Bolshoi simulation will have a big impact on the field.”

The standard explanation for how the universe evolved after the Big Bang is known as the Lambda Cold Dark Matter model, and it is the theoretical basis for the Bolshoi simulation. According to this model, gravity acted initially on slight density fluctuations present shortly after the Big Bang to pull together the first clumps of dark matter. These grew into larger and larger clumps through the hierarchical merging of smaller progenitors.

Although the nature of dark matter remains a mystery, it accounts for about 82 percent of the matter in the universe. As a result, the evolution of structure in the universe has been driven by the gravitational interactions of dark matter. The ordinary matter that forms stars and planets has fallen into the “gravitational wells” created by clumps of dark matter, giving rise to galaxies in the centers of dark matter halos.

A principal purpose of the Bolshoi simulation is to compute and model the evolution of dark matter halos.  A comparison of the Bolshoi predictions with galaxy observations from the Sloan Digital Sky Survey showed very good agreement, according to Primack.

Simulating one billion light-years on a side

The Bolshoi simulation focused on a representative section of the universe, computing the evolution of a cubic volume measuring about one billion light-years on a side and following the interactions of 8.6 billion particles of dark matter. It took 6 million CPU-hours to run the full computation on the Pleiades supercomputer, recently ranked as the seventh fastest supercomputer in the world.

A variant of the Bolshoi simulation, known as BigBolshoi or MultiDark, was run on the same supercomputer with the same number of particles, but this time in a volume 64 times larger. BigBolshoi was run to predict the properties and distribution of galaxy clusters and other very large structures in the universe, as well as to help with dark energy projects such as the Baryon Oscillation Spectroscopic Survey (BOSS).

Another variant, called MiniBolshoi, is currently being run on the Pleiades supercomputer. MiniBolshoi focuses on a smaller portion of the universe and provides even higher resolution than Bolshoi. The Bolshoi simulation and its two variants will be made publicly available to astrophysical researchers worldwide in phases via the MultiDark Database, hosted by the Potsdam Astrophysics Institute in Germany and supported by grants from Spain and Germany.

Primack, Klypin, and their collaborators are continuing to analyze the results of the Bolshoi simulation and submit papers for publication. Among their findings are results showing that the simulation correctly predicts the number of galaxies as bright as the Milky Way that have satellite galaxies as bright as the Milky Way’s major satellites, the Large and Small Magellanic Clouds.

Ref.: A. Klypin, S. Trujillo-Gomez, J. Primack, Halos and galaxies in the standard cosmological model: results from the Bolshoi simulation, arXiv.org, 2011; [arXiv:1002.3660v4]

Ref.: Sebastian Trujillo-Gomez, et al., Galaxies in LCDM with Halo Abundance Matching: luminosity-velocity relation, baryonic mass-velocity relation, velocity function and clustering, arXiv.org, 2011; [arXiv:1005.1289v3]