The general solution of hydrostatic equilibrium equations for a two-component fluid of ions and electrons without a local electroneutrality constraint is found in the framework of Newtonian gravity theory. In agreement with the Poincar\'e theorem on analyticity and in the context of Dyson's argument, the general solution is demonstrated to possess a fixed (essential) singularity in the gravitational constant $G$ at $ G = 0 $. The regular component of the general solution can be determined by perturbation theory in $G$ starting from a locally neutral solution. The non-perturbative component obtained using the method of Wentzel, Kramers and Brillouin is exponentially small in the inner layers of the star and grows rapidly in the outward direction. Near the surface of the star, both components are comparable in magnitude, and their non-linear interplay determines the properties of an electro- or ionosphere. The stellar charge varies within the limits of $- 0.1 $ to $150$ C per solar mass. The properties of electro- and ionospheres are exponentially sensitive to variations of the fluid densities in the central regions of the star. The general solutions of two exactly solvable stellar models without a local electroneutrality constraint are also presented.

If strange matter is absolutely stable, the ordinary nuclei decay to strangelets, while neutron stars convert into strange stars. Lifetimes of the ordinary nuclei are constrained experimentally to be above $\sim 10^{33}$ years, while lifetimes of the metastable neutron stars depend on the neutron star masses and can exceed the age of the Universe. As a consequence, the neutron stars and the strange stars can coexist in the Universe. We point out that numerical simulations of the conversion of neutron stars to strange stars, performed by M. Herzog and F. K. Roepke in Phys. Rev. D 84, 083002 (2011) [arXiv:1109.0539], are focused on a region in the parameter space of strange matter, in which low-mass neutron stars and strange stars are coexistent, whereas massive neutron stars are unstable and short lived on an astronomical timescale.

Phase transition from hadronic matter to quark-gluon matter is discussed for various regimes of temperature and baryon number density. For small and medium densities, the phase transition is accurately described in the framework of the Field Correlation Method, whereas at high density predictions are less certain and leave room for the phenomenological models. We study formation of multiquark states (MQS) at zero temperature and high density. Relevant MQS components of the nuclear matter can be described using a previously developed formalism of the quark compound bags (QCB). Partial-wave analysis of nucleon-nucleon scattering indicates the existence of 6QS which manifest themselves as poles of $P$-matrix. In the framework of the QCB model, we formulate a self-consistent system of coupled equations for the nucleon and 6QS propagators in nuclear matter and the G-matrix. The approach provides a link between high-density nuclear matter with the MQS components and the cumulative effect observed in reactions on the nuclei, which requires the admixture of MQS in the wave functions of nuclei kinematically. 6QS determine the natural scale of the density for a possible phase transition into the MQS phase of nuclear matter. Such a phase transition can lead to dynamic instability of newly born protoneutron stars and dramatically affect the dynamics of supernovae. Numerical simulations show that the phase transition may be a good remedy for the triggering supernova explosions in the spherically symmetric supernova models. A specific signature of the phase transition is an additional neutrino peak in the neutrino light curve. For a Galactic core-collapse supernova, such a peak could be resolved by the present neutrino detectors. The possibility of extracting the parameters of the phase of transition from observation of the neutrino signal is discussed also.

Using elementary geometric tools, we apply essentially the same methods to derive expressions for the rotation angle of the swing plane of Foucault's pendulum and the rotation angle of the spin of a relativistic particle moving in a circular orbit (Thomas precession effect).

Theories beyond the standard model include a number of new particles some of which might be light and weakly coupled to ordinary matter. Such particles affect equation of state of nuclear matter and can shift admissible masses of neutron stars to higher values. The internal structure of neutron stars is modified provided the ratio between coupling strength and mass squared of a weakly interacting light boson is above $g^2/\mu^2 \sim 25 ~\mathrm{GeV}^{-2}$. We provide limits on the couplings with the strange sector, which cannot be achieved from laboratory experiments analysis. When the couplings to the first family of quarks is considered the limits imposed by the neutron stars are not more stringent than the existing laboratory ones. The observations on neutron stars give evidence that equation of state of the $\beta$-equilibrated nuclear matter is stiffer than expected from many-body theory of nuclei and nuclear matter. A weakly interacting light vector boson coupled predominantly to the second family of the quarks can produce the required stiffening.