List of cosmological computation software explained

This List of Cosmological Computation Software catalogs the tools and programs used by scientists in cosmological research.

In the past few decades, the accelerating technological evolution has profoundly enhanced astronomical instrumentation, enabling more precise observations and expanding the breadth and depth of data collection by several orders of magnitude. Simultaneously, the exponential growth in computational power has enabled the creation of computer simulations that reveal details with unprecedented resolution and accuracy. For performing computer simulations of the cosmos and analyzing data from both cosmological experiments and simulations, many advanced methods and computational software codes are developed every year. These codes are widely used by researchers all across the globe, in all various fields and topics of cosmology.

The computational software used in cosmology can be classified into the following major classes:

Newtonian cosmological simulation codes

GADGET

GADGET, named "GAlaxies with Dark matter and Gas intEracT" is a code written in C++ for cosmological N-body/Smoothed-particle hydrodynamics (SPH) simulations on massively parallel computers with distributed memory.[15] Its first version was developed by German astrophysicist, Volker Springel and was published in 2000.[16] It was followed by two more official public versions, with GADGET-2[17] [18] released in 2005 and GADGET-4[19] [20] released in 2020, which is the most recent public version of the software suite currently. GADGET is capable to address a wide array of astrophysically interesting problems, e.g. the dynamics of the gaseous intergalactic medium, star formation and its regulation by feedback processes, colliding and merging galaxies, as well as the formation of large-scale structure in the Universe.

AREPO

AREPO[21] [22] is a massively parallel code for gravitational N-body systems, hydrodynamics and magnetohydrodynamics (MHD). It is named after the enigmatic word AREPO in the Latin palindromic sentence "sator arepo tenet opera rotas", the Sator Square. The first version of AREPO was written and published by Volker Springel in 2010, with further development by Rüdiger Pakmor and contributions by many other authors. The Arepo code utilizes an unstructured Voronoi-mesh and was designed to blend the benefits of finite-volume hydrodynamics and SPH. Primarily optimized for cosmological simulations, especially galaxy formation, Arepo supports a high dynamic range in space and time.[23]

GIZMO

GIZMO[24] is a flexible, massively parallel, multi-physics simulation code, written in ANSI C by Philip F. Hopkins. The code offers diverse methods to solve fluid equations. It also introduces novel methods, which optimize the resolution of simulations and minimize common errors found in previous methods that limited the accuracy of prior solvers. Originating from GADGET (hence the name "GIZMO", a play on words), the code maintains compatibility in naming/use conventions as well as input/output, making it user-friendly for those familiar with GADGET.[25]

StePS

StePS,[26] [27] which stands for "STEreographically Projected cosmological Simulations" is a freely available code that implements a novel N-body simulation method that models an infinite universe within a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Unlike traditional methods, which use unrealistic periodic boundary conditions for numerical simplicity, StePS offers a more observation-aligned approach. This technique enables detailed simulations of an infinite universe using less memory and provides results that are more in line with the observed universe geometry and topology.

Relativistic cosmological simulation codes

CosmoGRaPH

CosmoGRaPH (Cosmological General Relativity And (Perfect fluid | Particle) Hydrodynamics) is a C++ code used to explore cosmological problems in a fully general relativistic setting. It was developed by James Mertens and Chi Tian and was published in 2016. The code implements various novel methods for numerically solving the Einstein field equations, including an N-body solver, full AMR capabilities via SAMRAI, and raytracing.

Cosmological Boltzmann codes

CMBFAST

CMBFAST is a computer code, developed by Uroš Seljak and Matias Zaldarriaga (based on a Boltzmann code written by Edmund Bertschinger, Chung-Pei Ma and Paul Bode) for computing the power spectrum of the cosmic microwave background anisotropy. It is the first efficient program to do so, reducing the time taken to compute the anisotropy from several days to a few minutes by using a novel semi-analytic line-of-sight approach.

CAMB

Code for Anisotropies in the Microwave Background by Antony Lewis and Anthony Challinor. The code was originally based on CMBFAST. Later several developments are made to make it a faster and more accurate and compatible with the present research. The code is written in an object oriented manner to make it more user friendly.

CMBEASY

CMBEASY is a software package written by Michael Doran, Georg Robbers and Christian M. Müller. The code is based on the CMBFAST package. CMBEASY is fully object oriented C++. This considerably simplifies manipulations and extensions of the CMBFAST code. In addition, a powerful Spline class can be used to easily store and visualize data. Many features of the CMBEASY package are also accessible via a graphical user interface. This may be helpful for gaining intuition, as well as for instruction purposes.

CLASS

The purpose of the Cosmic Linear Anisotropy Solving System is to simulate the evolution of linear perturbations in the universe and to compute CMB and large scale structure observables. CLASS is written in plain C to achieve high performance, yet its modular structure emulates the architecture and philosophy of classes in object-oriented languages for enhanced readability and modularity. The name "CLASS" also derives from its object-oriented style, mimicking the notion of a class.

Parameter estimation packages

AnalizeThis

AnalizeThis is a parameter estimation package used by cosmologists. It comes with the CMBEASY package. The code is written in C++ and uses the global metropolis algorithm for estimation of cosmological parameters. The code was developed by Michael Doran, for parameter estimation using WMAP-5 likelihood. However, the code was not updated after 2008 for the new CMB experiments. Hence this package is currently not in use by the CMB research community. The package comes up with a nice GUI.

CosmoMC

CosmoMC is a Fortran 2003 Markov chain Monte Carlo (MCMC) engine for exploring cosmological parameter space. The code does brute force (but accurate) theoretical matter power spectrum and Cl calculations using CAMB. CosmoMC uses a simple local Metropolis algorithm along with an optimized fast-slow sampling method. This fast-slow sampling method provides faster convergence for the cases with many nuisance parameters like Planck. CosmoMC package also provides subroutines for post processing and plotting of the data.

CosmoMC was written by Antony Lewis in 2002 and later several versions are developed to keep the code up-to date with different cosmological experiments. It is presently the most used cosmological parameter estimation code.

SCoPE

SCoPE/Slick Cosmological Parameter Estimator is a newly developed cosmological MCMC package written by Santanu Das in C language. Apart from standard global metropolis algorithm the code uses three unique technique named as 'delayed rejection' that increases the acceptance rate of a chain, 'pre-fetching' that helps an individual chain to run on parallel CPUs and 'inter-chain covariance update' that prevents clustering of the chains allowing faster and better mixing of the chains. The code is capable of faster computation of cosmological parameters from WMAP and Planck data.

Other packages

Likelihood software packages

Different cosmology experiments, in particular the CMB experiments like WMAP and Planck measures the temperature fluctuations in the CMB sky and then measure the CMB power spectrum from the observed skymap. But for parameter estimation the χ² is required. Therefore, all these CMB experiments comes up with its own likelihood software.

See also

Notes and References

  1. Seljak. Uros. Zaldarriaga, Matias . A Line of Sight Approach to Cosmic Microwave Background Anisotropies. Astrophysical Journal. 1996. 469. 437–444. 10.1086/177793 . 1996ApJ...469..437S. astro-ph/9603033 . 3015599.
  2. Zaldarriaga. Matias. Uros Seljak . Edmund Bertschinger . Integral Solution for the Microwave Background Anisotropies in Non-flat Universes. Astrophysical Journal. 1998. 494. 2. 491–502. 10.1086/305223. 1998ApJ...494..491Z. astro-ph/9704265 . 15966880.
  3. Web site: Seljak. U., & Zaldarriaga, M.. CMBFAST. 2014-02-22. 2016-04-06. https://web.archive.org/web/20160406121934/http://lambda.gsfc.nasa.gov/toolbox/tb_cmbfast_ov.cfm. dead.
  4. Lewis. Antony. Challinor, Anthony . CAMB: Code for Anisotropies in the Microwave Background. Astrophysics Source Code Library . ascl:1102.026 . 2011ascl.soft02026L. . 2011.
  5. Web site: Doran. Michael. CMBEASY. 2014-02-22. https://web.archive.org/web/20140118063521/http://www.thphys.uni-heidelberg.de/~robbers/cmbeasy/. 2014-01-18. dead.
  6. Doran. Michael. CMBEASY:: an Object Oriented Code for the Cosmic Microwave Background. Journal of Cosmology and Astroparticle Physics. 27 Apr 2006. 0510. 10. 011. 10.1088/1475-7516/2005/10/011 . astro-ph/0302138 . 2005JCAP...10..011D . 5451633. Submitted manuscript.
  7. Blas. D.. J. Lesgourgues . T. Tram . CLASS II: Approximation schemes. Journal of Cosmology and Astroparticle Physics. 2011. 1107. 7. 034. 10.1088/1475-7516/2011/07/034. 1104.2933. 2011JCAP...07..034B. 53490516.
  8. Lesgourgues. J. CLASS I: Overview. 1104.2932. astro-ph.IM. 2011.
  9. Web site: CLASS, the Cosmic Linear Anisotropy Solving System . Lesgourgues . J. . https://web.archive.org/web/20230715215459/http://class-code.net/ . July 15, 2023 . dead.
  10. Das . Santanu . Phan . Anh . Cosmic Microwave Background Anisotropy numerical solution (CMBAns). Part I. An introduction to Cl calculation . Journal of Cosmology and Astroparticle Physics . 5 May 2020 . 2020 . 5 . 006 . 10.1088/1475-7516/2020/05/006 . 1910.00725 . 2020JCAP...05..006D . 203626962 .
  11. Lewis. Antony. Sarah Bridle . Cosmological parameters from CMB and other data: a Monte-Carlo approach. Physical Review D . 103511. 2002. 10.1103/PhysRevD.66.103511. astro-ph/0205436 . 2002PhRvD..66j3511L. 66. 10. 55316758.
  12. Lewis. Antony. Efficient sampling of fast and slow cosmological parameters. Physical Review D. 2013. 87. 10. 103529. 10.1103/PhysRevD.87.103529. 1304.4473 . 2013PhRvD..87j3529L . 119259816.
  13. Doran. Michael. Christian M. Mueller . Analyze This! A Cosmological Constraint Package for CMBEASY. Journal of Cosmology and Astroparticle Physics. 2004. 0409. 3. 003. 10.1088/1475-7516/2004/09/003 . astro-ph/0311311 . 2004JCAP...09..003D . 119333027.
  14. Das. Santanu. Tarun Souradeep . SCoPE: An efficient method of Cosmological Parameter Estimation. Journal of Cosmology and Astroparticle Physics. 2014. 1407. 18. 018. 10.1088/1475-7516/2014/07/018 . 1403.1271 . 2014JCAP...07..018D . 119233297.
  15. Web site: Cosmological simulations with GADGET .
  16. 10.1016/S1384-1076(01)00042-2 . GADGET: A code for collisionless and gasdynamical cosmological simulations . 2001 . Springel . Volker . Yoshida . Naoki . White . Simon D.M. . New Astronomy . 6 . 2 . 79–117 . astro-ph/0003162 . 2001NewA....6...79S . 5879269 .
  17. 10.1111/j.1365-2966.2005.09655.x . The cosmological simulation code gadget-2 . 2005 . Springel . Volker . Monthly Notices of the Royal Astronomical Society . 364 . 4 . 1105–1134 . astro-ph/0505010 . 2005MNRAS.364.1105S . 16378825 .
  18. Web site: ASCL.net - GADGET-2: A Code for Cosmological Simulations of Structure Formation . 2023-10-25 . ascl.net.
  19. 10.1093/mnras/stab1855 . Simulating cosmic structure formation with the gadget-4 code . 2021 . Springel . Volker . Pakmor . Rüdiger . Zier . Oliver . Reinecke . Martin . Monthly Notices of the Royal Astronomical Society . 506 . 2 . 2871–2949 . 2010.03567 .
  20. Web site: ASCL.net - GADGET-4: Parallel cosmological N-body and SPH code . 2023-10-25 . ascl.net.
  21. 10.1111/j.1365-2966.2009.15715.x . E pur si muove: Galilean-invariant cosmological hydrodynamical simulations on a moving mesh . 2010 . Springel . Volker . Monthly Notices of the Royal Astronomical Society . 401 . 2 . 791–851 . 0901.4107 . 2010MNRAS.401..791S . 119241866 .
  22. Web site: ASCL.net - AREPO: Cosmological magnetohydrodynamical moving-mesh simulation code . 2023-10-25 . ascl.net.
  23. Web site: Arepo documentation — Arepo 1.0 documentation .
  24. Web site: ASCL.net - GIZMO: Multi-method magneto-hydrodynamics+gravity code . 2023-10-25 . ascl.net.
  25. Web site: GIZMO Phil Hopkins' Research Group . 2023-10-25 . www.tapir.caltech.edu.
  26. 10.1016/j.ascom.2019.100303 . StePS: A multi-GPU cosmological N-body Code for compactified simulations . 2019 . Rácz . G. . Szapudi . I. . Dobos . L. . Csabai . I. . Szalay . A.S. . Astronomy and Computing . 28 . 100303 . 1811.05903 . 2019A&C....2800303R . 84842307 .
  27. Web site: ASCL.net - StePS: Stereographically Projected Cosmological Simulations . 2023-10-25 . ascl.net.
  28. Das . Santanu . SIToolBox: A package for Bayesian estimation of the isotropy violation in the CMB sky . MNRAS . November 2019 . 489 . 4 . 5889–5899 . 10.1093/mnras/stz2542 . free . 1810.09470 .
  29. Web site: C3 Research — CMB.
  30. Eriksen. H. K.. Jewell. J. B.. Dickinson. C.. Banday. A. J.. Górski. K. M.. Lawrence. C. R.. March 2008. Joint Bayesian Component Separation and CMB Power Spectrum Estimation. The Astrophysical Journal. 676. 1. 10–32. 10.1086/525277. 0709.1058. 2008ApJ...676...10E. 8662228.