The Moving Hotspots model for kHz QPOs in accreting neutron stars

3D MHD simulation of accretion onto neutron stars have shown in the last few years that the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move around the magnetic pole, while instabilities produce other hotspots that appear randomly and move along the equator or slightly above. The angular velocities of the two hotspots are different, the equatorial one being higher and both close to the Keplerian velocity in the inner region. Modeling of the lightcurves of these hotspots with Monte Carlo simulations show that the signatures produced in power specra by them, if observed, are QPOs plus low frequency components. Their frequencies, general behavior and features describe correctly most of the properties of kHz QPOs, if we assume the funnel flow hotspots as the origin of the lower kHz QPO and instabilities as the origin of the upper kHz QPO.Comment: 4 pages, 2 figures, to appear in the Proceedings of the ASTRONS 2010 conferenc

No Time for Dead Time: Use the Fourier Amplitude Differences to Normalize Dead-time-affected Periodograms

Dead time affects many of the instruments used in X-ray astronomy, by producing a strong distortion in power density spectra. This can make it difficult to model the aperiodic variability of the source or look for quasi-periodic oscillations. Whereas in some instruments a simple a priori correction for dead-time-affected power spectra is possible, this is not the case for others such as NuSTAR, where the dead time is non-constant and long (~2.5 ms). Bachetti et al. 2015 suggested the cospectrum obtained from light curves of independent detectors within the same instrument as a possible way out, but this solution has always only been a partial one: the measured rms was still affected by dead time, because the width of the power distribution of the cospectrum was modulated by dead time in a frequency-dependent way. In this Letter we suggest a new, powerful method to normalize cospectra and, with some caveats, even power density spectra. Our approach uses the difference of the Fourier amplitudes from two independent detectors to characterize and filter out the effect of dead time. This method is crucially important for the accurate modelling of periodograms derived from instruments affected by dead time on board current missions like NuSTAR and ASTROSAT, but also future missions such as IXPEComment: 8 pages, 5 figures, Published on ApJL on 2018 January 3

3D MHD Simulations of accreting neutron stars: evidence of QPO emission from the surface

3D Magnetohydrodynamic simulations show that when matter accretes onto neutron stars, in particular if the misalignment angle is small, it does not constantly fall at a fixed spot. Instead, the location at which matter reaches the star moves. These moving hot spots can be produced both during stable accretion, where matter falls near the magnetic poles of the star, and unstable accretion, characterized by the presence of several tongues of matter which fall on the star near the equator, due to Rayleigh-Taylor instabilities. Precise modeling with Monte Carlo simulations shows that those movements could be observed as high frequency Quasi Periodic Oscillations. We performed a number of new simulation runs with a much wider set of parameters, focusing on neutron stars with a small misalignment angle. In most cases we observe oscillations whose frequency is correlated with the mass accretion rate $\dot{M}$. Moreover, in some cases double QPOs appear, each of them showing the same correlation with $\dot{M}$.Comment: 2 pages, 1 figure, to appear in the Proceedings of the Bologna x-ray conference 2009, uses aipproc.cls, aip-6s.clo,

A search for hyperluminous X-ray sources in the XMM-Newton source catalog

We present a new method to identify luminous off-nuclear X-ray sources in the outskirts of galaxies from large public redshift surveys, distinguishing them from foreground and background interlopers. Using the 3XMM-DR5 catalog of X-ray sources and the SDSS DR12 spectroscopic sample of galaxies, with the help of this off-nuclear cross-matching technique, we selected 98 sources with inferred X-ray luminosities in the range $10^{41} < L_{\rm X} < 10^{44}\,{\rm erg\,s}^{-1}$, compatible with hyperluminous X-ray objects (HLX). To validate the method, we verify that it allowed us to recover known HLX candidates such as ESO 243$-$49 HLX$-$1 and M82 X$-$1. From a statistical study, we conservatively estimate that up to $71 \pm 11$ of these sources may be fore- or background sources, statistically leaving at least 16 that are likely to be HLXs, thus providing support for the existence of the HLX population. We identify two good HLX candidates and using other publicly available datasets, in particular the VLA FIRST in radio, UKIDSS in the near-infrared, GALEX in the ultra-violet and CFHT Megacam archive in the optical, we present evidence that these objects are unlikely to be foreground or background X-ray objects of conventional types, e.g. active galactic nuclei, BL Lac objects, Galactic X-ray binaries or nearby stars. However, additional dedicated X-ray and optical observations are needed to confirm their association with the assumed host galaxies and thus secure their HLX classification.Comment: 20 pages, 6 figures, accepted to Ap

Timing of the accreting millisecond pulsar IGR~J17511--3057

{Timing analysis of Accretion-powered Millisecond Pulsars (AMPs) is a powerful tool to probe the physics of compact objects. The recently discovered \newigrj is the 12 discovered out of the 13 AMPs known. The Rossi XTE satellite provided an extensive coverage of the 25 days-long observation of the source outburst.} {Our goal is to investigate the complex interaction between the neutron star magnetic field and the accretion disk, determining the angular momentum exchange between them. The presence of a millisecond coherent flux modulation allows us to investigate such interaction from the study of pulse arrival times. In order to separate the neutron star proper spin frequency variations from other effects, a precise set of orbital ephemeris is mandatory.} {Using timing techniques, we analysed the pulse phase delays fitting differential corrections to the orbital parameters. To remove the effects of pulse phase fluctuations we applied the timing technique already successfully applied to the case of an another AMP, XTE J1807-294.} {We report a precise set of orbital ephemeris. We demonstrate that the companion star is a main sequence star. We find pulse phase delays fluctuations on the first harmonic with a characteristic amplitude of about 0.05, similar to what also observed in the case of the AMP XTE J1814-338. For the second time an AMP shows a third harmonic detected during the entire outburst. The first harmonic phase delays show a puzzling behaviour, while the second harmonic phase delays show a clear spin-up. Also the third harmonic shows a spin-up, although not highly significant (3$\sigma$ c.l.). The presence of a fourth harmonic is also reported. In the hypothesis that the second harmonic is a good tracer of the spin frequency of the neutron star, we find a mean spin frequency derivative for this source of \np{1.65(18)}{-13} Hz s$^{-1}$.} (continue ...)Comment: 9 pages, 12 figures, A&A accepted on 23/10/201

Magnetar-like activity from the central compact object in the SNR RCW103

The 6.67 hr periodicity and the variable X-ray flux of the central compact object (CCO) at the center of the SNR RCW 103, named 1E 161348-5055, have been always difficult to interpret within the standard scenarios of an isolated neutron star or a binary system. On 2016 June 22, the Burst Alert Telescope (BAT) onboard Swift detected a magnetar-like short X-ray burst from the direction of 1E 161348-5055, also coincident with a large long-term X-ray outburst. Here we report on Chandra, NuSTAR, and Swift (BAT and XRT) observations of this peculiar source during its 2016 outburst peak. In particular, we study the properties of this magnetar-like burst, we discover a hard X-ray tail in the CCO spectrum during outburst, and we study its long-term outburst history (from 1999 to July 2016). We find the emission properties of 1E 161348-5055 consistent with it being a magnetar. However in this scenario, the 6.67 hr periodicity can only be interpreted as the rotation period of this strongly magnetized neutron star, which therefore represents the slowest pulsar ever detected, by orders of magnitude. We briefly discuss the viable slow-down scenarios, favoring a picture involving a period of fall-back accretion after the supernova explosion, similarly to what is invoked (although in a different regime) to explain the "anti-magnetar" scenario for other CCOs.Comment: 6 pages, 3 figures. To be published in the Astrophysical Journal Letters; replaced to match the version accepted for publication on 2016 August 1

An XMM-Newton and NuSTAR study of IGR J18214-1318: a non-pulsating high-mass X-ray binary with a neutron star

IGR J18214-1318, a Galactic source discovered by the International Gamma-Ray Astrophysics Laboratory, is a high-mass X-ray binary (HMXB) with a supergiant O-type stellar donor. We report on the XMM-Newton and NuSTAR observations that were undertaken to determine the nature of the compact object in this system. This source exhibits high levels of aperiodic variability, but no periodic pulsations are detected with a 90% confidence upper limit of 2% fractional rms between 0.00003-88 Hz, a frequency range that includes the typical pulse periods of neutron stars (NSs) in HMXBs (0.1-10$^3$ s). Although the lack of pulsations prevents us from definitively identifying the compact object in IGR J18214-1318, the presence of an exponential cutoff with e-folding energy $\lesssim30$ keV in its 0.3-79 keV spectrum strongly suggests that the compact object is an NS. The X-ray spectrum also shows a Fe K$\alpha$ emission line and a soft excess, which can be accounted for by either a partial-covering absorber with $N_{\mathrm{H}}\approx10^{23}$ cm$^{-2}$ which could be due to the inhomogeneous supergiant wind, or a blackbody component with $kT=1.74^{+0.04}_{-0.05}$ keV and $R_{BB}\approx0.3$ km, which may originate from NS hot spots. Although neither explanation for the soft excess can be excluded, the former is more consistent with the properties observed in other supergiant HMXBs. We compare IGR J18214-1318 to other HMXBs that lack pulsations or have long pulsation periods beyond the range covered by our observations.Comment: 15 pages, 12 figures, 4 table

On the Statistical Properties of Cospectra

In recent years, the cross-spectrum has received considerable attention as a means of characterizing the variability of astronomical sources as a function of wavelength. The cospectrum has only recently been understood as a means of mitigating instrumental effects dependent on temporal frequency in astronomical detectors, as well as a method of characterizing the coherent variability in two wavelength ranges on different timescales. In this paper, we lay out the statistical foundations of the cospectrum, starting with the simplest case of detecting a periodic signal in the presence of white noise, under the assumption that the same source is observed simultaneously in independent detectors in the same energy range. This case is especially relevant for detecting faint X-ray pulsars in detectors heavily affected by instrumental effects, including NuSTAR, Astrosat, and IXPE, which allow for even sampling and where the cospectrum can act as an effective way to mitigate dead time. We show that the statistical distributions of both single and averaged cospectra differ considerably from those for standard periodograms. While a single cospectrum follows a Laplace distribution exactly, averaged cospectra are approximated by a Gaussian distribution only for more than ~30 averaged segments, dependent on the number of trials. We provide an instructive example of a quasi-periodic oscillation in NuSTAR and show that applying standard periodogram statistics leads to underestimated tail probabilities for period detection. We also demonstrate the application of these distributions to a NuSTAR observation of the X-ray pulsar Hercules X-1

A multi-wavelength pipeline for pulsar searches

Pulsar studies in the recent years have shown, more than others, to have benefited from a multi-wavelength approach. The INAF - Astronomical Observatory in Cagliari (INAF-OAC) is a growing facility with a young group devoted to pulsar and fast transients studies across the electromagnetic spectrum. Taking advantage of this expertise we have worked to provide a suite of multi-wavelength software and databases for the observations of pulsars and compact Galactic objects at the Sardinia Radio Telescope (SRT). In turn, radio pulsar observations at SRT will be made available, in a processed format, to gamma-ray searches using AGILE and Fermi gamma-ray satellite and, in a near future, they will be complementary to polarimetric X-ray observations with IXPE.Comment: Accepted for publications in Rendiconti Lincei as Proceedings of "A Decade of AGILE: Results, Challenges and Prospects of Gamma-Ray Astrophysics

MHD Simulations of Magnetospheric Accretion, Ejection and Plasma-field Interaction

We review recent axisymmetric and three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetospheric accretion, plasma-field interaction and outflows from the disk-magnetosphere boundary.Comment: 11 pages, 8 figures, conference proceedings: "Physics at the Magnetospheric Boundary", Geneva, Switzerland, 25-28 June, 201