Draft version May 21, 2020 Typeset using LATEX twocolumn style in AASTeX63
NS 1987A in SN 1987A Dany Page,1 Mikhail V. Beznogov,1 Ivan´ Garibay,1 James M. Lattimer,2 Madappa Prakash,3 and Hans-Thomas Janka4
1Instituto de Astronom´ıa,Universidad Nacional Aut´onoma de M´exico, Ciudad de M´exico, CDMX 04510, Mexico 2Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA 3Department of Physics and Astronomy, Ohio University, Athens, OH 45701, USA 4Max-Planck-Institut f¨urAstrophysik, P.O. Box 1317, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany
(Received xx; Revised xx; Accepted xx)
ABSTRACT The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object’s predicted position. The infrared excess could be due to the decay of isotopes like 44Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spindown of a pulsar, or thermal emission from an embedded, cooling neutron star (NS 1987A). It is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 year old neutron star. Furthermore, models of cooling neutron stars within the Minimal Cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively 1 limited crustal n- S0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations don’t overlap, then pulsar spindown or accretion might be more likely, but the pulsar’s period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.
Keywords: stars: neutron - supernovae: individual (SN 1987A)
1. INTRODUCTION The chance to witness the birth of a neutron star arose when neutrinos were detected from the core-collapse su- arXiv:2004.06078v2 [astro-ph.HE] 20 May 2020 pernova SN 1987A on Feb 23, 1987 (Hirata et al. 1987; [email protected] Bionta et al. 1987). The further opportunity to observe [email protected] the early evolution of a neutron star has, however, been elusive. Owing to the dust and the ring surrounding the [email protected] supernova remnant (SNR) of SN 1987A, observational [email protected] searches in radio and X-rays for a neutron star remnant [email protected] have been unsuccessful (e.g., Alp et al. 2018; Esposito et al. 2018). Recently, however, Cigan et al.(2019) have [email protected] presented high angular resolution images of dust and molecules in SN 1987A ejecta obtained from the At- 2 Page et al. acama Large Millimeter/submillimeter Array (ALMA) (2016). The best-fit model, L15-1, had also been con- and concluded that the presence of a compact source sidered (Abell´anet al. 2017) with respect to the dis- in the remnant is strongly indicated. They observed a tribution of molecular CO 2-1 and SiO 5-4 emission in localized blob of warm dust with temperature 33 K the ejecta of SN 1987A, and well-fits the size, shape ' compared to that of the surrounding dust, (17–22) K. and clumpy character of its apparent ring geometry. ' The observed SiO and CO gas temperatures correspond By orienting the Fe+Si ejecta of this model with re- to a luminosity of the dust blob of Lcs = (40–90)L , spect to the ring plane and observer direction to obtain 33 1 where L = 3.826 10 erg s− is the solar luminosity the asymmetry seen in SN 1987A, the supernova kick × (the subscript ‘cs’ stands for ‘condensed source’), which turned out to have a northern component. The main requires a compact source to have an estimated power of reason for this is a big Fe+Si mass located below the this magnitude or slightly lower, if the source is embed- ring plane of SN 1987A, south of the connecting line to ded in the dust blob. However, if the compact source is the observer. 3D explosion simulations show that the not located within the blob, but heats it from afar, the supernova kick vector and the bulk mass of iron-group source must have a power somewhat greater than Lcs . and intermediate-mass elements should lie in opposite In proposing a compact object in the remnant of SN hemispheres (Wongwathanarat et al. 2013), compatible 1987A as the most likely explanation for the observed with observations (Katsuda et al. 2018; Holland-Ashford excess dust blob emission, Cigan et al.(2019) had to rule et al. 2017). The same 3D supernova model, L15-1, also out alternative possibilities, the foremost of which was allowed for a reasonable match of the redshift of the 44Ti heating by the decay of 44Ti synthesized by SN 1987A. emission and of the 56Co 847 keV line profile observed While the warm, extended dust is expected to be ra- in SN 1987A (Jerkstrand et al. 2020). The resulting dioactively heated by 44Ti, this is unlikely for the con- angle between the supernova kick vector and observer centrated, warmer blob. The existence of a single, high- direction should be about 30 degrees, and most likely density clump of 44Ti seems implausible, and even if it less than about 90 degrees (Jerkstrand et al. 2020). The were formed, its heating would not be strongly localized model L15-1 had a kick velocity only around 300 km/s, because of optically thin conditions to the 44Ti γ-rays while the actual kick may have been 500 km/s or more, as discussed in Section2. but the additional expense of fine-tuning was considered Cigan et al.(2019) also noted that there is an offset not worthwhile. between the location of the brightest pixel of the warm That the explosion was significantly asymmetric is not blob and the center of the SNR at the original posi- in doubt, since radioactive 56Co debris is seen mostly tion of the progenitor star. This displacement, which moving away from us (McCray & Fransson 2016). A could be associated with the supernova kick imparted NuSTAR observation (Boggs et al. 2015) also shows that to the compact object, is between about 20 mas and 85 44Ti is considerably redshifted, suggesting a kick veloc- mas, depending on how the center of the explosion is ity component along the line of sight towards us of sev- determined, e.g., by fitting the geometrical center of the eral hundred km/s. With a transverse component of 315 GHz emission seen by ALMA (Cigan et al. 2019), 160 km/s the compact source in SN1987A would have a the ATCA radio ring continuum (Potter et al. 2009), space velocity near the peak of the distribution observed or the ring hot spots on HST images (Alp et al. 2018). for young pulsars, whereas a transverse component of If connected with a supernova kick, the velocity com- nearly 700 km/s would place it in the high-velocity tail. ponent transverse to the line of sight is between 160 Explosion models of Utrobin et al.(2019) for state- km/s and nearly 700 km/s (for the 51.4 kpc distance to of-the-art progenitor models of SN 1987A indicate the the LMC). In addition, the kick, judging from the dis- baryon mass, MB, of its compact remnant to be (1.35– tribution of iron-group and intermediate-mass elements 1.66) M , while Ertl et al.(2020) predict (1 .48–1.56) M
(Larsson et al. 2016), should have a northerly compo- for single-star progenitors and (1.38–1.75) M for bi- nent in the sky (Janka et al. 2017), which matches the nary progenitor1. These baryon masses translate to a orientation of the dust blob with respect to the origi- gravitational mass M (1.22–1.62) M using the EOS- ' nal position of the progenitor. This offers evidence that the compact source is nearby or even surrounded by the 1 We have dropped results with the progenitors B15 in Utrobin blob. et al.(2019) and W15 in Ertl et al.(2020) (which is the same Specifically, Janka et al.(2017) analysed the geome- model) because its He core mass is too small to explain the light try of Fe and Si in a set of 3D supernova simulations curve peak, its pre-SN luminosity is too small, and it ejects too little O. Additionally, we omitted cases from Utrobin et al.(2019) for matching the shape and mass of the Fe+Si distri- with too-little ejected nickel or too-high explosion energies. bution of SN 1987A as determined by Larsson et al. NS 1987A in SN 1987A 3 independent relation (Lattimer & Prakash 2001) the pulsar wind or accretion explanation might be pre- ferred since only a fraction of a source’s power would MB M β − (1.2 0.1) , (1) be required. Either could easily be large enough (with M ' ± 2 β − plausible P and B, or accretion rates) while the expected where β = GM/Rc2 and R 11.5 1 km is the typical L would be quite insufficient. But it should be noted ' ± th neutron star radius2. These values are well below the that past observations (Alp et al. 2018) have set upper measured masses, M > 2 M , of several pulsars (PSR limits to the total (bolometric) emission of any kind of
J1614-2230, Demorest∼ et al. 2010; PSR J0348+0432, compact source of about 138 L in the presence of dust,
Antoniadis et al. 2013; and PSR J0740+6620, Cro- and 22 L without dust. Even with dust, this is 1000 martie et al. 2020), as well as an inferred upper limit times smaller than the Crab’s luminosity. An otherwise to the neutron star maximum mass Mmax < (2.2– hidden pulsar or accretion source can thus have a lumi- ∼ 2.3) M (Margalit & Metzger 2017) from GW170817, nosity at most 1.5 Lcs – 3.5 Lcs , so it must be unnatu- which strongly suggests that a black hole remnant in rally close to the blob and again raises the prospect of SN 1987A is unlikely. The identification of the com- fine-tuning that would disfavor these hypotheses. Sec- pact remnant with an intermediate-mass neutron star tion3 considers the possibility that the required energy is supported by the observations of neutrinos from SN stems from the spin-down of a young pulsar. 1987A and the association of their inferred total energy Our preferred hypothesis is, instead, that the power (Lattimer & Yahil 1989) with the binding energy of a source of the blob is NS 1987A, a central compact object proto-neutron star. If about 1/6 of the total energy was (CCO: Pavlov et al. 2002), defined to be a young neutron radiated as the observed electron anti-neutrinos, the im- star in a SNR whose luminosity, Lth, is predominantly plied (2.9 1.2) 1053 ergs binding energy suggests the due to surface thermal emission. Section4 examines the ± × neutron star gravitational mass to be (1.38 0.43) M cooling of a star following the “Minimal Cooling” sce- ± using Eq. (1). We assume in this paper that the com- nario (Page et al. 2004, 2009, 2011) which assumes the pact remnant produced by SN 1987A is most likely a lack of rapid neutrino cooling due to a direct Urca pro- neutron star, hereafter called NS 1987A, which is also cess (Lattimer et al. 1991), emphasizing the importance possibly a pulsar. of the envelope’s chemical composition. The question of If the neutron star is enclosed within the blob, the whether or not light elements can survive in the enve- most natural explanation is that the blob is heated by lope during the hot, early stages of a neutron star’s life its thermal emission Lth . As we show in Section4, the is addressed in Section5. Section6 considers the case expected Lth of a 30 year old neutron star is within a of a neutron star that has enhanced neutrino cooling, factor 3 of the inferred excess blob luminosity. possibly because it is relatively massive. In Section7, An alternative explanation, that the entire spindown a comparison of the cooling trajectories of NS 1987A power of a pulsar heats the surrounding blob, is disfa- and the neutron star, Cas A, in the Cassiopeia A SNR vored by the fine-tuning of the rotational period P and is made. Section8 contains a discussion and conclu- magnetic field B of the young neutron star that would sions. Essential details of the equation of state (EOS) be required. Both P and B could have values up to 2 models used in this work are given in AppendixA. The orders of magnitude higher or lower than what is nec- neutrino cooling processes considered are summarized essary. The same argument would apply for another in AppendixB, and neutron superfluid gaps used in the alternative, an accretion power source from a neutron inner crust are described in AppendixC. star or black hole. Here, a reasonable upper limit is 44 39 1 2. BLOB LUMINOSITY AND TI DECAY Lacc < 1.3(M/M ) 10 erg s− , which is 10 times the · 44 Eddington luminosity. But the lower limit is arbitrar- Radioactive decay of Ti might offer a possible expla- nation of the blob luminosity, which is L = (1.5–3.5) ily small, so the large range of possible accretion rates cs × 1035 erg s 1 (Cigan et al. 2019), if γ-rays and positrons implies fine-tuning is required to obtain Lacc Lcs . − ∼ 44 44 However, it is also possible that the blob and the produced through the decay channel of Ti Sc 44 → neutron star’s locations are disjoint, a situation well Ca were efficiently thermalized in the blob medium. → 44 44 known from the Crab Nebula, where the brightest part The decays of Ti to Sc proceed by electron capture, 44 44 of the pulsar wind and the pulsar are spatially separated and the transition of Sc to Ca is almost exclusively + (Weisskopf et al. 2000; Gomez et al. 2012). In this case, by β decays. A corresponding upper limit to the lumi- nosity is obtained for complete thermalization and given by 2 All radii quoted in this paper are circumferential radii, i.e., no dN (t) N (t) red-shift factor is applied. L (t) Ti E = Ti E . (2) blob ≤ dt decay τ decay 0
4 Page et al.
44 Above, NTi(t) = NTi,0 exp( t/τ0) is the time-dependent the observer (i.e., redshifted, consistent with the Ti 44 − number of Ti nuclei in the blob with NTi,0 being the decay-line measurements by Boggs et al. 2015) and along initial number and τ0 85 yr their decay time. The en- the symmetry axis of the ejecta that is inclined to the 44 ≈ ergy release per Ti decay, E 2.9 MeV includes east by about 14◦ relative to the symmetry axis of the decay ≈ the 0.068 MeV, 0.078 MeV, and 1.157 MeV γ-photons ring. This inferred apparent position of the 56Ni clumps from 44Sc and 44Ca de-excitation, as well as the energy and 44Ti heating north-east of the line of sight adds 2 + (2m c + E + ) from e e−-annihilation of the emitted relevance to the question of whether the blob emission e h e i positron, which possesses an average kinetic energy of is powered by 44Ti decay. E + 0.6 MeV (Cameron & Singh 1999). The blob However, in the following we will argue that this pos- h e i ∼ luminosity measured about t 30 yr after the explo- sibility is disfavored by the fact that the blob is not m ≈ sion thus leads to a constraint on the initial 44Ti mass massive enough to enable efficient thermalization of the contained by the blob: γ photons produced by radioactive decays. The optical depth of the blob is of the order of blob tm/τ0 Lblob(tm) τ0 mTi MTi,0 = NTi,0 mTi e , (3) ≥ Edecay Rblob τγ dr κγ ρblob , (8) 44 ∼ 0 where mTi is the mass of a Ti atom. Inserting num- Z bers, including Lblob(tm) = Lcs , where Rblob and ρblob are the blob radius and aver- blob 6 age density, respectively. For radiation particles i, the MTi,0 (4.5 10.5) 10− M , (4) ≥ − × 1 opacity is κi = Ye mb− σi, with Ye being the electron 44 which accounts for 3%–7% of the total mass of Ti faction, mb the average baryon mass, and σi the re- ∼ 4 ejected in SN 1987A ( 1.5 10− M ) as deduced from action cross section. The cross section is given by the ∼ × X-ray observations and light-curve analysis (Boggs et al. Klein-Nishina formula for standard Compton scattering 2015; Jerkstrand et al. 2011). of the 68 keV, 78 keV, and 1157 keV γ photons, which Mixing processes in 3D SN explosion models yield an are σ68 σ78 = 0.8σT and σ1157 = 0.3σT, respectively, 44 ≈ 25 2 initial Ti mass fraction relative to the freshly pro- where σ 6.65 10− cm is the Thomson scatter- T ≈ × duced 56Ni mass fraction that varies within the inter- ing cross section of electromagnetic radiation with elec- val 0.001–0.004 on large scales (Wongwathanarat et al. trons. The 44Sc decay positrons possess a broad energy 2017). Therefore, one estimates the initial mass ratio of spectrum with a peak around 600 keV and a maximum radioactive titanium and nickel in the blob to be energy of 1474 keV. Their typical annihilation cross sec- tion with electrons is of the order σT, but the γ-rays 3 blob blob 3 . 10− M /M 4 10− , (5) . Ti,0 Ni,0 . × emitted by the annihilation have a Klein-Nishina cross section that is smaller, namely in the range of (0.29– and the corresponding mass of 56Ni in the blob to be ∼ 0.43)σT. A safe upper limit of the γ opacity of the blob blob blob blob is therefore obtained with σγ 0.8σT, and is given by 250 MTi,0 . MNi,0 . 1000 MTi,0 . (6) ≤ 2 With Eq. (4), one thus derives Ye cm κ 0.16 . (9) γ . 0.5 g 3 blob (1.1 2.6) 10− M . MNi,0 − × 3 . (4.5 10.5) 10− M , (7) The 5σ emission knot, the blob, in Cigan et al.(2019) − × has an angular diameter of θ 20 mas. With a dis- blob ≈ which means that the blob should carry between 1.6% tance of 51.4 kpc this converts to a blob radius and 14.6% of the total 56Ni mass (0.072 M ) produced by SN 1987A. θblob R 7.66 1015 cm . (10) The lower limit of this range is in the ballpark of the blob ≈ × 20 mas 3 radioactive nickel mass, about 10− M , that has been estimated to be connected to a feature of the Hα line The average blob density can be expressed as profile during the so-called Bochum event (Utrobin et al. 3 M 3 M blob 1995). Wang et al.(2002) discussed an association of the ρ blob = Ni,0 , (11) 56 blob 3 3 blob Ni clumps of the Bochum event with the distribution ≈ 4π Rblob 4π RblobXNi,0 of 44Ti that powers the late-time ejecta. They argued blob 56 that both share the same orientation with respect to the where XNi,0 is the initial mass fraction of Ni in the line of sight, namely a location north of it and away from clump of ejecta forming the blob. Using Eqs. (9)–(11) NS 1987A in SN 1987A 5 in Eq. (8) leads to 8 Spin- τγ κγ Rblob ρblob down ∼ age = blob M κ 1014 4 yrs 2 Ni,0 γ 10 G 10 . 1.3 10− 2 1 × 0.01 M 0.16 cm g− ! Crab 10 6 2 12 12 -1 yrs − G 10 blob 1 θblob ) (X ) , (12) 1 erg s Ni,0 − 40 × 20 mas 10 14 108 yrs blob (s s -1 where the normalization of MNi,0 is guided by the upper ˙ P limit of the mass interval given in Eq. (7) and yields an erg s 10 10 yrs 16 10 35 J1210-5226 10 10 G upper limit of τγ . For efficient γ thermalization it is necessary that τ 1, which requires, at least, that Log J0821-4300 J1852+0040 γ & 18 1.5-3.5×10 -1 M blob blob 2 Ni,0 erg s XNi,0 . 1.3 10− , (13) 30 × 0.01 M 20 10 ! 108 G if all other factors in Eq. (12) are chosen to be unity. 22 blob blob 1 With X = M M − , this equation leads to a 3 2 1 0 1 Ni,0 Ni,0 blob lower limit of the blob mass: Log10 P (s)
Mblob & 0.77 M . (14) Figure 1. Measured period derivatives P˙ vs periods P for 2,256 pulsars from the ATNF Pulsar Catalogue (Manch- Therefore, the condition τγ & 1 can only be fulfilled with an implausibly large mass of the blob, nearly as much ester et al. 2005, and its on-line version at https://www. as the entire mass of oxygen ejected in SN 1987A. atnf.csiro.au/research/pulsar/psrcat). Red dots show pul- sars with a SNR association. Shown are lines of constant For all reasonable values of Mblob, we therefore con- spin-down timescale τsd ≡ P/(2P˙ ) (dashed blue), inferred clude that τ 1 and thermalization of the γ’s released 19 ˙ 1/2 γ polar magnetic field strength Bp ≡ 3.2 × 10 (P P ) G by 44Ti decay in the blob is highly inefficient. Under (Lyne & Graham-Smith 2012) (violet dots), and spin-down such circumstances, with only partial thermalization of power W ≡ |IΩΩ˙ | (red dots), where Ω ≡ 2π/P and an as- the γ photons being possible, a 44Ti mass much higher sumed moment of inertia I = 1.5 × 1045 g cm2 (Haensel et al. 2007). The dark red lines delimit the inferred range than our estimate of Eq. (4) would be needed to explain 35 1 Lcs ∼ (40–90) L = (1.5–3.5) × 10 erg s− corresponding L . Again, this is implausible and it would require ex- cs to NS 1987A; the adjacent thin dashed red line shows the treme fine-tuning if several tens of percent of the ejected 35 1 upper limit of 138 L = 5.2 × 10 erg s− for the luminosity 44 Ti, or even all of it, were concentrated in a single blob of any compact source in the SNR (Alp et al. 2018). The instead of being widely distributed in association with locations of the Crab pulsar and the only three CCOs with the expelled iron-group elements (for corresponding 3D measured P and P˙ , PSR J1852+0040 in Kes 79 (Halpern & explosion simulations in comparison to the gaseous rem- Gotthelf 2010), PSR 0821-4300 in Puppis A, and PSR J1210- nant of Cas A, see Wongwathanarat et al. 2017). For 5226 in PKS 1209-52 (Gotthelf et al. 2013), are indicated. all these reasons 44Ti-decay heating cannot provide a convincing explanation of the entire blob luminosity. 3. THE SPIN-DOWN POWER OF YOUNG Nevertheless, a small fraction of the blob luminosity PULSARS AND CENTRAL COMPACT of (40–90) L might still come from 44Ti-decay heat- OBJECTS ing (Cigan et al. 2019) or irradiation by external sources. An obvious explanation for the source of energy pow- If such mechanisms were responsible for heating the blob ering the dust blob is the pulsar’s spin-down power to the 22 K of its surroundings, they would produce ˙ ∼ W IΩΩ , where I is the moment of inertia and Ω only a fraction of (22/33)4 20 % of the blob luminosity ≡ | | ˙ ≈ is the rotational frequency. Fig.1 shows the P -P di- (blackbody emission is assumed for this estimate). An agram containing known pulsars to put this hypoth- over-dense cloud of dust, which is not implausible for esis in context. While the needed luminosity Lcs is the blob due to the strong gravitational attraction of an within the range of values for young pulsars associ- enclosed compact object, might be heated by radioactiv- ated with a SNR, that range is very large, extending ity even somewhat above the 22 K of the surroundings. 31 1 from 2 10 erg s− in the case of PSR J1210-5226 For a 25 K radioactive “floor value” of the blob tem- ∼ × 39 1 up to 10 erg s− for the Crab pulsar. Moreover, 44 ' perature, for example, the Ti-decay would account for most young pulsars with “canonical” polar magnetic (25/33)4 33 % of the blob luminosity. ≈ 6 Page et al.
12 13 field strengths Bp 10 − G have W Lcs , while physical conditions that control a young neutron star’s ∼ 14 magnetar-type pulsars with Bp > 10 G exhibit values luminosity around this age. substantially below Lcs . The three CCOs with mea- Subsequent to its birth in a core-collapse SN, the sured P and P˙ , which must also be young pulsars and compact remnant enters the proto-neutron star phase are displayed in Fig.1, also have W L . in which hot and dense matter temporarily traps neu- cs The large possible range for W suggests fine-tuning is trinos (Burrows & Lattimer 1986). This phase lasts a involved if a pulsar is embedded in the blob and power- minute or so due to eventual energy and lepton number ing its excess luminosity. It must lie between the double loss from neutrino leakage after which matter reaches its red lines in Fig.1. And it may be even more unlikely beta-equilibrated chemical composition. After this, the that the pulsar is located outside the blob, since in this neo-neutron star (Beznogov et al. 2020) has a rapidly case W must be larger than Lcs and the pulsar must lie decreasing Lth for about a year until reaching a plateau 35 1 above the double red lines in Fig.1. But there is an of (0.5–5) 10 erg s− (12–120 L ) that lasts for a few × observed (Alp et al. 2018) upper limit to the power of decades. This range agrees well with Lcs . a hidden pulsar in the 1987A SNR that is only 138L , In the neo-neutron star phase, the surface temper- i.e., (1.5–3.5) L , suggesting the pulsar is located in- ature T is initially controlled mainly by the matter ' cs s credibly close to the blob. directly underneath it. With progressing time, how- The emerging family of CCOs comprises about a ever, conduction allows deeper and deeper layers to play dozen objects, all young by definition as they are a role. The early cooling, after the first day, is en- found in SNRs (see the on-line listing at http://www. tirely due to neutrino emission by plasmon decay (γ∗ → iasf-milano.inaf.it/∼deluca/cco/main.htm), and, except νν¯) and results in a virtually universal temperature- for the three mentioned above, generally presenting density (T -ρ) profile in the outer layers. When the lo- no evidence for a significant magnetic field (see, e.g., cal T (ρ) has dropped below the local plasma temper- Gotthelf & Halpern 2008; de Luca 2008, 2017 for re- ature TP (ρ), the plasmon density is exponentially re- views). The ATNF Pulsar Catalogue (Manchester et al. duced, and, consequently, plasmon decay ceases as well 2005) lists 32 pulsars with spin-down timescales τsd (Adams et al. 1963). Hence, the T (ρ) profile follows ˙ ≡ P/(2P ) < 10, 000 yrs. Of these, 21 are “regular” pul- from the electron plasma temperature TP, e ~ΩP, e/kB 14 2 2 ≡ sars with Bp < 10 G, and 11 are “magnetars”, i.e., where ΩP, e 4πe ne/me∗ and ne is the net electron den- 14 ≡ having Bp > 10 G. Since 12 CCOs are known, this sity. The effective mass of ultra-relativistic electrons is indicates that a young neutron star has a significant me∗ = ~kF (e)/c, where the electron Fermi momentum 12/(21 + 11 + 12) 27% chance of being observed as a k (e) (3π2n )1/3. Thus, T (ρ) is directly determined ' F ≡ e CCO. (Of course, the actual probability of being a CCO by the ne(ρ) profile. The standard assumption is that, could be different because of selection biases.) Indeed, being formed by cooling from very high temperatures, NS 1987A had been proposed to be a CCO before the the crust of a neo-neutron star consists of “catalyzed class was even given a name (Muslimov & Page 1995; matter” in which the optimum nucleus minimizes the Page et al. 1998; Geppert et al. 1999). energy density (Baym et al. 1971) at each baryon num- The reasonably large probability of being observed as ber density. At low-enough densities, this implies the 56 a CCO and the fine-tuning in P and Bp required to ex- presence of Fe. As a consequence, ne(ρ) is a univer- plain the blob’s luminosity whether or not it encloses the sal function, and, from the above discussion and before neutron star imply the energy source powering the dust thermal diffusion takes over in a very young star, so is blob is most likely of a different nature than spin-down. T (ρ) (see Beznogov et al. 2020 for more details). We therefore favor the idea that NS 1987A is a CCO. Once plasmon decay shuts off, crustal neutrino emis- In the next section, it is shown that Lth, except under sion continues from electron-ion bremsstrahlung, which certain circumstances, is naturally the same magnitude is, however, relatively inefficient (Gnedin et al. 2001). as Lcs. Hence, subsequent cooling proceeds very slowly with further evolution mainly determined by heat transport. 4. THE THERMAL LUMINOSITY OF A 30 YEAR The core maintains more powerful neutrino emission and OLD NEUTRON STAR is therefore colder than the crust, resulting in the crustal thermal energy diffusing inward, which leads to the de- The early thermal evolution of a neutron star has been crease of Ts on the thermal diffusion time scale. specifically studied in several works (see, e.g., Lattimer Thermal diffusion through a layer of width l acts on et al. 1994; Gnedin et al. 2001; Shternin & Yakovlev a timescale (Brown et al. 1988) 2008) and in most cases the luminosity at age 30 years is within a factor 0.1–1 of L . Here we examine the τ l2 C /K , (15) cs th ≈ V NS 1987A in SN 1987A 7 where CV is the specific heat per unit volume and K To quantify the above considerations, we have per- is the thermal conductivity of the layer. The width l formed extensive numerical simulations using the code of any crustal layer scales with the total crust thickness NSCool (Page 1989, 2016) that we previously used in ∆R, which is primarily a function of the star’s mass M Page et al.(2004, 2009, 2011). Our results are pre- and radius R (Lattimer et al. 1994). The specific heat sented in Fig.2. Each panel corresponds to different poses a major uncertainty due to the presence of dripped assumptions about the core’s physics (i.e., the dense neutrons, which contribute negligibly if they are in a su- matter EOS and core superfluidity) or the stellar mass, perfluid state (i.e., if the local T (ρ) < Tc(ρ), the local whereas the various cooling curves in each panel corre- superfluid transition critical temperature), whereas they spond to different assumptions about the crust’s physics will largely dominate if they are in the normal state. Un- (superfluidity) and the envelope’s composition. certainty in knowledge of Tc(ρ) directly translates into Although the early evolution is dominated by the crust a large uncertainty in CV and has a significant impact and envelope, the neutron star core also plays a role on the thermal relaxation of the crust as shown below. through two different effects. First, the core’s EOS de- In contrast, the thermal conductivity is dominated by termines the stellar mass M, radius R and crust thick- well-understood electron scattering so that K(ρ) in it- ness ∆R. The latter determines the crustal length scales self has little uncertainty once the chemical composition l entering Eq. (15). Simulations are restricted to EOSs of the medium is determined3. Catalyzed matter in the satisfying radius constraints described in AppendixA. outer layers implies the presence of 56Fe up to densi- From the four EOSs we consider in AppendixA, only 7 3 ties 10 g cm− and increasingly more neutron-rich results for MS-A1 and APR that have the largest and ∼ nuclei as density increases. However, an important as- smallest R (and ∆R), respectively, are displayed. The pect of the problem is the possibility that some amount other two EOSs lead to intermediate results. of light elements (H/He/C/O) may have been accreted Second, neutrino cooling of the core determines the after the crust originally equilibrates. These light ele- final crust temperature once its thermal relaxation is ments have smaller Z and electron-ion scattering rates, completed. Core cooling generally follows one of two and therefore larger thermal conductivities, than heavier paths, “standard” (modified Urca or MU) or “enhanced” elements, which lead to a smaller thermal gradient in the (direct Urca or DU) cooling (see AppendixB). In this outer layers and higher surface temperature Ts. Thicker section, we only consider standard cooling, and, in Sec- light-element layers correspond to larger Ts (Potekhin tion6, we consider enhanced cooling which results in et al. 1997). We will call this outer layer the envelope,4 colder cores and possibly shorter crust relaxation time where a strong temperature gradient toward the surface since CV decreases and K increases (Eq. (15)). If re- may be present and the chemical composition may have stricted to standard or slow neutrino cooling, i.e., ad- been altered from that of the original catalyzed matter. hering to the Minimal Cooling paradigm, core neutrino 3 We denote by ρ = 10L g cm− the maximum density emission is mostly controlled by nucleon pairing, either L reached by light elements, there being catalyzed matter p-superconductivity or n-superfluidity. The main effects at higher densities. of pairing are to suppress the dominant MU process and/or trigger the efficient Cooper pair breaking and formation processes (see, e.g., Page et al. 2014). 3 Electron scattering by impurities and lattice defaults is the major source of uncertainty but only at temperatures much lower than Panels (a) and (b) of Fig.2 explore the consequences those present in the crust of a young neutron star. The reader of variations in ∆R due to the core’s EOS under the is referred to Page & Reddy(2012) for a review on the thermal assumption of no core pairing. EOS-induced changes properties of a neutron star’s crust. are seen to be small. Panels (b) and (c) explore the 4 Traditionally, see, e.g., Page et al.(2006), the envelope encom- suppression of the MU emission by p-superconductivity, passes the outer layers at densities below a boundary value ρb. These layers are treated separately for numerical convenience and resulting in higher luminosities after crust relaxation provide the outer boundary condition as a “Tb-Ts”relationship and the consequent smaller temporal extension of the (Gudmundsson et al. 1983) between the boundary temperature crust relaxation phase. The assumed p-pairing model, Tb and the surface Ts. Our present definition of the envelope is thus an extension to include all layers with possible non-catalyzed CCDK from Chen et al.(1993a), has extensive p- 11 3 chemical composition extending to a density ρc = 10 g cm− superconductivity covering almost the entire stellar core, so that ρ ρc always. We thus split the outer layers in an L thus maximizing the suppression of neutrino cooling. Fi- outer envelope≤ at densities ρ ρ included into the outer bound- ≤ b ary condition, and an inner envelope with ρb ρ ρc possibly nally, panels (d) through (f) display the effect of adding ≤ ≤ 3 contaminated by light elements (see Beznogov et al. 2020 for a core n- P2 superfluidity with a small gap as needed similar extension). (Page et al. 2011; Shternin et al. 2011) to explain the observed Lth and dLth /dt of the neutron star (Cas A) in 8 Page et al. (a) (b) (c) 36.0 3 3 3 R =13.1 km n P2: none R =11.4 km n P2: none R =11.4 km n P2: none
1 1 1 R =1.23 km p S0: none R =0.83 km p S0: none R =0.83 km p S0: CCDK
11 11 11 35.5 NS 1987A NS 1987A NS 1987A ]
1 10 10 10 9 9 9 8 8
[erg s 8 35.0 th
L H H H 10 log 34.5 F M =1.4 M M =1.4 M M =1.4 M N EOS : MS-A1 EOS : APR EOS : APR 34.0 (d) (e) (f) 36.0 3 3 3 R =11.5 km n P2: CasA R =11.4 km n P2: CasA R =11.3 km n P2: CasA
1 1 1 R =1.04 km p S0: CCDK R =0.83 km p S0: CCDK R =0.66 km p S0: CCDK 11 11 35.5 11 NS 1987A NS 1987A NS 1987A 10 ]
1 10