Search

The Neutron Star Interior Composition Explorer (NICER) has extensively monitored the August 2019 outburst of the 401 Hz millisecond X-ray pulsar SAX J1808.4-3658. In this Letter, we report on the detection of a bright helium-fueled Type I X-ray burst. With a bolometric peak flux of $(2.3\pm0.1)\times 10^{-7}$ erg/cm^2/s, this was the brightest X-ray burst among all bursting sources observed with NICER to date. The burst shows a remarkable two-stage evolution in flux, emission lines at $1.0$ keV and $6.7$ keV, and burst oscillations at the known pulsar spin frequency, with $\approx4$\% fractional sinusoidal amplitude. We interpret the burst flux evolution as the detection of the local Eddington limits associated with the hydrogen and helium layers of the neutron star envelope. The emission lines are likely associated with Fe, due to reprocessing of the burst emission in the accretion disk., Comment: 10 pages, 5 figures, submitted to ApJ Letters

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 (NICER) 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 $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{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., Comment: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State; 76 pages, 24 figures, 7 tables, 8 figure sets (available in the online journal or from the authors)

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 of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ 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 $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density., Comment: 49 pages, 16 figures, part of The Astrophysical Journal Letters focus issue on the Neutron Star Interior Composition Explorer

We report on a coherent timing analysis of the 163 Hz accreting millisecond X-ray pulsar IGR J17062-6143. Using data collected with the Neutron Star Interior Composition Explorer and XMM-Newton, we investigated the pulsar evolution over a timespan of four years. We obtained a unique phase-coherent timing solution for the stellar spin, finding the source to be spinning up at a rate of $(3.77\pm0.09)\times 10^{-15}$ Hz/s. We further find that the $0.4-6$ keV pulse fraction varies gradually between 0.5% and 2.5% following a sinusoidal oscillation with a $1210\pm40$ day period. Finally, we supplemented this analysis with an archival Rossi X-ray Timing Explorer observation, and obtained a phase coherent model for the binary orbit spanning 12 years, yielding an orbital period derivative measurement of $(8.4\pm2.0) \times 10^{-12}$ s/s. This large orbital period derivative is inconsistent with a binary evolution that is dominated by gravitational wave emission, and is suggestive of highly non-conservative mass transfer in the binary system., 16 pages, 7 figures, 5 tables. Accepted for publication in ApJ

We report the detection of 376.05 Hz (2.66 ms) coherent X-ray pulsations in NICER observations of a transient outburst of the low-mass X-ray binary IGR J17494-3030 in 2020 October/November. The system is an accreting millisecond X-ray pulsar in a 75 minute ultracompact binary. The mass donor is most likely a $\simeq 0.02 M_\odot$ finite-entropy white dwarf composed of He or C/O. The fractional rms pulsed amplitude is 7.4%, and the soft (1-3 keV) X-ray pulse profile contains a significant second harmonic. The pulsed amplitude and pulse phase lag (relative to our mean timing model) are energy-dependent, each having a local maximum at 4 keV and 1.5 keV, respectively. We also recovered the X-ray pulsations in archival 2012 XMM-Newton observations, allowing us to measure a long-term pulsar spin-down rate of $\dot\nu = -2.1(7)\times10^{-14}$ Hz/s and to infer a pulsar surface dipole magnetic field strength of $\simeq 10^9$ G. We show that the mass transfer in the binary is likely non-conservative, and we discuss various scenarios for mass loss from the system., Comment: 10 pages, 3 figures, and 1 table; accepted for publication in ApJL

Accretion disks around neutron stars regularly undergo sudden strong irradiation by Type I X-ray bursts powered by unstable thermonuclear burning on the stellar surface. We investigate the impact on the disk during one of the first X-ray burst observations with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. The burst is seen from Aql X-1 during the hard spectral state. In addition to thermal emission from the neutron star, the burst spectrum exhibits an excess of soft X-ray photons below 1 keV, where NICER's sensitivity peaks. We interpret the excess as a combination of reprocessing by the strongly photoionized disk and enhancement of the pre-burst persistent flux, possibly due to Poynting Robertson drag or coronal reprocessing. This is the first such detection for a short sub-Eddington burst. As these bursts are observed frequently, NICER will be able to study how X-ray bursts affect the disk and corona for a range of accreting neutron star systems and disk states., 6 pages, 5 figures, published by ApJ Letters

We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of mHz X-ray brightness oscillations from the "clocked burster." NICER observed the source in the periods 2017 June 20 - 29, July 11 - 13, and September 9 - 15, for a total useful exposure of 34 ks. Two consecutive dwells obtained on 2017 September 9 revealed highly significant oscillations at a frequency of 8 mHz. The fractional, sinusoidal modulation amplitude increases from 0.7 % at 1 keV to approximately 2 % at 6 keV. Similar oscillations were also detected at lower significance in three additional dwells. The oscillation frequency and amplitude are consistent with those of mHz QPOs reported in other accreting neutron star systems. A thermonuclear X-ray burst was also observed on 2017 June 22. The burst properties and X-ray colors are both consistent with the source being in a soft spectral state during these observations, findings that are confirmed by ongoing monitoring with MAXI and SWIFT-BAT. Assuming that the mHz oscillations are associated with black body emission from the neutron star surface, modeling of the phase-resolved spectra shows that the oscillation is consistent with being produced by modulation of the temperature component of this emission. In this interpretation, the black body normalization, proportional to the emitting surface area, is consistent with being constant through the oscillation cycle. We place the observations in the context of the current theory of marginally stable burning and briefly discuss the potential for constraining neutron star properties using mHz oscillations., 26 pages, 10 figures, accepted for publication in the Astrophysical Journal

Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin-frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation (QPO) near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L~1E+38 erg/s, the pulse profiles transition from single-peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turn-over to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ~1E+13 G., 35 pages, 7 figures (main text), 5 figures (appendix); Accepted for publication in ApJ on June 19, 2018

We present Neutron Star Interior Composition Explorer (NICER) observations of the neutron star low-mass X-ray binary Serpens X-1 during the early mission phase in 2017. With the high spectral sensitivity and low-energy X-ray passband of NICER, we are able to detect the Fe L line complex in addition to the signature broad, asymmetric Fe K line. We confirm the presence of these lines by comparing the NICER data to archival observations with XMM-Newton/RGS and NuSTAR. Both features originate close to the innermost stable circular orbit (ISCO). When modeling the lines with the relativistic line model RELLINE, we find the Fe L blend requires an inner disk radius of $1.4_{-0.1}^{+0.2}$ $R_{\mathrm{ISCO}}$ and Fe K is at $1.03_{-0.03}^{+0.13}$ $R_{\mathrm{ISCO}}$ (errors quoted at 90%). This corresponds to a position of $17.3_{-1.2}^{+2.5}$ km and $12.7_{-0.4}^{+1.6}$ km for a canonical neutron star mass ($M_{\mathrm{NS}}=1.4\ M_{\odot}$) and dimensionless spin value of $a=0$. Additionally, we employ a new version of the RELXILL model tailored for neutron stars and determine that these features arise from a dense disk and supersolar Fe abundance., Comment: 6 pages, 4 figures, 1 table, accepted for publication in ApJL

We report the discovery by the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) of the accreting millisecond X-ray pulsar IGR J17591-2342. Coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation were detected. This implies an orbital period of ∼8.8 h and a projected semi-major axis of ∼1.23 lt-s. With the binary mass function, we estimate a minimum companion mass of 0.42 M, obtained assuming a neutron star mass of 1.4[subscript ⊙] and an inclination angle lower than 60°, as suggested by the absence of eclipses or dips in the light curve of the source. The broad-band energy spectrum, obtained by combining NuSTAR, swift and INTEGRAL observations, is dominated by Comptonisation of soft thermal seed photons with a temperature of ∼0.7 keV by electrons heated to 21 keV. We also detect black-body-like thermal direct emission that is compatible with an emission region of a few kilometers and a temperature compatible with the seed source of Comptonisation. A weak Gaussian line centred on the iron Kα complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, which was measured when the outburst faded. Key words: accretion, accretion disks / stars: low-mass / pulsars: general / stars: neutron / X-rays: binaries, United States. National Aeronautics and Space Administration

We present results of recent Neutron Star Interior Composition Explorer observations of the accreting millisecond X-ray pulsar IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38 minute orbital period. NICER observed the source for approximately 26 ksec over a 5.3 day span in 2017 August, and again for 14 and 11 ksec in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hz in Rossi X-ray Timing Explorer data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z^2 method finds a best-fitting circular orbit with a period of 2278.21 s (37.97 min), a projected semi-major axis of 0.00390 lt-sec, and a barycentric pulsar frequency of 163.6561105 Hz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 e-8} solar masses, presently the smallest known for a stellar binary. The minimum donor mass ranges from about 0.005 - 0.007 solar masses, for a neutron star mass from 1.2 - 2 solar masses. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175 - 0.0155 solar masses and 19 deg < i < 27.5 deg, where the lower and upper bounds correspond to 1.4 and 2 solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 +- 0.11 % (0.3 - 3.2 keV)., 19 pages, 7 figures, published in the Astrophysical Journal Letters

The Neutron Star Interior Composition Explorer (NICER) on the International Space Station (ISS) observed strong photospheric expansion of the neutron star in 4U 1820-30 during a Type I X-ray burst. A thermonuclear helium flash in the star's envelope powered a burst that reached the Eddington limit. Radiation pressure pushed the photosphere out to ~200 km, while the blackbody temperature dropped to 0.45 keV. Previous observations of similar bursts were performed with instruments that are sensitive only above 3 keV, and the burst signal was weak at low temperatures. NICER's 0.2-12 keV passband enables the first complete detailed observation of strong expansion bursts. The strong expansion lasted only 0.6 s, and was followed by moderate expansion with a 20 km apparent radius, before the photosphere finally settled back down at 3 s after the burst onset. In addition to thermal emission from the neutron star, the NICER spectra reveal a second component that is well fit by optically thick Comptonization. During the strong expansion, this component is six times brighter than prior to the burst, and it accounts for 71% of the flux. In the moderate expansion phase, the Comptonization flux drops, while the thermal component brightens, and the total flux remains constant at the Eddington limit. We speculate that the thermal emission is reprocessed in the accretion environment to form the Comptonization component, and that changes in the covering fraction of the star explain the evolution of the relative contributions to the total flux., 6 pages, 4 figures, 1 table, published by ApJ Letters

We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar Swift J1756.9$-$2508 during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combining our NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we also studied the long-term spin and orbital evolution of this source. We find that the binary system is well described by a constant orbital period model, with an upper limit on the orbital period derivative of $|\dot P_b| < 7.4\times 10^{-13}$ s s$^{-1}$. Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times, finding R.A. = $17^{\rm h}56^{\rm m}57.18^{\rm s}\pm0.08^{\rm s}$ and Decl. = $-25^{\circ}06'27.8''\pm3.5''$, while simultaneously measuring the long-term spin frequency derivative as $\dot��= -7.3\times 10^{-16}$ Hz s$^{-1}$. We briefly discuss the implications of these measurements in the context of the wider population of accreting millisecond pulsars., 9 pages, 4 figures, 2 tables. Accepted for publication in ApJ

We report on the detection of a kilohertz quasi-periodic oscillation (QPO) with the Neutron Star Interior Composition Explorer (NICER). Analyzing approximately 165 ks of NICER exposure on the X-ray burster 4U 0614+09, we detect multiple instances of a single-peak upper kHz QPO, with centroid frequencies that range from 400 Hz to 750 Hz. We resolve the kHz QPO as a function of energy, and measure, for the first time, the QPO amplitude below 2 keV. We find the fractional amplitude at 1 keV is on the order of 2% rms, and discuss the implications for the QPO emission process in the context of Comptonization models., Comment: proof correction applied. 7 pages, 3 figures, 1 table, accepted for publication in ApJ Letters

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 Stoke's parameters, I(E), U(E), and Q(E). The equations for Stoke's 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(E')*R(E', E), where e(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 psi'_0(E'), respectively, and the model-predicted spectra for these relations are obtained by a convolution with the "modulated" energy response function, m(E')*e(E')R(E, E'), where m(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., 27 pages, 8 figures, accepted for publication in The Astrophysical Journal

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/GSC 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 5x10^10 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 ~100 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., Comment: 16 pages, 10 figures, 6 tables, accepted for publication in ApJ

We report for the first time below 1.5 keV, the detection of a secondary peak in an Eddington-limited thermonuclear X-ray burst observed by the Neutron Star Interior Composition Explorer (NICER) from the low-mass X-ray binary 4U 1608-52. Our time-resolved spectroscopy of the burst is consistent with a model consisting of a varying-temperature blackbody, and an evolving persistent flux contribution, likely attributed to the accretion process. The dip in the burst intensity before the secondary peak is also visible in the bolometric flux. Prior to the dip, the blackbody temperature reached a maximum of $\approx3$ keV. Our analysis suggests that the dip and secondary peak are not related to photospheric expansion, varying circumstellar absorption, or scattering. Instead, we discuss the observation in the context of hydrodynamical instabilities, thermonuclear flame spreading models, and re-burning in the cooling tail of the burst., 8 pages, 6 figures, Accepted for publication in ApJ

We present evidence for three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). Our timing analysis reveals three sudden spin-down events of magnitudes $\Delta\nu = -23, -30,$ and $-43 \,\mu$Hz (fractional amplitudes $\Delta\nu/\nu = -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfluid components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse., Comment: ApJ, accepted

The discovery of quasi-periodic brightness oscillations (QPOs) in the X-ray emission accompanying the giant flares of the soft gamma-ray repeaters SGR 1806-20 and SGR 1900+14 has led to intense speculation about their nature and what they might reveal about the interiors of neutron stars. Here we take a fresh look at the giant flare data for SGR 1806-20, and in particular we analyze short segments of the post-peak emission using a Bayesian procedure that has not previously been applied to these data. We find at best weak evidence that any QPO persists for more than $\sim 1$ second; instead, almost all the data are consistent with a picture in which there are numerous independently-excited modes that decay within a few tenths of a second. This has interesting implications for the rapidity of decay of the QPO modes, which could occur by the previously-suggested mechanism of coupling to the MHD continuum. The strongest QPOs favor certain rotational phases, which might suggest special regions of the crust or of the magnetosphere. We also find several previously unreported QPOs in these data, which may help in tracking down their origin., Comment: 20 pages, 5 figures, published in The Astrophysical Journal

We report the discovery of kHz quasi periodic oscillations (QPOs) in three Rossi X-ray Timing Explorer observations of the low mass X-ray binary (LMXB) XTE J1701-407. In one of the observations we detect a kHz QPO with a characteristic frequency of 1153 +/- 5 Hz, while in the other two observations we detect twin QPOs at characteristic frequencies of 740 +/- 5 Hz, 1112 +/- 17 Hz and 740 +/- 11 Hz, 1098 +/- 5 Hz. All detections happen when XTE J1701-407 was in its high intensity soft state, and their single trial significance are in the 3.1-7.5 sigma range. The frequency difference in the centroid frequencies of the twin kHz QPOs (385 +/- 13 Hz) is one of the largest seen till date. The 3-30 keV fractional rms amplitude of the upper kHz QPO varies between ~18 % and ~30 %. XTE J1701-407, with a persistent luminosity close to 1 % of the Eddington limit, is among the small group of low luminosity kHz QPO sources and has the highest rms for the upper kHz QPO detected in any source. The X-ray spectral and variability characteristics of this source indicate its atoll source nature., Comment: 10 pages, 9 figures, accepted for publication in MNRAS

The Polarimeter for Relativistic Astrophysical X-ray Sources (PRAXyS) is one of three Small Explorer (SMEX) missions selected by NASA for Phase A study, with a launch date in 2020. The PRAXyS Observatory exploits grazing incidence X-ray mirrors and Time Projection Chamber Polarimeters capable of measuring the linear polarization of cosmic X-ray sources in the 2-10 keV band. PRAXyS combines well-characterized instruments with spacecraft rotation to ensure low systematic errors. The PRAXyS payload is developed at the Goddard Space Flight Center with the Johns Hopkins University Applied Physics Laboratory, University of Iowa, and RIKEN (JAXA) collaborating on the Polarimeter Assembly. The LEOStar-2 spacecraft bus is developed by Orbital ATK, which also supplies the extendable optical bench that enables the Observatory to be compatible with a Pegasus class launch vehicle. A nine month primary mission will provide sensitive observations of multiple black hole and neutron star sources, where theory predicts polarization is a strong diagnostic, as well as exploratory observations of other high energy sources. The primary mission data will be released to the community rapidly and a Guest Observer extended mission will be vigorously proposed.

We analyze the latest \emph{Suzaku} observation of the bright neutron star low-mass X-ray binary Serpens X-1 taken in 2013 October and 2014 April. The observation was taken using the burst mode and only suffered mild pile-up effects. A broad iron line is clearly detected in the X-ray spectrum. We test different models and find that the iron line is asymmetric and best interpreted by relativistic reflection. The relativistically broadened iron line is generally believed to originate from the innermost regions of the accretion disk, where strong gravity causes a series of special and general relativistic effects. The iron line profile indicates an inner radius of $\sim8$ $R_{\rm G}$, which gives an upper limit on the size of the neutron star. The asymmetric iron line has been observed in a number of previous observations, which gives several inner radius measurements at different flux states. We find that the inner radius of Serpens X-1 does not evolve significantly over the range of $L/L_{\rm Edd}\sim0.4-0.6$, and the lack of flux dependence of the inner radius implies that the accretion disk may be truncated outside the innermost stable circular orbit by the boundary layer rather than the stellar magnetic field., 8 pages, 6 figures, accepted for publication in ApJ

We present an analysis of the high frequency X-ray variability of NGC 4051 (M_BH ~ 1.7*10^6 M_sun) based on a series of XMM-Newton observations taken in 2009 with a total exposure of ~570 ks (EPIC pn). These data reveal the form of the power spectrum over frequencies from 10^-4 Hz, below the previously detected power spectral break, to > 10^-2 Hz, above the frequency of the innermost stable circular orbit (ISCO) around the black hole (nu_ISCO ~ 10^-3 - 10^-2 Hz, depending on the black hole spin parameter j). This is equivalent to probing frequencies >1 kHz in a stellar mass (M_BH ~ 10 Msun) black hole binary system. The power spectrum is a featureless power law over the region of the expected ISCO frequency, suggesting no strong enhancement or change in the variability at the fastest orbital period in the system. Despite the huge amplitude of the flux variations between the observations (peak-to-peak factor of >50) the power spectrum appears to be stationary in shape, and varies in amplitude at all observed frequencies following the previously established linear rms-flux relation. The rms-flux relation is offset in flux by a small and energy-dependent amount. The simplest interpretation of the offset is in terms of a very soft spectral component that is practically constant (compared to the primary source of variability). One possible origin for this emission is a circumnuclear shock energised by a radiatively driven outflow from the central regions, and emitting via inverse-Compton scattering of the central engine's optical-UV continuum (as inferred from a separate analysis of the energy spectrum). A comparison with the power spectrum of a long XMM-Newton observation taken in 2001 gives only weak evidence for non-stationarity in power spectral shape or amplitude. [Abridged]

The discovery of two stable peaks at frequencies with a ratio of 3:2 in the power spectrum of X-ray emission from the brightest X-ray source in galaxy M82 suggests that, if the relationship between frequency and mass that holds for stellar-mass black holes can be extended to intermediate masses, the black hole believed to be the source of the emission has a mass approximately 400 times that of the Sun. Stable oscillations in the X-ray emissions of stellar-mass black holes correlate closely with black hole mass. Dheeraj Pasham et al. use this relationship to calculate the mass of the brightest X-ray source in the nearby starburst galaxy M82, thought to be an intermediate-mass black hole. Previous estimates were based on scaling relations with large intrinsic uncertainties that even allowed for a stellar-mass black hole. The new figure is around 428 solar masses — or alternatively about 415 plus or minus 63 solar masses based on a relativistic precession model. These findings confirm that the millihertz quasi-periodic oscillations of M82 X-1 are from an intermediate-mass black hole. M82 X-1, the brightest X-ray source in the galaxy M82, has been thought to be an intermediate-mass black hole (100 to 10,000 solar masses) because of its extremely high luminosity and variability characteristics1,2,3,4,5,6, although some models suggest that its mass may be only about 20 solar masses3,7. The previous mass estimates were based on scaling relations that use low-frequency characteristic timescales which have large intrinsic uncertainties8,9. For stellar-mass black holes, we know that the high-frequency quasi-periodic oscillations (100–450 hertz) in the X-ray emission that occur in a 3:2 frequency ratio are stable and scale in frequency inversely with black hole mass with a reasonably small dispersion10,11,12,13,14,15. The discovery of such stable oscillations thus potentially offers an alternative and less ambiguous means of mass determination for intermediate-mass black holes, but has hitherto not been realized. Here we report stable, twin-peak (3:2 frequency ratio) X-ray quasi-periodic oscillations from M82 X-1 at frequencies of 3.32 ± 0.06 hertz and 5.07 ± 0.06 hertz. Assuming that we can extrapolate the inverse-mass scaling that holds for stellar-mass black holes, we estimate the black hole mass of M82 X-1 to be 428 ± 105 solar masses. In addition, we can estimate the mass using the relativistic precession model, from which we get a value of 415 ± 63 solar masses.

We report the discovery of burst oscillations at the spin frequency in ten thermonuclear bursts from the accreting millisecond X-ray pulsar (AMXP) IGR J17511-3057. The burst oscillation properties are, like those from the persistent AMXPs SAX J1808.4-3658 and XTE J1814-338, anomalous compared to burst oscillations from intermittent pulsars or non-pulsing LMXBs. Like SAX J1808.4-3658 they show frequency drifts in the rising phase rather than the tail. There is also evidence for harmonic content. Where IGR J17511-3057 is unusual compared to the other two persistent pulsars is that oscillations are not detected throughout all bursts. As accretion rate drops the bursts get brighter and their rise/decay time scales become shorter, while the oscillation amplitude falls below the detection threshold: first in the burst peak and then also in the rise. None of the bursts from IGR J17511-3057 show evidence for photospheric radius expansion (which might be expected to suppress oscillation amplitude) which allow us to set an upper limit to the distance of 6.9 kpc. We discuss the implications of our results for models of the burst oscillation mechanism.

We report the discovery of X-ray eclipses in the recently discovered accreting millisecond X-ray pulsar Swift J1749.4-2807. This is the first detection of X-ray eclipses in a system of this type and should enable a precise neutron star mass measurement once the companion star is identified and studied. We present a combined pulse and eclipse timing solution that enables tight constraints on the orbital parameters and inclination and shows that the companion mass is in the range 0.6-0.8 M_sun for a likely range of neutron star masses, and that it is larger than a main sequence star of the same mass. We observed two individual eclipse egresses and a single ingress. Our timing model shows that the eclipse features are symmetric about the time of 90 deg longitude from the ascending node, as expected. Our eclipse timing solution gives an eclipse duration (from the mid-points of ingress to egress) of 2172 +/- 13 s. This represents 6.85% of the 8.82 hr orbital period. This system also presents a potential measurement of "Shapiro" delay due to General Relativity; through this technique alone, we set an upper limit to the companion mass of 2.2 M_sun., Comment: Accepted to ApJ Letters

We report the discovery with XMM-Newton of correlated spectral and timing behavior in the ultraluminous X-ray source (ULX) NGC 5408 X-1. An ~100 ksec pointing with XMM/Newton obtained in January, 2008 reveals a strong 10 mHz QPO in the > 1 keV flux, as well as flat-topped, band limited noise breaking to a power law. The energy spectrum is again dominated by two components, a 0.16 keV thermal disk and a power-law with an index of ~2.5. These new measurements, combined with results from our previous January 2006 pointing in which we first detected QPOs, show for the first time in a ULX a pattern of spectral and temporal correlations strongly analogous to that seen in Galactic black hole sources, but at much higher X-ray luminosity and longer characteristic time-scales. We find that the QPO frequency is proportional to the inferred disk flux, while the QPO and broad-band noise amplitude (root mean squared, rms) are inversely proportional to the disk flux. Assuming that QPO frequency scales inversely with black hole mass at a given power-law spectral index we derive mass estimates using the observed QPO frequency - spectral index relations from five stellar-mass black hole systems with dynamical mass constraints. The results from all sources are consistent with a mass range for NGC 5408 X-1 from 1000 - 9000 solar masses. Moreover, the recent relation from Gierlinski et al. that relates black hole mass to the strength of variability at high frequencies (above the break in the power spectrum) is also indicative of such a high mass for NGC 5408 X-1. We argue that these new findings strongly support the conclusion that NGC 5408 X-1 harbors an intermediate mass black hole., Accepted for Publication in the Astrophysical Journal

RX J0806.3+1527 is a candidate double-degenerate binary with possibly the shortest known orbital period. The source shows an ≈100% X-ray intensity modulation at the putative orbital period of 321.5 s, but there still is no direct detection of orbital motion at this period. Here we report the results of high-resolution X-ray spectroscopic observations of the source, obtained with the Chandra Low Energy Transmission Grating (LETG) spectrometer. We obtained a total of ≈80 ks of exposure in two pointings separated by about 5 months. We detect the source in both the zeroth-order image and the dispersed spectrum. Detection of the 321.5 s pulsations in the zeroth-order image allows us to extract dispersed spectra as a function of the pulse phase of the X-ray modulation. The bulk of the source flux appears in the 20-60 A band. We find evidence for emission features in the 25-50 A range. The strongest feature (detected at better than 3 σ confidence) is an ≈0.6 A-wide (FWHM) emission line centered at 27 A. Weaker features at 31.6 and 46.4 A are also suggested, but they require confirmation with deeper exposures. The continuum is adequately described by an absorbed blackbody model with a temperature of 65 ± 5 eV and column depth of (2 ± 0.5) × 1020 cm−2. Interestingly, the spectrum does not show strong emission lines of typically abundant elements such as carbon, nitrogen, oxygen, and neon. We briefly discuss possible identifications for the emission features.

We report the discovery of burst oscillations at 414.7 Hz during a thermonuclear X-ray burst from the low mass X-ray binary (LMXB) 4U 0614+091 with the Burst Alert Telescope (BAT) onboard SWIFT. In a search of the BAT archive, we found two burst triggers consistent with the position of 4U 0614+091. We searched both bursts for high frequency timing signatures, and found a significant detection at 414.7 Hz during a 5 s interval in the cooling tail of the brighter burst. This result establishes the spin frequency of the neutron star in 4U 0614+091 as 415 Hz. The oscillation had an average amplitude (rms) of 14%, These results are consistent with those known for burst oscillations seen in other LMXBs. The inferred ratio of the frequency difference between the twin kHz QPOs, and the spin frequency in this source is strongly inconsistent with either 0.5 or 1, and tends to support the recent suggestions by Yin et al., and Mendez & Belloni, that the kHz QPO frequency difference may not have a strong connection to the neutron star spin frequency.

Millisecond period brightness oscillations during the intensity rise of thermonuclear X-ray bursts are likely caused by an azimuthally asymmetric, expanding burning region on the stellar surface. The time evolution of the oscillation amplitude during the intensity rise encodes information on how the thermonuclear flames spread across the stellar surface. This process depends on properties of the accreted burning layer, surface fluid motions, and the surface magnetic field structure, and thus can provide insight into these stellar properties. We present two examples of bursts from different sources that show a decrease in oscillation amplitude during the intensity rise. Using theoretical modeling, we demonstrate that the observed amplitude evolution of these bursts is not well described by a uniformly expanding circular burning region. We further show that by including in our model the salient aspects of the Coriolis force (as described by Spitkovsky, Levin, and Ushomirsky) we can qualitatively reproduce the observed evolution curves. Our modeling shows that the evolutionary structure of burst oscillation amplitude is sensitive to the nature of flame spreading, while the actual amplitude values can be very useful to constrain some source parameters., Comment: 8 pages, 4 figures, accepted for publication in ApJ Letters

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 to 60 degree is consistent with the lack of dips and eclipses from Ser X-1., Comment: 13 pages, 2 tables, 2 figures, published in ApJ Letters

The recent discovery of high frequency oscillations during giant flares from the Soft Gamma Repeaters SGR 1806-20 and SGR 1900+14 may be the first direct detection of vibrations in a neutron star crust. If this interpretation is correct it offers a novel means of testing the neutron star equation of state, crustal breaking strain, and magnetic field configuration. We review the observational data on the magnetar oscillations, including new timing analysis of the SGR 1806-20 giant flare using data from the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) and the Rossi X-ray Timing Explorer (RXTE). We discuss the implications for the study of neutron star structure and crust thickness, and outline areas for future investigation.

We report the finding of an unusual, weak precursor to a thermonuclear X-ray burst from the accreting millisecond pulsar SAX J1808.4-3658. The burst in question was observed on Oct. 19, 2002 with the Rossi X-Ray Timing Explorer (RXTE) proportional counter array (PCA). The precursor began approximately 1 s prior to the onset of a strong radius expansion burst, lasted for about 0.4 s, and exhibited strong oscillations at the 401 Hz spin frequency. Oscillations are not detected in the approximately 0.5 s interval between the precursor and the main burst. The estimated peak photon flux and energy fluence of the precursor are about 1/25, and 1/500 that of the main burst, respectively. From joint spectral and temporal modeling, we find that an expanding burning region with a relatively low temperature on the spinning neutron star surface can explain the oscillations, as well as the faintness of the precursor with respect to the main part of the burst. We dicuss some of the implications of our findings for the ignition and spreading of thermonuclear flames on neutron stars., 10 pages, 4 figures, accepted for publication in ApJ

We report the discovery with XMM-Newton of quasiperiodic variability in the 0.2 - 10 keV X-ray flux from the ultraluminous X-ray source NGC 5408 X-1. The average power spectrum of all EPIC-pn data reveals a strong 20 mHz QPO with an average amplitude (rms) of 9%, and a coherence of Q~6. In a 33 ksec time interval when the 20 mHz QPO is strongest we also find evidence for a 2nd QPO peak at 15 mHz. This is the first indication for a close pair of QPOs in a ULX source. Interestingly, the frequency ratio of this QPO pair is inconsistent with 3:2 at the 3sigma level, but is consistent with a 4:3 ratio. A powerlaw noise component with slope near 1.5 is also present below 0.1 Hz with evidence for a break to a flatter slope at about 3 mHz. The source shows substantial broadband variability, with a total amplitude (rms) of about 30% in the 0.1 - 100 mHz frequency band, and there is strong energy dependence to the variability. We discuss the implications of these findings in the context of models for ULXs, and their implications for the object's mass., Comment: 23 pages, 7 figures. Accepted for publication in the Astrophysical Journal

Superbursts are hours-long X-ray flares attributed to the thermonuclear runaway burning of carbon-rich material in the envelope of accreting neutron stars. By studying the details of the X-ray light curve, properties of carbon combustion can be determined. In particular, we show that the shape of the rise of the light curve is set by the the slope of the temperature profile left behind by the carbon flame. We analyse RXTE/PCA observations of 4U 1636-536 and separate the direct neutron star emission from evolving photoionized reflection and persistent spectral components. This procedure results in the highest quality light curve ever produced for the superburst rise and peak, and interesting behaviour is found in the tail. The rising light curve between 100 and 1000 seconds is inconsistent with the idea that the fuel burned locally and instantaneously everywhere, as assumed in some previous models. By fitting improved cooling models, we measure for the first time the radial temperature profile of the superbursting layer. We find $d\ln T/d\ln P=1/4$. Furthermore, 20% of the fuel may be left unburned. This gives a new constraint on models of carbon burning and propagation in superbursts., Comment: 9 pages, 5 figures, accepted for publication in MNRAS

Type I X-ray bursts are thermonuclear flashes observed from the surfaces of accreting neutron stars (NSs) in Low Mass X-ray Binaries. Oscillations have been observed during the rise and/or decay of some of these X-ray bursts. Those seen during the rise can be well explained by a spreading hot spot model, but large amplitude oscillations in the decay phase remain mysterious because of the absence of a clear-cut source of asymmetry. To date there have not been any quantitative studies that consistently track the oscillation amplitude both during the rise and decay (cooling tail) of bursts. Here we compute the light curves and amplitudes of oscillations in X-ray burst models that realistically account for both flame spreading and subsequent cooling. We present results for several such "cooling wake" models, a "canonical" cooling model where each patch on the NS surface heats and cools identically, or with a latitude-dependent cooling timescale set by the local effective gravity, and an "asymmetric" model where parts of the star cool at significantly different rates. We show that while the canonical cooling models can generate oscillations in the tails of bursts, they cannot easily produce the highest observed modulation amplitudes. Alternatively, a simple phenomenological model with asymmetric cooling can achieve higher amplitudes consistent with the observations., Comment: 8 pages, 7 figures, Accepted for publication in ApJ, Additional calculations and discussion compared to v1

We analyse archival Rossi X-Ray Timing Explorer (RXTE) proportional counter array (PCA) data of thermonuclear X-ray bursts from the 2002 outburst of the accreting millisecond pulsar SAX J1808.4-3658. We present evidence of nonmonotonic variations of oscillation frequency during burst rise, and correlations among the time evolution of the oscillation frequency, amplitude, and the inferred burning region area. We also discuss that the amplitude and burning region area evolutions are consistent with thermonuclear flame spreading on the neutron star surface. Based on this discussion, we infer that for the 2002 Oct. 15 thermonuclear burst, the ignition likely occured in the mid-latitudes, the burning region took ~ 0.2 s to nearly encircle the equatorial region of the neutron star, and after that the lower amplitude oscillation originated from the remaining asymmetry of the burning front in the same hemisphere where the burst ignited. Our observational findings and theoretical discussion indicate that studies of the evolution of burst oscillation properties during burst rise can provide a powerful tool to understand thermonuclear flame spreading on neutron star surfaces under extreme physical conditions., 12 pages, 1 table, 2 figures, accepted for publication in ApJ Letters

The accreting millisecond pulsar XTE J1814-338 exhibits oscillations at the known spin frequency during Type I X-ray bursts. The properties of the burst oscillations reflect the nature of the thermal asymmetry on the stellar surface. We present an analysis of the variability of the burst oscillations of this source, focusing on three characteristics: fractional amplitude, harmonic content and frequency. Fractional amplitude and harmonic content constrain the size, shape and position of the emitting region, whilst variations in frequency indicate motion of the emitting region on the neutron star surface. We examine both long-term variability over the course of the outburst, and short-term variability during the bursts. For most of the bursts, fractional amplitude is consistent with that of the accretion pulsations, implying a low degree of fuel spread. There is however a population of bursts whose fractional amplitudes are substantially lower, implying a higher degree of fuel spread, possibly forced by the explosive burning front of a precursor burst. For the first harmonic, substantial differences between the burst and accretion pulsations suggest that hotspot geometry is not the only mechanism giving rise to harmonic content in the latter. Fractional amplitude variability during the bursts is low; we cannot rule out the hypothesis that the fractional amplitude remains constant for bursts that do not exhibit photospheric radius expansion (PRE). There are no significant variations in frequency in any of the bursts except for the one burst that exhibits PRE. This burst exhibits a highly significant but small ($\approx 0.1$Hz) drop in frequency in the burst rise. The timescale of the frequency shift is slower than simple burning layer expansion models predict, suggesting that other mechanisms may be at work., Comment: 20 pages, 20 figures, accepted for publication in ApJ. Uses emulateapj.cls

We report the discovery of complex high frequency variability during the August 27, 1998 giant flare from SGR 1900+14 using the Rossi X-ray Timing Explorer (RXTE). We detect an 84 Hz oscillation (QPO) during a 1 s interval beginning approximately 1 min after the initial hard spike. The modulation 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 December 27, 2004 flare from SGR 1806-20, suggesting 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 magnetic field estimates from pulse timing., Comment: 13 Pages, 5 figures, AASTeX, accepted for publication in the Astrophysical Journal Letters

Detailed modeling of the millisecond brightness oscillations during thermonuclear bursts from low mass X-ray binaries can provide important information about neutron star structure. Until now the implementation of this idea has not been entirely successful, largely because of the negligible harmonic content in burst oscillation lightcurves. However, the recent discovery of non-sinusoidal burst oscillation lightcurves from the accreting millisecond pulsar XTE J1814-338 has changed this situation. We, therefore, for the first time, make use of this opportunity to constrain neutron star parameters. In our detailed study of the lightcurves of 22 bursts, we fit the burst oscillation lightcurves with fully general relativistic models that include light-bending and frame-dragging for lightcurve calculation, and compute numerically the structure of neutron stars using realistic equations of state. We find that for our model and parameter grid values, at the 90% confidence level, Rc^2/GM > 4.2 for the neutron star in XTE J1814-338. We also find that the photons from the thermonuclear flash come out through the layers of accreted matter under conditions consistent with Thomson scattering, and show that the secondary companion is a hydrogen burning main sequence star, with possible bloating (probably due to X-ray heating)., 32 pages, 14 figures, 4 tables, thorough revision, accepted for publication in ApJ

We report on Chandra ACIS observations of the ultracompact AM CVn binary ES Cet. This object has a 10.3 minute binary period and is the most compact of the confirmed AM CVn systems. We have, for the first time, unambiguously detected the X-ray counterpart to ES Cet. In a 20 ks ACIS-S image a point-like X-ray source is found within 1'' of the cataloged optical position. The mean count rate in ACIS-S is 0.013 s-1, and there is no strong evidence for variability. We folded the X-ray data using the optical ephemeris of Warner & Woudt, but did not detect any significant modulation. If a ~100% modulation similar to those seen in the ultracompact candidates V407 Vul and RX J0806.3+1527 were present, then we would have detected it. The upper limit (3 σ) to any modulation at the putative orbital period is ~15% (rms). We extract the first X-ray spectrum from ES Cet, and find that it is not well described by simple continuum models. We find suggestive evidence for discrete spectral components at ~470 and 890 eV that can be modeled as Gaussian emission lines. In comparison with recent X-ray detections of nitrogen and neon in another AM CVn system (GP Com), it appears possible that these features may represent emission lines from these same elements; however, deeper spectroscopy will be required to confirm this. Our best spectral model includes a blackbody continuum with kT = 0.8 keV along with the Gaussian lines. The unabsorbed 0.2-5 keV X-ray flux was ~7 × 10-14 ergs cm-2 s-1. The observed luminosity in the hard component is a small fraction of the total expected accretion luminosity, most of which should be radiated below 0.1 keV at the expected mass accretion rate for this orbital period. We discuss the implications of our results for the nature of ES Cet.

V407 Vul (RX J1914.4+2456) is a candidate double-degenerate binary with a putative 1.756 mHz (9.5 minute) orbital frequency. In a previous timing study using archival ROSAT and ASCA data, we reported evidence for an increase of this frequency at a rate consistent with expectations for gravitational radiation from a detached ultracompact binary system. Here we report the results of new Chandra timing observations, which confirm the previous indications of spin-up of the X-ray frequency and provide much tighter constraints on the frequency derivative, . We obtained with Chandra a total of 90 ks of exposure in two epochs separated in time by 11.5 months. The total time span of the archival ROSAT, ASCA, and new Chandra data is now ≈10.5 yr. This more than doubles the interval spanned by the ROSAT and ASCA data alone, providing much greater sensitivity to a frequency derivative. With the addition of the Chandra data, an increasing frequency is evident, and the mean is (7.0 ± 0.8) × 10-18 Hz s-1. Although a long-term spin-up trend is confirmed, there is excess variance in the phase timing residuals, perhaps indicating shorter timescale torque fluctuations or phase instability associated with the source of the X-ray flux. Power spectral searches for periods longward of the 9.5 minute period do not find any significant modulations; however, the sensitivity of searches in this frequency range are somewhat compromised by the dithering of the Chandra attitude. The observed spin-up is of a magnitude consistent with that expected from gravitational radiation decay; however, the factor of ≈3 variations in flux combined with the timing noise could conceivably result from accretion-induced spin-up of a white dwarf. Continued monitoring to explore correlations of torque with X-ray flux could provide a further test of this hypothesis.

We report on X-ray spectroscopic observations with XMM/Newton of the ultracompact, double white dwarf binary GP Com. With the Reflection Grating Spectrometers (RGSs), we detect the Lyα and Lyβ lines of hydrogen-like nitrogen (N VII) and neon (Ne X), as well as the helium-like triplets (N VI and Ne IX) of these same elements. All the emission lines are unresolved. These are the first detections of X-ray emission lines from a double-degenerate, AM CVn system. We detect the resonance (r) and intercombination (i) lines of the N VI triplet, but not the forbidden (f) line. The implied line ratios for N VI, R = f/i < 0.3 and G = (f + i)/r ≈ 1, combined with the strong resonance line are consistent with formation in a dense, collision-dominated plasma. Both the RGS and EPIC/MOS spectra are well fitted by emission from an optically thin thermal plasma with an emission measure proportional to (kT/6.5 keV)0.8 (model cevmkl in XSPEC). Helium, nitrogen, oxygen, and neon are required to adequately model the spectrum; however, the inclusion of sulphur and iron further improves the fit, suggesting that these elements may also be present at low abundance. We confirm in the X-rays the underabundance of both carbon and oxygen relative to nitrogen, first deduced from optical spectroscopy by Marsh et al. The average X-ray luminosity of ≈3 × 1030 ergs s-1 implies a mass accretion rate ≈ 9 × 10-13 M☉ yr-1. The implied temperature and density of the emitting plasma, combined with the presence of narrow emission lines and the low- value, are consistent with production of the X-ray emission in an optically thin boundary layer just above the surface of the white dwarf.

We report the discovery with the EPIC CCD cameras onboard XMM-Newton of a 54 mHz quasiperiodic oscillation (QPO) in the > 2 keV X-ray flux from an ultraluminous X-ray source (ULX) in the starburst galaxy M82. This is the first detection of a QPO in the X-ray flux from an extra-Galactic ULX, and confirms that the source is a compact object. Based on the QPO strength and previous Chandra observations it appears likely that the QPO is associated with the most luminous object in the central region of M82, CXOM82 J095550.2+694047, however, XMM imaging alone is not sufficient to unambiguously confirm this. The other plausible candidate is CXOM82 J095551.1+694045, however, the QPO luminosity is comparable to the peak luminosity of this object in Chandra observations, which argues against it being the source of the QPO. The QPO had a centroid frequency of 54.3 +- 0.9 mHz, a coherence of 5, and an amplitude (rms) in the 2 - 10 keV band of 8.5%. Below 0.2 Hz the power spectrum can be described by a power-law with index ~1, and amplitude (rms) of 13.5%. The X-ray spectrum requires a curving continuum, with a disk-blackbody (diskbb) at T = 3.1 keV providing an acceptable fit. A broad Fe line centered at 6.55 keV is required in all fits, but the equivalent width (EW) is sensitive to the continuum model. There is no evidence of a reflection component. The implied bolometric luminosity is 4 - 5 x 10^{40} ergs/sec. Archival Rossi X-ray Timing Explorer (RXTE) pointings at M82 also show evidence for QPOs in the 50 - 100 mHz frequency range. We discuss the implications of our findings for models of ULXs.

Oscillating Black Hole The massive black holes that reside in the centers of galaxies can occasionally capture and tidally disrupt stars that wander too close. One such tidal disruption event was detected last year by the Swift satellite. Follow-up x-ray observations analyzed by Reis et al. (p. 949 , published online 2 August; see the Perspective by McKinney ) show quasi-periodic oscillations that suggest that an accretion disk formed around the black hole shortly after the tidal disruption event. This type of oscillation is commonly seen in the x-ray light from the much lighter black holes that result from the gravitational collapse of stars, but has been seen only once in a massive black hole residing in the center of a galaxy.

We analysed the RXTE archival data of 1E 1048.1-5937 covering a time span of more than one year. The spin down rate of this source decreases by 30 percent during the observation. We could not resolve the X-ray flux variations because of contamination by Eta Carinae. We find that the level of pulse frequency fluctuations of 1E 1048.1-5937 is consistent with typical noise levels of accretion powered pulsars. Recent RXTE observations of 1E 2259+586 have shown a constant spin down with a very low upper limit on timing noise. We used the RXTE archival X-ray observations of 1E 2259+586 to show that the intrinsic X-ray luminosity times series is also stable, with an rms fractional variation of less than 15 percent. The source could have been in a quiet phase of accretion with a constant X-ray luminosity and spin down rate.

We present a study of the distribution of asymptotic burst oscillation periods for a sample of 26 bursts from 4U 1636-53 observed with the Rossi X-ray Timing Explorer (RXTE). The burst sample includes both archival and proprietary data and spans more than 4.5 years. We also present new optical light curves of V801 Arae, the optical counterpart of 4U 1636-53, obtained during 1998-2001. We use these optical data to refine the binary period measured by Augusteijn et al. (1998) to 3.7931206(152) hours. We show that a subset of ~70% of the bursts form a tightly clustered distribution of asymptotic periods consistent with a period stability of ~ 1/10,000. The tightness of this distribution, made up of bursts spanning more than 4 years in time, suggests that the underlying period is highly stable, with a time to change the period of ~30,000 years, comparable to similar numbers derived for X-ray pulsars. We investigate the period and orbital phase data for our burst sample and show that it is consistent with binary motion of the neutron star with v sin i < 55 km/sec at 90% confidence. We use this limit as well as previous radial velocity data to constrain the binary geometry and component masses in 4U 1636-53., AASTEX, 35 pages, 9 figures. Submitted to the Astrophysical Journal

We present a detailed observational and theoretical study of a ~3 hr long X-ray burst (the ``super burst'') observed by the Rossi X-ray Timing Explorer (RXTE) from the low mass X-ray binary (LMXB) 4U 1820-30. This is the longest X-ray burst ever observed from this source, and perhaps one of the longest ever observed in great detail from any source. We show that the super burst is thermonuclear in origin. The level of the accretion driven flux as well as the total energy release of ~1.5 x 10^{42} ergs indicate that helium could not be the energy source for the super burst. We outline the physics relevant to carbon production and burning on helium accreting neutron stars and present calculations of the thermal evolution and stability of a carbon layer and show that this process is the most likely explanation for the super burst. We show that for large columns of accreted carbon fuel, a substantial fraction of the energy released in the carbon burning layer is radiated away as neutrinos, and the heat that is conducted from the burning layer in large part flows inward, only to be released on timescales longer than the observed burst. Thus the energy released possibly exceeds that observed in X-rays by more than a factor of ten. Spectral analysis during the super burst reveals the presence of a broad emission line between 5.8 - 6.4 keV and an edge at 8 - 9 keV likely due to reflection of the burst flux from the inner accretion disk in 4U 1820-30. We believe this is the first time such a signature has been unambiguously detected in the spectrum of an X-ray burst.

We have modeled X-ray burst oscillations observed with the Rossi X-Ray Timing Explorer from two low-mass X-ray binaries (LMXB): 4U 1636-53 with a frequency of 580 Hz and 4U 1728-34 at a frequency of 363 Hz. We have computed least-squares fits to the oscillations observed during the rising phase of bursts using a model that includes emission from either a single circular hot spot or a pair of circular antipodal hot spots on the surface of a neutron star. We model the spreading of the thermonuclear hot spots by assuming that the hot spot angular size grows linearly with time. We calculate the flux as a function of rotational phase from the hot spots and take into account photon deflection in the relativistic gravitational field of the neutron star, assuming the exterior spacetime is the Schwarzschild metric. We find acceptable fits with our model in a χ2 sense and use these to place constraints on the compactness of the neutron stars in these sources. For 4U 1636-53, in which detection of a 290 Hz subharmonic supports the two-spot model, we find that the compactness (i.e., mass/radius ratio) is constrained to be M/R < 0.163 at 90% confidence (G = c = 1). This requires a relatively stiff equation of state (EOS) for the stellar interior. For example, if the neutron star has a mass of 1.4 M☉, then its radius must be greater than 12.8 km. Fits using a single hot spot model are not as highly constraining. We discuss the implications of our findings for recent efforts to calculate the EOS of dense nucleon matter and the structure of neutron stars.

We present the results of {\it Chandra} observations of two non-accreting millisecond pulsars, PSRs J1640$+$2224 (J1640) and J1709$+$2313 (J1709), with low inferred magnetic fields and spin-down rates in order to constrain their surface temperatures, obtain limits on the amplitude of unstable $r$-modes in them, and make comparisons with similar limits obtained for a sample of accreting low-mass X-ray binary (LMXB) neutron stars. We detect both pulsars in the X-ray band for the first time. They are faint, with inferred soft X-ray fluxes ($0.3-3$ keV) of $\approx$ $6\times10^{-15}$ and $3\times 10^{-15}$ erg cm$^{-2}$ s$^{-1}$ for J1640 and J1709, respectively. Spectral analysis assuming hydrogen atmosphere emission gives global effective temperature upper limits ($90\%$ confidence) of $3.3 - 4.3 \times 10^5$ K for J1640 and $3.6 - 4.7 \times 10^5$ K for J1709, where the low end of the range corresponds to canonical neutron stars ($M=1.4 M_{\odot}$), and the upper end corresponds to higher-mass stars ($M=2.21 M_{\odot}$). Under the assumption that $r$-mode heating provides the thermal support, we obtain dimensionless $r$-mode amplitude upper limits of $3.2 - 4.8 \times 10^{-8}$ and $1.8 - 2.8 \times 10^{-7}$ for J1640 and J1709, respectively, where again the low end of the range corresponds to lower-mass, canonical neutron stars ($M=1.4 M_{\odot}$). These limits are about an order of magnitude lower than those we derived previously for a sample of LMXBs, except for the accreting millisecond X-ray pulsar (AMXP) SAX J1808.4$-$3658, which has a comparable amplitude limit to J1640 and J1709., 9 pages, 4 figures, published in ApJ