English
 
Help Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

New pseudospectral code for the construction of initial data

MPS-Authors
/persons/resource/persons192099

Dietrich,  Tim
Multi-messenger Astrophysics of Compact Binaries, AEI-Golm, MPI for Gravitational Physics, Max Planck Society;
Astrophysical and Cosmological Relativity, AEI-Golm, MPI for Gravitational Physics, Max Planck Society;

External Resource
No external resources are shared
Fulltext (restricted access)
There are currently no full texts shared for your IP range.
Fulltext (public)

2109.14511.pdf
(Preprint), 2MB

Supplementary Material (public)
There is no public supplementary material available
Citation

Rashti, A., Fabbri, F. M., Brügmann, B., Chaurasia, S. V., Dietrich, T., Ujevic, M., et al. (2022). New pseudospectral code for the construction of initial data. Physical Review D, 105: 104027. doi:10.1103/PhysRevD.105.104027.


Cite as: https://hdl.handle.net/21.11116/0000-000A-97CB-7
Abstract
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 pseudo-spectral 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 $\sim 0.8$.