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

arXiv commentPRD version

References in corpus (57)

• Observation of Gravitational Waves from a Binary Black Hole Merger
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs
• GW170817: Measurements of Neutron Star Radii and Equation of State
• GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run
• Properties of the binary neutron star merger GW170817
• 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
• Dark Energy after GW170817 and GRB170817A
• Tests of General Relativity with the Binary Black Hole Signals from the LIGO-Virgo Catalog GWTC-1
• Tests of General Relativity with Binary Black Holes from the second LIGO-Virgo Gravitational-Wave Transient Catalog
• Observation of gravitational waves from two neutron star-black hole coalescences
• Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817
• Tests of General Relativity with GW170817
• A simple construction of initial data for multiple black holes
• A single-domain spectral method for black hole puncture data
• Conformal ``thin sandwich'' data for the initial-value problem of general relativity
• A multidomain spectral method for solving elliptic equations
• Extrinsic Curvature and the Einstein Constraints
• Spectral Methods for Numerical Relativity
• Spin Flips and Precession in Black-Hole-Binary Mergers
• Excision boundary conditions for black hole initial data
• Gravitational waves from spinning black hole-neutron star binaries: dependence on black hole spins and on neutron star equations of state
• Numerical relativity simulations of binary neutron stars
• 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
• The spin parameter of uniformly rotating compact stars
• Gravitational waveforms from binary neutron star mergers with high-order WENO schemes in numerical relativity
• A new public code for initial data of unequal-mass, spinning compact-object binaries
• Accurate and realistic initial data for black hole-neutron star binaries
• Excision without excision: the relativistic turducken
• Constructing quasi-equilibrium initial data for binary neutron stars with arbitrary spins
• A new numerical method to construct binary neutron star initial data
• Relativistic black hole-neutron star binaries in quasiequilibrium: effects of the black hole excision boundary condition
• Quasiequilibrium black hole-neutron star binaries in general relativity
• Puncture Evolution of Schwarzschild Black Holes
• Binary Neutron Stars with Arbitrary Spins in Numerical Relativity
• A double-domain spectral method for black hole excision data
• 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
• A new code for equilibriums and quasiequilibrium initial data of compact objects
• Relativistic hydrodynamics in the presence of puncture black holes
• The initial value problem as it relates to numerical relativity
• Puncture Initial Data for Black-Hole Binaries with High Spins and High Boosts
• Quasiequilibrium states of black hole-neutron star binaries in the moving-puncture framework
• Constraints as evolutionary systems
• Hyperbolic Relaxation Method for Elliptic Equations
• Alternatives to standard puncture initial data for binary black hole evolution
• 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
• Unified discontinuous Galerkin scheme for a large class of elliptic equations
• Black-Hole Neutron Star Simulations with the BAM code: First Tests and Simulations
• Numerical investigations of the asymptotics of solutions to the evolutionary form of the constraints
• Toward computing gravitational initial data without elliptic solvers
• Black Hole - Neutron Star Binary Mergers: The Imprint of Tidal Deformations and Debris
• Reducing orbital eccentricity in initial data of black hole--neutron star binaries in the puncture framework

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