We have recently developed a neutron star model fulfilling global and not local charge neutrality, both in the static and in the uniformly rotating cases. The model is described by the coupled Einstein-Maxwell-Thomas- Fermi (EMTF) equations, in which all fundamental interactions are accounted for in the framework of general relativity and relativistic mean field theory. Uniform rotation is introduced following the Hartle's formalism. We show that the use of realistic parameters of rotating neutron stars obtained from numerical integration of the self-consistent axisymmetric general relativistic equations of equilibrium leads to values of the magnetic field and radiation efficiency of pulsars very different from estimates based on fiducial parameters assuming a neutron star mass, M = 1.4 Msun, radius R = 10 km, and moment of inertia, I = 10^45 g cm^2. In addition, we compare and contrast the magnetic field inferred from the traditional Newtonian rotating magnetic dipole model with respect to the one obtained from its general relativistic analog which takes into due account the effect of the finite size of the source. We apply these considerations to the specific high-magnetic field pulsars class and show that, indeed, all these sources can be described as canonical pulsars driven by the rotational energy of the neutron star, and with magnetic fields lower than the quantum critical field for any value of the neutron star mass.

In our previous treatment of neutron stars, we have developed the model fulfilling global and not local charge neutrality. In order to implement such a model, we have shown the essential role by the Thomas-Fermi equations, duly generalized to the case of electromagnetic field equations in a general relativistic framework, forming a coupled system of equations that we have denominated Einstein-Maxwell-Thomas-Fermi (EMTF) equations. From the microphysical point of view, the weak interactions are accounted for by requesting the \beta\ stability of the system, and the strong interactions by using the \sigma-\omega-\rho\ nuclear model, where \sigma, \omega\ and \rho\ are the mediator massive vector mesons. Here we examine the equilibrium configurations of slowly rotating neutron stars by using the Hartle formalism in the case of the EMTF equations indicated above. We integrate these equations of equilibrium for different central densities and circular angular velocities and compute the mass, polar and equatorial radii, angular momentum, eccentricity, moment of inertia, as well as quadrupole moment of the configurations. Both the Keplerian mass-shedding limit and the axisymmetric secular instability are used to construct the new mass-radius relation. We compute the maximum and minimum masses and rotation frequencies of neutron stars. We compare and contrast all the results for the global and local charge neutrality cases.

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 sigma, omega, and rho 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 parametrization 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 $\sim \hbar/(m_e c)$ of opposite charge, as well as a neutral crust at lower densities. Inside the core there is a Coulomb potential well of depth $\sim m_\pi c^2/e$. The constancy of the Klein potentials in the transition from the core to the crust, impose the presence of an overcritical electric field $\sim (m_\pi/m_e)^2 E_c$, the critical field being $E_c=m^2_e c^3/(e \hbar)$. The electron chemical potential and the density decrease, in the boundary interface, until values $\mu^{\rm crust}_e < \mu^{\rm core}_e$ and $\rho_{\rm crust}<\rho_{\rm 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 $\rho_{\rm crust}=\rho_{\rm drip}\sim 4.3\times 10^{11}$ g/cm$^3$ separates neutron stars with and without inner crust. We present here the novel neutron star mass-radius for the case $\rho_{\rm crust}=\rho_{\rm drip}$ and compare and contrast it with the one obtained from the Tolman-Oppenheimer-Volkoff treatment.