Neutron star equilibrium configurations within a fully relativistic theory with strong, weak, electromagnetic, and gravitational interactions

We formulate the equations of equilibrium of neutron stars taking into account strong, weak, electromagnetic, and gravitational interactions within the framework of general relativity. The nuclear interactions are described by the exchange of the σ, ω, and ρ virtual mesons. The equilibrium conditions are given by our recently developed theoretical framework based on the Einstein–Maxwell–Thomas–Fermi equations along with the constancy of the general relativistic Fermi energies of particles, the “Klein potentials”, throughout the configuration. The equations are solved numerically in the case of zero temperatures and for selected parameterizations of the nuclear models. The solutions lead to a new structure of the star: a positively charged core at supranuclear densities surrounded by an electronic distribution of thickness ∼ℏ/(mec)∼102ℏ/(mπc) of opposite charge, as well as a neutral crust at lower densities. Inside the core there is a Coulomb potential well of depth ∼mπc2/e. The constancy of the Klein potentials in the transition from the core to the crust, imposes the presence of an overcritical electric field ∼(mπ/me)2Ec, the critical field being . The electron chemical potential and the density decrease, in the boundary interface, until values and ρcrust<ρcore. For each central density, an entire family of core–crust interface boundaries and, correspondingly, an entire family of crusts with different mass and thickness, exist. The configuration with separates neutron stars with and without inner crust. We present here the novel neutron star mass–radius for the especial case ρcrust=ρdrip and compare and contrast it with the one obtained from the traditional Tolman–Oppenheimer–Volkoff treatment.