Elliptica: a new pseudo-spectral code for the construction of initial data

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$.

PRD version

References in corpus (21)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Dark Energy after GW170817: dead ends and the road ahead
• Constraints on a phenomenologically parameterized neutron-star equation of state
• Strong constraints on cosmological gravity from GW170817 and GRB 170817A
• Observation of gravitational waves from two neutron star-black hole coalescences
• Spin Flips and Precession in Black-Hole-Binary Mergers
• Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement
• Initial data for black hole-neutron star binaries: a flexible, high-accuracy spectral method
• A new public code for initial data of unequal-mass, spinning compact-object binaries
• Excision without excision: the relativistic turducken
• Quasiequilibrium black hole-neutron star binaries in general relativity
• Relativistic black hole-neutron star binaries in quasiequilibrium: effects of the black hole excision boundary condition
• Puncture Evolution of Schwarzschild Black Holes
• Constructing Binary Neutron Star Initial Data with High Spins, High Compactness, and High Mass-Ratios
• Quasiequilibrium sequences of black-hole--neutron-star binaries in general relativity
• Relativistic hydrodynamics in the presence of puncture black holes
• The initial value problem as it relates to numerical relativity
• A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code
• Initial data for black hole-neutron star binaries, with rotating stars
• Black-Hole Neutron Star Simulations with the BAM code: First Tests and Simulations
• Unified discontinuous Galerkin scheme for a large class of elliptic equations

Cited by in corpus (8)

• The Science of the Einstein Telescope
• Quasi-equilibrium configurations of binary systems of dark matter admixed neutron stars
• General-Relativistic Hydrodynamics Simulation of a Neutron Star - Sub-Solar-Mass Black Hole Merger
• A high-order shock capturing discontinuous Galerkin-finite-difference hybrid method for GRMHD
• Simulating binary black hole mergers using discontinuous Galerkin methods
• Multi-domain spectral method for self-force calculations
• On the stability of covariant BSSN formulation
• The Parallel Compact Object CALculator: An Efficient General Relativistic Initial Data Solver for Compact Objects