Asteroseismology of Vibration Powered Neutron Stars

Chapter from the book"Astrophysics", p.287-308. Edited by Ibrahim Kucuk, ISBN 978-953-51-0473-5, InTech, March 3, 2012.


Introduction
There is a general recognition today that basic features of asteroseismology of non-convective final stage (FS) stars, such as white dwarfs, neutron stars and and strange quark star, can be properly understood working from the model of vibrating solid star, rather than the liquid star model lying at the base of asteroseismology of convective main-sequence (MS) stars. In accord with this, most of current investigations of an FS-star vibrations rests on principles of solid-mechanical theory of continuous media, contrary to the study of the MS-star vibrations which are treated in terms of fluidmechanical theory. This means that super dense matter of FS-stars (whose gravitational pressure is counterbalanced by degeneracy pressure of constituting Fermi-matter), possesses elasticity and viscosity generic to solid state of condensed matter, whereas a fairly dilute matter of the MS (whose internal pressure of self-gravity is opposed by radiative pressure) possesses property of fluidity which is generic to the liquid state of a highly conducting condensed matter. This feature of the MS-star matter plays crucial role in generation of their magnetic fields in the dynamo processes involving macroscopic flows which are supported by energy of nuclear reactions in the central, reactive zone, of these stars. In the meantime, in the FS stars, like white dwarfs and neutron stars, there are no nuclear energy sources to support convection. The prevailing today view, therefore, is that a highly stable to spontaneous decay dipolar magnetic fields of neutron stars are fossil. The extremely large intensity of magnetic fields of degenerate solid stars is attributed to amplification of fossil magnetic field in the magnetic-flux-conserving process of core-collapse supernova.
Even still before discovery of neutron stars, it has been realized that, for absence of nuclear sources of energy, the radiative activity of these of compact objects should be powered by energy of stored in either rotation or vibrations and that the key role in maintaining the neutron star radiation should play an ultra strong magnetic field. As is commonly know today, the neutron star capability of accommodating such a field is central to understanding pulsating character of magneto-dipole radiation of radio pulsars whose radiative power is provided, as is commonly believed, by the energy of rigid-body rotation. The discovery of soft gamma-ray repeaters and their identification with magnetars [1] -quaking neutron stars endowed with ultra strong magnetic fields experiencing decay -has stimulated enhanced interest in the study of models of quake-induced magneto-mechanical seismic vibrations of neutron star and resultant electromagnetic radiation. Of particular interest in this study are torsional magneto-mechanical vibrations about axis of magnetic dipole moment of the star driven by forces of magnetic-field-dependent stresses in a perfectly conducting matter and in a permanently magnetized non-conducting matter [2]. Most, if not all, reported up to now computations of frequency spectra of poloidal and toroidal Alfvén vibration modes in pulsars and magnetars rest on tacitly adopted assumption about constant-in-time magnetic field in which a perfectly conducting neutron star matter undergoes Lorentz-force-driven oscillations [3,4,5,6,7,8,9,10,11,12] (see, also, references therein). A special place in the study of the above Alfvén modes of pure shear magneto-mechanical seismic vibrations (a-modes) occupies a homogeneous model of a solid star with the uniform density ρ and frozen-in poloidal static magnetic field of both homogeneous and inhomogeneous internal and dipolar external configuration, and we start section 2 with a brief outline of this model. In section 3, a mathematical background of quaking neutron star model is briefly outlined with emphasis on the loss of vibration energy caused by depletion of internal magnetic field pressure and resulting vibration-energy powered magneto-dipole radiation. The decreasing of magnetic field pressure in the star is presumed to be caused by coupling between vibrating star and outgoing material which is expelled by quake, but mechanisms of star-envelope interaction resulting in the decay of magnetic field, during the time of vibrational relaxation, are not considered. It is shown that physically meaningful inferences regarding radiative activity of quaking neutron star can be made even when detailed mechanisms of depletion of magnetic field pressure in the process of vibrations triggered by quakes are not exactly known. This statement is demonstrated by a set of representative examples of magnetic field decay. The basic results of the model of vibration powered neutron star are summarized in section 4 with emphasis of its relevance to astrophysics of magnetars.

Lorentz-force-driven torsion vibrations of a neutron star
To gain better understanding of the basic physics behind the interconnection between seismic and radiative activity of quaking neutron star, we start with a brief outline a fiducial model of a solid star with frozen-in homogeneous internal magnetic field and dipolar external magnetic field, as is shown in Fig.1. In the last equation, B is the field intensity [in Gauss] and b(r) stands for the dimensionless vector-function of spatial distribution of the field. The above form of B(r) has been utilized in recent works [10,11,12] in which the discrete spectra of frequencies of node-free torsional Alfvén oscillations has been computed in analytic form on the basis of equations of linear magneto-solid mechanics Here u = u(r, t) is the field of material displacement, the fundamental dynamical variable of solid-mechanical theory of elastically deformable (non-flowing) material continua. These equations describe Lorentz-force-driven non-compressional vibrations of a perfectly conducting elastic matter of a non-convective solid star with Ampére form of fluctuating current density δj. Equation for δB describing coupling between fluctuating field of material displacements u and background magnetic field B pervading stellar material is the mathematical form of Alfvén theorem about frozen-in lines of magnetic field in perfectly conducting matter. The adopted for Alfvén vibrations of non-convective solid stars terminology [which is not new of course, see for instance [13], namely, equations of magneto-solid mechanics and/or solid-magnetics is used as a solid-mechanical counterpart of well-known terms like equations of magneto-fluid mechanics, magnetohydrodynamics (MHD) and describing flowing magento-active plasma in terms of the velocity of fluctuating flow δv =u and fluctuating magnetic field δB. The MHD approach is normally utilized in astrophysics of convective main-sequence (MS) liquid stars such, for instance, as rapidly oscillating Ap (roAp) stars, chemically peculiar magnetic stars exhibiting high-frequency oscillations that have been and still remain the subject of extensive investigation [14,15]. Inserting (2) in (3) the former equation of solid-magnetics takes the form The analogy between perfectly conducting medium pervaded by magnetic field (magneto-active plasma) and elastic solid, regarded as a material continuum, is strengthened by the following tensor representation of the last equation This form is identical in appearance to canonical equation of solid-mechanics ρü i = ∇ k σ ik , where σ ik = 2µ u ik + [κ − (2/3)µ] u jj δ ik is the Hookean tensor of mechanical stresses and u ik = (1/2)[∇ i u k + ∇ k u i ] is the tensor of shear deformations in an isotropic elastic continuous matter with shear modulus µ and bulk modulus κ (having physical dimension of pressure). The most prominent manifest of physical similarity between these two material continua is their capability of transmitting non-compressional perturbation by transverse waves. Unlike incompressible liquid, an elastic solid can respond to impulsive non-compressional load by transverse waves of shear material displacements traveling with the speed c t = µ/ρ. The unique feature of an incompressible perfectly conducting and magnetized continuous matter (in liquid or solid aggregated state) is the capability of transmitting perturbation by transverse magneto-mechanical, Alfvén, wave in which mechanical displacements of material and fluctuations of magnetic field undergo coupled oscillations traveling with the speed v A = B/ √ 4πρ = 2P B /ρ (where P B = B 2 /8π is magnetic field pressure) along magnetic axis. It is stated, therefore, that magnetic field pervading perfectly conducting medium imparts to it a supplementary portion of solid-mechanical elasticity [16,17]. This suggests that hydromagnetic Alfvén vibrations of a spherical mass of a perfectly conducting matter with frozen-in magnetic field can be specified in a manner of eigenstates of elastic vibrations of a solid sphere. As for the general asteroseismology of compact objects is concerned, the above equations seems to be appropriate not only for neutron stars but also white dwarfs [18] and quark stars. The superdense material of these latter yet hypothetical compact stars is too expected to be in solid state [19,20]. Equation (4) serves as a basis of our further analysis. The studied in above works regime of node-free torsion vibrations under the action of Lorentz restoring force is of some interest in that the rate of differentially rotational material displacementṡ has one and the same form as in torsion elastic mode of node-free vibrations under the action of Hooke's force of mechanical shear stresses.
Hereafter P ℓ (cos θ) stands for Legendre polynomial of degree ℓ specifying the overtone of toroidal a-mode. Fig.2 shows quadrupole and octupole overtones of such vibrations. The time-dependent amplitude α(t) describes temporal evolution of above vibrations; the governing equation for α(t) is obtained form equation (4). The prime purpose of above works was to get some insight into difference between spectra of discrete frequencies of toroidal a-modes in neutron star models having one and the same mass M and radius R, but different shapes of constant-in-time poloidal magnetic fields. By use of the energy method, it was found that each specific form of spatial configuration of static magnetic field about axis of which the neutron star matter undergoes nodeless torsional oscillations is uniquely reflected in the discrete frequency spectra by form of dependence of frequency upon overtone ℓ of nodeless vibration. It worth noting that first computation of discrete spectra of frequencies of toroidal Alfvén stellar vibrations in the standing wave-regime, has been reported by Chandrasekhar [23]. The extensive review of other earlier computations of discrete frequency spectra of a-modes, ω ℓ = ω A s ℓ = Bκ ℓ , is given in well-known review of Ledoux and Walraven [14].
The assumption about constant in time undisturbed magnetic field means that the internal magnetic field pressure, P B , the velocity v A of Alfvén wave in the star bulk and, hence, the frequency Alfvén oscillations remain constant in the process of vibrations whose amplitude α(t) subjects to standard equation of undamped harmonic oscillator [10,11]. The allow for viscosity of stellar material leads to exponential damping of amplitude, but the frequency ν A and, hence, the period P A preserve one and the same values as in the case of non-viscous vibrations [12]. In Fig.3 these latter quantities are plotted as functions of intensity B of undisturbed poloidal magnetic field in the neutron star models with indicated mass M and radius R. The practical usefulness of chosen logarithmic scale in this figure is that it shows absolute vales of ν A and P A for global vibrations of typical in mass and radius neutron stars.
In what follows, we relax the assumption about constant-in-time magnetic field and examine the impact of its decay on the vibration energy and period. A brief analysis of such a case has been given in the context of magnetic white dwarfs [21]. In this paper we present a highly extensive consideration of this problem in the context of neutron stars with emphasis on its relevance to the post-quake radiation of magnetars.

Vibration-energy powered magneto-dipole emission of neutron star
In the following we consider a model of a neutron star with time-dependent intensity of homogeneous poloidal magnetic field which can be conveniently represented in the form On account of this the equation of solid-magnetics, (4), takes the form Inserting here the following separable form of fluctuating material displacements Scalar product of (12) with the time-independent field of instantaneous displacements a(r) followed by integration over the star volume leads to equation for amplitude α(t) having the form of equation of oscillator with depending on time spring constant The solution of equation of non-isochronal (non-uniform in duration) and non-stationary vibrations with time-dependent frequency [α(t) + ω 2 (t)α(t) = 0 where ω 2 (t) = K(t)/M] is non-trivial and fairly formidable task [22]. But solution of such an equation, however, is not a prime purpose of this work. The main subject is the impact of depletion of magnetic-field-pressure on the total energy of Alfvén vibrations and the discrete spectrum of frequency of the toroidal a-mode It is to be stated clearly from the onset that it is not our goal here to speculate about possible mechanisms of neutron star demagnetization and advocate conceivable laws of magnetic field decay. The main purpose is to gain some insight into the effect of arbitrary law of magnetic field decay in quaking neutron star on period of Lorentz-force-driven torsional seismic vibrations (whose quadrupole and octupole overtones are pictured in Fig.2) and radiative activity of the star brought about by such vibrations. In the reminder of the paper we focus on the case of torsional Alfvén vibrations in quadrupole (ℓ = 2) overtone. In so doing we omit index ℓ putting ω(t) = ω ℓ=2 (t).

Magnetic-field-decay induced loss of vibration energy
The total energy stored in quake-induced Alfvén seismic vibrations of the star is given by Perhaps most striking consequence of the magnetic-field-pressure depletion during the post-quake vibrational relaxation of neutron star is that it leads to the loss of vibration energy at a rate proportional to the rate of magnetic field decay In the model under consideration, a fairly rapid decay of magnetic field during the time of postquake vibrational relaxation of the star is though of as caused, to a large extent, by coupling of the vibrating star with material expelled by quake. In other words, escaping material removing a part of magnetic flux density from the star is considered to be a most plausible reason of depletion of internal magnetic field pressure in the star. It seems quite likely that, contrary to viscous dissipation, the loss of vibration energy due to decay of magnetic field must be accompanied by coherent (non-thermal) electromagnetic radiation. Adhering to this supposition in the next section special consideration is given to the conversion of the energy of Lorentz-force-driven seismic vibrations into the energy of magneto-dipole emission whose flux oscillates with frequency of torsional Alfvén magneto-mechanical vibrations of the final stage solid stars, like magnetic white dwarf and neutron stars. Fig.4 replicates seismic torsional Alfvén vibrations of the star which are accompanied by oscillations of lines of dipolar magnetic field defining the beam direction of outburst X-ray emission.

Conversion of vibration energy into power of magneto-dipole radiation
The point of departure in the study of vibration-energy powered magneto-dipole emission of the star (whose radiation power, P, is given by Larmor's formula) is the equation Consider a model of quaking neutron star whose torsional magneto-mechanical 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 magneto-mechanical vibrations of stellar matter which are described by equation for α(t). This means that δµ(t) and α(t) must obey equations of similar form, namely Figure 4: Schematic view of quadrupole overtone of seismic torsional Alfvén vibrations of neutron star with homogeneous internal and dipolar external field whose lines, defining direction of outburst beam, oscillate with the frequency of this seismic a-mode.
It is easy to see that equations (23) and (24) can be reconciled if Then, from (23), it follows δμ = −ω 2 µα. Given this and equating we arrive at the equation of time evolution of magnetic field which yields the following law of its decay The lifetime of magnetic field τ is regarded as a parameter whose value is established from below given relations between the period P and its time derivativeṖ which are taken from observations. Knowing from observations P (t) andṖ (t) and estimating τ one can get information about the magnitude of total magnetic moment and the strength of undisturbed magnetic field. It is worth noting that in the model of vibration-energy powered magneto-dipole emission under consideration, the equation of magnetic field evolution is obtained in similar fashion as equation for the angular velocity Ω does in the standard model of rotation-energy powered emission of neutron star. As is shown in the next section, the substantial physical difference between models of rotationenergy and vibration-energy powered pulsating emission of neutron stars is that in the model of quaking neutron star vibrating in toroidal a-mode, the elongation of period of pulses is attributed to magnetic field decay, whereas in canonical Pacini-Gold model of radio-pulsar the lengthening of period of pulses is ascribed to the slow down of the neutron star rotation [39,40,49].

Lengthening of vibration period
The immediate consequence of above line of argument is the magnetic-field-decay induced lengthening of vibration period It follows that lifetime τ is determined by This inference of the model under consideration is demonstrated in Fig.5. The difference between periods evaluated at successive moments of time t 1 = 0 and t 2 = t is given by The practical usefulness of these general relations is that they can be used as a guide in search for fingerprints of Alfvén seismic vibrations in data on oscillating emission from quaking neutron star.
All the above shows that the magnetic-field-decay induced loss of vibration energy is substantially different from the vibration energy dissipation caused by shear viscosity of matter resulting in heating of stellar material [26,27]. As was noted, the characteristic feature of this latter mechanism of vibration energy conversion into the heat (i.e., into the energy of non-coherent electromagnetic emission responsible for the formation of photosphere of the star) is that the frequency and, hence, period of vibrations are the same as in the case of viscous-free vibrations [12]. However, it is no longer so in the case under consideration. It follows from above that depletion of magnetic field pressure resulting in the loss of total energy of Alfvén vibrations of the star causes its vibration period to lengthen at a rate proportional to the rate of magnetic field decay.

Numerical analysis
To get an idea of the magnitude of characteristic parameter of vibrations providing energy supply of magneto-dipole radiation of a neutron star, in Table 1 we present results of numerical computations of the fundamental frequency of neutron star oscillations in quadrupole toroidal a-mode and time of decay of magnetic field τ as functions of increasing magnetic field. As was emphasized, the most striking feature of considered model of vibration powered radiation is the lengthening of periods of pulsating emission caused by decay of internal magnetic field. This suggests that this model is relevant to electromagnetic activity of magnetars -neutron stars endowed with magnetic field of extremely high intensity the radiative activity of which is ultimately related to the magnetic field decay. Such a view is substantiated by estimates of Alfvén frequency presented in the table. For magnetic fields of typical rotation powered radio pulsars, B ∼ 10 12 G, the computed frequency ν A is much smaller than the detected frequency of pulses whose origin is attributed to lighthouse effect. In the meantime, for neutron stars with magnetic fields B ∼ 10 14 G the estimates of ν A are in the realm of observed frequencies of high-energy pulsating emission of soft gamma repeaters (SGRs), anomalous X-ray pulsars (AXPs) and sources exhibiting similar features. According to common belief, these are magnetars -highly magnetized neutron stars whose radiative activity is related with magnetic field decay. The amplitude of vibrations is estimated as [28]   The luminosity powered by neutron star vibrations in quadrupole toroidal a-mode is given by The presented in Fig.8 computations of power of magneto-dipole radiation of a neutron stars (of one and the same mass and radius but different values of magnetic fields) exhibit oscillating character of luminosity. The frequency of these oscillations equal to that of torsional Alfvén seismic vibrations of neutron star. All above suggests that developed theory of vibration-energy powered emission of neutron star is relevant to electromagnetic activity of magnetars -neutron stars endowed with magnetic field of extremely high intensity the radiative activity of which is ultimately related to the magnetic field decay. This subclass of highly magnetized compact objects is commonly associated with soft gamma repeaters and anomalous X-ray pulsars [29,30,31,32] -young isolated and seismically active neutron stars [33,34]. The magnetar quakes are exhibited by short-duration thermonuclear gammaray flash followed by rapidly oscillating X-ray flare of several-hundred-seconds duration. During this latter stage of quake-induced radiation of magnetar, a long-periodic (2-12 sec) modulation of brightness was observed lasting for about 3-4 minutes. Such long period of pulsed emission might be expected from old rotation-energy powered pulsars (due to slow down of rotation), but not from magnetars which are young neutron stars, as follows from their association with pretty young supernovae. Taking this into account (and also the fact that energy release during X-ray flare is much larger than the energy of rigid-body rotation with so long periods) it has been suggested [2] that detected long-periodic pulsating emission of magnetars is powered by the energy of torsional magneto-elastic vibrations triggered by quake. In view of key role of ultra-strong magnetic field it is quite likely that quasi-periodic oscillations (QPOs) of outburst flux from SGR 1806-20 and SGR 1900+14 detected in [35,36,37] are produced by torsional seismic vibrations predominately sustained by Lorentz force [10,11,12] which are accompanied, as was argued above, by monotonic decay of background magnetic field. If so, the predicted elongation of QPOs period of oscillating outburst emission from quaking magnetars should be traced in existing and future observations.

Comparison with the rotation powered neutron star
For a sake of comparison, in the considered model of vibration powered radiation, the equation of magnetic field evolution is obtained in similar fashion as that for the angular velocity Ω(t) does in the standard model of rotation powered neutron star which rests on equations [38] One of the basic postulates of the model of rotation-energy powered emission is that the time evolution of total magnetic moment of the star is governed by the equation It follows where θ is angle of inclination of µ to Ω(t). The total magnetic moment of non-rotating neutron star is parametrized by equation of uniformly magnetized (along the polar axis) sphere This parametrization presumes that in the rotation powered neutron star, the frozen-in the star magnetic field operates like a passive promoter of magneto-dipole radiation, that is, intensity of the internal magnetic field remain constant in the process of radiation. As a result, the equation of energy conversion from rotation to magnetic dipole radiation is reduced to equation of slow down of rotationΩ where θ is angle of inclination of µ to Ω(t). From the last equation it follows and, hence, the lifetime τ is related with P (t) andṖ (t) as Equating two independent estimates for τ given by equations (54) and (57) we arrive at widely utilized analytic estimate of magnetic field on the neutron star pole: For a neutron star of mass M = M ⊙ , and radius R = 13 km, one has B = 3.2 10 19 P (t)Ṗ (t), G.
The outlined treatment of rotation-powered magneto-dipole radiation of a neutron star emphasizes kinematic nature of variation of magnetic moment of pulsar whose magnetic field is regarded as independent of time. Thus, the substantial physical difference between vibration-powered and rotation-powered neutron star models is that in the former the elongation of pulse period is attributed to magnetic field decay, whereas in the latter the period lengthening is ascribed to slow down of rotation [39,40,41].

Comment on magnetic field decay in quaking neutron star
The considered law of magnetic field decay cannot be, of course, regarded as universal because it reflects a quite concrete line of argument regarding the fluctuations of magnetic moment of the star. The interrelation between quake-induced oscillations of total magnetic moment of the star, δµ(t), and the amplitude, α(t), of its seismic magneto-mechanical oscillations can be consistently interpreted with the aid of the function of dipole demagnetization f(B(t)) which is defined by the following condition of self-consistency in α of right and left hand sides of equation (22), namely The vector-function of dipole demagnetization f(B(t)) depending on decaying magnetic field reflects temporal changes of electromagnetic properties of neutron star matter as well as evolution of magnetic-field-promoted coupling between neutron star and its environment. This means that specific form of this phenomenological function should be motivated by heuristic arguments taking into account these factors. With this form of δμ(t), equation (22) is transformed to magnetic field decay of the form In the above considered case, this function is given by The practical usefulness of the dipole demagnetization function, f(B(t)), consists in that it provides economic way of studying a vast variety of heuristically motivated laws of magnetic field decay, B = B(t), whose inferences can ultimately be tested by observations. In this subsection, with no discussing any specific physical mechanism which could be responsible for magnetic field decay, we consider a set of representative examples of demagnetization function f(B(t)) some of which have been regarded before, though in a somewhat different context [42][43][44][45][46][47][48][49]. Here we stress again that the model under consideration deals with magnetic field decay in the course of vibrations triggered by starquake, not with long-term secular decay which has been the subject of these latter investigations (see also references therein).   1. As a first representative example, a model of quaking neutron star whose function of dipole demagnetization has the following form is been considered. In such demagnetization model, the temporal evolution of the magnetic field and the lengthening of the vibration period obey the following laws where B(0) is the intensity of magnetic field before quake. 2. For a quaking neutron star model whose function of demagnetization is given by the resultant equation of magnetic field evolution and the vibration period elongation read which lead to a fairly non-trivial logarithmic law of magnetic field decay and vibration period lengthening For each of above examples, the magnetic field and the resultant lengthening of vibration period, computed as functions of fractional time t/τ , are shown in Fig. 9. This third example is interesting in that computed in this model ratioν(t)/ν(0), pictured in Fig.10, as a function of t, is similar to that which exhibit data on post-glitch emission of PSR J1846-0258 [50]. These examples show that the period elongation is the common effect of magnetic field decay. The interrelations between periods and its derivatives substantially depend on specific form of the magnetic field decay law although physical processes responsible for magnetic field decay remain uncertain. Practical significance of above heuristic line of argument is that it leads to meaningful conclusion (regarding elongation of periods of oscillating magneto-dipole emission) even when detailed mechanisms of magnetic field decay in the course of vibrations are not exactly known. ν(x)/ν(0) Figure 10: The rate of frequencyν(x) normalized to ν(0) as a function of x = t/τ for logarithmic law of magnetic filed decay.

Summary
It is generally realized today that the standard model of inclined rotator, lying at the base of our understanding of radio pulsars, faces serious difficulties in explaining the long-periodic (2 < P < 12 s) pulsed radiation of soft gamma repeaters (SGRs) and anamalous X-ray pulsars (AXPs).
Observations show that persistent X-ray luminosity of these sources (10 34 < L X < 10 36 erg s −1 ) is appreciably (10-100 times) larger than expected from neutron star deriving radiation power from energy of rotation with frequency of detected pulses. It is believed that this discrepancy can be resolved assuming that AXP/SGR-like sources are magnetars -young, isolated and seismically active neutron stars whose energy supply of pulsating high-energy radiation comes not from rotation (as is the case of radio pulsars) but from different process involving decay of ultra strong magnetic field, 10 14 < B < 10 16 G. Adhering to this attitude we have presented the model of quaking neutron star deriving radiation power from the energy of torsional Lorentz-force-driven oscillations. It is appropriate to remind early works of the infancy of neutron star era [51,52] in which it has been pointed out for the first time that vibrating neutron star should operate like Hertzian magnetic dipole deriving radiative power of magneto-dipole emission from the energy of magneto-mechanical vibrations [53]. What is newly disclosed here is that the main prerequisite of the energy conversion from vibrations into radiation is the decay of magnetic field in the star. Since the magnetic field decay is one of the most conspicuous features distinguishing magnetars from rotation powered pulsars, it seems meaningful to expect that at least some of AXP/SGR -like sources are magnetars deriving power of pulsating magnetic dipole radiation from the energy of quake-induced torsion seismic vibrations.