GR-Athena++: Puncture evolutions on vertex-centered oct-tree adaptive mesh refinement

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ∼1.2 × 10^4 CPUs and excellent weak scaling is shown up to ∼10^5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.

Files

Metadata

Work Title GR-Athena++: Puncture evolutions on vertex-centered oct-tree adaptive mesh refinement
Access
Open Access
Creators
  1. Boris Daszuta
  2. Francesco Zappa
  3. William Cook
  4. David Radice
  5. Sebastiano Bernuzzi
  6. Viktoriya Morozova
License In Copyright (Rights Reserved)
Work Type Article
Publisher
  1. The Astrophysical Journal Supplement Series
Publication Date November 11, 2021
Publisher Identifier (DOI)
  1. https://doi.org/10.3847/1538-4365/ac157b
Deposited March 14, 2025

Versions

Analytics

Collections

This resource is currently not in any collection.

Work History

Version 1
published

  • Created

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added 2101.08289-1.pdf

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added Creator Boris Daszuta

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added Creator Francesco Zappa

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added Creator William Cook

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added Creator David Radice

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added Creator Sebastiano Bernuzzi

    March 14, 2025 20:18 by Researcher Metadata Database

  • Added Creator Viktoriya Morozova

    March 14, 2025 20:18 by Researcher Metadata Database

  • Published

    March 14, 2025 20:18 by Researcher Metadata Database

  • Updated

    March 14, 2025 22:04 by [unknown user]

  • Updated Publisher, Description Show Changes

    April 15, 2025 20:57 by avs5190

    Publisher
    • Astrophysical Journal, Supplement Series
    • The Astrophysical Journal Supplement Series
    Description
    • Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ∼1.2 × 10<sup>4</sup> CPUs and excellent weak scaling is shown up to ∼10<sup>5</sup> CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.
    • Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ∼1.2 × 10^4 CPUs and excellent weak scaling is shown up to ∼10^5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.