Bolshoi Cosmological Simulation

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search

The Bolshoi simulation, a computer model of the universe run in 2010 on the Pleiades supercomputer at the NASA Ames Research Center, was the most accurate cosmological simulation to that date of the evolution of the large-scale structure of the universe.[1] The Bolshoi simulation used the now-standard ΛCDM (Lambda-CDM) model of the universe and the WMAP five-year and seven-year cosmological parameters from NASA's Wilkinson Microwave Anisotropy Probe team.[2] "The principal purpose of the Bolshoi simulation is to compute and model the evolution of dark matter halos, thereby rendering the invisible visible for astronomers to study, and to predict visible structure that astronomers can seek to observe."[3] “Bolshoi” is a Russian word meaning “big.”

The first two of a series of research papers describing Bolshoi and its implications were published in 2011 in the Astrophysical Journal.[4][5] The first data release of Bolshoi outputs has been made publicly available to the world's astronomers and astrophysicists.[6] The data include output from the Bolshoi simulation and from the BigBolshoi, or MultiDark, simulation of a volume 64 times that of Bolshoi.[7] The Bolshoi-Planck simulation, with the same resolution as Bolshoi, was run in 2013 on the Pleiades supercomputer using the Planck satellite team's cosmological parameters released in March 2013. The Bolshoi-Planck simulation is currently being analyzed in preparation for publication and distribution of its results in 2014.[8][9]

Bolshoi simulations continue to be developed as of 2018.


Joel R. Primack's team at the University of California, Santa Cruz, partnered with Anatoly Klypin's group at New Mexico State University, in Las Cruces[4][5] to run and analyze the Bolshoi simulations. Further analysis and comparison with observations by Risa Wechsler's group at Stanford and others are reflected in the papers based on the Bolshoi simulations.[10]


A successful large-scale simulation of the evolution of galaxies, with results consistent with what is actually seen by astronomers in the night sky, provides evidence that the theoretical underpinnings of the models employed, i.e., the supercomputer implementations ΛCDM, are sound bases for understanding galactic dynamics and the history of the universe, and opens avenues to further research. The Bolshoi Simulation isn't the first large-scale simulation of the universe, but it is the first to rival the extraordinary precision of modern astrophysical observations.[1]

The previous largest and most successful simulation of galactic evolution was the Millennium Simulation Project, led by Volker Springel.[11] Although the success of that project stimulated more than 400 research papers, the Millennium simulations used early WMAP cosmological parameters that have since become obsolete. As a result, they led to some predictions, for example about the distribution of galaxies, that do not match very well with observations. The Bolshoi simulations use the latest cosmological parameters, are higher in resolution, and have been analyzed in greater detail.[10]


The Bolshoi simulation follows the evolving distribution of a statistical ensemble of 8.6 billion particles of dark matter, each of which represents about 200 million solar masses, in a cube of 3-dimensional space about 1 billion light years on edge. Dark matter and dark energy dominate the evolution of the cosmos in this model. The dynamics are modeled with the ΛCDM theory and Albert Einstein's general theory of relativity, with the model including cold dark matter (CDM) and the Λ cosmological constant term simulating the cosmic acceleration referred to as dark energy.

The first 100 million years (Myr) or so of the evolution of the universe after the Big Bang can be derived analytically.[12] The Bolshoi simulation was started at redshift z=80, corresponding to about 20 Myr after the Big Bang. Initial parameters were calculated with linear theory as implemented by the CAMB[13] tools,[14] part of the WMAP website.[15] The tools provide the initial conditions, including a statistical distribution of positions and velocities of the particles in the ensemble, for the much more demanding Bolshoi simulation of the next approximately 13.8 billion years. The experimental volume thus represents a random region of the universe, so comparisons with observations must be statistical.

Key Cosmological Parameters σ8 and ΩM from Observations Compared with Simulations
Two key cosmological parameters, σ8 and ΩM, with values and 1-σ uncertainties from observations and values used in three cosmological simulations. The parameter σ8 represents the amplitude of the fluctuation spectrum on the scale of clusters of galaxies, and the parameter ΩM is the dark + ordinary matter fraction of the cosmic density. The observations represented by the shapes on the figure are from X-ray and gravitational lensing studies of clusters of galaxies. The observations with error bars are from cosmic microwave background data combined with other data from the Wilkinson Microwave Anisotropy Probe (WMAP) five-year (2009), seven-year (2011), and nine-year (2013) publications and the Planck (2013) data release. The simulations are the Millennium I, II, and XXL simulations (which all used the same cosmological parameters consistent with the WMAP first-year data release 2003), and the Bolshoi (2011) and Bolshoi-Planck (2014) simulations.

The Bolshoi simulation employs a version of an adaptive mesh refinement (AMR) algorithm called an adaptive refinement tree (ART), in which a cube in space with more than a predefined density of matter is recursively divided into a mesh of smaller cubes. The subdivision continues to a limiting level, chosen to avoid using too much supercomputer time. Neighboring cubes are not permitted to vary by too many levels, in the case of Bolshoi by more than one level of subdivision, to avoid large discontinuities. The AMR/ART method is well suited to model the increasingly inhomogeneous distribution of matter that evolves as the simulation proceeds. “Once constructed, the mesh, rather than being destroyed at each time step, is promptly adjusted to the evolving particle distribution.”[16] As the Bolshoi simulation ran, the position and velocity of each of the 8.6 billion particles representing dark matter was recorded in 180 snapshots roughly evenly spaced over the simulated 13.8-billion-year run on the Pleiades supercomputer.[4] Each snapshot was then analyzed to find all the dark matter halos and the properties of each (particle membership, location, density distribution, rotation, shape, etc.). All this data was then used to determine the entire growth and merging history of every halo. These results are used in turn to predict where galaxies will form and how they will evolve. How well these predictions correspond to observations provides a measure of the success of the simulation. Other checks were also made.[5]


The Bolshoi simulation is considered to have produced the best approximation to reality so far obtained for so large a volume of space, about 1 billion light years across. “Bolshoi produces a model universe that bears a striking and uncanny resemblance to the real thing. Starting with initial conditions based on the known distribution of matter shortly after the Big Bang, and using Einstein’s general theory of relativity as the ‘rules’ of the simulation, Bolshoi predicts a modern-day universe with galaxies lining up into hundred-million-light-year-long filaments that surround immense voids, forming a cosmic foam-like structure that precisely matches the cosmic web as revealed by deep galaxy studies such as the Sloan Digital Sky Survey. To achieve such a close match, Bolshoi is clearly giving cosmologists a fairly accurate picture of how the universe actually evolved.”[17]


This research was supported by grants from NASA and NSF to Joel Primack and Anatoly Klypin, including massive grants of supercomputer time on the NASA Advanced Supercomputing (NAS) supercomputer Pleiades at NASA Ames Research Center. Hosting of the Bolshoi outputs and analyses at Leibniz Institute for Astrophysics Potsdam (AIP) is partially supported by the MultiDark grant from the Spanish MICINN Programme.[18]

In popular culture[edit]

A visualization from the Bolshoi simulation was narrated in the National Geographic TV special Inside the Milky Way.[7][19] The Icelandic singer-songwriter Björk used footage from the Bolshoi cosmological simulation in the performance of her musical number “Dark Matter” in her Biophilia concert.[20]


  1. ^ a b Primack, Joel R. (1 Oct 2012). "The Cosmological Supercomputer. How the Bolshoi simulation evolves the universe all over again". IEEE Spectrum. IEEE Spectrum. Retrieved 31 Dec 2013. CS1 maint: discouraged parameter (link)
  2. ^ Hayes, Brian. "A Box of Universe". American Scientist. Sigma Xi, The Scientific Research Society. Archived from the original on 10 August 2014. Retrieved 11 Jan 2014. CS1 maint: discouraged parameter (link)
  3. ^ Primack, J.; Bell, T. (July 2012). "Supercomputer modeling is transforming cosmology from a purely observational science into an experimental science" (PDF). University of California High-Performance AstroComputing Center. Sky & Telescope. Retrieved 31 Dec 2013. CS1 maint: discouraged parameter (link)
  4. ^ a b c Klypin, Anatoly A.; Trujillo-Gomez, Sebastian; Primack, Joel (20 Oct 2011). "Dark Matter Halos in the Standard Cosmological Model: Results from the Bolshoi Simulation" (PDF). The Astrophysical Journal. 740 (2): 102. arXiv:1002.3660. Bibcode:2011ApJ...740..102K. doi:10.1088/0004-637X/740/2/102. S2CID 16517863. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  5. ^ a b c Trujillo-Gomez, Sebastian; Klypin, Anatoly; Primack, Joel; Romanowsky, Aaron J. (23 Sep 2011). "Galaxies in ΛCDM with Halo Abundance Matching: Luminosity-Velocity Relation, Baryonic Mass-Velocity Relation, Velocity Function, and Clustering" (PDF). The Astrophysical Journal. 742 (1): 16. arXiv:1005.1289. Bibcode:2011ApJ...742...16T. doi:10.1088/0004-637X/742/1/16. S2CID 53004003. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  6. ^ Kristin Riebe; Adrian M. Partl; Harry Enke; Jaime Forero-Romero; Stefan Gottloeber; Anatoly Klypin; Gerard Lemson; Francisco Prada; Joel R. Primack; Matthias Steinmetz; Victor Turchaninov (August 2013). "The MultiDark Database: Release of the Bolshoi and MultiDark Cosmological Simulations". Astronomische Nachrichten. 334 (7): 691–708. arXiv:1109.0003. Bibcode:2013AN....334..691R. doi:10.1002/asna.201211900. S2CID 16512696. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  7. ^ a b "Introduction: The Bolshoi Simulation at UC-HiPACC". Bolshoi Cosmological Simulations. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  8. ^ Primack, Joel. "Computing the Universe". Los Alamos National Laboratory. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  9. ^ Primack, Joel. "Bolshoi-Planck Cosmological Simulation. Anatoly Klypin & Joel Primack" (PDF). University of California High-Performance AstroComputing Center. p. 25. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  10. ^ a b High-Performance AstroComputing Center, University of California. "Publications". Bolshoi Cosmic Simulation. UC-HiPACC. Retrieved 3 Jan 2014. CS1 maint: discouraged parameter (link)
  11. ^ Boylan-Kolchin, Michael; Volker Springel; Simon D. M. White; Adrian Jenkins; Gerard Lemson (5 Jun 2009). "Resolving cosmic structure formation with the Millennium-II Simulation". Monthly Notices of the Royal Astronomical Society. 398 (3): 1150–1164. arXiv:0903.3041. Bibcode:2009MNRAS.398.1150B. doi:10.1111/j.1365-2966.2009.15191.x. S2CID 9703617.
  12. ^ Loeb, Abraham (2010). How Did the First Stars and Galaxies Form?. Princeton Frontiers in Physics. Princeton, NJ: Princeton University Press. ISBN 9781400834068. Retrieved 3 Jan 2014. CS1 maint: discouraged parameter (link)
  13. ^ Lewis, Antony; Challinor, Anthony. "Code for Anisotropies in the Microwave Background". Antony Lewis. Retrieved 3 Jan 2014. CS1 maint: discouraged parameter (link)
  14. ^ "CAMB Web Interface". National Aeronautics and Space Administration. Goddard Space Flight Center. Retrieved 3 Jan 2014. CS1 maint: discouraged parameter (link)
  15. ^ National Aeronautics and Space Administration, Goddard Space Flight Center. "Wilkinson Microwave Anisotropy Probe". NASA. Retrieved 3 Jan 2014. CS1 maint: discouraged parameter (link)
  16. ^ Kravtsov, Andrey V.; Klypin, Anatoly A.; Khokhlov, Alexei M. (1997). "Adaptive Refinement Tree − A New High-Resolution N-Body Code for Cosmological Simulations" (PDF). The Astrophysical Journal Supplement Series. 111 (1): 73. arXiv:astro-ph/9701195. Bibcode:1997ApJS..111...73K. doi:10.1086/313015. S2CID 14416883.
  17. ^ "Introduction to Interview with Joel Primack". Sky and Telescope. 19 Mar 2012. Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  18. ^ "What was the support for the Bolshoi simulation suite?". Frequently Asked Questions. UC-HiPACC. Retrieved 12 Jan 2014. CS1 maint: discouraged parameter (link)
  19. ^ "Excerpts from "Inside the Milky Way", produced by National Geographic TV". Retrieved 1 Jan 2014. CS1 maint: discouraged parameter (link)
  20. ^ Björk. "Björk - Dark Matter @ Bestival 2011". YouTube. Retrieved 3 Jan 2014. CS1 maint: discouraged parameter (link)

References for figure[edit]

  • Mantz, A., Allen, S. W., Ebeling, H., & Rapetti, D. 2008, MNRAS, 387, 1179
  • Henry, J. P., Evrard, A. E., Hoekstra, H., Babul, A., & Mahdavi, A. 2009, ApJ,691, 1307
  • Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 692, 1060
  • Rozo, E., Rykoff, E. S., Evrard, A., et al. 2009, ApJ, 699, 768

External links[edit]