The problem of computing the pulse profiles from thermally emitting spots on the surface of a neutron star in general relativity is reconsidered. We show that it is possible to extend Beloborodov (2002) approach to include (multiple) spots of finite size in different positions on the star surface. Results for the pulse profiles are expressed by comparatively simple analytical formulas which involve only elementary functions.
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 (M7), 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 EPIC cameras onboard the XMM-Newton satellite allow 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 M7, but seen at greater distances and thus undergoing higher interstellar absorptions. Identifications in more than 170 astronomical catalogues and visual screening allowed 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 CVs and AGN, 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.
Context: Accretion-powered X-ray pulsars show highly energy-dependent and complex pulse-profile morphologies. Significant deviations from the average pulse profile can appear, in particular close to the cyclotron line energies. These deviations can be described as energy-dependent phase lags, that is, as energy-dependent shifts of main features in the pulse profile. Aims: Using a numerical study we explore the effect of cyclotron resonant scattering on observable, energy-resolved pulse profiles. Methods: We generated the observable emission as a function of spin phase, using Monte Carlo simulations for cyclotron resonant scattering and a numerical ray-tracing routine accounting for general relativistic light-bending effects on the intrinsic emission from the accretion columns. Results: We find strong changes in the pulse profile coincident with the cyclotron line energies. Features in the pulse profile vary strongly with respect to the average pulse profile with the observing geometry and shift and smear out in energy additionally when assuming a non-static plasma. Conclusions: We demonstrate how phase lags at the cyclotron energies arise as a consequence of the effects of angular redistribution of X-rays by cyclotron resonance scattering in a strong magnetic field combined with relativistic effects. We also show that phase lags are strongly dependent on the accretion geometry. These intrinsic effects will in principle allow us to constrain a systems accretion geometry.
The isolated neutron star (INS) 2XMM J104608.7-594306 is one of the only two to be discovered through their thermal emission since the ROSAT era. In a first dedicated XMM-Newton observation of the source, we found intriguing evidence of a very fast spin period. We re-observed 2XMM J104608.7-594306 with XMM-Newton to better characterise the spectral energy distribution of the source, confirm the candidate spin period, and possibly constrain the pulsar spin-down. Statistically acceptable spectral fits and meaningful physical parameters for the source are only obtained when the purely thermal spectrum is modified by at least one line in absorption. The implied distance is consistent with a location in (or in front of) the Carina nebula, and radiation radii are compatible with emission originating on most of the surface. Non-thermal X-ray emission is ruled out at levels above 0.5% of the source luminosity. Unfortunately, the second XMM-Newton observation proved inconclusive in terms of confirming (discarding) the fast candidate spin, providing an upper limit on the pulsed fraction of the source that is very close to the limiting sensitivity for detecting the modulation found previously. In the absence of an unambiguous period determination and an estimate of the magnetic field, the nature of the source remains open to interpretation. Its likely association with the Carina cluster and its overall spectral properties (only shared by a handful of other peculiar INSs) disfavour a standard evolutionary path, or one in which the source was previously recycled by accretion in a binary system. The INS 2XMM J104608.7-594306 may be similar to Calvera (1RXS J141256.0+792204), a neutron star for which the scenario of an evolved `anti-magnetar has been discussed. A better age estimate and deeper radio and gamma-ray limits are required to further constrain the evolutionary state of the neutron star.
The fastest-spinning neutron stars in low-mass X-ray binaries, despite having undergone millions of years of accretion, have been observed to spin well below the Keplerian break-up frequency. We simulate the spin evolution of synthetic populations of accreting neutron stars in order to assess whether gravitational waves can explain this behaviour and provide the distribution of spins that is observed. We model both persistent and transient accretion and consider two gravitational-wave-production mechanisms that could be present in these systems: thermal mountains and unstable $r$-modes. We consider the case of no gravitational-wave emission and observe that this does not match well with observation. We find evidence for gravitational waves being able to provide the observed spin distribution; the most promising mechanisms being a permanent quadrupole, thermal mountains and unstable $r$-modes. However, based on the resultant distributions alone it is difficult to distinguish between the competing mechanisms.
Previous XMM-Newton observations of the thermally emitting isolated neutron star RX J1605.3+3249 provided a candidate for a shallow periodic signal and evidence of a fast spin down, which suggested a high dipolar magnetic field and an evolution from a magnetar. We obtained a large programme with XMM-Newton to confirm its candidate timing solution, understand the energy-dependent amplitude of the modulation, and investigate the spectral features of the source. We performed extensive high-resolution and broadband periodicity searches in the new observations, using the combined photons of the three EPIC cameras and allowing for moderate changes of pulsed fraction and the optimal energy range for detection. A deep $4sigma$ upper limit of $1.33(6)%$ for modulations in the relevant frequency range conservatively rules out the candidate period previously reported. Blind searches revealed no other periodic signal above the $1.5%$ level $(3sigma$) in any of the four new observations. While theoretical models fall short at physically describing the complex energy distribution of the source, best-fit X-ray spectral parameters are obtained for a fully or partially ionized neutron star hydrogen atmosphere model with $B=10^{13}$ G, modified by a broad Gaussian absorption line at $385pm10$ eV. The deep limits from the timing analysis disfavour equally well-fit double temperature blackbody models where both the star surface and small hotspots contribute to the X-ray flux of the source. We identified a low significance ($1sigma$) temporal trend on the parameters of the source in the analysis of RGS data dating back to 2002, which may be explained by unaccounted calibration issues and spectral model uncertainties. The new dataset also shows no evidence of the previously reported narrow absorption feature at $sim570$ eV, whose possible transient nature disfavours an atmospheric origin.