Published December 19, 2024 | Version Published
Journal Article Open

Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method

Abstract

We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods.

Copyright and License

© 2024 The Author(s). Published by IOP Publishing Ltd.

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 license. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

Acknowledgement

Charm++/Converse [144] was developed by the Parallel Programming Laboratory in the Department of Computer Science at the University of Illinois at Urbana-Champaign. The figures in this article were produced with matplotlib [145146], TikZ [147] and ParaView [148149]. Computations were performed with the Wheeler cluster at Caltech and the mbot cluster at Cornell. This work was supported in part by the Sherman Fairchild Foundation and by NSF Grants Nos. PHY-2309211, PHY-2309231, and OAC-2209656 at Caltech, and NSF Grants Nos. PHY-2207342, PHY-2407742 and OAC-2209655 at Cornell. FF gratefully acknowledges support from the Department of Energy, Office of Science, Office of Nuclear Physics, under Contract Number DE-AC02-05CH11231, from NASA through Grant 80NSSC22K0719, and from the NSF through Grant AST-2107932. MD gratefully acknowledges support from the NSF through Grant PHY-2110287 and support from NASA through Grant 80NSSC22K0719. GL and MSB acknowledge support from NSF Award PHY-2208014, the Dan Black Family Trust and Nicholas and Lee Begovich. IL acknowledges support from the Department of Energy under Award Number DE-SC0023101. ERM acknowledges support by the National Science Foundation under Grant Nos. AST-2307394 and PHY-2309210, the NSF Frontera supercomputer under Grant AST21006, and Delta at the National Center for Supercomputing Applications (NCSA) through allocation PHY210074 from the Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support (ACCESS) program, which is supported by National Science Foundation Grants #2138259, #2138286, #2138307, #2137603, and #2138296. ERM further acknowledges support on Perlmutter through NERSC under Grant m4575. PK acknowledges support of the Department of Atomic Energy, Government of India, under Project No. RTI4001, and by the Ashok and Gita Vaish Early Career Faculty Fellowship at the International Centre for Theoretical Sciences.

Funding

This work was supported in part by the Sherman Fairchild Foundation and by NSF Grants Nos. PHY-2309211, PHY-2309231, and OAC-2209656 at Caltech, and NSF Grants Nos. PHY-2207342, PHY-2407742 and OAC-2209655 at Cornell. FF gratefully acknowledges support from the Department of Energy, Office of Science, Office of Nuclear Physics, under Contract Number DE-AC02-05CH11231, from NASA through Grant 80NSSC22K0719, and from the NSF through Grant AST-2107932. MD gratefully acknowledges support from the NSF through Grant PHY-2110287 and support from NASA through Grant 80NSSC22K0719. GL and MSB acknowledge support from NSF Award PHY-2208014, the Dan Black Family Trust and Nicholas and Lee Begovich. IL acknowledges support from the Department of Energy under Award Number DE-SC0023101. ERM acknowledges support by the National Science Foundation under Grant Nos. AST-2307394 and PHY-2309210, the NSF Frontera supercomputer under Grant AST21006, and Delta at the National Center for Supercomputing Applications (NCSA) through allocation PHY210074 from the Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support (ACCESS) program, which is supported by National Science Foundation Grants #2138259, #2138286, #2138307, #2137603, and #2138296. ERM further acknowledges support on Perlmutter through NERSC under Grant m4575. PK acknowledges support of the Department of Atomic Energy, Government of India, under Project No. RTI4001, and by the Ashok and Gita Vaish Early Career Faculty Fellowship at the International Centre for Theoretical Sciences.

Data Availability

The data that support the findings of this study are openly available at the following URL/DOI: github.com/sxs-collaboration/spectre.

Files

Deppe_2024_Class._Quantum_Grav._41_245002.pdf

Files (2.2 MB)

Name Size Download all
md5:de06393445ccd7fc4cbce92f07b9faf2
2.2 MB Preview Download

Additional details

Related works

Is new version of
Discussion Paper: arXiv:2406.19038 (arXiv)
Is supplemented by
Dataset: https://github.com/sxs-collaboration/spectre (URL)

Funding

Sherman Fairchild Foundation
National Science Foundation
PHY-2309211
National Science Foundation
PHY-2309231
National Science Foundation
OAC-2209656
National Science Foundation
PHY-2207342
National Science Foundation
PHY-2407742
National Science Foundation
OAC-2209655
United States Department of Energy
DE-AC02-05CH11231
National Aeronautics and Space Administration
80NSSC22K0719
National Science Foundation
AST-2107932
National Science Foundation
PHY-2110287
National Science Foundation
PHY-2208014
United States Department of Energy
DE-SC0023101
National Science Foundation
AST-2307394
National Science Foundation
PHY-2309210
National Science Foundation
AST21006
National Science Foundation
PHY210074
National Science Foundation
OAC-2138259
National Science Foundation
OAC-2138286
National Science Foundation
OAC-2138307
National Science Foundation
OAC-2137603
National Science Foundation
OAC-2138296
Natural Sciences and Engineering Research Council
m4575
Department of Atomic Energy
RTI4001
International Centre for Theoretical Sciences

Dates

Accepted
2024-10-18
Accepted
Available
2024-11-11
Published

Caltech Custom Metadata

Caltech groups
Astronomy Department, TAPIR, Walter Burke Institute for Theoretical Physics
Publication Status
Published