Timing of high-count rate sources with the NuSTAR Small Explorer Mission requires specialized analysis techniques. NuSTAR was primarily designed for spectroscopic observations of sources with relatively low count-rates rather than for timing analysis
of bright objects. The instrumental dead time per event is relatively long (~2.5 msec), and varies by a few percent event-to-event. The most obvious effect is a distortion of the white noise level in the power density spectrum (PDS) that cannot be modeled easily with the standard techniques due to the variable nature of the dead time. In this paper, we show that it is possible to exploit the presence of two completely independent focal planes and use the cross power density spectrum to obtain a good proxy of the white noise-subtracted PDS. Thereafter, one can use a Monte Carlo approach to estimate the remaining effects of dead time, namely a frequency-dependent modulation of the variance and a frequency-independent drop of the sensitivity to variability. In this way, most of the standard timing analysis can be performed, albeit with a sacrifice in signal to noise relative to what would be achieved using more standard techniques. We apply this technique to NuSTAR observations of the black hole binaries GX 339-4, Cyg X-1 and GRS 1915+105.
We present results for two Ultraluminous X-ray Sources (ULXs), IC 342 X-1 and IC 342 X-2, using two epochs of XMM-Newton and NuSTAR observations separated by $sim$7 days. We observe little spectral or flux variability above 1 keV between epochs, with
unabsorbed 0.3--30 keV luminosities being $1.04^{+0.08}_{-0.06} times 10^{40}$ erg s$^{-1}$ for IC 342 X-1 and $7.40pm0.20 times 10^{39}$ erg s$^{-1}$ for IC 342 X-2, so that both were observed in a similar, luminous state. Both sources have a high absorbing column in excess of the Galactic value. Neither source has a spectrum consistent with a black hole binary in low/hard state, and both ULXs exhibit strong curvature in their broadband X-ray spectra. This curvature rules out models that invoke a simple reflection-dominated spectrum with a broadened iron line and no cutoff in the illuminating power-law continuum. X-ray spectrum of IC 342 X-1 can be characterized by a soft disk-like black body component at low energies and a cool, optically thick Comptonization continuum at high energies, but unique physical interpretation of the spectral components remains challenging. The broadband spectrum of IC 342 X-2 can be fit by either a hot (3.8 keV) accretion disk, or a Comptonized continuum with no indication of a seed photon population. Although the seed photon component may be masked by soft excess emission unlikely to be associated with the binary system, combined with the high absorption column, it is more plausible that the broadband X-ray emission arises from a simple thin blackbody disk component. Secure identification of the origin of the spectral components in these sources will likely require broadband spectral variability studies.
We present the results of NuSTAR and XMM-Newton observations of the two ultraluminous X-ray sources (ULX) NGC 1313 X-1 and X-2. The combined spectral bandpass of the two satellites enables us to produce the first spectrum of X-1 between 0.3 and 30 ke
V, while X-2 is not significantly detected by NuSTAR above 10 keV. The NuSTAR data demonstrate that X-1 has a clear cutoff above 10 keV, whose presence was only marginally detectable with previous X-ray observations. This cutoff rules out the interpretation of X-1 as a black hole in a standard low/hard state, and it is deeper than predicted for the downturn of a broadened iron line in a reflection-dominated regime. The cutoff differs from the prediction of a single-temperature Comptonization model. Further, a cold disk-like black body component at ~0.3 keV is required by the data, confirming previous measurements by XMM-Newton only. We observe a spectral transition in X-2, from a state with high luminosity and strong variability to a lower-luminosity state with no detectable variability, and we link this behavior to a transition from a super-Eddington to a sub-Eddington regime.
We investigate the quality factor and RMS amplitude of the lower kHz QPOs from XTE J1701-462, a unique X-ray source which was observed in both the so-called Z and atoll states. Correcting for the frequency drift of the QPO, we show that, as in all so
urces for which such a correction can be applied, the quality factor and RMS amplitude drops sharply above above a critical frequency. For XTE J1701-462 this frequency is estimated to be ~800 Hz, where the quality factor reaches a maximum of ~200 (e.g. a value consistent with the one observed from more classical systems, such as 4U~1636-536). Such a drop has been interpreted as the signature of the innermost stable circular orbit, and that interpretation is consistent with the observations we report here. The kHz QPOs in the Z state are much less coherent and lower amplitude than they are in the atoll state. We argue that the change of the QPO properties between the two source states is related to the change of the scale height of the accretion disk; a prediction of the toy model proposed by barret et al. (2007). As a by-product of our analysis, we also increased the significance of the upper kHz QPO detected in the atoll phase up to 4.8 sigma (single trial significance), and show that the frequency separation (266.5+/-13.1 Hz) is comparable with the one measured from simultaneous twin QPOs the Z phase.
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.
