AREPO

AREPO
Name	AREPO
Title	AREPO
Developer	Volker Springel and collaborators
Released	2009
Programming language	C, C++
Platform	Linux, macOS
License	Proprietary / Academic

AREPO.

AREPO is a computational code for astrophysical and cosmological simulations that uses a moving unstructured mesh to solve the equations of hydrodynamics coupled to gravity, magnetohydrodynamics, and radiative processes. It was developed to bridge techniques used in grid-based Eulerian codes and particle-based Lagrangian methods, and it has been employed by research groups studying galaxy formation, large-scale structure, galaxy clusters, star formation, and feedback. The code has influenced work at institutions and collaborations including the Max Planck Society, Harvard University, Princeton University, the Kavli Institute, and the Center for Computational Astrophysics.

Overview

AREPO implements a finite-volume Godunov scheme on a Voronoi tessellation that moves with the flow, combining ideas from the Tessellation literature, the Voronoi diagram formalism, and shock-capturing schemes developed in the context of the Riemann problem. It couples hydrodynamics with self-gravity via tree and particle-mesh methods similar to those used in N-body simulation toolkits and incorporates subgrid physics modules inspired by implementations from projects like GADGET-2, ENZO, RAMSES, FLASH (software), and ZEUS (code). The code architecture supports adaptive resolution that follows mass elements, enabling studies comparable to work carried out with the Millennium simulation, the Illustris project, and the EAGLE (simulation) collaboration. AREPO's numerical approach is related to techniques used by researchers at the Max Planck Institute for Astrophysics, Harvard-Smithsonian Center for Astrophysics, and the Space Telescope Science Institute.

Development and Algorithms

AREPO originated from work by Volker Springel and colleagues at institutions including the Max Planck Society and the Heidelberg Institute for Theoretical Studies. Its core algorithm constructs a moving Voronoi mesh from a set of mesh-generating points, then solves local Riemann problems across faces using approximate solvers akin to those implemented in HLLC and other shock-capturing frameworks. Gravity is computed via a hierarchical multipole Barnes–Hut tree and particle-mesh methods that use Fourier techniques similar to those in the Particle-mesh method and hybrid algorithms employed in the Gadget series. AREPO supports magnetohydrodynamics via constrained transport and divergence-cleaning strategies comparable to approaches in Athena (code) and PLUTO (code), and integrates cooling, chemistry, and radiative transfer modules following methods used by teams behind CLOUDY, CHIANTI (database), and radiative transfer codes like RADMC-3D. Parallelization uses domain decomposition, MPI, and load-balancing strategies derived from best practices at institutions such as Lawrence Livermore National Laboratory, Argonne National Laboratory, and the Oak Ridge National Laboratory.

Applications and Simulations

AREPO has been applied to cosmological volumes, zoom simulations of galaxy halos, isolated galaxy models, cluster physics, interstellar medium turbulence, and protostellar disk evolution. Notable projects that have used or been inspired by AREPO methods include the Illustris project, the IllustrisTNG suite, zoom-in campaigns related to the Aquarius Project, and comparisons with outputs from the Millennium simulation and Bolshoi simulation. Teams at the Harvard Center for Astrophysics, Princeton University, University of California, Berkeley, Cambridge University, University of Oxford, and the Kavli Institute for Cosmology have employed AREPO to study feedback from active galactic nucleuss, star formation histories compared against the Sloan Digital Sky Survey, chemical enrichment linked to yields from Type Ia supernovae and Type II supernovae, and baryon cycling in agreement or tension with observations from Hubble Space Telescope, Chandra X-ray Observatory, ALMA, and the Very Large Telescope.

Performance and Validation

Performance studies of AREPO compare favorably with both moving-particle codes like GADGET-3 and fixed-grid codes such as ENZO and RAMSES for problems involving strong shocks, shear flows, and multiphase media. Validation tests include standard suites: the Sod shock tube, the Kelvin–Helmholtz instability, the Rayleigh–Taylor instability, Sedov blast wave tests similar to those used in FLASH (software) comparisons, and cosmological structure formation benchmarks traceable to outputs from the Millennium simulation and the Coyote Universe. Parallel scaling analyses have been reported on high-performance computing systems at centers including Max Planck Computing and Data Facility, NERSC, PRACE facilities, and leadership-class machines at Argonne National Laboratory and Oak Ridge National Laboratory. Comparisons against observational datasets from Planck (spacecraft), WMAP, and redshift surveys like 2dF Galaxy Redshift Survey and SDSS have been used to calibrate subgrid prescriptions.

Licensing and Availability

AREPO's distribution has been controlled by its original authors and host institutions, with access typically granted to academic collaborators and through individual agreements, similar in practice to the managed-release models adopted by codes developed at the Max Planck Institute for Astrophysics and research groups at Heidelberg University. Source access and licensing terms have been discussed at conferences organized by societies such as the American Astronomical Society, International Astronomical Union, and workshops at facilities like the Swinburne University Centre for Astrophysics and Supercomputing. Users often coordinate with groups at Heidelberg Institute for Theoretical Studies and the Max Planck Society for collaboration and support.

Criticism and Limitations

Critiques of AREPO include concerns about the proprietary access model compared with fully open-source projects like ENZO, FLASH (software), GIZMO, and PLUTO (code), debates over numerical diffusion and angular momentum conservation relative to smoothed particle hydrodynamics implementations exemplified by GADGET-2 and grid approaches like RAMSES, and the complexity of coupling detailed radiative transfer and chemistry modules as done in community codes such as CLOUDY and MOCASSIN. Limitations cited in literature involve mesh regularization choices, overheads for mesh construction in extreme adaptive scenarios, and challenges reproducing certain observational scaling relations without calibrated subgrid models—issues addressed in follow-up work by groups at Harvard University, Princeton University, University of Cambridge, and the Max Planck Institute for Astrophysics.

Category:Astrophysical simulation software