[Paper Review] Spectral methods in general relativity -- toward the simulation of 3D-gravitational collapse of neutron stars

S. Bonazzola, J. Frieben+2 more|ArXiv.org|Apr 12, 1996

Pulsars and Gravitational Waves Research30 references11 citations

TL;DR

This paper presents a spectral method in a general spherical-like coordinate system for solving Einstein's equations in general relativity, enabling highly accurate simulations of neutron star collapse. The approach achieves orders of magnitude higher accuracy than finite difference methods, with successful applications to spherically symmetric collapse and magnetized, rotating neutron stars, and lays the groundwork for 3D non-axisymmetric simulations using a tensor-based scheme.

ABSTRACT

Several applications of spectral methods to problems related to the relativistic astrophysics of compact objects are presented. Based on a proper definition of the analytical properties of regular tensorial functions we have developed a spectral method in a general sphericallike coordinate system. The applications include the investigation of spherically symmetric neutron star collapse as well as the solution of the coupled 2D-Einstein-Maxwell equations for magnetized, rapidly rotating neutron stars. In both cases the resulting codes are efficient and give results typically several orders of magnitude more accurate than equivalent codes based on finite difference schemes. We further report the current status of a 3D-code aiming at the simulation of non-axisymmetric neutron star collapse where we have chosen a tensor based numerical scheme.

Motivation & Objective

• To develop a high-accuracy spectral method for solving the Einstein equations in general relativity applicable to compact relativistic objects.
• To address the numerical challenges in simulating 3D gravitational collapse of neutron stars with improved precision over finite difference methods.
• To extend spectral techniques to non-axisymmetric, 3D configurations of neutron stars, including magnetic and rotational effects.
• To establish a robust numerical framework using tensorial spectral methods for solving the coupled Einstein-Maxwell equations in 2D and 3D.
• To achieve computational efficiency and accuracy suitable for simulating complex astrophysical phenomena such as neutron star collapse and gravitational wave emission.

Proposed method

• The method employs spectral collocation techniques in a general spherically-like coordinate system to discretize the Einstein field equations.
• It is based on a rigorous mathematical treatment of regular tensorial functions in curved spacetime, ensuring stability and convergence.
• The spectral approach uses global basis functions (e.g., Chebyshev or Legendre polynomials) to represent spatial derivatives with high-order accuracy.
• A tensor-based numerical scheme is implemented to handle the full 3D system of equations for non-axisymmetric collapse.
• The method solves the coupled 2D-Einstein-Maxwell equations for magnetized, rapidly rotating neutron stars using spectral decomposition in angular and radial coordinates.
• The code is validated on spherically symmetric cases and extended to 3D using a pseudo-spectral approach with adaptive grid resolution.

Experimental results

Research questions

• RQ1 Can spectral methods achieve significantly higher accuracy than finite difference schemes in simulating relativistic neutron star collapse?
• RQ2 How can spectral methods be generalized to handle tensorial fields in non-trivial, spherically-like coordinate systems relevant to neutron stars?
• RQ3 What is the feasibility and performance of extending spectral methods to 3D non-axisymmetric gravitational collapse of neutron stars?
• RQ4 How accurately can spectral methods solve the coupled Einstein-Maxwell equations for magnetized, rotating neutron stars?
• RQ5 Can spectral methods maintain stability and convergence in complex, nonlinear systems like 3D neutron star collapse with strong gravity and magnetic fields?

Key findings

• The spectral method achieves results several orders of magnitude more accurate than equivalent finite difference codes for spherically symmetric neutron star collapse.
• The method successfully solves the 2D-Einstein-Maxwell system for magnetized, rapidly rotating neutron stars with high precision and stability.
• The 3D code framework is under development and shows promise for simulating non-axisymmetric collapse using a tensor-based spectral approach.
• The use of spectral methods in a general spherically-like coordinate system enables accurate representation of tensorial fields in curved spacetime.
• The method demonstrates superior convergence rates and reduced numerical dissipation compared to finite difference schemes.
• The approach is computationally efficient for problems with smooth solutions, such as equilibrium and quasi-equilibrium neutron star configurations.