arxiv: v1 [astro-ph.he] 17 Nov 2011

Transcription

1 arxiv: v1 [astro-ph.he] 17 Nov 2011 The influence of magnetic field geometry on magnetars X-ray spectra D Viganò 1, N Parkins 2,3, S Zane 3, R Turolla 4,3, J A Pons 1 and J A Miralles 1 1 Departament de Fisica Aplicada, Universitat d Alacant, Ctra. San Vicent del Raspeig, s/n San Vicent del Raspeig, Spain 2 University of Liverpool, Liverpool L69 3BX, UK 3 Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey RH5 6NT, UK 4 Università di Padova, Dipartimento di Fisica, via F. Marzolo 8, I Padova, Italy daniele.vigano@ua.es Abstract. Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars or magnetars is one of the most valuable tools to gain insight into the physical processes occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering of the thermal radiation coming from the star surface. Motivated by the increased quality of the observational data, much theoretical work has been devoted to develop Monte Carlo (MC) codes that incorporate the effects of resonant Compton scattering (RCS) in the modeling of radiative transfer of photons through the magnetosphere. The two key ingredients in this simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The MF geometry is expected to be complex, but up to now only mathematically simple solutions (self-similar solutions) have been employed. In this work, we discuss the effects of new, more realistic, MF geometries on synthetic spectra. We use new force-free solutions [14] in a previously developed MC code [9] to assess the influence of MF geometry on the emerging spectra. Our main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters of the magnetospheric plasma, but the MF geometry plays an important role on the angledependence of the spectra. 1. Introduction Neutron stars (NSs) can be powered by different sources of energy: rotation, accretion, residual heat and magnetic field. Up to now, more than twenty X-ray pulsating sources are identified as isolated, slowly rotating (P 2 10 s) neutron stars, whose large X-ray luminosity L x erg/s cannot usually be explained in terms of rotational energy, at variance with rotation-powered pulsars (RPPs). Historically, they are classified as Anomalous X-ray Pulsars (AXPs) or Soft Gamma Repeaters (SGRs), and they are recognized as magnetars. In these objects, MF decay and evolution is the dominant energy source and crustal deformations produced by the huge internal magnetic stresses are the likely cause of bursts and episodic giant flares. In most cases, the spin-down inferred external magnetic field turns out to be super-strong. Recently, increasing evidence gathered in favor of a unified scenario for isolated neutron stars,

2 their different observed manifestations being mainly due to a different topology and strength of the internal MF, and to evolution ([11]). The typical soft X-ray (< 10 kev) spectrum of magnetars is given by a thermal component (kt kev), plus a hard tail phenomenologically described by a power law with a photon index Γ 3 4 [7]. This hard tail is believed to be produced by resonant compton (up)scattering (RCS) of photons emitted by the star surface onto mildly relativistic particles flowing in magnetosphere. In order to take into account for this effect, several codes have been developedin recent years [4, 8, 9, 5]. Thetwobasicinputsaretheexternal MFandaprescription for the particle space and velocity distribution. With these ingredients, and with a given energy distribution for the incoming surface photons (the seed spectrum), the MC simulations of RCS provide the emerging, reprocessed spectrum, which must be compared with observational data. Magnetars also exhibit a high-energy ( kev) power-law tail, the origin of which is still under debate. One of the possible mechanisms is again RCS, whose efficiency at high energy is very sensitive to MF and plasma velocity. In this paper we discuss the effect of the geometry of the MF, which is obtained from new calculations of magnetostatic, force-free equilibria. 2. Force-free twisted magnetosphere The shape of the internal MF of NSs is largely unknown, but it is thought to be complex, with toroidal and poloidal components of similar strength, and with significant deviations from a dipolar geometry. In the crust, a large toroidal field is a necessary ingredient to explain the observed thermal and bursting properties of magnetars, even in the case in which they show a low value of the inferred external dipolar field [13]. The complex internal topology must be smoothly matched to an external solution, so that the internal MF geometry at the star surface has a strong effect on the external MF. In this framework, the RPP phenomenolgy is expected to be better explained by the presence of a NS with a simple, nearly potential, dipolar field, while magnetars activity is compatible with the presence of a complex external field, as a result of the transfer of magnetic helicity from the internal to the external field. A non-potential MF needs supporting currents, and charges can be either extracted from the surface or produced by one-photon pair creation in the strong MF [1]. In the inner magnetospheric region MF lines corotate with the star, with footprints connected with the crustal field, which evolves on long timescales ( years). Furthermore, resistive processes in the magnetosphere act on a typical timescale of years [2], much longer than the typical response of the tenuous plasma, whose Alfvén velocity is near to the velocity of light. For these reasons, a reasonable approximation is to build stationary, force-free magnetospheric solutions, assuming that the footprints are anchored in the crust and neglecting resistivity. Then, the ideal MHD force-free equation is a simple balance between electromagnetic forces, provided that inertial, pressure, and gravity terms are negligible ρ e E + 1 c J B = 0, (1) where we have introduced the rotationally induced electric field E = v rot B, the related charge density ρ e and the current density J = c B/4π. The electric field gives corrections of order (ρω/c) 2. Neglecting such contributions is safe if we consider the region closest to a magnetar surface (r 100r ), and remembering that the light cylinder for slow rotators lies typically at a few thousands stellar radii (r l = c/ω and magnetar periods are in the 1 10 s range). Therefore, condition (1) reduces simply to J B: the particles can flow only along MF lines. Even if the pulsar equation simplifies in this limit, it still contains an arbitrary function, which basically describes the toroidal field as a function of the poloidal one (see [14] for the mathematical formulation). In the literature, only the simplest (semi-)analytical solutions have been derived and later employed in applications. In

3 Figure 1. Numerical model 1 (left) and model 2 (right). First panel: B (grey logarithmic scale) with scattering surfaces for photons of 1, 3 and 5 kev (here β = 0.5, µ = 0.5; see text); second panel: current intensity J /c (grey linear scale), in units of G/r. particular, in the context of magnetars, the self-similar twisted dipole family of solutions is the most popular choice. Proposed for solar corona studies [6], and later extended to the magnetar scenario [12], in self-similar models the radial dependence of the MF components is described by a simple power law, B i r (p+2). These solutions are semi-analytical and described by only one parameter, the radial index p (see [10] for globally twisted multipolar fields). However, globally twisted solutions necessary evolve in time and magnetic configurations in which the twist affects only a limited bundle of field lines are likely required by observations [2]. As an alternative to self-similar models, in [14] we presented a numerical method to build general solutions. The force-free solutions are obtained by iteration starting from an arbitrary, non-force-free poloidal plus toroidal field. By the introduction of a fictitious velocity field proportional to the Lorentz force, the code is designed to dissipate only the part of the current not parallel to the MF. The boundaryconditions are an arbitrary, fixed radial MF at the surface, B r (r,θ), andtherequirementofacontinuousmatchingwithapotential dipoleat largedistances r out 100r. This method is able to produce general configurations, whose features may be radically different from the self-similar solutions. In this work, we try to reproduce some j- bundle like structure, aiming at providing input for MC simulations [9]. In particular, we focus on two models, whose MF and current distribution are shown in Fig. 1. The firstone is obtained fixing a dipolar form for radial MF, while the second one has a concentration of field lines near the south pole. To get these profiles, we have chosen the φ component of the vector potential as: model 1: model 2: sinθ A φ (θ,r ) = B 0 2 sinθ A φ (θ,r ) = B 0 2 [ ] 1+4e (θ 0.9π) (2) (3) from which we recover B r = ( A φˆφ)r (B 0 is the value of B r at north pole). Both models have high helicity, but it is not homogeneously distributed like in self-similar models. Model 1 is characterized by an equatorial asymmetry, with current and twist concentrated in a closed bundle near the equatorial region and a j bundle near the southern semi-axis. Model 2 is more extreme: its helicity is 25 times larger than for model 1. Both are supposed to be more realistic than the self-similar models, since currents are more concentrated near the axis.

4 3. Resonant Compton scattering The MC code we employ is presented in [8, 9] and we refer to these works for further details. Here we just summarize some basic concepts. It is assumed that currents flow in a non-potential magnetosphere and that the electron density is large enough for the medium to be optically thick to resonant scattering: photons produced by the cooling star surface can be then up-scattered by RCS, populating the hard tail of the spectrum. In the stellar frame, the resonant energy corresponds to the Doppler-shifted gyration frequency of the electron ω D (B,β) = 1 eb γ(1 βµ) m e c (4) where B = B is the local intensity of MF, cβ the electron velocity, γ = 1/ 1 β 2 the Lorentz factor, and µ the cosine of the angle between the photon direction and the electron momentum. The non-resonant contributions are negligible, since the resonant cross-section is up to five orders of magnitude larger than the Thomson cross-section. Figs. 1 show the scattering surfaces in the two models, for photons of 1, 3 and 5 kev, assuming β = 0.5 and µ = 0.5. The surfaces are far from being spherically symmetric. In model 2, which has the strongest helicity, the scattering surfaces of soft X-ray photons lie tens of stellar radii away from the surface, especially near the axis. The scattering optical depth for a given seed photon energy depends on the local particle density n e, the velocity of the scatterer and the local value of the MF. Since photons can be pushed to energies of some hundreds kevs, non-conservative scattering and the fully relativistic QED cross sections are used [9]. The velocity distribution is assumed to be a one-dimensional (along the field line), relativistic Maxwellian distribution with a plasma temperature T e and a bulk velocity β bulk, as described in [8]. The strong simplification in this approach is that the velocity distribution does not depend on position, which reduces the microphysical inputs to two parameters: the temperature of the plasma T e and β bulk. Only recently, a more consistent treatment of the charges space and velocity distribution was attempted [3]. These new results show that the velocity distribution is non-trivial and strongly depends on position, which opens a new range of possibilities, so far unexplored. Given a seed spectrum for the photons emitted from the surface (assumed as a blackbody with temperature T ), the MC code follows the photon propagation through the scattering magnetosphere. When a photon enters in a parameter region where no more resonant scatterings are possible, the photon escapes to infinity, where its energy and direction are stored. The sky at infinity is divided into patches, so that viewing effects can be accounted for (e.g. if line-of-sight is along the angles θ s, φ s, only photons collected in the patch which contains that pair of angles are analyzed). The angle-averaged spectrum is obtained simply by averaging over all patches. 4. Results and discussion In previous works [15], self-similar magnetospheric models were used to obtain synthetic spectra that were satisfactorily fitted to X-ray observations, giving typical values of β bulk [0.1,0.7] and φ [0.4,2.0]. In this paper we explore the dependence of the spectrum on the magnetosphere model by comparing results obtained with self-similar models and with our new numerical solutions for force-free magnetospheres. We present results comparing models with a magnetic field intensity at the pole of B p = G. The seed spectrum is assumed to be a kt = 0.4 kev Planckian isotropic distribution for both ordinary and extraordinary photons. In Figs. 2 and 3 we plot the angle-averaged spectra, i.e., the distribution of all photons escaped to infinity, independently on the final direction. Fig. 2 shows the effects produced by changes in the bulk velocity, keeping fixed the MF configuration (model 1), and the electron temperature (kt e = 20 kev). The main fact is the the spectrum in the region E > 10 kev becomes harder as

5 Figure 2. Angle-averaged spectra for kt e = 20 kev, model 1, changing β bulk. Figure 3. Angle-averaged spectra for kt e = 20 kev, β bulk = 0.9, changing MF model. β bulk is increased, while it is depleted of photons of energies in the 1 10 kev range. A similar, but less pronounced effect, is obtained by increasing the electron temperature [9]. In Fig. 3 we compare spectra obtained with fixed values of kt e = 20 kev and β bulk = 0.9, varying only the MF topology. The three lines correspond to model 1 (solid), model 2 (dashes) and a self-similar model (dash-dotted line) with φ = 1.36, which approximately has the same total helicity of model1. Wehave chosen ahighvalueof β bulk toshow acasewheretheeffects are larger, but our conclusions do not change qualitatively for lower values of β bulk. The high-energy cutoff produced by electron recoil is also evident. We see that the major differences arising from changes in MF topology appear in the hard tail of the spectra, which is clearly harder for the self-similar model. However, the thermal part of the spectrum is more depleted, especially for model 2. However, in the angle-averaged spectrum, the differences due to the MF topology are partially smeared out, and the comparison between synthetic spectra as seen from different viewing angles (first three panles in Fig. 4) is more interesting. The comptonization degree reflects the inhomogeneous particle density distribution and therefore the particular current distribution produces important differences between different viewing angles. Since neutron stars rotate, the study of pulse profiles and phase-resolved spectra can trace the geometric features of the scattering region (tens of stellar radii). In the lower right panel we show the light curves for the three MF models, in bands kev and kev, for an oblique rotator with a certain line of sight (ξ = χ = 30 in [8]). For self-similar models (lower left panel), the thermal part of the spectrum shows a smooth dependence with the viewing angle, which in turn translates into relatively regular light curves. In models 1 and 2, differences in the spectra induced by viewing angles are more pronounced. In these cases, the spectrum is much more irregular and asymmetries are larger. In particular, model 1 (upper left) shows a softer spectrum when seen from northern latitudes, with important spectral differences. The pulsed fraction of model 1 is very high in the hard X-ray band, while it is comparable with the self-similar model in the soft range. Model 2 (upper right) has a more symmetric distribution of currents, and the comptonization degree at different angles depends on the energy band in a non trivial way. This results in large pulsed fraction for both energy ranges and in notable differences between their pulse profiles. The high energy tail, seen at different colatitudes, can vary by one order of magnitude or more. Investigating the geometry can help to recognize the different components seen in the hard tails via pulse phase spectroscopy. On the other hand, the macrophysical approach used in this work to obtain the MF topology should be accompanied by the corresponding microphysical

6 Figure 4. Synthetic spectra for model 1 (upper left panel), model 2 (upper right), self-similar model(lower left) with the parameters indicated in the figures. Legends indicate the four viewing angles here considered: aligned with north or south pole, just above and just below the equator. The seed blackbody is shown for comparison (solid line). Fourth panel shows the pulse profiles obtained for model 1 (solid line), model 2 (dotted), self-similar (dashed), in the energy ranges kev (thin lines) and kev (thick lines), with ξ = χ = 30 (see text). description to determine the velocity distribution of particles, here simply taken as a model parameter. A fully consistent description of both MF geometry and the spatial and velocity distribution of particles is needed to advance our interpretation of magnetar spectra. References [1] Beloborodov A M and Thompson C 2007 ApJ [2] Beloborodov A M 2009 ApJ [3] Beloborodov A M 2011 High-Energy Emission from Pulsars and their Systems ed Torres D F and Rea N p 299 [4] Fernández R and Thompson C 2007 ApJ [5] Fernández R and Davis S W 2011 ApJ [6] Low B C and Lou Y Q 1990 ApJ [7] Mereghetti S 2008 A&ARv [8] Nobili L, Turolla R and Zane S 2008 MNRAS [9] Nobili L, Turolla R and Zane S 2008 MNRAS [10] Pavan L, Turolla R, Zane S and Nobili L 2009 MNRAS [11] Pons J A and Perna R 2011 ApJ [12] Thompson C, Lyutikov M and Kulkarni S R 2002 ApJ [13] Turolla R, Zane S, Pons J A, Esposito P and Rea N 2011 ApJ [14] Viganò D, Pons J and Miralles J A 2011 A&A 533 A125 [15] Zane S, Rea N, Turolla R and Nobili L 2009 MNRAS