Relativistic models of magnetars: the twisted torus magnetic field configuration

We find general relativistic solutions of equilibrium magnetic field configurations in magnetars, extending previous results of Colaiuda et al. Our method is based on the solution of the relativistic Grad–Shafranov equation, to which Maxwell's equations can be reduced. We obtain equilibrium solutions with the toroidal magnetic field component confined into a finite region inside the star, and the poloidal component extending to the exterior. These so-called twisted torus configurations have been found to be the final outcome of dynamical simulations in the framework of Newtonian gravity, and appear to be more stable than other configurations. The solutions include higher-order multipoles, which are coupled to the dominant dipolar field. We use arguments of minimal energy to constrain the ratio of the toroidal to the poloidal field.     

Early evolution of newly born magnetars with a strong toroidal field

We present a state-of-the-art scenario for newly born magnetars as strong sources of gravitational waves (GWs) in the early days after formation. We address several aspects of the astrophysics of rapidly rotating, ultra-magnetized neutron stars (NSs), including early cooling before transition to superfluidity, the effects of the magnetic field on the equilibrium shape of NSs, the internal dynamical state of a fully degenerate, oblique rotator and the strength of the electromagnetic torque on the newly born NS. We show that our scenario is consistent with recent studies of supernova remnant surrounding Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) in the Galaxy that constrains the electromagnetic energy input from the central NS to be  ≤ 10⁵¹  erg. We further show that if this condition is met, then the GW signal from such sources is potentially detectable with the forthcoming generation of GW detectors up to Virgo cluster distances where an event rate ∼1 yr⁻¹ can be estimated. Finally, we point out that the decay of an internal magnetic field in the 10¹⁶ G range couples strongly with the NS cooling at very early stages, thus significantly slowing down both processes: the field can remain this strong for at least 10³ yr, during which the core temperature stays higher than several times 10⁸ K.     

Magnetic deformation of magnetars for the giant flares of the soft gamma-ray repeaters

We present one possible mechanism for the giant flares of the soft gamma-ray repeaters (SGRs) within the framework of the magnetar (superstrongly magnetized neutron star) model, motivated by the positive period increase associated with the August 27 event from SGR 1900+14. From second-order perturbation analysis of the equilibrium of the magnetic polytrope, we find that there exist different equilibrium states separated by the energy of the giant flares and the shift in the moment of inertia to cause the period increase. This suggests that, if we assume that global reconfiguration of the internal magnetic field of     suddenly occurs, the positive period increase     as well as the energy ≳10⁴⁴ erg of the giant flares may be explained. The moment of inertia can increase with a release of energy, because the star shape deformed by the magnetic field can be prolate rather than oblate. In this mechanism, since oscillation of the neutron star will be excited, a ∼ ms-period pulsation of the burst profile and an emission of gravitational waves are expected. The gravitational waves could be detected by planned interferometers such as LIGO, VIRGO and LCGT.     

Magnetic fields in axisymmetric neutron stars

We derive general equations for axisymmetric Newtonian magnetohydrodynamics and use these as the basis of a code for calculating equilibrium configurations of rotating magnetized neutron stars in a stationary state. We investigate the field configurations that result from our formalism, which include purely poloidal, purely toroidal and mixed fields. For the mixed-field formalism, the toroidal component appears to be bounded at less than 7 per cent. We calculate distortions induced both by magnetic fields and by rotation. From our non-linear work, we are able to look at the realm of validity of perturbative work: we find for our results that perturbative-regime formulae for magnetic distortions agree to within 10 per cent of the non-linear results if the ellipticity is less than 0.15 or the average field strength is less than 10¹⁷ G. We also consider how magnetized equilibrium structures vary for different polytropic indices.     

Topology of magnetars external field – I. Axially symmetric fields

There is an increasing theoretical and observational evidence that the external magnetic field of magnetars may contain a toroidal component, likely of the same order of the poloidal one. Such 'twisted magnetospheres' are threaded by currents flowing along the closed field lines which can efficiently interact with soft thermal photons via resonant cyclotron scatterings (RCS). Actually, RCS spectral models proved quite successful in explaining the persistent ∼1–10 keV emission from the magnetar candidates, the soft γ-ray repeaters (SGRs) and the anomalous X-ray pulsars (AXPs). Moreover, it has been proposed that, in the presence of highly relativistic electrons, the same process can give rise to the observed hard X-ray spectral tails extending up to  ∼200 keV  . Spectral calculations have been restricted up to now to the case of a globally twisted dipolar magnetosphere, although there are indications that the twist may be confined only to a portion of the magnetosphere, and/or that the large-scale field is more complex than a simple dipole. In this paper, we investigate multipolar, force–free magnetospheres of ultramagnetized neutron stars. We first discuss a general method to generate multipolar solutions of the Grad-Schlüter-Shafranov (GSS) equation, and analyse in detail dipolar, quadrupolar and octupolar fields. The spectra and lightcurves for these multipolar, globally twisted fields are then computed using a Monte Carlo code and compared with those of a purely dipolar configuration. Finally, the phase-resolved spectra and energy-dependent lightcurves obtained with a simple model of a locally sheared field are confronted with the International Gamma-Ray Astrophysics Laboratory ( INTEGRAL ) observations of the AXPs 1RXS J1708−4009 and 4U 0142+61. Results support a picture in which the field in these two sources is not globally twisted.     

Heating and cooling of magnetars with accreted envelopes

We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (   B ≳ 10¹⁴  G) non-accreted and accreted outermost envelopes composed of different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.     

The birthrate of magnetars

Magnetars, neutron stars with ultrastrong magnetic fields  ( B ∼ 10¹⁴−10¹⁵G)  , manifest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray pulsars. This study estimates the birthrate of magnetars to be ∼0.22 per century with a Galactic population comprising ∼17 objects. A population synthesis was carried out based on the five anomalous X-ray pulsars detected in the ROSAT All Sky Survey by comparing their number to that of massive OB stars in a well-defined volume. Additionally, the group of seven X-ray dim isolated neutron stars detected in the same survey were found to have a birthrate of ∼2 per century with a Galactic population of ∼22 000 objects.     

The XMM‐Newton view of magnetars

Anomalous X‐ray Pulsars and Soft Gamma‐ray Repeaters are believed to be magnetars: isolated neutron stars powered by the decay of extremely high magnetic fields. We review some of the main results obtained with XMM‐Newton and discuss the prospects for future observations of this small but extremely interesting class of objects. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Constraints on the magnetic field geometry of magnetars

We study the effect of the magnetic field geometry on the oscillation spectra of strongly magnetized stars. We construct a configuration of magnetic field where a toroidal component is added to the standard poloidal one. We consider a star with a type I superconductor core so that both components of the magnetic field are expelled from the core and confined in the crust. Our results show that the toroidal contribution does not influence significantly the torsional oscillations of the crust. On the contrary, the confinement of the magnetic field in the crust drastically affects the torsional oscillation spectrum. A comparison with estimations for the magnetic field strength, from observations, excludes the possibility that magnetars will have a magnetic field solely confined in the crust, that is, our results suggest that the magnetic field in whatever geometry has to permeate the whole star.     

On the quasi-periodic oscillations in magnetars

We study torsional Alfvén oscillations of magnetars, that is neutron stars with a strong magnetic field. We consider the poloidal and toroidal components of the magnetic field and a wide range of equilibrium stellar models. We use a new coordinate system  ( X , Y )  , where     and     and a ₁ is the radial component of the magnetic field. In this coordinate system, the one+two-dimensional evolution equation describing the quasi-periodic oscillations (QPOs), see Sotani et al., is reduced to a one+one-dimensional equation where the perturbations propagate only along the y -axis. We solve the one+one-dimensional equation for different boundary conditions and the open magnetic field lines, that is magnetic field lines that reach the surface and there match up with the exterior dipole magnetic field as well as closed magnetic lines, i.e. magnetic lines that never reach the stellar surface. For the open field lines, we find two families of QPO frequencies: a family of 'lower' QPO frequencies which is located near the x -axis and a family of 'upper' frequencies located near the y -axis. According to Levin, the fundamental frequencies of these two families can be interpreted as the turning point of the continuous spectrum. We find that the upper frequencies are multiples of the lower ones by a constant equalling  2 n + 1  . For the closed lines, the corresponding factor is   n + 1  . By using these relations, we can explain both the lower and the higher observed frequencies in SGR 1806−20 and SGR 1900+14.     

Progenitors with enhanced rotation and the origin of magnetars

Among the dozen known magnetar candidates, there are no binary objects. Given that the fraction of binary neutron stars is estimated to be about 3–10 per cent, it is reasonable to address the question of solitarity of magnetars, to estimate theoretically the fraction of binary objects among them, and to identify the most probable companions. We present population synthesis calculations of massive binary systems. In this study, we adopt the hypothesis that magnetic field of a magnetar is generated at the protoneutron star stage due to a dynamo mechanism, so rapid rotation of the core of a progenitor star is essential. Our goal is to estimate the number of neutron stars originated from progenitors with enhanced rotation. In our calculations, the fraction of neutron stars originating from such progenitors is about 8–9 per cent. This should be considered as an upper limit to the fraction of magnetars, as some of the progenitors can lose momentum. Most of these objects are isolated due to coalescences of components prior to neutron star formation, or due to system disruption after the second supernova explosion. The fraction of such neutron stars in surviving binaries is about 1 per cent or lower. Their most numerous companions are black holes.     

Multidimensional thermal structure of magnetized neutron star envelopes

Recently launched X-ray telescopes have discovered several candidate isolated neutron stars. The thermal radiation from these objects may potentially constrain our understanding of nuclear physics in a realm inaccessible to terrestrial experiments. To translate the observed fluxes from neutron stars into constraints, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We describe models of the thermal structure of the envelopes of neutron stars with magnetic fields up to 10¹⁴ G. Unlike earlier work, we infer the properties of envelope models in two dimensions and precisely account for the quantization of the electron phase-space. Both dipole and uniformly magnetized envelopes are considered.     

Almost analytic models of ultramagnetized neutron star envelopes

Recent ROSAT measurements show that the X-ray emission from isolated neutron stars is modulated at the stellar rotation period. To interpret these measurements, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We present nearly analytic models of the thermal structure of the envelopes of ultramagnetized neutron stars. Specifically, we examine the limit in which only the ground Landau level is filled. We use the models to estimate the amplitude of modulation expected from non-uniformities in the surface temperatures of strongly magnetized neutron stars. In addition, we estimate cooling rates for stars with fields B  ∼ 10¹⁵ − 10¹⁶ G, which are relevant to models that invoke 'magnetars' to account for soft γ-ray emission from some repeating sources.