Deformed neutron stars

J. T. Quartuccio,P. H. R. S. Moraes,J. D. V. Arbañil

2024-02-29

Abstract:We present solutions for non-spherically symmetric neutron stars. We begin by deriving the Tolman-Oppenheimer-Volkoff equations from a parameterized metric that takes into account the deformation of the star due to differences in equatorial and polar pressures, expressed in terms of a parameter D, which is the ratio between polar and equatorial radius. The stellar structure is solved using the GM1 equation of state and the Tolman-Oppenheimer-Volkoff equations for deformed objects are numerically integrated using the fourth-order Runge-Kutta method for different values of the parameter D. We show that larger values of D > 1, that describe prolate neutron stars, yield smaller values of mass and radius, while for smaller values of D < 1, describing oblate neutron stars, larger values for mass and radius are attained. From the confrontation of our model theoretical predictions with recent observational data on pulsars, it is possible to constrain the values of the parameter D.

General Relativity and Quantum Cosmology

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to study the structure and properties of non - spherically symmetric neutron stars (deformed neutron stars, DNSs). Specifically, by introducing a deformation parameter \(D\) that describes the pressure difference between the polar and equatorial regions, the authors derived the general relativity equations applicable to such non - spherically symmetric celestial bodies - the so - called “\(D\)-TOV equations”. These equations are used for numerical simulations of neutron stars with different \(D\) values and to explore how their mass, radius and other physical properties change with \(D\).

### Main problems
1. **Influence of non - sphericity**: Traditional models usually assume that neutron stars are spherically symmetric celestial bodies, but this may not be the case in reality. This paper aims to consider the shape deformation of neutron stars due to internal pressure differences, especially by introducing the deformation parameter \(D\) to quantify this non - sphericity.

2. **Comparison between theory and observational data**: By performing numerical simulations of neutron stars with different \(D\) values, the authors hope to compare theoretical predictions with existing pulsar observational data, thereby constraining the reasonable value range of the deformation parameter \(D\).

3. **Exploration of the equation of state under extreme conditions**: To study the behavior of neutron stars under extremely high - density conditions (such as the secondary compact object in the GW190814 event), which helps to understand the equation of state (EoS) of matter under extreme conditions and verify whether there is a more “stiff” EoS that can explain the observed mass scale.

### Key formulas
- **Metric**:
\[
ds^2=-e^{2\phi(r)}dt^2+(1 - \frac{2m(r)}{r})^{-D}dr^2+r^2(d\theta^2+\sin^2\theta d\phi^2)
\]
where \(D\) is the deformation parameter, representing the ratio of the polar radius to the equatorial radius.

- **Energy - momentum tensor**:
\[
T_{\mu\nu}=\begin{pmatrix}
\rho&0&0&0\\
0&p_r&0&0\\
0&0&p_z&0\\
0&0&0&p_{rz}
\end{pmatrix}
\]

- **\(D\)-TOV equations**:
\[
\frac{dp_r}{dr}=-(\rho + p_r)\left[\frac{8\pi r^2p_r-(1 - \frac{2m}{r})^D + 1}{2r(1 - \frac{2m}{r})^D}\right]
\]
\[
\frac{dp_z}{dz}=-D(\rho + p_z)\left[\frac{8\pi(\frac{z}{D})^2p_z-(1 - \frac{2mD}{z})^D + 1}{2z(1 - \frac{2mD}{z})^D}\right]
\]

Through these formulas, the authors can systematically study the internal structure and physical properties of non - spherically symmetric neutron stars, thus providing a new perspective for understanding compact celestial bodies.