It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below $\sim 10^{-2}$ of the Eddington limit ($L_{\rm Edd}$). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work, we use NuSTAR and NICER observations obtained during the 2018 accretion outburst of the neutron star LMXB 4U 1608-52, to study changes in the reflection spectrum. We find that the broad Fe-K$\alpha$ line and Compton hump, clearly seen during the peak of the outburst when the X-ray luminosity is $\sim 10^{37}$ erg/s ($\sim 0.05$ $L_{\rm Edd}$), disappear during the decay of the outburst when the source luminosity drops to $\sim 4.5 \times 10^{35}$ erg/s ($\sim 0.002$ $L_{\rm Edd}$). We show that this non-detection of the reflection features cannot be explained by the lower signal-to-noise at lower flux, but is instead caused by physical changes in the accretion flow. Simulating synthetic NuSTAR observations on a grid of inner disk radius, disk ionisation, and reflection fraction, we find that the disappearance of the reflection features can be explained by either increased disk ionisation ($\log \xi \geq 4.1$) or a much decreased reflection fraction. A changing disk truncation alone, however, cannot account for the lack of reprocessed Fe-K$\alpha$ emission. The required increase in ionisation parameter could occur if the inner accretion flow evaporates from a thin disk into a geometrically thicker flow, such as the commonly assumed formation of an radiatively inefficient accretion flow at lower mass accretion rates.

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.

NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.

We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.

We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.

Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.

We report on high-resolution X-ray spectroscopy of the ultracompact X-ray binary pulsar 4U 1626-67 with Chandra/HETGS acquired in 2010, two years after the pulsar experienced a torque reversal. The well-known strong Ne and O emission lines with Keplerian profiles are shown to arise at the inner edge of the magnetically-channeled accretion disk. We exclude a photoionization model for these lines based on the absence of sharp radiative recombination continua. Instead, we show that the lines arise from a collisional plasma in the inner-disk atmosphere, with $T\simeq 10^7$ K and $n_e \sim 10^{17}$ cm^(-3). We suggest that the lines are powered by X-ray heating of the optically-thick disk inner edge at normal incidence. Comparison of the line profiles in HETGS observations from 2000, 2003, and 2010 show that the inner disk radius decreased by a factor of two after the pulsar went from spin-down to spin-up, as predicted by magnetic accretion torque theory. The inner disk is well inside the corotation radius during spin-up, and slightly beyond the corotation radius during spin-down. Based on the disk radius and accretion torque measured during steady spin-up, the pulsar's X-ray luminosity is $2\times 10^{36}$ erg/s, yielding a source distance of 3.5(+0.2-0.3) kpc. The mass accretion rate is an order of magnitude larger than expected from gravitational radiation reaction, possibly due to X-ray heating of the donor. The line profiles also indicate a binary inclination of 39(+20-10) degrees, consistent with a 0.02 Msun donor star. Our emission measure analysis favors a He white dwarf or a highly-evolved H-poor main sequence remnant for the donor star, rather than a C-O or O-Ne white dwarf. The measured Ne/O ratio is 0.46+-0.14 by number. In an appendix, we show how to express the emission measure of a H-depleted collisional plasma without reference to a H number density.

In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.

We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.

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.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

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.

The Neutron Star Interior Composition Explorer (NICER) has observed seven thermonuclear X-ray bursts from the Low Mass X-ray Binary (LMXB) neutron star 4U 1728-34 from the start of the mission's operations until February of 2019. Three of these bursts show oscillations in their decaying tail with frequencies that are within 1 Hz of the previously detected burst oscillations from this source. Two of these burst oscillations have unusual properties: They have large fractional rms amplitudes of $ 48 \pm 9 \%$ and $ 46 \pm 9 \%$, and they are detected only at photon energies above 6 keV. By contrast, the third detected burst oscillation is compatible with previous observations of this source, with a fractional rms amplitude of $7.7 \pm 1.5\%$ rms in the 0.3 to 6.2 keV energy band. We discuss the implications of these large-amplitude burst oscillations, finding they are difficult to explain with the current theoretical models for X-ray burst tail oscillations.

We report on X-ray and radio observations of the ultra-compact X-ray binary 4U 1543-624 taken in August 2017 during an enhanced accretion episode. We obtained NICER monitoring of the source over a $\sim10$ day period during which target-of-opportunity observations were also conducted with Swift, INTEGRAL, and ATCA. Emission lines were measured in the NICER X-ray spectrum at $\sim0.64$ keV and $\sim6.4$ keV that correspond to O and Fe, respectively. By modeling these line components, we are able to track changes in the accretion disk throughout this period. The innermost accretion flow appears to move inwards from hundreds of gravitational radii ($R_{g}=GM/c^{2}$) at the beginning of the outburst to $<8.7$ $R_{g}$ at peak intensity. We do not detect the source in radio, but are able to place a $3\sigma$ upper limit on the flux density at $27$ $\mu$Jy beam$^{-1}$. Comparing the radio and X-ray luminosities, we find that the source lies significantly away from the range typical of black holes in the ${L}_{{r}}$-${L}_{{x}}$ plane, suggesting a neutron star (NS) primary. This adds to the evidence that NSs do not follow a single track in the ${L}_{{r}}$-${L}_{{x}}$ plane, limiting its use in distinguishing between different classes of NSs based on radio and X-ray observations alone.

We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of narrow emission and absorption lines during photospheric radius expansion (PRE) X-ray bursts from the ultracompact binary 4U 1820-30. NICER observed the source in 2017 August accumulating about 60 ks of exposure. Five thermonuclear X-ray bursts were detected of which four showed clear signs of PRE. We extracted spectra during the PRE phases and fit each to a model that includes a comptonized component to describe the accretion-driven emission, and a black body for the burst thermal radiation. The temperature and spherical emitting radius of the fitted black body are used to assess the strength of PRE in each burst. The two strongest PRE bursts (burst pair 1) had black body temperatures of approximately 0.6 keV and emitting radii of 100 km (at a distance of 8.4 kpc). The other two bursts (burst pair 2) had higher temperatures (~0.67 keV) and smaller radii (75 km). All of the PRE bursts show evidence of narrow line emission near 1 keV. By co-adding the PRE phase spectra of burst pairs 1 and, separately, 2 we find, in both co-added spectra, significant, narrow, spectral features near 1.0 (emission), 1.7 and 3.0 keV (both in absorption). Remarkably, all the fitted line centroids in the co-added spectrum of burst pair 1 appear systematically blue-shifted by a factor of $1.046 \pm 0.006$ compared to the centroids of pair 2, strongly indicative of a gravitational red-shift, a wind-induced blue-shift, or more likely some combination of both effects. The observed shifts are consistent with this scenario in that the stronger PRE bursts in pair 1 reach larger photospheric radii, and thus have weaker gravitational red-shifts, and they generate faster outflows, yielding higher blue-shifts. We discuss possible elemental identifications for the observed features in the context of recent burst-driven wind models.

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.

We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of 7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3-10 keV flux increased by 5% while the 30-60 keV flux decreased significantly by 35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical NPEX continuum with an added broad Gaussian emission component around the spectral peak at 20 keV. Taken together with the observed Pdot value obtained from Fermi/GBM, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase of the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field.

We report on the Neutron Star Interior Composition Explorer (NICER) monitoring campaign of the 468 Hz accreting millisecond X-ray pulsar IGR J17379-3747. From a detailed spectral and timing analysis of the coherent pulsations we find that they show a strong energy dependence, with soft thermal emission lagging about 640 microseconds behind the hard, Comptonized emission. Additionally, we observe uncommonly large pulse fractions, with measured amplitudes in excess of 20% sinusoidal fractional amplitude across the NICER passband and fluctuations of up to ~70%. Based on a phase-resolved spectral analysis, we suggest that these extreme properties might be explained if the source has an unusually favorable viewing geometry with a large magnetic misalignment angle. Due to these large pulse fractions, we were able to detect pulsations down to quiescent luminosities (~5 x 10^33 erg s^-1). We discuss these low-luminosity pulsations in the context of transitional millisecond pulsars.

The unknown state of matter at ultra-high density, large proton/neutron number asymmetry, and low temperature is a major long-standing problem in modern physics. Neutron stars provide the only known setting in the Universe where matter in this regime can stably exist. Valuable information about the interior structure of neutron stars can be extracted via sensitive observations of their exteriors. There are several complementary techniques that require different combinations of high time resolution, superb spectral resolution, and high spatial resolution. In the upcoming decade and beyond, measurements of the masses and radii of an ensemble of neutron stars using these techniques, based on data from multiple proposed next-generation X-ray telescopes, can produce definitive empirical constraints on the allowed dense matter equation of state.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

The immense tidal forces of massive black holes can rip apart stars that come too close to them. As the resulting stellar debris spirals inwards, it heats up and emits x-rays when near the black hole. Here, we report the discovery of an exceptionally stable 131-second x-ray quasi-periodicity from a black hole after it disrupted a star. Using a black hole mass indicated from host galaxy scaling relations implies that, (1) this periodicity originates from very close to the black hole's event horizon, and (2) the black hole is rapidly spinning. Our findings suggest that other disruption events with similar highly sensitive observations likely also exhibit quasi-periodicities that encode information about the fundamental properties of their black holes.

We report on 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, detecting coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation. This implies an orbital period of ~8.8 hours and a projected semi-major axis of ~1.23 lt-s. From the binary mass function, we estimate a minimum companion mass of 0.42 solar masses, obtained assuming a neutron star mass of 1.4 solar masses and an inclination angle lower than 60 degrees, as suggested by the absence of eclipses or dips in the light-curve of the source. The broad-band energy spectrum 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 compatible with an emission region of a few kilometers and temperature compatible with the seed source of Comptonisation. A weak Gaussian line centered on the iron K-alpha; complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, measured during the outburst fading.

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.

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.

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).

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.

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.

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\nu = -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.

We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration.

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.

We present the analysis of seven \emph{Chandra} High Energy Transmission Grating Spectrometer and six simultaneous \emph{RXTE} Proportional Counter Array observations of the persistent neutron star (NS) low-mass X-ray binary GX 13+1 on its normal and horizontal branches. Across nearly 10 years, GX 13+1 is consistently found to be accreting at $50-70$\% Eddington, and all observations exhibit multiple narrow, blueshifted absorption features, the signature of a disk wind, despite the association of normal and horizontal branches with jet activity. A single absorber with standard abundances cannot account for all seven major disk wind features, indicating that multiple absorption zones may be present. Two or three absorbers can produce all of the absorption features at their observed broadened widths and reveal that multiple kinematic components produce the accretion disk wind signature. Assuming the most ionized absorber reflects the physical conditions closest to the NS, we estimate a wind launching radius of $7\times10^{10}$ cm, for an electron density of $10^{12}$ cm$^{-3}$. This is consistent with the Compton radius and also with a thermally driven wind. Because of the source's high Eddington fraction, radiation pressure likely facilitates the wind launching.

We report on a spectral-timing analysis of the neutron star low-mass X-ray binary Aql~X-1 with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. Aql~X-1 was observed with NICER during a dim outburst in 2017 July, collecting approximately $50$ ks of good exposure. The spectral and timing properties of the source correspond to that of an (hard) extreme island state in the atoll classification. We find that the fractional amplitude of the low frequency ($<0.3$ Hz) band-limited noise shows a dramatic turnover as a function of energy: it peaks at 0.5 keV with nearly 5\% rms, drops to $12\%$ rms at 2 keV, and rises to $15\%$ rms at 10 keV. Through the analysis of covariance spectra, we demonstrate that band-limited noise exists in both the soft thermal emission and the power-law emission. Additionally, we measure hard time lags, indicating the thermal emission at $0.5$ keV leads the power-law emission at 10 keV on a timescale of $\sim100$ ms at $0.3$ Hz to $\sim10$ ms at $3$ Hz. Our results demonstrate that the thermal emission in the hard state is intrinsically variable, and driving the modulation of the higher energy power-law. Interpreting the thermal spectrum as disk emission, we find our results are consistent with the disk propagation model proposed for accretion onto black holes.

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.

We report on the kinematics of two interacting CMEs observed on 13 and 14 June 2012. Both CMEs originated from the same active region NOAA 11504. After their launches which were separated by several hours, they were observed to interact at a distance of 100 Rs from the Sun. The interaction led to a moderate geomagnetic storm at the Earth with Dst index of approximately, -86 nT. The kinematics of the two CMEs is estimated using data from the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) onboard the Solar Terrestrial Relations Observatory (STEREO). Assuming a head-on collision scenario, we find that the collision is inelastic in nature. Further, the signatures of their interaction are examined using the in situ observations obtained by Wind and the Advance Composition Explorer (ACE) spacecraft. It is also found that this interaction event led to the strongest sudden storm commencement (SSC) (approximately 150 nT) of the present Solar Cycle 24. The SSC was of long duration, approximately 20 hours. The role of interacting CMEs in enhancing the geoeffectiveness is examined.

We present a comprehensive X-ray study of the neutron star low-mass X-ray binary IGR J17062-6143, which has been accreting at low luminosities since its discovery in $2006$. Analysing NuSTAR, XMM-Newton and Swift observations, we investigate the very faint nature of this source through three approaches: modelling the relativistic reflection spectrum to constrain the accretion geometry, performing high-resolution X-ray spectroscopy to search for an outflow, and searching for the recently reported millisecond X-ray pulsations. We find a strongly truncated accretion disk at $77^{+22}_{-18}$ gravitational radii ($\sim 164$ km) assuming a high inclination, although a low inclination and a disk extending to the neutron star cannot be excluded. The high-resolution spectroscopy reveals evidence for oxygen-rich circumbinary material, possibly resulting from a blueshifted, collisionally-ionised outflow. Finally, we do not detect any pulsations. We discuss these results in the broader context of possible explanations for the persistent faint nature of weakly accreting neutron stars. The results are consistent with both an ultra-compact binary orbit and a magnetically truncated accretion flow, although both cannot be unambigiously inferred. We also discuss the nature of the donor star and conclude that it is likely a CO or O-Ne-Mg white dwarf, consistent with recent multi-wavelength modelling.

The Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X) probes strong gravity for stellar mass to supermassive black holes and ultradense matter with unprecedented effective area, high time-resolution, and good spectral resolution, while providing a powerful time-domain X-ray observatory.

A millisecond pulsar is a neutron star that has been substantially spun up by accretion from a binary companion. A previously unrecognized factor governing the spin evolution of such pulsars is the crucial effect of non-steady or transient accretion. We numerically compute the evolution of accreting neutron stars through a series of outburst and quiescent phases considering the drastic variation of the accretion rate and the standard disk-magnetosphere interaction. We find that, for the same long-term average accretion rate, X-ray transients can spin up pulsars to rates several times higher than can persistent accretors, even when the spin down due to electromagnetic radiation during quiescence is included. We also compute an analytical expression for the equilibrium spin frequency in transients, by taking spin equilibrium to mean that no net angular momentum is transferred to the neutron star in each outburst cycle. We find that the equilibrium spin rate for transients, which depends on the peak accretion rate during outbursts, can be much higher than that for persistent sources. This explains our numerical finding. This finding implies that any meaningful study of neutron star spin and magnetic field distributions requires the inclusion of the transient accretion effect, since most accreting neutron star sources are transients. Our finding also implies the existence of a submillisecond pulsar population, which is not observed. This may point to the need for a competing spin-down mechanism for the fastest-rotating accreting pulsars, such as gravitational radiation.

Recent work on jets and disk winds in low-mass X-ray binaries (LMXBs) suggests that they are to a large extent mutually exclusive, with jets observed in spectrally hard states and disk winds observed in spectrally soft states. In this paper we use existing literature on jets and disk winds in the luminous neutron star (NS) LMXB GX 13+1, in combination with archival Rossi X-ray Timing Explorer data, to show that this source is likely able to produce jets and disk winds simultaneously. We find that jets and disk winds occur in the same location on the source's track in its X-ray color-color diagram. A further study of literature on other luminous LMXBs reveals that this behavior is more common, with indications for simultaneous jets and disk winds in the black hole LMXBs V404 Cyg and GRS 1915+105 and the NS LMXBs Sco X-1 and Cir X-1. For the three sources for which we have the necessary spectral information, we find that the simultaneous jets/winds all occur in their spectrally hardest states. Our findings indicate that in LMXBs with luminosities above a few tens of percent of the Eddington luminosity, jets and disk winds are not mutually exclusive, and that the presence of disk winds does not necessarily result in jet suppression.

We present a new test of hypothesis in which we seek the probability of the null conditioned on the data, where the null is a simplification undertaken to counter the intractability of the more complex model, that the simpler null model is nested within. With the more complex model rendered intractable, the null model uses a simplifying assumption that capacitates the learning of an unknown parameter vector given the data. Bayes factors are shown to be known only up to a ratio of unknown data-dependent constants--a problem that cannot be cured using prescriptions similar to those suggested to solve the problem caused to Bayes factor computation, by non-informative priors. Thus, a new test is needed in which we can circumvent Bayes factor computation. In this test, we undertake generation of data from the model in which the null hypothesis is true and can achieve support in the measured data for the null by comparing the marginalised posterior of the model parameter given the measured data, to that given such generated data. However, such a ratio of marginalised posteriors can confound interpretation of comparison of support in one measured data for a null, with that in another data set for a different null. Given an application in which such comparison is undertaken, we alternatively define support in a measured data set for a null by identifying the model parameters that are less consistent with the measured data than is minimally possible given the generated data, and realising that the higher the number of such parameter values, less is the support in the measured data for the null. Then, the probability of the null conditional on the data is given within an MCMC-based scheme, by marginalising the posterior given the measured data, over parameter values that are as, or more consistent with the measured data, than with the generated data.

Due to observational challenges our knowledge of low-level accretion flows around neutron stars is limited. We present Nustar, Swift and Chandra observations of the low-mass X-ray binary IGR J17062-6143, which has been persistently accreting at ~0.1 per cent of the Eddington limit since 2006. Our simultaneous Nustar/Swift observations show that the 0.5-79 keV spectrum can be described by a combination of a power law with a photon index of Gamma~2, a black body with a temperature of kT_bb~0.5 keV (presumably arising from the neutron star surface), and disk reflection. Modeling the reflection spectrum suggests that the inner accretion disk was located at R_in>100 GM/c2 (>225 km) from the neutron star. The apparent truncation may be due to evaporation of the inner disk into a radiatively-inefficient accretion flow, or due to the pressure of the neutron star magnetic field. Our Chandra gratings data reveal possible narrow emission lines near 1 keV that can be modeled as reflection or collisionally-ionized gas, and possible low-energy absorption features that could point to the presence of an outflow. We consider a scenario in which this neutron star has been able to sustain its low accretion rate through magnetic inhibition of the accretion flow, which gives some constraints on its magnetic field strength and spin period. In this configuration, IGR J17062-6143 could exhibit a strong radio jet as well as a (propeller-driven) wind-like outflow.

LOFT-P is a concept for a NASA Astrophysics Probe-Class (<$1B) X-ray timing mission, based on the LOFT concept originally proposed to ESAs M3 and M4 calls. LOFT-P requires very large collecting area (>6 m^2, >10x RXTE), high time resolution, good spectral resolution, broad-band spectral coverage (2-30 keV), highly flexible scheduling, and an ability to detect and respond promptly to time-critical targets of opportunity. It addresses science questions such as: What is the equation of state of ultra dense matter? What are the effects of strong gravity on matter spiraling into black holes? It would be optimized for sub-millisecond timing to study phenomena at the natural timescales of neutron star surfaces and black hole event horizons and to measure mass and spin of black holes. These measurements are synergistic to imaging and high-resolution spectroscopy instruments, addressing much smaller distance scales than are possible without very long baseline X-ray interferometry, and using complementary techniques to address the geometry and dynamics of emission regions. A sky monitor (2-50 keV) acts as a trigger for pointed observations, providing high duty cycle, high time resolution monitoring of the X-ray sky with ~20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger studies. A probe-class mission concept would employ lightweight collimator technology and large-area solid-state detectors, technologies which have been recently greatly advanced during the ESA M3 study. Given the large community interested in LOFT (>800 supporters, the scientific productivity of this mission is expected to be very high, similar to or greater than RXTE (~2000 refereed publications). We describe the results of a study, recently completed by the MSFC Advanced Concepts Office, that demonstrates that LOFT-P is feasible within a NASA probe-class mission budget.

A detailed investigation was carried out to understand a prolonged (~44 hours) weakly southward interplanetary magnetic field (IMF-Bz) condition during May 02--04, 1998. In-situ observations, during the period, showed the passage of an expanding magnetic cloud embedded in an interplanetary coronal mass ejection (ICME), followed up by a shock and an interplanetary discontinuity driven by another ICME. It is the arrival of the ICMEs and the upfront shocks that cause the prolonged southward IMF-Bz condition. The magnetic configuration of the source regions of the IMF associated with the ICME interval were also examined, which showed open magnetic field structures, emanating from a small active region on the north of the heliospheric current sheet (HCS). The structures remained constantly to the north of the HCS, both on April 29 and May 01, suggesting no change in their polarity. The draping of these outward directed radial field lines around the propagating CMEs in the shocked plasma explained the observed polarity changes of the IMF-Bz at 1 AU. In addition, multiple enhancements were also detected in the geomagnetic field variations, which showed a distinct one-to-one correspondence with the density pulses observed at 1 AU, during 0700--1700 UT on May 03. The spectral analysis of both the variations showed the same discrete frequencies of 0.48, 0.65 and 0.75 mHz, demonstrating that the solar wind density enhancements can cause detectable global geomagnetic disturbances. The observations, thus, provide a deeper insight into the possible causes and geomagnetic consequences of a prolonged weakly southward IMF-Bz condition.

We present results from multi-wavelength simultaneous X-ray and radio observations of the black hole X-ray binary V404 Cyg in quiescence. Our coverage with NuSTAR provides the very first opportunity to study the X-ray spectrum of V404 Cyg at energies above 10 keV. The unabsorbed broad-band (0.3--30 keV) quiescent luminosity of the source is 8.9$\times$10$^{32}$ erg s$^{-1}$ for a distance of 2.4 kpc. The source shows clear variability on short time scales (an hour to a couple of hours) in radio, soft X-ray and hard X-ray bands in the form of multiple flares. The broad-band X-ray spectra obtained from XMM-Newton and NuSTAR can be characterized with a power-law model having photon index $\Gamma$=2.12$\pm$0.07 (90\% confidence errors); however, residuals at high energies indicate spectral curvature significant at a 3$\sigma$ confidence level with e-folding energy of the cutoff to be 20$^{+20}_{-7}$ keV. Such curvature can be explained using synchrotron emission from the base of a jet outflow. Radio observations using the VLA reveal that the spectral index evolves on very fast time-scales (as short as 10 min.), switching between optically thick and thin synchrotron emission, possibly due to instabilities in the compact jet or stochastic instabilities in accretion rate. We explore different scenarios to explain this very fast variability.

We present NuSTAR and Swift observations of the neutron star Aquila X-1 during the peak of its July 2014 outburst. The spectrum is soft with strong evidence for a broad Fe K\alpha line. Modeled with a relativistically broadened reflection model, we find that the inner disk is truncated with an inner radius of 15+/-3 R_G. The disk is likely truncated by either the boundary layer and/or a magnetic field. Associating the truncated inner disk with pressure from a magnetic field gives an upper limit of B<5+/-2x10^8G. Although the radius is truncated far from the stellar surface, material is still reaching the neutron star surface as evidenced by the X-ray burst present in the t NuSTAR observation.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

Observational evidence has been accumulating that thermonuclear X-ray bursts ignited on the surface of neutron stars influence the surrounding accretion flow. Here, we exploit the excellent sensitivity of NuSTAR up to 79 keV to analyze the impact of an X-ray burst on the accretion emission of the neutron star LMXB 4U 1608-52. The ~200 s long X-ray burst occurred during a hard X-ray spectral state, and had a peak intensity of ~30-50 per cent of the Eddington limit with no signs of photospheric radius expansion. Spectral analysis suggests that the accretion emission was enhanced up to a factor of ~5 during the X-ray burst. We also applied a linear unsupervised decomposition method, namely non-negative matrix factorization (NMF), to study this X-ray burst. We find that the NMF performs well in characterizing the evolution of the burst emission and is a promising technique to study changes in the underlying accretion emission in more detail than is possible through conventional spectral fitting. For the burst of 4U 1608-52, the NMF suggests a possible softening of the accretion spectrum during the X-ray burst, which could potentially be ascribed to cooling of a corona. Finally, we report a small (~3 per cent) but significant rise in the accretion emission ~0.5 h before the X-ray burst, although it is unclear whether this was related to the X-ray burst ignition.

The X-ray binary 4U 1543$-$624 has been provisionally identified as an ultracompact system with an orbital period of $\simeq$18~min. We have carried out time-resolved optical imaging of the binary to verify the ultra-short orbital period. Using 140\,min of high-cadence $r'$-band photometry we recover the previously-seen sinusoidal modulation and determine a period $P=18.20\pm0.09$\,min. In addition, we also see a 7.0$\times 10^{-4}$\,mag\,min$^{-1}$ linear decay, likely related to variations in the source's accretion activity. Assuming that the sinusoidal modulation arises from X-ray heating of the inner face of the companion star, we estimate a distance of 6.0--6.7\,kpc and an inclination angle of 34$^{\circ}$--61$^{\circ}$ (90\% confidence) for the binary. Given the stability of the modulation we can confirm that the modulation is orbital in origin and 4U 1543$-$624 is an ultracompact X-ray binary.

We report on NuSTAR and Swift observations of a soft state of the neutron star low-mass X-ray binary GS 1826-24, commonly known as the "clocked" burster. The transition to the soft state was recorded in 2014 June through an increase of the 2-20 keV source intensity measured by MAXI, simultaneous with a decrease of the 15-50 keV intensity measured by Swift/BAT. The episode lasted approximately two months, after which the source returned to its usual hard state. We analyze the broad-band spectrum measured by Swift/XRT and NuSTAR, and estimate the accretion rate during the soft episode to be about 13% of Eddington, within the range of previous observations. However, the best fit spectral model, adopting the double Comptonization used previously, exhibits significantly softer components. We detect seven type-I X-ray bursts, all significantly weaker (and with shorter rise and decay times) than observed previously. The burst profiles and recurrence times vary significantly, ruling out the regular bursts that are typical for this source. One burst exhibited photospheric radius expansion, and we estimate the source distance at about (5.7 / xi_b^1/2) kpc, where xi_b parameterizes the possible anisotropy of the burst emission. Interpreting the soft state as a transition from an optically thin inner flow to an optically thick flow passing through a boundary layer, as is commonly observed in similar systems, is contradicted by the lower optical depth measured for the double Comptonization model we find for this soft state. The effect of a change in disk geometry on the burst behavior remains unclear.

We present spectral and timing analysis of NuSTAR observations of the accreting X-ray pulsar 2RXP J130159.6-635806. The source was serendipitously observed during a campaign focused on the gamma-ray binary PSR B1259-63 and was later targeted for a dedicated observation. The spectrum has a typical shape for accreting X-ray pulsars, consisting of a simple power law with an exponential cutoff starting at ~7 keV with a folding energy of E_fold=~18 keV. There is also an indication of the presence of a 6.4 keV iron line in the spectrum at the ~3 sigma significance level. NuSTAR measurements of the pulsation period reveal that the pulsar has undergone a strong and steady spin-up for the last 20 years. The pulsed fraction is estimated to be ~80%, and is constant with energy up to 40 keV. The power density spectrum shows a break towards higher frequencies relative to the current spin period. This, together with steady persistent luminosity, points to a long-term mass accretion rate high enough to bring the pulsar out of spin equilibrium.

Studying the reflection of X-rays off the inner edge of the accretion disk in a neutron star low-mass X-ray binary, allows us to investigate the accretion geometry and to constrain the radius of the neutron star. We report on a NuSTAR observation of 4U 1608-52 obtained during a faint outburst in 2014 when the neutron star, which has a known spin frequency of 620 Hz, was accreting at ~1-2% of the Eddington limit. The 3-79 keV continuum emission was dominated by a Gamma~2 power law, with a ~1-2% contribution from a kTbb~0.3-0.6 keV black body component. The high-quality NuSTAR spectrum reveals the hallmarks of disk reflection; a broad iron line peaking near 7~keV and a Compton back-scattering hump around ~20-30 keV. Modeling the disk reflection spectrum points to a binary inclination of i~30-40 degrees and a small `coronal' height of h<8.5 GM/c2. Furthermore, our spectral analysis suggests that the inner disk radius extended to Rin~7-10 GM/c2, close to the innermost stable circular obit. This constrains the neutron star radius to R<21 km and the redshift from the stellar surface to z>0.12, for a mass of M=1.5 Msun and a spin parameter of a=0.29.

We report on a 10 ks simultaneous Chandra/HETG-NuSTAR observation of the Bursting Pulsar, GRO J1744-28, during its third detected outburst since discovery and after nearly 18 years of quiescence. The source is detected up to 60 keV with an Eddington persistent flux level. Seven bursts, followed by dips, are seen with Chandra, three of which are also detected with NuSTAR. Timing analysis reveals a slight increase in the persistent emission pulsed fraction with energy (from 10% to 15%) up to 10 keV, above which it remains constant. The 0.5-70 keV spectra of the persistent and dip emission are the same within errors, and well described by a blackbody (BB), a power-law with an exponential rolloff, a 10 keV feature, and a 6.7 keV emission feature, all modified by neutral absorption. Assuming that the BB emission originates in an accretion disc, we estimate its inner (magnetospheric) radius to be about 4x10^7 cm, which translates to a surface dipole field B~9x10^10 G. The Chandra/HETG spectrum resolves the 6.7 keV feature into (quasi-)neutral and highly ionized Fe XXV and Fe XXVI emission lines. XSTAR modeling shows these lines to also emanate from a truncated accretion disk. The burst spectra, with a peak flux more than an order of magnitude higher than Eddington, are well fit with a power-law with an exponential rolloff and a 10~keV feature, with similar fit values compared to the persistent and dip spectra. The burst spectra lack a thermal component and any Fe features. Anisotropic (beamed) burst emission would explain both the lack of the BB and any Fe components.

As of 2014 August, the Nuclear Spectroscopic Telescope Array (NuSTAR) had observed ~30 X-ray binaries either as part of the planned program, as targets of opportunity, or for instrument calibration. The main science goals for the observations include probing the inner part of the accretion disk and constraining black hole spins via reflection components, providing the first observations of hard X-ray emission from quiescent Low Mass X-ray Binaries (LMXBs), measuring cyclotron lines from accreting pulsars, and studying type I X-ray bursts from neutron stars. Here, we describe the science objectives in more depth and give an overview of the NuSTAR observations that have been carried out to achieve the objectives. These include observation of four "IGR" High Mass X-ray Binaries (HMXBs) discovered by INTEGRAL. We also summarize the results that have been obtained and their implications. Among the IGR HMXBs, we focus on the discovery of a cyclotron line in the spectrum of IGR J17544-2619.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of thermonuclear X-ray bursts on accreting neutron stars. For a summary, we refer to the paper.

Tracking the spectral evolution of transiently accreting neutron stars between outburst and quiescence probes relatively poorly understood accretion regimes. Such studies are challenging because they require frequent monitoring of sources with luminosities below the thresholds of current all-sky X-ray monitors. We present the analysis of over 30 observations of the neutron star low-mass X-ray binary SAX J1750.8-2900 taken across four years with the X-ray telescope aboard Swift. We find spectral softening with decreasing luminosity both on long ($\sim$1 year) and short ($\sim$days to week) timescales. As the luminosity decreases from $4\times10^{36}$ erg s$^{-1}$ to $ \sim1\times10^{35} $ erg s$^{-1}$ (0.5-10 keV), the power law photon index increases from from 1.4 to 2.9. Although not statistically required, our spectral fits allow an additional soft component that displays a decreasing temperature as the luminosity decreases from $4 \times 10^{36} $ to $6 \times 10^{34}$ erg s$^{-1}$. Spectral softening exhibited by SAX J1750.8-2900 is consistent both with accretion emission whose spectral shape steepens with decreasing luminosity and also with being dominated by a changing soft component, possibly associated with accretion onto the neutron star surface, as the luminosity declines.

We present NuSTAR spectral and timing studies of the Supergiant Fast X-ray Transient (SFXT) IGR J17544-2619. The spectrum is well-described by a ~1 keV blackbody and a hard continuum component, as expected from an accreting X-ray pulsar. We detect a cyclotron line at 17 keV, confirming that the compact object in IGR J17544-2619 is indeed a neutron star. This is the first measurement of the magnetic field in a SFXT. The inferred magnetic field strength, B = (1.45 +/- 0.03) * 10^12 G * (1+z) is typical of neutron stars in X-ray binaries, and rules out a magnetar nature for the compact object. We do not find any significant pulsations in the source on time scales of 1-2000 s.

We report on two NuSTAR observations of GRS 1741.9-2853, a faint neutron star 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\times 10^{34}$ erg~s$^{-1}$ on 2013 July 31, and $5\times 10^{35}$ 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, a H mass fraction $X=0.7$ in the atmosphere, and a neutron star mass $M=1.4 M_{\odot}$, we determine a new lower limit on the distance for this source of $6.3 \pm 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 neutron star, 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 $\sigma$), a narrow absorption line at $5.46\pm0.10$ keV during the PRE phase of the burst, reminiscent of the detection by Waki et al. (1984). We propose that the line, if real, is formed in the wind above the photosphere of the neutron star by a resonant K$\alpha$ 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_{\odot}$.

The low-mass X-ray binary Cen X-4 is the brightest and closest (<1.2 kpc) quiescent neutron star transient. Previous 0.5-10 keV X-ray observations of Cen X-4 in quiescence identified two spectral components: soft thermal emission from the neutron star atmosphere and a hard power-law tail of unknown origin. We report here on a simultaneous observation of Cen X-4 with NuSTAR (3-79 keV) and XMM-Newton (0.3-10 keV) in 2013 January, providing the first sensitive hard X-ray spectrum of a quiescent neutron star transient. The 0.3-79 keV luminosity was 1.1 x 10^(33) erg/s (for D=1kpc), with around 60 percent in the thermal component. We clearly detect a cutoff of the hard spectral tail above 10 keV, the first time such a feature has been detected in this source class. We show that thermal Comptonization and synchrotron shock origins for the hard X-ray emission are ruled out on physical grounds. However, the hard X-ray spectrum is well fit by a thermal bremsstrahlung model with an 18 keV electron temperature, which can be understood as arising either in a hot layer above the neutron star atmosphere or in a radiatively-inefficient accretion flow (RIAF). The power-law cutoff energy may be set by the degree of Compton cooling of the bremsstrahlung electrons by thermal seed photons from the neutron star surface. Lower thermal luminosities should lead to higher (possibly undetectable) cutoff energies. We compare Cen~X-4's behavior with PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859, which have shown transitions between LMXB and radio pulsar modes at a similar X-ray luminosity.

Ultraluminous X-ray sources (ULX) are off-nuclear point sources in nearby galaxies whose X-ray luminosity exceeds the theoretical maximum for spherical infall (the Eddington limit) onto stellar-mass black holes. Their luminosity ranges from $10^{40}$ erg s$^{-1} < L_X$(0.5 - 10 keV) $<10^{40}$ erg s$^{-1}$. Since higher masses imply less extreme ratios of the luminosity to the isotropic Eddington limit theoretical models have focused on black hole rather than neutron star systems. The most challenging sources to explain are those at the luminous end ($L_X$ > $10^{40}$ erg s$^{-1}$), which require black hole masses MBH >50 solar masses and/or significant departures from the standard thin disk accretion that powers bright Galactic X-ray binaries. Here we report broadband X-ray observations of the nuclear region of the galaxy M82, which contains two bright ULXs. The observations reveal pulsations of average period 1.37 s with a 2.5-day sinusoidal modulation. The pulsations result from the rotation of a magnetized neutron star, and the modulation arises from its binary orbit. The pulsed flux alone corresponds to $L_X$(3 - 30 keV) = $4.9 \times 10^{39}$ erg s$^{-1}$. The pulsating source is spatially coincident with a variable ULX which can reach $L_X$ (0.3 - 10 keV) = $1.8 \times 10^{40}$ erg s$^{-1}$. This association implies a luminosity ~100 times the Eddington limit for a 1.4 solar mass object, or more than ten times brighter than any known accreting pulsar. This finding implies that neutron stars may not be rare in the ULX population, and it challenges physical models for the accretion of matter onto magnetized compact objects.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We propose a method to estimate the location of the Sun in the disk of the Milky Way using a method based on the Hellinger distance and construct confidence sets on our estimate of the unknown location using a bootstrap based method. Assuming the Galactic disk to be two-dimensional, the sought solar location then reduces to the radial distance separating the Sun from the Galactic center and the angular separation of the Galactic center to Sun line, from a pre-fixed line on the disk. On astronomical scales, the unknown solar location is equivalent to the location of us earthlings who observe the velocities of a sample of stars in the neighborhood of the Sun. This unknown location is estimated by undertaking pairwise comparisons of the estimated density of the observed set of velocities of the sampled stars, with densities estimated using synthetic stellar velocity data sets generated at chosen locations in the Milky Way disk according to four base astrophysical models. The "match" between the pair of estimated densities is parameterized by the affinity measure based on the familiar Hellinger distance. We perform a novel cross-validation procedure to establish a desirable "consistency" property of the proposed method.

Understanding of the kinematic evolution of Coronal Mass Ejections (CMEs) in the heliosphere is important to estimate their arrival time at the Earth. It is found that kinematics of CMEs can change when they interact or collide with each other as they propagate in the heliosphere. In this paper, we analyze the collision and post-interaction characteristics of two Earth-directed CMEs, launched successively on 2012 November 9 and 10, using white light imaging observations from STEREO/SECCHI and in situ observations taken from WIND spacecraft. We tracked two density enhanced features associated with leading and trailing edge of November 9 CME and one density enhanced feature associated with leading edge of November 10 CME by constructing J-maps. We found that the leading edge of November 10 CME interacted with the trailing edge of November 9 CME. We also estimated the kinematics of these features of the CMEs and found a significant change in their dynamics after interaction. In in situ observations, we identified distinct structures associated with interacted CMEs and also noticed their heating and compression as signatures of CME-CME interaction. Our analysis shows an improvement in arrival time prediction of CMEs using their post-collision dynamics than using pre-collision dynamics. Estimating the true masses and speeds of these colliding CMEs, we investigated the nature of observed collision which is found to be close to perfectly inelastic. The investigation also places in perspective the geomagnetic consequences of the two CMEs and their interaction in terms of occurrence of geomagnetic storm and triggering of magnetospheric substorms.

We report NuSTAR observations of the millisecond pulsar - low mass X-ray binary (LMXB) transition system PSR J1023+0038 from June and October 2013, before and after the formation of an accretion disk around the neutron star. Between June 10-12, a few days to two weeks before the radio disappearance of the pulsar, the 3-79 keV X-ray spectrum was well fit by a simple power law with a photon index of Gamma=1.17 +/-0.08 (at 90% confidence) with a 3-79 keV luminosity of 7.4+/-0.4 x 10^32 erg/s. Significant orbital modulation was observed with a modulation fraction of 36+/-10%. During the October 19-21 observation, the spectrum is described by a softer power law (Gamma=1.66+/-0.06) with an average luminosity of 5.8+/-0.2 x 10^33 erg/s and a peak luminosity of ~1.2 x 10^34 erg/s observed during a flare. No significant orbital modulation was detected. The spectral observations are consistent with previous and current multi-wavelength observations and show the hard X-ray power law extending to 79 keV without a spectral break. Sharp edged, flat bottomed `dips' are observed with widths between 30-1000 s and ingress and egress time-scales of 30-60 s. No change in hardness ratio was observed during the dips. Consecutive dip separations are log-normal in distribution with a typical separation of approximately 400 s. These dips are distinct from dipping activity observed in LMXBs. We compare and contrast these dips to observations of dips and state changes in the similar transition systems PSR J1824-2452I and XSS J1227.0-4859 and discuss possible interpretations based on the transitions in the inner disk.

Modeling the amplitudes and shapes of the X-ray pulsations observed from hot, rotating neutron stars provides a direct method for measuring neutron-star properties. This technique constitutes an important part of the science case for the forthcoming NICER and proposed LOFT X-ray missions. In this paper, we determine the number of distinct observables that can be derived from pulse profile modeling and show that using only bolometric pulse profiles is insufficient for breaking the degeneracy between inferred neutron-star radius and mass. However, we also show that for moderately spinning (300-800 Hz) neutron stars, analysis of pulse profiles in two different energy bands provides additional constraints that allow a unique determination of the neutron-star properties. Using the fractional amplitudes of the fundamental and the second harmonic of the pulse profile in addition to the amplitude and phase difference of the spectral color oscillations, we quantify the signal-to-noise ratio necessary to achieve a specified measurement precision for neutron star radius. We find that accumulating 10^6 counts in a pulse profile is sufficient to achieve a <5% uncertainty in the neutron star radius, which is the level of accuracy required to determine the equation of state of neutron-star matter. Finally, we formally derive the background limits that can be tolerated in the measurements of the various pulsation amplitudes as a function of the system parameters.

GRO J1008$-$57 is a high-mass X-ray binary for which several claims of a cyclotron resonance scattering feature near 80 keV have been reported. We use NuSTAR, Suzaku, and Swift data from its giant outburst of November 2012 to confirm the existence of the 80 keV feature and perform the most sensitive search to date for cyclotron scattering features at lower energies. We find evidence for a 78$^{+3}_{-2}$ keV line in the NuSTAR and Suzaku data at $> 4\sigma$ significance, confirming the detection using Suzaku alone by Yamamoto et al. (2014). A search of both the time-integrated and phase-resolved data rules out a fundamental at lower energies with optical depth larger than 5% of the 78 keV line. These results indicate that GRO J1008$-$57 has a magnetic field of $6.7\times10^{12} (1+z)$ G, the highest of known accreting pulsars.

We present a spectral analysis of three simultaneous NuSTAR and Swift/XRT observations of the transient Be-neutron star binary KS 1947+300 taken during its outburst in 2013/2014. These broad-band observations were supported by Swift/XRT monitoring snap-shots every 3 days, which we use to study the evolution of the spectrum over the outburst. We find strong changes of the power-law photon index, which shows a weak trend of softening with increasing X-ray flux. The neutron star shows very strong pulsations with a period of P ~ 18.8 s. The 0.8-79 keV broad-band spectrum can be described by a power-law with an exponential cutoff and a black-body component at low energies. During the second observation we detect a cyclotron resonant scattering feature at 12.5 keV, which is absent in the phase-averaged spectra of observations 1 and 3. Pulse phase-resolved spectroscopy reveals that the strength of the feature changes strongly with pulse phase and is most prominent during the broad minimum of the pulse profile. At the same phases the line also becomes visible in the first and third observation at the same energy. This discovery implies that KS 1947+300 has a magnetic field strength of B ~ 1.1e12 (1+z)G, which is at the lower end of known cyclotron line sources.

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four magnetars, two rotation-powered pulsars and the cataclysmic variable AE Aquarii. NuSTAR also discovered 3.76-s pulsations from the transient source SGR J1745-29 recently found by Swift very close to the Galactic Center, clearly identifying the source as a transient magnetar. For magnetar 1E 1841-045, we show that the spectrum is well fit by an absorbed blackbody plus broken power-law model with a hard power-law photon index of ~1.3. This is consistent with previous results by INTEGRAL and RXTE. We also find an interesting double-peaked pulse profile in the 25-35 keV band. For AE Aquarii, we show that the spectrum can be described by a multi-temperature thermal model or a thermal plus non-thermal model; a multi-temperature thermal model without a non-thermal component cannot be ruled out. Furthermore, we do not see a spiky pulse profile in the hard X-ray band, as previously reported based on Suzaku observations. For other magnetars and rotation-powered pulsars observed with NuSTAR, data analysis results will be soon available.

In this paper we focus on a type of inverse problem in which the data is expressed as an unknown function of the sought and unknown model function (or its discretised representation as a model parameter vector). In particular, we deal with situations in which training data is not available. Then we cannot model the unknown functional relationship between data and the unknown model function (or parameter vector) with a Gaussian Process of appropriate dimensionality. A Bayesian method based on state space modelling is advanced instead. Within this framework, the likelihood is expressed in terms of the probability density function ($pdf$) of the state space variable and the sought model parameter vector is embedded within the domain of this $pdf$. As the measurable vector lives only inside an identified sub-volume of the system state space, the $pdf$ of the state space variable is projected onto the space of the measurables, and it is in terms of the projected state space density that the likelihood is written; the final form of the likelihood is achieved after convolution with the distribution of measurement errors. Application motivated vague priors are invoked and the posterior probability density of the model parameter vectors, given the data is computed. Inference is performed by taking posterior samples with adaptive MCMC. The method is illustrated on synthetic as well as real galactic data.

We report on an observation of the neutron star low-mass X-ray binary Serpens X-1, made with NuSTAR. The extraordinary sensitivity afforded by NuSTAR facilitated the detection of a clear, robust, relativistic Fe K emission line from the inner disk. A relativistic profile is required over a single Gaussian line from any charge state of Fe at the 5-sigma level of confidence, and any two Gaussians of equal width at the same confidence. The Compton back-scattering "hump" peaking in the 10-20 keV band is detected for the first time in a neutron star X-ray binary. Fits with relativistically-blurred disk reflection models suggest that the disk likely extends close to the innermost stable circular orbit (ISCO) or stellar surface. The best-fit blurred reflection models constrain the gravitational redshift from the stellar surface to be z > 0.16. The data are broadly compatible with the disk extending to the ISCO; in that case, z > 0.22 and R < 12.6 km (assuming M = 1.4 Msun and a=0, where a = cJ/GM^2). If the star is as large or larger than its ISCO, or if the effective reflecting disk leaks across the ISCO to the surface, the redshift constraints become measurements. We discuss our results in the context of efforts to measure fundamental properties of neutron stars, and models for accretion onto compact objects.

We report on the results of spectral fits made to a NuSTAR observation of the black hole GRS 1915+105 in a "plateau" state. This state is of special interest because it is similar to the "low/hard" state seen in other black holes, especially in that compact, steady jets are launched in this phase. The 3-79 keV bandpass of NuSTAR, and its ability to obtain moderate-resolution spectra free from distortions such as photon pile-up, are extremely well suited to studies of disk reflection in X-ray binaries. In only 15 ks of net exposure, an extraordinarily sensitive spectrum of GRS 1915+105 was measured across the full bandpass. Ionized reflection from a disk around a rapidly-spinning black hole is clearly required to fit the spectra; even hybrid Comptonization models including ionized reflection from a disk around a Schwarzschild black hole proved inadequate. A spin parameter of a = 0.98 +/- 0.01 (1-sigma statistical error) is measured via the best-fit model; low spins are ruled out at a high level of confidence. This result suggests that jets can be launched from a disk extending to the innermost stable circular orbit. A very steep inner disk emissivity profile is also measured, consistent with models of compact coronae above Kerr black holes. These results support an emerging association between the hard X-ray corona and the base of the relativistic jet.

The NuSTAR hard X-ray telescope observed the transient Be/X-ray binary GS 0834-430 during its 2012 outburst - the first active state of this system observed in the past 19 years. We performed timing and spectral analysis, and measured the X-ray spectrum between 3-79keV with high statistical significance. We find the phase-averaged spectrum to be consistent with that observed in many other magnetized accreting pulsars. We fail to detect cyclotron resonance scattering features that would allow us to constrain the pulsar's magnetic field in either phase-averaged or phase-resolved spectra. Timing analysis shows a clearly detected pulse period of ~12.29s in all energy bands. The pulse profiles show a strong, energy-dependent hard phase lag of up to 0.3 cycles in phase, or about 4s. Such dramatic energy-dependent lags in the pulse profile have never before been reported in high-mass X-ray binary (HMXB) pulsars. Previously reported lags have been significantly smaller in phase and restricted to low-energies (E<10keV). We investigate the possible mechanisms that might produce this energy-dependent pulse phase shift. We find the most likely explanation for this effect to be a complex beam geometry.

Using both wavelet and Fourier analysis, a study has been undertaken of the changes in the quasi-periodic variations in solar photospheric fields in the build-up to one of the deepest solar minima experienced in the past 100 years. This unusual and deep solar minimum occurred between solar cycles 23 and 24. The study, carried out using ground based synoptic magnetograms spanning the period 1975.14 to 2009.86, covered solar cycles 21, 22 and 23. A hemispheric asymmetry in periodicities of the photospheric fields was seen only at latitudes above $\pm45{^{\circ}}$ when the data was divided, based on a wavelet analysis, into two parts: one prior to 1996 and the other after 1996. Furthermore, the hemispheric asymmetry was observed to be confined to the latitude range 45${^{\circ}}$ to 60${^{\circ}}$. This can be attributed to the variations in polar surges that primarily depend on both the emergence of surface magnetic flux and varying solar surface flows. The observed asymmetry when coupled with the fact that both solar fields above $\pm45{^{\circ}}$ and micro-turbulence levels in the inner-heliosphere have been decreasing since the early to mid nineties (Janardhan et al., 2011) suggests that around this time active changes occurred in the solar dynamo that governs the underlying basic processes in the sun. These changes in turn probably initiated the build-up to the very deep solar minimum at the end of the cycle 23. The decline in fields above $\pm45{^{\circ}}$ for well over a solar cycle, would imply that weak polar fields have been generated in the past two successive solar cycles \textit{viz.} cycles 22 and 23. A continuation of this declining trend beyond 22 years, if it occurs, will have serious implications on our current understanding of the solar dynamo.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

We report on Chandra X-ray and Gemini-North near-infrared K-band observations of the eclipsing accretion-powered millisecond X-ray pulsar SWIFT J1749.4-2807 in quiescence. Using the Chandra observation we derive a source position of Right Ascencion: 17:49:31.73 and Declination:-28:08:05.08. The position is accurate to 0.6" (90 per cent confidence). We find one source at a magnitude K=18.44+-0.03 with a position fully consistent with the accurate Chandra X-ray localization and a second source at K=19.2+-0.1 that falls close to the edge of the error circle in the deep K-band images. The presence of a few weaker sources as suggested by previous H-band observations presented in the literature cannot the ruled out. There is marginal evidence that the brighter of the these two sources is variable. Follow-up spectroscopy of this potential counterpart will show if this source is the true counterpart to SWIFT J1749.4-2807. If so, baring the presence of complicating effects such as heating of the mass-donor star, it would allow for the mass of the neutron star to be measured through the measurement of periodic radial velocity variations.

Thermonuclear bursts from slowly accreting neutron stars (NSs) have proven difficult to detect, yet they are potential probes of the thermal properties of the neutron star interior. During the first year of a systematic all-sky search for X-ray bursts using the Gamma-ray Burst Monitor (GBM) aboard the Fermi Gamma-ray Space Telescope we have detected 15 thermonuclear bursts from the NS low-mass X-ray binary 4U 0614+09, when it was accreting at nearly 1% of the Eddington limit. We measured an average burst recurrence time of 12+/-3 d (68% confidence interval) between March 2010 and March 2011, classified all bursts as normal duration bursts and placed a lower limit on the recurrence time of long/intermediate bursts of 62 d (95% confidence level). We discuss how observations of thermonuclear bursts in the hard X-ray band compare to pointed soft X-ray observations, and quantify such bandpass effects on measurements of burst radiated energy and duration. We put our results for 4U 0614+09 in the context of other bursters and briefly discuss the constraints on ignition models. Interestingly, we find that the burst energies in 4U 0614+09 are on average between those of normal duration bursts and those measured in long/intermediate bursts. Such a continuous distribution in burst energy provides a new observational link between normal and long/intermediate bursts. We suggest that the apparent bimodal distribution that defined normal and long/intermediate duration bursts during the last decade could be due to an observational bias towards detecting only the longest and most energetic bursts from slowly accreting NSs.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We present a comprehensive study of the thermonuclear bursts and millihertz quasi-periodic oscillations (mHz QPOs) from the neutron star (NS) transient and 11 Hz X-ray pulsar IGR J17480-2446, located in the globular cluster Terzan 5. The increase in burst rate that we found during its 2010 outburst, when persistent luminosity rose from 0.1 to 0.5 times the Eddington limit, is in qualitative agreement with thermonuclear burning theory yet opposite to all previous observations of thermonuclear bursts. Thermonuclear bursts gradually evolved into a mHz QPO when the accretion rate increased, and vice versa. The mHz QPOs from IGR J17480-2446 resemble those previously observed in other accreting NSs, yet they feature lower frequencies (by a factor ~3) and occur when the persistent luminosity is higher (by a factor 4-25). We find four distinct bursting regimes and a steep (close to inverse cubic) decrease of the burst recurrence time with increasing persistent luminosity. We compare these findings to nuclear burning models and find evidence for a transition between the pure helium and mixed hydrogen/helium ignition regimes when the persistent luminosity was about 0.3 times the Eddington limit. We also point out important discrepancies between the observed bursts and theory, which predicts brighter and less frequent bursts, and suggest that an additional source of heat in the NS envelope is required to reconcile the observed and expected burst properties. We discuss the impact of NS magnetic field and spin on the expected nuclear burning regimes, in the context of this particular pulsar.

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

The accreting millisecond pulsar SAX J1808.4-3658 has shown a peculiar orbital evolution in the past with an orbital expansion much faster than expected from standard binary evolutionary scenarios. Previous limits on the pulsar spin frequency derivative during transient accretion outbursts were smaller than predicted by standard magnetic accretion torque theory, while the spin evolution between outbursts was consistent with magnetic dipole spin-down. In this paper we present the results of a coherent timing analysis of the 2011 outburst observed by the Rossi X-ray Timing Explorer and extend our previous long-term measurements of the orbital and spin evolution over a baseline of thirteen years. We find that the expansion of the 2 hr orbit is accelerating at a rate 1.6E-20 s/s^2 and we interpret this as the effect of short-term angular momentum exchange between the mass donor and the orbit. The gravitational quadrupole coupling due to variations in the oblateness of the companion can be a viable mechanism for explaining the observations. No significant spin frequency derivatives are detected during the 2011 outburst (<4E-13 Hz/s) and the long term spin down remains stable over thirteen years with a rate of approximately -1E-15 Hz/s.

A non-proprietary Mission Concept available for presentation to NASA, providing High Area, High Resolution Imaging Spectroscopy and Timing with Arcmin Angular Resolution Submitted in response to NASA 2011 RFI NNH11ZDA018L 'Concepts for the Next NASA X-ray Astronomy Mission'

AXTAR is a NASA MIDEX mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument is a collimated, thick Si pixel detector with 2-50 keV coverage and over 3 square meters effective area. For timing observations of accreting neutron stars and black holes, AXTAR provides at least a factor of five improvement in sensitivity over the RXTE PCA. AXTAR also carries a sensitive sky monitor that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and design choices that face a next generation timing mission. We then describe the technical concept for AXTAR and summarize a preliminary mission design study at the NASA/MSFC Advanced Concepts Office.

The neutron star transient and 11 Hz X-ray pulsar IGR J17480-2446, recently discovered in the globular cluster Terzan 5, showed unprecedented bursting activity during its 2010 October-November outburst. We analyzed all X-ray bursts detected with the Rossi X-ray Timing Explorer and find strong evidence that they all have a thermonuclear origin, despite the fact that many do not show the canonical spectral softening along the decay imprinted on type I X-ray bursts by the cooling of the neutron star photosphere. We show that the persistent-to-burst power ratio is fully consistent with the accretion-to-thermonuclear efficiency ratio along the whole outburst, as is typical for type I X-ray bursts. The burst energy, peak luminosity and daily-averaged spectral profiles all evolve smoothly throughout the outburst, in parallel with the persistent (non-burst) luminosity. We also find that the peak burst to persistent luminosity ratio determines whether or not cooling is present in the bursts from IGR J17480-2446, and argue that the apparent lack of cooling is due to the "non-cooling" bursts having both a lower peak temperature and a higher non-burst (persistent) emission. We conclude that the detection of cooling along the decay is a sufficient, but not a necessary condition to identify an X-ray burst as thermonuclear. Finally, we compare these findings with X-ray bursts from other rapidly accreting neutron stars.

The accretion-powered millisecond pulsar IGR J00291+5934 underwent two ~10 d long outbursts during 2008, separated by 30 d in quiescence. Such a short quiescent period between outbursts has never been seen before from a neutron star X-ray transient. X-ray pulsations at the 599 Hz spin frequency are detected throughout both outbursts. For the first time, we derive a pulse phase model that connects two outbursts, providing a long baseline for spin frequency measurement. Comparison with the frequency measured during the 2004 outburst of this source gives a spin-down during quiescence of -4(1)x10^-15 Hz/s, approximately an order of magnitude larger than the long-term spin-down observed in the 401 Hz accretion-powered pulsar SAX J1808.4-3658. If this spin-down is due to magnetic dipole radiation, it requires a 2x10^8 G field strength, and its high spin-down luminosity may be detectable with the Fermi Large Area Telescope. Alternatively, this large spin-down could be produced by gravitational wave emission from a fractional mass quadrupole moment of Q/I = 1x10^{-9}. The rapid succession of the outbursts also provides a unique test of models for accretion in low-mass X-ray binaries. Disk instability models generally predict that an outburst will leave the accretion disk too depleted to fuel a second outburst after such a brief quiescence. We suggest a modification in which the outburst is shut off by the onset of a propeller effect before the disk is depleted. This model can explain the short quiescence and the unusually slow rise of the light curve of the second 2008 outburst.

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

X-ray observations of EXO 0748-676 during thermonuclear bursts revealed a set of narrow absorption lines that potentially originate from the stellar photosphere. The identification of these lines with particular atomic transitions led to the measurement of the surface gravitational redshift of the neutron star and to constraints on its mass and radius. However, the recent detection of 552 Hz oscillations at 15% rms amplitude revealed the spin frequency of the neutron star and brought into question the consistency of such a rapid spin with the narrow width of the absorption lines. Here, we calculate the amplitudes of burst oscillations and the width of absorption lines emerging from the surface of a rapidly rotating neutron star for a wide range of model parameters. We show that no combination of neutron-star and geometric parameters can simultaneously reproduce the narrowness of the absorption lines, the high amplitude of the oscillations, and the observed flux at the time the oscillations were detected. We, therefore, conclude that the observed absorption lines are unlikely to originate from the surface of this neutron star.

We present the results of a systematic search in approximately 14 years of Rossi X-ray Timing Explorer All-Sky Monitor data for evidence of periodicities not reported by Wen et al. (2006). Two variations of the commonly used Fourier analysis search method have been employed to achieve significant improvements in sensitivity. The use of these methods and the accumulation of additional data have resulted in the detection of the signatures of the orbital periods of eight low-mass X-ray binary systems and of ten high-mass X-ray binaries not listed in the tables of Wen et al.

We report the detection of 15 X-ray bursts with RXTE and Swift observations of the peculiar X-ray binary Circinus X-1 during its May 2010 X-ray re-brightening. These are the first X-ray bursts observed from the source after the initial discovery by Tennant and collaborators, twenty-five years ago. By studying their spectral evolution, we firmly identify nine of the bursts as type I (thermonuclear) X-ray bursts. We obtain an arcsecond location of the bursts that confirms once and for all the identification of Cir X-1 as a type I X-ray burst source, and therefore as a low magnetic field accreting neutron star. The first five bursts observed by RXTE are weak and show approximately symmetric light curves, without detectable signs of cooling along the burst decay. We discuss their possible nature. Finally, we explore a scenario to explain why Cir X-1 shows thermonuclear bursts now but not in the past, when it was extensively observed and accreting at a similar rate.

The Advanced X-ray Timing Array (AXTAR) is a mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument, the Large Area Timing Array (LATA) is a collimated instrument with 2-50 keV coverage and over 3 square meters effective area. The LATA is made up of an array of supermodules that house 2-mm thick silicon pixel detectors. AXTAR will provide a significant improvement in effective area (a factor of 7 at 4 keV and a factor of 36 at 30 keV) over the RXTE PCA. AXTAR will also carry a sensitive Sky Monitor (SM) that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and technical concept for AXTAR and present results from a preliminary mission design study.

We observed the ultra-compact binary candidate 4U 0614+091 for a total of 200 ksec with the high-energy transmission gratings onboard the \chandra X-ray Observatory. The source is found at various intensity levels with spectral variations present. X-ray luminosities vary between 2.0$\times10^{36}$ \ergsec and 3.5$\times10^{36}$ \ergsec. Continuum variations are present at all times and spectra can be well fit with a powerlaw component, a high kT blackbody component, and a broad line component near oxygen. The spectra require adjustments to the Ne K edge and in some occasions also to the Mg K edge. The Ne K edge appears variable in terms of optical depths and morphology. The edge reveals average blue- and red-shifted values implying Doppler velocities of the order of 3500 \kms. The data show that Ne K exhibits excess column densities of up to several 10$^{18}$ cm$^{-2}$. The variability proves that the excess is intrinsic to the source. The correponding disk velocities also imply an outer disk radius of the order of $< 10^9$ cm consistent with an ultra-compact binary nature. We also detect a prominent soft emission line complex near the \oviii L$\alpha$ position which appears extremely broad and relativistic effects from near the innermost disk have to be included. Gravitationally broadened line fits also provide nearly edge-on angles of inclination between 86 and 89$^{\circ}$. The emissions appear consistent with an ionized disk with ionization parameters of the order of 10$^4$ at radii of a few 10$^7$ cm. The line wavelengths with respect to \oviiia\ are found variably blue-shifted indicating more complex inner disk dynamics.

We report the detection of pulsations at 552 Hz in the rising phase of two type-I (thermonuclear) X-ray bursts observed from the accreting neutron star EXO 0748-676 in 2007 January and December, by the Rossi X-ray Timing Explorer. The fractional amplitude was 15% (rms). The dynamic power density spectrum for each burst revealed an increase in frequency of approx. 1-2 Hz while the oscillation was present. The frequency drift, the high significance of the detections and the almost identical signal frequencies measured in two bursts separated by 11 months, confirms this signal as a burst oscillation similar to those found in 13 other sources to date. We thus conclude that the spin frequency in EXO 0748-676 is within a few Hz of 552 Hz, rather than 45 Hz as was suggested from an earlier signal detection by Villarreal & Strohmayer (2004). Consequently, Doppler broadening must significantly affect spectral features arising from the neutron star surface, so that the narrow absorption features previously reported from an XMM-Newton spectrum could not have arisen there. The origin of both the previously reported 45 Hz oscillation and the X-ray absorption lines is now uncertain.

The X-ray source 4U1820-30 in the globular cluster NGC 6624 is known as the most compact binary among the identified X-ray binaries. Having an orbital period of 685.0 s, the source consists of a neutron star primary and likely 0.06--0.08 Msun white dwarf secondary. Here we report on far-ultraviolet (FUV) observations of this X-ray binary, made with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. From our Fourier spectral analysis of the FUV timing data, we obtain a period of 693.5+/-1.3 s, which is significantly different from the orbital period. The light curve folded at this period can be described by a sinusoid, with a fractional semiamplitude of 6.3% and the phase zero (maximum of the sinusoid) at MJD 50886.015384+/-0.000043 (TDB). While the discovered FUV period may be consistent with a hierarchical triple system model that was previously considered for 4U 1820-30, we suggest that it could instead be the indication of superhump modulation, which arises from an eccentric accretion disk in the binary. The X-ray and FUV periods would be the orbital and superhump periods, respectively, indicating a 1% superhump excess and a white-dwarf/neutron-star mass ratio around 0.06. Considering 4U 1820-30 as a superhump source, we discuss the implications.

We present a new Bayesian non-parametric deprojection algorithm DOPING (Deprojection of Observed Photometry using and INverse Gambit), that is designed to extract 3-D luminosity density distributions $\rho$ from observed surface brightness maps $I$, in generalised geometries, while taking into account changes in intrinsic shape with radius, using a penalised likelihood approach and an MCMC optimiser. We provide the most likely solution to the integral equation that represents deprojection of the measured $I$ to $\rho$. In order to keep the solution modular, we choose to express $\rho$ as a function of the line-of-sight (LOS) coordinate $z$. We calculate the extent of the system along the ${\bf z}$-axis, for a given point on the image that lies within an identified isophotal annulus. The extent along the LOS is binned and density is held a constant over each such $z$-bin. The code begins with a seed density and at the beginning of an iterative step, the trial $\rho$ is updated. Comparison of the projection of the current choice of $\rho$ and the observed $I$ defines the likelihood function (which is supplemented by Laplacian regularisation), the maximal region of which is sought by the optimiser (Metropolis Hastings). The algorithm is successfully tested on a set of test galaxies, the morphology of which ranges from an elliptical galaxy with varying eccentricity to an infinitesimally thin disk galaxy marked by an abruptly varying eccentricity profile. Applications are made to faint dwarf elliptical galaxy Ic~3019 and another dwarf elliptical that is characterised by a central spheroidal nuclear component superimposed upon a more extended flattened component. The result of deprojection of the X-ray image of triaxial cluster A1413 is also presented.

The low mass X-ray binary 4U 2129+47 was discovered during a previous X-ray outburst phase and was classified as an accretion disk corona source. A 1% delay between two mid-eclipse epochs measured ~22 days apart was reported from two XMM-Newton observations taken in 2005, providing support to the previous suggestion that 4U 2129+47 might be in a hierarchical triple system. In this work we present timing and spectral analysis of three recent XMM-Newton observations of 4U 2129+47, carried out between November 2007 and January 2008. We found that absent the two 2005 XMM-Newton observations, all other observations are consistent with a linear ephemeris with a constant period of 18857.63s; however, we confirm the time delay reported for the two 2005 XMM-Newton observations. Compared to a Chandra observation taken in 2000, these new observations also confirm the disappearance of the sinusoidal modulation of the lightcurve as reported from two 2005 XMM-Newton observations. We further show that, compared to the Chandra observation, all of the XMM-Newton observations have 40% lower 0.5--2 keV absorbed fluxes, and the most recent XMM-Newton observations have a combined 2--6 keV flux that is nearly 80% lower. Taken as a whole, the timing results support the hypothesis that the system is in a hierarchical triple system (with a third body period of at least 175 days). The spectral results raise the question of whether the drop in soft X-ray flux is solely attributable to the loss of the hard X-ray tail (which might be related to the loss of sinusoidal orbital modulation), or is indicative of further cooling of the quiescent neutron star after cessation of residual, low-level accretion.

We present the results of an X-ray mass analysis of the early-type galaxy NGC 4636, using Chandra data. We have compared the X-ray mass density profile with that derived from a dynamical analysis of the system's globular clusters (GCs). Given the observed interaction between the central active galactic nucleus and the X-ray emitting gas in NGC 4636, we would expect to see a discrepancy in the masses recovered by the two methods. Such a discrepancy exists within the central ~10kpc, which we interpret as the result of non-thermal pressure support or a local inflow. However, over the radial range ~10-30kpc, the mass profiles agree within the 1-sigma errors, indicating that even in this highly disturbed system, agreement can be sought at an acceptable level of significance over intermediate radii, with both methods also indicating the need for a dark matter halo. However, at radii larger than 30kpc, the X-ray mass exceeds the dynamical mass, by a factor of 4-5 at the largest disagreement. A Fully Bayesian Significance Test finds no statistical reason to reject our assumption of velocity isotropy, and an analysis of X-ray mass profiles in different directions from the galaxy centre suggests that local disturbances at large radius are not the cause of the discrepancy. We instead attribute the discrepancy to the paucity of GC kinematics at large radius, coupled with not knowing the overall state of the gas at the radius where we are reaching the group regime (>30kpc), or a combination of the two.

X-ray emission from the enigmatic Rotating RAdio Transients (RRATs) offers a vital clue to understanding these objects and how they relate to the greater neutron star population. An X-ray counterpart to J1819-1458 is known, and its properties are similar to those of other middle-aged (0.1 Myr) neutron stars. We have searched for X-ray emission with Chandra/ACIS at the positions of two RRATs with arcsecond (or better) localisation, J0847-4316 and J1846-0257. Despite deep searches (especially for J1847-0257) we did not detect any emission with 0.3-8 keV count-rate limits of 1 counts/ks and 0.068 counts/ks, respectively, at 3sigma confidence. Assuming thermal emission similar to that seen from J1819-1458 (a blackbody with radius of approximately 20 km), we derive effective temperature limits of 77 eV and 91 eV for the nominal values of the distances and column densities to both sources, although both of those quantities are highly uncertain and correlated. If we instead fix the temperature of the emission (a blackbody with kT=0.14 keV), we derive unabsorbed luminosity limits in the 0.3-8 keV range of 1e32 erg/s and 3e32 erg/s. These limits are considerably below the luminosity of J1819-1458 (4e33 erg/s), suggesting that RRATs J0847-4316 and J1846-0257 have cooled beyond the point of visibility (plausible given the differences in characteristic age). However, as we have not detected X-ray emission, it may also be that the emission from RRATs J0847-4316 and J1846-0257 has a different character from that of J1819-1458. The two non-detections may prove a counterpoint to J1819-1458, but more detections are certainly needed before we can begin to derive general X-ray emission properties for the RRAT populations.

The Rossi X-ray Timing Explorer has observed five outbursts from the transient 2.5 ms accretion-powered pulsar SAX J1808.4-3658 during 1998-2008. We present a pulse timing study of the most recent outburst and compare it with the previous timing solutions. The spin frequency of the source continues to decrease at a rate of (-5.5+/-1.2)x10^-18 Hz/s, which is consistent with the previously determined spin derivative. The spin-down occurs mostly during quiescence, and it is most likely due to the magnetic dipole torque from a B = 1.5x10^8 G dipolar field at the neutron star surface. We also find that the 2 hr binary orbital period is increasing at a rate of (3.80+/-0.06)x10^-12 s/s, also consistent with previous measurements. It remains uncertain whether this orbital change reflects secular evolution or short-term variability.

The dark matter halos in $\Lambda$CDM cosmological simulations are triaxial and highly flattened. In many cases, these triaxial equilibria are also tumbling slowly, typically about their short axes, with periods of order a Hubble time. Halos may therefore exert a slowly changing external torque on spiral galaxies that can affect their dynamical evolution in interesting ways. We examine the effect of the external torques exerted by a tumbling quadrupolar tidal field on the evolution of spiral galaxies using N-body simulations with realistic, disk galaxy models. We measure the amplitude of the external quadrupole moments of dark halos in cosmological simulations and use these to force disk galaxy models in a series of N-body experiments for a range of pattern speeds. We find that the torques are strong enough to induce long lived transient warps in disks similar to those observed in real spirals and also induce the bar instability at later times in some galaxy models that are otherwise stable for long periods of time in isolation. We also observe forced spiral structure near the edge of the disk where normally self-gravity is too weak to be responsible for such structure. This overlooked influence of dark halos may well be responsible for many of the peculiar aspects of disk galaxy dynamics.

We report the discovery of a 4.5 um counterpart to the anomalous X-ray pulsar (magnetar) 1E 2259+586 with the Spitzer Space Telescope. The mid-infrared flux density is 6.3+/-1.0 uJy at 4.5 um and <20 uJy (at 95% confidence) at 8 um, or 0.02% of the 2-10 keV X-ray flux (corrected for extinction). Combining our Spitzer measurements with previously published near-infrared data, we show that the overall infrared emission from 1E 2259+586 is qualitatively similar to that from the magnetar 4U 0142+61. Therefore the passive X-ray-heated dust disk model originally developed for 4U 0142+61 might also apply to 1E 2259+586. However, the IR data from this source can also be fitted by a simple power-law spectrum as might be expected from magnetospheric emission.

The state of a non-relativistic gravitational dynamical system is known at any time $t$ if the dynamical rule, i.e. Newton's equations of motion, can be solved; this requires specification of the gravitational potential. The evolution of a bunch of phase space coordinates ${\bf w}$ is deterministic, though generally non-linear. We discuss the novel Bayesian non-parametric algorithm CHASSIS that gives phase space $pdf$ $f({\bf w})$ and potential $\Phi({\bf x})$ of a relaxed gravitational system. CHASSIS is undemanding in terms of input requirements in that it is viable given incomplete, single-component velocity information of system members. Here ${\bf x}$ is the 3-D spatial coordinate and ${\bf w}={\bf x+v}$ where ${\bf v}$ is the 3-D velocity vector. CHASSIS works with a 2-integral $f=f(E, L)$ where energy $E=\Phi + v^2/2, \: v^2 = \sum_{i=1}^{3}{v_i^2}$ and the angular momentum is $L = |{\bf r}\times{\bf v}|$, where ${\bf r}$ is the spherical spatial vector. Also, we assume spherical symmetry. CHASSIS obtains the $f(\cdot)$ from which the kinematic data is most likely to have been drawn, in the best choice for $\Phi(\cdot)$, using an MCMC optimiser (Metropolis-Hastings). The likelihood function ${\cal{L}}$ is defined in terms of the projections of $f(\cdot)$ into the space of observables and the maximum in ${\cal{L}}$ is sought by the optimiser.

Soft phase lags, in which X-ray pulses in lower energy bands arrive later than pulses in higher energy bands, have been observed in nearly all accretion-powered millisecond pulsars, but their origin remains an open question. In a study of the 2.5 ms accretion-powered pulsar SAX J1808.4-3658, we report that the magnitude of these lags is strongly dependent on the accretion rate. During the brightest stage of the outbursts from this source, the lags increase in magnitude as the accretion rate drops; when the outbursts enter their dimmer flaring-tail stage, the relationship reverses. We evaluate this complex dependence in the context of two theoretical models for the lags, one relying on the scattering of photons by the accretion disk and the other invoking a two-component model for the photon emission. In both cases, the turnover suggests that we are observing the source transitioning into the "propeller" accretion regime.

We present a coherent timing analysis of the 2003 outburst of the accreting millisecond pulsar XTE J1807-294. We find an upper limit for the spin frequency derivative of 5E-14 Hz/s. The sinusoidal fractional amplitudes of the pulsations are the highest observed among the accreting millisecond pulsars and can reach values of up to 27% (2.5-30 keV). The pulse arrival time residuals of the fundamental follow a linear anti-correlation with the fractional amplitudes that suggests hot spot motion over the surface of the neutron star both in longitude and latitude. An anti-correlation between residuals and X-ray flux suggests an influence of accretion rate on pulse phase, and casts doubts on the use of standard timing techniques to measure spin frequencies and torques on the neutron star.

We develop a novel way of finding total mass density profiles in Sersic ellipticals, to about 3 times the major axis effective radius, using no other information other than what is typically available for distant galaxies, namely the observed surface brightness distribution and the central velocity dispersion $\sigma_0$. The luminosity density profile of the observed galaxy is extracted by deprojecting the measured brightness distribution and scaling it by a fiduciary, step-function shaped, $raw$ mass-to-light ratio profile ($M/L$). The resulting raw, discontinuous, total, 3-D mass density profile is then smoothed according to a proposed smoothing prescription. The parameters of this raw $M/L$ are characterised by implementing the observables in a model-based study. The complete characterisation of the formalism is provided as a function of the measurements of the brightness distribution and $\sigma_0$. The formalism, thus specified, is demonstrated to yield the mass density profiles of a suite of test galaxies and is successfully applied to extract the gravitational mass distribution in NGC 3379 and NGC 4499, out to about 3 effective radii.

We present ongoing Rossi X-ray Timing Explorer (RXTE) monitoring observations of the 377.3 Hz accretion-powered pulsar, HETE J1900.1-2455 Activity continues in this system more than 3 years after discovery, at a mean luminosity of 4.4e36 erg/s (for d=5 kpc), although pulsations were present only within the first 70 days. X-ray variability has increased each year, notably with a brief interval of nondetection in 2007, during which the luminosity dropped to below 1e-3 of the mean level. A deep search of data from the intervals of nondetection in 2005 revealed evidence for extremely weak pulsations at an amplitude of 0.29% rms, a factor of ten less than the largest amplitude seen early in the outburst. X-ray burst activity continued through 2008, with bursts typically featuring strong radius expansion. Spectral analysis of the most intense burst detected by RXTE early in the outburst revealed unusual variations in the inferred photospheric radius, as well as significant deviations from a blackbody. We obtained much better fits instead with a comptonisation model.

AXTAR is an X-ray observatory mission concept, currently under study in the U.S., that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument is a collimated, thick Si pixel detector with 2-50 keV coverage and 8 square meters collecting area. For timing observations of accreting neutron stars and black holes, AXTAR provides at least an order of magnitude improvement in sensitivity over both RXTE and Constellation-X. AXTAR also carries a sensitive sky monitor that acts as a trigger for pointed observations of X-ray transients and also provides continuous monitoring of the X-ray sky with 20 times the sensitivity of the RXTE ASM. AXTAR builds on detector and electronics technology previously developed for other applications and thus combines high technical readiness and well understood cost.

The spin frequency distribution of accreting millisecond X-ray pulsars cuts off sharply above 730 Hz, well below the breakup spin rate for most neutron star equations of state. I review several different ideas for explaining this cutoff. There is currently considerable interest in the idea that gravitational radiation from rapidly rotating pulsars might act to limit spin up by accretion, possibly allowing eventual direct detection with gravitational wave interferometers. I describe how long-term X-ray timing of fast accreting millisecond pulsars like the 599 Hz source IGR J00291+5934 can test the gravitational wave model for the spin frequency limit.

The radio sky is poorly sampled for rapidly varying transients because of the narrow field-of-view of most imaging radio telescopes at cm and shorter wavelengths. The emergence of sensitive long wavelength observations with intrinsically larger fields-of-view are changing this situation, as partly illustrated by our ongoing meter-wavelength monitoring observations and archival studies of the Galactic Center. In this search, we discovered a transient, bursting, radio source in the direction of the Galactic Center, GCRT J1745-3009, with extremely unusual properties. Its flux and rapid variability imply a brightness temperature >10^12 K if it is at a distance >70 pc, implying that it is a coherent emitter. I will discuss the discovery of the source and the subsequent re-detections, as well as searches for counterparts at other wavelengths, and several proposed models.

We present an optical/near-infrared search for a counterpart to the perplexing radio transient GCRT J1745-3009, a source located ~1 degree from the Galactic Center. Motivated by some similarities to radio bursts from nearby ultracool dwarfs, and by a distance upper limit of 70 pc for the emission to not violate the 1e12 K brightness temperature limit for incoherent radiation, we searched for a nearby star at the position of GCRT J1745-3009. We found only a single marginal candidate, limiting the presence of any late-type star to >1 kpc (spectral types earlier than M9), >200 pc (spectral types L and T0-T4), and >100 pc (spectral types T4-T7), thus severely restricting the possible local counterparts to GCRT J1745-3009. We also exclude any white dwarf within 1 kpc or a supergiant star out to the distance of the Galactic Center as possible counterparts. This implies that GCRT J1745-3009 likely requires a coherent emission process, although whether or not it reflects a new class of sources is unclear.

We have assembled a sample of 1187 thermonuclear (type-I) X-ray bursts from 48 accreting neutron stars by the Rossi X-ray Timing Explorer, spanning more than ten years. The sample contains examples of two of the three theoretical ignition regimes and likely examples of the third. We present a detailed analysis of the variation of the burst profiles, energetics, recurrence times, presence of photospheric radius expansion, and presence of burst oscillations, as a function of accretion rate. We estimated the distance for 35 sources exhibiting radius-expansion bursts, and found that the peak flux of such bursts varies typically by 13%. We classified sources into two main groups based on the burst properties: both long and short bursts (indicating mixed H/He accretion), and consistently short bursts (primarily He accretion). The decrease in burst rate observed for both groups at >0.06 Mdot_Edd (>~2E37 erg/s) is associated with a transition in the persistent spectral state and (as has been suggested previously) may be related to the increasing role of steady He-burning. We found examples of bursts separated by <30 min, including burst triplets and even quadruplets. We describe the oscillation amplitudes for 13 of the 16 burst oscillation sources, as well as the stages and properties of the bursts in which the oscillations are detected. The burst properties are correlated with the burst oscillation frequency; sources at <400 Hz generally have consistently short bursts, while the more rapidly-spinning systems have both long and short bursts. This correlation suggests either that shear-mediated mixing dominates the burst properties, or that the nature of the mass donor (and hence the evolutionary history) has an influence on the long-term spin evolution.

{The determination of cluster masses is a complex problem that would be aided by information about the cluster shape and orientation (along the line-of-sight).} {It is in this context, that we have developed a scheme for identifying the intrinsic morphology and inclination of a cluster, by looking for the signature of the true cluster characteristics in the inter-comparison of the different deprojected emissivity profiles (that all project to the same X-ray brightness distribution) and by using SZe data when available.} {We deproject the cluster X-ray surface brightness profile under the assumptions of four different geometry and inclination configurations for the observed system; these 4 configurations correspond to four extreme geometry+inclination scenarios. The deprojection in question is performed by the non-parametric algorithm DOPING. The formalism is tested with model systems and then is applied to a sample of 24 clusters. While the shape determination is possible by implementing the X-ray brightness alone, the estimation of the inclination is usually markedly improved by the usage of SZe data that is available for the considered sample.}{We spot 8 prolate systems, 1 oblate and 15 of the clusters in our sample as triaxial. In fact, for systems identified as triaxial, we are able to discern how the three semi-axes lengths compare with each other. This, when compounded by the information about the line-of-sight extent, allows us to constrain the inclination quite tightly and offers the two intrinsic axial ratios of the triaxial systems.}{}

{} {To quantify the amount of chaos that exists in the local phase space.} {A sample of orbits from four different models of the Solar neighbourhood phase space are analysed by a new chaos identification (and quantification) technique. While three of the used models bear the signature of the perturbation due to both the Galactic bar and the spiral pattern, the last of the models is a bar only one. We explore the models by inter-comparing the corresponding values of chaos strength that is induced at the various energy levels .}{(1) We find that of all the viable models that have been demonstrated to successfully reproduce the local phase space structure, i.e. those that include the bar as well as the spiral, bear strong chaoticity, though the model that implies the highest degree of chaos is the one in which overlap of the major resonances of the bar and the spiral occurs. The bar only model is found to display regularity. (2) We advance chaos to be primarily responsible for the splitting of the Hyades-Pleiades mode (the larger mode) of the local velocity distribution}{}

We present the density structure of dark matter in the outer parts (to about 7 effective radii) of the elliptical galaxy NGC 4636, from the radial velocities of 174 globular clusters, using the non-parametric, Bayesian algorithm CHASSIS. We find the galaxy to be rich in dark matter, with $R$-band mass-to-light ($M/L$) ratios rising to about 30, at nearly 4$R_e$; the $K$-band $M/L$ at about 3$R_e$ is found to be nearly 10. The result does not depend on applying the method to the red and blue globular clusters separately. This estimate of $M/L$ is higher than the previous analysis from the same kinematic data. We also find that the dark matter distribution is highly concentrated towards the inner halo.

SWIFT J1756.9-2508 is one of the few accreting millisecond pulsars (AMPs) discovered to date. We report here the results of our analysis of its aperiodic X-ray variability, as measured with the Rossi X-ray Timing Explorer during the 2007 outburst of the source. We detect strong (~35%) flat-topped broadband noise throughout the outburst with low characteristic frequencies (~0.1 Hz). This makes SWIFT J1756.9-2508 similar to the rest of AMPs and to other low luminosity accreting neutron stars when they are in their hard states, and enables us to classify this AMP as an atoll source in the extreme island state. We also find a hard tail in its energy spectrum extending up to 100 keV, fully consistent with such source and state classification.

We attempt to understand the state of the local phase space by comparing simulated 2-D velocity distributions to the distribution that is constructed for the solar neighbourhood, from measurements of stellar radial and transverse velocities. The joint perurbation of the central bar in the Galaxy and the spiral pattern is found to be a must, in order to produce successful models of the local phase space. The existence of chaos is found to be an important ingredient in the formation of the observed phase space structure.

We present a new non-parametric deprojection algorithm DOPING (Deprojection of Observed Photometry using and INverse Gambit), that is designed to extract the three dimensional luminosity density distribution $\rho$, from the observed surface brightness profile of an astrophysical system such as a galaxy or a galaxy cluster, in a generalised geometry, while taking into account changes in the intrinsic shape of the system. The observable is the 2-D surface brightness distribution of the system. While the deprojection schemes presented hitherto have always worked within the limits of an assumed intrinsic geometry, in DOPING, geometry and inclination can be provided as inputs. The $\rho$ that is most likely to project to the observed brightness data is sought; the maximisation of the likelihood is performed with the Metropolis algorithm. Unless the likelihood function is maximised, $\rho$ is tweaked in shape and amplitude, while maintaining positivity, but otherwise the luminosity distribution is allowed to be completely free-form. Tests and applications of the algorithm are discussed.

We present a 7 yr timing study of the 2.5 ms X-ray pulsar SAX J1808.4-3658, an X-ray transient with a recurrence time of ~2 yr, using data from the Rossi X-ray Timing Explorer covering 4 transient outbursts (1998-2005). We verify that the 401 Hz pulsation traces the spin frequency fundamental and not a harmonic. Substantial pulse shape variability, both stochastic and systematic, was observed during each outburst. Analysis of the systematic pulse shape changes suggests that, as an outburst dims, the X-ray "hot spot" on the pulsar surface drifts longitudinally and a second hot spot may appear. The overall pulse shape variability limits the ability to measure spin frequency evolution within a given X-ray outburst (and calls previous nudot measurements of this source into question), with typical upper limits of |nudot| < 2.5x10^{-14} Hz/s (2 sigma). However, combining data from all the outbursts shows with high (6 sigma) significance that the pulsar is undergoing long-term spin down at a rate nudot = (-5.6+/-2.0)x10^{-16} Hz/s, with most of the spin evolution occurring during X-ray quiescence. We discuss the possible contributions of magnetic propeller torques, magnetic dipole radiation, and gravitational radiation to the measured spin down, setting an upper limit of B < 1.5x10^8 G for the pulsar's surface dipole magnetic field and and Q/I < 5x10^{-9} for the fractional mass quadrupole moment. We also measured an orbital period derivative of Pdot = (3.5+/-0.2)x10^{-12} s/s. This surprising large Pdot is reminiscent of the large and quasi-cyclic orbital period variation observed in the so-called "black widow" millisecond radio pulsars, supporting speculation that SAX J1808.4-3658 may turn on as a radio pulsar during quiescence. In an appendix we derive an improved (0.15 arcsec) source position from optical data.

We report on the discovery by the Swift Gamma-Ray Burst Explorer of the eighth known transient accretion-powered millisecond pulsar, SWIFT J1756.9-2508, as part of routine observations with the Swift Burst Alert Telescope hard X-ray transient monitor. The pulsar was subsequently observed by both the X-Ray Telescope on Swift and the Rossi X-Ray Timing Explorer Proportional Counter Array. It has a spin frequency of 182 Hz (5.5 ms) and an orbital period of 54.7 minutes. The minimum companion mass is between 0.0067 and 0.0086 solar masses, depending on the mass of the neutron star, and the upper limit on the mass is 0.030 solar masses (95% confidence level). Such a low mass is inconsistent with brown dwarf models, and comparison with white dwarf models suggests that the companion is a He-dominated donor whose thermal cooling has been at least modestly slowed by irradiation from the accretion flux. No X-ray bursts, dips, eclipses or quasi-periodic oscillations were detected. The current outburst lasted approximately 13 days and no earlier outbursts were found in archival data.

We report on our analysis of 300 ks of Chandra observations of the neutron star soft X-ray transient 1H1905+000 in quiescence. We do not detect the source down to a 95% confidence unabsorbed flux upper limit of 2E-16 erg cm-2 s-1 in the 0.5-10 keV energy range for an assumed Gamma=2 power law spectral model. A limit of 1.4E-16 erg cm-2 s-1 is derived if we assume that the spectrum of 1H1905+000 in quiescence is described well with a black body of temperature of 0.2 keV. For the upper limit to the source distance of 10 kpc this yields a 0.5-10 keV luminosity limit of 2.4E30 / 1.7E30 erg/s for the abovementioned power law or black body spectrum, respectively. This luminosity limit is lower than the luminosity of A0620-00, the weakest black hole soft X-ray transient in quiescence reported so far. Together with the uncertainties in relating the mass transfer and mass accretion rates we come to the conclusion that the claim that there is evidence for the presence of a black hole event horizon on the basis of a lower quiescent luminosity for black holes than for neutron stars is unproven. We also briefly discuss the implications of the low quiescent luminosity of 1H1905+000 for the neutron star equation of state. Using deep Magellan images of the field of 1H1905+000 obtained at excellent observing conditions we do not detect the quiescent counterpart at the position of the outburst optical counterpart down to a magnitude limit of i'>25.3. This can be converted to a limit on the absolute magnitude of the counterpart of I>9.6 which implies that the counterpart can only be a brown or a white dwarf.

We examine the idea that dynamical parameters can be estimated by identifying locations in the solar neighbourhood where simulated velocity distributions match the observed local distribution. Here, the dynamical influence of both the Galactic bar and the outer spiral pattern are taken into account. The Milky Way disc is stirred by analytical potentials that are chosen to represent the two perturbations, the ratio of pattern speeds of which is explored, rather than held constant. The velocity structure of the final configuration is presented as heliocentric velocity distributions at different locations. These model velocity distributions are compared to the observed distribution in terms of a goodness-of-fit parameter that has been formulated here. We monitor the spatial distribution of the maximal value of this parameter, in order to constrain the solar position from a model. Efficiency of a model is based on a study of this distribution as well as on other independent dynamical considerations. We reject the bar only and spiral only models and arrive at the following bar parameters from the bar+spiral simulations: bar pattern speed of 57.4^{+2.8}_{-3.3} km/s/kpc and a bar angle in [0^\circ$, 30^\circ$]. However, extracting information in this way is no longer viable when the dynamical influence of the spiral pattern does not succumb to that of the bar. Orbital analysis indicates that even though the basic bimodality in the local velocity distribution can be attributed to scattering off the Outer Lindblad Resonance of the bar, it is the interaction of irregular orbits and orbits of other resonant families, that is responsible for the other moving groups; it is realised that such interaction increases with the warmth of the background disk.

Using archival X-ray data, we find that the catalog location of the X-ray binary Scutum X-1 (Sct X-1) is incorrect, and that the correct location is that of the X-ray source AX J183528-0737, which is 15' to the west. Our identification is made on the basis of the 112-s pulse period for this object detected in an XMM-Newton observation, as well as spatial coincidence between AX J183528-0737 and previous X-ray observations. Based on the XMM-Newton data and archival RXTE data, we confirm secular spin-down over 17 years with period derivative Pdot~3.9e-9 s/s, but do not detect a previously reported X-ray iron fluorescence line. We identify a bright (Ks=6.55), red (J-Ks=5.51), optical and infrared counterpart to AX J183528-0737 from 2MASS, a number of mid-IR surveys, and deep optical observations, which we use to constrain the extinction to and distance of Sct X-1. From these data, as well as limited near-IR spectroscopy, we conclude that Sct X-1 is most likely a binary system comprised of a late-type giant or supergiant and a neutron star.

Mass modelling of early-type systems is a thorny issue. In this paper, we present a new mass modelling formalism for ellipticals that invokes no observations other than the central velocity dispersion ($\sigma_0$) and photometry. The essence of the method lies in choosing a local mass-to-light ratio ($M/L$) profile for a galaxy, with which the deprojected luminosity density distribution (along the major axis coordinate $x$) is scaled. The resulting discontinuous mass density profile is then smoothed, according to a laid out prescription; the local $M/L$ profile that stems from this smoothed mass density, is found to be significantly different from the raw $M/L$ distribution. A suite of model galaxies (both Sersic and cored in nature) is used for extensive experimentation in order to characterise this raw $M/L$ profile and in each case, the mass density recovered from this mass modelling technique is compared to the known mass distribution. We opt to work with a raw $M/L$ profile that is a simple two-stepped function of $x$, with a low inner and higher outer value of $M/L-\Upsilon_{in}$ and $\Upsilon_{out}$, respectively. The only constraint that we have on this profile is in the centre of the galaxy, via $\sigma_0$. This value of $\sigma_0$ is implemented in the virial theorem to obtain an estimate of the central $M/L$ ratio of the galaxy. The fallibility of the virial mass estimate is taken care of, by allowing for a range in the values of $\Upsilon_{in}$ that can be used for a given galaxy model. Moreover, our experiments indicate that $\Upsilon_{out}$ is uniquely known, for a given $\Upsilon_{in}$; for cored galaxies, this functional form is found uniquely dependant on the core radius. The jump radius of the raw $M/L$ profile is chosen to be thrice the effective radius of the galaxy.

[abridged] We report results from four recent observations of the ultracompact LMXB pulsar 4U 1626-67. All the observations obtained high-resolution X-ray spectra of the system, two from the Chandra X-ray Observatory using the HETGS, and two from the XMM-Newton Observatory using the RGS as well as the EPIC PN and MOS. These data allow us to study in detail the prominent Ne and O emission line complexes which make 4U 1626-67 unique among LMXBs. The observations were spaced over a period of 3 years for a total observing time of 238 ks, allowing us to monitor the line regions as well as the overall source flux, continuum spectrum, and timing properties. The structure of the emission lines and the ratios of the components of the helium-like Ne IX and O VII triplets support the hypothesis that they are formed in the high-density environment of the accretion disk. We do not find any significant changes in the line widths or ratios over this time period, though we note that the line equivalent widths decrease. We are able to place constraints on the strengths of the Ne K, Fe L, and O K photoelectric absorption edges, and find that the data do not require an overabundance of Ne or O in the system relative to the expected ISM values. We find that the pulsar is still spinning down, and note that the pulse profile has changed significantly from what was found prior to the torque reversal in 1990, suggesting that this event may be linked to a change in the geometry of the accretion column. The flux of 4U 1626-67 continues to decrease, in keeping with the trend of the last approximately 30 years over which it has been observed. Taking into consideration current theory on disk stability, we expect that 4U 1626-67 will enter a period of quiescence in 2-15 years.

We describe observations of the seventh accretion-powered millisecond pulsar, HETE J1900.1-2455 made with the Rossi X-ray Timing Explorer during the year of activity that followed its discovery in 2005 June. We detected intermittent pulsations at a peak fractional amplitude of 3%, but only in the first two months of the outburst. On three occasions during this time we observed an abrupt increase in the pulse amplitude, approximately coincident with the time of a thermonuclear burst, followed by a steady decrease on a timescale of approx. 10 d. HETE J1900.1-2455 has shown the longest active period by far for any transient accretion-powered millisecond pulsar, comparable instead to the outburst cycles for other transient X-ray binaries. Since the last detection of pulsations, HETE J1900.1-2455 has been indistinguishable from a low-accretion rate, non-pulsing LMXB; we hypothesize that other, presently active LMXBs may have also been detectable initially as millisecond X-ray pulsars.

We report on our search for the optical/infrared counterparts to the central compact objects in four young supernova remnants: Pup A, PKS 1209-52, RCW 103, and Cas A. The X-ray point sources in these supernova remnants are excellent targets for probing the existence of supernova fallback disks, since irradiation of a disk by a central X-ray source should lead to an infrared excess. We used ground-based optical and near-infrared imaging and Spitzer Space Telescope mid-infrared imaging to search for optical/infrared counterparts at the X-ray point source positions measured by the Chandra X-Ray Observatory. We did not detect any counterparts, and hence find no evidence for fallback disks around any of these sources. In PKS 1209-52, we are able to exclude a nearby optical/infrared candidate counterpart. In RCW 103, a blend of 3 faint stars at the X-ray source position prevents us from deriving useful limits. For the other targets, the upper limits on the infrared/X-ray flux ratio are as deep as (1.0--1.7)$\times 10^{-4}$. Comparing these limits to the ratio of $\approx 6\times10^{-5}$ measured for 4U 0142+61 (a young pulsar recently found with an X-ray irradiated dust disk), we conclude that the non-detection of any disks around young neutron stars studied here are consistent with their relatively low X-ray luminosities, although we note that a similar dust disk around the neutron star in Pup A should be detectable by deeper infrared observations.

We present high-resolution spectroscopy of the neon K-shell and iron L-shell interstellar absorption edges in nine X-ray binaries using the High Energy Transmission Grating Spectrometer (HETGS) onboard the Chandra X-ray Observatory. We found that the iron absorption is well fit by an experimental determination of the cross-section for metallic iron, although with a slight wavelength shift of ~20 mA. The neon edge region is best fit by a model that includes the neutral neon edge and three Gaussian absorption lines. We identify these lines as due to the 1s-2p transitions from Ne II, Ne III, and Ne IX. As we found in our oxygen edge study, the theoretical predictions for neutral and low-ionization lines all require shifts of ~20 mA to match our data. Combined with our earlier oxygen edge study, we find that a best fit O/Ne ratio of 5.4+/-1.6, consistent with standard interstellar abundances. Our best fit Fe/Ne ratio of 0.20+/-0.03 is significantly lower than the interstellar value. We attribute this difference to iron depletion into dust grains in the interstellar medium. We make the first measurement of the neon ionization fraction in the ISM. We find Ne II/Ne I ~ 0.3 and Ne III/Ne I ~ 0.07. These values are larger than is expected given the measured ionization of interstellar helium. For Ne IX, our results confirm the detection of the hot ionized interstellar medium of the Galaxy.

Pulsars are rotating, magnetized neutron stars that are born in supernova explosions following the collapse of the cores of massive stars. If some of the explosion ejecta fails to escape, it may fall back onto the neutron star or it may possess sufficient angular momentum to form a disk. Such 'fallback' is both a general prediction of current supernova models and, if the material pushes the neutron star over its stability limit, a possible mode of black hole formation. Fallback disks could dramatically affect the early evolution of pulsars, yet there are few observational constraints on whether significant fallback occurs or even the actual existence of such disks. Here we report the discovery of mid-infrared emission from a cool disk around an isolated young X-ray pulsar. The disk does not power the pulsar's X-ray emission but is passively illuminated by these X-rays. The estimated mass of the disk is of order 10 Earth masses, and its lifetime (at least a million years) significantly exceeds the spin-down age of the pulsar, supporting a supernova fallback origin. The disk resembles protoplanetary disks seen around ordinary young stars, suggesting the possibility of planet formation around young neutron stars.

We present the high-resolution Chandra X-ray Observatory persistent (non-dip) spectrum of 4U 1916-05 which revealed narrow absorption lines from hydrogenic neon, magnesium, silicon, and sulfur, in addition to the previous identified hydrogenic and helium-like iron absorption lines. This makes 4U 1916-05 only the second of the classical X-ray dipper systems to show narrow absorption lines from elements other than iron. We propose two possible explanations for the small measured line widths (>~ 500-2000 km s^{-1}), compared to the expected Keplerian velocities (> 1000 km s^{-1}) of the accretion disk in this 50-min orbital period system, and lack of wavelength shifts (>~ 250 km s^{-1}). First, the ionized absorber may be stationary. Alternatively, the line properties may measure the relative size of the emission region. From this hypothesis, we find that the emission region is constrained to be >~ 0.25 times the radial extent of the absorber. Our results also imply that the ionized absorber spans a range of ionization parameters.

In this Paper we report on our analysis of a ~25 ksec. Chandra X-ray observation of the neutron star soft X-ray transient (SXT) 1H1905+000 in quiescence. Furthermore, we discuss our findings of the analysis of optical photometric observations which we obtained using the Magellan telescope and photometric and spectroscopic observations which we obtained using the Very Large Telescope at Paranal. The X-ray counterpart of 1H1905+000 was not detected in our Chandra data, with a 95 per cent confidence limit to the source count rate of 1.2x10^-4 counts s^-1. For different spectral models this yields an upper limit on the luminosity of 1.8x10^31 erg s^-1 (for an upper limit on the distance of 10 kpc.) This luminosity limit makes 1H1905+000 the faintest neutron star SXT in quiescence observed to date. The neutron star luminosity is so low that it is similar to the lowest luminosities derived for black hole SXTs in quiescence. This low luminosity for a neutron star SXT challanges the hypothesis presented in the literature that black hole SXTs in quiescence have lower luminosities than neutron star SXTs as a result of the presence of a black hole event horizon. Furthermore, the limit on the neutron star luminosity obtained less than 20 years after the outburst has ceased, constrains the thermal conductivity of the neutron star crust. Finally, the neutron star core must be so cold that unless the time averaged mass accretion rate is lower than 2x10^-12 M_sun yr^-1, core cooling has to proceed via enhanced neutrino emission processes. We derive a limit on the absolute I-band magnitude of the quiescent counterpart of M_I>7.8 assuming the source is at 10 kpc. This is in line with 1H1905+000 being an ultra-compact X-ray binary, as has been proposed based on the low outburst V-band absolute magnitude.

We report the detection of probable optical counterparts for two Intermediate Mass Binary Pulsar (IMBP) systems, PSR J1528-3146 and PSR J1757-5322. Recent radio pulsar surveys have uncovered a handful of these systems with putative massive white dwarf companions, thought to have an evolutionary history different from that of the more numerous class of Low Mass Binary Pulsars (LMBPs) with He white dwarf companions. The study of IMBP companions via optical observations offers us several new diagnostics: the evolution of main sequence stars near the white-dwarf-neutron star boundary, the physics of white dwarfs close to the Chandrasekhar limit, and insights into the recycling process by which old pulsars are spun up to high rotation frequencies. We were unsuccessful in our attempt to detect optical counterparts of PSR J1141-6545, PSR J1157-5112, PSR J1435-6100, and PSR J1454-5846.

The observed radial and transverse velocities of individual stars in M15 are implemented as inputs to a fully non-parametric code (CHASSIS) in order to estimate the equilibrium stellar distribution function and the three-dimensional mass density profile. In particular, the paper explores the possibility of the existence of a central black hole in M15 via several runs that utilize the radial velocity data set which offers kinematic measurements closer to the centre of the cluster than the proper motion data. These runs are distinguished from each other in the choice of the initial seed for the cluster characteristics; however, the profiles identified by the algorithm at the end of each run concur with each other, within error bars, thus confirming the robustness of CHASSIS. The recovered density profiles are noted to exhibit unequivocal flattening, inner to about 0.0525pc. Also, the enclosed mass profile is very close to being a power-law function of radius inside 0.1pc and is not horizontal. Simplistically speaking, these trends negate the possibility of the central mass to be concentrated in a black hole, the lower bound on the radius of the sphere of influence of which would be $\gtrsim$0.041pc, had it existed. However, proper analysis suggests that the mass enclosed within the inner 0.01pc could be in the form of a black hole of mass $\sim{10^3}$M$_{\odot}$, under two different scenarios, which are discussed. The line-of-sight velocity dispersion is visually found to be very similar to the observed dispersion profile. The enclosed mass and velocity dispersion profiles calculated from runs done with the proper motion data are found to be consistent with the profiles obtained with the radial velocity data.

We analyzed 123 thermonuclear (type-I) X-ray bursts observed by the Rossi X-ray Timing Explorer from the low-mass X-ray binary 4U 1636-536. All but two of the 40 radius-exansion bursts in this sample reached peak fluxes which were normally distributed about a mean of 6.4e-8 ergs/cm^2/s, with a standard deviation of 7.6%. The remaining two radius-expansion bursts reached peak fluxes a factor of 1.69+/-0.13 lower than this mean value; as a consequence, the overall variation in the peak flux of the radius-expansion bursts was a factor of ~2. This variation is comparable to the range of the Eddington limit between material with solar H-fraction (X=0.7) and pure He. Such a variation may arise if, for the bright radius-expansion bursts, most of the accreted H is eliminated either by steady hot CNO burning or expelled in a radiatively-driven wind. However, steady burning cannot exhaust the accreted H for solar composition material within the typical ~2 hr burst recurrence time, nor can it result in sufficient elemental stratification to allow selective ejection of the H only. An additional stratification mechanism appears to be required to separate the accreted elements and thus allow preferential ejection of the hydrogen. We found no evidence for a gap in the peak flux distribution between the radius-expansion and non-radius expansion bursts, previously observed in smaller samples. Assuming that the faint radius-expansion bursts reached the Eddington limit for H-rich material (X~0.7), and the brighter bursts the limit for pure He (X=0), we estimate the distance to 4U 1636-536 (for a canonical neutron star with M_NS=1.4M_sun, R_NS=10 km) to be 6.0+/-0.5 kpc, or for M_NS=2M_sun at most 7.1 kpc. (Abstract abridged)

In this Paper we report on our analysis of three Chandra observations of the accretion-powered millisecond X-ray pulsar IGR J00291+5934 obtained during the late stages of the 2004 outburst. We also report the serendipitous detection of the source in quiescence by ROSAT during MJD 48830-48839. The detected 0.3-10 keV source count rates varied significantly between the Chandra observations from (7.2+-1.2)x10^-3, (6.8+-0.9)x10^-3, and (1.4+-0.1)x10^-2 counts per second for the 1st, 2nd, and 3rd Chandra observation, on MJD 53371.88, 53383.99, and 53407.57, respectively. The count rate for the 3rd observation is 2.0+-0.4 times as high as that of the average of the first two observations. The unabsorbed 0.5-10 keV source flux for the best-fit power-law model to the source spectrum was (7.9+-2.5)x10^-14, (7.3+-2.0)x10^-14, and (1.17+-0.22)x10^-13 erg cm^-2 s^-1 for the 1st, 2nd, and 3rd Chandra observation, respectively. We find that this source flux is consistent with that found by ROSAT [~(5.4+-2.4)x10^-14 erg cm^-2 s^-1]. Under the assumption that the interstellar extinction, N_H, does not vary between the observations, we find that the blackbody temperature during the 2nd Chandra observation is significantly higher than that during the 1st and 3rd observation. Furthermore, the effective temperature of the neutron star derived from fitting an absorbed blackbody or neutron star atmosphere model to the data is rather high in comparison with many other neutron star soft X-ray transients in quiescence, even during the 1st and 3rd observation. If we assume that the source quiescent luminosity is similar to that measured for two other accretion powered millisecond pulsars in quiescence, the distance to IGR J00291+5934 is 2.6-3.6 kpc.

We report the precise optical and X-ray localization of the 3.2 ms accretion-powered X-ray pulsar XTE J1814-338 with data from the Chandra X-Ray Observatory as well as optical observations conducted during the 2003 June discovery outburst. Optical imaging of the field during the outburst of this soft X-ray transient reveals an R = 18 star at the X-ray position. This star is absent (R > 20) from an archival 1989 image of the field and brightened during the 2003 outburst, and we therefore identify it as the optical counterpart of XTE J1814-338. The best source position derived from optical astrometry is R.A. = 18h13m39.s04, Dec.= -33d46m22.3s (J2000). The featureless X-ray spectrum of the pulsar in outburst is best fit by an absorbed power-law (with photon index = 1.41 +- 0.06) plus blackbody (with kT = 0.95 +- 0.13 keV) model, where the blackbody component contributes approximately 10% of the source flux. The optical broad-band spectrum shows evidence for an excess of infrared emission with respect to an X-ray heated accretion disk model, suggesting a significant contribution from the secondary or from a synchrotron-emitting region. A follow-up observation performed when XTE J1814-338 was in quiescence reveals no counterpart to a limiting magnitude of R = 23.3. This suggests that the secondary is an M3 V or later-type star, and therefore very unlikely to be responsible for the soft excess, making synchroton emission a more reasonable candidate.

We present high-resolution spectroscopy of the neutron star/low-mass X-ray binaries (LMXBs) 4U 1850-087 and 4U 0513-40 as part of our continuing study of known and candidate ultracompact binaries. The LMXB 4U 1850-087 is one of four systems in which we had previously inferred an unusual Ne/O ratio in the absorption along the line of sight, most likely from material local to the binaries. However, our recent Chandra X-ray Observatory LETGS spectrum of 4U 1850-087 finds a Ne/O ratio by number of 0.22+/-0.05, smaller than previously measured and consistent with the expected interstellar value. We propose that variations in the Ne/O ratio due to source variability, as previously observed in these sources, can explain the difference between the low- and high-resolution spectral results for 4U 1850-087. Our XMM-Newton RGS observation of 4U 0513-40 also shows no unusual abundance ratios in the absorption along the line of sight. We also present spectral results from a third candidate ultracompact binary, 4U 1822-000, whose spectrum is well fit by an absorbed power-law + blackbody model with absorption consistent with the expected interstellar value. Finally, we present the non-detection of a fourth candidate ultracompact binary, 4U 1905+000, with an upper limit on the source luminosity of < 1 x 10^{32} erg s^{-1}. Using archival data, we show that the source has entered an extended quiescent state.

We report on observations of the sixth accretion-powered millisecond pulsar, IGR J00291+5934, with the Rossi X-Ray Timing Explorer. The source is a faint, recurrent X-ray transient initially identified by INTEGRAL. The 599 Hz (1.67 ms) pulsation had a fractional rms amplitude of 8% in the 2-20 keV range, and its shape was approximately sinusoidal. The pulses show an energy-dependent phase delay, with the 6-9 keV pulses arriving up to 85 us earlier than those at lower energies. No X-ray bursts, dips, or eclipses were detected. The neutron star is in a circular 2.46 hr orbit with a very low-mass donor, most likely a brown dwarf. The binary parameters of the system are similar to those of the first known accreting millisecond pulsar, SAX J1808.4-3658. Assuming that the mass transfer is driven by gravitational radiation and that the 2004 outburst fluence is typical, the 3-yr recurrence time implies a distance of at least 4 kpc.

We report the discovery of the orbital period of the ultracompact low-mass X-ray binary (LMXB) 4U 1543-624 using time-resolved optical photometry taken with the 6.5-m Clay (Magellan II) telescope in Chile. The light curve in the Sloan r' band clearly shows a periodic, sinusoidal modulation at 18.2+-0.1 min with a fractional semiamplitude of 8%, which we identify as the binary period. This is the second shortest orbital period among all the known LMXBs, and it verifies the earlier suggestion of 4U 1543-624 as an ultracompact binary based on X-ray spectroscopic properties. The sinusoidal shape of the optical modulation suggests that it arises from X-ray heating of the mass donor in a relatively low-inclination binary, although it could also be a superhump oscillation in which case the orbital period is slightly shorter. If the donor is a C-O white dwarf as previously suggested, its likely mass and radius are around 0.03 M_sun and 0.03 R_sun, respectively. For conservative mass transfer onto a neutron star and driven by gravitational radiation, this implies an X-ray luminosity of 6.5X10^36 erg/s and a source distance of 7 kpc. We also discuss optical photometry of another LMXB, the candidate ultracompact binary 4U 1822-000. We detected significant optical variability on a time scale of about 90 min, but it is not yet clear whether this was due to a periodic modulation.

Despite considerable evidence verifying that millisecond pulsars are spun up through sustained accretion in low-mass X-ray binaries (LMXBs), it has proven surprisingly difficult to actually detect millisecond X-ray pulsars in LMXBs. There are only 5 accretion-powered millisecond X-ray pulsars known among more than 80 LMXBs containing neutron stars, but there are another 11 "nuclear powered" millisecond pulsars which reveal their spin only during brief, thermonuclear X-ray bursts. In addition, 2 of the accretion-powered pulsars also exhibit X-ray burst oscillations, and their unusual properties, along with the absence of persistent pulsations in most LMXBs, suggest that the magnetic fields in many LMXBs may be hidden by accreted material. Interestingly, the nuclear-powered pulsars offer a statistically unbiased probe of the spin distribution of recycled pulsars and show that this distribution cuts off sharply above 730 Hz, well below the breakup spin rate for most neutron star equations of state. This indicates that some mechanism acts to halt or balance spin-up due to accretion and that submillisecond pulsars must be very rare (and are possibly nonexistent). It is unclear what provides the necessary angular momentum sink, although gravitational radiation is an attractive possibility.

We propose a new non-parametric algorithm that can be implemented to study and characterize stellar clusters. The scheme attempts to simultaneously recover the stellar distribution function and the cluster potential by using projected radii and velocity information about the cluster members. The pair of these functions that is most consistent with the input data is detected by the Metropolis algorithm. In this work, the cluster characteristics recovered by CHASSIS are calibrated against the N-body realizations of two clusters, namely Hyades and Arches. The cluster mass and line-of-sight projected velocity dispersion profiles are correctly reproduced by the algorithm when the cluster obeys the assumption used in the code, namely isotropy in phase space. The results recovered by the code are shown to be insensitive to the choice of the initial parameters. The results are also not influenced by increasing the number of input data points as long as this number exceeds a minimum value which is moderately low for an input data set that obeys the assumptions of isotropy and sphericity.

(Abbrev.) We present high-resolution spectroscopy of the oxygen K-shell interstellar absorption edge in 7 X-ray binaries using the HETGS onboard Chandra. Using the brightest sources as templates, we found a best-fit model of 2 absorption edges and 5 Gaussian absorption lines. All of these features can be explained by the recent predictions of K-shell absorption from neutral and ionized atomic oxygen. We identify the K alpha and K beta absorption lines from neutral oxygen, as well as the S=3/2 absorption edge. The expected S=1/2 edge is not detected in these data due to overlap with instrumental features. We also identify the K alpha absorption lines from singly and doubly ionized oxygen. The OI K alpha absorption line is used as a benchmark with which to adjust the absolute wavelength scale for theoretical predictions of the absorption cross-sections. We find that shifts of 30-50 mA are required, consistent with differences previously noticed from comparisons of the theory with laboratory measurements. Significant oxygen features from dust or molecular components, as suggested in previous studies, are not required by our HETGS spectra. With these spectra, we can begin to measure the large-scale properties of the ISM. We place a limit on the velocity dispersion of the neutral lines of <200 km s^{-1}, consistent with measurements at other wavelengths. We also make the first measurement of the oxygen ionization fractions in the ISM. We constrain the interstellar ratio of OII/OI to ~0.1 and the ratio of OIII/OI to <0.1.

Recently we have made measurements of thermonuclear burst energetics and recurrence times which are unprecedented in their precision, largely thanks to the sensitivity of the Rossi X-ray Timing Explorer. In the "Clocked Burster", GS 1826-24, hydrogen burns during the burst via the rapid-proton (rp) process, which has received particular attention in recent years through theoretical and modelling studies. The burst energies and the measured variation of alpha (the ratio of persistent to burst flux) with accretion rate strongly suggests solar metallicity in the neutron star atmosphere, although this is not consistent with the corresponding variation of the recurrence time. Possible explanations include extra heating between the bursts, or a change in the fraction of the neutron star over which accretion takes place. I also present results from 4U 1746-37, which exhibits regular burst trains which are interrupted by "out of phase" bursts.

This paper describes test particle simulations of the response of the outer parts of Galactic disks to barring and spiral structure. Simulations are conducted for cold Mestel disks and warm quasi-exponential disks with completely flat rotation curves, subjected to pure quadrupoles and logarithmic spirals. Even though the starting velocity distributions are smooth, the end-points of the bar simulations show bimodality and multi-peaked structures at locations near the outer Lindblad resonance (OLR), although spirality can make this smoother. The growth of a bar may cause the disk isophotes to become boxy at the OLR, as stars accummulate particularly along the minor axis. The growth of a bar is also accompanied by substantial heating of the disk stars near the OLR. For the growth of a $10^{10} M_{\rm \odot}$ bar, the radial velocity dispersion is typically quadrupled for initially cold disks (initial \sigmau $\sim 10$ \kms), and typically doubled for disks with final \sigmau $\sim 45$ \kms. Simulations performed of the growth and dissolution of bars give very similar results, demonstrating that the heat once given to disk stars is very difficult to remove.

Previously, observations with the Rossi X-ray Timing Explorer showed that millisecond oscillations occur preferentially in thermonuclear X-ray bursts with photospheric radius expansion from sources rotating near 600 Hz, while they occur with equal likelihood in X-ray bursts with and without radius expansion for sources rotating near 300 Hz. With a larger sample of data than in previous studies, we find that the detectability of the oscillations is not directly determined by the properties of the X-ray bursts. Instead, we find that (1) the oscillations are observed almost exclusively when the accretion rate onto the neutron star is high, but that (2) radius expansion is only observed at high accretion rates from the 600 Hz sources, whereas it occurs only at low accretion rates in the 300 Hz sources. The persistent millisecond pulsars provide the only apparent exceptions to these trends. The first result might be explained if the oscillation amplitudes are attenuated at low accretion rates by an extended electron corona. The second result indicates that the rotation period of the neutron star determines how the burst properties vary with accretion rate, possibly through the differences in the effective surface gravity or the strength of the Coriolis force.

(abridged) An outburst of more than 80 individual bursts, similar to those seen from Soft Gamma Repeaters (SGRs), was detected from the Anomalous X-ray Pulsar (AXP) 1E 2259+586 in 2002 June. Coincident with this burst activity were gross changes in the pulsed flux, persistent flux, energy spectrum, pulse profile and spin down of the underlying X-ray source. We present RXTE and XMM-Newton observations of 1E 2259+586 that show the evolution of the aforementioned source parameters during and following this episode. Specifically, we observe an X-ray flux increase by more than an order of magnitude having two distinct components. The first component is linked to the burst activity and decays within ~2 days during which the energy spectrum is considerably harder than during the quiescent state of the source. The second component decays over the year following the glitch according to a power law in time with an exponent -0.22 +/- 0.01. The pulsed fraction decreased initially to ~15% RMS, and the pulse profile changed significantly during the outburst. A glitch was observed in 1E 2259+586 that preceded the observed burst activity. A fraction of the glitch (~19%) recovered, although the recovery was not purely exponential. An exponential rise of ~20% of the frequency jump with a time scale of ~14 days results in a significantly better fit to the data, however, contamination from a systematic drift in the phase of the pulse profile cannot be excluded. The long-term post-glitch spin-down rate decreased in magnitude relative to the pre-glitch value. A comparison with SGR outburst properties, a physical interpretation of the results, and implications on the number of magnetar candidates in our Galaxy are discussed.

We analyze 24 type I X-ray bursts from GS 1826-24 observed by the Rossi X-ray Timing Explorer between 1997 November and 2002 July. The bursts observed between 1997-98 were consistent with a stable recurrence time of 5.74 +/- 0.13 hr. The persistent intensity of GS 1826-24 increased by 36% between 1997-2000, by which time the burst interval had decreased to 4.10 +/- 0.08 hr. In 2002 July the recurrence time was shorter again, at 3.56 +/- 0.03 hr. The bursts within each epoch had remarkably identical lightcurves over the full approx. 150 s burst duration; both the initial decay timescale from the peak, and the burst fluence, increased slightly with the rise in persistent flux. The decrease in the burst recurrence time was proportional to Mdot^(-1.05+/-0.02) (where Mdot is assumed to be linearly proportional to the X-ray flux), so that the ratio alpha between the integrated persistent and burst fluxes was inversely correlated with Mdot. The average value of alpha was 41.7 +/- 1.6. Both the alpha value, and the long burst durations indicate that the hydrogen is burning during the burst via the rapid-proton (rp) process. The variation in alpha with Mdot implies that hydrogen is burning stably between bursts, requiring solar metallicity (Z ~ 0.02) in the accreted layer. We show that solar metallicity ignition models naturally reproduce the observed burst energies, but do not match the observed variations in recurrence time and burst fluence. Low metallicity models (Z ~ 0.001) reproduce the observed trends in recurrence time and fluence, but are ruled out by the variation in alpha. We discuss possible explanations, including extra heating between bursts, or that the fraction of the neutron star covered by the accreted fuel increases with Mdot.

We present high-resolution spectroscopy of the neutron-star/low-mass X-ray binaries 2S 0918-549 and 4U 1543-624 with the High Energy Transmission Grating Spectrometer onboard the Chandra X-ray Observatory and the Reflection Grating Spectrometer onboard XMM-Newton. Previous low-resolution spectra of both sources showed a broad line-like feature at 0.7 keV that was originally attributed to unresolved line emission. We recently showed that this feature could also be due to excess neutral Ne absorption, and this is confirmed by the new high-resolution Chandra spectra. The Chandra spectra are each well fit by an absorbed power-law + blackbody model with a modified Ne/O number ratio of 0.52+/-0.12 for 2S 0918-549 and 1.5+/-0.3 for 4U 1543-624, compared to the interstellar-medium value of 0.18. The XMM spectrum of 2S 0918-549 is best fit by an absorbed power-law model with a Ne/O number ratio of 0.46+/-0.03, consistent with the Chandra result. On the other hand, the XMM spectrum of 4U 1543-624 is softer and less luminous than the Chandra spectrum and has a best-fit Ne/O number ratio of 0.54+/-0.03. The difference between the measured abundances and the expected interstellar ratio, as well as the variation of the column densities of O and Ne in 4U 1543-624, supports the suggestion that there is absorption local to these binaries. We propose that the variations in the O and Ne column densities of 4U 1543-624 are caused by changes in the ionization structure of the local absorbing material. It is important to understand the effect of ionization on the measured absorption columns before the abundance of the local material can be determined. This work supports our earlier suggestion that 2S 0918-549 and 4U 1543-624 are ultracompact binaries with Ne-rich companions.

The relativistic plasma flows onto neutron stars that are accreting material from stellar companions can be used to probe strong-field gravity as well as the physical conditions in the supranuclear-density interiors of neutron stars. Plasma inhomogeneities orbiting a few kilometres above the stars are observable as X-ray brightness fluctuations on the millisecond dynamical timescale of the flows. Two frequencies in the kilohertz range dominate these fluctuations: the twin kilohertz quasi-periodic oscillations (kHz QPOs). Competing models for the origins of these oscillations (based on orbital motions) all predict that they should be related to the stellar spin frequency, but tests have been difficult because the spins were not unambiguously known. Here we report the detection of kHz QPOs from a pulsar whose spin frequency is known. Our measurements establish a clear link between kHz QPOs and stellar spin, but one not predicted by any current model. A new approach to understanding kHz QPOs is now required. We suggest that a resonance between the spin and general relativistic orbital and epicyclic frequencies could provide the observed relation between QPOs and spin.

Millisecond pulsars are neutron stars (NSs) that are thought to have been spun-up by mass accretion from a stellar companion. It is unknown whether there is a natural brake for this process, or if it continues until the centrifugal breakup limit is reached at submillisecond periods. Many NSs that are accreting from a companion exhibit thermonuclear X-ray bursts that last tens of seconds, caused by unstable nuclear burning on their surfaces. Millisecond brightness oscillations during bursts from ten NSs (as distinct from other rapid X-ray variability that is also observed) are thought to measure the stellar spin, but direct proof of a rotational origin has been lacking. Here, we report the detection of burst oscillations at the known spin frequency of an accreting millisecond pulsar, and we show that these oscillations always have the same rotational phase. This firmly establishes burst oscillations as nuclear-powered pulsations tracing the spin of accreting NSs, corroborating earlier evidence. The distribution of spin frequencies of the 11 nuclear-powered pulsars cuts off well below the breakup frequency for most NS models, supporting theoretical predictions that gravitational radiation losses can limit accretion torques in spinning up millisecond pulsars.

We have obtained optical and near-infrared images of the field of the accreting millisecond X-ray pulsar XTE J1751-305. There are no stars in the 0.7'' error circle (0.7'' is the overall uncertainty arising from tying the optical and X-ray images and from the intrinsic uncertainty in the Chandra X-ray astrometric solution). We derive limiting magnitudes for the counterpart of R>23.1, I>21.6, Z>20.6, J>19.6, K>19.2. We compare these upper limits with the magnitudes one would expect for simple models for the possible donor stars and the accretion disk subject to the reddening observed in X-rays for XTE J1751-305 and when put at the distance of the Galactic Centre (8.5 kpc). We conclude that our non-detection does not constrain any of the models for the accretion disk or possible donor stars. Deep, near-infrared images obtained during quiescence will, however, constrain possible models for the donor stars in this ultra-compact system.

We report a major outburst from the Anomalous X-ray Pulsar 1E 2259+586, in which over 80 X-ray bursts were detected in four hours using the Rossi X-ray Timing Explorer. The bursts range in duration from 2 ms to 3 s and have fluences in the 2-10 keV band that range from 3 x 10^{-11} to 5 x 10^{-9} erg cm^{-2}. We simultaneously observed increases of the pulsed and persistent X-ray emission by over an order of magnitude relative to quiescent levels. Both decayed significantly during the course of our 14 ks observation. Correlated spectral hardening was also observed, with the spectrum softening during the observation. In addition, we observed a pulse profile change, in which the amplitudes of the two peaks in the pulse profile were swapped. The profile relaxed back to its pre-outburst morphology after ~6 days. The pulsar also underwent a sudden spin-up (df/f = 4 x 10^{-6}), followed by a large (factor of \~2) increase in spin-down rate which persisted for >18 days. We also observed, using the Gemini-North telescope, an infrared enhancement, in which the K_s (2.15 um) flux increased, relative to that measured in a observation made in 2000, by a factor of ~3, three days post-outburst. The IR counterpart then faded by a factor of ~2 one week later. In addition, we report an upper limit of 50 uJy on radio emission at 1.4 GHz two days post-outburst. The X-ray properties of this outburst are like those seen only in Soft Gamma Repeaters. This conclusively unifies Anomalous X-ray Pulsars and Soft Gamma Repeaters, as predicted uniquely by the magnetar model.

We investigate the limitations of thermonuclear X-ray bursts as a distance indicator for the weakly-magnetized accreting neutron star 4U 1728-34. We measured the unabsorbed peak flux of 81 bursts in public data from the Rossi X-Ray Timing Explorer (RXTE). The distribution of peak fluxes was bimodal: 66 bursts exhibited photospheric radius expansion and were distributed about a mean bolometric flux of 9.2e-8 erg/cm^2/s, while the remaining (non-radius expansion) bursts reached 4.5e-8 erg/cm^2/s, on average. The peak fluxes of the radius-expansion bursts were not constant, exhibiting a standard deviation of 9.4% and a total variation of 46%. These bursts showed significant correlations between their peak flux and the X-ray colors of the persistent emission immediately prior to the burst. We also found evidence for quasi-periodic variation of the peak fluxes of radius-expansion bursts, with a time scale of approximately 40 d. The persistent flux observed with RXTE/ASM over 5.8 yr exhibited quasi-periodic variability on a similar time scale. We suggest that these variations may have a common origin in reflection from a warped accretion disk. Once the systematic variation of the peak burst fluxes is subtracted, the residual scatter is only approximately 3%, roughly consistent with the measurement uncertainties. The narrowness of this distribution strongly suggests that i) the radiation from the neutron star atmosphere during radius-expansion episodes is nearly spherically symmetric, and ii) the radius-expansion bursts reach a common peak flux which may be interpreted as a standard candle intensity.Adopting the minimum peak flux for the radius-expansion bursts as the Eddington flux limit, we derive a distance for the source of 4.4-4.8 kpc.

We examine the energy-resolved pulse profiles of 51 flux oscillations observed during the decline of thermonuclear X-ray bursts from accreting weakly- magnetized neutron stars with the Rossi X-ray Timing Explorer. We find that the fractional rms amplitudes of the oscillations increase as a function of energy by 0.25%/keV to 0.9%/keV between 3-20 keV, and are as large as 20% in the 13-20 keV band. We also show that the pulses observed in the higher energy bands generally lag behind those in lower energy bands by 0.002 cycles/keV to 0.007 cycles/keV between 3-20 keV. This amounts to total delays of 0.03-0.12 cycles between the lowest and highest energy bands, or time delays that range from 100-200 mus. We then model the oscillations as flux variations arising from temperature patterns on the surfaces of rapidly rotating neutron stars. In this framework, we find that the increase in the pulse amplitude with photon energy can be explained if the cooler regions on the neutron star emit in the lower energy bands, reducing the flux variations there. On the other hand, the Doppler shifts caused by the rapid rotation of the neutron star should cause the hard pulses to precede the soft pulses by about 0.05 cycles (100 mus), in contrast to the observations. This suggests that the photons originating from the stellar surface are reprocessed by a hot corona of electrons before they reach the observer.

We present an analysis of the first high-resolution spectra measured from an accretion-driven millisecond X-ray pulsar in outburst. We observed XTE J1751-305 with XMM-Newton on 2002 April 7 for approximately 35 ksec. Using a simple absorbed blackbody plus power-law model, we measure an unabsorbed flux of (6.6 +/- 0.1) * 10^(-10) erg/cm^2/s (0.5--10.0 keV). A hard power-law component (Gamma = 1.44 +/- 0.01) contributes 83% of the unabsorbed flux in the 0.5-10.0 keV band, but a blackbody component (kT = 1.05 +/- 0.01 keV) is required. We find no clear evidence for narrow or broad emission or absorption lines in the time-averaged spectra, and the sensitivity of this observation has allowed us to set constraining upper-limits on the strength of important features. The lack of line features is at odds with spectra measured from some other X-ray binaries which share some similarities with XTE J1751-305. We discuss the implications of these findings on the accretion flow geometry in XTE J1751-305.

We present the high-resolution spectrum of the accretion-powered millisecond pulsar XTE J0929-314 during its 2002 outburst, measured using the Low Energy Transmission Grating Spectrometer onboard the Chandra X-ray Observatory. The 1.5-25.3 A (0.5-8.3 keV) Chandra spectrum is well fit by a power-law + blackbody model with photon index Gamma=1.55+/-0.03, blackbody temperature kT_bb=0.65+/-0.03 keV, and blackbody normalization R_bb,km/d_10kpc=7.6+/-0.8. No emission or absorption features are found in the high-resolution spectrum, with a 3 sigma equivalent width upper limit of <0.007 A at 1.5 A and <0.12 A at 24 A. The neutral absorption edge depths are consistent with the estimated interstellar absorption along the line of sight to the source. We found no orbital modulation of the 2-10 keV X-ray flux, to a 3 sigma upper limit of 1.1%, which implies an upper limit on the binary inclination angle of i<85deg for a Roche-lobe-filling companion. We also present the broadband spectrum measured over the course of the outburst by the Rossi X-ray Timing Explorer. The XTE spectrum of XTE J0929-314 is also well fit with a power-law + blackbody model, with average values of Gamma=1.76+/-0.03, kT_bb=0.66+/-0.06 keV, and R_bb,km/d_10kpc=5.9+/-1.3 in the 2-50 keV energy range. The blackbody flux remained constant over the course of the outburst, while the power-law flux was strongly correlated to the (decreasing) flux of the source. We find that the difference in power-law photon indices measured from Chandra and XTE spectra can be explained by a change in the power-law photon index at low energies.

We present a comprehensive observational and theoretical analysis of the amplitudes and profiles of oscillations that occur during thermonuclear X-ray bursts from weakly-magnetized neutron stars in low mass X-ray binaries. Our sample contains 59 oscillations from six sources observed with the Rossi X-ray Timing Explorer. The oscillations that we examined occurred primarily during the decaying portions of bursts, and lasted for several seconds each. We find that the oscillations are as large as 15% during the declines of the bursts, and they appear and disappear due to intrinsic variations in their fractional amplitudes. However, the maxima in the amplitudes are not related to the underlying flux in the burst. We derive folded profiles for each oscillation train to study the pulse morphologies. The mean rms amplitudes of the oscillations are 5%, although the eclipsing source MXB 1659-298 routinely produces 10% oscillations in weak bursts. We also produce combined profiles from all of the oscillations from each source. Using these pulse profiles, we place upper limits on the fractional amplitudes of harmonic and half-frequency signals of 0.3% and 0.6%, respectively (95% confidence). We then compare the pulse morphologies to theoretical profiles from models with one or two antipodal bright regions on the surface of a rotating neutron star. We find that if one bright region is present on the star, it must either lie near the rotational pole or cover nearly half the neutron star in order to be consistent with the observed lack of harmonic signals. If an antipodal pattern is present, the hot regions must form very near the rotational equator. We discuss how these geometric constraints challenge current models for the production of brightness variations on the surface of a neutron star. (abridged)

We report our discovery of the likely near-infrared counterpart to the anomalous X-ray pulsar (AXP) 1E1048.1-5937, using observations from the 6.5-m Baade (Magellan-I) telescope in Chile. We derived a precise position for the X-ray source using archival data from the Chandra X-Ray Observatory. This position is inconsistent with a position reported earlier from XMM-Newton, but we show that the originally reported XMM-Newton position suffered from attitude reconstruction problems. Only two of the infrared objects in a 17 arcsec x 17 arcsec field containing the target have unusual colors, and one of these has colors consistent with those of the identified counterparts of two other AXPs. The latter object is also the only source detected within the 0.6 arcsec Chandra error circle, and we identify it as the counterpart to 1E1048.1-5937. This is the first AXP counterpart detected in multiple infrared bands, with magnitudes J=21.7(3), H=20.8(3), and K=19.4(3). There is marginal evidence for spectral flattening at longer wavelengths.

Observations of the 297 s X-ray pulsar 1E1145.1-6141 with the Rossi X-Ray Timing Explorer have revealed its 14.4 d eccentric orbit around its B supergiant companion. The best-fit orbital elements are: P_orb = 14.365(2) d, a_x sin i=99.4(18) light-s, e=0.20(3). No eclipses are detected, indicating that the binary inclination is less than 55 deg.

We used optical spectroscopy of the neutron star X-ray transient XTE J2123-058 to measure the rotational broadening of the binary companion's stellar absorption lines and determined that the companion's projected rotational velocity is vsini = 121 (+21)(-29) km/s. This value is considerably larger than the measurements of vsini obtained previously for three other neutron star systems where the values are between 34 and 65 km/s. The larger value is likely due to the combination of high binary inclination and short (6 hr) orbital period for XTE J2123-058. Along with previously measured parameters, the rotational broadening measurement allowed us to determine the binary parameters, including the ratio of the neutron star mass to the companion mass, q = M_1/M_2 = 2.7 (+1.3)(-0.9), the neutron star mass, M_1 = 1.46 (+0.30)(-0.39) solar masses, and the companion mass, M_2 = 0.53 (+0.28)(-0.39) solar masses, assuming a Roche lobe filling companion synchronously rotating with the binary orbit. These values are consistent with a previous measurement where the H alpha emission line was used to determine the semiamplitude of the neutron star's radial velocity curve (K_1). It is a significant result that the two methods give consistent values. We also report the first measurement of the XTE J2123-058 companion's radius. Assuming synchronous rotation, we obtain R_2 = 0.62 (+0.11)(-0.15) solar radii, which is in-line with that expected for a late K-type star on or close to the main sequence. Finally, we report the first precise spectroscopic determination of the binary orbital period (P_orb = 21442.3 +/- 1.8 seconds).

We analyze the frequency evolution of millisecond oscillations observed during type I X-ray bursts with the Rossi X-ray Timing Explorer in order to establish the stability of the mechanism underlying the oscillations. Our sample contains 68 pulse trains detected in a search of 159 bursts from 8 accreting neutron stars. As a first step, we confirm that the oscillations usually drift upward in frequency by about 1% toward an apparent saturation frequency. Previously noted anomalies, such as drifts toward lower frequencies as the oscillations disappear (``spin-down'' episodes) and instances of two signals present simultaneously at frequencies separated by a few Hz, occur in 5% of oscillations. Having verified the generally accepted description of burst oscillations, we proceed to study the coherence of the oscillations during individual bursts, and the dispersion in the asymptotic frequencies in bursts observed over five years. On short time scales, we find that 30% of the oscillation trains do not appear to evolve smoothly in phase. This suggests either that two signals are present simultaneously with a frequency difference too small to resolve (< 1 Hz), that the frequency evolution is discontinuous, or that discrete phase jumps occur. On time scales of years, the maximum frequencies of the oscillations exhibit fractional dispersions of < 0.004. In the case of 4U 1636-536, this dispersion is uncorrelated with the known orbital phase, which indicates that a mechanism besides orbital Doppler shifts prevents the oscillations from appearing perfectly stable. We interpret these results under the assumption that the oscillations originate from anisotropies in the emission from the surfaces of these rotating neutron stars. (Abridged)

We have identified the third known accretion-powered millisecond pulsar, XTE J0929-314, with the Rossi X-Ray Timing Explorer. The source is a faint, high-Galactic-latitude X-ray transient (d >~ 5 kpc) that was in outburst during 2002 April-June. The 185 Hz (5.4 ms) pulsation had a fractional rms amplitude of 3-7% and was generally broad and sinusoidal, although occasionally double-peaked. The hard X-ray pulses arrived up to 770 microseconds earlier than the soft X-ray pulses. The pulsar was spinning down at an average rate of -(9.2 +/- 0.4) * 10^-14 Hz/s; the spin-down torque may arise from magnetic coupling to the accretion disk, a magnetohydrodynamic wind, or gravitational radiation from the rapidly spinning pulsar. The pulsations were modulated by a 43.6 min ultracompact binary orbit, yielding the smallest measured mass function (2.7 * 10^-7 M_sun) of any stellar binary. The binary parameters imply an approximately 0.01 M_sun white dwarf donor and a moderately high inclination. We note that all three known accreting millisecond pulsars are X-ray transients in very close binaries with extremely low mass transfer rates. This is an important clue to the physics governing whether or not persistent millisecond pulsations are detected in low-mass X-ray binaries.

We demonstrate that measurable vertical structure can be excited in the accretion disc of a close binary system by a dipolar magnetic field centred on the secondary star. We present the first high resolution hydrodynamic simulations to show the initial development of a uniform warp in a tidally truncated accretion disc. The warp precesses retrogradely with respect to the inertial frame. The amplitude depends on the phase of the warp with respect to the binary frame. A warped disc is the best available explanation for negative superhumps.

We have measured the precise position of the 38-s eclipsing X-ray pulsar OAO 1657-415 with the Chandra X-Ray Observatory: RA = 17h00m48.90s, Dec = -41d39m21.6s, equninox J2000, error radius = 0.5 arcsec. Based on the previously measured pulsar mass function and X-ray eclipse duration, this 10.4-d high-mass X-ray binary is believed to contain a B supergiant companion. Deep optical imaging of the field did not detect any stars at the Chandra source position, setting a limit of V>23. However, near-IR imaging revealed a relatively bright star (J=14.1, H=11.9, K_s=10.7) coincident with the Chandra position, and we identify this star as the IR counterpart of OAO 1657-415. The IR colors and magnitudes and the optical non-detections for this star are all consistent with a highly reddened B supergiant (A_V= 20.4 +/- 1.3) at a distance of 6.4 +/- 1.5 kpc. This implies an X-ray luminosity of 3e36 erg/s (2-10 keV). IR spectroscopy can verify the spectral type of the companion and measure its radial velocity curve, yielding a neutron star mass measurement.

Low-mass X-ray binaries containing weakly magnetized neutron stars may be divided into two classes, Z and atoll sources, based in large part upon the patterns which they trace in plots of two X-ray colors. In this paper, we examine color-color diagrams of nine atoll sources and five Z sources using data from the Rossi X-ray Timing Explorer. The five-year span of data we have examined is significantly longer than those of previous studies. We find that the previous clear distinction between atoll and Z sources may be an artifact of incomplete sampling, as those atoll sources which are sampled over a wide dynamic range in intensity (greater than a factor of 80) trace three-branched color-color patterns similar to the tracks for which Z sources are named. This suggests that the two classes are more intimately related than previously realized, which also has been suggested based upon recent timing studies. However, detailed differences in the spectral behavior of the two classes argue against any simple unification.

We report on a 25 ks observation of the 8.7 s anomalous X-ray pulsar 4U~0142+61 with the High Energy Transmission Grating Spectrometer (HETGS) on the Chandra X-ray Observatory. The continuum spectrum is consistent with previous measurements and is well fit by an absorbed power-law + blackbody with photon index Gamma=3.3+/-0.4 and kT=0.418+/-0.013 keV. No evidence was found for emission or absorption lines, with an upper limit of ~50 eV on the equivalent width of broad features in the 2.5-13 A (0.95-5.0 keV) range and an upper limit of ~10 eV on the equivalent width of narrow features in the 4.1-17.7 A (0.7-3.0 keV) range. If the source is a magnetar, then the absence of a proton cyclotron line strongly constrains magnetar atmosphere models and hence the magnetic field strength of the neutron star. We find no strong features that are indicative of cyclotron absorption for magnetic field strengths of (1.9-9.8)x10^{14} G. This is still consistent with the dipole field strength of B=1.3x10^{14} G (at the polar cap) estimated from the pulsar's spindown.

We report on our search for the optical counterparts of the Southern Hemisphere anomalous X-ray pulsar 1E1048.1-5937 and the radio-quiet neutron stars in supernova remnants Puppis A, RCW 103, and PKS 1209-52. The observations were carried out with the new MIT/CfA MagIC camera on the Magellan-I 6.5 m telescope in Chile. We present deep multiband optical images of the X-ray error circles for each of these targets and discuss the resulting candidates and limits.

We report on a long-term monitoring campaign of 1E 1841-045, the 12-s anomalous X-ray pulsar and magnetar candidate at the center of the supernova remnant Kes 73. We have obtained approximately monthly observations of the pulsar with the Rossi X-ray Timing Explorer spanning over two years, during which time 1E 1841-045 is found to be rotating with sufficient stability to derive a phase-connected timing solution. A linear ephemeris is consistent with observations of the pulse period made over the last 15 yrs with the Ginga, ASCA, XTE, & SAX observatories. Phase residuals suggest the presence of "timing noise", as is typically observed from young radio pulsars. These results confirm a rapid, constant spin-down for the pulsar, which continues to maintain a steady flux; this is inconsistent with most accretion scenarios.

We report on long-term monitoring of anomalous X-ray pulsars (AXPs) using the Rossi X-ray Timing Explorer. Using phase-coherent timing, we find a wide variety of behaviors among the sources, ranging from high stability (in 1E 2259.1+586 and 4U 0142+61), to instabilities so severe that phase-coherent timing is not possible (in 1E 1048.1-5937). We note a correlation in which timing stability in AXPs decreases with increasing spin-down rate. The timing stability of soft gamma repeaters in quiescence is consistent with this trend, which is similar to one seen in radio pulsars. We consider high signal-to-noise ratio average pulse profiles as a function of energy for each AXP, and find a variety of behaviors. We find no large variability in pulse morphology nor in pulsed flux as a function of time.

We report on a 25 ks observation of the anomalous X-ray pulsar 4U 0142+61 with the High Energy Transmission Grating Spectrometer (HETGS) on the Chandra X-ray Observatory. The continuum spectrum is consistent with previous measurements and is well fit by an absorbed power-law + blackbody model with photon index Gamma=3.3+/-0.4 and blackbody temperature kT=0.418+/-0.013 keV. The pulsar frequency was 0.1150966+/-0.0000017 Hz and the pulse fractions were between 8.7% and 21%, which are also consistent with past measurements. No evidence was found for emission or absorption lines with an upper limit of ~50 eV on the equivalent width of a broad feature in the 2.5-13 A (0.95-4.96 keV) range. The absence of a proton cyclotron line strongly constrains magnetar atmosphere models and hence the magnetic field strength of the neutron star. For the energy range given above, the allowed magnetic field strengths of 4U 0142+61 are B<1.9x 10^14 G and B>9.8x 10^14 G.

We report on a 39 ks observation of the 7.7-s low-mass X-ray binary pulsar 4U1626-67 with the High Energy Transmission Grating Spectrometer (HETGS) on the Chandra X-Ray Observatory. This ultracompact system consists of a disk-accreting magnetic neutron star and a very low mass, hydrogen-depleted companion in a 42-min binary. We have resolved the previously reported Ne/O emission line complex near 1 keV into Doppler pairs of broadened (2500 km/s FWHM) lines from highly ionized Ne and O. In most cases, the blue and red line components are of comparable strength, with blueshifts of 1550-2610 km/s and redshifts of 770-1900 km/s. The lines appear to originate in hot (10^6 K), dense material just below the X-ray-heated skin of the outer Keplerian accretion disk, or else possibly in a disk wind driven from the pulsar's magnetopause. The observed photoelectric absorption edges of Ne and O appear nearly an order of magnitude stronger than expected from interstellar material and are likely formed in cool, metal-rich material local to the source. Based on the inferred local abundance ratios, we argue that the mass donor in this binary is probably the 0.02 M_sun chemically fractionated core of a C-O-Ne or O-Ne-Mg white dwarf which has previously crystallized.

We present multiband optical/infrared photometry of V4580 Sgr, the optical counterpart of the accretion-powered millisecond pulsar SAX J1808.4-3658, taken during the 1998 X-ray outburst of the system. The optical flux is consistent with emission from an X-ray-heated accretion disk. Self-consistent modeling of the X-ray and optical emission during the outburst yields best-fit extinction Av=0.68 and inclination cos(i)=0.65, assuming a distance of 2.5 kpc. The allowed inclination range requires that the pulsar's stellar companion has extremely low mass, 0.05-0.10 Msun. Some of the infrared observations are not consistent with disk emission and are too bright to be from either the disk or the companion, even in the presence of X-ray heating.

We present a detailed study of the correlations between the burst properties and the inferred mass accretion rate for the X-ray transient MXB 1659-298. The bursts which exhibited oscillations were observed when the source was at relatively high mass accretion rate, similar to what has been seen for other sources. However, due to the limited number of observations at lower mass accretion rates, no bursts were observed at such accretion rates and it is still possible that when MXB 1659-298 accretes at such low mass accretion rates, bursts can occur which might still exhibit burst oscillations. No clear correlations were found between the different burst properties and accretion rate, in contrast to what has been found for KS 1731-260 and 4U 1728-34, but similar to what has been reported for Aql X-1. However, this lack of correlation for MXB 1659-298 and Aql X-1 might be due to the limited range of mass accretion rate observed for those sources compared to KS 1731-260 and 4U 1728-34.

We measured the radial velocity curve of the companion of the neutron star X-ray transient XTE J2123-058. Its semi-amplitude (K_2) of 298.5 +/- 6.9 km/s is the highest value that has been measured for any neutron star LMXB. The high value for K_2 is, in part, due to the high binary inclination of the system but may also indicate a high neutron star mass. The mass function (f_2) of 0.684 +/- 0.047 solar masses, along with our constraints on the companion's spectral type (K5V-K9V) and previous constraints on the inclination, gives a likely range of neutron star masses from 1.2 to 1.8 solar masses. We also derive a source distance of 8.5 +/- 2.5 kpc, indicating that XTE J2123-058 is unusually far, 5.0 +/- 1.5 kpc, from the Galactic plane. Our measurement of the systemic radial velocity is -94.5 +/- 5.5 km/s, which is significantly different from what would be observed if this object corotates with the disk of the Galaxy.

There are three low-mass X-ray binaries (4U 0614+091, 2S 0918-549, and 4U 1543-624) for which broad line emission near 0.7 keV was previously reported. A recent high-resolution observation of 4U 0614+091 with the Chandra/LETGS found evidence for an unusually high Ne/O abundance ratio along the line of sight but failed to detect the previously reported 0.7 keV feature. We have made a search of archival ASCA spectra and identified a fourth source with the 0.7 keV feature, the 20-min ultracompact binary 4U 1850-087. In all four of these sources, the 0.7 keV residual is eliminated from the ASCA spectra by allowing excess photoelectric absorption due to a non-solar relative abundance of neon, just as in the LETGS spectrum of 4U 0614+091. The optical properties of these systems suggest that all four are ultracompact (P_orb<80 min) binaries. We propose that there is excess neon local to each of these sources, as also found in the ultracompact binary pulsar 4U 1626-67. We suggest that the mass donor in these systems is a low-mass, neon-rich degenerate dwarf and that the binaries are all ultracompact.

We have detected large-amplitude 0.3-1.2 mHz quasi-periodic oscillations (QPOs) from the low-mass X-ray binary pulsar 4U 1626-67/KZ TrA, using UV photometry from HST and ground-based optical photometry. These 1 mHz QPOs, which have coherence (Nu/DeltaNu)=8, are entirely distinct from the 130 mHz pulsar spin frequency, a previously known 48 mHz QPO, and the 42 min binary period (independently confirmed here). Unlike the 48 mHz and 130 mHz oscillations which are present in both the optical/UV and the X-ray emission, the 1 mHz QPOs are not detected in simultaneous observations with RXTE. The rms amplitude of the mHz QPO decreases from 15% in the far UV to 3% in the optical, while the upper limit on a corresponding X-ray QPO is as low as 0.8%. We suggest that the mHz oscillations are due to warping of the inner accretion disk. We also report the detection of coherent upper and lower sidebands of the 130 mHz optical pulsations, with unequal amplitude and a spacing of 1.93 mHz around the main pulsation. The origin of these sidebands remains unclear.

We report on long-term monitoring of the anomalous X-ray pulsar (AXP) 1E 1048.1-5937 using the Rossi X-ray Timing Explorer (RXTE). The timing behavior of this pulsar is different from that of other AXPs being monitored with RXTE. In particular, we show that the pulsar shows significant deviations from simple spin-down such that phase-coherent timing has not been possible over time spans longer than a few months. We find that the deviations from simple spin down are not consistent with single "glitch" type events, nor are they consistent with radiative precession. We show that in spite of the rotational irregularities, the pulsar exhibits neither pulse profile changes nor large pulsed flux variations. We discuss the implications of our results for AXP models. In the context of the magnetar model, we suggest that 1E 1048.1-5937 may be a transition object between the soft gamma-ray repeater and AXP populations, and the AXP most likely to one day undergo an outburst.

We report on long-term monitoring of the anomalous X-ray pulsar 1E 1048.1-5937 using the Rossi X-ray Timing Explorer. This pulsar's timing behavior is different from that of other AXPs. In particular, the pulsar shows significant deviations from simple spin-down such that phase-coherent timing has not been possible over time spans longer than a few months. We show that in spite of the rotational irregularities, the pulsar exhibits neither pulse profile changes nor large pulsed flux variations. We discuss the implications of our results for AXP models. We suggest that 1E 1048.1-5937 may be a transition object between the soft gamma-ray repeater and AXP populations, and the AXP most likely to one day undergo an outburst.

We report on the highly variable and erratic long-term X-ray luminosity behavior of the only known accretion-driven millisecond X-ray pulsar SAX J1808.4-3658 during its 2000 outburst, as observed with RXTE. The maximum observed luminosity is ~2.5 x 10^35 erg s^-1 (3-25 keV; for a distance of 2.5 kpc), which is a factor of ~10 lower than that observed during the 1996 and 1998 outbursts. Due to solar constraints, the source could not be observed for several months with RXTE before 21 Jan. 2000. Therefore, the exact moment of the outburst onset is unknown and the peak luminosity could have been significantly higher. On some occasions SAX J1808.4-3658 was observed with luminosities of ~10^35 erg s^-1, but on other occasions it could not be detected with RXTE resulting in upper limits of a few times 10^33 erg s^-1 (3-25 keV). The non-detections of the source during its 2000 outburst obtained with BeppoSAX demonstrate that its luminosity was at times <10^32 erg s^-1 (0.5-10 keV). However, only a few days after these BeppoSAX observations, we detected the source again with RXTE at high luminosities, giving a factor of >1000 of luminosity swings in this system on time scales of days. The last detection of SAX J1808.4-3658 with RXTE was on 2000 May 13, almost 4 months after the first detection. Due to the lack of sensitivity and observations during the 1996 and 1998 outbursts, it cannot be excluded that after those outbursts the source remained active for months and that the source behavior during the 2000 outburst is not unique. Long duration activity at low luminosities has been observed in other transients (both neutron stars and black holes), although not with such extreme variability which might point to a different origin for this behavior for the pulsar.

We use archival data from the Rossi X-Ray Timing Explorer to examine 125 type I X-ray bursts from the 9 weakly magnetic accreting neutron stars where millisecond oscillations have been detected during some bursts. We find that oscillations from the 6 "fast" (approximately 600 Hz) sources are almost always observed during radius expansion bursts, whereas oscillations from the 3 "slow" (about 300 Hz) sources are about equally likely to be found in bursts both with and without photospheric radius expansion. This strongly suggests that the distinction between these two source groups cannot be an observational selection effect, but must instead arise from some physical mechanism.

BeppoSAX has revealed a population of faint neutron star X-ray transients in the Galactic bulge. King conjectured that these neutron stars are accreting from brown dwarfs (BD) with a time-averaged mass transfer rate <\dot M>\approx 10^{-11} M_sun/year. We show that the orbital parameters of the 401 Hz accreting millisecond pulsar SAX J1808.4-3658 support this hypothesis. A main-sequence mass donor requires a nearly face-on inclination and a higher <\dot M> than observed. However, the range of allowed inclinations is substantially relaxed and the predicted mass transfer rate is consistent with that observed if a hot 0.05 M_sun BD is the donor. We then explain the BD radius required (0.13 R_\sun) to fill the Roche lobe. All transiently accreting neutron stars have a minimum luminosity in quiescence set by the time-averaged mass transfer rate onto the neutron star. We show that the constant heating of the BD by this emission can maintain the higher entropy implied by a 0.13 R_\sun radius. Our work bolsters the case that SAX J1808.4-3658 is a progenitor to compact ms radio pulsar binaries (e.g. like those found in 47 Tuc). The very low <\dot M>'s implies that the progenitors to these radio pulsars are long-lived (Gyrs) transient systems, rather than short-lived (Myrs) Eddington-limited accretors. Hence, the accreting progenitor population to millisecond radio pulsars in 47 Tuc could still be present and found in quiescence with Chandra.(abridged)

We report the first extensive set of optical photometric observations of the counterpart to SAX J1808.4-3658 (V4580 Sagittarius) in quiescence. The source was detected at V~21, 5 magnitudes fainter than at the peak of its 1998 outburst. However, a comparable ~6% semi-amplitude 2hr modulation of its flux is revealed. This has the same phasing and approximately sinusoidal modulation as seen during outburst and with photometric minimum when the pulsar is behind the companion. The lack of a double-humped morphology rules out an ellipsoidal origin, implying that the bulk of the optical flux does not arise from the companion. Moreover, applying crude modelling to the disc and X-ray irradiated face of the donor shows that the internal energy release of a remnant disc (with mass transfer driven by gravitational radiation) is sufficient to explain most of the optical emission, and with the modulation due to the varying contribution of the heated star's face. We note that this model is also consistent with the much lower X-ray to optical flux ratio in quiescence versus outburst, and with the phasing of the optical modulation.

We estimate the mass of the central black hole in our Galaxy from stellar kinematical data published by Ghez et al. (1998) and Genzel et al. (2000). For this we develop a method, related to Merritt (1993), for non-parametrically reconstructing the mass profile and the stellar distribution function in the central region of the Galaxy from discrete kinematic data, including velocity errors. Models are searched using the Metropolis algorithm. We assume that the mass distribution is spherical and the stellar velocity distribution is isotropic, and devise a test of this assumption. From proper motions we obtain an enclosed mass of $2.0\pm{0.7}\times10^6{\rm M}_{\odot}$ within the inner $0.0044\rm pc$, from radial velocities we obtain a mass of $2.2^{+1.6}_{-1.0}\times10^6{\rm M}_{\odot}$ within 0.046pc and from three-dimensional velocities we obtain $1.8^{+0.4}_{-0.3}\times10^6{\rm M}_{\odot}$ within 0.046pc.

The energy source powering the X-ray emission from anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) is still uncertain. In one scenario, the presence of an ultramagnetized neutron star, or ``magnetar'', with B on the order of 10^{14} - 10^{15} G is invoked. To investigate this hypothesis, we have analyzed archival ASCA data for several known AXPs and SGRs, and fitted them with a model in which all or part of the X-ray flux originates as thermal emission from a magnetar. Our magnetar spectral model includes the effects of the anisotropy of the heat flow through an ultramagnetized neutron star envelope, reprocessing by a light element atmosphere, and general relativistic corrections to the observed spectrum. We obtain good fits to the data with radii for the emitting areas which are generally consistent with those expected for neutron stars, in contrast to blackbody (BB) fits, which imply much smaller radii. Furthermore, the inclusion of atmospheric effects results in inferred temperatures which are lower than those implied by BB fits, but however still too high to be accounted by thermal cooling alone. An extra source of heating (possibly due to magnetic field decay) is needed. Despite the harder tail in the spectrum produced by reprocessing of the outgoing flux through the atmosphere, spectral fits still require a considerable fraction of the flux to be in a power-law component.

We report the discovery of a highly coherent oscillation in a type-I X-ray burst observed from 4U 1916-053 by the Rossi X-ray Timing Explorer (RXTE). The oscillation was most strongly detected approx. 1 s after the burst onset at a frequency of 269.3 Hz, and it increased in frequency over the following 4 seconds of the burst decay to a maximum of around 272 Hz. The total measured drift of 3.58 +/- 0.41 Hz (1 sigma) represents the largest fractional change in frequency (1.32 +/- 0.15 %) yet observed in any burst oscillation. If the asymptotic frequency of the oscillation is interpreted in terms of a decoupled surface burning layer, the implied neutron star spin period is around 3.7 ms. However, the expansion of the burning layer required to explain frequency drift during the burst is around 80 m, substantially larger than expected theoretically (assuming rigid rotation). The oscillation was not present in the persistent emission before the burst, nor in the initial rise. When detected its amplitude was 6-12% (RMS) with a roughly sinusoidal profile. The burst containing the oscillation showed no evidence for photospheric radius expansion, while at least 5 of the other 9 bursts observed from the source by RXTE during 1996 and 1998 did. No comparable oscillations were detected in the other bursts. A pair of kilohertz quasi-periodic oscillations (QPOs) has been previously reported from this source with a mean separation of 348 +/- 12 Hz. 4U 1916-053 is the first example of a source where the burst oscillation frequency is significantly smaller than the frequency separation of the kHz QPOs.

The spectacular first light observation by the Chandra X-Ray Observatory revealed an X-ray point source near the center of the 300 yr old Cas A supernova remnant. We present an analysis of the public X-ray spectral and timing data. No coherent pulsations were detected in the Chandra/HRC data. The 3-sigma upper limit on the pulsed fraction is <35% for P>20 ms. The Chandra/ACIS spectrum of the point source may be fit with an ideal blackbody (kT=0.5 keV), or with BB models modified by the presence of a NS atmosphere (kT=0.25-0.35 keV), but the temperature is higher and the inferred emitting area lower than expected for a 300 yr old NS according to standard cooling models. The spectrum may also be fit with a power law model (photon index 2.8-3.6). Both the spectral properties and the timing limits of the point source are inconsistent with a young Crab-like pulsar, but are quite similar to the properties of the anomalous X-ray pulsars. The spectral parameters are also very similar to those of the other radio-quiet X-ray point sources in the supernova remnants Pup A, RCW 103, and PKS 1209-52. Current limits on an optical counterpart for the Cas A point source rule out models that invoke fallback accretion onto a compact object if fallback disk properties are similar to those in quiescent low-mass X-ray binaries. However, the optical limits are marginally consistent with plausible alternative assumptions for a fallback disk. In this case, accreting NS models can explain the X-ray data, but an accreting BH model is not promising.

We report the detection of a sudden spin-up of the 11 s anomalous X-ray pulsar 1RXS J170849.0-4000910 in regular timing observations made with the Rossi X-ray Timing Explorer. The event, which occurred between MJD 51446 (1999 September 25) and 51472 (1999 October 21), is well characterized by an increase in the rotational frequency of magnitude (Delta nu / nu) = (6.2 +/- 0.3) x 10^-7 and an increase in the rate of spin down (Delta \dot\nu / \dot\nu) = (1.38 +/- 0.25) x 10^-2. These values are very similar to those of glitches observed in the Vela radio pulsar and other young radio pulsars. The event therefore suggests that the internal structure of this anomalous X-ray pulsar is similar to those of the radio pulsars. In particular, it implies that the fractional moment of inertia in neutron superfluid that is not corotating with the crust is >= 1 %. The detection of a glitch in this anomalous X-ray pulsar constrains models for the origin of glitches in neutron stars. Most notably, it challenges models that preclude glitches in long-period pulsars, and, under the magnetar hypothesis, suggests that large glitches can occur in hot neutron stars. The glitch is consistent with the predictions of the magnetar model for anomalous X-ray pulsars, but accretion-powered scenarios cannot be excluded using our observations alone.

The new transient X-ray pulsar XTE J0111.2-7317 was observed with Advanced Satellite for Cosmology and Astrophysics (ASCA) on 1998 November 18, a few days after its discovery with the Proportional Counter Array onboard the Rossi X-ray Timing Explorer. The source was detected at a flux level of 3.6x10^-10 erg cm^-2 s^-1 in the 0.7--10.0 keV band, which corresponds to the X-ray luminosity of 1.8x10^38 erg s^-1, if a distance of 65 kpc for this pulsar in the Small Magellanic Cloud is assumed. Nearly sinusoidal pulsations with a period of 30.9497 +/- 0.0004 s were unambiguously detected during the ASCA observation. The pulsed fraction is low and slightly energy dependent with average value of \~27%. The energy spectrum shows a large soft excess below ~2 keV when fitted to a simple power-law type model. The soft excess is eliminated if the spectrum is fitted to an ``inversely broken power-law'' model, in which photon indices below and above a break energy of 1.5 keV are 2.3 and 0.8, respectively. The soft excess can also be described by a blackbody or a thermal bremsstrahlung when the spectrum above ~2 keV is modeled by a power-law. In these models, however, the thermal soft component requires a very large emission zone, and hence it is difficult to explain the observed pulsations at energies below 2 keV. A bright state of the source enables us to identify a weak iron line feature at 6.4 keV with an equivalent width of 50 +/- 14 eV. Pulse phase resolved spectroscopy revealed a slight hardening of the spectrum and marginal indication of an increase in the iron line strength during the pulse maximum.

A deep all-sky imaging hard x-ray survey and wide-field monitor is needed to extend soft (ROSAT) and medium (ABRIXAS2) x-ray surveys into the 10-100 keV band (and beyond) at comparable sensitivity (~0.05 mCrab). This would enable discovery and study of >3000 obscured AGN, which probably dominate the hard x-ray background; detailed study of spectra and variability of accreting black holes and a census of BHs in the Galaxy; Gamma-ray bursts and associated massive star formation (PopIII) at very high redshift and Soft Gamma-ray Repeaters throughout the Local Group; and a full galactic survey for obscured supernova remnants. The Energetic X-ray Imaging Survey Telescope (EXIST) is a proposed array of 8 x 1m^2 coded aperture telescopes fixed on the International Space Station (ISS) with 160deg x 40deg field of view which images the full sky each 90 min orbit. EXIST has been included in the most recent NASA Strategic Plan as a candidate mission for the next decade. An overview of the science goals and mission concept is presented.

The USA Experiment is a new X-ray timing experiment with large collecting area and microsecond time resolution capable of conducting a broad program of studies of galactic X-ray binaries. USA is one of nine experiments aboard the Advanced Research and Global Observation Satellite which was launched February 23, 1999. USA is a collimated proportional counter X-ray telescope with about 1000 cm$^2$ of effective area per detector with two detectors sensitive to photons in the energy range 1-15 keV. A unique feature of USA is that photon events are time tagged by reference to an onboard GPS receiver allowing precise absolute time and location determination. We will present an overview of the USA instrument, capabilities, and scientific observing plan as well as the current status of the instrument.

We report on long-term X-ray timing of two anomalous X-ray pulsars, 1RXS J170849.0-400910 and 1E 2259+586, using the Rossi X-ray Timing Explorer. In monthly observations made over 1.4 yr and 2.6 yr for the two pulsars, respectively, we have obtained phase-coherent timing solutions which imply that these objects have been rotating with great stability throughout the course of our observations. For 1RXS J170849.0-400910, we find a rotation frequency of 0.0909169331(5) Hz and frequency derivative -15.687(4) x 10^(-14) Hz/s, for epoch MJD 51215.931. For 1E 2259+586, we find a rotation frequency of 0.1432880613(2)Hz, and frequency derivative -1.0026(7) x 10^(-14) Hz/s, for epoch MJD 51195.583. RMS phase residuals from these simple models are only about 0.01 cycles for both sources. We show that the frequency derivative for 1E 2259+586 is inconsistent with that inferred from incoherent frequency observations made over the last 20 yr. Our observations are consistent with the magnetar hypothesis and make binary accretion scenarios appear unlikely.

The recent discovery of the first known accretion-powered millisecond pulsar with the Rossi X-Ray Timing Explorer provides the first direct probe of the interaction of an accretion disk with the magnetic field of a weakly magnetic (B<10^10 G) neutron star. We demonstrate that the presence of coherent pulsations from a weakly magnetic neutron star over a wide range of accretion rates places strong constraints on models of the disk-magnetosphere interaction. We argue that the simple Mdot^(3/7) scaling law for the Keplerian frequency at the magnetic interaction radius, widely used to model disk accretion onto magnetic stars, is not consistent with observations of SAX J1808.4-3658 for most proposed equations of state for stable neutron stars. We show that the usually neglected effects of multipole magnetic moments, radiation drag forces, and general relativity must be considered when modeling such weakly magnetic systems. Using only very general assumptions, we obtain a robust estimate of mu~(1-10)x10^26 G cm^3 for the dipole magnetic moment of SAX J1808.4-3658, implying a surface dipole field of ~10^8-10^9 G at the stellar equator. We therefore infer that after the end of its accretion phase, this source will become a normal millisecond radio pulsar. Finally, we compare the physical properties of this pulsar to those of the non-pulsing, weakly magnetic neutron stars in low-mass X-ray binaries and argue that the absence of coherent pulsations from the latter does not necessarily imply that these neutron stars have significantly different magnetic field strengths from SAX J1808.4-3658.

We present observations of a series of 10 outbursts of pulsed hard X-ray flux from the transient 10.6 mHz accreting pulsar GS 1843-02, using the Burst and Transient Source Experiment on the Compton Gamma Ray Observatory. These outbursts occurred regularly every 242 days, coincident with the ephemeris of the periodic transient GRO J1849-03 (Zhang et al. 1996), which has recently been identified with the SAS 3 source 2S 1845-024 (Soffitta et al. 1998). Our pulsed detection provides the first clear identification of GS 1843-02 with 2S 1845-024. We present a pulse timing analysis which shows that the 2S 1845-024 outbursts occur near the periastron passage of the neutron star's highly eccentric (e = 0.88+-0.01) 242.18+-0.01 day period binary orbit about a high mass (M > 7 solar masses) companion. The orbit and transient outburst pattern strongly suggest the pulsar is in a binary system with a Be star. Our observations show a long-term spin-up trend, with most of the spin-up occurring during the outbursts. From the measured spin-up rates and inferred luminosities we conclude that an accretion disk is present during the outbursts.

Typical radio pulsars are magnetized neutron stars that are born rapidly rotating and slow down as they age on time scales of 10 to 100 million years. However, millisecond radio pulsars spin very rapidly even though many are billions of years old. The most compelling explanation is that they have been "spun up" by the transfer of angular momentum during accretion of material from a companion star in so-called low-mass X-ray binary systems, LMXBs. (LMXBs consist of a neutron star or black hole accreting from a companion less than one solar mass.) The recent detection of coherent X-ray pulsations with a millisecond period from a suspected LMXB system appears to confirm this link. Here we report observations showing that the orbital period of this binary system is two hours, which establishes it as an LMXB. We also find an apparent modulation of the X-ray flux at the orbital period (at the two per cent level), with a broad minimum when the pulsar is behind this low-mass companion star. This system seems closely related to the "black widow" millisecond radio pulsars, which are evaporating their companions through irradiation. It may appear as an eclipsing radio pulsar during periods of X-ray quiescence.

Enigmatic transitions between spin-up and spin-down have been observed in several X-ray pulsars accreting matter via an accretion disk. In these transitions, the torque changes sign but remains at nearly the same magnitude. It has been noted previously that alternating prograde and retrograde disk flows would explain many features of the torque reversals, although it has been unclear how a stable retrograde disk could be formed. We suggest that the reversals may be related to the disk at times being warped to such an extent that the inner region becomes tilted by more than 90 degrees. This region would thus become retrograde, leading to a negative torque. Accretion disk models can show such behavior, if account is taken of a warping instability due to irradiation. The resulting `flip-overs' of the inner parts of the disk can reproduce most characteristics of the observations, although it remains unclear what sets the timescale on which the phenomenon occurs. If this model were correct, it would have a number of ramifications, for instance that in the spin-down state the X-ray source would mostly be observed through the accretion disk.

There are now about fifty known radio pulsars in binary systems, including at least five in double neutron star binaries. In some cases, the stellar masses can be directly determined from measurements of relativistic orbital effects. In others, only an indirect or statistical estimate of the masses is possible. We review the general problem of mass measurement in radio pulsar binaries, and critically discuss all current estimates of the masses of radio pulsars and their companions. We find that significant constraints exist on the masses of twenty-one radio pulsars, and on five neutron star companions of radio pulsars. All the measurements are consistent with a remarkably narrow underlying gaussian mass distribution, m = 1.35 +- 0.04 solar masses. There is no evidence that extensive mass accretion (delta m >~ 0.1 solar mass) has occurred in these systems. We also show that the observed inclinations of millisecond pulsar binaries are consistent with a random distribution, and thus find no evidence for either alignment or counteralignment of millisecond pulsar magnetic fields.

We report on high-speed optical photometry of the radio positions of two young rotation-powered pulsars. No pulsations were detected from the optical counterpart proposed by Caraveo et al. (1994) for PSR B1509-58, with a 2 sigma upper limit on the pulsed fraction of 12 percent, significantly lower than that measured in the five known optical pulsars. Given its low pulsed fraction, high optical luminosity, and significant (8 percent) chance coincidence probability, we suggest that this candidate is not associated with the pulsar. We also find that the still-unidentified optical counterpart of PSR B1706-44 has R>18 and lies within 3 arcsec of an R=16.6 star.

GRO J2058+42, a transient 198 second x-ray pulsar, was discovered by the Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory (CGRO), during a "giant" outburst in 1995 September-October. The total flux peaked at about 300 mCrab (20-50 keV) as measured by Earth occultation. The pulse period decreased from 198 s to 196 s during the 46-day outburst. The pulse shape evolved over the course of the outburst and exhibited energy dependent variations. BATSE observed five additional weak outbursts from GRO J2058+42, each with two week duration and peak pulsed flux of about 15 mCrab (20-50 keV), that were spaced by about 110 days. An observation of the 1996 November outburst by the Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA) localized the source to within a 4\arcmin radius error circle (90 % confidence) centered on R.A. = 20h 59m.0, Decl. = 41 deg 43 arcmin (J2000). Additional shorter outbursts with peak pulsed fluxes of about 8 mCrab were detected by BATSE halfway between the first four 15 mCrab outbursts. The RXTE All-Sky Monitor detected all 8 weak outbursts with approximately equal durations and intensities. GRO J2058+42 is most likely a Be/X-ray binary that appears to outburst at periastron and apastron. No optical counterpart has been identified to date and no x-ray source was present in the error circle in archival ROSAT observations.

We present infrared spectroscopy of the low-mass X-ray binary GX 1+4/V2116 Oph. This symbiotic binary consists of a 2-min accretion-powered pulsar and an M5 III red giant. A strong He I 1.083 micron emission line with a pronounced P Cygni profile was observed. From the blue edge of this feature, we infer an outflow velocity of 250(50) km/s. This is an order of magnitude faster than a typical red giant wind, and we suggest that radiation from the accretion disk or the neutron star may contribute to the acceleration of the outflow. We infer a wind mass loss rate of around 10^-6 Msun/yr. Accretion from such a strong stellar wind provides a plausible alternative to Roche lobe overflow for supplying the accretion disk which powers the X-ray source. The H I Paschen beta and He I 1.083 micron lines showed no evidence for the dramatic changes previously reported in some optical lines, and no evidence for pulsations at the 2-min pulsar period.

The low-mass X-ray binary pulsar 4U 1626-67 shows 0.048 Hz quasi-periodic oscillations (QPOs) and red noise variability as well as coherent pulsations at the 0.130 Hz neutron star spin frequency. Power density spectra of observations made with the Rossi X-ray Timing Explorer show significant sidebands separated from the pulsar spin frequency (and its harmonics) by the QPO frequency. These show that the instantaneous amplitude of the coherent pulsations is modulated by the amplitude of the QPOs. This phenomenon is expected in models such as the magnetospheric beat frequency model where the QPOs originate near the polar caps of the neutron star. In the 4--8 keV energy range, however, the lower-frequency sidebands are significantly stronger than their higher-frequency complements; this is inconsistent with the magnetospheric beat frequency model. We suggest that the 0.048 Hz QPOs are instead produced by a structure orbiting the neutron star at the QPO frequency. This structure crosses the line of sight once per orbit and attenuates the pulsar beam, producing the symmetric (amplitude modulation) sidebands. It also reprocesses the pulsar beam at the beat frequencies between the neutron star spin frequency and the QPOs, producing the excess variability observed in the lower-frequency sidebands. Quite independently, we find no evidence that the red noise variability modulates the amplitude of the coherent pulsations. This is also in contrast to the expectations of the magnetospheric beat frequency model and differs from the behavior in some high-mass X-ray binary pulsars.

Dramatic torque reversals between spin up and spin down have been observed in half of the persistent X-ray pulsars monitored by the BATSE all-sky monitor on CGRO. Theoretical models developed to explain early pulsar timing data can explain spin down torques via a disk-magnetosphere interaction if the star nearly corotates with the inner accretion disk. To produce the observed BATSE torque reversals, however, these equilibrium models require the disk to alternate between two mass accretion rates, with $\dot M_{\pm}$ producing accretion torques of similar magnitude, but always of opposite sign. Moreover, in at least one pulsar (GX 1+4) undergoing secular spin down the neutron star spins down faster during brief ($\sim 20$ day) hard X-ray flares -- this is opposite the correlation expected from standard theory, assuming BATSE pulsed flux increases with mass accretion rate. The $10$ day to 10 yr intervals between torque reversals in these systems are much longer than any characteristic magnetic or viscous time scale near the inner disk boundary and are more suggestive of a global disk phenomenon. We discuss possible explanations of the observed torque behavior. Despite the preferred sense of rotation defined by the binary orbit, the BATSE observations are surprisingly consistent with an earlier suggestion by Makishima \etal (1988) for GX~1+4: the disks in these systems somehow alternate between episodes of prograde and retrograde rotation. We are unaware of any mechanism that could produce a stable retrograde disk in a binary undergoing Roche-lobe overflow, but such flip-flop behavior does occur in numerical simulations of wind-fed systems. One possibility is that the disks in some of these binaries are fed by an X-ray excited wind.

We summarize five years of continuous monitoring of accretion-powered pulsars with the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory. Our 20-70 keV observations have determined or refined the orbital parameters of 13 binaries, discovered 5 new transient accreting pulsars, measured the pulsed flux history during outbursts of 12 transients (GRO J1744-28, 4U 0115+634, GRO J1750-27, GS 0834-430, 2S 1417-624, GRO J1948+32, EXO 2030+375, GRO J1008-57, A 0535+26, GRO J2058+42, 4U 1145-619 and A 1118-616), and also measured the accretion torque history of during outbursts of 6 of those transients whose orbital parameters were also known. We have also continuously measured the pulsed flux and spin frequency for eight persistently accreting pulsars (Her X-1, Cen X-3, Vela X-1, OAO 1657-415, GX 301-2, 4U 1626-67, 4U 1538-52, and GX 1+4). Because of their continuity and uniformity over a long baseline, BATSE observations have provided new insights into the long-term behavior of accreting magnetic stars. We have found that all accreting pulsars show stochastic variations in their spin frequencies and luminosities, including those displaying secular spin-up or spin-down on long time scales, blurring the conventional distinction between disk-fed and wind-fed binaries. Pulsed flux and accretion torque are strongly correlated in outbursts of transient accreting pulsars, but uncorrelated, or even anticorrelated, in persistent sources.

We present multiwavelength observations of this S-type symbiotic LMXB which consists of a 2-min X-ray pulsar accreting from an M6 III giant. This is the only symbiotic system definitely known to contain a neutron star. The steady interstellar extinction toward the binary (Av=5) contrasts the variable N_H inferred from X-ray measurements, most likely evidence for a stellar wind. The mass donor is probably near the tip of the first-ascent red giant branch, in which case the system is 3-6 kpc distant and has an X-ray luminosity of 10^37 erg/s. It is also possible, though less likely, that the donor star is just beginning its ascent of the asymptotic giant branch, in which case the system is 12-15 kpc distant and has an X-ray luminosity of 10^38 erg/s. However, our measured Av argues against such a large distance. We show that the dense (10^9 cm^-3) emission-line nebula enshrouding the binary is powered by UV radiation from an accretion disk. The emission-line spectrum constrains the temperature and inner radius of the disk (and thus the pulsar's magnetic field strength), and we discuss this in the context of the accretion torque reversals observed in the pulsar. We also show that the binary period must be >100 d and is most likely >260 d, making GX 1+4 the only known LMXB with Porb>10 d. If the mass donor fills its Roche lobe, the mass transfer must be highly super-Eddington, requiring much mass loss from the binary. We discuss the alternative that the disk forms from the slow, dense stellar wind expected from the red giant.

Rapid UBVRI photometry of this ultracompact LMXB pulsar has detected 7.67 s optical pulsations in all five bands. The optical pulsations, which are at the same frequency as the X-ray pulsations, are understood as reprocessing of pulsed X-rays in the accretion disk or on the companion surface. The optical pulsed fraction is 6%, independent of wavelength, indicating that the optical emission is dominated by X-ray reprocessing. A weaker (1.5%) sideband, shifted down 0.395(15) mHz from the main optical pulsation, is also present. This is consistent with a previously reported sideband and corroborating the 42-min binary period proposed earlier by Middleditch et al. (1981). A 0.048 Hz optical QPO, corresponding to a known X-ray feature, was also detected in some of the observations. This is the first measurement of an optical QPO in an X-ray binary pulsar. I discuss constraints on the nature of the mass donor and show that mass transfer via a radiatively-driven wind is inconsistent with the data. I also review the basic theory of X-ray-heated accretion disks and show that such models provide a good fit to the optical photometry. If the X-ray albedo of LMXB accretion disks is as high as recently reported, then the optical data imply a distance of 8 kpc and an X-ray luminosity of 10^37 erg/s.

Over five years of daily hard X-ray (>20 keV) monitoring of the 2-min accretion-powered pulsar GX 1+4 with the Compton Gamma Ray Observatory/BATSE large-area detectors has found nearly continuous rapid spin-down, interrupted by a bright 200-d spin-up episode. During spin-down, the torque becomes more negative as the luminosity increases (assuming that the 20-60 keV pulsed flux traces bolometric luminosity), the opposite of what is predicted by standard accretion torque theory. No changes in the shape of the 20-100 keV pulsed energy spectrum were detected, so that a very drastic change in the spectrum below 20 keV or the pulsed fraction would be required to make the 20-60 keV pulsed flux a poor luminosity tracer. These are the first observations which flatly contradict standard magnetic disk accretion theory, and they may have important implications for understanding the spin evolution of X-ray binaries, cataclysmic variables, and protostars. We briefly discuss the possibility that GX 1+4 may be accreting from a retrograde disk during spin-down, as previously suggested.

We present results from a deep search for radio emission from the X-ray pulsar systems GX 1+4, GS 0834-430 and 4U 0115+63 which have variously been suggested to possess radio jets and to be good candidates for propellor ejection mechanisms. None of these sources is detected at their optical positions, to 3 sigma limits of a few hundred micro-Jy. This places upper limits on their radio luminosities and thus on the internal energy and numbers of any relativistic electrons which are three to four orders of magnitude below those of radio-bright X-ray binaries such as SS 433 and Cyg X-3. Spectral and structural information on the proposed `radio lobes' of GX 1+4 make their association with the source unlikely. The lack of detected radio emission from any X-ray pulsar system is discussed statistically, and it is found that X-ray pulsations and radio emission from X-ray binaries are strongly anti-correlated.