Reconstruction of the neutron star equation of state from $w$-quasinormal modes spectra with a piecewise polytropic meshing and refinement method

TL;DR

This paper introduces a novel piecewise polytropic meshing and refinement method to reconstruct the neutron star equation of state from quasinormal mode spectra, enabling accurate predictions of physical parameters and constraints from gravitational wave data.

Contribution

The paper presents a new inverse method combining quasinormal modes and piecewise polytropic parametrization for neutron star EoS reconstruction, improving accuracy and efficiency.

Findings

Reconstructed EoS closely matches original with less than 2.5% difference in physical parameters.

Method effectively predicts tidal deformability and rotation parameters.

Constraints from GW170817 are incorporated into the reconstruction process.

Abstract

In this paper we present a new approach to the inverse problem for relativistic stars using quasinormal modes and the piecewise polytropic parametrization of the equation of state. The algorithm is a piecewise polytropic meshing and refinement method that reconstructs the neutron star equation of state from experimental data of the mass and the $wI$-quasinormal modes. We present an algorithm able to numerically calculate axial quasinormal modes of neutron stars in an efficient way. We use an initial mesh of $27440$ equations of state in a $4$-volume of piecewise polytropic parameters that contains most of the candidate equations of state used today. The refinement process drives us to the reconstruction of the equation of state with a certain precision. Using the reconstructed equation of state, we calculate predictions for tidal deformability and slow rotation parameters (moment of inertia and quadrupole moment, for example). In order to check the method with an explicit example, we use as input data a few (five) configurations of a given equation of state. We reconstruct the equation of state in a quite good approximation, and then we compare the curves of physical parameters from the original equation of state and the reconstructed one. We obtain a relative difference for all of the parameters smaller than $2.5\%$ except for the tidal deformability, for which we obtain a relative difference smaller than $6.5\%$. We also study constraints from GW170817 event for the reconstructed equation of state.