[en] 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 M. Moreover, in some cases double QPOs appear, each of them showing the same correlation with M.

[en] The search for astronomical pulsed signals within noisy data in the radio band is usually performed through an initial Fourier analysis to find “candidate” frequencies and then refined through the folding of the time series using trial frequencies close to the candidate. In order to establish the significance of the pulsed profiles found at these trial frequencies, pulsed profiles are evaluated with a χ ² test to establish how much they depart from a null hypothesis where the signal is consistent with a flat distribution of noisy measurements. In high-energy astronomy, the χ ² statistic has widely been replaced by the $Z_n^2$ statistic and H-test, as they are more sensitive to extra information, such as the harmonic content of the pulsed profile. The $Z_n^2$ statistic and H-test were originally developed for use with “event data” composed of arrival times of single photons, leaving it unclear how these methods could be used in radio astronomy. In this paper, we present a version of the $Z_n^2$ statistic and H-test for pulse profiles with Gaussian uncertainties appropriate for radio or even optical pulse profiles. We show how these statistical indicators provide better sensitivity to low-significance pulsar candidates with respect to the usual χ ² method and a straightforward way to discriminate between pulse profile shapes. Moreover, they provide an additional tool for radio frequency interference rejection.

[en] Neutron stars are thought to be born rapidly rotating and then exhibit a phase of rotation-powered pulsations as they slow down to 1–10 s periods. The significant population of millisecond pulsars observed in our Galaxy is explained by the recycling concept: during an epoch of accretion from a donor star in a binary system, the neutron star is spun up to millisecond periods. However, only a few pulsars are observed during this recycling process, with relatively high rotational frequencies. Here we report the detection of an X-ray pulsar with $P_{\mathrm{s}\mathrm{p}\mathrm{i}\mathrm{n}}=1.20\mathrm{s}$ in the globular cluster B091D in the Andromeda galaxy, the slowest pulsar ever found in a globular cluster. This bright (up to 30% of the Eddington luminosity) spinning-up pulsar, persistent over the 12 years of observations, must have started accreting less than 1 Myr ago and has not yet had time to accelerate to hundreds of Hertz. The neutron star in this unique wide binary with an orbital period $P_{\mathrm{o}\mathrm{r}\mathrm{b}}=30.5\mathrm{h}\mathrm{r}$ in a 12 Gyr old, metal-rich star cluster accretes from a low-mass, slightly evolved post-main-sequence companion. We argue that we are witnessing a binary formed at a relatively recent epoch by getting a ∼0.8 $M_{\odot <!--\; ???\; -->}$ star in a dynamical interaction—a viable scenario in a massive, dense globular cluster like B091D with high global and specific stellar encounter rates. This intensively accreting non-recycled X-ray pulsar therefore provides a long-sought missing piece in the standard pulsar recycling picture.

[en] The recent discovery by Bachetti et al. of a pulsar in M82 that can reach luminosities of up to 10⁴⁰ erg s⁻¹, 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 (L_X > 10³⁹ erg s⁻¹) is Γ = 1.33 ± 0.15. For the disk blackbody model, the average temperature is T_in = 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 $E_{\mathrm{C}}={14}_{-<!--\; ???\; -->3}^{+5}$ 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 10³⁹ erg s⁻¹, greater than 10 times its Eddington limit.

[en] We present an analysis of two type-I X-ray bursts observed by NuSTAR originating from the very faint transient neutron star (NS) low-mass X-ray binary GRS 1741.9–2853 during a period of outburst in 2020 May. We show that the persistent emission can be modeled as an absorbed, Comptonized blackbody in addition to Fe Kα emission, which can be attributed to relativistic disk reflection. We measure a persistent bolometric, unabsorbed luminosity of $L_{\mathrm{b}\mathrm{o}\mathrm{l}}={7.03}_{-<!--\; ???\; -->0.05}^{+0.04}\times <!--\; ??\; -->{10}^{36}\mathrm{e}\mathrm{r}\mathrm{g}{\mathrm{s}}^{-<!--\; ???\; -->1}$, assuming a distance of 7 kpc, corresponding to an Eddington ratio of 4.5%. This persistent luminosity combined with light-curve analysis leads us to infer that the bursts were the result of pure He burning rather than mixed H/He burning. Time-resolved spectroscopy reveals that the bolometric flux of the first burst exhibits a double-peaked structure, placing the source within a small population of accreting NSs that exhibit multiple-peaked type-I X-ray bursts. We find that the second, brighter burst shows evidence for photospheric radius expansion (PRE) and that at its peak, this PRE event had an unabsorbed bolometric flux of $F_{\mathrm{p}\mathrm{e}\mathrm{a}\mathrm{k}}={2.94}_{-<!--\; ???\; -->0.26}^{+0.28}\times <!--\; ??\; -->{10}^{-<!--\; ???\; -->8}\mathrm{e}\mathrm{r}\mathrm{g}{\mathrm{c}\mathrm{m}}^{-<!--\; ???\; -->2}{\mathrm{s}}^{-<!--\; ???\; -->1}$. This yields a new distance estimate of d = 9.0 ± 0.5 kpc, assuming that this corresponds to the Eddington limit for pure He burning on the surface of a canonical NS. Additionally, we performed a detailed timing analysis that failed to find evidence for quasi-periodic oscillations or burst oscillations, and we place an upper limit of 16% on the rms variability around 589 Hz, the frequency at which oscillations have previously been reported.

[en] 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³ 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 $\lesssim <!--\; ???\; -->30\mathrm{k}\mathrm{e}\mathrm{V}$ 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α emission line and a soft excess, which can be accounted for by either a partial-covering absorber with $N_{\mathrm{H}}\approx <!--\; ???\; -->{10}^{23}$ cm⁻², which could be due to the inhomogeneous supergiant wind, or a blackbody component with $kT={1.74}_{-<!--\; ???\; -->0.05}^{+0.04}$ keV and $R_{BB}\approx <!--\; ???\; -->0.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.

[en] We report on two NuSTAR observations of GRS 1741.9-2853, a faint neutron star (NS) low-mass X-ray binary burster located 10' away from the Galactic center. NuSTAR detected the source serendipitously as it was emerging from quiescence: its luminosity was 6 × 10³⁴ erg s^–1 on 2013 July 31 and 5 × 10³⁵ erg s^–1 in a second observation on 2013 August 3. A bright, 800 s long, H-triggered mixed H/He thermonuclear Type I burst with mild photospheric radius expansion (PRE) was present during the second observation. Assuming that the luminosity during the PRE was at the Eddington level, an H mass fraction X = 0.7 in the atmosphere, and an NS mass M = 1.4 M _☉, we determine a new lower limit on the distance for this source of 6.3 ± 0.5 kpc. Combining with previous upper limits, this places GRS 1741.9-2853 at a distance of 7 kpc. Energy independent (achromatic) variability is observed during the cooling of the NS, which could result from the disturbance of the inner accretion disk by the burst. The large dynamic range of this burst reveals a long power-law decay tail. We also detect, at a 95.6% confidence level (1.7σ), a narrow absorption line at 5.46 ± 0.10 keV during the PRE phase of the burst, reminiscent of the detection by Waki et al. We propose that the line, if real, is formed in the wind above the photosphere of the NS by a resonant Kα transition from H-like Cr gravitationally redshifted by a factor 1 + z = 1.09, corresponding to a radius range of 29.0-41.4 km for a mass range of 1.4-2.0 M _☉

[en] We report on the first hard X-ray detection of the Geminga pulsar above 10 keV using a 150 ks observation with the Nuclear Spectroscopic Telescope Array (NuSTAR) observatory. The double-peaked pulse profile of non-thermal emission seen in the soft X-ray band persists at higher energies. Broadband phase-integrated spectra over the 0.2-20 keV band with NuSTAR and archival XMM-Newton data do not fit to a conventional two-component model of a blackbody plus power law, but instead exhibit spectral hardening above ∼5 keV. We find that two spectral models fit the data well: (1) a blackbody (kT ₁ ∼ 42 eV) with a broken power law (Γ₁ ∼ 2.0, Γ₂ ∼ 1.4 and E _break ∼ 3.4 keV) and (2) two blackbody components (kT ₁ ∼ 44 eV and kT ₂ ∼ 195 eV) with a power-law component (Γ ∼ 1.7). In both cases, the extrapolation of the Rayleigh-Jeans tail of the thermal component is consistent with the UV data, while the non-thermal component overpredicts the near-infrared data, requiring a spectral flattening at E ∼ 0.05-0.5 keV. While strong phase variation of the power-law index is present below ∼5 keV, our phase-resolved spectroscopy with NuSTAR indicates that another hard non-thermal component with Γ ∼ 1.3 emerges above ∼5 keV. The spectral hardening in non-thermal X-ray emission as well as spectral flattening between the optical and X-ray bands argue against the conjecture that a single power law may account for multi-wavelength non-thermal spectra of middle-aged pulsars.

[en] We present the results of NuSTAR and XMM-Newton observations of the two ultraluminous X-ray sources: 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 keV, 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 blackbody 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.

[en] 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 ∼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.06^+0.08×10⁴⁰ erg s^–1 for IC 342 X-1 and 7.40 ± 0.20 × 10³⁹ 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 blackbody 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