Modelling of the surface emission of the low magnetic field magnetar SGR 0418+5729

We perform a detailed modelling of the post-outburst surface emission of the low magnetic field magnetar SGR 0418+5729. The dipolar magnetic field of this source, B=6×1012G estimated from its spin-down rate, is in the observed range of magnetic fields for normal pulsars. The source is further characterized by a high pulse fraction and a single-peak profile. Using synthetic temperature distribution profiles, and fully accounting for the general-relativistic effects of light deflection and gravitational redshift, we generate synthetic X-ray spectra and pulse profiles that we fit to the observations. We find that asymmetric and symmetric surface temperature distributions can reproduce equally well the observed pulse profiles and spectra of SGR 0418. None the less, the modelling allows us to place constraints on the system geometry (i.e. the angles ψ and ξ that the rotation axis makes with the line of sight and the dipolar axis, respectively), as well as on the spot size and temperature contrast on the neutron star surface. After performing an analysis iterating between the pulse profile and spectra, as done in similar previous works, we further employed, for the first time in this context, a Markov-Chain Monte Carlo approach to extract constraints on the model parameters from the pulse profiles and spectra, simultaneously. We find that, to reproduce the observed spectrum and flux modulation: (a) the angles must be restricted to 65° ≲ ψ + ξ ≲ 125° or 235° ≲ ψ + ξ ≲ 295°; (b) the temperature contrast between the poles and the equator must be at least a factor of ∼6, and (c) the size of the hottest region ranges between 0.2 and 0.7 km (including uncertainties on the source distance). Lastly, we interpret our findings within the context of internal and external heating models.

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.

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.

Chandra and Suzaku observations of the Be/X-ray star HD110432

We present an analysis of a pointed 141 ks Chandra high-resolution transmission gratings observation of the Be X-ray emitting star HD110432, a prominent member of the γ Cas analogs. This observation represents the first high-resolution spectrum taken for this source as well as the longest uninterrupted observation of any γ Cas analog. The Chandra light curve shows a high variability but its analysis fails to detect any coherent periodicity up to a frequency of 0.05 Hz. Hardness ratio versus intensity analyses demonstrate that the relative contributions of the [1.5-3] Å, [3-6] Å, and [6-16] Å energy bands to the total flux change rapidly in the short term. The analysis of the Chandra High Energy Transmission Grating (HETG) spectrum shows that, to correctly describe the spectrum, three model components are needed. Two of those components are optically thin thermal plasmas of different temperatures (kT ≈ 8-9 and 0.2-0.3 keV, respectively) described by the models vmekal or bvapec. The Fe abundance in each of these two components appears equal within the errors and is slightly subsolar with Z ≈ 0.75 Z ☉. The bvapec model better describes the Fe L transitions, although it cannot fit well the Na XI Lyα line at 10.02 Å, which appears to be overabundant. Two different models seem to describe well the third component. One possibility is a third hot optically thin thermal plasma at kT = 16-21 keV with an Fe abundance Z ≈ 0.3 Z ☉, definitely smaller than for the other two thermal components. Furthermore, the bvapec model describes well the Fe K shell transitions because it accounts for the turbulence broadening of the Fe XXV and Fe XXVI lines with a v turb ≈ 1200 km s–1. These two lines, contributed mainly by the hot thermal plasma, are significantly wider than the Fe Kα line whose FWHM < 5 mÅ is not resolved by Chandra. Alternatively, the third component can be described by a power law with a photon index of Γ = 1.56. In either case, the Chandra HETG spectrum establishes that each one of these components must be modified by distinct absorption columns. The analysis of a noncontemporaneous 25 ks Suzaku observation shows the presence of a hard tail extending up to at least 33 keV. The Suzaku spectrum is described with the sum of two components: an optically thin thermal plasma at kT ≈ 9 keV and Z ≈ 0.74 Z ☉, and a very hot second plasma with kT ≈ 33 keV or, alternatively, a power law with photon index of Γ = 1.58. In either case, each one of the two components must be affected by different absorption columns. Therefore, the kT = 8-9 keV component is definitely needed while the nature of the harder emission cannot be unambiguously established with the present data sets. The analysis of the Si XIII and S XV He-like triplets present in the Chandra spectrum points to a very dense (ne ~ 1013 cm–3) plasma located either close to the stellar surface (r < 3R *) of the Be star or, alternatively, very close (r ~ 1.5R WD) to the surface of a (hypothetical) white dwarf companion. We argue, however, that the available data support the first scenario.

A Be-type star with a black-hole companion

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars1, 2, 3 (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the approximately 80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions2, 3, 4. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (ref. 5; also known as HD 215227), although that conclusion was based on a single radial velocity curve of the Be star, a mistaken spectral classification6 and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of MWC 656. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe ii profiles from the equatorial disk, and a refined Be classification (to that of a B1.5–B2 III star), indicates that a black hole of 3.8 to 6.9 solar masses orbits MWC 656, the candidate counterpart of the γ-ray source AGL J2241+4454 (refs 5, 6). The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 × 10−7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect in conventional X-ray surveys.

Population synthesis of isolated neutron stars with magneto-rotational evolution

We revisit the population synthesis of isolated radio-pulsars incorporating recent advances on the evolution of the magnetic field and the angle between the magnetic and rotational axes from new simulations of the magneto-thermal evolution and magnetosphere models, respectively. An interesting novelty in our approach is that we do not assume the existence of a death line. We discuss regions in parameter space that are more consistent with the observational data. In particular, we find that any broad distribution of birth spin periods with P0 ≲ 0.5 s can fit the data, and that if the alignment angle is allowed to vary consistently with the torque model, realistic magnetospheric models are favoured compared to models with classical magneto-dipolar radiation losses. Assuming that the initial magnetic field is given by a lognormal distribution, our optimal model has mean strength 〈log B0[G]〉 ≈ 13.0–13.2 with width σ(log B0) = 0.6–0.7. However, there are strong correlations between parameters. This degeneracy in the parameter space can be broken by an independent estimate of the pulsar birth rate or by future studies correlating this information with the population in other observational bands (X-rays and γ-rays).

Ion balance in waters through inductively coupled plasma optical emission spectrometry

Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) has been employed to carry out the determination of both major anions and cations in water samples. The anion quantification has been performed by means of a new automatic accessory. In this device chloride has been determined by continuously adding a silver nitrate solution. As a result solid silver chloride particles are formed and retained on a nylon filter inserted in the line. The emission intensity is read at a silver characteristic wavelength. By plotting the drop in silver signal versus the chloride concentration, a straight line is obtained. As regards bicarbonate, this anion has been on-line transformed into carbon dioxide with the help of a 2.0 mol L−1 nitric acid stream. Carbon signal is linearly related with bicarbonate concentration. Finally, information about sulfate concentration has been achieved by means of the measurement of sulfur emission intensity. All the steps have been simultaneously and automatically performed. With this setup detection limits have been 1.0, 0.4 and 0.09 mg L−1 for chloride, bicarbonate and sulfate, respectively. Furthermore, it affords good precision with RSD below 6 %. Cation (Ca, Mg, Na and K) concentration, in turn, has been obtained by simultaneously reading the emission intensity at characteristic wavelengths. The obtained limits of detection have been 8 × 10−3, 2 × 10−3, 8 × 10−4 and 10−2 mg L−1 for sodium, potassium, magnesium and calcium, respectively. As regards sample throughput, about 30 samples h−1 can be analysed. Validation results have revealed that the obtained concentrations for these anions are not significantly different as compared to the data provided by conventional methods. Finally, by considering the data for anions and cations, precise ion balances have been obtained for well and mineral water samples.