Extremely High-Order Convergence in Simulations of Relativistic Stars

Abstract

We provide a road towards obtaining gravitational waveforms from inspiraling material binaries with an accuracy viable for third-generation gravitational wave detectors, without necessarily advancing computational hardware or massively-parallel software infrastructure. We demonstrate a proof-of-principle 1+1-dimensional numerical implementation that exhibits up to 7th-order convergence for highly dynamic barotropic stars in curved spacetime, and numerical errors up to 6 orders of magnitude smaller than a standard method. Aside from high-order interpolation errors (Runge's phenomenon), there are no obvious fundamental obstacles to obtaining convergence of even higher order. The implementation uses a novel surface-tracking method, where the surface is evolved and high-order accurate boundary conditions are imposed there. Computational memory does not need to be allocated to fluid variables in the vacuum region of spacetime. We anticipate the application of this new method to full $3+1$-dimensional simulations of the inspiral phase of compact binary systems with at least one material body. The additional challenge of a deformable surface must be addressed in multiple spatial dimensions, but it is also an opportunity to input more precise surface tension physics.

Pdf ~1 Mb. 13 page body, 2 pages of appendices. Video content available at https://youtu.be/rBpz-aVNSLo . v2: added references and footnote, minor style change

References in corpus (12)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Advanced Virgo: a 2nd generation interferometric gravitational wave detector
• Advanced LIGO
• Exploring the Sensitivity of Next Generation Gravitational Wave Detectors
• Dynamical Tides in General Relativity: Effective Action and Effective-One-Body Hamiltonian
• Sixth post-Newtonian local-in-time dynamics of binary systems
• Binary dynamics at the fifth and fifth-and-a-half post-Newtonian orders
• Sub-radian-accuracy gravitational waveforms of coalescing binary neutron stars in numerical relativity
• Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves
• Smooth equations of state for high-accuracy simulations of neutron star binaries
• Increasing the Accuracy of Binary Neutron Star Simulations with an improved Vacuum Treatment
• Hamilton-Jacobi hydrodynamics of pulsating relativistic stars