Source term method for binary neutron stars initial data

DOI10.1088/1361-6382/ABFC29zbMATH Open1481.85001arXiv2010.08733OpenAlexW3093470146MaRDI QIDQ5163391FDOQ5163391

Authors: Bing Jyun Tsao, Roland Haas, Antonios Tsokaros 

Publication date: 3 November 2021

Published in: Classical and Quantum Gravity (Search for Journal in Brave)

Abstract: The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: a) the boundary is not known a-priori, but constitutes part of the solution of the problem; b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver COCAL, we demonstrate agreement to approximately $1\mathrm{\%}$. Our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

Full work available at URL: https://arxiv.org/abs/2010.08733

Recommendations

• A new numerical method to construct binary neutron star initial data
• A multi-domain spectral method for initial data of arbitrary binaries in general relativity
• Initial data for high-compactness black hole-neutron star binaries
• Helically reduced wave equations and binary neutron stars
• Models of helically symmetric binary systems

zbMATH Keywords

initial databinary neutron starsPoisson solverssource term method

Mathematics Subject Classification ID

Iterative numerical methods for linear systems (65F10) Galactic and stellar dynamics (85A05) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Euler-Poisson-Darboux equations (35Q05) Boundary value problems on manifolds (58J32) Two-body problems (70F05) Propagation of singularities; initial value problems on manifolds (58J47)

Cites Work

• Level set methods and dynamic implicit surfaces
• Two methods for discretizing a delta function supported on a level set
• Finite difference methods for approximating Heaviside functions
• A source term method for Poisson problems on irregular domains
• A new numerical method to construct binary neutron star initial data

Cited In (5)

• A source term method for Poisson problems with a discontinuous diffusion coefficient
• A new numerical method to construct binary neutron star initial data
• Methods for relativistic self-gravitating fluids: from binary neutron stars to black hole-disks and magnetized rotating neutron stars
• Accurate evolutions of inspiralling neutron-star binaries: assessment of the truncation error
• Realistic binary neutron star initial data with \texttt{Elliptica}

This page was built for publication: Source term method for binary neutron stars initial data

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5163391)