A numerical rotating neutron star solver is used to study the temporal evolution of accreting neutron stars using a multi-polytrope model for the nuclear equation of state named ACB5. The solver is based on a quadrupole expansion of the metric, but confirms the results of previous works, revealing the possibility of an abrupt transition of a neutron star from a purely hadronic branch to a third-family branch of stable hybrid stars, passing through an unstable intermediate branch. The accretion is described through a sequence of stationary rotating {stellar} configurations which lose angular momentum through magnetic dipole emission while, at the same time, gaining angular momentum through mass accretion. The model has several free parameters which are inferred from observations. The mass accretion scenario is studied in dependence on the effectiveness of angular momentum transfer which determines at which spin frequency the neutron star will become unstable against gravitational collapse to the corresponding hybrid star on the stable third-family branch. It is conceivable that the neutrino burst which accompanies the deconfinement transition may trigger a pulsar kick which results in the eccentric orbit. A consequence of the present model is the prediction of a correlation between the spin frequency of the millisecond pulsar in the eccentric orbit and its mass at birth.

Observational evidence implying the accelerated expansion of the universe has been the motivation to develop various classes of modified gravity theories. One of them uses the so-called "screening mechanism", which is successful in reproducing the observed gravitational behavior in large scales as well as being in agreement with tests of general relativity in the solar system. In this work, we investigate an example of scalar-tensor theories with screening mechanism, namely the profile of a Chameleon field around inhomogeneous astrophysical objects. According to [Khoury and Weltman(2004)], one can define two kinds of approaches applicable to the thin shell and thick shell regimes, that allow for a solution to the Chameleon equation of motion. For sufficiently large objects, the scalar field can be assumed to propagate from a thin shell of the object instead of the whole body, which simplifies the problem. On the other hand, this solution is not practical in small objects. We find that in inhomogeneous objects this is not trivial and at least one more factor, which turns out to be the density, can change the way of approaching this problem.

We present a systematic investigation of the cosmological constant effects in compact stars interiors in the framework of Einstein's gravity. Consideration of a cosmological constant $\Lambda$ in compact stars is motivated by the mechanism of acceleration of the observable universe, where $\Lambda$ is usually related to the dark energy. In particular, we consider compact star mass twins, hybrid neutron stars that populate both the second and third branch of the mass-radius diagram. For those models, the need of consideration of excluded volume effects in the equation of state, resulting from the finite size volume of nucleons, leads to a stiffening of matter causing compact stars to acquire higher mass and radius values. We demonstrate that certain values of the cosmological constant can also modify the compact star properties but in an opposite way. In addition, we find that the inclusion of $\Lambda$ can have a similar effect to the existence of pasta phases at the hadron-quark interface.