Using the spin and flux evolution of the accreting pulsar 4U 1626$-$67 across the 2008 torque reversal, we determine the fastness parameter dependence of the dimensionless torque acting on the pulsar. We find that the dimensionless torque is qualitatively different from the existing models: it is concave-up across the torque equilibrium whereas the existing torque models predict a concave-down (convex) relation with the fastness parameter. We show that the dimensionless torque has a cubic dependence on the fastness parameter near the torque equilibrium. We also find that the torque can not attain large values away from the equilibrium, either in the positive or the negative side, but saturates at limited values. The spin-down torque can attain a 2.5 times larger magnitude at the saturation limit than the spin-up torque. From the evolution of the frequency of quasi-periodic oscillations of 4U 1626$-$67 across the torque reversal of 1990, we determine the critical fastness parameter corresponding to torque equilibrium to be $\omega_{\rm c} \simeq 0.75$ within the framework of the beat frequency model and boundary region model for reasonable values of the model parameters. We find that the disc magnetosphere interaction becomes unstable when the inner radius approaches the corotation radius as predicted by some models, though with a longer timescale. We also find that there is an unstable regime that is triggered when the fastness parameter is 0.8 times the critical fastness parameter ($\omega=0.6$ for $\omega_{\rm c} \simeq 0.75$) possibly associated with an instability observed in numerical simulations.

The glitch activity of young pulsars arises from the exchange of angular momentum between the crust and the interior of the star. Recently, it was inferred that the moment of inertia of the crust of a neutron star is not sufficient to explain the observed glitches. Such estimates are presumed in the Einstein's general relativity in describing the hydrostatic equilibrium of neutron stars. The crust of the neutron star has a space-time curvature of 14 orders of magnitude larger than that probed in solar system tests. This makes gravity the weakest constrained physics input in the crust related processes. We calculate the ratio of the crustal to the total moment of inertia of neutron stars in the scalar-tensor theory of gravity and the non-perturbative $f({\cal R})={\cal R}+ a {\cal R}^2$ gravity. We find for the former that the crust to core ratio of the moment of inertia does not change significantly from what is inferred in general relativity. For the latter we find that the ratio increases significantly from what is inferred in general relativity in the case of high mass objects. Our results suggest that the glitch activity of pulsars may be used to probe gravity models, although the gravity models explored in this work are not appropriate candidates.

Neutron stars (NSs) are thought to be excellent laboratories for determining the equation of state (EoS) of cold dense matter. Their strong gravity suggests that they can also be used to constrain gravity models. The mass and radius (M-R) of a NS both depend on the choice of EoS and gravity, meaning that NSs cannot be simultaneously good laboratories for both of these questions. A measurement of M-R would constrain the less well known physics input. The assumption that M-R measurements can be used to constrain EoS-presumes general relativity (GR) is the ultimate model of gravity in the classical regime. We calculate the radial profile of compactness and curvature (square root of the full contraction of the Weyl tensor) within a NS and determine the domain not probed by the Solar System tests of GR. We find that, except for a tiny sphere of radius less than a millimeter at the center, the curvature is several orders of magnitude above the values present in Solar System tests. The compactness is beyond the solar surface value for r>10 m, and increases by 5 orders of magnitude towards the surface. With the density being only an order of magnitude higher than that probed by nuclear scattering experiments, our results suggest that the employment of GR as the theory of gravity describing the hydrostatic equilibrium of NSs is a rather remarkable extrapolation from the regime of tested validity, as opposed to that of EoS models. Our larger ignorance of gravity within NSs suggests that a measurement of M-R constrains gravity rather than EoS, and given that EoS has yet to be determined by nucleon scattering experiments, M-R measurements cannot tightly constrain the gravity models either. Near the surface the curvature and compactness attain their largest values, while EoS in this region is fairly well known. This renders the crust as the best site to look for deviations from GR.