For finding the mass-radius curve of neutron stars, we can solve TOV Equations which are constraint equations got by solving Einsteins equations. The boundary conditions are $m(r=0)=0$ and $\rho(r=0)=\rho_c$. Then we put a physical condition that for $r\rightarrow r_*$(radius of star), $\rho\rightarrow 0$, and by this we can get the limiting radius $r_*$ and mass $m_*$ and plots like
and for density
I understood and took it from MIT OCW
One can also do the same process for Modified Gravities and more involved $P(\rho)$ matter equations
- I have seen many places where the mass-radius curve looks like
This curve is from the paper "Neutron stars in Einstein-Λ gravity: the cosmological constant effects G. H. Bordbar1, , S. H. Hendi1, and B. Eslam Panah"
- Confusion: Why isn't the mass of neutron star cumulative as radius increases, i.e. it must be $0$ at $r=0$ and reach its maximum at $r=r_*$? Maybe I am not getting the right interpretation of these curves. So, I would be grateful for some explanation or references.