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.

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.

The increasing richness of data related to cold dense matter, from laboratory experiments to neutron-star observations, requires a framework for constraining the properties of such matter that makes use of all relevant information. Here, we present a rigorous but practical Bayesian approach that can include diverse evidence, such as nuclear data and the inferred masses, radii, tidal deformabilities, moments of inertia, and gravitational binding energies of neutron stars. We emphasize that the full posterior probability distributions of measurements should be used rather than, as is common, imposing a cut on the maximum mass or other quantities. Our method can be used with any parameterization of the equation of state (EOS). We use both a spectral parameterization and a piecewise polytropic parameterization with variable transition densities to illustrate the implications of current measurements and show how future measurements in many domains could improve our understanding of cold catalyzed matter. We find that different types of measurements will play distinct roles in constraining the EOS in different density ranges. For example, better symmetry energy measurements will have a major influence on our understanding of matter somewhat below nuclear saturation density but little influence above that density. In contrast, precise radius measurements or multiple tidal deformability measurements of the quality of those from GW170817 or better will improve our knowledge of the EOS over a broader density range.

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.

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.

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.

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

This erratum addresses an error in the text of Lo et al. (2013) and more importantly corrects a computational error that was made in all previously published relativistic computations of absolute pulse waveform fluxes produced by hot spots on rotating stars that followed the radiation forward in time from the star to the observer. These computations used an incorrect expression for the differential emitting area on the stellar surface, as measured in the comoving frame. This erratum explains this error, lists all the published papers that contain results affected by this error, and summarizes in a general way why and how the results in Lo et al. (2013) were affected by this error, in a way that makes clear why and how the results in the other papers were affected. The idealized radiation spectra, beaming patterns, and spot shapes assumed in many of the previously reported computations of pulse waveforms very likely introduced systematic errors that are much larger than the relatively small errors introduced by using an incorrect expression for the differential emitting area. This erratum also corrects an error in the caption of Figure 5, an error in the text on page 19, and an incomplete description of our results on page 24.

Precise and reliable measurements of the masses and radii of neutron stars with a variety of masses would provide valuable guidance for improving models of the properties of cold matter with densities above the saturation density of nuclear matter. Several different approaches for measuring the masses and radii of neutron stars have been tried or proposed, including analyzing the X-ray fluxes and spectra of the emission from neutron stars in quiescent low-mass X-ray binary systems and thermonuclear burst sources; fitting the energy-dependent X-ray waveforms of rotation-powered millisecond pulsars, burst oscillations with millisecond periods, and accretion-powered millisecond pulsars; and modeling the gravitational radiation waveforms of coalescing double neutron star and neutron star -- black hole binary systems. We describe the strengths and weaknesses of these approaches, most of which currently have substantial systematic errors, and discuss the prospects for decreasing the systematic errors in each method.

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]

We describe sophisticated new Bayesian analysis methods that make it possible to estimate quickly the masses and radii of rapidly rotating, oblate neutron stars by fitting oblate-star waveform models to energy-resolved observations of the X-ray oscillations produced by a hot spot on such stars. We conclude that models that take the oblate shape of the star into account should be used for stars with large radii and rotation rates $>300$ Hz. We find that a 25% variation of the temperature of the hot spot with latitude does not significantly bias estimates of the mass $M$ and equatorial radius $R_{\rm eq}$ derived by fitting a model that assumes a uniform-temperature spot. Our results show that fits of oblate-star waveform models to waveform data can simultaneously determine $M$ and $R_{\rm eq}$ with uncertainties $\lesssim\,$7% if (1) the star's rotation rate is $\gtrsim\,$$600$ Hz; (2) the spot center and observer's sightline are both within $30^\circ$ of the star's rotational equator; (3) the oscillations have a fractional rms amplitude $\gtrsim\,$10%; and (4) $\gtrsim$$10^7$ counts are collected from the star. This is a realistic fractional amplitude, and this many counts could be obtained from a single star by the accepted NICER and proposed LOFT and AXTAR space missions by combining data from many X-ray bursts. These uncertainties are small enough to improve substantially our understanding of cold, ultradense matter.

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.

Simultaneous, precise measurements of the mass $M$ and radius $R$ of neutron stars can yield uniquely valuable information about the still uncertain properties of cold matter at several times the density of nuclear matter. One method that could be used to measure $M$ and $R$ is to analyze the energy-dependent waveforms of the X-ray flux oscillations seen during some thermonuclear bursts from some neutron stars. These oscillations are thought to be produced by X-ray emission from hotter regions on the surface of the star that are rotating at or near the spin frequency of the star. Here we explore how well $M$ and $R$ could be determined by generating, and analyzing using Bayesian techniques, synthetic energy-resolved X-ray data that we produce assuming a future space mission having 2--30 keV energy coverage and an effective area of 10 m$^2$, such as the proposed \textit{LOFT} or \textit{AXTAR} missions. We find that if the hot spot is within 10$^\circ$ of the rotation equator, both $M$ and $R$ can usually be determined with an uncertainty of about 10% if there are $10^6$ total counts from the spot, whereas waveforms from spots within 20$^\circ$ of the rotation pole provide no useful constraints. These constraints can usually be achieved even if the burst oscillations vary with time and data from multiple bursts must be used to obtain 10$^6$ counts from the hot spot. This is therefore a promising method to constrain $M$ and $R$ tightly enough to discriminate strongly between competing models of cold, high-density matter.

Precise and accurate measurements of neutron star masses and radii would provide valuable information about the still uncertain properties of cold matter at supranuclear densities. One promising approach to making such measurements involves analysis of the X-ray flux oscillations often seen during thermonuclear (type 1) X-ray bursts. These oscillations are almost certainly produced by emission from hotter regions on the stellar surface modulated by the rotation of the star. One consequence of the rotation is that the oscillation should appear earlier at higher photon energies than at lower energies. Ford (1999) found compelling evidence for such a hard lead in the tail oscillations of one type 1 burst from Aql X-1. We have therefore analyzed individually the oscillations observed in the tails of the four type 1 bursts from 4U 1636-536 that, when averaged, provided the strongest evidence for a soft lead in the analysis by Muno et al. (2003). We have also analyzed the oscillation observed during the superburst from this star. We find that the data from these five bursts, treated both individually and jointly, are fully consistent with a rotating hot spot model. Unfortunately, the uncertainties in these data are too large to provide interesting constraints on the mass and radius of this star.

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 show results of global 3D MHD simulations of accretion onto stars with superposition of the dipole and octupole fields, where we vary the ratio between components. Simulations show that if octupolar field strongly dominates at the disc-magnetosphere boundary, then matter flows into the ring-like octupolar poles, forming ring-shape spots at the surface of the star above and below equator. The light-curves are complex and may have two peaks per period. In case where the dipole field dominates, matter accretes in two ordered funnel streams towards poles, however the polar spots are meridionally-elongated due to the action of the octupolar component. In the case when the fields are of similar strengths, both, polar and belt-like spots are present. In many cases the light-curves show the evidence of complex fields, excluding the cases of small inclinations angles, where sinusoidal light-curve %is observed and `hides' the information about the field complexity. We also propose new mechanisms of phase shift in stars with complex magnetic fields. We suggest that the phase shifts can be connected with: (1) temporal variation of the star's intrinsic magnetic field and subsequent redistribution of main magnetic poles; (2) variation of the accretion rate, which causes the disc to interact with the magnetic fields associated with different magnetic moments. We use our model to demonstrate these phase shift mechanisms, and we discuss possible applications of these mechanisms to accreting millisecond pulsars and young stars.

X-ray burst spectra have long been used to estimate neutron star masses and radii. These estimates assumed that burst spectra are accurately described by the model atmosphere spectra developed over the last three decades. We compared RXTE data from a superburst with these spectra and found that the spectra predicted by previously published model atmospheres are strongly inconsistent with these high-precision measurements. In contrast, a simple Bose-Einstein spectrum is fully consistent with the data, as are recently published model atmosphere spectra. We discuss the implications of our results for determinations of neutron star masses and radii via constraints on their surface gravity and redshift, as originally suggested by Majczyna and Madej.

The magnetic field of the classical T Tauri star V2129 Oph can be modeled approximately by superposing slightly tilted dipole and octupole moments, with polar magnetic field strengths of 0.35kG and 1.2kG respectively (Donati et al. 2007). Here we construct a numerical model of V2129 Oph incorporating this result and simulate accretion onto the star. Simulations show that the disk is truncated by the dipole component and matter flows towards the star in two funnel streams. Closer to the star, the flow is redirected by the octupolar component, with some of the matter flowing towards the high-latitude poles, and the rest into the octupolar belts. The shape and position of the spots differ from those in a pure dipole case, where crescent-shaped spots are observed at the intermediate latitudes. Simulations show that if the disk is truncated at the distance of 6.2 R_* which is comparable with the co-rotation radius, 6.8 R_*, then the high-latitude polar spots dominate, but the accretion rate obtained from the simulations is about an order of magnitude lower than the observed one. The accretion rate matches the observed one if the disk is disrupted much closer to the star, at 3.4 R_*. However, the octupolar belt spots strongly dominate. Better match has been obtained in experiments with a dipole field twice as strong. The torque on the star from the disk-magnetosphere interaction is small, and the time-scale of spin evolution, 2 x10^7-10^9 years is longer than the 2x10^6 years age of V2129 Oph. The external magnetic flux of the star is strongly influenced by the disk: the field lines connecting the disk and the star inflate and form magnetic towers above and below the disk. The potential (vacuum) approximation is still valid inside the Alfv\'en (magnetospheric) surface where the magnetic stress dominates over the matter stress.

It has been known for nearly three decades that the energy spectra of thermonuclear X-ray bursts are often well-fit by Planck functions with temperatures so high that they imply a super-Eddington radiative flux at the emitting surface, even during portions of bursts when there is no evidence of photospheric radius expansion. This apparent inconsistency is usually set aside by assuming that the flux is actually sub-Eddington and that the fitted temperature is so high because the spectrum has been distorted by the energy-dependent opacity of the atmosphere. Here we show that the spectra predicted by currently available conventional atmosphere models appear incompatible with the highest-precision measurements of burst spectra made using the Rossi X-ray Timing Explorer, such as during the 4U 1820-30 superburst and a long burst from GX 17+2. In contrast, these measurements are well-fit by Bose-Einstein spectra with high temperatures and modest chemical potentials. Such spectra are very similar to Planck spectra. They imply surface radiative fluxes more than a factor of three larger than the Eddington flux. We find that segments of many other bursts from many sources are well-fit by similar Bose-Einstein spectra, suggesting that the radiative flux at the emitting surface also exceeds the Eddington flux during these segments. We suggest that burst spectra can closely approximate Bose-Einstein spectra and have fluxes that exceed the Eddington flux because they are formed by Comptonization in an extended, low-density radiating gas supported by the outward radiation force and confined by a tangled magnetic field.

We investigate further a model of the accreting millisecond X-ray pulsars we proposed earlier. In this model, the X-ray-emitting regions of these pulsars are near their spin axes but move. This is to be expected if the magnetic poles of these stars are close to their spin axes, so that accreting gas is channeled there. As the accretion rate and the structure of the inner disk vary, gas is channeled along different field lines to different locations on the stellar surface, causing the X-ray-emitting areas to move. We show that this "nearly aligned moving spot model" can explain many properties of the accreting millisecond X-ray pulsars, including their generally low oscillation amplitudes and nearly sinusoidal waveforms; the variability of their pulse amplitudes, shapes, and phases; the correlations in this variability; and the similarity of the accretion- and nuclear-powered pulse shapes and phases in some. It may also explain why accretion-powered millisecond pulsars are difficult to detect, why some are intermittent, and why all detected so far are transients. This model can be tested by comparing with observations the waveform changes it predicts, including the changes with accretion rate.

We have shown previously that many of the properties of persistent accretion-powered millisecond pulsars can be understood if their X-ray emitting areas are near their spin axes and move as the accretion rate and structure of the inner disk vary. Here we show that this "nearly aligned moving spot model" may also explain the intermittent accretion-powered pulsations that have been detected in three weakly magnetic accreting neutron stars. We show that movement of the emitting area from very close to the spin axis to about 10 degrees away can increase the fractional rms amplitude from less than about 0.5 percent, which is usually undetectable with current instruments, to a few percent, which is easily detectable. The second harmonic of the spin frequency usually would not be detected, in agreement with observations. The model produces intermittently detectable oscillations for a range of emitting area sizes and beaming patterns, stellar masses and radii, and viewing directions. Intermittent oscillations are more likely in stars that are more compact. In addition to explaining the sudden appearance of accretion-powered millisecond oscillations in some neutron stars with millisecond spin periods, the model explains why accretion-powered millisecond oscillations are relatively rare and predicts that the persistent accretion-powered millisecond oscillations of other stars may become undetectable for brief intervals. It suggests why millisecond oscillations are frequently detected during the X-ray bursts of some neutron stars but not others and suggests mechanisms that could explain the occasional temporal association of intermittent accretion-powered oscillations with thermonuclear X-ray bursts.

Detection of paired kilohertz quasi-periodic oscillations (kHz QPOs) in the X-ray emission of a compact object is compelling evidence that the object is an accreting neutron star. In many neutron stars, the stellar spin rate is equal or roughly equal to Delta-nu, the frequency separation of the QPO pair, or to 2Delta-nu. Hence, if the mechanism that produces the kilohertz QPOs is similar in all stars, measurement of Delta-nu can provide an estimate of the star's spin rate. The involvement of the stellar spin in producing Delta-nu indicates that the magnetic fields of these stars are dynamically important. We focus here on the implications of the paired kHz QPOs recently discovered in the low-mass X-ray binary (LMXB) system Cir X-1 (Boutloukos et al. 2006). The kHz QPOs discovered in Cir X-1 are generally similar to those seen in other stars, establishing that the compact object in the Cir X-1 system is a neutron star. However, the frequency nu-u of its upper kHz QPO is up to a factor of three smaller than is typical, and Delta-nu varies by about a factor 2 (167 Hz, the largest variation so far observed). Periodic oscillations have not yet been detected from Cir X-1, so its spin rate has not yet been measured directly. The low values of nu-u and the large variation of Delta-nu challenge current models of the generation of kHz QPOs. Improving our understanding of Cir X-1 will improve our knowledge of the spin rates and magnetic fields of all neutron stars.

Using the Rossi X-ray Timing Explorer (RossiXTE), astronomers have discovered that disk-accreting neutron stars with weak magnetic fields produce three distinct types of high-frequency X-ray oscillations. These oscillations are powered by release of the binding energy of matter falling into the strong gravitational field of the star or by the sudden nuclear burning of matter that has accumulated in the outermost layers of the star. The frequencies of the oscillations reflect the orbital frequencies of gas deep in the gravitational field of the star and/or the spin frequency of the star. These oscillations can therefore be used to explore fundamental physics, such as strong-field gravity and the properties of matter under extreme conditions, and important astrophysical questions, such as the formation and evolution of millisecond pulsars. Observations using RossiXTE have shown that some two dozen neutron stars in low-mass X-ray binary systems have the spin rates and magnetic fields required to become millisecond radio-emitting pulsars when accretion ceases, but that few have spin rates above about 600 Hz. The properties of these stars show that the paucity of spin rates greater than 600 Hz is due in part to the magnetic braking component of the accretion torque and to the limited amount of angular momentum that can be accreted in such systems. Further study will show whether braking by gravitational radiation is also a factor. Analysis of the kilohertz oscillations has provided the first evidence for the existence of the innermost stable circular orbit around dense relativistic stars that is predicted by strong-field general relativity. It has also greatly narrowed the possible descriptions of ultradense matter.

We have analyzed observations of the black hole transient {GX 339$-$4} made with the {\it Rossi} X-ray Timing Explorer (RXTE) and the Burst and Transient Source Experiment (BATSE) on board the Compton Gamma-ray Observatory (CGRO). We have found a nearly linear relation between the peak flux during the low/hard (LH) state that occurs at the beginning of an outburst and the time since the flux peak of the latest LH state identified in the previous outburst. Assuming that the rate at which mass accumulates in the accretion disk between these peaks is constant and that any mass that remains in the disk after an outburst has a negligible effect on the next outburst, this nearly linear relation suggests that the peak flux during the LH state that occurs at the beginning of an outburst is related to the mass in the disk, and thus that the entire disk is probably involved in powering these LH states. We have also found a positive correlation between the peak luminosities of the LH state in the three recent outbursts of {GX 339$-$4} and the peak luminosities of the following HS state. This correlation is similar to the correlations reported previously for {Aql X$-$1}, {4U 1705$-$44}, and {XTE J1550$-$564}, providing further support that the accretion flow that powers the LH state is related to the accretion flow that powers the following HS state. Although the luminosity of the LH-to-HS transition varies by up to an order of magnitude, the neutron stars and the black holes are distinguishable in the state transition luminosity. We discuss the implications for the mass determination of the compact stars in the ultraluminous X-ray sources (ULXs).

Quasi-periodic oscillations (QPOs) have been discovered in the X-ray emission of many neutron stars and black holes. The QPOs with frequencies greater than about 300 Hz are thought to be produced near the surfaces of neutron stars and the event horizons of black holes. I first summarize some of the most important properties of the QPOs seen in neutron star and black hole systems. I then review some of the models that have been proposed and compare them with observational data. Finally, I describe how these QPOs can be used to determine the properties of dense matter and strong gravitational fields.

Observations using the Rossi X-ray Timing Explorer have discovered dozens of accreting neutron stars with millisecond spin periods in low-mass binary star systems. Eighteen are millisecond X-ray pulsars powered by accretion or nuclear burning or both. These stars have magnetic fields strong enough for them to become millisecond rotation-powered (radio) pulsars when accretion ceases. Few, if any, accretion- or rotation-powered pulsars have spin rates higher than 750 Hz. There is strong evidence that the spin-up of some accreting neutron stars is limited by magnetic spin equilibrium whereas the spin-up of others is halted when accretion ends. Further study will show whether the spin rates of some accretion- or rotation-powered pulsars are or were limited by emission of gravitational radiation.

Motivated by the report by Cottam et al. (2002) of iron resonance scattering lines in the spectra of thermonuclear bursts from EXO 0748-676, we have investigated the information about neutron star structure and the geometry of the emission region that can be obtained by analyzing the profiles of atomic lines formed at the surface of the star. We have calculated the detailed profiles of such lines, taking into account the star's spin and the full effects of special and general relativity, including light-bending and frame-dragging. We discuss the line shapes produced by rotational Doppler broadening and magnetic splitting of atomic lines, for the spin rates and magnetic fields expected in neutron stars in low-mass X-ray binary systems. We show that narrow lines are possible even for rapidly spinning stars, if the emission region or the line of sight are close to the spin axis. For most neutron stars in low-mass systems, magnetic splitting is too small to obscure the effects of special and general relativity. We show that the ratio of the star's mass to its equatorial radius can be determined to within 5% using atomic line profiles, even if the lines are broad and skewed. This is the precision required to constrain strongly the equation of state of neutron star matter. We show further that if the radius and latitude of emission are known to ~ 5-10% accuracy, then frame-dragging has a potentially detectable effect on the profiles of atomic lines formed at the stellar surface.

We discuss millisecond period brightness oscillations and surface atomic spectral lines observed during type I X-ray bursts from a neutron star in a low mass X-ray binary system. We show that modeling of these phenomena can constrain models of the dense cold matter at the cores of neutron stars. We demonstrate that, even for a broad and asymmetric spectral line, the stellar radius-to-mass ratio can be inferred to better than 5%. We also fit our theoretical models to the burst oscillation data of the low mass X-ray binary XTE J1814-338, and find that the 90% confidence lower limit of the neutron star's dimensionless radius-to-mass ratio is 4.2.

This short encyclopedia article, reviewing current information on neutron stars, is intended for a broad scientific audience.

Observations made using the Rossi X-ray Timing Explorer have shown that accreting weak-field neutron stars in low-mass X-ray binary systems (LMXBs) produce three distinct types of millisecond X-ray oscillations that can be used to determine the spin rates and estimate the magnetic fields of these stars. These oscillations show that more than two dozen neutron stars in LMXBs have spin rates and magnetic fields in the range that will make them radio-emitting millisecond pulsars when accretion ceases, supporting the hypothesis that neutron stars in LMXBs are the progenitors of the rotation-powered millisecond pulsars. The spins of the 16 known accreting millisecond pulsars in LMXBs are consistent with spin-up to accretion spin equilibrium for magnetic fields ranging from 3 x 10^7 G to 3 x 10^8 G and time-averaged accretion rates ranging from 3 x 10^{-3} MdotE to MdotE, provided these stars have been accreting long enough. The P-Pdot distribution of rotation-powered millisecond pulsars indicates that their initial spins are set by spin-up to spin equilibrium at accretion rates ranging from 10^{-3} MdotE to MdotE or that they never reached accretion spin equilibrium.

KHz QPOs have now been detected in more than twenty accreting neutron stars in low-mass binary systems. Two kHz QPOs are usually detected in each star. Burst oscillations and two kHz QPOs have recently been detected in the 401 Hz accretion-powered X-ray pulsar SAX J1808.4-3658. In this star the frequency of the burst oscillation is approximately equal to the star's spin frequency nu_spin whereas the frequency separation of the two kHz QPOs is approximately nu_spin/2. If as expected the frequency of the burst oscillations in other stars is also approximately nu_spin, the frequency separation is approximately nu_spin in some stars but approximately nu_spin/2 in others. A frequency separation approximately equal to nu_spin/2 is unexplained in all existing models of the kHz QPOs. Here we propose a modified version of the sonic-point beat-frequency model that can explain within a single framework why the frequency separation is close to nu_spin in some stars but close to nu_spin/2 in others. As in the original sonic-point model, the frequency nu_QPO2 of the upper kHz QPO is close to the orbital frequency nu_orb at the radius r_{sp} of the sonic point in the disk flow. We show that magnetic and radiation fields rotating with the star will preferentially excite vertical motions in the disk at the "spin-resonance'' radius r_{sr} where nu_orb - nu_spin is equal to the vertical epicyclic frequency, producing vertical motions in the disk that modulate the X-ray flux at approximately nu_QPO2 - nu_spin or approximately nu_QPO2 - nu_spin/2, depending on whether the disk flow at r_{sr} is smooth or clumped. This sonic-point and spin-resonance model can also explain quantitatively the decrease of the kHz QPO frequency separation with increasing accretion rate that is observed in many sources.

The instability of r-mode oscillations in rapidly rotating neutron stars has attracted attention as a potential mechanism for producing high frequency, almost periodic gravitational waves. The analyses carried so far have shown the existence of these modes and have considered damping by shear and bulk viscosity. However, the magnetohydrodynamic coupling of the modes with a stellar magnetic field and its role in the damping of the instability has not been fully investigated yet. Following our introductory paper (Rezzolla, Lamb and Shapiro 2000), we here discuss in more detail the existence of secular higher-order kinematical effects which will produce toroidal fluid drifts. We also define the sets of equations that account for the time evolution of the magnetic fields produced by these secular velocity fields and show that the magnetic fields produced can reach equipartition in less than a year. The full numerical calculations as well as the evaluation of the impact of strong magnetic fields on the onset and evolution of the r-mode instability will be presented in a companion paper.

The evolution of the r-mode instability is likely to be accompanied by secular kinematic effects which will produce differential rotation with large scale drifts of fluid elements, mostly in the azimuthal direction. As first discussed by Rezzolla, Lamb and Shapiro 2000, the interaction of these secular velocity fields with a pre-existing neutron star magnetic field could result in the generation of intense and large scale toroidal fields. Following their derivation in the companion paper, we here discuss the numerical solution of the evolution equations for the magnetic field. The values of the magnetic fields obtained in this way are used to estimate the conditions under which the r-mode instability might be prevented or suppressed. We also assess the impact of the generation of large magnetic fields on the gravitational wave detectability of r-mode unstable neutron stars. Our results indicate that the signal to noise ratio in the detection of gravitational waves from the r-mode instability might be considerably decreased if the latter develops in neutron stars with initial magnetic fields larger than 10^10 G.

We investigate the suggestion that the frequencies of the kilohertz and other QPOs in the emission from neutron stars in low-mass X-ray binaries are generated by geodesic motion of gas clumps around the star. First we assume, following previous work, that the dominant frequencies produced by such motions are the azimuthal frequency \nuk, the apsidal precession frequency \nuap, and the first overtone of the nodal precession frequency \nunp. We investigate whether geodesics can be found for which these frequencies agree with the observed frequencies. Correcting calculational errors made when the geodesic precession hypothesis was first proposed, we find that \nuk and \nuap disagree qualitatively with the frequencies of the kilohertz QPOs for infinitesimally or even moderately eccentric geodesics. These frequencies are similar to the frequencies of the kilohertz QPOs only for highly eccentric geodesics, with apastron to periastron ratios \sim3--4; for these geodesics, 2\nunp differs qualitatively from the frequencies of the low-frequency QPOs. Next we investigate whether these frequencies would be the dominant frequencies produced by orbiting clumps. We find that they are not and that the dominant frequencies are instead harmonics and sidebands of the radial epicyclic frequency \nur. Finally, we show that gas dynamical constraints restrict the radii of orbiting clumps to be \ll 10^{-3} of the stellar radius; the fractional modulation that could be produced by clumps is therefore \ll 10^{-6}. We conclude that there are significant difficulties in interpreting the kilohertz and low-frequency QPOs as consequences of geodesic motion of gas clumps.

Previous work on the sonic-point beat-frequency (SPBF) model of the kilohertz quasi-periodic oscillations (QPOs) observed in the X-ray flux from neutron stars in low-mass binary systems has shown that it naturally explains many properties of these QPOs. These include the existence of just two principal QPOs in a given source, the commensurability of the frequency separation \Dnu\ of the two kilohertz QPOs and the spin frequency \nus inferred from burst oscillations, and the high frequencies, coherence, and amplitudes of these QPOs. Here we show that the SPBF model predicts that \Dnu is less than but close to \nus, consistent with the observed differences between \Dnu and \nus. It also explains naturally the decrease in \Dnu with increasing QPO frequency seen in some sources and the plateau in the QPO frequency--X-ray flux observed in 4U 1820-30. The model fits well the QPO frequency behavior observed in Sco X-1, 4U 1608-52, 4U 1728-34, and 4U 1820-30 (chisqdof = 0.4-2.1), giving masses ranging from 1.59 to 2.0 Msun and spin rates ranging from 279 to 364 Hz. In the SPBF model, the kilohertz QPOs are effects of strong-field gravity. Thus, if the model is validated, the kilohertz QPOs can be used not only to determine the properties of neutron stars but also to explore quantitatively general relativistic effects in the strong-field regime.

We present an analysis of ~390 ksec of data of the Z source GX 340+0 taken during 24 observations with the Rossi X-ray Timing Explorer satellite. We report the discovery of a new broad component in the power spectra. The frequency of this component varied between 9 and 14 Hz, and remained close to half that of the horizontal branch quasi-periodic oscillations (HBO). Its rms amplitude was consistent with being constant around ~5%, while its FWHM increased with frequency from 7 to 18 Hz. If this sub-HBO component is the fundamental frequency, then the HBO and its second harmonic are the second and fourth harmonic component, while the third harmonic was not detected. This is similar to what was recently found for the black hole candidate XTE J1550-564. The profiles of both the horizontal- and the normal branch quasi-periodic oscillation peaks were asymmetric when they were strongest. We describe this in terms of a shoulder component at the high frequency side of the quasi-periodic oscillation peak, whose rms amplitudes were approximately constant at ~4% and \~3%, respectively. The peak separation between the twin kHz quasi-periodic oscillations was consistent with being constant at 339+-8 Hz but a trend similar to that seen in, e.g. Sco X-1 could not be excluded. We discuss our results within the framework of the various models which have been proposed for the kHz QPOs and low frequency peaks.

We show that r-mode oscillations distort the magnetic fields of neutron stars and that their occurrence is likely to be limited by this interaction. If the field is gtrsim 10^{16} (Omega/Omega_B) G, where Omega and Omega_B are the angular velocities of the star and at which mass shedding occurs, r-mode oscillations cannot occur. Much weaker fields will prevent gravitational radiation from exciting r-mode oscillations or damp them on a relatively short timescale by extracting energy from the modes faster than gravitational wave emission can pump energy into them. For example, a 10^{10} G poloidal magnetic field that threads the star's superconducting core is likely to prevent the ell=2 mode from being excited unless Omega exceeds 0.35 Omega_B. If Omega is larger than 0.35 Omega_B initially, the ell=2 mode may be excited but is likely to decay rapidly once Omega falls below 0.35 Omega_B, which happens in lesssim 15^d if the saturation amplitude is gtrsim 0.1. The r-mode oscillations may play an important role in determining the structure of neutron star magnetic fields.

We have measured the phase-delay and rms amplitude spectra of CygX-2 and ScoX-1. Using EXOSAT data from the normal branch of CygX-2 we confirm earlier (Ginga) results, showing that at an energy near 6keV there is both a minimum in the QPO rms amplitude spectrum and a 150 degree phase jump in the quasi-periodic oscillations (QPO) phase-delay spectrum. Surprisingly, using EXOSAT and Ginga data, we find no evidence for a phase jump of this kind in the phase-delay spectrum of ScoX-1 on either the normal or flaring branch. Upper limits (90% confidence) of 42 degrees can be set on any phase jump in the energy range 2-10keV on the normal branch, and 88 degrees on the flaring branch. The QPO rms amplitude spectrum of ScoX-1 increases steeply with energy on both the normal and flaring branches. These results suggest that the X-ray spectrum pivots about an energy of <~2keV or >~10keV or that normal branch QPO of ScoX-1 does not have a pivoting spectrum. We discuss the implications of these results in terms of the radiation-hydrodynamic model for normal branch QPO.

The great success of the Rossi X-Ray Timing Explorer (RXTE) has given us a new probe to study strong gravitational fields and to measure the physical properties of black holes and neutron stars. Here, we describe a "next-generation" x-ray timing mission, the Relativistic Astrophysics Explorer (RAE), designed to fit within the envelope of a "medium-sized" mission. The main instruments will be a narrow-field x-ray detector array with an area of at least ten times that of RXTE, and a wide-field x-ray monitor with good sensitivity and few arcminute position resolution. We describe the design of the instruments and the science which will be possible with a factor of ten increase in collecting area.

We study the horizontal branch oscillations (HBO) and the two simultaneous kilohertz QPOs comparing their properties in five Z sources with the predictions of the magnetospheric beat-frequency and Lense-Thirring precession models. We find that the variation of the HBO frequency with accretion rate predicted by the magnetospheric beat-frequency model for a purely dipolar stellar magnetic field and a radiation-pressure-dominated inner accretion disk agrees well with the observed variation. This model implies that the neutron stars in the Z sources are near magnetic spin equilibrium, that their magnetic field strengths are ~10^9-10^10 G, and that the critical fastness parameter for these sources is >0.8. If the frequency of the upper kilohertz QPO is an orbital frequency in the accretion disk, the magnetospheric beat-frequency model requires that a small fraction of the gas in the disk does not couple strongly to the stellar magnetic field at 3-4 stellar radii but instead drifts slowly inward in nearly circular orbits until it is within a few kilometers of the neutron star surface. The Lense-Thirring precession model is consistent with the observed magnitudes of the HBO frequencies only if the moments of inertia of the neutron stars in the Z sources are 4-5 times larger than the largest values predicted by realistic neutron-star equations of state. The trend of the correlation between the HBO and upper kHz QPO frequencies is not consistent with this model. We argue that the change in the slope of the correlation between the frequency of the horizontal branch oscillation and the frequency of the upper kilohertz QPO, when the latter is greater than 850 Hz, is directly related to the varying frequency separation of the kilohertz QPOs.

Strong, coherent, quasi-periodic brightness oscillations (QPOs) with frequencies ranging from about 300 Hz to 1200 Hz have been discovered with the Rossi X-ray Timing Explorer in the X-ray emission from some fifteen neutron stars in low-mass binary systems. Two simultaneous kilohertz QPOs differing in frequency by 250 to 350 Hertz have been detected in twelve of the fifteen sources. Here we propose a model for these QPOs. In this model the X-ray source is a neutron star with a surface magnetic field of 10^7 to 10^10 G and a spin frequency of a few hundred Hertz, accreting gas via a Keplerian disk. The frequency of the higher-frequency QPO in a kilohertz QPO pair is the Keplerian frequency at a radius near the sonic point at the inner edge of the Keplerian flow whereas the frequency of the lower-frequency QPO is approximately the difference between the Keplerian frequency at a radius near the sonic point and the stellar spin frequency. This model explains naturally many properties of the kilohertz QPOs, including their frequencies, amplitudes, and coherence. We show that if the frequency of the higher-frequency QPO in a pair is an orbital frequency, as in the sonic-point model, the frequencies of these QPOs place interesting upper bounds on the masses and radii of the neutron stars in the kilohertz QPO sources and provide new constraints on the equation of state of matter at high densities. Further observations of these QPOs may provide compelling evidence for the existence of a marginally stable orbit, confirming a key prediction of general relativity in the strong-field regime.

It has recently been suggested that gravitomagnetic precession of the inner part of the accretion disk, possibly driven by radiation torques, may be responsible for some of the quasi-periodic X-ray brightness oscillations (QPOs) at 20 - 300 Hz, observed in some low-mass binary systems containing accreting neutron stars and black hole candidates. We have explored warp modes of geometrically thin disks in the presence of gravitomagnetic and radiation torques. We have found a family of overdamped, low-frequency gravitomagnetic (LFGM) modes with precession frequencies lower than a certain critical frequency \omega_{crit}, which is ~ 1 Hz for a compact object of solar mass. A radiation warping torque can cause a few of the lowest-frequency LFGM modes to grow with time, but even a strong radiation warping torque has essentially no effect on the LFGM modes with frequencies > 10^{-4} Hz. We have also discovered a second family of high-frequency gravitomagnetic (HFGM) modes with precession frequencies that range from $\omega_{crit}$ up to slightly less than the gravitomagnetic precession frequency $\omega_{gm,i}$ of a particle at the inner edge of the disk, which is 30 Hz if the disk extends inward to the innermost stable circular orbit around a compact object of 2 solar masses with dimensionless angular momentum cJ/GM^2 = 0.2. The highest-frequency HFGM modes are very localized spiral corrugations of the inner disk and are weakly damped, with Q values 2-50. We discuss the implications of our results for the observability of Lense-Thirring precession in X-ray binaries.

We discovered two simultaneous kHz quasi-periodic oscillations (QPOs) in the bright low-mass X-ray binary and Z source GX 5-1 with the Rossi X-ray Timing Explorer. In the X-ray color-color and hardness-intensity diagram a clear Z track is traced out, which shifted between observations. The frequencies of the two kHz QPOs increased from ~215 Hz and ~500 Hz on the left part of the horizontal branch to ~700 Hz and ~890 Hz, respectively, on the upper part of the normal branch. With increasing frequency the FWHM and rms amplitude (8.6-60 keV) of the higher-frequency kHz QPO decreased from 300 to 30 Hz, and from 6.6 to 2.4%, respectively. The FWHM and amplitude of the lower-frequency kHz QPO (50-100 Hz and 3-4% rms) did not correlate with the position of the source on the Z track. The kHz QPO separation was consistent with being constant at 298+/-11 Hz. Simultaneously with the kHz QPOs horizontal branch oscillations (HBOs) were detected with frequencies between 18 and 56 Hz.

Pairs of quasi-periodic oscillations (QPOs) at kilohertz frequencies are a common phenomenon in several neutron-star low-mass X-ray binaries. The frequency separation of the QPO peaks in the pair appears to be constant in many sources and directly related to the neutron star spin frequency. However, in Sco X-1 and possibly in 4U 1608-52, the frequency separation of the QPOs decreases with increasing inferred mass accretion rate. We show that the currently available Rossi X-ray Timing Explorer data are consistent with the hypothesis that the frequency separations in all sources vary by amounts similar to the variation in Sco X-1. We discuss the implications for models of the kilohertz QPOs.

Quasi-periodic X-ray brightness oscillations (QPOs) with frequencies around a kilohertz have now been discovered in more than a dozen neutron stars in low-mass X-ray binary systems using the Rossi X-ray Timing Explorer. There is strong evidence that the frequencies of the kilohertz oscillations are the orbital frequencies of accreting gas in nearly circular orbits around these stars. Some stars that produce kilohertz QPOs may have spin frequencies greater than 400 Hertz. For spin rates this high, first-order analytic treatments of the effects of the star's rotation on its structure and the spacetime are inaccurate. Here we use the results of a large number of fully relativistic, self-consistent numerical calculations of the stellar structure of rapidly rotating neutron stars and the interior and exterior spacetime to investigate the constraints on the properties of such stars that can be derived if stable circular orbits of various frequencies are observed. We have computed the equatorial radius of the star, the radius of the innermost stable circular orbit, and the frequency of the highest-frequency stable circular orbit as functions of the stellar spin rate, for spin rates up to the maximum possible and for several illustrative equations of state. Our calculations show that the upper bounds on the stiffness of neutron star matter implied by a given orbital frequency are typically significantly stricter for stars with spin frequencies more than 400 Hertz than for slowly rotating stars.

We have discovered two simultaneous kHz quasi-periodic oscillations (QPOs) in the Z source GX 340+0 with the Rossi X-ray Timing Explorer. The X-ray hardness-intensity and color-color diagram each show a full Z-track, with an extra limb branching off the flaring branch of the Z. Both peaks moved to higher frequencies when the mass accretion rate increased. The two peaks moved from 247 +/- 6 and 567 +/- 39 Hz at the left end of the horizontal branch to 625 +/- 18 and 820 +/- 19 Hz at its right end. The higher frequency peak's rms amplitude (5-60 keV) and FWHM decreased from ~5% and 383 +/- 135 Hz to ~2%, and 145 +/- 62 Hz, respectively. The rms amplitude and FWHM of the lower peak were consistent with being constant near 2.5 % and 100 Hz. The kHz QPO separation was consistent with being constant at 325 +/- 10 Hz. Simultaneous with the kHz QPOs we detected the horizontal branch oscillations (HBO) and its second harmonic, at frequencies between 20 and 50 Hz, and 38 and 69 Hz, respectively. The normal branch oscillations were only detected on the upper and middle normal branch, and became undetectable on the lower normal branch. The HBO frequencies do not fall within the range predicted for Lense-Thirring (LT) precession, unless either the ratio of the neutron star moment of inertia to neutron star mass is at least 4, 10^45 gcm^2/M_sun, the frequencies of the HBO are in fact the sub-harmonic oscillations, or the observed kHz peak difference is half the spin frequency and not the spin frequency. During a 1.2 day gap between two observations, the Z-track in the hardness-intensity diagram moved to higher count rates by about 3.5%. Comparing data before and after this shift, we find that the HBO properties are determined by position on the Z-track and not directly by count rate or X-ray colors.

One of the most dramatic discoveries made with the Rossi X-Ray Timing Explorer is that many accreting neutron stars in low-mass binary systems produce strong, remarkably coherent, high-frequency X-ray brightness oscillations. The 325-1200 Hz quasi-periodic oscillations (QPOs) observed in the accretion-powered emission are thought to be produced by gas orbiting very close to the neutron star, whereas the 360-600 Hz brightness oscillations seen during thermonuclear X-ray bursts are produced by one or two hot spots rotating with the star and have frequencies equal to the stellar spin frequency or its first overtone. The oscillations constrain the masses and radii of these neutron stars, which are thought to be the progenitors of the millisecond pulsars. Modeling indicates that the stars have spin frequencies 250-350 Hz and magnetic fields 1e7-5e9 G.

One of the most dramatic discoveries made so far with the Rossi X-Ray Timing Explorer is that many accreting neutron stars with weak magnetic fields generate strong, remarkably coherent, high-frequency X-ray brightness oscillations. The 325-1200 Hz quasi-periodic oscillations (QPOs) observed in the accretion-powered emission are almost certainly produced by gas orbiting very close to the stellar surface and have frequencies related to the orbital frequencies of the gas. The 360-600 Hz brightness oscillations seen during thermonuclear X-ray bursts are produced by one or two hotter regions on the stellar surface and have frequencies equal to the stellar spin frequency or its first overtone. Measurements of these oscillations are providing tight upper bounds on the masses and radii of neutron stars, and important new constraints on the equation of state of neutron star matter.

The dramatic discovery with the Rossi X-Ray Timing Explorer satellite of remarkably coherent $\sim$300--1200 Hz oscillations in the X-ray brightness of some sixteen neutron stars in low-mass binary systems has spurred theoretical modeling of these oscillations and investigation of their implications for the neutron stars and accretion flows in these systems. High-frequency oscillations are observed both during thermonuclear X-ray bursts and during intervals of accretion-powered emission and appear to be a characteristic feature of disk-accreting neutron stars with weak magnetic fields. In this review we focus on the high-frequency quasi-periodic oscillations (QPOs) seen in the accretion-powered emission. We first summarize the key properties of these kilohertz QPOs and then describe briefly the models that have been proposed to explain them. The existing evidence strongly favors beat-frequency models. We mention several of the difficulties encountered in applying the magnetospheric beat-frequency model to the kilohertz QPOs. The most fully developed and successful model is the sonic-point beat-frequency model. We describe the work on this model in some detail. We then discuss observations that could help to distinguish between models. We conclude by noting some of the ways in which study of the kilohertz QPOs may advance our understanding of dense matter and strong gravitational fields.

The frequencies of the highest-frequency kilohertz QPOs recently discovered in some sixteen neutron stars in low-mass X-ray binary systems are most likely the orbital frequencies of gas in Keplerian orbit around these neutron stars. If so, these QPOs provide tight upper bounds on the masses and radii of these neutron stars and important new constraints on the equation of state of neutron star matter. If the frequency of a kilohertz QPO can be established as the orbital frequency of gas at the innermost stable circular orbit, this would confirm one of the key predictions of general relativity in the strong-gravity regime. If the spin frequency of the neutron star can also be determined, the frequency of the QPO would fix the mass of the neutron star for each assumed equation of state. Here we show how to derive mass and radius bounds, using the kilohertz QPOs, for nonrotating and slowly rotating stars, and discuss how these bounds are affected by rapid stellar rotation and radial radiation forces. We also describe observational results that would be strong evidence for the presence of an innermost stable circular orbit. No such strong evidence is present in current data, but future prospects are excellent.

We discovered a single kHz quasi-periodic oscillation (QPO) near 1150 Hz in the Rossi X-ray Timing Explorer X-ray light curve of the low mass X-ray binary and atoll source 4U 1735-44. The rms amplitude of this peak was 2-3%, and the FWHM 6-40 Hz. There are indications that the kHz QPO frequency decreased from 1160 Hz to 1145 Hz when the count rate increased, which would be quite different from what is observed in other atoll sources for which kHz QPOs have been discovered. In the X-ray color-color diagram and hardness-intensity diagram the source traced out the curved branch (the so-called banana branch) which has been found by previous instruments. The kHz QPO was only detected when the source was at the lowest count rates during our observations, i.e. on the lower part of the banana branch. When 4U 1735-44 was at higher count rates, i.e. on the upper part of the banana branch and at higher inferred mass accretion rate with respect to that on the lower part of the banana branch, the QPO was not detected. Besides the kHz QPO we discovered a low frequency QPO with a frequency near 67 Hz, together with a complex broad peaked noise component below 30 Hz. This 67 Hz QPO may be related to the magnetospheric beat-frequency QPO, which is observed on the horizontal branch of Z sources. This idea is supported by the (peaked) noise found in both 4U 1735-44 and Z sources at frequencies just below the QPO frequency.

The frequencies of the highest-frequency kilohertz QPOs recently discovered in some sixteen neutron stars in low-mass X-ray binary systems are most likely orbital frequencies. If so, these QPOs provide tight upper bounds on the masses and radii of these neutron stars and interesting new constraints on the equation of state of neutron star matter. If the frequency of a kilohertz QPO can be established as the orbital frequency of the innermost stable circular orbit, this would confirm one of the key predictions of general relativity in the strong-field regime. If the spin frequency of the neutron star can also be determined, the frequency of the QPO would fix the mass of the neutron star for each assumed equation of state. Here we describe how bounds on the stellar mass and radius can be derived and how these bounds are affected by the stellar spin. We also discuss detection of the innermost stable circular orbit.

Using a new technique to improve the sensitivity to weak Quasi-Periodic Oscillations (QPO) we discovered a new QPO peak at about 1100 Hz in the March 1996 outburst observations of 4U 1608-52, simultaneous with the ~ 600-900 Hz peak previously reported from these data. The frequency separation between the upper and the lower QPO peak varied significantly from $232.7 \pm 11.5$ Hz on March 3, to $293.1 \pm 6.6$ Hz on March 6. This is the first case of a variable kHz peak separation in an atoll source.

Using a new technique to improve the sensitivity to weak Quasi-Periodic Oscillations (QPO) we discovered a new QPO peak at about 1100 Hz in the March 1996 outburst observations of 4U 1608-52, simultaneous with the ~ 600-900 Hz peak previously reported from these data. The frequency separation between the upper and the lower QPO peak varied significantly from $232.7 \pm 11.5$ Hz on March 3, to $293.1 \pm 6.6$ Hz on March 6. This is the first case of a variable kHz peak separation in an atoll source. We discuss to what extent this result could be accommodated in beat-frequency models such as proposed for the kHz QPOs. We measured the rms fractional amplitude of both QPOs as a function of energy, and we found that the relation is steeper for the lower than for the upper frequency peak. This is the first source where such difference between the energy spectrum of the two kHz QPOs could be measured.

The discovery of high-amplitude brightness oscillations at the spin frequency or its first overtone in six neutron stars in low-mass X-ray binaries during type~1 X-ray bursts provides a powerful new way to constrain the compactness of these stars, and hence to constrain the equation of state of the dense matter in all neutron stars. Here we present the results of general relativistic calculations of the maximum fractional rms amplitudes that can be observed during bursts. In particular, we determine the dependence of the amplitude on the compactness of the star, the angular dependence of the emission from the surface, the rotational velocity at the stellar surface, and whether there are one or two emitting poles. We show that if two poles are emitting, as is strongly indicated by independent evidence in 4U 1636-536 and KS 1731-26, the resulting limits on the compactness of the star can be extremely restrictive. We also discuss the expected amplitudes of X-ray color oscillations and the observational signatures necessary to derive convincing constraints on neutron star compactness from the amplitudes of burst oscillations.

During observations with RXTE we discovered two simultaneous kHz QPOs near 500 Hz and 860 Hz in the Z source Cyg X-2. In the X-ray color-color diagram and hardness-intensity diagram (HID) a clear Z track was traced out. A single kHz QPO peak was detected at the left end of the horizontal branch (HB) of the Z track, with a frequency of 731+/-20 Hz and an amplitude of 4.7+0.8-0.6 % rms in the energy band 5.0-60 keV. Further to the right on the HB, at somewhat higher count rates, an additional peak at 532+/-43 Hz was detected with an rms amplitude of 3.0+1.0-0.7 %. When the source moved down the HB, thus when the inferred mass accretion rate increased, the frequency of the higher-frequency QPO increased to 839+/-13 Hz, and its amplitude decreased to 3.5+0.4-0.3 % rms. The higher-frequency QPO was also detected on the upper normal branch (NB) with an rms amplitude of 1.8+0.6-0.4 % and a frequency of 1007+/-15 Hz; its peak width did not show a clear correlation with inferred mass accretion rate. The lower-frequency QPO was most of the time undetectable, with typical upper limits of 2 % rms, no conclusion on how this QPO behaved with mass accretion rate can be drawn. If the peak separation between the QPOs is the neutron star spin frequency (as required in some beat-frequency models) then the neutron star spin period is 2.9+/-0.2 ms (346+/-29 Hz).

We observed the low-mass X-ray binary and Z source GX 17+2 with the Rossi X-ray Timing Explorer during 6-8 Feb 1997, 1-4 Apr 1997, and 26-27 Jul 1997. The X-ray color-color diagram shows a clear Z track. Two simultaneous kHz quasi-periodic oscillations (QPOs) are present in each observation, whose frequencies are well correlated with the position of the source on the Z-track. At the left end of the horizontal branch (HB) only the higher frequency peak is observed, at 645+/-9 Hz, with an rms amplitude of 5.7+/-0.5%, and a FWHM of 183+/-35 Hz. When the source moves down the Z track to the upper normal branch the frequency of the kHz QPO increases to 1087+/-12 Hz, and the rms amplitude and FWHM decrease by a factor 2. Further down the Z the QPO becomes undetectable, with rms upper limits of typically 2.0%. Halfway down the HB a second QPO appears in the power spectra with a frequency of 480+/-23 Hz. The frequency of this QPO also increases when the source moves along the Z track, up to 781+/-11 Hz halfway down the normal branch, while the rms amplitude and FWHM stay approximately constant at 2.5% and 70 Hz. The QPO frequency difference is constant at 293.5+/-7.5 Hz. Simultaneously with the kHz QPOs we detect HB QPOs (HBOs). The simultaneous presence of HBOs and kHz QPOs excludes the magnetospheric beat-frequency model as the explanation for at least one of these two phenomena.

Compton scattering of photons by nonrelativistic particles is thought to play an important role in forming the radiation spectrum of many astrophysical systems. Here we derive the time-dependent photon kinetic equation that describes spontaneous and induced Compton scattering as well as absorption and emission by static and moving media, the corresponding radiative transfer equation, and their zeroth and first moments, in both the system frame and in the frame comoving with the medium. We show that it is necessary to use the correct relativistic differential scattering cross section in order to obtain a photon kinetic equation that is correct to first order in epsilon/m_e, T_e/m_e, and V, where epsilon is the photon energy, T_e and m_e are the electron temperature and rest mass, and V is the electron bulk velocity in units of the speed of light. We also demonstrate that the terms in the radiative transfer equation that are second-order in V usually should be retained, because if the radiation energy density is sufficiently large compared to the radiation flux, the effects of bulk Comptonization described by the terms that are second-order in V are at least as important as the effects described by the terms that are first-order in V, even when V is small. Our equations are valid for systems of arbitrary optical depth and can therefore be used in both the free-streaming and the diffusion regimes. We demonstrate that Comptonization by the electron bulk motion occurs whether or not the radiation field is isotropic or the bulk flow converges and that it is more important than thermal Comptonization if V^2 > 3 T_e/m_e.

We report three RXTE/PCA observations of the low-mass X-ray binary 4U 0614+09. They show strong (~ 30% rms) band-limited noise with a cut-off frequency varying between 0.7 and 15 Hz in correlation with the X-ray flux. We observe two non-simultaneous 11-15% (rms) kHz peaks near 728 and 629 Hz in the power spectra of two of our observations when the X-ray flux is ~ 10^{-9} erg cm^{-2} s{-1} (2-10 keV), but find no QPO (< 6% rms) when the X-ray flux is half that. We suggest that count rate may not be a good measure for $\dot M$ even in sources as intrinsically weak as 4U 0614+09, and that QPO frequency and noise cutoff frequency track $\dot M$ more closely than count rate. The QPO increases in rms amplitude from ($11 \pm 1.3$) % at 3 keV to ($37 \pm 12$) % at 23 keV; the fractional amplitude of the band-limited noise is energy-independent. This suggests different sites of origin for these two phenomena. The spectrum of the oscillating flux roughly corresponds to a black body with temperature ($1.56 \pm 0.2$) keV and radius ($ 500 \pm 200$) m (other models fit as well), which might indicate the oscillations originate at a small region on the neutron star surface.

Using the Rossi X-ray Timing Explorer we have measured 27 microsecond time delays in 830 Hz quasi-periodic oscillations (QPO) between 4-6 keV and 11-17 keV in 4U 1608-52, with high-energy photons lagging low energy photons, and found upper limits to the time delays of 45 microseconds between 2-6.5 keV and 6.5-67 keV in 730 Hz QPO in 4U 0614+091 and 30 microseconds between 8.7-12.4 keV and 12.4-67 keV in 870 Hz QPO in 4U 1636-53. We also find that the cross coherence function between QPO at different energies is > 0.85 with 95% confidence in 4U 1608-52 and 4U 1636-53. If Compton upscattering of low-energy X-rays in a region with an optical depth of a few is responsible for the delays, then the Compton upscattering region is between a few kilometers and a few tens of kilometers in size.

Quasi-periodic brightness oscillations (QPOs) with frequencies in the range 300-1200 Hz have been detected from at least nine neutron star low-mass X-ray binaries. Here we summarize the sonic-point model for these brightness oscillations, which we present in detail elsewhere. If the sonic-point interpretation of kilohertz QPOs is confirmed, measurements of kilohertz QPO frequencies in low-mass X-ray binaries will provide new bounds on the masses and radii of neutron stars in these systems and new constraints on the equation of state of matter at high densities.

We present results of RXTE observations of the low-mass X-ray binary and atoll source 4U 1608$-$52 made over 9 days during the decline of an X-ray intensity outburst in March 1996. A fast-timing analysis shows a strong and narrow quasi periodic oscillation (QPO) peak at frequencies between 850 and 890 Hz on March 3 and 6, and a broad peak around 690 Hz on March 9. Observations on March 12 show no significant signal. On March 3, the X-ray spectrum of the QPO is quite hard; its strength increases steadily from 5 \% at $\sim$2 to $\sim$20 \% at $\sim$12~keV. The QPO frequency varies between 850 and 890 Hz on that day, and the peak widens and its rms decreases with centroid frequency in a way very similar to the well-known horizontal branch oscillations (HBO) in Z-sources. We apply the HBO beat frequency model to atoll sources, and suggest that, whereas the model could produce QPOs at the observed frequencies, the lack of correlation we observe between QPO properties and X-ray count rate is hard to reconcile with this model.

A simple, physically consistent model has been proposed that seeks to explain in a unified way the X-ray spectra and rapid X-ray variability of the so-called Z sources and other accreting neutron stars in low-mass systems. Here we summarize the results of detailed numerical calculations of the X-ray spectra of the Z sources predicted by this model. Our computations show that in the Z sources, photons are produced primarily by electron cyclotron emission in the neutron star magnetosphere. Comptonization of these photons by the hot central corona and radial inflow produces X-ray spectra, color-color tracks, and countrate variations like those observed in the Z sources.