This paper describes the design and implementation of stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality ...comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build allows for easy extensions and incorporation of its methods into data analysis workflows and pipelines. We aim for the library to be a platform for the implementation of future spectral-timing techniques. We describe the overall vision and framework, core functionality, extensions, and connections to high-level command-line and graphical interfaces. The code is well tested, with a test coverage of currently 95%, and is accompanied by extensive Application Program Interface (API) documentation and a set of step-by-step tutorials.
Aims. In light of recent discoveries of pulsating ultraluminous X-ray sources (ULXs) and recently introduced theoretical schemes that propose neutron stars (NSs) as the central engines of ULXs, we ...revisit the spectra of eighteen well known ULXs, in search of indications that favour this newly emerging hypothesis. Methods. We examine the spectra from high-quality XMM-Newton and NuSTAR observations. We use a combination of elementary black body and multicolour disk black body (MCD) models, to diagnose the predictions of classic and novel theoretical models of accretion onto NSs. We re-interpret the well established spectral characteristics of ULXs in terms of accretion onto lowly or highly magnetised NSs, and explore the resulting parameter space for consistency. Results. We confirm the previously noted presence of the low-energy (≲6 keV) spectral rollover and argue that it could be interpreted as due to thermal emission. The spectra are well described by a double thermal model consisting of a “hot” (≳1 keV) and a “cool” (≲0.7 keV) multicolour black body (MCB). Under the assumption that the “cool” MCD emission originates in a disk truncated at the neutron star magnetosphere, we find that all ULXs in our sample are consistent with accretion onto a highly magnetised (B ≳ 1012 G) neutron star. We note a strong correlation between the strength of the magnetic field, the temperature of the “hot” thermal component and the total unabsorbed luminosity. Examination of the NuSTAR data supports this interpretation and also confirms the presence of a weak, high-energy (≳15 keV) tail, most likely the result of modification of the MCB emission by inverse Compton scattering. We also note that the apparent high-energy tail, may simply be the result of mismodelling of MCB emission with an atypical temperature (T) versus radius (r) gradient, using a standard MCD model with a fixed gradient of T ~ r-0.75. Conclusions. We have offered a new and robust physical interpretation for the dual-thermal spectra of ULXs. We find that the best-fit derived parameters of our model, are in excellent agreement with recent theoretical predictions that favour super-critically accreting NSs as the engines of a large fraction of ULXs. Nevertheless, the considerable degeneracy between models and the lack of unequivocal evidence cannot rule out other equally plausible interpretations. Deeper broadband observations and time-resolved spectroscopy are warranted to further explore this newly emerging framework.
Abstract
Over the past few decades, the measurement precision of some pulsar timing experiments has advanced from ∼10
μ
s to ∼10 ns, revealing many subtle phenomena. Such high precision demands both ...careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present
PINT
(
P
INT
I
s
N
ot
T
empo3
), a high-precision
Python
pulsar timing data analysis package, which is hosted on GitHub and available on the
Python
Package Index (PyPI) as
pint-pulsar
.
PINT
is well tested, validated, object oriented, and modular, enabling interactive data analysis and providing an extensible and flexible development platform for timing applications. It utilizes well-debugged public
Python
packages (e.g., the N
um
P
y
and A
stropy
libraries) and modern software development schemes (e.g., version control and efficient development with
git
and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility.
PINT
is developed and implemented without referring to, copying, or transcribing the code from other traditional pulsar timing software packages (e.g.,
Tempo
/
Tempo2
) and therefore provides a robust tool for cross-checking timing analyses and simulating pulse arrival times. In this paper, we describe the design, use, and validation of
PINT
, and we compare timing results between it and
Tempo
and
Tempo2
.
Abstract
The Gas Pixel Detector (GPD) is an X-ray polarimeter to fly onboard IXPE and other missions. To correctly measure the source polarization, the response of IXPE’s GPDs to unpolarized ...radiation has to be calibrated and corrected. In this paper, we describe the way such response is measured with laboratory sources and the algorithm to apply such correction to the observations of celestial sources. The latter allows to correct the response to polarization of single photons, therefore allowing great flexibility in all the subsequent analysis. Our correction approach is tested against both monochromatic and nonmonochromatic laboratory sources and with simulations, finding that it correctly retrieves the polarization up to the statistical limits of the planned IXPE observations.
Abstract
Imaging X-ray Polarimetry Explorer (IXPE) is a Small Explorer mission that was launched at the end of 2021 to measure the polarization of X-ray emission from tens of astronomical sources. ...Its focal-plane detectors are based on the Gas Pixel Detector, which measures the polarization by imaging photoelectron tracks in a gas mixture and reconstructing their initial directions. The quality of the single track, and then the capability of correctly determining the original direction of the photoelectron, depends on many factors, e.g., whether the photoelectron is emitted at low or high inclination with respect to the collection plane or the occurrence of a large Coulomb scattering close to the generation point. The reconstruction algorithm used by IXPE to obtain the photoelectron emission direction also calculates several properties of the shape of the tracks that characterize the process. In this paper we compare several such properties and identify the best one to weight each track on the basis of the reconstruction accuracy. We demonstrate that significant improvement in sensitivity can be achieved with this approach and for this reason it will be the baseline for IXPE data analysis.
ABSTRACT The recent discovery by Bachetti et al. of a pulsar in M82 that can reach luminosities of up to 1040 erg s−1, a factor of ∼100 times the Eddington luminosity for a 1.4 M compact object, ...poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We analyze 19 ACIS observations where the point-spread function (PSF) of the pulsar is not contaminated by nearby sources. We fit the Chandra spectra of the pulsar with a power-law model and a disk blackbody model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where four observations have a pile-up fraction of >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (LX > 1039 erg s−1) is Γ = 1.33 0.15. For the disk blackbody model, the average temperature is Tin = 3.24 0.65 keV, the spectral shape being consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cutoff, where Γ = 0.6 0.3 and keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 1039 erg s−1, greater than 10 times its Eddington limit.
The first pulsating ultraluminous X-ray source (PULX) to be identified is M82 X-2. After the discovery in 2014, NuSTAR observed the M82 field 15 times throughout 2015 and 2016. In this paper, we ...report the results of pulsation searches in all of these data sets and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin-down rate between 2014 and 2016 of ∼−6 × 10−11 Hz s−1, which is in contrast to the strong spin-up seen during the 2014 observations, representing the first detection of spin-down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 data set. We find a glitch superimposed on the very strong and variable spin-up already reported-the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs.
Abstract We observed the nearby and relatively understudied ultraluminous X-ray source (ULX) NGC 4190 ULX-1 jointly with Neutron Star Interior Composition Explorer (NICER) and NuSTAR to investigate ...its broadband spectrum, timing properties, and spectral variation over time. We found NGC 4190 ULX-1 to have a hard spectrum characterized by two thermal components (with temperatures ∼0.25 and ∼1.6 keV) and a high-energy excess typical of the ULX population although the spectrum turns over at an unusually low energy. While no pulsations were detected (with pulsed fraction 3 σ upper limits of 16% for NICER and 35% for NuSTAR), the source shows significant stochastic variability, and the covariance spectrum indicates the presence of a high-energy cutoff power-law component, potentially indicative of an accretion column. Additionally, when fitting archival XMM-Newton data with a similar model, we find that the luminosity–temperature evolution of the hot thermal component follows the behavior of a super-Eddington slim disk though the expected spectral broadening for such a disk is not seen, suggesting that the inner accretion disk may be truncated by a magnetic field. Therefore, despite the lack of detected pulsations, there is tantalizing evidence for NGC 4190 ULX-1 being a candidate neutron star accretor although further broadband observations will be required to confirm this behavior.
Abstract
We present results for the first observed outburst from the transient X-ray binary source MAXI J0637–430. This study is based on eight observations from the Nuclear Spectroscopic Telescope ...Array (NuSTAR) and six observations from the Neil Gehrels Swift Observatory X-Ray Telescope (Swift/XRT) collected from 2019 November 19 to 2020 April 26 as the 3–79 keV source flux declined from 8.2 × 10
−10
to 1.4 × 10
−12
erg cm
−2
s
−1
. We see the source transition from a soft state with a strong disk-blackbody component to a hard state dominated by a power-law or thermal Comptonization component. NuSTAR provides the first reported coverage of MAXI J0637–430 above 10 keV, and these broadband spectra show that a two-component model does not provide an adequate description of the soft-state spectrum. As such, we test whether blackbody emission from the plunging region could explain the excess emission. As an alternative, we test a reflection model that includes a physical Comptonization continuum. Finally, we also test a spectral component based on reflection of a blackbody illumination spectrum, which can be interpreted as a simple approximation to the reflection produced by returning disk radiation due to the bending of light by the strong gravity of the black hole. We discuss the physical implications of each scenario and demonstrate the value of constraining the source distance.
The Nuclear Spectroscopic Telescope Array (NuSTAR) mission is the first focusing X-ray telescope in the hard X-ray (3-79 keV) band. Among the phenomena that can be studied in this energy band, some ...require high time resolution and stability: rotation-powered and accreting millisecond pulsars, fast variability from black holes and neutron stars, X-ray bursts, and more. Moreover, a good alignment of the timestamps of X-ray photons to UTC is key for multi-instrument studies of fast astrophysical processes. In this paper, we describe the timing calibration of the NuSTAR mission. In particular, we present a method to correct the temperature-dependent frequency response of the on-board temperature-compensated crystal oscillator. Together with measurements of the spacecraft clock offsets obtained during downlinks passes, this allows a precise characterization of the behavior of the oscillator. The calibrated NuSTAR event timestamps for a typical observation are shown to be accurate to a precision of ∼65 s.