Bolshoi Cosmological Simulation

Last updated

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.”

Contents

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.

Contributors

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 University and others are reflected in the papers based on the Bolshoi simulations. [10]

Rationale

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]

Methods

The Bolshoi simulation follows the evolving distribution of a statistical ensemble of 8.6 billion particles of dark matter, each of which represents about 100 million solar masses, [1] 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.

Two key cosmological parameters, s8 and OM, with values and 1-s uncertainties from observations and values used in three cosmological simulations. The parameter s8 represents the amplitude of the fluctuation spectrum on the scale of clusters of galaxies, and the parameter OM 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. Key Cosmological Parameters s8 and OM from Observations Compared with Simulations.jpg
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]

Results

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] The Bolshoi simulation found that the Sheth–Tormen approximation overpredicts the abundance of halos by a factor of for redshifts . [4]

Support

This research was supported by grants from NASA and the National Science Foundation (U.S.) 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]

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]

Related Research Articles

<span class="mw-page-title-main">Physical cosmology</span> Branch of cosmology which studies mathematical models of the universe

Physical cosmology is a branch of cosmology concerned with the study of cosmological models. A cosmological model, or simply cosmology, provides a description of the largest-scale structures and dynamics of the universe and allows study of fundamental questions about its origin, structure, evolution, and ultimate fate. Cosmology as a science originated with the Copernican principle, which implies that celestial bodies obey identical physical laws to those on Earth, and Newtonian mechanics, which first allowed those physical laws to be understood.

<span class="mw-page-title-main">Dark matter</span> Concept in cosmology

In astronomy, dark matter is a hypothetical form of matter that appears not to interact with light or the electromagnetic field. Dark matter is implied by gravitational effects which cannot be explained by general relativity unless more matter is present than can be seen. Such effects occur in the context of formation and evolution of galaxies, gravitational lensing, the observable universe's current structure, mass position in galactic collisions, the motion of galaxies within galaxy clusters, and cosmic microwave background anisotropies.

<span class="mw-page-title-main">Galaxy formation and evolution</span>

The study of galaxy formation and evolution is concerned with the processes that formed a heterogeneous universe from a homogeneous beginning, the formation of the first galaxies, the way galaxies change over time, and the processes that have generated the variety of structures observed in nearby galaxies. Galaxy formation is hypothesized to occur from structure formation theories, as a result of tiny quantum fluctuations in the aftermath of the Big Bang. The simplest model in general agreement with observed phenomena is the Lambda-CDM model—that is, clustering and merging allows galaxies to accumulate mass, determining both their shape and structure. Hydrodynamics simulation, which simulates both baryons and dark matter, is widely used to study galaxy formation and evolution.

In cosmology and physics, cold dark matter (CDM) is a hypothetical type of dark matter. According to the current standard model of cosmology, Lambda-CDM model, approximately 27% of the universe is dark matter and 68% is dark energy, with only a small fraction being the ordinary baryonic matter that composes stars, planets, and living organisms. Cold refers to the fact that the dark matter moves slowly compared to the speed of light, giving it a vanishing equation of state. Dark indicates that it interacts very weakly with ordinary matter and electromagnetic radiation. Proposed candidates for CDM include weakly interacting massive particles, primordial black holes, and axions.

<span class="mw-page-title-main">Observable universe</span> All of space observable from the Earth at the present

The observable universe is a ball-shaped region of the universe consisting of all matter that can be observed from Earth or its space-based telescopes and exploratory probes at the present time; the electromagnetic radiation from these objects has had time to reach the Solar System and Earth since the beginning of the cosmological expansion. Initially, it was estimated that there may be 2 trillion galaxies in the observable universe. That number was reduced in 2021 to several hundred billion based on data from New Horizons. Assuming the universe is isotropic, the distance to the edge of the observable universe is roughly the same in every direction. That is, the observable universe is a spherical region centered on the observer. Every location in the universe has its own observable universe, which may or may not overlap with the one centered on Earth.

<span class="mw-page-title-main">Age of the universe</span> Time elapsed since the Big Bang

In physical cosmology, the age of the universe is the time elapsed since the Big Bang. Astronomers have derived two different measurements of the age of the universe: a measurement based on direct observations of an early state of the universe, which indicate an age of 13.787±0.020 billion years as interpreted with the Lambda-CDM concordance model as of 2021; and a measurement based on the observations of the local, modern universe, which suggest a younger age. The uncertainty of the first kind of measurement has been narrowed down to 20 million years, based on a number of studies that all show similar figures for the age. These studies include researches of the microwave background radiation by the Planck spacecraft, the Wilkinson Microwave Anisotropy Probe and other space probes. Measurements of the cosmic background radiation give the cooling time of the universe since the Big Bang, and measurements of the expansion rate of the universe can be used to calculate its approximate age by extrapolating backwards in time. The range of the estimate is also within the range of the estimate for the oldest observed star in the universe.

The cuspy halo problem is a discrepancy between the inferred dark matter density profiles of low-mass galaxies and the density profiles predicted by cosmological N-body simulations. Nearly all simulations form dark matter halos which have "cuspy" dark matter distributions, with density increasing steeply at small radii, while the rotation curves of most observed dwarf galaxies suggest that they have flat central dark matter density profiles ("cores").

<span class="mw-page-title-main">Lambda-CDM model</span> Model of Big Bang cosmology

The Lambda-CDM, Lambda cold dark matter, or ΛCDM model is a mathematical model of the Big Bang theory with three major components:

  1. a cosmological constant, denoted by lambda (Λ), associated with dark energy
  2. the postulated cold dark matter, denoted by CDM
  3. ordinary matter
<span class="mw-page-title-main">Dark matter halo</span> Theoretical cosmological structure

In modern models of physical cosmology, a dark matter halo is a basic unit of cosmological structure. It is a hypothetical region that has decoupled from cosmic expansion and contains gravitationally bound matter. A single dark matter halo may contain multiple virialized clumps of dark matter bound together by gravity, known as subhalos. Modern cosmological models, such as ΛCDM, propose that dark matter halos and subhalos may contain galaxies. The dark matter halo of a galaxy envelops the galactic disc and extends well beyond the edge of the visible galaxy. Thought to consist of dark matter, halos have not been observed directly. Their existence is inferred through observations of their effects on the motions of stars and gas in galaxies and gravitational lensing. Dark matter halos play a key role in current models of galaxy formation and evolution. Theories that attempt to explain the nature of dark matter halos with varying degrees of success include cold dark matter (CDM), warm dark matter, and massive compact halo objects (MACHOs).

<span class="mw-page-title-main">Structure formation</span> Formation of galaxies, galaxy clusters and larger structures from small early density fluctuations

In physical cosmology, structure formation describes the creation of galaxies, galaxy clusters, and larger structures starting from small fluctuations in mass density resulting from processes that created matter. The universe, as is now known from observations of the cosmic microwave background radiation, began in a hot, dense, nearly uniform state approximately 13.8 billion years ago. However, looking at the night sky today, structures on all scales can be seen, from stars and planets to galaxies. On even larger scales, galaxy clusters and sheet-like structures of galaxies are separated by enormous voids containing few galaxies. Structure formation models gravitational instability of small ripples in mass density to predict these shapes, confirming the consistency of the physical model.

The Millennium Run, or Millennium Simulation is a computer N-body simulation used to investigate how the distribution of matter in the Universe has evolved over time, in particular, how the observed population of galaxies was formed. It is used by scientists working in physical cosmology to compare observations with theoretical predictions.

Joel R. Primack is an American physicist. He is a professor of physics and astrophysics at the University of California, Santa Cruz, and is a member of the Santa Cruz Institute for Particle Physics.

The Navarro–Frenk–White (NFW) profile is a spatial mass distribution of dark matter fitted to dark matter halos identified in N-body simulations by Julio Navarro, Carlos Frenk and Simon White. The NFW profile is one of the most commonly used model profiles for dark matter halos.

<span class="mw-page-title-main">Simon White</span> British astronomer

Simon David Manton White, FRS, is a British-German astrophysicist. He was one of directors at the Max Planck Institute for Astrophysics before his retirement in late 2019.

<span class="mw-page-title-main">Institute for Computational Cosmology</span> Research institute at Durham University

The Institute for Computational Cosmology (ICC) is a research institute at Durham University, England. It was founded in November 2002 as part of the Ogden Centre for Fundamental Physics, which also includes the Institute for Particle Physics Phenomenology (IPPP). The ICC's primary mission is to advance fundamental knowledge in cosmology. Topics of active research include: the nature of dark matter and dark energy, the evolution of cosmic structure, the formation of galaxies, and the determination of fundamental parameters.

The University of California High-Performance AstroComputing Center (UC-HiPACC) based at the University of California at Santa Cruz (UCSC) is a consortium of nine University of California campuses and three Department of Energy laboratories. The consortium's goal is to support and facilitate original research and education in computational astrophysics and to engage in public outreach and education.

<span class="mw-page-title-main">Illustris project</span> Computer-simulated universes

The Illustris project is an ongoing series of astrophysical simulations run by an international collaboration of scientists. The aim is to study the processes of galaxy formation and evolution in the universe with a comprehensive physical model. Early results were described in a number of publications following widespread press coverage. The project publicly released all data produced by the simulations in April, 2015. Key developers of the Illustris simulation have been Volker Springel and Mark Vogelsberger. The Illustris simulation framework and galaxy formation model has been used for a wide range of spin-off projects, starting with Auriga and IllustrisTNG followed by Thesan (2021), MillenniumTNG (2022) and TNG-Cluster (2023).

<span class="mw-page-title-main">Ben Moore (astrophysicist)</span> American professor of astrophysics

Ben Moore is an English professor of astrophysics, author, musician, and director of the Center for Theoretical Astrophysics and Cosmology at the University of Zürich. His research is focussed on cosmology, gravity, astroparticle physics, and planet formation. He has authored in excess of 200 scientific papers on the origin of planets and galaxies, as well as dark matter and dark energy. In his research, he simulates the universe using custom-built supercomputers.

<span class="mw-page-title-main">UniverseMachine</span> Computer simulated universes

The UniverseMachine is a project carrying out astrophysical supercomputer simulations of various models of possible universes, created by astronomer Peter Behroozi and his research team at the Steward Observatory and the University of Arizona. Numerous universes with different physical characteristics may be simulated in order to develop insights into the possible beginning and evolution of our universe. A major objective is to better understand the role of dark matter in the development of the universe. According to Behroozi, "On the computer, we can create many different universes and compare them to the actual one, and that lets us infer which rules lead to the one we see."

The Sheth–Tormen approximation is a halo mass function.

References

  1. 1 2 3 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.
  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 January 2014.
  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.
  4. 1 2 3 4 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 January 2014.
  5. 1 2 3 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 January 2014.
  6. Riebe, Kristin; Partl, Adrian M.; Enke, Harry; Forero-Romero, Jaime; Gottloeber, Stefan; Klypin, Anatoly; Lemson, Gerard; Prada, Francisco; Primack, Joel R.; Steinmetz, Matthias; Turchaninov, Victor (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 January 2014.
  7. 1 2 "Introduction: The Bolshoi Simulation at UC-HiPACC". Bolshoi Cosmological Simulations. Retrieved 1 January 2014.
  8. Primack, Joel. "Computing the Universe". Los Alamos National Laboratory . Retrieved 1 January 2014.
  9. Primack, Joel. "Bolshoi-Planck Cosmological Simulation. Anatoly Klypin & Joel Primack" (PDF). University of California High-Performance AstroComputing Center. p. 25. Retrieved 1 January 2014.
  10. 1 2 High-Performance AstroComputing Center, University of California. "Publications". Bolshoi Cosmic Simulation. University of California-HiPACC. Retrieved 3 January 2014.
  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, New Jersey: Princeton University Press. ISBN   9781400834068 . Retrieved 3 January 2014.
  13. Lewis, Antony; Challinor, Anthony. "Code for Anisotropies in the Microwave Background". Antony Lewis. Retrieved 3 January 2014.
  14. "CAMB Web Interface". National Aeronautics and Space Administration. Goddard Space Flight Center. Retrieved 3 January 2014.
  15. National Aeronautics and Space Administration, Goddard Space Flight Center. "Wilkinson Microwave Anisotropy Probe". NASA . Retrieved 3 January 2014.
  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 January 2014.
  18. "What was the support for the Bolshoi simulation suite?". Frequently Asked Questions. University of California-HiPACC. Retrieved 12 January 2014.
  19. "Excerpts from "Inside the Milky Way", produced by National Geographic TV" . Retrieved 1 January 2014.
  20. Björk. "Björk – Dark Matter @ Bestival 2011". YouTube . Retrieved 3 January 2014.

References for figure