# Enhanced cooling of neutron stars via Cooper-pairing neutrino emission

M. E. Gusakov 1 Ioffe Physical Technical Institute, Politekhnicheskaya 26, 194021 St. Petersburg, Russia, ; ; 1    A. D. Kaminker 1 Ioffe Physical Technical Institute, Politekhnicheskaya 26, 194021 St. Petersburg, Russia, ; ; 1    D.G. Yakovlev 1 Ioffe Physical Technical Institute, Politekhnicheskaya 26, 194021 St. Petersburg, Russia, ; ; 1    O. Y. Gnedin 2Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA, 2
Received xx xxxxx 2004 / Accepted xx xxxxx 2004
###### Key Words.:
Stars: neutron – dense matter
offprints: M.E. Gusakov

We simulate cooling of superfluid neutron stars with nucleon cores where direct Urca process is forbidden. We adopt density dependent critical temperatures and of singlet-state proton and triplet-state neutron pairing in a stellar core and consider a strong proton pairing (with maximum K) and a moderate neutron pairing ( K). When the internal stellar temperature falls below , the neutrino luminosity due to Cooper pairing of neutrons behaves , just as that produced by modified Urca process (in a non-superfluid star) but is higher by about two orders of magnitude. In this case the Cooper-pairing neutrino emission acts like an enhanced cooling agent. By tuning the density dependence we can explain observations of cooling isolated neutron stars in the scenario in which direct Urca process or similar process in kaon/pion condensed or quark matter are absent.

## 1 Introduction

Thanks to Chandra and XMM-Newton missions, there is a great progress in observations of thermal radiation emergent from the surfaces of isolated (cooling) middle-aged neutron stars (e.g., Pavlov & Zavlin pz03 ). A comparison of these data with theoretical models of cooling neutron stars gives a method to constrain (still poorly known) fundamental properties of supranuclear matter in neutron-star cores, such as the composition and equation of state of the matter and its superfluid properties.

So far, the observations can be explained by a number of vastly different theoretical models (e.g., Page page98a ; page98b , Tsuruta et al. tsurutaetal02 , Khodel et al. khodeletal04 , Blaschke et al. blaschke , Yakovlev & Pethick yp04 , and references therein). Particularly, one can employ the simplest models of neutron stars with the cores composed of nucleons (or nucleons/hyperons), or containing pion condensates, kaon condensates or quarks. The simplest model of a non-superfluid nucleon core which cools via modified Urca process of neutrino emission (without any powerful direct Urca process) cannot explain the observations: some neutron stars (e.g., PSR B1055–52) turn out to be much warmer, while others (e.g., the Vela pulsar) are much colder than those given by this model. Warmer stars can be explained (Kaminker et al. khy01 ) assuming a strong proton superfluidity in the core: such a superfluidity suppresses modified Urca process and slows down the cooling. However colder stars require some cooling mechanism which is faster than the modified Urca process.

Explanations of observations of colder stars presented in the literature invoke usually either a powerful direct Urca process in nucleon (or nucleon/hyperon) matter or similar processes in kaon-condensed, pion-condensed, quark matter in the inner cores of massive neutron stars.

In this paper we present a new scenario of neutron star cooling. We adopt the simplest model equation of state of supranuclear matter in neutron star cores (Douchin & Haensel dh01 ) involving only nucleons, electrons and muons. This equation of state forbids direct Urca process in all stable neutron stars. We will show that the enhanced cooling required to explain colder isolated neutron stars can be produced by neutrino emission due to a moderately strong triplet-state pairing of neutrons. This new interpretation is possible only for a specific density dependence of the critical temperature of neutron pairing.

In the next section we outline the observational basis; the cooling scenario is given afterwards.

## 2 Observations

Table 1 summarize observations of isolated (cooling) middle-aged ( yr) neutron stars, whose thermal surface radiation has been detected (or constrained). We present the estimated stellar ages and effective surface temperatures (as detected by a distant observer).

Two young objects, RX J0822–4300 and 1E 1207.4–5209, are radio-quiet neutron stars in supernova remnants; RX J1856.4–3754 and RX J0720.4–3125 are also radio-quiet neutron stars. Other objects — the Crab and the Vela pulsars, PSR B1706–44, PSR J0538+2817, Geminga, and PSR B1055–52 — are observed as radio pulsars.

RX J0205+6449 and the Crab pulsar are associated with historical supernovae and their ages are certain. For RX J0822–4300, we take the age of the host supernova remnant, Puppis A. As can be deduced, e.g., from a discussion in Arendt et al.  (adp91 ), its age ranges from 2 to 5 kyr; the central value is kyr (Winkler et al. winkler88 ). For 1E 1207.4–5209, we also adopt the age of the host supernova remnant (G296.5+10). According to Roger et al. (roger88 ), it is kyr. For the Vela pulsar, we take the age interval from the standard characteristic spindown age of the pulsar to the characteristic age corrected due to the pulsar glitching behaviour (Lyne et al. lyne96 ). The age of PSR J0538+2817, kyr, was estimated by Kramer et al. (kramer03 ) from the measurements of the pulsar proper motion relative to the center of the host supernova remnant, S147. The age of RX J1856.4–3754 has been revised recently by Walter & Lattimer (wl02 ) from the kinematical reasons. Following these authors we take the central value kyr and choose such an error-bar of to clearly distinguish the revised value from the value kyr reported previously by Walter (walter01 ) on the basis of less accurate parallax measurement. The characteristic age of RX J0720.4–3125 has been estimated by Zane et al. (zane02 ), Kaplan et al.  (kaplanetal02 ) and Cropper et al. (cropper ) from X-ray measurements of the neutron-star spindown rate. We adopt the central value kyr with an uncertainty by a factor of 2. The ages of three other pulsars, PSR B1706–44, Geminga, and PSR B1055–52, are the characteristic pulsar ages assuming an uncertainty by a factor of 2.

For two youngest sources, RX J0205+6449 and the Crab pulsar, no thermal emission has been detected, but the upper limits on surface temperature have been established (Slane et al. slane02 , Weisskopf et al. weisskopf04 ). The surface temperatures of the next five sources, RX J0822–4300, 1E 1207.4–5209, Vela, PSR B1706–44, and PSR J0538+2817, have been obtained using hydrogen atmosphere models (see references in Table 1). Such models are more consistent with other information on these sources (e.g., Pavlov et al. pavlovetal02 ) than the blackbody model. On the contrary, for the Geminga and PSR B1055–52 we present the values of inferred using the blackbody spectrum because this spectrum is more consistent for these sources.

Let us notice that from Table 1 we have excluded PSR B0656+14 which was considered earlier (e.g., Yakovlev et al. yakovlevetal02 ). A combined analysis of new X-ray and optical observations of the source (with the improved distance from new parallax measurements of Brisken et al. briskenetal03 ) leads either to unrealistically small values of the neutron star radius (in the blackbody model) or to unreasonably small distance to the star (in the hydrogen atmosphere model); see, e.g., Zavlin & Pavlov (zp02 ). This makes current interpretations of the data unreliable.

The surface temperature of RX J1856.4–3754 is still rather uncertain. A wide scatter of , obtained by different authors, takes place because X-ray and optical observations are not described by one blackbody model. This can be explained, for instance, by the presence of hot spots on the neutron star surface. Thus, we adopt the upper limit MK, which agrees with the value of obtained either with the “Si-ash” atmosphere model of Pons et al. (ponsetal02 ) or with the model of condensed surface layers of Burwitz et al. (burwitzetal03 ). It agrees also with the model of nonuniform surface temperature distribution suggested by Pavlov & Zavlin (pz03 ). In the latter case, the mean surface temperature MK is below our upper limit of .

Finally, for RX J0720.4–3125 is taken from Motch et al. (motchetal03 ) who have interpreted the observed spectrum with a model of a hydrogen atmosphere of finite depth.

For PSR J0538–4300, PSR B1055-52, and RX J0720.4–3125, the authors cited in Table 1 have not reported any error bars of . We adopt 20% uncertainties which seem to be appropriate for these sources.

## 3 Physics input and calculations

We will simulate cooling of neutron stars using our generally relativistic cooling code described by Gnedin et al. (gyp01 ). We adopt a moderately stiff equation of state of neutron star interiors proposed by Douchin & Haensel (dh01 ). According to this equation of state, neutron star cores (regions of the densities g cm) consist of neutrons, with the admixture of protons, electrons and muons. All constituents exist everywhere in a core, except for muons which appear at g cm. The most massive stable star has the (gravitational) mass , the central density g cm, and the (circumferential) radius km. The central densities and masses of eight neutron star models (with from 1.111  to 1.994 ) are presented in the right panel of Fig. 1.

All physics input is standard. The effects of muons are included as described by Bejger et al. (byg03 ). We assume no envelope of light elements on stellar surfaces (Sect. 5). The code calculates the cooling curves, which give the dependence of the effective surface stellar temperature on stellar age . Let us remind that neutron stars are born hot in supernova explosions (with internal temperatures K) but gradually cool down via neutrino emission from the entire stellar body and via heat diffusion to the surface and thermal surface emission of photons. Qualitatively, one can distinguish three cooling stages. At the first (‘non-isothermal’) stage ( yr) the main cooling mechanism is neutrino emission but the stellar interior stays highly non-isothermal. At the second (‘neutrino’) stage ( yr) the cooling goes mainly via neutrino emission from isothermal interiors. At the third (‘photon’) stage ( yr) a star cools predominantly through the surface photon emission. Figure 1: Left: Density dependence of critical temperature of model p1 for proton superfluidity and model nt1 for neutron superfluidity in a neutron-star core; vertical dot-and-dash line indicates the central density of a maximum-mass neutron star. Right: Observations (Table 1) compared with theoretical cooling curves of eight neutron stars (1–8) with different masses. All solid curves refer to neutron stars with model superfluidities from the left panel. The dotted curve 7 is for a non-superfluid star. Insert table gives masses and central densities of stars 1–8.

The new element of our present studies is the equation of state of Douchin & Haensel (dh01 ). We have chosen it because it forbids the powerful direct Urca process of neutrino emission (Lattimer et al. lpph91 ) in all stable neutron stars (). In this case, a non-superfluid neutron star of any mass will have almost the same (universal) cooling curve (the dotted curve in the right panel of Fig. 1). At the neutrino cooling stage, this curve is determined by the neutrino emission due to the modified Urca process. The curve is almost independent of the equation of state of neutron star cores (Page & Applegate pa92 ) as long as the direct Urca process is forbidden. As has been indicated by many authors (see, e.g., Yakovlev & Pethick yp04 and references therein) and seen from Fig. 1, this universal cooling model is certainly unable to explain the data. For instance, it gives much lower than that of PSR B1055–52, but much higher than that of the Vela pulsar. We will show that all the data can be explained assuming superfluidity of neutron-star cores.

It is well known that neutrons and protons in stellar cores can be in superfluid state. Proton superfluidity is caused by singlet-state proton pairing, while neutron superfluidity is produced by triplet-state neutron pairing. These superfluidities can be specified by density dependent critical temperatures for protons and neutrons, and . Results of calculations of these temperatures from microscopic theories show a large scatter of critical temperatures depending on a nucleon-nucleon interaction model and a many-body theory employed. In particular, recently Schwenk & Friman (sf04 ) and Zuo et al. (zuo ) have obtained weak neutron and proton pairing in neutron star cores but many other calculations give much stronger superfluidity (e.g., Lombardo & Schulze ls01 ; also see references in Yakovlev et al. yls99 ). In this situation it is reasonable to consider and as unknown functions of (consistent with predictions of microscopic theories) which can hopefully be constrained by comparing theoretical cooling curves with observations.

Superfluidity of neutrons and/or protons in neutron-star cores affects the heat capacity of nucleons and reduces neutrino reactions (Urca and nucleon-nucleon bremsstrahlung processes) involving superfluid nucleons (as reviewed, e.g., by Yakovlev et al. yls99 ). Moreover, superfluidity initiates an additional neutrino emission mechanism associated with Cooper pairing of nucleons (Flowers et al. frs76 ). All these effects of superfluidity are incorporated into our cooling code.

In our calculations we adopt one model of strong superfluidity of protons (with the maximum of about K) and several models of moderate superfluidity of neutrons (with K) in a neutron-star core. These models are phenomenological but consistent with the results of microscopic theories. A pair of models: proton superfluidity p1 and neutron superfluidty nt1 is plotted in the left panel of Fig. 1.

The strong proton superfluidity is required to slow down the cooling of low-mass stars, , whose central densities are g cm. This scenario was suggested by Kaminker et al. (khy01 ). In a low-mass star, one has  K everywhere in the core. The proton superfluidity occurs at the early cooling stage ( yr) and suppresses modified Urca processes of neutrino emission as well as neutrino generation in proton-proton and proton-neutron collisions. Neutrino emission due to Cooper pairing of protons is switched on too early and becomes inefficient in middle-aged neutron stars we are interested in. In contrast, the adopted neutron superfluidity is too weak in low-mass stars (the left panel of Fig. 1) to appear at the neutrino cooling stage. This superfluidity does not suppress the neutrino emission in neutron-neutron collisions which becomes the leading mechanism of neutrino cooling. It is much weaker than the modified Urca process (which would be leading in non-superfluid stars). As a consequence, the cooling curves of low-mass stars go noticeably higher than the universal cooling curve of non-superfluid stars. Actually, these cooling curves also merge into one almost universal curve, which is independent of the equation of state in a stellar core and the exact behaviour of the (Kaminker et al. kyg02 ). This upper curve 1 allows one to explain observations of neutron stars hottest for their age (RX J0822–4300, 1E 1207.4–5209, PSR B1055–52, RX J0720.4–3125) as cooling low-mass neutron stars.

Now we come to observations of neutron stars coldest for their age (first of all, PSR J0205+6449, the Vela pulsar, and Geminga). It has been widely proposed to interpret these objects as rather massive neutron stars with the neutrino emission enhanced by direct Urca process in nucleon cores (or by similar processes in pion-condensed, kaon-condensed or quark cores). We will show that coldest objects can be explained without invoking these mechanisms by tuning the model of moderate neutron superfluidity at g cm. Let us consider the most massive neutron star (, curve 8) in Fig. 1. Its central density is higher than the density at which neutron superfluidity nt1 dies out. When the internal temperature of the star becomes lower than the maximum critical temperature of the neutron superfluidity, the neutrino emission due to Cooper pairing of neutrons switches on and becomes a powerful neutrino emission mechanism, which can be about two orders of magnitude more efficient than the modified Urca process in a non-superfluid star (see Sect. 4). This emission produces enhanced cooling (attributed to direct Urca or similar processes in previous calculations). The enhancement is not too strong (e.g., the direct Urca process in a nucleon stellar core would further enhance the neutrino luminosity by about 4–5 orders of magnitude). However, even this not too strong enhancement is sufficient to explain observations of the coldest neutron stars (particularly, PSR J0205+6449, the Vela and Geminga pulsars). Evidently, all neutron stars with g cm (in our model) will cool nearly as fast as the star in Fig. 1.

Therefore, we come to three distinct classes of cooling neutron stars (similar to those described by Kaminker et al. kyg02 for the case of enhanced cooling due to direct Urca process). The first class contains low-mass, very slowly cooling stars (curve 1 in the right panel of Fig. 1). Another class contains high-mass stars with enhanced cooling (curve 8). Finally, there is a class of medium-mass neutron stars (curves 2–6) which show intermediate cooling. Their cooling curves fill in the space between the upper curve for low-mass stars and the lower curve for high-mass stars. These curves explain observations of PSR B1706–44, PSR J0538+2817, and RX J1856.4–3754.

## 4 Cooper-pairing neutrino emission as a fast-cooling agent

Let us give a simple explanation of the computer results on enhanced neutrino emission due to Cooper pairing of neutrons. We start from the expression for the neutrino emissivity due this process (e.g., Eq. (236) in Yakovlev et al. ykhg01 ). It can be written as

 QCP(ρ,T)=q(ρ,T)F(τ), (1)

where

 q(ρ,T) ≈ 1.17×1021(m∗NmN)(pFmNc) (2) ×T79NνaNergcm−3s−1,

K is the internal stellar temperature, is the bare nucleon ( or ) mass, is the nucleon effective mass in dense matter, is the nucleon Fermi momentum, is a dimensionless constant combined of squared weak-interaction constants of vector and axial-vector nucleon currents, is the number of neutrino flavors, and is a function of . The constant depends on nucleon species and pairing type, while depends on pairing type. We have for the triplet-state neutron pairing under discussion. This value can be renormalized by many-body effects (for instance, the renormalization of the axial-vector constant was considered by Carter & Prakash cp02 ). However, theoretical cooling curves are not too sensitive to the exact value of , and we use the non-renormalized value. The analytic fit expression for is presented, for instance, by Yakovlev et al. (ykhg01 ). Let us remind that just after superfluidity onset (immediately after falls below ) and at . Thus, the emissivity is exponentially suppressed at . Figure 2: A sketch of neutrino luminosities produced by the modified Urca process (LMurca) and Cooper pairing process (LCP) as well as of the photon luminosity Lγ of a neutron star versus internal temperature T for three models of neutron superfluidity in the stellar core with Tmaxcn=108, 3×108 and 109 K.

For our qualitative analysis in this section we employ the simplest dependence of the neutron critical temperature on distance from the stellar center:

 Tcn(r)=Tcm{1−(r−rm)2(Δrm)2} (3)

at (with the maximum at ), and at .

Neglecting, for simplicity, general relativistic effects and assuming an isothermal stellar core at a temperature , the neutrino luminosity due to Cooper pairing of neutrons can be written as

 LCP=4π∫r2r1r2QCPdr. (4)

Here, and restrict the superfluid layer, where and the neutrino process in question is allowed. To be specific, let us assume that the widest superfluid layer (which is realized at and extends from to ) entirely falls in the neutron star core.

The factor in the emissivity , Eq. (1), is a more rapidly varying function of than . Thus we can set and (with ) in all functions under the integral but in . A simple replacement of integration variable leads then to

 LCP = 8πr2mΔrmq(ρm,T)τmℓ(τm), (5) ℓ(τ) = 12∫1τdτ′F(τ′)τ′3/2√τ′−τ, (6)

where . The integration can be done numerically; the appropriate analytic fit (for triplet-state neutron pairing) is

 ℓ(τ) = (1−τ)3/2[3.844(1−τ)+3.142τ2 (7) + 13.99τ(1−τ)+25.4τ2.5(1−τ)2((τ−0.2493)2+0.03694)0.7].

Evidently, the luminosity vanishes in a hot star, where and neutron superfluidity is absent. It switches on as falls below ; it grows almost linearly while decreases to ; afterwards, it reaches maximum at (with ) and then decreases. At the increasing and maximum-luminosity stage, is collected from a superfluid spherical stellar layer in the vicinity of the maximum critical temperature, . This creates a splash of neutrino emission associated with Cooper pairing of neutrons.

For typical values of the parameters, the maximum value of can be one-two orders of magnitude higher than the neutrino luminosity of a non-superfluid star (with forbidden direct Urca process). This is demonstrated in Fig. 2 using a toy model of cooling neutron stars described by Yakovlev & Haensel (yh03 ) – there is no need to employ accurate models in this section. The parameters of the neutron-star model presented at the figure are: , km, g cm, km. The three superfluidity models () are selfsimilar and differ by the values of , , and K. Three solid lines exhibit the Cooper-pairing neutrino luminosity calculated from Eqs.  (4)–(7) for three models of neutron superfluidity. Since and , the Cooper-pairing luminosity is more competitive at weaker superfluidity (lower ). However, at K this luminosity becomes lower than the photon thermal luminosity of the star (Fig. 2) which makes it insignificant for stellar cooling. It is worth to notice that, for realistic parameters, is much smaller than the neutrino luminosity due to the direct Urca process in a non-superfluid star (if the direct Urca process is open).

The decreasing part of is even more fascinating. We have as , resulting in the scaling relation

 LCP∝ΔrmT8/Tcm, (8)

which becomes sufficiently accurate at . This neutrino emission is actually produced from two thin spherical shells (near and ), where is just below . The widths of these shells are proportional to , which explains the power-law (instead of the exponential decrease of the emissivity in a local element of superfluid matter). Therefore, the decreasing part of the Cooper-pairing neutrino luminosity has the same temperature dependence as all slow neutrino emission mechanisms (modified Urca, nucleon-nucleon bremsstrahlung) in non-superfluid cores of stars with forbidden direct Urca processes. In other words, the superfluidity suppresses the neutrino emission available in non-superfluid stars but initiates the Cooper-pairing neutrino emission in such a way that it acts as a new nonsuppressed neutrino cooling mechanism. Moreover, the new emission can be more intense than that in a non-superfluid star and provide enhanced cooling. This important feature appears in realistic models of cooling neutron stars with density dependent critical temperatures (and does not appear in the models with density-independent ). In particular, it implies that once takes on leadership in competition with just after the superfluidity onset, it will not lose it during subsequent evolution (especially because is actually suppressed by superfluidity, which is not taken into account in Fig. 2). This is clearly seen from Fig.  2.

Let us add that at we can obtain a better formula for than Eq. (4), without employing the specific profile, Eq. (3). It is sufficient to start from Eq. (4) and notice that the main contribution into comes from two thin shells, at and , where . In each shell, the gradient can be taken constant. Then we get

 LCP=8π[r21H1q(ρ1,T)+r22H2q(ρ2,T)]ℓ(0), (9)

where and are characteristic widths of our shells, , , and . Strictly speaking, , , , , , and depend slightly on , but this dependence can be regarded as parametric. It is easy to verify that if is given by Eq. (3) at and , our new expression for coincides with Eq. (5). Equation (5) is expected to be useful just after the superfluidity onset, at (where the parabolic dependence may be a good approximation), while Eq. (9) is more exact at lower . Both equations enable one to incorporate the Cooper-pairing neutrino emission in simplified cooling models (like a toy model of Yakovlev & Haensel yh03 ), useful for understanding main features of neutron star cooling without complicated cooling codes.

The above analysis is valid as long as vanishes in the stellar interior. If not, there is a minimum value of , and will become exponentially suppressed at .

## 5 Testing the cooling scenario and discussion

After clarifying the efficiency of the Cooper-pairing neutrino emission let us return to the cooling scenario described in Sect. 3. As we have already mentioned, the scenario is rather insensitive to a specific model of proton superfluidity (required to raise the surface temperature of low-mass stars for explaining observations of the sources hottest for their ages). The only serious constraint on the proton pairing is that should be high ( K) in the cores of low-mass stars. Figure 3: Left: One model p1 for proton superfluidity and five models nt1–nt5 for neutron superfluidity in a neutron-star core. Right: Cooling curves of low-mass (1.111 M⊙) and high-mass (1.994 M⊙) stars with model p1 proton superfluidity and one of the models of neutron superfluidity from the left panel, compared with observations. Cooling of the low-mass star is insensitive to selected models of neutron superfluidity (except for model nt2 at t>300 kyr). Insert shows the comparison of cooling curves of the high-mass star with observations of the Vela pulsar in more details.

However, the constraints on the neutron critical temperature in a stellar core should be really strong. This is illustrated in Fig. 3. The left panel displays the critical temperatures of our proton superfluidity model (p1) and five neutron superfluidity models (nt1–nt5), including our basic model nt1 used in Sect. 3. The right panel shows cooling curves of a low-mass () star and a high-mass () star. Any curve is calculated for model p1 of the proton superfluidity and one model of the neutron superfluidity from the left panel of Fig. 3. Any observational point between an upper curve and a lower curve can be explained by a given superfluid model. The constraints on the neutron superfluidity are as follows.

(1) The neutron superfluidity should be weak in low-mass stars. In our case (for the equation of state of Douchin & Haensel dh01 ) this means that K at g cm. Under this condition the neutron superfluidity does not affect the cooling (at least at the neutrino cooling stage) of low-mass stars () and does not violate our interpretation of the sources hottest for their age (first of all, RX J0822–4300 and PSR B1055–52). Accordingly, all five cooling curves (for superfluids nt1–nt5) of low-mass stars merge in one upper (solid) cooling curve in Fig. 3. The only exclusion is provided by model nt2 with highest pre-peak among models nt1–nt5. In a low-mass star this superfluidity occurs at kyr. The Cooper-pairing neutrino emission and reduced heat capacity of neutrons noticeably accelerate the cooling at this late stage (the upper short-dashed curve).

(2) The neutron superfluidity should be moderate at g cm, with the peak maximum K (model nt1 in Fig. 3, the solid curve). In this case it switches on just in time to initiate the enhanced cooling in a high-mass star. Its level is sufficient to explain observations of neutron stars coldest for their ages (first of all, PSR J0205+6449 and the Vela pulsar). The asymptotic neutrino-cooling regime given by the scaling expression (8) is realized at kyr. If were slightly higher than K (model nt3, K, the dotted curve), the Cooper-pairing neutrino emission will start operating in a younger massive star but becomes less efficient at kyr, which is less favorable for explaining the observations of the Vela pulsar. This cooling behaviour is naturally explained by the scaling (8). If were slightly lower than K (model nt4, K, the long-dashed curve), the Cooper-pairing neutrino emission will start operating too late which would violate the interpretation of the observations of PSR J0205+6449.

(3) The results are also sensitive to the width of the peak of the curve. For instance, retaining the peak maximum of K but making the peak narrower (model nt5, the dot-dashed curve) will reduce the neutrino emissivity due to neutron pairing, raise the temperature of the massive star and complicate the interpretation of the Vela pulsar (again, in agreement with the scaling (8)). However, the cooling curves are rather insensitive to the exact position of the maximum. We can slightly shift the maximum to higher or lower (confining the peak within the kernel of a massive star) but these shifts will not change the cooling curves of massive stars (such tests are not shown in Fig. 3). However, the shift of the maximum to g cm would cause the enhanced cooling of low-mass stars. The cooling curves of low-mass stars would become close to those of high-mass stars which would violate the interpretation of neutron stars hottest for their ages (see item (1)).

This discussion shows that the cooling curve of a massive neutron star implying model nt1 of neutron superfluidity is close to the lowest cooling curve (in the scenario, where the cooling is enhanced by Cooper-pairing neutrino emission). Observations of cold neutron stars, PSR J0205+6449 and the Vela pulsar, provide excellent tests for this scenario. If these pulsars were noticeably colder we would be unable to explain them within our scheme. Notice, that the upper limit of the surface temperature of PSR J0205+6449 was inferred from observations (Slane et al. slane02 ) using the blackbody spectrum of surface emission. If this pulsar has a hydrogen atmosphere, the upper limit on could be expected to be about twice lower than for the blackbody case. In that case we would be unable to explain this source within the proposed scenario.

Although we have used one equation of state of dense matter (Douchin & Haensel dh01 ) we would obtain similar results for other equations of state which forbid direct Urca processes (and other similar processes of fast neutrino cooling) in neutron star cores. Taking different equations of state would lead to attributing different masses to the same sources (Fig. 1); similar problem has been discussed by Kaminker et al. (kyg02 ).

In addition, we could take an equation of state in the stellar core which opens direct Urca process at highest densities (in the central kernels of most massive stable neutron stars; similar to the equation of state of Akmal & Pandharipande ap98 ). Applying the same model of nucleon superfluidty as in Fig. 1, we would get five types of cooling neutron stars (instead of three). Three types would be the same as those mentioned in Sect. 3: low-mass, very slowly cooling stars; massive stars whose cooling is enhanced by Cooper-pairing neutrino emission; and of medium-mass stars whose cooling is intermediate. In addition, we would have: most massive neutron stars demonstrating a very fast cooling via the direct Urca process; and stars whose cooling is intermediate between that enhanced by the Cooper-pairing neutrino emission and by the direct Urca process. The transition from the Cooper-pairing neutrino cooling to direct-Urca cooling with increasing mass will be very sharp and the number of intermediate-cooling sources will be small. The maximum-mass neutron stars would be extremely cold ( K at kyr), about the same as discussed, e.g., by Kaminker et al. (kyg02 ). A discovery of such stars would definitely indicate the operation of the direct Urca process in their cores. An indirect evidence of their existence is provided by a non-detection of neutron stars in some supernova remnants (Kaplan et al. kaplan ).

Note that the cooling of neutron stars can also be affected by the singlet-state superfluidity of neutrons in inner stellar crusts, by the presence of surface layers of light (accreted) elements, and by stellar magnetic fields (e.g., Potekhin et al. pycg03 , Geppert et al. geppert ). These effects can be especially important in low-mass stars. We have neglected them in the present paper since we have mainly focused on enhanced cooling of massive stars but we will consider them in a future publication.

## 6 Conclusions

We have proposed a new scenario of cooling of isolated neutron stars. We have shown that the present observational data on thermal emission from isolated middle-aged neutron stars can be explained assuming that neutron star cores are composed of neutrons, protons and electrons (and possibly muons) with forbidden direct Urca process of neutrino emission. In our scenario, an enhanced neutrino emission, which is required for interpretation of neutron stars coldest for their age, is provided by neutrino process associated with Cooper pairing of neutrons. We have shown that the neutrino luminosity due to this process (at internal temperatures ) behaves as . In this way it “mimics” the neutrino luminosity produced either by modified Urca processes or nucleon-nucleon bremsstrahlung processes in non-superfluid stars, but it can be one-two orders of magnitude higher. The proposed cooling scenario imposes very stringent constraints on the density dependence of neutron-pairing temperature . The constraints result from the comparison of the cooling theory with two most important “testing sources”, PSR J0205+6449 and the Vela pulsar (Sect. 5). This scenario is the first one in which a moderate superfluidity and associated neutrino emission are helpful for explaining the data (cf. with previous cooling scenarios, where a moderate superfluidity has violated interpretation of observations, e.g., Kaminker et al. kyg02 ).

Our interpretation implies the presence of a strong proton superfluidity and a moderate neutron superfluidity in neutron star cores (Sect. 3). We need the proton superfluidity to explain observations of neutron stars hottest for their age, and the neutron superfluidity to explain observations of stars coldest for their age. However, as has been demonstrated by Gusakov et al. (gkyg04 ), cooling curves are not too sensitive to exchanging neutron and proton superfluidities () in neutron-star cores. Therefore, we would also be able to explain observational data in the scenario with a strong neutron superfluidity and a moderate proton superfluidity in stellar cores.

We need a strong superfluidity to suppress modified Urca process in low-mass stars, rise the surface temperature of these stars and explain observations of neutron stars hotter for their age. In fact, we can rise the temperature of low-mass middle-aged neutron stars by assuming the presence of surface layers of light (accreted) elements. The mass of light elements may decrease with time, e.g., due to diffusive nuclear burning (Chang & Bildsten cb03 ), which opens additional freedom to regulate the cooling. In this way, the presence of a strong (proton or neutron) superfluidity in a neutron star core is not vitally important for our interpretation. We will show this in a future publication. However, the presence of a moderate superfluidity (of neutrons or protons) with a tuned density dependence of critical temperature (Sect. 5) is crucial for this scenario, where this tuned dependence is combined with the remarkable simplicity of the equation of state of neutron-star cores (nucleon composition with forbidden direct Urca process). We hope that this scenario can be taken into consideration along with many other scenarios (reviewed or proposed, e.g., by Page page98a ; page98b , Tsuruta et al. tsurutaetal02 , Khodel et al. khodeletal04 , Blaschke et al. blaschke , Yakovlev & Pethick yp04 ). The correct scenario should be selected in future observations of neutron stars combined with new advanced theoretical results.

After this paper was prepared for submission we became aware of the paper of Page et al. (pageetal04 ). These authors give a detailed consideration of enhanced cooling via neutrino emission due to Cooper pairing of neutrons in neutron-star cores composed of nucleons with forbidden direct Urca process. The idea to enhance the cooling by Cooper-pairing neutrino emission is the same as in our paper, but its realization is different. Particularly, Page et al. (pageetal04 ) use a set of superfluidity models obtained from microscopic theories. Their main models for neutron superfluidity in a stellar core (for, instance, model (a) in their Fig. 9) have too high peak temperatures K and too high at the pre-peak densities to explain the observations of PSR J0205+6449 and the Vela pulsar and to obtain a pronounced dependence of cooling curves on neutron star mass. In contrast, our models are phenomenological but tuning them we obtain a noticeable dependence of the cooling on . It enables us to explain all the data by one model of nucleon superfluidity (even neglecting the effect of accreted envelopes).

###### Acknowledgements.
We are grateful to A.Y. Potekhin for helpful discussions. This work has been supported partly by the RFBR, grants 02-02-17668 and 03-07-90200, the RLSS, grant 1115.2003.2, and by the INTAS, grant YSF 03-55-2397.
{"ft_lang_label":"__label__en","ft_lang_prob":0.86063445,"math_prob":0.8598341,"size":41734,"snap":"2023-40-2023-50","text_gpt3_token_len":10919,"char_repetition_ratio":0.16101126,"word_repetition_ratio":0.047249854,"special_character_ratio":0.26053098,"punctuation_ratio":0.18071856,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.9512607,"pos_list":[0,1,2,3,4,5,6],"im_url_duplicate_count":[null,1,null,1,null,1,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2023-10-04T16:13:11Z\",\"WARC-Record-ID\":\"<urn:uuid:02564997-884d-4266-901d-a72b18a36e72>\",\"Content-Length\":\"541025\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:3a384a6f-e71f-4e55-903e-59f0628e8a20>\",\"WARC-Concurrent-To\":\"<urn:uuid:ed58d4bb-4291-4fac-b8d1-43ff87c15b0a>\",\"WARC-IP-Address\":\"172.67.158.169\",\"WARC-Target-URI\":\"https://www.arxiv-vanity.com/papers/astro-ph/0404002/\",\"WARC-Payload-Digest\":\"sha1:NPTJR7N2HLGW5U6LI4VUEGYTCOKZR4EJ\",\"WARC-Block-Digest\":\"sha1:IWQH2ZH36IJ7CRHNSSUYNODO3MT4QG5A\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2023/CC-MAIN-2023-40/CC-MAIN-2023-40_segments_1695233511386.54_warc_CC-MAIN-20231004152134-20231004182134-00755.warc.gz\"}"}