New pseudospectral code for the construction of initial data
- Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudospectral code ELLIPTICA, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, ELLIPTICA deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of ELLIPTICA and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolveNumerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudospectral code ELLIPTICA, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, ELLIPTICA deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of ELLIPTICA and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude ∼0.8.…
Author details: | Alireza RashtiORCiD, Francesco Maria FabbriORCiD, Bernd BrügmannORCiD, Swami Vivekanandji ChaurasiaORCiD, Tim DietrichORCiDGND, Maximiliano Ujevic, Wolfgang TichyORCiD |
---|---|
DOI: | https://doi.org/10.1103/PhysRevD.105.104027 |
ISSN: | 2470-0010 |
ISSN: | 2470-0029 |
Title of parent work (English): | Physical review D |
Publisher: | American Physical Society |
Place of publishing: | College Park |
Publication type: | Article |
Language: | English |
Date of first publication: | 2022/05/16 |
Publication year: | 2022 |
Release date: | 2024/06/13 |
Volume: | 105 |
Issue: | 10 |
Article number: | 104027 |
Number of pages: | 18 |
Funding institution: | NSF [PHY-1707227, PHY-2011729]; DFG [BR 2176/5-1]; Coordenacao de; Aperfeicoamento de Pessoal de Nivel Superior-Brasil (CAPES); [88887.571346/2020-00]; Swedish Research council (VR) [2016-06012];; North German Supercomputing Alliance (HLRN) [bbp00049]; HighPerformance; Computing Center Stuttgart (HLRS) [GWanalysis 44189]; Leibniz; Supercomputing Centre (LRZ) [pn29ba] |
Organizational units: | Mathematisch-Naturwissenschaftliche Fakultät / Institut für Physik und Astronomie |
DDC classification: | 5 Naturwissenschaften und Mathematik / 52 Astronomie / 520 Astronomie und zugeordnete Wissenschaften |
5 Naturwissenschaften und Mathematik / 53 Physik / 530 Physik | |
Peer review: | Referiert |
License (German): | Keine öffentliche Lizenz: Unter Urheberrechtsschutz |