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.

Fallback disks are expected to form around young neutron stars. The presence of these disks can be revealed by their blackbody spectrum in the infrared, optical and UV bands. We present a re-reduction of the archival optical and infrared data of the Vela pulsar, together with the existing infrared and UV spectrum of Vela, and model their unpulsed components with the blackbody spectrum of a supernova debris disk. We invoke the quiescent disk solution of Sunyaev and Shakura for the description of the disk in the propeller stage and find the inner radius of the disk to be inside the light cylinder radius. We perform a high resolution X-ray analysis with XMM-Newton and find a narrow absorption feature at 0.57 keV which can be interpreted as K$_{\alpha}$ line of He-like oxygen (OVII). The strength of the line indicates an element over-abundance in our line of sight exceeding the amounts that would be expected from ISM. The spectral feature may originate from the pulsar wind nebula and may be partly caused by the reprocessed X-ray radiation by the fallback disk. We discuss the lower-than-3 braking index of Vela as partially due to the contribution of the propeller torques. Our results suggest that the pulsar mechanism can work simultaneously with the propeller processes and that the debris disks can survive the radiation pressure for at least $\sim 10^4$ years. As Vela is a relatively close object, and a prototypical pulsar, the presence of a disk, if confirmed, may indicate the ubiquity of debris disks around young neutron stars.