We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition ...Explorer X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius Req are, respectively, 1.34 − 0.16 + 0.15 M and 12.71 − 1.19 + 1.14 km , while the compactness GM R eq c 2 = 0.156 − 0.010 + 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.
We derive a generalization of forward fitting for X-ray spectroscopy to include linear polarization of X-ray sources, appropriate for the anticipated next generation of space-based photoelectric ...polarimeters. We show that the inclusion of polarization sensitivity requires joint fitting to three observed spectra, one for each of the Stokes parameters, I(E), U(E), and Q(E). The equations for Stokes' I(E) (the total intensity spectrum) are identical to the familiar case with no polarization sensitivity, and for which the model-predicted spectrum is obtained by a convolution of the source spectrum, F(E′), with the familiar energy response function, ϵ(E′) R(E′, E), where ϵ(E′) and R(E′, E) are the effective area and energy redistribution matrix, respectively. In addition to the energy spectrum, the two new relations for U(E) and Q(E) include the source polarization fraction and position angle versus energy, a(E′), and , respectively, and the model-predicted spectra for these relations are obtained by a convolution with the "modulated" energy response function, (E′) ϵ (E′)R(E, E′), where (E′) is the energy-dependent modulation fraction that quantifies a polarimeter's angular response to 100% polarized radiation. We present results of simulations with response parameters appropriate for the proposed PRAXyS Small Explorer observatory to illustrate the procedures and methods, and we discuss some aspects of photoelectric polarimeters with relevance to understanding their calibration and operation.
Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in ...their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star's equatorial circumferential radius Re and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are km and (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of Re and M that are consistent with those found in this work. We show that our measurements of Re and M for PSR J0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.
ABSTRACT
In 2019 November, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6−7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of ...the system with NICER, which led to the discovery of X-ray pulsations with a period of 9.3 s. We analysed simultaneous X-ray data obtained with NuSTAR and NICER, allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst, the maximum broad-band X-ray luminosity of the system reached (1–2) × 1039 erg s−1, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 M⊙ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of 3 × 1012 G. We conclude that RX J0209.6−7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6−7427, the NS appears to have only a moderate magnetic field strength.
Supermassive black holes (SMBHs; mass is greater than or approximately 10 5 times that of the Sun) are known to exist at the center of most galaxies with sufficient stellar mass. In the local ...universe, it is possible to infer their properties from the surrounding stars or gas. However, at high redshifts we require active, continuous accretion to infer the presence of the SMBHs, which often comes in the form of long-term accretion in active galactic nuclei. SMBHs can also capture and tidally disrupt stars orbiting nearby, resulting in bright flares from otherwise quiescent black holes. Here, we report on a ∼200-second x-ray quasi-periodicity around a previously dormant SMBH located in the center of a galaxy at redshift z = 0.3534. This result may open the possibility of probing general relativity beyond our local universe.
We report on the first simultaneous Neutron Star Interior Composition Explore (NICER) and Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the neutron star (NS) low-mass X-ray binary 4U ...1735−44, obtained in 2018 August. The source was at a luminosity of ∼1.8 (D/5.6 kpc)2 × 1037 erg s−1 in the 0.4-30 keV band. We account for the continuum emission with two different continuum descriptions that have been used to model the source previously. Despite the choice in continuum model, the combined passband reveals a broad Fe K line indicative of reflection in the spectrum. In order to account for the reflection spectrum we utilize a modified version of the reflection model relxill that is tailored for thermal emission from accreting NSs. Alternatively, we also use the reflection convolution model of rfxconv to model the reflected emission that would arise from a Comptonized thermal component for comparison. We determine that the innermost region of the accretion disk extends close to the innermost stable circular orbit (RISCO) at the 90% confidence level regardless of reflection model. Moreover, the current flux calibration of NICER is within 5% of the NuSTAR/FPMA(B).
We report the discovery of complex high-frequency variability during the 1998 August 27 giant flare from SGR 1900+14 using the Rossi X-Ray Timing Explorer (RXTE). We detect an -84 Hz quasi-periodic ...oscillation (QPO) during a 1 s interval beginning approximately 1 minute after the initial hard spike. The amplitude is energy-dependent, reaching a maximum of 26% (rms) for photons above 30 keV and is not detected below 11 keV, with a 90% confidence upper limit of 14% (rms). Remarkably, additional QPOs are detected in the average power spectrum of data segments centered on the rotational phase at which the 84 Hz signal was detected. Two signals, at 53.5 and 155.1 Hz, are strongly detected, while a third feature at 28 Hz is found with lower significance. These QPOs are not detected at other rotational phases. The phenomenology seen in the SGR 1900+14 flare is similar to that of QPOs recently reported by Israel et al. from the 2004 December 27 flare from SGR 1806-20, suggesting that they may have a common origin, perhaps torsional vibrations of the neutron star crust. Indeed, an association of the four frequencies (in increasing order) found in SGR 1900+14 with l = 2, 4, 7, and 13 toroidal modes appears plausible. We discuss our findings in the context of this model and show that if the stars have similar masses, then the magnetic field in SGR 1806-20 must be about twice as large as in SGR 1900+14, broadly consistent with estimates from pulse timing. We briefly discuss how mode identifications could lead to constraints on the nuclear equation of state.
We report an analysis of archival RXTE data from the 2004 December hyperflare from SGR 1806-20. In addition to the -90 Hz QPO first discovered by Israel et al., we report the detection of higher ...frequency oscillations at -150, 625, and 1840 Hz. We also find evidence of oscillations at -720, and 2384 Hz, but with lower significances. The 150 Hz QPO has a width (FWHM) of about 17 Hz, an average amplitude (rms) of 6.8% and is associated with the strongest peak in the pulse profile. The 625 Hz oscillation was detected in an average power spectrum from nine successive cycles beginning approximately 180 s after the initial hard spike. It has a width (FWHM) of -2 Hz and an average amplitude (rms) during this interval of 8.5%. We find a strong detection of the 625 Hz oscillation in a pair of successive rotation cycles beginning about 230 s after the start of the flare. In these cycles we also detect the 1840 Hz QPO. When the 625 Hz QPO is detected we also confirm the simultaneous presence of 30 and 92 Hz QPOs. The centroid frequency of the 625 Hz QPO detected with RXTE is within 1 Hz of the -626 Hz oscillation recently found in RHESSI data by Watts & Strohmayer. We argue that these new findings provide further evidence for a connection of these oscillations with global oscillation modes of neutron stars, in particular, the high-frequency signals may represent toroidal modes with at least one radial node in the crust. We discuss their implications in the context of this model and for the depth of neutron star crusts.
Thermonuclear X-ray bursts from accreting neutron stars power brief but strong irradiation of their surroundings, providing a unique way to study accretion physics. We analyze MAXI/Gas Slit Camera ...and Swift/XRT spectra of a day-long flash observed from IGR J17062-6143 in 2015. It is a rare case of recurring bursts at a low accretion luminosity of 0.15% Eddington. Spectra from MAXI, Chandra, and NuSTAR observations taken between the 2015 burst and the previous one in 2012 are used to determine the accretion column. We find it to be consistent with the burst ignition column of 5 × 1010 g cm−2, which indicates that it is likely powered by burning in a deep helium layer. The burst flux is observed for over a day, and decays as a straight power law: F ∝ t−1.15. The burst and persistent spectra are well described by thermal emission from the neutron star, Comptonization of this emission in a hot optically thin medium surrounding the star, and reflection off the photoionized accretion disk. At the burst peak, the Comptonized component disappears, when the burst may dissipate the Comptonizing gas, and it returns in the burst tail. The reflection signal suggests that the inner disk is truncated at ∼102 gravitational radii before the burst, but may move closer to the star during the burst. At the end of the burst, the flux drops below the burst cooling trend for 2 days, before returning to the pre-burst level.
We report on an analysis of XMM-Newton data from the neutron star low-mass X-ray binary (LMXB) Serpens X-1 (Ser X-1). Spectral analysis of EPIC PN data indicates that the previously known broad iron ...K alpha emission line from this source has a significantly skewed structure with a moderately extended red wing. The asymmetric shape of the line is well described with the laor and diskline models in XSPEC and strongly supports an inner accretion disk origin of the line. To our knowledge, this is the first strong evidence of a relativistic line in a neutron star LMXB. This finding suggests that the broad lines seen in other neutron star LMXBs likely originate from the inner disk as well. Detailed study of such lines opens up a new way to probe neutron star parameters and their strong gravitational fields. The red wing of the iron line from Ser X-1 is not as broad as that observed from some black hole systems. This is not unreasonable for a neutron star system, as the accretion disk has to terminate at or before the hard stellar surface. Finally, the inferred source inclination angle in the approximate range 40 degree -60 degree is consistent with the lack of dips and eclipses from Ser X-1.