List of cosmological computation software

Last updated

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

Contents

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[ according to whom? ]:

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]

RAMSES

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]

PKDGRAV3

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. [28]


Relativistic cosmological simulation codes

gevolution

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

A snapshot of AnalyzeThis (CMBEASY) GUI package. The plot shows the marginalize probability distribution from a MCMC chain. Snapshot3-2-.png
A snapshot of AnalyzeThis (CMBEASY) GUI package. The plot shows the marginalize probability distribution from a MCMC chain.

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

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

Related Research Articles

<span class="mw-page-title-main">Cosmic microwave background</span> Trace radiation from the early universe

The cosmic microwave background is microwave radiation that fills all space in the observable universe. It is a remnant that provides an important source of data on the primordial universe. With a standard optical telescope, the background space between stars and galaxies is almost completely dark. However, a sufficiently sensitive radio telescope detects a faint background glow that is almost uniform and is not associated with any star, galaxy, or other object. This glow is strongest in the microwave region of the radio spectrum. The accidental discovery of the CMB in 1965 by American radio astronomers Arno Penzias and Robert Wilson was the culmination of work initiated in the 1940s.

In physics, quintessence is a hypothetical form of dark energy, more precisely a scalar field, postulated as an explanation of the observation of an accelerating rate of expansion of the universe. The first example of this scenario was proposed by Ratra and Peebles (1988) and Wetterich (1988). The concept was expanded to more general types of time-varying dark energy, and the term "quintessence" was first introduced in a 1998 paper by Robert R. Caldwell, Rahul Dave and Paul Steinhardt. It has been proposed by some physicists to be a fifth fundamental force. Quintessence differs from the cosmological constant explanation of dark energy in that it is dynamic; that is, it changes over time, unlike the cosmological constant which, by definition, does not change. Quintessence can be either attractive or repulsive depending on the ratio of its kinetic and potential energy. Those working with this postulate believe that quintessence became repulsive about ten billion years ago, about 3.5 billion years after the Big Bang.

In modern physical cosmology, the cosmological principle is the notion that the spatial distribution of matter in the universe is uniformly isotropic and homogeneous when viewed on a large enough scale, since the forces are expected to act equally throughout the universe on a large scale, and should, therefore, produce no observable inequalities in the large-scale structuring over the course of evolution of the matter field that was initially laid down by the Big Bang.

The Sunyaev–Zeldovich effect is the spectral distortion of the cosmic microwave background (CMB) through inverse Compton scattering by high-energy electrons in galaxy clusters, in which the low-energy CMB photons receive an average energy boost during collision with the high-energy cluster electrons. Observed distortions of the cosmic microwave background spectrum are used to detect the disturbance of density in the universe. Using the Sunyaev–Zeldovich effect, dense clusters of galaxies have been observed.

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">Very Small Array</span> Radio telescope in the Canary Islands

The Very Small Array (VSA) was a 14-element interferometric radio telescope operating between 26 and 36 GHz that is used to study the cosmic microwave background radiation. It was a collaboration between the University of Cambridge, University of Manchester and the Instituto de Astrofisica de Canarias (Tenerife), and was located at the Observatorio del Teide on Tenerife. The array was built at the Mullard Radio Astronomy Observatory by the Cavendish Astrophysics Group and Jodrell Bank Observatory, and was funded by PPARC. The design was strongly based on the Cosmic Anisotropy Telescope.

In physical cosmology, CMBFAST is a computer code, written by Uroš Seljak and Matias Zaldarriaga, for computing the anisotropy of the cosmic microwave background. It was 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.

<span class="mw-page-title-main">GADGET</span> Computer software for cosmological simulations

GADGET is free software for cosmological N-body/SPH simulations written by Volker Springel at the Max Planck Institute for Astrophysics. The name is an acronym of "GAlaxies with Dark matter and Gas intEracT". It is released under the GNU GPL. It can be used to study for example galaxy formation and dark matter.

<span class="mw-page-title-main">Archeops</span>

Archeops was a balloon-borne instrument dedicated to measuring the Cosmic microwave background (CMB) temperature anisotropies. The study of this radiation is essential to obtain precise information on the evolution of the Universe: density, Hubble constant, age of the Universe, etc. To achieve this goal, measurements were done with devices cooled down at 100mK temperature placed at the focus of a warm telescope. To avoid atmospheric disturbance the whole apparatus is placed on a gondola below a helium balloon that reaches 40 km altitude.

<span class="mw-page-title-main">CMB cold spot</span> Region in space

The CMB Cold Spot or WMAP Cold Spot is a region of the sky seen in microwaves that has been found to be unusually large and cold relative to the expected properties of the cosmic microwave background radiation (CMBR). The "Cold Spot" is approximately 70 µK (0.00007 K) colder than the average CMB temperature, whereas the root mean square of typical temperature variations is only 18 µK. At some points, the "cold spot" is 140 µK colder than the average CMB temperature.

Conformal cyclic cosmology (CCC) is a cosmological model in the framework of general relativity and proposed by theoretical physicist Roger Penrose. In CCC, the universe iterates through infinite cycles, with the future timelike infinity of each previous iteration being identified with the Big Bang singularity of the next. Penrose popularized this theory in his 2010 book Cycles of Time: An Extraordinary New View of the Universe.

<span class="mw-page-title-main">Uroš Seljak</span> Slovenian cosmologist

Uroš Seljak is a Slovenian cosmologist and a professor of astronomy and physics at University of California, Berkeley. He is particularly well-known for his research in cosmology and approximate Bayesian statistical methods.

In cosmology, decoupling is a period in the development of the universe when different types of particles fall out of thermal equilibrium with each other. This occurs as a result of the expansion of the universe, as their interaction rates decrease up to this critical point. The two verified instances of decoupling since the Big Bang which are most often discussed are photon decoupling and neutrino decoupling, as these led to the cosmic microwave background and cosmic neutrino background, respectively.

Minicharged particles are a proposed type of subatomic particle. They are charged, but with a tiny fraction of the charge of the electron. They weakly interact with matter. Minicharged particles are not part of the Standard Model. One proposal to detect them involved photons tunneling through an opaque barrier in the presence of a perpendicular magnetic field, the rationale being that a pair of oppositely charged minicharged particles are produced that curve in opposite directions, and recombine on the other side of the barrier reproducing the photon again.

<span class="mw-page-title-main">Cosmology Large Angular Scale Surveyor</span> Microwave telescope array in Chile

The Cosmology Large Angular Scale Surveyor (CLASS) is an array of microwave telescopes at a high-altitude site in the Atacama Desert of Chile as part of the Parque Astronómico de Atacama. The CLASS experiment aims to improve our understanding of cosmic dawn when the first stars turned on, test the theory of cosmic inflation, and distinguish between inflationary models of the very early universe by making precise measurements of the polarization of the Cosmic Microwave Background (CMB) over 65% of the sky at multiple frequencies in the microwave region of the electromagnetic spectrum.

The Illustris project is an ongoing series of astrophysical simulations run by an international collaboration of scientists. The aim was 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">Atacama B-Mode Search</span>

The Atacama B-Mode Search (ABS) was an experiment to test the theory of cosmic inflation and distinguish between inflationary models of the very early universe by making precise measurements of the polarization of the Cosmic Microwave Background (CMB). ABS was located at a high-altitude site in the Atacama Desert of Chile as part of the Parque Astronómico de Atacama. ABS began observations in February 2012 and completed observations in October 2014.

The "axis of evil" is a name given to the apparent correlation between the plane of the Solar System and aspects of the cosmic microwave background (CMB). It gives the plane of the Solar System and hence the location of Earth a greater significance than might be expected by chance – a result which has been claimed to be evidence of a departure from the Copernican principle as assumed in the concordance model.

In cosmological inflation, within the slow-roll paradigm, the Lyth argument places a theoretical upper bound on the amount of gravitational waves produced during inflation, given the amount of departure from the homogeneity of the cosmic microwave background (CMB).

<span class="mw-page-title-main">Cosmic microwave background spectral distortions</span> Fluctuations in the energy spectrum of the microwave background

CMB spectral distortions are tiny departures of the average cosmic microwave background (CMB) frequency spectrum from the predictions given by a perfect black body. They can be produced by a number of standard and non-standard processes occurring at the early stages of cosmic history, and therefore allow us to probe the standard picture of cosmology. Importantly, the CMB frequency spectrum and its distortions should not be confused with the CMB anisotropy power spectrum, which relates to spatial fluctuations of the CMB temperature in different directions of the sky.