The plasma emission model of RBS1774
Key Words.:X-rays – stars: pulsars: individual RBS1774 – radiation mechanisms: non-thermal
In the present paper we construct a self-consistent theory, interpreting the observational properties of RBS1774. It is well known that the distribution function of relativistic particles is one-dimensional at the pulsar surface. However, cyclotron instability causes an appearance of transverse momenta of relativistic electrons, which as a result, start to radiate in the synchrotron regime. We study the process of the quasi-linear diffusion developed by means of the cyclotron instability on the basis of the Vlasov’s kinetic equation. This mechanism provides generation of measured optical and X-ray emission on the light cylinder lengthscales. A different approach of the synchrotron theory is considered, giving the spectral energy distribution that is in a good agreement with the XMM-Newton observational data. We also provide the possible explanation of the spectral feature at keV, in the framework of the model.
RBS1774 (1RXS J214303.7+065419) has been the most recent XDIN (Xray dim isolated neutron star) to be found (Zampieri et al., 2001). Its X-ray spectrum is well reproduced by an absorbed blackbody with a temperature eV and with a total column density of . Application of more sophisticated, and physically motivated models for the surface emission (atmospheric models) result in worse agreement with the data (Zane et al., 2005). According to Schwope et al. (2009), a fit to the X-ray spectra extracted from RGS spectrographs onboard XMM-Newton yields that the best result is obtained when the two-temperature blackbody model is used. But the same model applied to the X-ray spectra extracted from three EPIC detectors does not improve the fit compared to the simple blackbody model. However, the formation of a non-uniform distribution of the surface temperature is more likely artificial and needs to be examined by convincing theory.
Alternatively, the observational properties of RBS1774 can be explained in the framework of the plasma emission model first developed by Machabeli & Usov (1979) and Lominadze et al. (1983). According to these works, in the electron-positron plasma of a pulsar magnetosphere the waves excited by the cyclotron resonance interact with particles, leading to the appearance of pitch angles, which obviously causes synchrotron radiation. We suppose that the X-ray emission from this object is generated by the synchrotron mechanism. According to the standard theory of the synchrotron emission (Bekefi & Barrett, 1977; Ginzburg, 1981) the typical synchrotron spectrum is a power-law, when the present model suggests different spectral distribution. The main reason for this is that we take into account the mechanism of creation of the pitch angles, consequently restricting their values. Contrary to this, in the standard theory of the synchrotron radiation, it is assumed that along the line of sight the magnetic field is chaotic, leading to the broad interval (from to ) of the pitch angles. The present model gives successful fit for the observed X-ray spectrum, when the originally excited cyclotron modes enter the same domain as the measured optical emission of RBS1774. We suppose that the observed spectral feature at keV in the X-ray spectrum of RBS1774 is caused by wave damping process developed near the light cylinder due to the cyclotron instability.
In this paper, we describe the emission model (Sec. 2), derive theoretical X-ray spectrum of RBS1774 and fit with XMM-Newton observations (Sec. 3), explain the possible nature of the observed spectral feature at keV (Sec.4), and discuss our results (Sec. 5).
2 Emission model
The distribution function of relativistic particles is one dimensional at the pulsar surface, because any transverse momenta () of relativistic electrons are lost in a very short time(s) via synchrotron emission in very strong magnetic fields. For typical pulsars the plasma consists of the following components: the bulk of plasma with an average Lorentz-factor , a tail on the distribution function with , and the primary beam with (see Fig. 1). However, plasma with an anisotropic distribution function becomes unstable, which can lead to a wave excitation in the pulsar magnetosphere. The generation of waves is possible during the further motion of the relativistic particles along the dipolar magnetic field lines if the condition of cyclotron resonance is fulfilled (Kazbegi et al., 1991):
where is the drift velocity of the particles due to curvature of the field lines, is the radius of curvature of the field lines and is the cyclotron frequency. During the generation of waves by resonant particles, one also has a simultaneous feedback of these waves on the electrons (Vedenov et al., 1961). This mechanism is described by quasi-linear diffusion, leading to the diffusion of particles as along as across the magnetic field lines. Therefore, resonant electrons acquire transverse momenta (pitch angles) and, as a result, start to radiate through the synchrotron mechanism.
where is the force responsible for conserving the adiabatic invariant const, - is the radiation deceleration force produced by synchrotron emission, and is the reaction force of the curvature radiation. They can be written in the form:
where , and is the pitch angle.
Now let us compare the transverse components of the forces and . If we consider the case we will have :
where is the magnetic field at the star surface, is the star radius and is the distance from the pulsar. For the typical parameter values of pulsars . Then taking into account that the equation for the diffusion across the magnetic field can be written in the form (Chkheidze et al., 2010)
is the diffusion coefficient and is the density of electric energy in the waves.
The transversal quasi-linear diffusion increases the pitch-angle, whereas force F resists this process, leading to the stationary state (). Then the solution of Eq. (7) is
To evaluate , we use the quantity
where is the frequency of original waves, excited during the cyclotron resonance and can be estimated from Eq. (1) as follows . Consequently, we will get
The mean value of the pitch-angle (i. e. the assumption done at the beginning of our computations proves to be true). As a result of the appearance of the pitch angles, the synchrotron emission is generated.
3 X-ray spectrum
Let us consider the synchrotron emission of the set of electrons. The number of emitting particles in the elementary volume is , with momenta from the intervals and , and with the velocities that lie inside the solid angle near the direction of . If we write the parallel distribution function of the emitting particles as , then the emission flux of the set of electrons will be (Ginzburg, 1981)
where is the Stokes parameter, which is additive in this case, as the observed synchrotron radiation wavelength is much less than the value of - the average distance between particles, where is the density of plasma component electrons. The integral (12) is easily reduced to (see Ginzburg (1981))
Here keV is the photon energy of the maximum of synchrotron spectrum of a single electron and is a Macdonald function. After substituting the mean value of the pitch-angle in the above expression for , we get
For the primary beam electrons with the Lorentz factor the emitted photon energy keV comes in the energy domain of the observed X-ray emission of RBS1774. Thus we suppose that the measured X-ray spectrum is the result of the synchrotron emission of primary beam electrons (the resonance occurs on the right slope of the distribution function of beam electrons (see Fig. 1)), switched on as the result of acquirement of pitch angles by particles during the quasi-linear stage of the cyclotron instability.
To find the synchrotron flux in our case, we need to know the one-dimensional distribution function of the emitting particles . Let us multiply both sides of Eq. (2) on and integrate it over . Using Eqs. (3), (4), (5) and the following expressions for the diffusion coefficients (Chkheidze et al., 2010)
And also taking into account that the distribution function vanishes at the boundaries of integration, Eq. (2) reduces to
Let us estimate the contribution of different terms on the righthand side of Eq. (16). The estimations show that the first term is much bigger than two other terms. Consequently, for the primary-beam electrons instead of Eq. (16), one gets
Considering the quasi-stationary case we find
For , a magnetic field inhomogeneity does not affect the process of wave excitation. The equation that describes the cyclotron noise level, in this case, has the form (Lominadze et al., 1983)
is the growth rate of the instability. Here can be found from the resonance condition (1)
Combining Eqs. (17) and (19) one finds
which reduces to
Taking into account that for the initial moment the major contribution of the lefthand side of the Eq. (24) comes from , the corresponding expression writes as
The distribution function f is proportional to , then one should neglect in comparison with . Consequently, the above equation reduces to
As we can see the function drastically depends on the form of the initial distribution of the primary beam electrons. Here we assume that the initial energy distribution in the beam has a Gaussian shape
where - is the half width of the distribution function and is the density of primary beam electrons, equal to the Goldreich-Julian density (Goldreich & Julian, 1969). Since , this distribution is very close to -function. Consequently, the electron distribution can be taken as monoenergetic.
In this case for the energy density of the waves we get
The effective value of the pitch angle depends on as follows
According to our emission model, the observed radiation comes from a region near the light cylinder radius, where the magnetic field lines are practically straight and parallel to each other (Osmanov et al., 2009), therefore, electrons with efficiently emit in the observer’s direction.
Using expression (18), (28) and (29), and replacing the integration variable by , from Eq. (13) we will get
The energy of the beam electrons vary in a small interval. In this case the integral (30) can be approximately expressed by the following function
We performed a spectral analysis by fitting the model spectrum (Eq. (31)) absorbed by cold interstellar matter, with the combined data extracted of the three EPIC X-ray cameras of the XMM-Newton Telescope. The resulting and the amount of interstellar matter , which appears to be close to the total Galactic absorption in the source direction ( Dickey & Lockman (1990)). The spectral feature at keV that is mostly described as an absorption edge or line (Zane et al., 2005; Schwope et al., 2009) is also evident in our case from inspection of Fig. 2. We find that adding an absorption edge at keV improves the fit, leading to a reduced . The best-fitting energy of the edge is keV, and the optical depth is (see Fig. 3). The fitting results are listed in Table 1. We suppose that existence of the absorption feature in spectra of RBS1774 is caused by wave damping at photon energies keV, which takes place near the light cylinder.
4 Possible nature of the spectral feature
During the farther motion in the pulsar magnetosphere, the X-ray emission of RBS1774 that is generated on the light cylinder lengthscales, might come in the cyclotron damping range (Khechinashvili & Melikidze, 1997):
The condition for the development of the cyclotron instability may be easily derived for the small angles of propagation with respect to the magnetic field. Representing the dispersion of the waves as
and neglecting the drift term, the resonance condition (32) may then be written as
where is the angle between the wave vector and the magnetic field. Taking into account that one finds from Eq. (34) the frequency of damped waves
If we assume that the resonant particles are the primary beam electrons, then the estimation shows that on the light cylinder lengthscales keV.
According to the generally accepted point of view, the X-ray spectrum of RBS1774 is purely thermal and is best represented by a Planckian shape. A fit with a pure blackbody component absorbed by cold interstellar matter gives (Schwope et al., 2009). Including a Gaussian absorption line at keV (as the largest discrepansies between model and data are around keV) improves the fit (parameters are listed in Tab.1). However, the nature of this spectral feature is not fully clarified as yet. The most likely interpretation is that it is due to proton cyclotron resonance, which implies ultrastrong magnetic field of G (Zane et al., 2005; Rea et al., 2007). Although, the required strong magnetic field is inconstistent with timing measurements giving G (Kaplan & van Kerkwijk, 2009).
We are not about to reject the existing thermal emission models, but in present paper we propose an alternative explanation of the observed X-ray spectrum of RBS1774. It is supposed that the emission of this source is generated by the synchrotron mechanism. The distribution function of relativistic particles is one dimensional at the pulsar surface, but plasma with an anisotropic distribution function is unstable which can lead to wave excitation. The main mechanism of wave generation in plasmas of the pulsar magnetosphere is the cyclotron instability, which develops on the light cylinder lengthscales. During the quasi-linear stage of the instability, a diffusion of particles arises along and across the magnetic field lines. Therefore, plasma particles acquire transverse momenta and, as a result, the synchrotron mechanism is switched on. If the resonant particles are the primary beam electrons with their synchrotron emission enter the same energy domain as the measured X-ray spectrum of RBS1774.
We construct a self-consistent theory interpreting the observations of RBS1774. Differently from the standard theory of the synchrotron emission (Ginzburg, 1981), which only provides a power-law spectrum with the spectral index less than 1, our approach gives the possibility to obtain different spectral energy distributions. In the standard theory of the synchrotron emission, it is supposed that the observed radiation is collected from a large spacial region in various parts of which, the magnetic field is oriented randomly. Thus, it is supposed that along the line of sight the magnetic field directions are chaotic and when finding emission flux, Eq. (13) is averaged over all directions of the magnetic field (which means integration over varying from to ). In our case the emission comes from a region of the pulsar magnetosphere where the magnetic field lines are practically straight and parallel to each other. And differently from standard theory, we take into account the mechanism of creation of the pitch angles. Thus we obtain a certain distribution function of the emitting particles from their perpendicular momenta (see Eq. (10)), which restricts the possible values of the pitch angles. In the framework of the model we obtain the following theoretical X-ray spectrum of RBS1774 and perform a spectral analysis by fitting data from the three EPIC detectors simultaneously. The fit with a model spectrum absorbed by cold interstellar matter yields (see parameters in Tab.1).
During the farther motion in the pulsar magnetosphere, the X-ray emission of RBS1774 comes in the cyclotron damping range (see Eq. (32)). If we assume that damping happens on the left slope of the distribution function of primary beam electrons (see Fig. 1), then the photon energy of damped waves will be keV. Taking into account the shape of the distribution function of beam electrons, we interpret the large residuals around keV (see Fig. 2) as an absorption edge. Including an absorption edge improves the fit leading to a reduced . The best-fitting energy of the edge is keV, and the optical depth is (see Tab.1). However, adding an absorption edge to the model spectrum does not produce a statistically significant improvement of the fitting. According to Schwope et al. (2009) if one uses the RGS X-ray spectra of RBS1774 in place of EPIC spectra, the resulting is changed just marginally when a Gaussian absorption line is included at keV. Thus, we conclude that the nature of the feature at keV is uncertain and might be related to calibration uncertainties of the CCDs and the RGS at those very soft X-ray energies. The same can be told about a feature at keV (the large residuals around keV are evident from inspection of Fig. 2 and 3). A feature of possible similar nature was detected in EPIC-pn spectra of the much brighter prototypical object RXJ1856.4-3754 and classified as remaining calibration problem by Haberl (2007). Consequently, more data are necessary to finally prove or disprove the existence of those features.
The author is grateful to George Machabeli for valuable discussions and Axel Schwope for providing the X-ray data.
- email: email@example.com
- email: firstname.lastname@example.org
- Arons, J. 1981, In: Proc. Varenna Summer School and Workshop on Plasma Astrophysics, ESA, p.273
- Bekefi George & Barrett Alan H., 1977, Electromagnetic vibrations, waves and radiation, The MIT Press, Cambridge, Massachusetts and London, England
- Burwitz, V., Zavlin, V. E., Neuhäuser, R., Predehl, P., Trümper, J. Brinkman, A. C., 2001, A&A, 379, L35
- Burwitz, V., Haberl, F. Neuhäuser, R., Predehl, P., Trümper, J. Zavlin, V. E., 2003, A&A, 399, 1109
- Chkheidze, N., Machabeli, G., 2007, A&A, 471, 599
- Chkheidze, N., Lomiashvili, D., 2008, NewA, 13, 12
- Chkheidze, N., Machabeli, G., Osmanov, Z., 2010, ApJ, (submitted)
- Dickey, J., & Lockman, F. 1990, ARA&A, 28, 215
- Drake, J. J., Marshall, H. L., Dreizler, S., Freeman, P. E., Fruscione, A., et al. 2002, ApJ, 572, 996
- Ginzburg, V. L., ”Teoreticheskaia Fizika i Astrofizika” ,Nauka, Moskva 1981
- Goldreich, P., Julian, W. H., 1969, ApJ, 157, 869
- Haberl, F., 2007, Ap&SS, 308, 181
- Ho, W. C. G., Kaplan, D. L., Chang, P., Adelsberg, M., Potekhin, A. Y., 2007, MNRAS 375, 281H
- Ho, W. C. G., 2007, MNRAS 380, 71H
- Kaplan, D. L., van Kerkwijk, M. H., 2009, ApJ, 692, 62
- Kazbegi A. Z., Machabeli G. Z., Melikidze G. I., 1991, MNRAS, 253, 377
- Kazbegi A. Z., Machabeli G. Z., Melikidze G. I., Shukre C., 1996, A&A, 309, 515
- Khechinashvili, D. G., Melikidze, G. I., 1997, A&A, 320, L45
- Lai, D., Salpeter, E. E., 1997, ApJ, 491, 270
- Lai, D., 2001, Rev. of Mod. Phys. 73, 629
- Lomiashvili D., Machabeli G., & Malov I., 2006, ApJ, 637, 1010
- Lominadze J.G., Machabeli G. Z., & Usov V. V., 1983, Ap&SS, 90, 19L
- Machabeli, G. Z., & Usov, V. V., 1979 AZh Pi’sma, 5, 445
- Machabeli, G. Z., Luo, Q., Vladimirov, S. V., & Melrose, D. B., 2002, Phys. Rev. E, 65, 036408
- Malov, I, F., Machabeli, G. Z., 2002, Astronomy Reports, Vol. 46, Issue 8, p.684
- Malov, I. F., Machabeli, G. Z., 2007, Ap&SS, 29M
- Motch, C., Zavlin, V. E., & Haberl, F. 2003, A&A, 408, 323
- Osmanov, Z., Shapakidze, D., & Machbeli, G., 2009, A&A, 509, 19
- Pacholczyk, A. G. 1970, Radio Astrophysics (San Francisco: W. H. Freeman)
- Pavlov, G. G., Zavlin, V. E., Trümper, J., Neuhäuser, R., 1996, ApJ, 472, L33
- Pavlov, G. G., 2000, Talk at the ITP/UCSB workshop ”Spin and Magnetism of Young Neutron Stars”
- Pavlov, G, G., Zavlin, V. E., Sanwal, D., in Neutron Stars and Supernova Remnants. Eds. W. Becher, H. Lesch, & J. Trümper, 2002, MPE Report 278, 273
- Pons, J. A., Walter, F. M., Lattimer, J. M., Prakash, M., Neuhaäuser, R., An, P., 2002, ApJ, 564, 981
- Ransom, S. M., Gaensler, B. M., Slane, P.O. 2002, ApJ, 570, L75
- Rea, N. er al. 2007, MNRAS, 379, 1484
- Schwope, A. D., Erben, T., Kohnert, J., Lamer, G., Steinmetz, M., et al. 2009, A&A, 499, 267S
- Sturrock P. A., 1971, ApJ, 164, 529
- Tiengo, A., Mereghetti, S., 2007, ApJ, 657, L101
- Treves, A., Turolla, R., Zane, S., & Colpi, M. 2000, PASP, 112, 297
- Turolla, R., Zane, S. Drake, J. J., 2004, ApJ, 603, 265
- Vedenov A.A., Velikhov E.P. & Sagdeev R.Z., 1961, Soviet Physics Uspekhi, Volume 4, Issue 2, 332
- Walter, F. M., Wolk, S. J., Neuhäuser, R., 1996, Nature 379,233
- Zampieri, L., et al. 2001, A&A, 378, L5
- Zane, S., Cropper,M., Turolla, R., et al. 2005, ApJ, 627, 397
- Zane, S., et al. 2008, ApJ, 682, 487
- Zhelezniakov, V. V., ”Volni v Kosmicheskoi Plazme” Nauka, Moskva 1977