Arxiv:Astro-Ph/0002114V2 14 Mar 2000

Home , Thermal runaway

arXiv:astro-ph/0002114v2 14 Mar 2000 eetydsoee -oeisaiiy(e neso & Andersson the (see on instability relies r-mode model gradi- discovered gravitational-wave temperature second recently generated The accretion a star ents. neutron of to deep due because the are Cutler crust Bild- in waves induced Ushomirsky, by gravitational deformation by the (proposed quadrupole refined (2000)), LMXB Bildsten recently the & and for (1998) In model that sten such torque. idea accretion first the on the the based balances being radiation mod- both Zhang gravitational two different, remaining & rather the are White for els while model to spin-up, magnetosphere (due induced standard the sur- accretion model on this first based is explain The (1997)) to proposed result. Three been prising 2000). have Klis, der models (van perhaps Hz different rapidly, 260-590 spin range LMXB narrow sug- in the in stars models Various neutron the see theorists. pro- that for stars, gest observations challenge neutron RXTE further recent a recent of rotating vide The review of a periods (1998)). modelling (for shorter Stergioulas the PSR1937+21 at reason in of pulsars the ms progress radio know 1.56 of not the lack do than apparent still the we for theoretical example, and/or For observations further modeling. by answered millisecond be to of accretion. to formation due the support spin-up for providing via (MSP) (2000)), model pulsars Klis re- standard der a the van (for for see stars neutron review, these spinning cent that rapidly suggest Low- contain strongly dozen systems a (LMXB) than Binaries neu- quasiperiodic more X-ray of in of Mass frequencies understanding observations kHz our Detailed at in phenomena era physics. new star a tron heralded 1995 in rpittpstuigL using typeset Preprint 2018 22, October version Draft ept hs dacssvrldffiutqetosremain questions difficult several advances these Despite (RXTE) Explorer Timing X-ray Rossi the of launch The siae ffcstepeitdiiilsi eid fhtyugneu We young hot star. of neutron periods cold headings: spin phase. relatively Subject initial spin-down and predicted driven pr old r-mode the the an brief affects to the in estimates due during crust dissipation radiated the of are r-m strength of that of the base mechan instability of damping the estimate the viscous at revised i) as a released are: iii) heat model Low-Mass to and due the the runaway in in thermal ingredients ii) stars main neutron The both pulsars. of observations current with ilmil ecnndt arwrneo spin-frequencies; of range narrow a to confined be mainly will epeetasml pneouinmdlta rdcsta ail ro rapidly that predicts that model evolution spin simple a present We eateto ahmtc,Uiest fSuhmtn Sou Southampton, of University Mathematics, of Department A 1. T E tl mltajv 04/03/99 v. emulateapj style X INTRODUCTION eateto hsc,AitteUiest fThessaloni of University Aristotle Physics, of Department -OERNWYADRPDYRTTN ETO STARS NEUTRON ROTATING RAPIDLY AND RUNAWAY R-MODE es atr–gaiain–sas eto tr:rtto s — rotation stars: – neutron stars: – gravitation – matter dense otsD oktsadNklo Stergioulas Nikolaos and Kokkotas D. Kostas oktsatoat.r [email protected] [email protected], amtsstna.k [email protected] [email protected], isAdrsnadDvdInJones Ian David and Andersson Nils rf eso coe 2 2018 22, October version Draft Acpe ac.1,2000) 10, March. (Accepted ABSTRACT 1 ttebs ftesldcuti eaieycl neutron cold relatively a in crust arises solid that the layer of boundary base dissipa- viscous the on the a at have that with would estimated associated crust They solid tion a oscillations. of r-mode presence the relevant the investigated that not (2000) effect are Ushomirsky the & r-modes Bildsten LMXB. the the that for of conclusion fraction the tiny phasize a for radiate gravitational lifetime. only detectable systems will of the source they a since as rota- waves, r-modes low LMXB the rather out galactic rules a would in conclusion to this r-mode Essentially, down the star rate. months tion neutron few accreting a an a spin in triggers and mode-heating runaway the thermal (via increases, weaker temperature up gets rel- viscosity the the star shear as at the the r-mode Since heat the temperatures). counteracts will evant that generat- r-modes viscosity to shear the the addition system, momentum in the angular the dissipate since, (1998), from that viable Spruit waves gravitational and not ing (1999) is as Levin idea However, by prime original out waves. be would pointed gravitational LMXB was Should detectable in for stars (1984). spin- Pa- neutron sources Wagoner a by the and r-modes suggestions reach happen, (1978) of this the would Pringle vein level & which the paloizou in significant in much a very was star equilibrium, to it accreting Originally, excited Stergioulas an were (1998)). & that Bildsten Kokkotas also an- thought Andersson, see first by (but was (1999) possibility detail bal- This in that torque. alyzed agent accretion the the provide ances r-modes unstable accreted from waves the dissipate star. neutron to the review) from momentum a angular for (2000) Kokkotas te eetrslswud(tfis ih)se oem- to seem sight) first (at would results recent Other gravitational that idea the reexamine we Letter this In hmtnS1 B,Uie Kingdom United 1BJ, SO17 thampton i hsaoii506 Greece 54006, Thessaloniki ki, P eas reytuho o h new the how on touch briefly also We -a iaisadmlieodradio millisecond and Binaries X-ray 1 = rnstars. tron sscutrc h -oegrowth, r-mode the counteract isms sneo icu onaylayer boundary viscous a of esence . 5 ics h rvttoa waves gravitational the discuss dsaoeaciia pnrate, spin critical a above odes − aigaceignurnstars neutron accreting tating s hsi nagreement in is This ms. 5 as oscillations tars: 2 R-mode runaway and rapidly rotating neutron stars star would greatly exceed that of the standard shear vis- From the fundamental relations, e.g. the formula for the cosity. Thus, Bildsten and Ushomirsky concluded that the gravitational-wave energy radiated via the current multi- r-mode instability would only be relevant for very high ro- poles, one can see that the gravitational-wave timescale tation rates, and could therefore not play a role in the should scale with M, R and P in the way manifested in LMXB. (1). Thus, we believe that Bildsten and Ushomirsky un- We have reassessed the effect of the viscous boundary derestimate the strength of radiation reaction significantly, layer (correcting an erring factor in the estimates of Bild- which motivates us to reassess the relevance of the viscous sten & Ushomirsky (2000)). Our new estimates show that boundary layer. the presence of the crust is important, but that the insta- We should, of course, emphasize at this point that our bility operates at significantly lower spin rates than sug- current understanding of the r-mode instability is based gested by Bildsten and Ushomirsky. Once we combine our on crude estimates of the various timescales. In order to estimates with the thermal runaway (now due to heating understand the role of the instability in an astrophysical caused mainly by the presence of the viscous boundary context we must improve our modelling of many aspects layer), that results as the star is spun up to the point of neutron star physics such as the effect of general rela- at which the instability sets in, we arrive at a model for tivity on the r-modes, cooling rates, viscosity coefficients, the spin-evolution of accreting neutron stars. Remarkably, magnetic fields, potential superfluidity, the formation of a this simple model agrees well with existing observations of solid crust etcetera (see Andersson & Kokkotas (2000) for rapidly rotating neutron stars, covering both the LMXB a description of recent progress in these various directions). and MSP populations. In the following we will mainly consider uniform density stars, i.e. use the gravitational-wave timescale given 2. DISSIPATION DUE TO A VISCOUS BOUNDARY LAYER by (1). In estimating the dissipation timescale tvbl due The r-mode instability follows after a tug of war be- to the presence of the crust, we need to evaluate tvbl ≈ tween (mainly current multipole) gravitational radiation 2E/(dE/dt) where E is the mode-energy, and dE/dt fol- − that drives the mode and various dissipation mechanisms lows from an integral over the surface area at the crust- that counteract the fluid motion. In the simplest model, core boundary (assumed to be located at radius Rb), cf. the mode is dominated by shear viscosity at low tempera- equation (3) of Bildsten & Ushomirsky (2000). To evalu- tures while bulk viscosity may suppress the mode at high ate this integral we use the standard result for the shear temperatures. At intermediate temperatures, the r-mode viscosity in a normal fluid. To incorporate our uniform sets an upper limit on the neutron star spin rate. In an density model, we make the reasonable assumption that interesting recent paper, Bildsten & Ushomirsky (2000) the density of the star is constant ( M/R3) inside radius ∼ estimate the strength of dissipation due to the solid crust R. Then it falls off rapidly in such a way that the base of an old neutron star, and find that the presence of a of the crust is located at a radius only slightly larger than boundary layer at the base of the crust leads to a very R. Hence, it makes sense to use Rb R. If we neglect the ≈ strong damping of the r-modes. small mass located outside radius R we can then immedi- While we agree with the main idea and the various as- ately compare the result for the viscous boundary layer to sumptions made by Bildsten and Ushomirsky, we would the timescales used by Andersson, Kokkotas & Stergioulas like to point out one important difference between their (1999). In the end, our estimate for the dissipation due to results and ones used previously in the literature. Their the presence of the viscous boundary layer is assumed timescale for gravitational radiation reaction dif- 2 1/2 fers significantly from, for example, the uniform density M 10 km T P tvbl 200 s . result derived by Kokkotas & Stergioulas (1999) (and sub- ≈ 1.4M⊙ R 108 K 1 ms sequently used by several authors, see Andersson & Kokko- (2) tas (2000)). This is surprising since the uniform density which is a factor of 2 larger than that of Bildsten and result, which can be written Ushomirsky. This difference arises simply because the mode-energy E is this factor larger for uniform density 4 6 1.4M⊙ 10 km P models. The star is assumed to have a uniform tempera- t 22 s , (1) gw ≈− M R 1 ms ture distribution, with core temperature T . In Figure 1 we show the instability window obtained (where the negative sign indicates that the mode is unsta- from our revised estimate. As is clear from this Figure, ble) has been shown to be close (within a factor of two) the presence of a viscous boundary layer in an old, rela- to the results for n = 1 polytropes. M, R, and P repre- tively cold neutron star is, indeed, important. However, sent the mass, radius and spin period of the star, respec- Bildsten and Ushomirsky’s conclusion that the r-mode in- tively. In contrast, Bildsten & Ushomirsky (2000) use the stability is irrelevant for the LMXB cannot be drawn from n = 1 polytrope result and argue that it corresponds to Figure 1. On the contrary, the Figure suggests that the tgw 146 s for a canonical neutron star rotating with instability may well be limiting the rotation of these sys- a period≈ − of 1 ms, i.e. assume that radiation reaction is tems. almost one order of magnitude weaker than in (1). This difference occurs because Bildsten and Ushomirsky have 3. THERMAL RUNAWAY IN RAPIDLY SPINNING NEUTRON STARS only rescaled the fiducial rotation frequency Ω0 = √πGρ¯ (whereρ ¯ represents the average density) in terms of which The fact that our revised instability curve for r-modes the n = 1 polytrope results of Owen et al. (1998) were ex- damped by dissipation in a viscous boundary layer agrees 6 pressed (tgw 3.26(Ω0/Ω) for a specific polytropic stel- well with the fastest observed neutron star spin frequen- lar model). Unfortunately,≈− this procedure is not correct. cies, cf. Figure 1, motivates us to speculate further on the Andersson et al. 3 relevance of the instability. We want to model how the main prediction: Once an accreting neutron star has been potential presence of an unstable r-mode affects the spin- spun up beyond (say) 5 ms it must remain in the rather evolution of rapidly spinning, accreting neutron stars. To narrow range of periods 1.5 5 ms until it has stopped do this we use the phenomenological two-parameter model accreting and magnetic dipole− braking eventually slows it devised by Owen et al. (1998), which is centered on evo- down. Since a given star can go through several r-mode lution equations for the rotation frequency Ω and the (di- cycles before accretion is halted one would expect most mensionless) r-mode amplitude α. Complete details of our neutron stars in LMXB and the MSP to be found in the particular version of this model will be given elsewhere. predicted range of rotation rates. As is clear from Figure 1, At the qualitative level, our results are not surprising. this prediction agrees well with the range of rotation peri- Accreting stars in the LMXB are expected to have core ods inferred from observed kHz quasiperiodic oscillations temperatures in the range 1 4 108 K (Brown & Bild- in LMXB. The observed range shown in Figure 1 corre- sten, 1998). For such temperatures− × the dissipation due to sponds to rotation frequencies in the range 260-590 Hz the viscous boundary layer gets weaker as the tempera- (cf. van der Klis (2000)). Our model also agrees with ture increases. Consequently, the situation here is essen- the observed data for MSP, which are mainly found in the tially identical to that considered by Levin (1999) (see also range 1.56 6 ms, see Figure 2. In other words, our pro- Spruit (1998) and Bildsten & Ushomirsky (2000)). After posed model− is in agreement with current observed data accreting and spinning up for something like 107 years, for rapidly rotating neutron stars. the star reaches the period at which the r-mode instabil- Finally, it is worthwhile discussing briefly the detectability sets in. For our particular estimates this corresponds ity of the gravitational waves that are radiated during the to a period of 1.5 ms (at a core temperature of 108 K). relatively short time when the r-mode is saturated and the It is notable that this value is close to the 1.56 ms pe- star spins down. As was argued by Levin (1999), the fact riod of PSR1937+21. Once the r-mode becomes unstable that the r-mode is active only for a small fraction of the (point A in Figure 1), viscous heating (now mainly due to lifetime of the system (something like 1 month out of the the energy released in the viscous boundary layer) rapidly 107 years it takes to complete one full cycle) means that heats the star up to a few times 109 K. The r-mode am- even though these sources would be supremely detectable plitude increases until it reaches a prescribed saturation from within our galaxy the event rate is far too low to make level (amplitude αs) at which unspecified nonlinear effects them relevant. However, it is interesting to note that the halt further growth (point B in Figure 1). Once the mode spin-evolution is rather similar to that of a hot young neu- has saturated, the neutron star rapidly spins down as ex- tron star once the r-mode has reached its saturation am- cess angular momentum is radiated as gravitational waves. plitude. This means that we can analyze the detectability When the star has spun down to the point where the mode of the emerging gravitational waves using the framework again becomes stable (point C in Figure 1), the amplitude of Owen et al. (1998). We then find that these events can starts to decay and the mode plays no further role in the be observed from rather distant galaxies. For a source in spin evolution of the star (point D in Figure 1) unless the the Virgo cluster (assumed to be at a distance of 15 Mpc) star is again spun up to the instability limit. Two exam- these gravitational waves could be detected with a signal ples of such r-mode cycles (corresponding to αs =0.1 and to noise ratio of a few using LIGO II. However, even at the 1, respectively) are shown in Figure 1. distance of the Virgo cluster these events would be quite The real surprise here concerns the quantitative pre- rare. By combining a birth rate for LMXB of 7 10−6 dictions of our model. As already mentioned, the model per year per galaxy with the fact that the the volume× of suggests that an accreting star will not spin up beyond space out to the Virgo cluster contains 103 galaxies, and 1.5 ms. This value obviously depends on the chosen stel- the possibility that each LMXB passes∼ through (say) four lar model, but it is independent of the r-mode saturation r-mode cycles during its lifetime we deduce that one can amplitude and only weakly dependent on the accretion only hope to see a few events per century in Virgo. In rate (through a slight change in core temperature). In order to see several events per year the detector must be fact, the accretion rate only affects the time it takes the sensitive enough to detect these gravitational waves from star to complete one full cycle. As soon as the mode be- (say) 150 Mpc. This would require a more advanced de- comes unstable the spin-evolution is dominated by gravi- tector configuration such as a narrow-banded LIGO II. We tational radiation and viscous heating. Once the star has will discuss this issue in more detail elsewhere. gone through the brief phase when the r-mode is active it 4. has spun down to a period in the range 2.8-4.8 ms (cor- ADDITIONAL REMARKS responding to 0.01 αs 1). Based on these results we Before concluding our discussion, we recall that the ini- ≤ ≤ propose the following spin-evolution scenario: An accret- tial excitement over the r-mode instability was related to ing neutron star will never spin up beyond (say) 1.5 ms. the fact that it provided an explanation for the relatively Once it has reached this level the r-mode instability sets slow inferred spin rates for young pulsars. In view of this, in and spins the star down to a several ms period. At this it is natural to digress somewhat and discuss how the pic- point the mode is again stable and continued accretion ture of the r-mode instability in hot, newly born neutron may resume to spin the star up. Since the star must ac- stars is affected by the possible formation of a solid crust. crete roughly 0.1M⊙ to reach the instability point, and the Hence, we consider the evolution of a neutron star just af- LMXB companions have masses in the range 0.1 0.4M⊙, ter its birth in a supernova explosion. At a first glance, we it can pass through several “r-mode cycles” during− its life- might expect to model its r-mode amplitude in the stan- time. dard way, cf. Owen et al. (1998), using the normal (crust- Let us confront this simple model with current obser- free) fluid viscous damping times for stellar temperatures vations. To do this we note that our model leads to one above the melting temperature of the crust (Tm), and the 4 R-mode runaway and rapidly rotating neutron stars

12 0.8 10 0.7 0.6 8 A B 0.5 6

PK/P 0.4 LMXB 4 0.3 C α =0.1 s number of recycled pulsars 2 0.2 α D s=1 0.1 0 0 2 4 6 8 10 12 14 16 0.0 P (ms) 7.0 7.5 8.0 8.5 9.0 9.5 10.0 10.5 r-mode cycle

log10 (T/1K) Fig. 2.— We compare the “r-mode cycle” predicted for ac- Fig. 1.— The r-mode instability window relevant for old neutron creting neutron stars to the observed rapidly spinning pulsars. stars. We show results for the simplest (crust-free) model, where The observed data is taken from the Princeton pulsar database gravitational radiation reaction is balanced by standard shear vis- (http://pulsar.princeton.edu/pulsar/catalog.shtml) as well as a sam- cosity at low temperatures and bulk viscosity at high temperatures ple of MSP recently discovered in 47Tuc (Camilo et al., 1999). (thin solid line). Also shown are Bildsten and Ushomirsky’s estimate for a star with a crust (dashed line) and our improved estimate of this situation (thick solid line). The r-mode is potentially unstable in the region above each curve. The illustrated results correspond to a canonical neutron star for which mass shedding at the equator sets in at the Kepler period PK ≈ 0.8 ms. We illustrate two typical r-mode cycles (for mode saturation amplitudes αs = 0.1 and 1), resulting from thermo-gravitational runaway after the onset of instability. Once accretion has spun the star up to the critical period (along the indicated spin-up line (thick vertical line)) the r-mode becomes unstable and the star evolves along the path A-B-C-D. Af- ter a month or so, the mode is stable and the star will cool down until it again reaches the spin-up line. For comparison with obser- vational data, we indicate the possible range of spin-periods inferred from current LMXB data (shaded box) as well as the observed periods and estimated upper limits of the temperature (cf. Andersson, Kokkotas & Stergioulas (1999)) of some of the most rapidly spinning MSP (short arrows). viscous boundary layer damping time for temperatures be- Such a scenario was described by Bildsten & Ushomirsky low Tm. However, the situation is a little more complicated (2000), who noted that the r-mode instability would then than this. Recall that the latent heat (i.e. the Coulomb not spin the star down beyond a much higher frequency. 48 binding energy) of a typical crust is Elat 10 ergs, while Using our estimated timescales the resultant spin period the r-mode energy is E 2α2(1 ms/P∼ )2 1051 ergs. would be 3 5 ms in this scenario. Which scenario ap- m ≈ × − Provided that the time taken for the star to cool to Tm plies depends sensitively on the early cooling of the star, is suﬃciently long, the energy in the r-mode (which grows the crustal formation temperature and perhaps most im- exponentially on a timescale t 20 s) will exceed E , portantly the initial amplitude of the r-mode following the gw ∼ lat preventing the formation of the crust, even when T

REFERENCES

Andersson, N., Kokkotas, K. D., & Stergioulas, N. 1999, ApJ, 516, Owen, B. J., Lindblom, L., Cutler, C., Schutz, B. F., Vecchio, A., & 307 Andersson, N. 1998, Phys. Rev. D, 58, 084020 Andersson N., & Kokkotas K.D. 2000, to appear in Int. J. Mod. Papaloizou J. & Pringle J.E. 1978, MNRAS, 184, 501 Phys. D Spruit H., 1999, Astron. Astrop., 341, L1 Bildsten, L. 1998, ApJ, 501, L89 Stergioulas N. 1998, Living Reviews in Relativity, 1998-8 Bildsten, L. & Ushomirsky G. 2000, ApJ, 529, 33 http://www.livingreviews.org/Articles/Volume1/1998-8stergio/ Brown E.F. and Bildsten L. 1998, ApJ, 496, 915 Ushomirsky G., Cutler C. & Bildsten, L. 2000, Deformations of Camilo F., Lorimer D.R., Freire P., Lyne A.G.& Manchester R.N., accreting neutron star crusts and gravitational wave emission 1999 to appear in ApJ, [astro-ph/9911234] [astro-ph/0001136] Kokkotas, K. D., & Stergioulas N. , 1999, Astron. Astrop., 341, 110 van der Klis, M. 2000, Millisecond oscillations in X-ray binaries to Levin, Y. 1999, ApJ, 517, 328 appear in Ann. Rev. Astron. Astrop. [astro-ph/0001167] Marshall F.E., Gotthelf E.V., Zhang W., Middleditch J. and Wang Wagoner R.V. 1984, ApJ, 278, 345 Q.D. 1998, ApJ499, 179 White N.E. and Zhang W. 1997, ApJ, 490, L87