A new algorithm to accelerate harmonic analysis of time series

The Lomb periodogram has been traditionally a tool that allows us to elucidate if a frequency turns out to be important for explaining the behaviour of a given time series. Many linear and nonlinear reiterative harmonic processes that are used for studying the spectral content of a time series take into account this periodogram in order to avoid including spurious frequencies in their models due to the leakage problem of energy from one frequency to others. However, the estimation of the periodogram requires long computation time that makes the harmonic analysis slower when we deal with certain time series. Here we propose an algorithm that accelerates the extraction of the most remarkable frequencies from the periodogram, avoiding its whole estimation of the harmonic process at each iteration. This algorithm allows the user to perform a specific analysis of a given scalar time series. As a result, we obtain a functional model made of (1) a trend component, (2) a linear combination of Fourier terms, and (3) the so-called mixed secular terms by reducing the computation time of the estimation of the periodogram.

Force-free twisted magnetospheres of neutron stars

Context. The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. Methods. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. Results. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. Conclusions. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

Pulsar timing irregularities and the imprint of magnetic field evolution

Context. The rotational evolution of isolated neutron stars is dominated by the magnetic field anchored to the solid crust of the star. Assuming that the core field evolves on much longer timescales, the crustal field evolves mainly though Ohmic dissipation and the Hall drift, and it may be subject to relatively rapid changes with remarkable effects on the observed timing properties. Aims. We investigate whether changes of the magnetic field structure and strength during the star evolution may have observable consequences in the braking index n. This is the most sensitive quantity to reflect small variations of the timing properties that are caused by magnetic field rearrangements. Methods. We performed axisymmetric, long-term simulations of the magneto-thermal evolution of neutron stars with state-of-the-art microphysical inputs to calculate the evolution of the braking index. Relatively rapid magnetic field modifications can be expected only in the crust of neutron stars, where we focus our study. Results. We find that the effect of the magnetic field evolution on the braking index can be divided into three qualitatively different stages depending on the age and the internal temperature: a first stage that may be different for standard pulsars (with n ~ 3) or low field neutron stars that accreted fallback matter during the supernova explosion (systematically n < 3); in a second stage, the evolution is governed by almost pure Ohmic field decay, and a braking index n > 3 is expected; in the third stage, at late times, when the interior temperature has dropped to very low values, Hall oscillatory modes in the neutron star crust result in braking indices of a high absolute value and both positive and negative signs. Conclusions. Current magneto-thermal evolution models predict a large contribution to the timing noise and, in particular, to the braking index, from temporal variations of the magnetic field. Models with strong (≳ 1014 G) multipolar or toroidal components, even with a weak (~1012 G) dipolar field are consistent with the observed trend of the timing properties.

A new code for the Hall-driven magnetic evolution of neutron stars

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

Unifying the observational diversity of isolated neutron stars via magneto-thermal evolution models

Observations of magnetars and some of the high magnetic field pulsars have shown that their thermal luminosity is systematically higher than that of classical radio-pulsars, thus confirming the idea that magnetic fields are involved in their X-ray emission. Here we present the results of 2D simulations of the fully coupled evolution of temperature and magnetic field in neutron stars, including the state-of-the-art kinetic coefficients and, for the first time, the important effect of the Hall term. After gathering and thoroughly re-analysing in a consistent way all the best available data on isolated, thermally emitting neutron stars, we compare our theoretical models to a data sample of 40 sources. We find that our evolutionary models can explain the phenomenological diversity of magnetars, high-B radio-pulsars, and isolated nearby neutron stars by only varying their initial magnetic field, mass and envelope composition. Nearly all sources appear to follow the expectations of the standard theoretical models. Finally, we discuss the expected outburst rates and the evolutionary links between different classes. Our results constitute a major step towards the grand unification of the isolated neutron star zoo.

The outburst decay of the low magnetic field magnetar SGR 0418+5729

We report on the long-term X-ray monitoring of the outburst decay of the low magnetic field magnetar SGR 0418+5729 using all the available X-ray data obtained with RXTE, Swift, Chandra, and XMM-Newton observations from the discovery of the source in 2009 June up to 2012 August. The timing analysis allowed us to obtain the first measurement of the period derivative of SGR 0418+5729: ˙ P = 4(1) × 10−15 s s−1, significant at a ∼3.5σ confidence level. This leads to a surface dipolar magnetic field of Bdip 6 × 1012 G. This measurement confirms SGR 0418+5729 as the lowest magnetic field magnetar. Following the flux and spectral evolution from the beginning of the outburst up to ∼1200 days, we observe a gradual cooling of the tiny hot spot responsible for the X-ray emission, from a temperature of ∼0.9 to 0.3 keV. Simultaneously, the X-ray flux decreased by about three orders of magnitude: from about 1.4 × 10−11 to 1.2 × 10−14 erg s−1 cm−2. Deep radio, millimeter, optical, and gamma-ray observations did not detect the source counterpart, implying stringent limits on its multi-band emission, as well as constraints on the presence of a fossil disk. By modeling the magneto-thermal secular evolution of SGR 0418+5729, we infer a realistic age of ∼550 kyr, and a dipolar magnetic field at birth of ∼1014 G. The outburst characteristics suggest the presence of a thin twisted bundle with a small heated spot at its base. The bundle untwisted in the first few months following the outburst, while the hot spot decreases in temperature and size. We estimate the outburst rate of low magnetic field magnetars to be about one per year per galaxy, and we briefly discuss the consequences of such a result in several other astrophysical contexts.

A strongly magnetized pulsar within the grasp of the milky way’s supermassive black hole

The center of our Galaxy hosts a supermassive black hole, Sagittarius (Sgr) A∗. Young, massive stars within 0.5 pc of Sgr A∗ are evidence of an episode of intense star formation near the black hole a few million years ago, which might have left behind a young neutron star traveling deep into Sgr A∗’s gravitational potential. On 2013 April 25, a short X-ray burst was observed from the direction of the Galactic center. With a series of observations with the Chandra and the Swift satellites, we pinpoint the associated magnetar at an angular distance of 2.4±0.3 arcsec from Sgr A∗, and refine the source spin period and its derivative (P = 3.7635537(2) s and ˙ P = 6.61(4) × 10−12 s s−1), confirmed by quasi simultaneous radio observations performed with the Green Bank Telescope and Parkes Radio Telescope, which also constrain a dispersion measure of DM = 1750 ± 50 pc cm−3, the highest ever observed for a radio pulsar. We have found that this X-ray source is a young magnetar at ≈0.07–2 pc from Sgr A∗. Simulations of its possible motion around Sgr A∗ show that it is likely (∼90% probability) in a bound orbit around the black hole. The radiation front produced by the past activity from the magnetar passing through the molecular clouds surrounding the Galactic center region might be responsible for a large fraction of the light echoes observed in the Fe fluorescence features.

The imprint of the crustal magnetic field on the thermal spectra and pulse profiles of isolated neutron stars

Isolated neutron stars (NSs) show a bewildering variety of astrophysical manifestations, presumably shaped by the magnetic field strength and topology at birth. Here, using state-of-the-art calculations of the coupled magnetic and thermal evolution of NSs, we compute the thermal spectra and pulse profiles expected for a variety of initial magnetic field configurations. In particular, we contrast models with purely poloidal magnetic fields to models dominated by a strong internal toroidal component. We find that, while the former displays double-peaked profiles and very low pulsed fractions, in the latter, the anisotropy in the surface temperature produced by the toroidal field often results in a single pulse profile, with pulsed fractions that can exceed the 50–60 per cent level even for perfectly isotropic local emission. We further use our theoretical results to generate simulated ‘observed’ spectra, and show that blackbody (BB) fits result in inferred radii that can be significantly smaller than the actual NS radius, even as low as ∼1–2 km for old NSs with strong internal toroidal fields and a high absorption column density along their line of sight. We compute the size of the inferred BB radius for a few representative magnetic field configurations, NS ages and magnitudes of the column density. Our theoretical results are of direct relevance to the interpretation of X-ray observations of isolated NSs, as well as to the constraints on the equation of state of dense matter through radius measurements.

Central compact objects and the hidden magnetic field scenario

Central compact objects (CCOs) are X-ray sources lying close to the centre of supernova remnants, with inferred values of the surface magnetic fields significantly lower (≲1011 G) than those of standard pulsars. In this paper, we revise the hidden magnetic field scenario, presenting the first 2D simulations of the submergence and re-emergence of the magnetic field in the crust of a neutron star. A post-supernova accretion stage of about 10−4–10−3 M⊙ over a vast region of the surface is required to bury the magnetic field into the inner crust. When accretion stops, the field re-emerges on a typical time-scale of 1–100 kyr, depending on the submergence conditions. After this stage, the surface magnetic field is restored close to its birth values. A possible observable consequence of the hidden magnetic field is the anisotropy of the surface temperature distribution, in agreement with observations of several of these sources. We conclude that the hidden magnetic field model is viable as an alternative to the antimagnetar scenario, and it could provide the missing link between CCOs and the other classes of isolated neutron stars.

A highly resistive layer within the crust of X-ray pulsars limits their spin periods

The lack of isolated X-ray pulsars with spin periods longer than 12 s raises the question of where the population of evolved high-magnetic-field neutron stars has gone. Unlike canonical radiopulsars, X-ray pulsars are not subject to physical limits to the emission mechanism nor observational biases against the detection of sources with longer periods. Here we show that a highly resistive layer in the innermost part of the crust of neutron stars naturally limits the spin period to a maximum value of about 10–20 s. This highly resistive layer is expected if the inner crust is amorphous and heterogeneous in nuclear charge, possibly owing to the existence of a nuclear ‘pasta’ phase. Our findings suggest that the maximum period of isolated X-ray pulsars may be the first observational evidence for an amorphous inner crust, whose properties can be further constrained by future X-ray timing missions combined with more detailed models.

Is SGR 0418+5729 indeed a waning magnetar?

SGR 0418+5729 is a transient soft gamma-ray repeater which underwent a major outburst in 2009 June, during which the emission of short bursts was observed. Its properties appeared quite typical of other sources of the same class until long-term X-ray monitoring failed to detect any period derivative. The present upper limit on P implies that the surface dipole field is Bp lsim 7.5 × 1012 G, well below those measured in other soft gamma-ray repeaters (SGRs) and in the Anomalous X-ray Pulsars (AXPs), a group of similar sources. Both SGRs and AXPs are currently believed to be powered by ultra-magnetized neutron stars (magnetars, Bp ≈ 1014-1015 G). SGR 0418+5729 hardly seems to fit in such a picture. We show that the magneto-rotational properties of SGR 0418+5729 can be reproduced if this is an aged magnetar, ≈1 Myr old, which experienced substantial field decay. The large initial toroidal component of the internal field required to match the observed properties of SGR 0418+5729 ensures that crustal fractures, and hence bursting activity, can still occur at the present time. The thermal spectrum observed during the outburst decay is compatible with the predictions of a resonant Compton scattering model (as in other SGRs/AXPs) if the field is low and the magnetospheric twist is moderate.

3XMM J185246.6+003317: another low magnetic field magnetar

We study the outburst of the newly discovered X-ray transient 3XMMJ185246.6+003317, re-analyzing all available XMM-Newton observations of the source to perform a phase-coherent timing analysis, and derive updated values of the period and period derivative. We find the source rotating at P = 11.55871346(6) s (90% confidence level; at epoch MJD 54728.7) but no evidence for a period derivative in the seven months of outburst decay spanned by the observations. This translates to a 3σ upper limit for the period derivative of ˙ P <1.4×10−13 s s−1, which, assuming the classical magneto-dipolar braking model, gives a limit on the dipolar magnetic field of Bdip < 4.1×1013 G. The X-ray outburst and spectral characteristics of 3XMM J185246.6+003317 confirm its identification as a magnetar, but the magnetic field upper limit we derive defines it as the third “low-B” magnetar discovered in the past 3 yr, after SGR 0418+5729 and Swift J1822.3−1606. We have also obtained an upper limit to the quiescent luminosity (<4×1033 erg s−1), in line with the expectations for an old magnetar. The discovery of this new low field magnetar reaffirms the prediction of about one outburst per year from the hidden population of aged magnetars.

A note on the periodic orbits of a self excited rigid body

The aim of the present paper is to study the periodic orbits of a perturbed self excited rigid body with a fixed point. For studying these periodic orbits we shall use averaging theory of first order.

Quiescent state and outburst evolution of SGR 0501+4516

We report on the quiescent state of the soft gamma repeater SGR 0501+4516 observed by XMM–Newton on 2009 August 30. The source exhibits an absorbed flux ∼75 times lower than that measured at the peak of the 2008 outburst, and a rather soft spectrum, with the same value of the blackbody temperature observed with ROSAT back in 1992. This new observation is put into the context of all existing X-ray data since its discovery in 2008 August, allowing us to complete the study of the timing and spectral evolution of the source from outburst until its quiescent state. The set of deep XMM–Newton observations performed during the few years time-scale of its outburst allows us to monitor the spectral characteristics of this magnetar as a function of its rotational period, and their evolution along these years. After the first ∼10 d, the initially hot and bright surface spot progressively cooled down during the decay. We discuss the behaviour of this magnetar in the context of its simulated secular evolution, inferring a plausible dipolar field at birth of 3 × 1014 G, and a current (magnetothermal) age of ∼10 kyr.

VdBH 222: a starburst cluster in the inner Milky Way

Context. It has been suggested that the compact open cluster VdBH 222 is a young massive distant object. Aims. We set out to characterise VdBH 222 using a comprehensive set of multi-wavelength observations. Methods. We obtained multi-band optical (UBVR) and near-infrared (JHKS) photometry of the cluster field, as well as multi-object and long-slit optical spectroscopy for a large sample of stars in the field. We applied classical photometric analysis, as well as more sophisticated methods using the CHORIZOS code, to determine the reddening to the cluster. We then plotted dereddened HR diagrams and determined cluster parameters via isochrone fitting. Results. We have identified a large population of luminous supergiants confirmed as cluster members via radial velocity measurements. We find nine red supergiants (plus one other candidate) and two yellow supergiants. We also identify a large population of OB stars. Ten of them are bright enough to be blue supergiants. The cluster lies behind ≈7.5 mag of extinction for the preferred value of RV = 2.9. Isochrone fitting allows for a narrow range of ages between 12 and 16 Ma. The cluster radial velocity is compatible with distances of ~6 and ~10 kpc. The shorter distance is inconsistent with the age range and Galactic structure. The longer distance implies an age ≈ 12 Ma and a location not far from the position where some Galactic models place the far end of the Galactic bar. Conclusions. VdBH 222 is a young massive cluster with a likely mass >20 000 M⊙. Its population of massive evolved stars is comparable to that of large associations, such as Per OB1. Its location in the inner Galaxy, presumably close to the end of the Galactic bar, adds to the increasing evidence for vigorous star formation in the inner regions of the Milky Way.