A search for thermally emitting isolated neutron stars in the 2XMMp catalogue

The relatively large number of nearby radio-quiet and thermally emitting isolated neutron stars (INSs) discovered in the ROSAT All-Sky Survey, dubbed the ""Magnificent Seven"", suggests that they belong to a formerly neglected major component of the overall INS population. So far, attempts to discover similar INSs beyond the solar vicinity failed to confirm any reliable candidate. The good positional accuracy and soft X-ray sensitivity of the EPIC cameras onboard the XMM-Newton satellite allow us to efficiently search for new thermally emitting INSs. We used the 2XMMp catalogue to select sources with no catalogued candidate counterparts and with X-ray spectra similar to those of the Magnificent Seven, but seen at greater distances and thus undergoing higher interstellar absorptions. Identifications in more than 170 astronomical catalogues and visual screening allowed us to select fewer than 30 good INS candidates. In order to rule out alternative identifications, we obtained deep ESO-VLT and SOAR optical imaging for the X-ray brightest candidates. We report here on the optical follow-up results of our search and discuss the possible nature of 8 of our candidates. A high X-ray-to-optical flux ratio together with a stable flux and soft X-ray spectrum make the brightest source of our sample, 2XMM J104608.7-594306, a newly discovered thermally emitting INS. The X-ray source 2XMM J010642.3+005032 has no evident optical counterpart and should be further investigated. The remaining X-ray sources are most probably identified with cataclysmic variables and active galactic nuclei, as inferred from the colours and flux ratios of their likely optical counterparts. Beyond the finding of new thermally emitting INSs, our study aims at constraining the space density of this Galactic population at great distances and at determining whether their apparently high density is a local anomaly or not.

Proper motions of thermally emitting isolated neutron stars measured with Chandra

The remarkable astrometric capabilities of the Chandra Observatory offer the possibility to measure proper motions of X-ray sources with an unprecedented accuracy in this wavelength range. We recently completed a proper motion survey of three of the seven thermally emitting radio-quiet isolated neutron stars (INSs) discovered in the ROSAT all-sky survey. These INSs (RXJ0420.0-5022, RXJ0806.4-4123 and RXJ1308.6+2127) either lack an optical counterpart or have one so faint that ground based or space born optical observations push the current possibilities of the instrumentation to the limit. Pairs of ACIS observations were acquired 3 to 5 years apart to measure the displacement of the sources on the X-ray sky using as a reference the background of extragalactic or remote Galactic X-ray sources. We derive 2 sigma upper limits of 123 mas yr(-1) and 86 mas yr(-1) on the proper motion of RXJ0420.0-5022 and RXJ0806.4-4123, respectively. RXJ1308.6+2127 exhibits a very significant displacement (similar to 9 sigma) yielding mu = 220 +/- 25 mas yr(-1), the second fastest measured among all ROSAT-discovered INSs. The source is probably moving away rapidly from the Galactic plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of similar to 740 (d/700 pc) km s(-1) might be the largest of all ROSAT INSs and its corresponding spatial velocity lies among the fastest recorded for neutron stars. RXJ1308.6+2127 is thus a middle-aged (age similar to 1 My) high velocity cooling neutron star. We investigate its possible origin in nearby OB associations or from a field OB star. In most cases, the flight time from birth place appears significantly shorter than the characteristic age derived from spin down rate. Overall, the distribution in transverse velocity of the ROSAT INSs is not statistically different from that of normal radio pulsars.

Maximum mass of neutron stars

We determine the structure of neutron stars within a Brueckner-Hartree-Fock approach based on realistic nucleon-nucleon, nucleon-hyperon, and hyperon-hyperon interactions. Our results indicate rather low maximum masses below 1.4 solar masses. This feature is insensitive to the nucleonic part of the EOS due to a strong compensation mechanism caused by the appearance of hyperons and represents thus strong evidence for the presence of nonbaryonic "quark" matter in the interior of heavy stars.

Maximum mass of neutron stars

We determine the structure of neutron stars within a Brueckner-Hartree-Fock approach based on realistic nucleon-nucleon, nucleon-hyperon, and hyperon-hyperon interactions. Our results indicate rather low maximum masses below 1.4 solar masses. This feature is insensitive to the nucleonic part of the EOS due to a strong compensation mechanism caused by the appearance of hyperons and represents thus strong evidence for the presence of nonbaryonic "quark" matter in the interior of heavy stars.

On the mass distribution of neutron stars

The distribution of masses for neutron stars is analysed using the Bayesian statistical inference, evaluating the likelihood of the proposed Gaussian peaks by using 54 measured points obtained in a variety of systems. The results strongly suggest the existence of a bimodal distribution of the masses, with the first peak around 1.37 M(circle dot) and a much wider second peak at 1.73 M(circle dot). The results support earlier views related to the different evolutionary histories of the members for the first two peaks, which produces a natural separation (even if no attempt to `label` the systems has been made here). They also accommodate the recent findings of similar to M(circle dot) masses quite naturally. Finally, we explore the existence of a subgroup around 1.25 M(circle dot), finding weak, if any, evidence for it. This recently claimed low-mass subgroup, possibly related to the O-Mg-Ne core collapse events, has a monotonically decreasing likelihood and does not stand out clearly from the rest of the sample.

Pion mass effects on axion emission from neutron stars through NN bremsstrahlung processes

The rates of axion emission by nucleon-nucleon bremsstrahlung are calculated with the inclusion of the full momentum contribution from a nuclear one pion exchange (OPE) potential. The contributions of the neutron-neutron (nn), proton-proton (pp) and neutron-proton (np) processes in both the non-degenerate and degenerate limits are explicitly given. We find that the finite-momentum corrections to the emissivities are quantitatively significant for the non-degenerate regime and temperature-dependent, and should affect the existing axion mass hounds. The trend of these nuclear effects is to diminish the emissivities. (C) 2009 Elsevier B.V. All rights reserved.

Evolution of magnetic field and spin period in accreting neutron stars

Based on the accretion-induced magnetic field decay model, in which a frozen field and an incompressible fluid are assumed, we obtain the following results: (1) an analytic relation between the magnetic field and spin period, if the fastness parameter of the accretion disk is neglected (The evolutionary tracks of accreting neutron stars in the P-B diagram in our model are different from the equilibrium period lines when the influence of the fastness parameter is taken into account.); (2) the theoretical minimum spin period of an accreting neutron star is max(1.1ms (DeltaM/M(circle dot))(-1)R(6)(-5/14) I(45)(M/M(circle dot))(-1/2),1.1ms (M/M(circle dot))(-1/2) R(6)(17/14)), independent of the accretion rate (X-ray luminosity) but dependent on the total accretion mass, DeltaM; however, the minimum magnetic field depends on the accretion rate; (3) the magnetic field strength decreases faster with time than does the period.

The saturated magnetic field of the neutron stars

Considering the ferromagnetic screening for the decay of the X-ray neutron star magnetic field in the binary accretion phase, the phase transition of ferromagnetic materials in the crust of neutron star induces the ferromagnetic screening saturation of the accreted crust, which results in the minimum surface magnetic field of the accreting neutron star, about 108 G, if the accreted matter has completely replaced the crust mass of the neutron star. The magnetic field evolution versus accreted mass is given as Bs ∝ ΔM-0.9, and the obtained magnetic field versus spin period relation is consistent with the distribution of the binary X-ray sources and recycled pulsars. The further thermal effect on the magnetic evolution is also studied.

Effects of vacuum structure on neutron stars

We study the equation of state for neutron matter using the Walecka model including quantum corrections for baryons and sigma mesons through a realignment of the vacuum. We next use this equation of state to calculate the radius, mass, and other properties of rotating neutron stars.

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.

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 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.

Spectral features in isolated neutron stars induced by inhomogeneous surface temperatures

The thermal X-ray spectra of several isolated neutron stars display deviations from a pure blackbody. The accurate physical interpretation of these spectral features bears profound implications for our understanding of the atmospheric composition, magnetic field strength and topology, and equation of state of dense matter. With specific details varying from source to source, common explanations for the features have ranged from atomic transitions in the magnetized atmospheres or condensed surface, to cyclotron lines generated in a hot ionized layer near the surface. Here, we quantitatively evaluate the X-ray spectral distortions induced by inhomogeneous temperature distributions of the neutron star surface. To this aim, we explore several surface temperature distributions, we simulate their corresponding general relativistic X-ray spectra (assuming an isotropic, blackbody emission), and fit the latter with a single blackbody model. We find that, in some cases, the presence of a spurious ‘spectral line’ is required at a high significance level in order to obtain statistically acceptable fits, with central energy and equivalent width similar to the values typically observed. We also perform a fit to a specific object, RX J0806.4−4123, finding several surface temperature distributions able to model the observed spectrum. The explored effect is unlikely to work in all sources with detected lines, but in some cases it can indeed be responsible for the appearance of such lines. Our results enforce the idea that surface temperature anisotropy can be an important factor that should be considered and explored also in combination with more sophisticated emission models like atmospheres.