When a rotating neutron star loses angular momentum, the progressive reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fern\'andez & Reisenegger (2005) for non-superfluid neutron stars and by Petrovich & Reisenegger (2010) for superfluid millisecond pulsars. Both studies found that pulsars reach a quasi-steady state in which the compression driving the matter out of beta equilibrium is balanced by the reactions trying to restore the equilibrium. We extend previous studies by considering the effect of density-dependence and anisotropy of the superfluid energy gaps, for the case in which the dominant reactions are the modified Urca processes, the protons are non-superconducting, and the neutron superfluidity is parametrized by models proposed in the literature. By comparing our predictions with the surface temperature of the millisecond pulsar PSR J0437-4715 and upper limits for twenty-one classical pulsars, we find the millisecond pulsar can be only explained by the models with the effectively largest energy gaps (type B models), the classical pulsars require with the gap models that vanish for some angle (type C) and two different envelope compositions. Thus, no single model for neutron superfluidity can simultaneously account for the thermal emission of all available observations of non-accreting neutron stars, possibly due to our neglect of proton superconductivity.

When a rotating neutron star loses angular momentum, the reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fern\'andez and Reisenegger for neutron stars of non-superfluid matter and by Petrovich and Reisenegger for superfluid matter, finding that the system in both cases reaches a quasi-steady state, corresponding to a partial equilibration between compression, due to the loss of angular momentum, and reactions that try to restore the equilibrium. However, Petrovich and Reisenegger assumes a constant value of the superfluid energy gap, whereas theoretical models predict density-dependent gap amplitudes, and therefore gaps that depend on the location in the star. In this work, we try to discriminate between several proposed gap models, comparing predicted surface temperatures to the value measured for the nearest millisecond pulsar, J0437-4715.