We are investigating universal relations between different normalisations of the moment of inertia and the compactness of neutron stars in slow rotation approximation. We study the relations in particular class of massive scalar-sensor theories with self-interaction, for which significant deviations from General Relativity are allowed for values of the parameters that are in agreement with the observations. Moment of inertia-compactness relations are examined for different normalisation of the moment of inertia. It is shown that for all studied cases the deviations from EOS universality are small for the examined equations of state. On the other hand the scalarization can lead to large deviations from the general relativistic universal relations for values of the parameters that are in agreement with the current observations that can be potentially used to set further test the scalar-tensor theories.

Binary pulsar observations and gravitational wave detections seriously constrained scalar-tensor theories with massless scalar field allowing only small deviations from general relativity. If we consider a nonzero mass of the scalar field, though, significant deviations from general relativity are allowed for values of the parameters that are in agreement with the observations. In the present paper we extend this idea and we study scalar-tensor theory with massive field with self-interaction term in the potential. The additional term suppresses the scalar field in the neutron star models in addition to the effect of the mass of the scalar field but still, large deviations from pure GR can be observed for values of the parameters that are in agreement with the observations.

In the scalar-tensor theories with a massive scalar field the coupling constants, and the coupling functions in general, which are observationally allowed, can differ significantly from those in the massless case. This fact naturally implies that the scalar-tensor neutron stars with a massive scalar field can have rather different structure and properties in comparison with their counterparts in the massless case and in general relativity. In the present paper we study slowly rotating neutron stars in scalar-tensor theories with a massive gravitational scalar. Two examples of scalar-tensor theories are examined - the first example is the massive Brans-Dicke theory and the second one is a massive scalar-tensor theory indistinguishable from general relativity in the weak field limit. In the later case we study the effect of the scalar field mass on the spontaneous scalarization of neutron stars. Our numerical results show that the inclusion of a mass term for the scalar field indeed changes the picture drastically compared to the massless case. It turns out that mass, radius and moment of inertia for neutron stars in massive scalar-tensor theories can differ drastically from the pure general relativistic solutions if sufficiently large masses of the scalar field are considered.