Vibration-Powered Radiation of Quaking Magnetar

Abstract

In juxtaposition with the standard model of rotation-powered pulsar, the model of vibration-powered magnetar undergoing quake-induced torsional Alfvén vibrations in its own ultrastrong magnetic field experiencing decay is considered. The presented line of argument suggests that the gradual decrease of frequencies (lengthening of periods) of long-periodic-pulsed radiation detected from a set of X-ray sources can be attributed to magnetic-field-decay-induced energy conversion from seismic vibrations to magnetodipole radiation of quaking magnetar.

1. Introduction

There is a common recognition today that the standard (lighthouse) model of inclined rotator, lying at the base of current understanding of radio pulsars, faces serious difficulties in explaining the long-periodic (2<𝑃<12 s) -pulsed radiation of soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs). The persistent X-ray luminosity of AXP/SGR sources (10³⁴<𝐿_𝑋<10³⁶ erg s⁻¹) is appreciably (10–100 times) larger than expected from a neutron star deriving radiation power from the energy of rotation with frequency of detected pulses. Such an understanding has come soon after the detection on March 5, 1979 of the first 0.2-second long gamma burst [1], which followed by a 200-second emission that showed a clear 8-second pulsation period [2], and association of this event to a supernova remnant known as N49 in the Large Magellanic Cloud [3]. This object is very young (only a few thousand years old), but the period of pulsating emission is typical of a much older neutron star. In works [4, 5] it has been proposed that discovered object, today designated SGR 0526-66, is a vibrating neutron star, that is, the detected for the first time long-periodic pulses owe their origin to the neutron star vibrations, rather than rotation as is the case with radio pulsars. During the following decades, the study of these objects has been guided by the idea [6, 7] that electromagnetic activity of magnetars, both AXP’s and SGR’s, is primarily determined by decay of ultrastrong magnetic field (10¹⁴<𝐵<10¹⁶ G) and that a highly intensive gamma bursts are manifestation of magnetar quakes [8–10].

In this paper we investigate in some detail the model of vibration-powered magnetar which is in line with the current treatment of quasiperiodic oscillations of outburst luminosity of soft gamma repeaters as being produced by Lorentz-force-driven torsional seismic vibrations triggered by quake. As an extension of this point of view, in this paper we focus on impact of the magnetic field decay on Alfvén vibrations and magnetodipole radiation generated by such vibrations. Before doing so, it seems appropriate to recall a seminal paper of Woltjer [11] who was first to observe that magnetic-flux-conserving core-collapse supernova can produce a neutron star with the above magnetic field intensity of typical magnetar. Based on this observation, Hoyle et al. [12] proposed that a strongly magnetized neutron star can generate magnetodipole radiation powered by energy of hydromagnetic, Alfvén, vibrations stored in the star after its birth in supernova event (see, also, [13]). Some peculiarities of this mechanism of vibration-powered radiation have been scrutinized in our recent work [14], devoted to radiative activity of pulsating magnetic white dwarfs, in which it was found that the necessary condition for the energy conversion from Alfvén vibrations into electromagnetic radiation is the decay of magnetic field. As was stressed, the magnetic field decay is one of the most conspicuous features distinguishing magnetars from normal rotation-powered pulsars. It seems not implausible, therefore, to expect that at least some of currently monitoring AXP/SGR-like sources are magnatars deriving power of pulsating magnetodipole radiation from energy of Alfvénic vibrations of highly conducting matter in the ultrastrong magnetic field experiencing decay.

In approaching Alfvén vibrations of neutron star in its own time-evolving magnetic field, we rely on the results of recent investigations [15–18] of both even parity poloidal and odd parity toroidal (according to Chandrasekhar terminology [19]) node-free Alfvén vibrations of magnetars in constant-in-time magnetic field. The extensive review of earlier investigations of standing-wave regime of Alfvénic stellar vibrations can be found in [20]. The spectral formula for discrete frequencies of both poloidal and toroidal 𝑎-modes in a neutron star with mass 𝑀, radius 𝑅, and magnetic field of typical magnetar, 𝐵₁₄=𝐵/10¹⁴ G, reads 𝜔_ℓ=𝜔_𝐴𝑠_ℓ,𝜔_𝐴=𝑣^𝐴𝑅=𝐵𝑅3𝑀,𝑀=43𝜋𝜌𝑅³,(1)𝜈_𝐴=𝜔^𝐴2𝜋=0.2055𝐵¹⁴𝑅^61/2,𝑀𝑀_⊙^−1/2,Hz,(2) where numerical factor 𝑠_ℓ is unique to each specific shape of magnetic field frozen in the neutron star of one and the same mass 𝑀 and radius 𝑅.

2. Alfvén Vibrations of Magnetar in Time-Evolving Magnetic Field

In above cited work it was shown that Lorentz-force-driven shear node-free vibrations of magnetar in its own magnetic field field 𝐁 can be properly described in terms of material displacements 𝐮 obeying equation of magneto-solid-mechanics 𝜌𝐮̈(𝐫,𝑡)=14𝜋[∇×[∇×[𝐮(𝐫,𝑡)×𝐁(𝐫,𝑡)]]]×𝐁(𝐫,𝑡),(3)𝐮̇(𝐫,𝑡)=[𝝎(𝐫,𝑡)×𝐫],𝝎(𝐫,𝑡)=𝐴_𝑡[∇𝜒(𝑟)]𝛼̇(𝑡).(4) The field 𝐮̇(𝐫,𝑡) is identical to that for torsion node-free vibrations restored by Hooke’s force of elastic stresses [18, 21] with 𝜒(𝐫)=𝐴_ℓ𝑓_ℓ(𝑟)𝑃_ℓ(cos𝜃), where 𝑓_ℓ(𝑟) is the nodeless function of distance from the star center and 𝑃_ℓ(cos𝜃) is Legendre polynomial of degree ℓ specifying the overtone of toroidal mode. In (4), the amplitude 𝛼(𝑡) is the basic dynamical variable describing time evolution of vibrations, which is different for each specific overtone; in what follows, we confine our analysis to solely one quadrupole overtone. The central to the subject of our study is the following representation of the time-evolving internal magnetic field: 𝐁(𝐫,𝑡)=𝐵(𝑡)𝐛(𝐫),(5) where 𝐵(𝑡) is the time-dependent intensity and 𝐛(𝐫) is dimensionless vector-function of the field distribution over the star volume. Scalar product of (1) with the separable form of material displacements 𝐮(𝐫,𝑡)=𝐚(𝐫)𝛼(𝑡)(6) followed by integration over the star volume leads to equation for amplitude 𝛼(𝑡) having the form of equation of oscillator with time-dependent spring constant ℳ𝛼̈(𝑡)+𝒦(𝑡)𝛼(𝑡)=0.(7) The total vibration energy and frequency are given by 𝐸_𝐴(𝑡)=ℳ𝛼̇²(𝑡)2+𝒦(𝐵(𝑡))𝛼²(𝑡)𝜔(𝑡)=2,𝒦(𝑡)ℳ=𝐵(𝑡)𝜅,𝜅=𝑅3𝑀𝑠.(8) It follows that 𝑑𝐸_𝐴(𝑡)𝑑𝑡=𝛼²(𝑡)2𝑑𝒦𝑑𝐵(𝐵)𝑑𝐵𝑑𝑡(𝑡).(9) This shows that the variation in time of magnetic field intensity in quaking magnetar causes variation in the vibration energy. In Section 3, we focus on conversion of energy of Lorentz-force-driven seismic vibrations of magnetar into the energy of magnetodipole radiation.

3. Vibration-Powered Radiation of Quaking Magnetar

The point of departure in the study of vibration energy-powered magnetodipole emission of the star (whose radiation power, 𝒫, is given by Larmor’s formula) is 𝑑𝐸_𝐴(𝑡)𝑑𝑡=−𝒫(𝑡),𝒫(𝑡)=32𝑐³𝛿𝝁̈²(𝑡).(10) We consider a model of a quaking neutron star whose torsional magnetomechanical oscillations are accompanied by fluctuations of total magnetic moment preserving its initial (in seismically quiescent state) direction: 𝝁=𝜇𝐧=constant. The total magnetic dipole moment should execute oscillations with frequency 𝜔(𝑡) equal to that for magnetomechanical vibrations of stellar matter, which are described by equation for 𝛼(𝑡). This means that 𝛿𝝁(𝑡) and 𝛼(𝑡) must obey equations of similar form, namely, 𝛼̈(𝑡)+𝜔𝛿₂𝝁̈(𝑡)+𝜔²(𝑡)𝛿𝝁(𝑡)=0,(𝑡)𝛼(𝑡)=0,𝜔²(𝑡)=𝐵²(𝑡)𝜅².(11) It is easy to see that (11) can be reconciled if 𝛿𝝁(𝑡)=𝝁𝛼(𝑡).(12) Given this, we arrive at the following law of magnetic field decay: 𝑑𝐵(𝑡)𝑑𝑡=−𝛾𝐵³(𝑡),𝛾=2𝜇²𝜅²3ℳ𝑐³𝐵(𝑡)=𝐵(0)=constant,(13)√1+𝑡/𝜏,𝜏⁻¹=2𝛾𝐵²(0).(14) The last equation shows that the lifetime of quake-induced vibrations in question substantially depends upon the intensity of initial (before quake) magnetic field 𝐵(0): the larger the 𝐵(0) the shorter the 𝜏. For neutron stars with one and the same mass 𝑀=1.4𝑀_⊙ and radius 𝑅=15km, and magnetic field of typical pulsar 𝐵(0)=10¹² G, we obtain 𝜏∼3×10⁷ years, whereas for magnetar with 𝐵(0)=10¹⁴ G, 𝜏∼7×10³ years.

The equation for vibration amplitude 𝜔𝛼̈(𝑡)+𝜔²(𝑡)𝛼(𝑡)=0,₂(𝑡)=𝜔²(0)1+𝑡/𝜏,𝜔(0)=𝜔^𝐴𝜅.(15) with help of substitution 𝑠=1+𝑡/𝜏 is transformed to 𝑠𝛼^(𝑠)+𝛽²𝛼(𝑠)=0,𝛽²=𝜔²(0)𝜏²=constant(16) permitting general solution [22]. The solution of this equation, obeying two conditions 𝛼(𝑡=0)=𝛼₀ and 𝛼(𝑡=𝜏)=0, can be represented in the form 𝛼(𝜂𝑡)==𝐽𝐶₁𝑠^1/2𝐽₁(𝑧(𝑡))−𝜂𝑌₁(𝑧(𝑡)),𝑧(𝑡)=2𝛽𝑠^1/2(𝑡),(𝑧(𝜏))𝑌₁(𝑧(𝜏)),𝐶=𝛼⁰𝐽₁(𝑧(0))−𝜂𝑌₁(𝑧(0))⁻¹,(17) where 𝐽₁(𝑧) and 𝑌₁(𝑧) are Bessel functions [23] and 𝛼²⁰=2𝐸_𝐴(0)𝑀𝜔²(0)=2𝐸_𝐴(0)𝐾(0),𝜔²(0)=𝐾(0)𝑀.(18) Here by 𝐸_𝐴(0), the average energy stored in torsional Alfvén vibrations of magnetar is understood. If all the detected energy 𝐸_burst of X-ray outburst goes in the quake-induced vibrations, 𝐸_burst=𝐸_𝐴, then the initial amplitude 𝛼₀ is determined unambiguously. The impact of magnetic field decay on frequency and amplitude of torsional Alfvén vibrations in quadrupole overtone is illustrated in Figure 1, where we plot 𝛼(𝑡) with pointed out parameters 𝛽 and 𝜂. The magnetic-field-decay-induced lengthening of period of pulsating radiation (equal to period of vibrations) is described by 𝑃(𝑡)=𝑃(0)1+𝑡𝜏^1/2,𝑃̇(𝑡)=12𝜏[1+𝑃(𝑡(/0𝜏))]^1/2𝜏=𝑃₂,(0)2𝑃(𝑡)𝑃̇(𝑡),𝑃(0)=𝜅𝐵2𝜋(0).(19) On comparing 𝜏 given by (14) and (19), one finds that interrelation between the equilibrium value of the total magnetic moment 𝜇 of a neutron star of mass 𝑀=1.4𝑀_⊙ and radius 𝑅=10km vibrating in quadrupole overtone of toroidal 𝑎-mode is given by 𝜇=𝐴(𝑀,𝑅)𝑃(𝑡)𝑃̇(𝑡),𝐴=3ℳ𝑐³8𝜋²𝜇=3.8×10₃₇,𝑃(𝑡)𝑃̇(𝑡)Gcm³.(20) For a sake of comparison, in the considered model of vibration-powered radiation, the equation of magnetic field evolution is obtained in a similar fashion as that for the angular velocity Ω(𝑡) in the standard model of rotation-powered neutron star which rests on 𝑑𝐸_𝑅𝑑𝑡=−32𝑐³𝐸𝛿𝝁̈²(𝑡),_𝑅(𝑡)=12𝐼Ω²(𝑡),𝐼=25𝑀𝑅²𝛿𝝁̈(𝑡𝝁)==𝜇𝛀𝐧(,𝑡)×𝜇=𝛀1(𝑡)×𝝁,,2𝐵𝑅³(21) which lead to Ω̇(𝑡)=−𝐾Ω³(𝑡),𝐾=2𝜇^2⟂3𝐼𝑐³Ω(𝑡)=Ω(0),𝜇_⟂=𝜇sin𝜃,√1+𝑡/𝜏,𝜏⁻¹=2𝐾Ω²(0),(22) where 𝜃 is angle of inclination of 𝝁 to Ω(𝑡). The time evolution of 𝑃(𝑡), 𝑃̇(𝑡) and expression for 𝜏 are too described by (19). It is these equations which lead to widely used exact analytic estimate of magnetic field on the neutron star pole: 𝐵=[3𝐼𝑐³/(2𝜋²𝑅⁶)]^1/2√𝑃(𝑡)𝑃̇(𝑡). For a neutron star of mass 𝑀=𝑀_⊙ and radius 𝑅=13km, one has 𝐵=3.210¹⁹√𝑃(𝑡)𝑃̇(𝑡) G. Thus, the substantial physical difference between vibration- and rotation-powered neutron star models is that in the former model the elongation of pulse period is attributed to magnetic field decay, whereas in the latter one the period lengthening is ascribed to the slow-down of rotation [24, 25].

Figure 1 

(Color online) The figure illustrates the effect of magnetic field decay on the vibration frequency and amplitude of quadrupole toroidal 𝑎-mode presented as functions of 𝑥=𝑡/𝜏.

4. Summary

The last two decades have seen a growing understanding that magnetic field decay is the key property distinguishing magnetars from pulsars. The magnetic fields frozen in the immobile matter of pulsars (most of which are fairly stable to starquakes) operate like a passive (unaltered in time) promoter of their persistent radiation along the axis of dipole magnetic moment of the star inclined to the rotation axis. The presented treatment of spindown effect emphasizes kinematic nature of time evolution of pulsar magnetic moment as brought about by rapid rotation with gradually decreasing frequency and showing that rotation does not affect intensity of its internal magnetic field. Unlike pulsars, magnetars (SGR-like sources) are isolated quaking neutron stars whose seismic and radiative activity is thought to be dominated by decay of magnetic field. Perhaps the most striking manifest of seismic vibrations of magnetar are quasiperiodic oscillations of outburst luminosity rapidly decreasing (in several hundred seconds) from about 𝐿_𝑋∼10⁴⁴ erg s⁻¹ in giant flare to about 𝐿_𝑋∼10³⁴ erg s⁻¹ in quiescent regime of long-periodic X-ray emission. As was emphasized, this lowest value of luminosity is well above of the rate of energy of rotation with period equal to that of detected emission. With all above in mind, we have set up a model of vibration-powered radiation of quaking magnetar in which the key role is attributed to magnetic field decay. The most striking outcome of presented line of argument is the gradual decrease of frequencies (lengthening of periods) of magnetodipole radiation (pulsating with the frequency of toroidal 𝑎-mode of seismic vibrations) owing its origin to magnetic field decay. Remarkably, this prediction of the model is consistent with data on gradual decrease of frequency of pulsed X-ray emission of such magnetars as 1E 2259.1+586 [26] and XTE J1810-197 [27]. It is hoped, therefore, that the theoretical results presented here can be efficiently utilized as a guide in observational quest for vibration-powered neutron stars.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant nos. 10935001 and 10973002), the National Basic Research Program of China (Grant no. 2009CB824800), and the John Templeton Foundation.