# Vibration-Powered Radiation of Quaking Magnetar

Authors: S. Bastrukov; J. W. Yu; R. X. Xu; I. Molodtsova
Journal: ISRN Astronomy and Astrophysics
 (2011)

---

## 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.


---

## Body

## 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<P<12 s) -pulsed radiation of soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs). The persistent X-ray luminosity of AXP/SGR sources (1034<LX<1036 erg s−1) 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 (1014<B<1016 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 a-modes in a neutron star with mass M, radius R, and magnetic field of typical magnetar, B14=B/1014 G, reads
(1)ωl=ωAsl,ωA=vAR=BR3M,M=4π3ρR3,(2)νA=ωA2π=0.2055B14R61/2,(MM⊙)-1/2,Hz,
where numerical factor sℓ is unique to each specific shape of magnetic field frozen in the neutron star of one and the same mass M and radius R.


## 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 fieldB can be properly described in terms of material displacements u obeying equation of magneto-solid-mechanics
(3)ρü(r,t)=14π[∇×[∇×[u(r,t)×B(r,t)]]]×B(r,t),(4)u̇(r,t)=[ω(r,t)×r],ω(r,t)=At[∇χ(r)]α̇(t).
The field u̇(r,t) is identical to that for torsion node-free vibrations restored by Hooke’s force of elastic stresses [18, 21] with χ(r)=Aℓfℓ(r)Pℓ(cosθ), where fℓ(r) is the nodeless function of distance from the star center and Pℓ(cosθ) is Legendre polynomial of degree ℓ specifying the overtone of toroidal mode. In (4), the amplitude α(t) 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)B(r,t)=B(t)b(r),
where B(t) is the time-dependent intensity and b(r) is dimensionless vector-function of the field distribution over the star volume. Scalar product of (1) with the separable form of material displacements
(6)u(r,t)=a(r)α(t)
followed by integration over the star volume leads to equation for amplitude α(t) having the form of equation of oscillator with time-dependent spring constant
(7)Mα̈(t)+K(t)α(t)=0.
The total vibration energy and frequency are given by
(8)EA(t)=Mα̇2(t)2+K(B(t))α2(t)2,ω(t)=K(t)M=B(t)κ,κ=R3Ms.
It follows that
(9)dEA(t)dt=α2(t)2dK(B)dBdB(t)dt.
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
(10)dEA(t)dt=-P(t),P(t)=23c3δμ̈2(t).
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: μ=μn=constant. The total magnetic dipole moment should execute oscillations with frequency ω(t) equal to that for magnetomechanical vibrations of stellar matter, which are described by equation for α(t). This means that δμ(t) and α(t) must obey equations of similar form, namely,
(11)δμ̈(t)+ω2(t)δμ(t)=0,α̈(t)+ω2(t)α(t)=0,ω2(t)=B2(t)κ2.
It is easy to see that (11) can be reconciled if
(12)δμ(t)=μα(t).
Given this, we arrive at the following law of magnetic field decay:
(13)dB(t)dt=-γB3(t),γ=2μ2κ23Mc3=constant,(14)B(t)=B(0)1+t/τ,τ-1=2γB2(0).
The last equation shows that the lifetime of quake-induced vibrations in question substantially depends upon the intensity of initial (before quake) magnetic field B(0): the larger the B(0) the shorter the τ. For neutron stars with one and the same mass M=1.4M⊙ and radius R=15km, and magnetic field of typical pulsar B(0)=1012 G, we obtain τ~3×107 years, whereas for magnetar with B(0)=1014 G, τ~7×103 years.The equation for vibration amplitude(15)α̈(t)+ω2(t)α(t)=0,ω2(t)=ω2(0)1+t/τ,ω(0)=ωAκ.
with help of substitution s=1+t/τ is transformed to
(16)sα′′(s)+β2α(s)=0,β2=ω2(0)τ2=constant
permitting general solution [22]. The solution of this equation, obeying two conditions α(t=0)=α0 and α(t=τ)=0, can be represented in the form
(17)α(t)=Cs1/2{J1(z(t))-ηY1(z(t))},z(t)=2βs1/2(t),η=J1(z(τ))Y1(z(τ)),C=α0[J1(z(0))-ηY1(z(0))]-1,
where J1(z) and Y1(z) are Bessel functions [23] and
(18)α02=2E¯A(0)Mω2(0)=2E¯A(0)K(0),ω2(0)=K(0)M.
Here by E¯A(0), the average energy stored in torsional Alfvén vibrations of magnetar is understood. If all the detected energy Eburst of X-ray outburst goes in the quake-induced vibrations, Eburst=EA, then the initial amplitude α0 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 α(t) with pointed out parameters β and η. The magnetic-field-decay-induced lengthening of period of pulsating radiation (equal to period of vibrations) is described by
(19)P(t)=P(0)[1+(tτ)]1/2,Ṗ(t)=12τP(0)[1+(t/τ)]1/2,τ=P2(0)2P(t)Ṗ(t),P(0)=2πκB(0).
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 M=1.4M⊙ and radius R=10km vibrating in quadrupole overtone of toroidal a-mode is given by
(20)μ=A(M,R)P(t)Ṗ(t),A=3Mc38π2,μ=3.8×1037P(t)Ṗ(t)Gcm3.
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 Ω(t) in the standard model of rotation-powered neutron star which rests on
(21)dERdt=-23c3δμ̈2(t),ER(t)=12IΩ2(t),I=25MR2,δμ̈(t)=[Ω(t)×[Ω(t)×μ]],μ=μn,μ=12BR3
which lead to
(22)Ω̇(t)=-KΩ3(t),K=2μ⊥23Ic3,μ⊥=μsin⁡θ,Ω(t)=Ω(0)1+t/τ,τ-1=2KΩ2(0),
where θ is angle of inclination of μ to Ω(t). The time evolution of P(t), Ṗ(t) 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: B=[3Ic3/(2π2R6)]1/2P(t)Ṗ(t). For a neutron star of mass M=M⊙ and radius R=13km, one has B=3.21019P(t)Ṗ(t) 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 toroidala-mode presented as functions of x=t/τ.


## 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 aboutLX~1044 erg s−1 in giant flare to about LX~1034 erg s−1 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 a-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.



---

*Source: 843825-2012-01-05.xml*