A simple model of cooling neutron stars with superfluid cores: comparison with observations

Cooling of neutron stars (NSs) with superfluid cores is simulated taking into account neutrino emission produced by Cooper pairing of nucleons. The critical temperatures of neutron and proton superfluidities, $T_{cn}$ and $T_{cp}$, are assumed to be constant over the NS core, and treated as free parameters. They are constrained using the surface temperatures $T_s$ of isolated NSs (RX J0822-43, PSR 1055-52, 1E 1207-52, Vela, Geminga, PSR 0656+14, RX J0002+62), obtained by interpretation of observed thermal radiation either with black body spectrum or with hydrogen atmosphere models.