The self-consistent general relativistic solution for a system of degenerate neutrons, protons and electrons in β-equilibrium

We present the self-consistent treatment of the simplest, nontrivial, self-gravitating system of degenerate neutrons, protons and electrons in β-equilibrium within relativistic quantum statistics and the Einstein-Maxwell equations. The impossibility of imposing the condition of local charge neutrality on such systems is proved, consequently overcoming the traditional Tolman-Oppenheimer-Volkoff treatment. We emphasize the crucial role of imposing the constancy of the generalized Fermi energies. A new approach based on the coupled system of the general relativistic Thomas-Fermi-Einstein-Maxwell equations is presented and solved. We obtain an explicit solution fulfilling global and not local charge neutrality by solving a sophisticated eigenvalue problem of the general relativistic Thomas-Fermi equation. The value of the Coulomb potential at the center of the configuration is eV(0)≃mπc2 and the system is intrinsically stable against Coulomb repulsion in the proton component. This approach is necessary, but not sufficient, when strong interactions are introduced.