arxiv:astro-ph/ v1 21 Sep 1999

Transcription

1 Mon. Not. R. Aston. So. 000, (0000) Pinted 18 Apil 2008 (MN LATEX style file v1.4) Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. A.N. Timokhin, 1,2,3 G.S. Bisnovatyi-Kogan, 2 H.C. Spuit 3 1 Stenbeg Astonomial Institute, Univesitetskij p. 13, Mosow, Russia 2 Spae Reseah Institute, Pofsoyuznaya 84/32, Mosow, Russia 3 Max-Plank-Institut fü Astophysik, Kal-Shwazshild St.1, D Gahing, Gemany axiv:asto-ph/ v1 21 Sep Apil INTRODUCTION ABSTRACT We genealize a fomula fo the Goldeih-Julian hage density (ρ GJ ), oiginally deived fo otating neuton sta, fo abitay osillations of a neuton sta with abitay magneti field onfiguation unde assumption of low uent density in the inne pats of the magnetosphee. As an appliation we onside tooidal osillation of a neuton sta with dipole magneti field and alulate enegy losses. Fo some osillation modes the longitudinal eleti field an not be aneled by putting haged patiles in the magnetosphee without a pesene of stong eleti uent (j ω ρ GJ ). It is shown that the enegy losses ae stongly affeted by plasma in the magnetosphee, and annot be desibed by vauum fomulas. Key wods: stas:neuton, osillations, magneti fields pulsas:geneal Osillations of neuton stas ould give an oppotunity to study thei intenal stutue. Osillations of a neuton sta (NS) may be obsevable in pulsas just afte a glith. The enegy elease duing the glith is estimated as E 1 M ust IΩ 2 Ω 2 M NS Ω, whee I is the moment of inetia of the NS, Ω is an angula veloity of otation, Ω is a jump of angula veloity, M NS is a mass of the NS, M ust is a mass of the ust. Fo the Vela pulsa E egs. Thee ae about 20 glithing pulsas, whih had suffeed a total of 45 glithes (Lyne 1996). The enegy elease in eah of these glithes is simila to that of the Vela pulsa. If seveal peent of this enegy goes into exitation of osillations it would be possible to obseve them. Anothe example how osillations of NS ould be obseved in pulsas is by obsevations of miostutue in single pulses. Boiakoff (1976) poposed that vibations of a neuton sta may ause the peiodiity of miopulses obseved in some pulsas. A diffiulty in this theoy is that no effetive mehanism has been poposed fo exitation of osillations (Hankins 1996). Osillating neuton stas have been poposed as a soue fo Galati Gamma-Ray Busts (GRBs) by Pahini & Rudeman (1974) and Tsygan (1975). This idea is futhe developed by Blaes et al. (1989), Smith & Epstein (1993), Fatuzzo & Melia (1993). Osillation indued had gamma adiation fom NS was poposed by Bisnovatyi-Kogan (1995) and Cheng & Ding (1997) fo explanation of the had delayed emission obseved fom some GRBs (Huley et al. 1994), even if the Gamma-Ray Bust itself is geneated by a diffeent mehanism (Bisnovatyi-Kogan et al. 1975). Reent GRBs obsevations make a Galati oigin of GRBs unlikely, but a Galati model may not be exluded ompletely. At least fo Soft Gamma Repeates, whih ae believed to be neuton stas, osillations of sta an play an essential ole (Dunan 1998). Eigenfequenies and eigenfuntions of neuton sta osillations have been omputed by many authos (MDemott et al. 1988; Caoll et al. 1986; Stohmaye 1991). These alulations have shown that typial peiods of neuton sta osillations ange fom 0.1 up to seveal tens of milliseonds. On the othe hand, obseved pulsa adiation in diffeent spetal egions is geneated mostly in thei magnetosphees (Tümpe & Beke 1998). If we ae looking fo osillation indued adiation fom pulsas we should investigate poesses in the magnetosphee. A omplete solution fo vauum eletomagneti fields (1) Cuent addess: Max-Plank-Institut fü Kenphysik, Saupfehekweg 1, D Heidelbeg, Gemany 0000 RAS

2 2 A.N.Timokhin et al. nea an osillating magnetized neuton sta was obtained by Muslimov and Tsygan (1986). But fo typial NS s magneti field stength (> 10 9 G) the eleti field aising fom sta s osillation will be stong enough to pull haged patiles fom the sufae in the magnetosphee. Hene, any ealisti model of the magnetosphee of osillating NS must take into aount the pesene of haged patiles. This will also affet the eletomagneti enegy losses of a pulsating magnetized sta. Thee wee alulations of enegy losses by osillating NS in vauum (MDemott et al. 1984; Muslimov & Tsygan 1986). As a fist step in the investigation of magnetosphee of an osillating neuton sta with nonzeo hage density we onside the inne pat of the magnetosphee, assuming low uent density thee. We genealize a fomula fo Goldeih-Julian hage density, deived oiginally fo otating NS, fo a NS osillating in any abitay mode. This allow us to take into aount plasma pesent in the magnetosphee of the NS and alulate the enegy losses. The plan of the pape is following. In the fist setion we intodue some definitions and give algoithm fo alulation of Goldeih-Julian hage density fo abitay osillation mode and magneti field stutue unde the assumption of low uent density, and disuss the limitations of this assumption. In the seond setion we apply this fomalism to tooidal osillations of a NS with a dipole magneti field. We alulate the Goldeih-Julian hage density and enegy losses fo some osillation modes and disuss fatos affeting enegy losses of osillating neuton sta. 2 GOLDREICH-JULIAN CHARGE DENSITY FOR OSCILLATING NS 2.1 Basi definitions and equations fo the Goldeih-Julian hage density We assume a neuton sta to be a magnetized onduting sphee of a adius. We ae inteested only in the pulsation modes with non vanishing amplitude ξ on the sufae. Conside a egion lose to the NS (nea zone), at distanes fom the NS sufae smalle than the wave length λ = 2π/ω, whee ω is osillation fequeny and is the speed of light. In the nea zone < λ one an neglet the displaement uent tem. Let us also assume that the physial uent density in the inne pat of the magnetosphee is small enough, i.e. magneti field to the fist ode in ξ/ an be onsideed as geneated only by volume uents inside the sta and by sufae uents on the NS s sufae. Hene, Maxwell s equations fo eleti and magneti fields to the fist ode of ξ ξ/ in the nea zone ae: E = 4πρ, (2) E = 1 tb, B = 0, B = 0, whee B and E ae magneti and eleti fields and ρ is the hage density. We use in the notation fo patial deivatives: α / α. On the unpetubed sufae of the NS E and B must satisfy the following bounday onditions: B ( ) = B 0, (3) (4) (5) (6) B θ,φ ( ) = B 0θ,0φ ± 4π J φ,θ, E θ,φ ( ) = 1 (v B0) θ,φ, E ( ) = 1 (v B0) + 4πσ, (7) (8) (9) whee subsipts, θ, φ denote veto omponents in a spheial oodinate system, v is the veloity of osillation of the NS sufae, B 0 is the sufae magneti field inside the NS, σ and J ae indued sufae hage and uent density. Hee (9) and (7) ae used fo detemination of the sufae hage and uent densities oespondingly. Using the popeties of solenoidal veto fields we wite B in the fom B = (Pe ) + (Te ), whee e = /. Sala funtions P and T an be expanded in spheial hamonis Y lm as P(, t) = P lm (,t) Y lm (θ, φ) (11) T(, t) = lm T lm (, t) Y lm (θ, φ). lm Substituting B in the fom (10) into the equation (5) and multiplying the esult by e we get ΩT = 0, whee Ω is the angula pat of the Laplaian, (10) (12) (13)

3 Ω 1 sin θ θ(sinθ θ ) + 1 sin 2 θ φφ. Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 3 Substitution of the expansion (12) into the equation (13) gives us T lm 0 fo eah l, m. Hene, T 0 (see Muslimov & Tsygan (1986)) and the magneti field an be expessed in tems of only one sala funtion P as (14) B = (Pe ). (15) Substituting B fom expession (15) in the Faaday s law (3) we get E = 1 ( tpe). (16) Integating this equation: E = 1 ( tpe) Ψ, (17) whee Ψ is an abitay sala funtion. Fom the theoy of patial diffeential equations (Elsgolts 1965; Bandt 1947) it is known that fo any veto field A thee exist a veto field E GJ pependiula to A if and only if A ( A) = 0. (18) Fom equation (5) it follows that magneti B satisfies equation (18). Let us assume that thee exists an eleti field E GJ with no omponent paallel to B satisfying the bounday onditions (8),(9). Evidently it satisfies Maxwell equations (2)-(5) and onsequently has the fom E GJ = 1 ( tpe) ΨGJ. (19) The eleti field E GJ satisfies equation E GJ B = 0. (20) Substituting E GJ fom the expession (19) in the equality (20) we have an equation fo Ψ GJ 1 ( tpe) B + ΨGJ B = 0. Let us all E GJ the Goldeih-Julian (G-J) eleti field. Intodue the genealized Goldeih-Julian hage density ρ GJ as that a hage density that is neessay to ompensate a longitudinal eleti field in the magnetosphee, aising fom the fist, vauum, tem in the expession (19). Fom the equation (2) it follows that ρ GJ = 1 4π EGJ. Using fomula (19) fo E GJ we have fo ρ GJ ρ GJ = 1 4π ΨGJ. In othe wods, if one uses epesentation (15) fo a magneti field in the nea zone, then the Goldeih-Julian eleti field is witten as a sum of two tems (19). The fist one is the vauum tem, and the seond one ( Ψ GJ) epesents the ontibution of the haged patiles in the magnetosphee of the NS. The potential Ψ GJ is a solution of the equation (21). In spheial oodinates, equation (21) beomes (using expession [15]): ΩP Ψ GJ θ P θ Ψ GJ 1 sin 2 θ φp φ Ψ GJ + 1 sin θ ( φp θ tp θ P φ tp) = 0. (24) Let us onside small osillations of a NS, ξ 1. We expand the funtion P in a seies in ξ and appoximate it by the sum of the fist two tems: (21) (22) (23) P(t,, θ, φ) P 0(,θ, φ) + δp(t,, θ, φ), (25) whee the funtion P 0(,θ, φ) is esponsible fo the unpetubed magneti field and δp(t,, θ, φ) is the fist ode tem in an expansion of P(t,, θ, φ) in ξ. We expand δp in of spheial hamonis δp = δ p lm (t,) Y lm (θ, φ). (26) l,m Subsequently we will need only the time-deivative of δp. Afte substitution of the expansion (26) in the equation (5) we find (see Appendix A)

4 4 A.N.Timokhin et al. ( ) l t δp = tδp lm (t)y lm (θ, φ). l,m (27) The oeffiients tδp lm (t) ae fixed by the bounday onditions. 2.2 Bounday onditions The bounday onditions fo the potential Ψ GJ ae obtained fom the bounday onditions fo eleti field (8). The tangential omponents of the G-J eleti field on the sufae of the NS ae Eθ GJ = = 1 ( 1 1 ΩP v φ + 1 ) sin θ φp v (28) = Eφ GJ = = 1 ( ) 1 1 ΩP v θ + θ P v. (29) = Fom the expession (19) the tangential omponents of the GJ eleti field outside the NS ae Eθ GJ = 1 ( ) 1 1 sin θ φ tp + θ Ψ GJ Eφ GJ = 1 ( 1 θ tp 1 ) sin θ φψ GJ. Equating these expessions on the sufae of the NS and eliminating Ψ GJ one obtains bounday onditions fo the θ and φ (30) (31) deivatives of Ψ GJ θ Ψ GJ = = 1 ( 1 ΩP v φ + 1 sin θ φp v + 1 sin θ φ tp) = (32) φ Ψ GJ = = 1 ( ) 1 sin θ ΩP v θ + sin θ θ P v + sin θ θ tp. (33) = Fom these bounday onditions we obtain the bounday ondition to be applied to δp lm (t). Diffeentiating equation (32) with espet to φ and equation (33) with espet to θ and equating esults we get an expession fo δp lm (t) to fist ode in ξ (see Appendix B) tδp lm (t) = 1 l(l + 1) 4π [ ] dω Ylm v ( ΩP 0) + ΩP 0 ( v ) + 2 ( ( P 0) ) v. (34) = Bounday ondition fo Ψ GJ an be obtained by integating equation (32) o equation (33). Fo onveniene we will use as bounday ondition the esult of integating equation (32) ove θ. Fo a petubation depending on time t as e iωt we get to fist ode in ξ Ψ GJ = = 1 ( 1 dθ ΩP0 v φ + 1 sin θ φp 0 v + 1 sin θ φ tδp) + e iωt F(φ), (35) = whee F is a funtion of φ. Fo eah osillation mode the oesponding veloity field is ontinuously diffeentiable. Fom the bounday ondition fo the eleti field (8) it follows that the tangential omponents of E GJ ae finite. The vauum tem on the the ight hand side in expession (19) fo the eleti field E GJ is also finite (with the natual assumption that P is ontinuously diffeentiable with espet to θ and φ expessions (34) and (27) ae finite fo = ). Let us onside the azimuthal omponent of G-J eleti field E GJ φ nea the poles. In expession (31) fo E GJ φ of expession (29), ae finite. Consequently the seond tem ( φψ GJ sin θ outside the NS both Eφ GJ and the fist (vauum tem) on the ight side ) is also finite, hene φ Ψ GJ {θ=0,π; = } = 0. Thus = in the bounday ondition (35) one must hoose F suh that Ψ GJ θ=0,π; = = C 1,2 e iωt, whee C 1,2 ae some onstants. Using the gauge feedom we hoose Ψ GJ θ=0; = = 0. (36) Using appoximation (25) fo P we wite an equation fo Ψ GJ in the ase of small osillations ΩP 0 Ψ GJ θ P 0 θ Ψ GJ 1 sin 2 θ φp 0 φ Ψ GJ + 1 sin θ ( φp 0 θ tδp θ P 0 φ tδp) = 0. (37) points with θ = 0(π), =

5 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 5 The fist-ode patial diffeential equation (37) togethe with the bounday onditions (35), (36) and tδp desibed by fomulas (27), (34) detemines Ψ GJ. The solution of this poblem povides Ψ GJ to fist ode in ξ. The Goldeih-Julian hage density ρ GJ an then be obtained fom the fomula (23). Any deviation of the hage density in the magnetosphee of NS fom ρ GJ will lead to a geneation of a longitudinal eleti field and to aeleation of haged patiles. As we will show below, thee does not exist a well behaved solution of equations (37), (35), (36) fo all osillation modes. Fo some modes ρ GJ has singulaities. We now examine ou assumption that the uent does not influene the magneti field in the nea zone to fist ode in ξ. Fo this to be the ase, we must have 4π j B(1), whee B (1) is the fist ode in ξ petubation of the magneti field. To ode of magnitude this means: j 1 ( ) ξ B 0 = B0ωξ ( ). ω The maximum value of the Goldeih-Julian hage density on the sufae of the NS is, by ode of the magnitude, ρ max GJ ( ) B 0 V, whee V is osillation veloity amplitude. Thus we must equie ( ) j ρ GJ( ). (40) ω Fo a hage-sepaated solution, when in any egion of the nea zone thee ae haged patiles of only one sign of hage, maximum uent density is ρ GJ. The absolute value of the Goldeih-Julian hage density deeases with ineasing of and in the nea zone /ω, onsequently ondition (40) fo hage-sepaated solution is satisfied automatially. Fo situation when in some egions of the nea zone thee ae haged patiles of diffeent sign of hage ondition (40) an be violated. We believe in existene of a egula foe-fee solution at least in the inne pat of the magnetosphee E B = 0 (41) fo osillating NS beause in pesene of longitudinal eleti field haged patiles ae aeleated to high enegies and thei uvatue adiation in stong magneti field podue e + e pai asade (Stuok 1971; Daughety & Hading 1982). The patiles podued in the asade anel the aeleating eleti field. The pesene of singulaities in the solution tell us that assumption of a negligible uent density in the sense of inequality (40) fails. Fo osillation modes whee ρ GJ has singulaities thee ae stong eleti uents in the nea zone, ( ) j ρ GJ( ) (42) ω in ode to satisfy ondition (41). This situation will be onsideed in subsequent papes. (38) (39) 2.3 Rotating NS Pulsa At the end of this setion we onside an impotant patiula ase - a otating neuton sta. In this ase the patial deivative t an be eplaed by Ω φ, whee Ω is the angula veloity. Thus tp = Ω φ P. (43) Equation (24) fo Ψ GJ takes the fom ΩP Ψ GJ θ P θ Ψ GJ 1 sin 2 θ φp φ Ψ GJ Ω 1 sin θ ( φp θ φ P θ P φφ P) = 0. (44) By diet substitution of Ψ ot GJ = Ω sin θ θp (45) it an be shown that the potential (45) is a solution of the equation (44), and satisfies the bounday onditions (32, 33). Substituting this potential into the expession fo the GJ eleti field (19) one gets E GJ ot = 1 (Ω ) B. (46) This esult was obtained by Goldeih and Julian (1969), and Mestel (1971).

6 6 A.N.Timokhin et al. 3 SMALL OSCILLATIONS OF MAGNETIZED NEUTRON STAR WITH DIPOLE MAGNETIC FIELD. TOROIDAL OSCILLATIONS 3.1 Geneal Fomulas Let us onside tooidal osillations of a NS with dipola magneti field. Fo simpliity we assume that the mode axis is aligned with the dipole moment. We an do this without loss of geneality beause any osillation mode with mode axis not aligned with the dipola moment µ an be epesented by seies of osillation modes with mode axis paallel to µ. In this ase the unpetubed magneti field is B = B 0 ( ) 3 os θ e B0 ( ) 3 sin θ eθ, (47) whee e and e θ ae unit oodinate vetos. This field is desibed by the sala funtion P dip 0 aoding to the fomulas (A7), (A8) : P dip 0 = B03 os θ. 2 The fist ode field petubation is desibed by the funtion δp dip. The time deivative of this funtion fo an osillation mode (l, m), aoding fomulas (27), (B7), is ( ) tδp dip lm = 1 2 B0w φ Y lm, (49) l(l + 1) whee w e iωt W. Equation (37) fo the potential Ψ lm GJ of tooidal mode (l, m) takes the fom 2os θ Ψ lm GJ + 1 ( ) sin θ θψ lm GJ m2 B 0w 3 Ylm = 0, (50) l(l + 1) whee Ψ lm GJ oesponds to one exited osillation mode (l, m), and fo geneal mode Ψ GJ = lm Ψlm GJ. The haateistis of the equation (50) ae t = C 0 (48) Ψ lm GJ m2 B 0w l(l + 1) ( ) 2 sin 4 θ sin θ ( φ ) 1/2 = C 1 = C 2 (51) Ylm dθ = sin 5 C3 θ The integal of the equation (50) is an abitay funtion of onstants C 0, C 1, C 2, C 3 ϕ(c 0, C 1, C 2, C 3) = 0, Expessing Ψ lm GJ we have fo the geneal solution of equation (50) ( ) Ψ lm GJ = m2 B 0w 2 ( ( ) ) 1/2 sin 4 Ylm θ l(l + 1) sin 5 θ dθ + Φ lm sin θ, φ, t (52), (53) whee Φ is an abitay funtion. In ode fo Ψ lm GJ epesented by the expession (53) to be a G-J potential it must satisfy the bounday onditions on the sufae of the NS (35) and (36), whih in the ase of tooidal osillations of a NS with dipola magneti field take the fom Ψ lm GJ = = B0w Ψ lm GJ θ=0,= = 0. ) (os θ θ Y lm m2 Y lm dθ + e iωt F(φ), (54) l(l + 1) sin θ Substituting Ψ lm GJ given by the expession (53) into the bounday onditions (54), (55) we get the bounday ondition fo the funtion Φ lm : ( ) ) 1/2 ) Φ (sin lm θ, φ, t = (os B0w θ θ Y lm m2 Y lm dθ l(l + 1) sin θ = B0w m 2 Ylm l(l + 1) sin4 θ sin 5 θ dθ + ( dθ) + B0w m 2 sin 4 Ylm θ + l(l + 1) sin 5 θ θ=0 (55)

7 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 7 + B0w [ ( ) os θ θ Y lm m2 Y lm dθ] l(l + 1) sin θ θ=0, (56) whee the two last tems ae added in ode to satisfy the seond bounday ondition (55). The funtion Φ lm depends on θ in the ombination ( ) 1/2 sin θ. To get Φ lm one has to expess the ight hand side of equation (56) in tems of sin θ and to eplae sin θ by ( ) 1/2 sin θ. Substituting the funtion Φ lm into the expession (53) we get the potential Ψ lm GJ fo small tooidal osillations of a NS with a dipola magneti field fo the mode (l, m). Using the fomula (23) we get the Goldeih-Julian hage density fo that mode. We have developed a set of pogams on the ompute algeba language MATHEMATICA 3.0 (Wolfam 1996) fo alulating the analytial expessions of ρ lm GJ aoding to the algoithm desibed in this setion. These pogams wee tested fo some (l, m) by ompaing esults obtained by hand and by ompute, and also by heking the ondition E GJ B = Main Results Goldeih-Julian hage density In the solution fo the whole egion the expession (56) is used, whee on the ight hand side os θ should be eplaed by ( 1 ( ) sin 2 θ ) 1/2 fo 0 θ < π (fist hemisphee) and by ( 1 ( ) 2 sin 2 θ ) 1/2 fo π θ π (seond hemisphee). 2 Beause of this thee ae two diffeent expessions fo Ψ GJ fo both hemisphees. If these expessions give diffeent esults fo θ = π fo > 2 thee is a disontinuity in the funtion ΨGJ at the equatoial plane, and onsequently ρgj(, π, φ) beomes 2 infinite. In suh situation it is impossible to anel the aeleating eleti field in the whole magnetosphee of NS without pesene of a stong eleti uent (42) even in a egion of losed magneti field lines. Examples of suh osillation modes ae modes (l, m) = (1, 1), (2,0). The oesponding veloity fields ae given in Fig.1. The dependene of the potential Ψ lm GJ fo these modes along dipola magneti field lines is shown in Fig.2,4. In the following we nomalize the G-J potential Ψ GJ and hage density ρ GJ by B 0W / and B 0W/ espetively. The jump of Ψ GJ(,π/2, φ) deeases, as the angle at whih oesponding magneti field line inteset the NS s sufae is ineasing. Ψ GJ is a ontinuous funtion of θ on the sufae of the NS, nevetheless GJ hage density diveges on the equato even on the sufae of NS (see Fig.3, 5). In these figues the hage density ρ lm GJ = is shown in spheial oodinates ( ρ lm GJ =,θ, φ). Hee the adial oodinate epesents absolute value of the GJ hage density as funtion of the pola angle θ and the azimuthal angle φ. Analytial expessions fo the GJ hage density fo the disussed modes ae given in Appendix C. At the equatoial plane ρ lm GJ is infinite. The mode (2,1) is a epesentative of anothe lass of osillation modes, in whih Ψ lm GJ is a ontinuous funtion of θ (see Fig.6) and ρ lm GJ is finite eveywhee. The veloity field and ρ 21 GJ ae shown in Figs.1, 7. Fo suh modes it is possible to anel longitudinal eleti field without geneating stong uent along magneti field lines. Conside now the diffeene between two lasses of tooidal osillation modes fo a NS with a dipole magneti field. The bounday ondition fo Ψ lm GJ is popotional to the veto podut of v and B. The axis z is dieted along the dipole magneti momentum. If the veloity field is symmeti elative to the equatoial plane (l m is an odd numbe), then the bounday ondition fo Ψ lm GJ is also symmeti elatively to the equatoial plane. Hene, afte expession of all tigonometi funtions in (56) though os θ and sin θ, os θ appeas in the ight hand side of (56) only in even powes. Consequently the funtion Φ lm will be ontinuous fo and ρ lm GJ will be finite. If veloity field is antisymmeti elative to the equatoial plane (l m is an even numbe o zeo), then os θ is ontained in the ight side of (56) in odd powes, and the funtion Φ lm may have disontinuity in the plane θ = π, and 2 ρlm GJ fo suh modes beomes infinite at the equatoial plane. An example of osillation mode with veloity field antisymmeti elative to the equatoial plane and fo whih Ψ lm GJ is ontinuous funtion of θ is the mode (3, 3). Gaphis fo the potential Ψ 33 GJ, Goldeih-Julian hage density ρ 33 GJ ae shown in Figs. 8, 9. So, fo a dipola magneti field the longitudinal eleti field geneated by the osillation modes with odd (l m) an be aneled by putting of haged patiles in the magnetosphee. Among the modes with even o zeo (l m) thee ae some, fo whih longitudinal eleti field annot be aneled by putting of haged patiles without pesene of stong eleti uent along eah magneti field line. The Goldeih-Julian hage density gives only a haateistial hage density in the magnetosphee. The patile density an be obtained by solving full system of MHD equations, but this has not been done even fo an aligned otato with a dipola magneti field (Mihel 1991). The ase of nonadial osillation is muh moe ompliated, beause of the absene of stationaity and axial symmety. Fo osillation modes with odd (l m) thee an be hage sepaation in the nea zone beause the foot points of magneti field lines have the same sign of the potential, i.e fo suh modes haged patiles annot flow fom one hemisphee to anothe. Fo these modes ondition (40) is automatially satisfied (see end of setion 2.2). Evidently fo modes with even o zeo (l m) thee an not be a hage-sepaated solution in the nea zone.

8 8 A.N.Timokhin et al Enegy Losses In pevious woks the alulation of eletomagneti enegy losses of an osillating neuton sta inludes only adiation of eletomagneti waves in the vauum. If thee is plasma in the magnetosphee, then the enegy will be lost by a tansfomation of the osillation enegy into kineti enegy of an outflowing plasma, as it was poposed fo pulsa enegy losses by Goldeih and Julian (1969). They obtained the enegy losses of a otating aligned NS though outflow of the haged patiles fom a egion of open field lines, and found it to be equal to the loss in vauum though adiation of eletomagneti waves by a pependiula otato. In this pape we onside eleti and magneti fields in the nea zone, hene we an not expliitly show the existene of plasma outflow as in the aligned otato by Goldeih and Julian, but a qualitatively the pitue is the following. Eleti uent aising fom the motion of the haged patiles in the magnetosphee geneate a magneti field. In the egion whee the value of this field beomes lage than the unpetubed NS magneti field the field lines beome open, on the Alfveni sufae. Plasma esaping fom the egion of open field lines podues an eletomagnetially diven stella wind. The wind auses a net eleti uent whih loses at infinity. This uent flows beneath the stella sufae between positive and negative emission egions. Beause it must oss magneti field lines thee, it exets a baking toque on the osillating NS and edues the osillational enegy. This pitue is simila to the one poposed by Rudeman & Sutheland (1975) fo pulsas, but now we must get the bounday of the open field line egion self-onsistently. In ou ase we an not appoximate the losed field line egion by a zone of ootation, athe we should self-onsistently detemine the last losed field line as the last field line lying inside the Alfven sufae. Conside egion nea the pole. Let us denote as θ 0 a pola angle by whih the last losed field line intesets the sufae of the NS. Maximal value of eleti field indued by NS osillation is of the ode of magnitude E B V = B ξω. Fom fomulas (57) and (46) the atio of eleti fields and, onsequently, eleti potentials and Goldeih-Julian hage densities fo otating and osillation neuton sta is Ω : ω ξ. So, fo usually assumed (MDemott et al. 1988) paametes eleti fields geneated by osillations of a NS ae simila to the ones of 1 se pulsas. Simila to pulsas above a pola ap an aeleation zone and a zone of e + e pai geneation above it will be builded. Eleton-positon pais anel aeleating eleti field, i.e. motion of haged patile above the aeleating zone will be not influened by eleti filed geneated by stella osillations. In the egion of open field lines, whee patiles esape fom the neuton sta and fom elativisti wind, the aveage time of ossing of the aeleation zone by the patile is muh less then the osillation peiod (see Disussion), so patiles patially do not etun bak to the sta. Hene, aveaged ove the osillation peiod, the enegy loss though the outflow of plasma fom open field line egion is ǫ lm 1 T T 0 dt 2π θ0 0 0 j lm (, θ, φ) A lm (θ, φ) 2 sin θ dθdφ, (58) whee A(θ,φ) is the wok done by the eleti field to move a unit hage to the point with oodinate (, θ, φ): A(θ,φ) θ 0 E GJ θ (, θ, φ) dθ. Beause of easons disussed at the end of setion 1 hage density in the inne pats of the magnetosphee must be appoximately equal to the Goldeih-Julian hage density ρ GJ. In the egion of plasma outflow haged patiles ae steaming along magneti field lines in the same dietion with ultaelativisti veloities. Chage density of outflowing plasma is equal to the Goldeih-Julian hage density. Hene, outflowing ultaelativisti eleton-positon plasma builds a net uent density (57) (59) j ρ GJ. (60) The uent density an diffe fom the values given by expession (60) due to diffeene in aveaged veloities of eleton and positon omponents of the plasma (disussion on this topi see in Lyubaskii (1992)), but on the ode of magnitude expession (60) gives good estimation fo the uent density (Rudeman & Sutheland 1975). Hene, fo the open field line egion ondition (40) is satisfied and in expession (60) we an use ρ GJ obtained by solving of equation (37) fo any osillation mode. The last losed field line is the line fo whih the kineti enegy density of outflowing plasma on equato (at the point (R l, π/2, φ)) beomes equal to the oesponding enegy density of the NS magneti field: ǫ lm (θ 0) 4πRl 2 B2 8π θ= π 2,=R l. Expessing the ight hand side of this equation in tems of angle θ 0, we get two equations (58), (61) fo a self-onsistent detemination of θ 0 and ǫ lm. (61)

9 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 9 We have solved these equations fo osillation modes (1,1), (2,0), and (3,0). Osillation peiods wee taken fom MDemott et al. (1988). The atio of enegy losses to the vauum ones as a funtion of the amplitude of dimensionless tansvesal displaement ξ fo modes (1,1), (2,0), and (3,0) ae given in Table 1. One an see that enegy losses though outflowing plasma ae muh less than in vauum fo small values of ξ. On the othe hand, fo osillation modes l 2, fo displaement amplitudes ξ lage than a value ξ, enegy losses exeed the vauum losses (see Fig. 10). Fo some osillation modes fo whih θ 0 1, we an lineaize equations (58), (61) in θ 0 and get an analytial estimation of θ 0. Fo tooidal osillation mode (l, m) the veloity amplitude V nea the pole is of the ode { θ m 1, m 0 V W (62) θ, m = 0 The eleti field is of the ode E GJ B V/. Fom fomula (59) we have { A BW θ m, m 0 θ 2, m = 0 The Goldeih-Julian hage density is ρ GJ BW θm Substituting expessions (63) and (64) into equations (60), (59) and (58) we get the enegy loss though the plasma outflow along open field lines: ( ) W 2 { θ ǫ plasm B m+2 0, m 0 θ0, 4 (65) m = 0 The angle θ 0 is detemined fom the equation (61). Afte substitution of ǫ plasm into equation (61) we get ( ) 1 W 3 m, m 0 θ 0 (66) ( )1 W 2, m = 0 We see that the angle θ 0 is small only fo modes with m < 3. This is beause of the small value of V in the pola egion fo modes with m 3. Fo suh modes a linea analysis is not possible and the ase of lage θ 0 needs additional investigation, whih will be given elsewhee. The enegy loss though plasma outflow is ( ) ξ m ω 3 m, m = 1,2 ǫ plasm B0 2 2 (67) ( ) ω 4 ξ 4, m = 0 The vauum enegy loss fo the osillation mode (l, m), aoding to MDemott et al. (1984) is ( ) ω 8, l = 1, m = 0 ǫ va B0 2 2 ξ 2 (68) ( ) ω 2l+2, l = 1, m = 1; l 2, m = 0... l Fom these fomulas we get the osillation amplitude ξ fo whih the plasma enegy loss is lage than the vauum loss ( ) ω 2, l = 1, m = 0 ξ > ξ ( ) ω l 1, l 2, m = 0; l 1, m = 1 (69) ( ) ω l/3 1, l 2, m = 2 Fom fomula (69) it follows that fo modes (l = 1, m = 0), (l 2, m = 0, 1) and (l 4, m = 2) nonvauum enegy loss exeeds the vauum loss even fo small displaement amplitudes ξ, beause fo these modes ξ 1. (63) (64) Aoding to MDemott et al. (1988) the dependene of the eigenfequeny on l is vey weak, at least fo small l, so we assumed the same osillation peiods fo modes (2,0) and (3,0)

10 10 A.N.Timokhin et al. Mode ǫ pl /ǫ va R lm op /λ (1, 1) ξ 2 0.9/ ξ (2, 0) 77 ξ 2 τ 2 0.2/ ξ (3, 0) ξ 2 τ / ξ Table 1. Leading tems in expansion ove ξ of atios of eletomagneti enegy loss though plasma outflow to the vauum one epsilon pl /ǫ va and distane in the equatoial plane of the last losed field line to the wave length R lm op /λ fo tooidal osillation modes (1, 1), (2, 0) and (3, 0). τ is peiod of osillations in milliseonds. Eletomagneti waves adiated by the osillating NS ae seened by the plasma. Fo small ξ the nonvauum enegy loss aoding to (69) is smalle than the vauum loss beause at the GJ density of patiles the aeleated plasma annot esape fom the magnetosphee. The atio of the equatoial adius fo the last losed field line, to the wave length λ ae given in Table 1. Fo mode (3, 0) this dependene is shown in Fig. 11. We see that the last losed field line fo small values of displaement amplitude lies deep inside the wave zone. In the egion limited by the last losed field lines and a sphee of adius λ (fomal bounday of the wave zone) in this situation onfined plasma waves ae exited suh as Alfven o magnetosoni wave by the osillation of the NS, see Fig. 12. This will lead to aumulation of enegy in this egion in the fom of plasma wave enegy. When the density of this enegy beomes lage than the magneti enegy, field lines beomes open. These poesses will lead to a deease of the equatoial adius of the last losed field lines, and hene, a inease of the enegy loss. A seond poess whih ould lead to lage enegy losses than obtained by the solution of the equation (61), is otation of the neuton sta. In this ase the equatoial adius of the last losed field line is the adius of the light ylinde R L = Ω/, and sufae aea fo plasma outflow may be onsideably ineased. It is not lea, howeve, how the inteation between otation an osillation will influene the total enegy losses. 4 DISCUSSION We have developed a fomalism desibing the magnetosphee of an osillating neuton sta in the same way as in the ase of otation. With this fomalism it is possible to apply theoetial models developed fo pulsas to osillating neuton stas. Fo example, one an investigate the aeleation mehanism poposed by Shalemann, Aons & Fawley (1978) o by Rudeman & Sutheland (1975) fo osillating neuton sta. The inetial fame dagging mehanism poposed by Muslimov & Tsygan (1992) does not wok hee, beause we onside a nonotating sta. Chaateisti time of pai asade fomation is estimated as h PPF/, whee h PPF is the hight of pai fomation font above the sufae. In Rudemann-Sutheland model minimum value of h PPF is the same as fo a NS otating with the angula fequeny Ω eff = ω ξ, and fo usually assumed paametes of NS osillation is of ode m. In Aons-Shalemann like model the hight of pai fomation font an be smalle then in a pulsa otating with the angula fequeny Ω eff (seveal km), beause in geneal fo nonadial osillation atio ρ GJ/B along magneti field lines inease faste then in magnetosphee of a otating NS, what leads to lage aeleating eleti fields. Beause the haateisti time of pai asade fomation is muh shote than the peiod of osillation, one an onside the pola ap aeleation zone in any given moment of time as stationay. Poblem of etun uent egion fo the ase of osillating NS as in the ase of pulsa emains open. Simila to pulsas etun uent an flow along the last losed magneti field line. The enegy losses of an osillating neuton sta ae stongly affeted by plasma pesent in its magnetosphee. This affets all pevious alulations of the eletomagneti damping ate of NS osillation. Ou alulations give low bounday on eletomagneti enegy losses of an osillating NS beause of easons disussed at the end of setion The enegy of sta osillation depend on displaement amplitude as ξ 2. Enegy losses though plasma outflow in ontast to enegy losses though adiation of eletomagneti waves in vauum have anothe dependene on the displaement amplitude. Fo estimation of osillation damping time one an use values given in MDemott et al. (1988), Table 6, multiplying by funtions fom Table 1 of this pape fo given displaement amplitude. Hee we have estited ouself to a simple ase of tooidal osillations. Staquakes should exite a whole spetum of osillations, inluding p- and g- modes. The magnetosphee stutue podued by these modes will be onsideed in a futue pape. The authos wish to thank A. I. Tsygan fo helpful disussions. This wok has been suppoted by Russian Foundation of Basi Reseah (RFFI) unde gant

11 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 11 Figue 1. Veloity field on a sphee fo tooidal modes (1,1), (2,0), (2,1) and (3,3) ae shown as pojetion on the meidional plane φ = 115. Ciles oesponding to the longitudes 0 (ight), 90 (left) and latitude 0 ae shown by dashed lines.

12 12 A.N.Timokhin et al. Figue 2. The potential Ψ 11 GJ along a dipola magneti field line as a funtion of the pola angle θ is shown fo 5 field lines with azimuthal angle φ = π/2, fo t = π n/(2ω), at whih the maximal absolute values of the potential ae eahed. The values of θ of the left and ight ends of the lines 1-5 detemine the pola angle at whih oesponding magneti field line osses the sufae of the neuton sta. The disontinuity is shown by the jump of ψ oesponding to the same field line in the opposite hemisphees. Figue 3. left: The hage density ρ 11 GJ = is shown in a spheial oodinate system ( ρ lm GJ =, θ, φ) fo t = π n/(2ω). The adial oodinate epesents the absolute value of the GJ hage density as funtion of the pola angle θ and the azimuthal angle φ. Positive values of hage density ae shown by a gay sufae and negative ones by white. The dipole magneti moment µ is dieted upwads along the z axis. Azimuthal angle is ounted fom the x-axis, as in Fig. 1. Note that fo θ = π/2, ρ 11 GJ is infinite. ight: Coss-setion of sufaes fom the left figue by the plane φ = π/2, whee the maximal absolute values of the hage density ρ 11 GJ = ae eahed. Ciles oespond to the values 0.01, 0.02, and Positive values ae shown solid, negative ones dashed. Note that fo θ = π/2, ρ 11 GJ is infinite.

13 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 13 Figue 4. The potential Ψ 20 GJ along a dipola magneti field line as a funtion of the pola angle θ is shown fo 5 field lines fo t = π n/(2ω), at whih the maximal absolute values of the potential ae eahed. The values of θ of the left and ight ends of the lines 1-5 detemine the pola angle at whih oesponding magneti field line osses the sufae of the neuton sta. The disontinuity is shown by the jump of ψ oesponding to the same field line in the opposite hemisphees. Figue 5. left: The hage density ρ 20 GJ = is shown in a spheial oodinate system ( ρ lm GJ =, θ, φ) fo t = π n/(2ω). The adial oodinate epesents the absolute value of the GJ hage density as funtion of the pola angle θ and the azimuthal angle φ. Positive values of hage density ae shown by a gay sufae and negative ones by white. The dipole magneti moment µ is dieted upwads along the z axis. Azimuthal angle is ounted fom the x-axis, as in Fig. 1. Note that fo θ = π/2, ρ 20 GJ is infinite. ight: Coss-setion of sufaes fom the left figue by a meidional plane. Ciles oespond to the values 0.1, 0.2, and 0.3. Positive values ae shown solid, negative ones dashed. Note that fo θ = π/2, ρ 20 GJ is infinite.

14 14 A.N.Timokhin et al. Figue 6. The potential Ψ 21 GJ along a dipola magneti field line as a funtion of the pola angle θ is shown fo 5 field lines with azimuthal angle φ = π/2, fo t = π n/(2ω), at whih the maximal absolute values of the potential ae eahed. The values of θ of the left and ight ends of the lines 1-5 detemine the pola angle at whih oesponding magneti field line osses the sufae of the neuton sta. Figue 7. left: The hage density ρ 21 GJ = is shown in a spheial oodinate system ( ρ lm GJ =, θ, φ) fo t = π n/(2ω). The adial oodinate epesents the absolute value of the GJ hage density as funtion of the pola angle θ and the azimuthal angle φ. Positive values of hage density ae shown by a gay sufae and negative ones by white. The dipole magneti moment µ is dieted upwads along the z axis. Azimuthal angle is ounted fom the x-axis, as in Fig. 1. ight: Coss-setion of sufaes fom the left figue by the plane φ = π/2, whee the maximal absolute values of the hage density ρ 21 GJ = ae eahed. Ciles oespond to the values 0.05, 0.1, and Positive values ae shown solid, negative ones dashed.

15 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 15 Figue 8. The potential Ψ 33 GJ along a dipola magneti field line as a funtion of the pola angle θ is shown fo 5 field lines with azimuthal angle φ = π/6, fo t = π n/(2ω), at whih the maximal absolute values of the potential ae eahed. The values of θ of the left and ight ends of the lines 1-5 detemine the pola angle at whih oesponding magneti field line osses the sufae of the neuton sta. Figue 9. left: The hage density ρ 33 GJ = is shown in a spheial oodinate system ( ρ lm GJ =, θ, φ) fo t = π n/(2ω). The adial oodinate epesents the absolute value of the GJ hage density as funtion of the pola angle θ and the azimuthal angle φ. Positive values of hage density ae shown by a gay sufae and negative ones by white. The dipole magneti moment µ is dieted upwads along the z axis. Azimuthal angle is ounted fom the x-axis, as in Fig. 1. ight: Coss-setion of sufaes fom the left figue by the plane φ = π/6, whee the maximal absolute values of the hage density ρ 33 GJ = ae eahed. Ciles oespond to the values 0.04, 0.08, and Positive values ae shown solid, negative ones dashed.

16 16 A.N.Timokhin et al. Figue 10. Ratio of the enegy loss though outflow of plasma fom the egion of open field lines to the loss in vauum is shown as a funtion of the dimensionless displaement amplitude ξ fo the mode (3,0). The solid line oesponds the osillation peiod T = 19ms (mode 3 t 0 in notations of MDemott et al. (1988)), the dotted line T = 1.7ms (mode 3 t 1 ), the dashed line T = 1.08ms (mode 3 t 2 ); see note on the page 9. Figue 11. Ratio of the equatoial adius of the last losed field line to the wavelength λ is shown as a funtion of the dimensionless displaement amplitude ξ fo the mode (3,0). This atio does not depend on the osillation fequeny ω.

17 yy yyy Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 17 last losed field line plasma outflow plasma waves plasma waves bounday of the wave zone plasma outflow Figue 12. Stutue of the magnetosphee of an osillating neuton sta (abitay saled). See explanations in text.

18 18 A.N.Timokhin et al. APPENDIX APPENDIX A: CALCULATION OF P AND δp Substituting B fom fomula (15) using expansion (25) into the equation (5) we have (P 0 + δp)e = 0. (A1) Beause P does not depend on t, equation (A1) is equivalent to the following two equations: (P 0e ) = 0 (δpe ) = 0. (A2) (A3) In spheial oodinate equation (A3) is (δpe ) = e θ 1 sin θ Substituting expansion of δp (26) we get ( 1 2 φ ΩδP + φ δp) + e φ 1 ( 1 2 θ ΩδP + θ δp) = 0. (A4) 1 2 l(l + 1)δ p lm(, t) + δ p lm (, t) = 0. (A5) The solution of equation (A5) whih vanishes at infinity is ( ) l δ p lm (,t) = δplm (t). The time-dependent oeffiients p lm (t) have to be detemined fom the bounday onditions. Substituting expession (A6) into expansion (26) and taking the deivative with espet to t we get expession (27). Similaly, fo the oeffiients P 0 lm in the expansion of the funtion P 0 we have P 0 lm () = p lm ( ) l Substituting the expansion of the funtion P 0 in spheial hamonis into the bounday ondition (6), multiplying this expession by Ylm, and integating it ove solid angle 4π, we have (ompae with fomula (16) in Muslimov & Tsygan (1986)) p lm = 2 B 0 Ylm dω, (A8) l(l + 1) 4π whee Ylm is invese spheial hamoni Y lm Y l m dω = δ ll δ mm and dω = sin θdθ dφ is the solid angle (A6) (A7) APPENDIX B: CALCULATION OF δp lm (t) Diffeentiating equation (32) with espet to φ and equation (33) with espet to θ, equating the esults and multiplying by 1/sin θ we have [ 1 sin θ θ(sin θ θ ( tp)) + 1 { sin 2 θ φφ( tp)] = 1 [ 1 = ΩP sin θ θ(sin θ v θ ) + 1 ] sin θ φ v φ [ 1 v sin θ θ(sin θ θ P) + 1 ] sin 2 θ φφp [ 1 θ ΩP v θ + 1 ] sin θ φ ΩP v φ [ θ P θ v + 1 ]} sin 2 θ φp φ v (B1) = Simplifying equation (B1) by witing P as a sum of the two tems (25), we get an expession fo time-dependent pat δp to fist ode in ξ [ ( ) ( 1 t ΩδP = = v ΩP 0 + v θ 1 1 θ ΩP 0 + v φ sin θ φ ΩP 0 ΩP 0 sin θ θ(sinθ v θ ) + 1 ) sin θ φ v φ ( θ P 0 θ v + 1 )] sin 2 θ φp 0 φ v (B2) = o in veto notation t ΩδP = = { v ( ΩP 0) + ΩP 0( v ) + 2 ( ( P 0) ) v } =, (B3)

19 Magnetosphee of Osillating Neuton Sta. Nonvauum Teatment. 19 whee v is the tangential pat of the veloity, v = e θ v θ + e φ v φ, and is the angula pat of the gadient: = e θ 1/ θ + e φ 1/( sin θ) φ. Substituting the expansion of δp in spheial hamonis (27), multiplying the esult by Ylm and integating it ove the solid angle 4π, we get an expession fo the oeffiients of tδp(t) (34) in the expansion of δp in spheial hamonis: 1 [ ] tδp lm (t) = dω Ylm v ( ΩP 0) + ΩP 0( v ) + 2 ( ( P 0) ) v. l(l + 1) = 4π Next we give expliit expessions fo oeffiients δp lm (t) fo tooidal and spheoidal osillations modes (see Unno et al. (1979)). Spheoidal sepaation of vaiables gives the following expessions fo omponents of the osillation veloity in spheial oodinates v = e iωt U() Y lm (θ, φ), v θ = e iωt V () θ Y lm (θ, φ), v φ = e iωt 1 V () sin θ φy lm (θ, φ), (B4) whee U and V ae adial and tansvesal veloity amplitude espetively. Substituting the veloity omponents (B4) into fomula (27), we get the oeffiients tδp lm (t) fo spheoidal osillations [ ( 1 tδp(t) = dω Ylm e iωt U Y lm ΩP 0 + θ Y lm θ P ) l(l + 1) sin 2 θ φy lm φ P 0 + e iωt V 4π ( l(l + 1)Y lm ΩP 0 + θ Y lm θ ΩP sin 2 θ φy lm φ ΩP 0 )] = The osillation veloity omponents fo tooidal osillations ae v = 0, v θ = e iωt W() 1 sin θ φy lm (θ, φ), v φ = e iωt W() θ Y lm (θ, φ), (B6) whee W is tansvesal veloity amplitude. With the use of fomula (B6), the oeffiients in the expansion of tδp fo tooidal osillations ae tδp lm = 1 [ dω Ylm e iωt W 1 l(l + 1) sin θ ( θy lm φ ΩP 0 φ Y lm θ ΩP 0)]. (B7) = 4π (B5) APPENDIX C: EXPRESSIONS FOR GOLDREICH-JULIAN CHARGE DENSITY FOR MODES (1,1), (2,0), (2,1), (3,3) ρ 11 GJ fo 0 θ < π e iωt B W sin(θ) sin(φ) ( 2 32 π ( 6 ) os(2 θ) os(2 θ) (os(θ) + 3 os(3 θ)) os(2 θ) 3 2 os(4 θ)) (C1) ρ 11 GJ fo π 2 < θ π e iωt B W sin(θ) sin(φ) ( 2 32 π ( 6 ) os(2 θ) os(2 θ) (os(θ) + 3 os(3 θ)) os(2 θ) os(4 θ)) (C2) ρ 20 GJ fo 0 θ < π 2 e 3 5 iωt 2 2 B W ( ( 4 ) os(2 θ) + 15 os(4 θ)) 128 π os(2 θ) (C3) ρ 20 GJ fo π 2 < θ π e 3 5 iωt 2 2 B W ( ( 4 ) os(2 θ) + 15 os(4 θ)) 128 π os(2 θ) (C4) ρ 21 GJ fo 0 θ π

20 20 A.N.Timokhin et al. e iωt B W ( os(2 θ)) sin(θ) sin(φ) 32 π (C5) ρ 33 GJ fo 0 θ π e iωt B W os(θ) sin(θ) 3 sin(3 φ) 64 π (C6) REFERENCES Bisnovatyi-Kogan G. S., 1995, ApJS 97, 185 Bisnovatyi-Kogan G. S., Imshennik V. S., Nadiozhin D. K., Chehetkin V. M., 1975, Ap&SS, 35, 23 Blaes O., Blandfod R., Goldeih P., Madau P., 1989, ApJ, 343, 839 Boiakoff V., 1976, ApJ, 208, L43 Bandt L., 1947, Veto and Tenso Analysis. John Wiley & Sons, In., New Yok Caoll B. W., Zweibel E. G., Hansen C. J., MDemott P. N., Savedoff M. P., Thomas J. N., Van Hon H. M., 1986, ApJ, 305, 767 Daughety J. K., Hading A. K., 1982, ApJ, 252, 337 Ding K. Y., Cheng K. S., 1997, MNRAS, 287, 671 Dunan R. C., 1998, ApJ, 498, L45 Elsgolts L. E., 1965, Diffeentsialnye uavneniya i vaiatsionnoe ishislenie. Nauka, Moskva (in Russian) Fatuzzo M., Melia F., 1993, ApJ, 407, 680 Goldeih P., Julian W. H., 1969, ApJ, 157, 869 Hankins T. H., 1996, in Astonomial Soiety of the Paifi Confeene Seies, Vol. 105, IAU Colloquium 160, Pulsas: Poblems and Pogess ed. S. Johnston, M. A. Walke & M. Bailes, 197 Huley K., et al., 1994, Nat, 372, 652 Lyne A. G., 1996, in Astonomial Soiety of the Paifi Confeene Seies, Vol. 105, IAU Colloquium 160, Pulsas: Poblems and Pogess ed. S. Johnston, M. A. Walke & M. Bailes, 73 Lyubaskii Yu. E., 1992, A&A, 261, 544 MDemott P. N., Van Hon H. M., Hansen C. J., 1988, ApJ, 325, 725 MDemott P. N., Savedoff M. P., Van Hon H. M., Zweibel E. G., Hansen C. J., 1984, ApJ, 281, 746 Mestel L., 1971, Natue Physial Siene, 233, 149 Mihel F. C., 1991, Theoy of neuton sta magnetosphees. Univesity of Chiago Pess, Chiago Muslimov A. G., Tsygan A. I., 1986, Ap&SS, 120, 27 Muslimov A. G., Tsygan A. I., 1992, MNRAS, 255, 61 Paini F., Rudeman M., 1974, Nat, 251, 399 Rudeman M. A., Sutheland P. G., 1975, ApJ, 196, 51 Shalemann E. T., Aons J., Fawley W. M., 1978, ApJ, 222, 297 Smith I. A., Epstein R. I., 1993, ApJ, 410, 315 Stohmaye T. E., 1991, ApJ, 372, 573 Stuok P. A., 1971, ApJ, 164, 529 Tümpe J., Beke W., 1998, Advanes in Spae Reseah, 21, 203 Tsygan A. I., 1975, A&A, 44, 21 Unno W., Osaki Y., Ando H., Shibahashi H., 1979, Nonadial Osillations of Stas. Univesity of Tokyo Pess, Tokyo Wolfam S., 1996, The Mathematia book, Thid Edition. Wolfam Media/Cambidge Univesity Pess, Champaign, Illinois