We develop a new model for X-ray emission from tidal disruption events (TDEs), applying stationary general relativistic ``slim disk'' accretion solutions to supermassive black holes (SMBHs) and then ray-tracing the photon trajectories from the image plane to the disk surface, including gravitational redshift, Doppler, and lensing effects self-consistently. We simultaneously and successfully fit the multi-epoch XMM-Newton X-ray spectra for two TDEs: ASASSN-14li and ASASSN-15oi. We test explanations for the observed, unexpectedly slow X-ray brightening of ASASSN-15oi, including delayed disk formation and variable obscuration by a reprocessing layer. We propose a new mechanism that better fits the data: a ``Slimming Disk'' scenario in which accretion onto an edge-on disk slows, reducing the disk height and exposing more X-rays from the inner disk to the sightline over time.For ASASSN-15oi, we constrain the SMBH mass to $4.0^{+2.5}_{-3.1} \times 10^6M_\odot$. For ASASSN-14li, the SMBH mass is $10^{+1}_{-7}\times 10^6M_\odot$ and the spin is $>0.3$. For both TDEs, our fitted masses are consistent with independent estimates; for ASASSN-14li, application of the external mass constraint narrows our spin constraint to $>0.85$. The mass accretion rate of ASASSN-14li decays slowly, as $\propto t^{-1.1}$, perhaps due to inefficient debris circularization. Over $\approx$1100 days, its SMBH has accreted $\Delta M \approx 0.17 M_\odot$, implying a progenitor star mass of $> 0.34 M_\odot$, i.e., no ``missing energy problem.'' For both TDEs, the hydrogen column density declines to the host galaxy plus Milky Way value after a few hundred days, suggesting a characteristic timescale for the depletion or removal of obscuring gas.

The Galactic Center black hole Sgr A* shows significant variability and flares in the submillimeter, infrared, and X-ray wavelengths. Owing to its exquisite resolution in the IR bands, the GRAVITY experiment for the first time spatially resolved the locations of three flares and showed that a bright region moves in ellipse-like trajectories close to but offset from the black hole over the course of each event. We present a model for plasmoids that form during reconnection events and orbit in the coronal region around a black hole to explain these observations. We utilize general-relativistic radiative transfer calculations that include effects from finite light travel time, plasmoid motion, particle acceleration, and synchrotron cooling and obtain a rich structure in the flare lightcurves. This model can naturally account for the observed motion of the bright regions observed by the GRAVITY experiment and the offset between the center of the centroid motion and the position of the black hole. It also explains why some flares may be double-peaked while others have only a single peak and uncovers a correlation between the structure in the lightcurve and the location of the flare. Finally, we make predictions for future observations of flares from the inner accretion flow of Sgr A* that will provide a test of this model.

We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.

The Event Horizon Telescope (EHT) is taking the first images of black holes resolved at horizon scales to measure their shadows and probe accretion physics. A promising avenue for testing the hypothesis that astrophysical black holes are described by the Kerr solution to Einstein's equations is to compare the size and shape of the shadow a black hole casts on the surrounding emission to the predictions of the Kerr metric. We develop here an efficient parametric framework to perform this test. We carry out ray-tracing simulations for several parametrized non-Kerr metrics to create a large data set of non-Kerr shadows that probe the allowed parameter space for the free parameters of each metric. We then perform principal components analysis (PCA) on this set of shadows and show that only a small number of components are needed to accurately reconstruct all shadows within the set. We further show that the amplitude of the PCA components are smoothly related to the free parameters in the metrics, and therefore, that these PCA components can be fit to EHT observations in order to place constraints on the free parameters of these metrics that will help quantify any potential deviations from the Kerr solution.

Interferometers, such as the Event Horizon Telescope (EHT), do not directly observe the images of sources but rather measure their Fourier components at discrete spatial frequencies up to a maximum value set by the longest baseline in the array. Construction of images from the Fourier components or analysis of them with high-resolution models requires careful treatment of fine source structure nominally beyond the array resolution. The primary EHT targets, Sgr A* and M87, are expected to have black-hole shadows with sharp edges and strongly filamentary emission from the surrounding plasma on scales much smaller than those probed by the currently largest baselines. We show that for aliasing not to affect images reconstructed with regularized maximum likelihood methods and model images that are directly compared to the data, the sampling of these images (i.e., their pixel spacing) needs to be significantly finer than the scale probed by the largest baseline in the array. Using GRMHD simulations of black-hole images, we estimate the maximum allowable pixel spacing to be approximately equal to (1/8)GM/c^2; for both of the primary EHT targets, this corresponds to an angular pixel size of <0.5 microarcseconds. With aliasing under control, we then advocate use of the second-order Butterworth filter with a cut-off scale equal to the maximum array baseline as optimal for visualizing the reconstructed images. In contrast to the traditional Gaussian filters, this Butterworth filter retains most of the power at the scales probed by the array while suppressing the fine image details for which no data exist.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The neutron star equation of state (EOS) is now being constrained from a diverse set of multi-messenger data, including gravitational waves from binary neutron star mergers, X-ray observations of the neutron star radius, and many types of laboratory nuclear experiments. These measurements are typically mapped to a common domain -- either to a corresponding radius or to a parametrized EOS using a Bayesian inference scheme -- for comparison with one another. We explore here the statistical biases that can arise when such multi-messenger data are mapped to a common domain for comparison. We find that placing Bayesian priors individually in each domain of measurement can transform to biased constraints in the domain of comparison. Using the first two binary neutron star mergers as an example, we show that a uniform prior in the tidal deformability can produce artificial evidence for large radii, which the data do not support. We present a new prescription for defining Bayesian priors in any domain of measurement, that will allow for minimally-biased constraints in the domain of comparison. Finally, using this new prescription, we provide a status update on multi-messenger EOS constraints on the neutron star radius.

Closure phases are critical in astronomical interferometry. However, their uncertainties are difficult to compute numerically. We provide a method to efficiently compute interferometric closure phase distributions in terms of an approximate distribution that is valid in the low signal-to-noise ratio regime. This is done by first showing that the true phase distribution is well approximated by the von Mises distribution, then performing a convolution of three von Mises distributions. The resulting approximation is superior than the normal distribution for all signal-to-noise ratios and, being fully analytic, allow for fast computations in statistical algorithms.

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.

Thermal X-ray emission from rotation-powered pulsars is believed to originate from localized "hotspots" on the stellar surface occurring where large-scale currents from the magnetosphere return to heat the atmosphere. Lightcurve modeling has primarily been limited to simple models, such as circular antipodal emitting regions with constant temperature. We calculate more realistic temperature distributions within the polar caps, taking advantage of recent advances in magnetospheric theory, and we consider their effect on the predicted lightcurves. The emitting regions are non-circular even for a pure dipole magnetic field, and the inclusion of an aligned magnetic quadrupole moment introduces a north-south asymmetry. As the aligned quadrupole moment is increased, one hotspot grows in size before becoming a thin ring surrounding the star. For the pure dipole case, moving to the more realistic model changes the lightcurves by $5-10\%$ for millisecond pulsars, helping to quantify the systematic uncertainty present in current dipolar models. Including the quadrupole gives considerable freedom in generating more complex lightcurves. We explore whether these simple dipole+quadrupole models can account for the qualitative features of the lightcurve of PSR J0437$-$4715.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.

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.

Observations of isolated neutron stars place constraints on the equation of state (EOS) of cold, neutron-rich matter, while nuclear physics experiments probe the EOS of hot, symmetric matter. Many dynamical phenomena, such as core-collapse supernovae, the formation and cooling of proto-neutron stars, and neutron star mergers, lie between these two regimes and depend on the EOS at finite temperatures for matter with varying proton fractions. In this paper, we introduce a new framework to accurately calculate the thermal pressure of neutron-proton-electron matter at arbitrary density, temperature, and proton fraction. This framework can be expressed using a set of five physically-motivated parameters that span a narrow range of values for realistic EOS and are able to capture the leading-order effects of degenerate matter on the thermal pressure. We base two of these parameters on a new approximation of the Dirac effective mass, with which we reproduce the thermal pressure to within <~30% for a variety of realistic EOS at densities of interest. Three additional parameters, based on the behavior of the symmetry energy near the nuclear saturation density, allow for the extrapolation of any cold EOS in beta-equilibrium to arbitrary proton fractions. Our model thus allows a user to extend any cold nucleonic EOS, including piecewise-polytropes, to arbitrary temperature and proton fraction, for use in calculations and numerical simulations of astrophysical phenomena. We find that our formalism is able to reproduce realistic finite-temperature EOS with errors of <~20% and offers a 1-3 orders-of-magnitude improvement over existing ideal-fluid models.

The Neutron-star Interior Composition ExploreR (NICER) is collecting data to measure the radii of neutron stars by observing the pulsed emission from their surfaces. The primary targets are isolated, rotation-powered pulsars, in which the surface polar caps are heated by bombardment from magnetospheric currents of electrons and positrons. We investigate various stopping mechanisms for the beams of particles that bombard the atmosphere and calculate the heat deposition, the atmospheric temperature profiles, and the energy spectra and beaming of the emerging radiation. We find that low-energy particles with {\gamma} $\sim 2-10$ deposit most of their energy in the upper regions of the atmosphere, at low optical depth, resulting in beaming patterns that are substantially different than those of deep-heated, radiative equilibrium models. Only particles with energies {\gamma} $\gtrsim 50$ penetrate to high optical depths and fulfill the conditions necessary for a deep-heating approximation. We discuss the implications of our work for modeling the pulse profiles from rotation-powered pulsars and for the inference of their radii with NICER observations.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

Synthesis imaging of the black hole in the center of the Milky Way, Sgr A*, with the Event Horizon Telescope (EHT) rests on the assumption of a stationary image. We explore the limitations of this assumption using high-cadence GRMHD simulations of Sgr A*. We employ analytic models that capture the basic characteristics of the images to understand the origin of the variability in the simulated visibility amplitudes. We find that, in all simulations, the visibility amplitudes for baselines oriented perpendicular to the spin axis of the black hole typically decrease smoothly over baseline lengths that are comparable to those of the EHT. On the other hand, the visibility amplitudes for baselines oriented parallel to the spin axis show significant structure with one or more minima. This suggests that fitting EHT observations with geometric models will lead to reasonably accurate determination of the orientation of the black-hole on the plane of the sky. However, in the disk-dominated models, the locations and depths of the minima in the visibility amplitudes depend primarily on the width and asymmetry of the crescent-like images and are highly variable. In the jet-dominated models, the locations of the minima are determined by the separation of the two image components but their depths depend primarily on the relative brightness of the two components and are also variable. This suggests that using time-independent models to infer additional black-hole parameters, such as the shadow size or the spin magnitude, will be severely affected by the variability of the accretion flow.

The Event Horizon Telescope is a millimeter VLBI array that aims to take the first pictures of the black holes in the center of the Milky Way and of the M87 galaxy, with horizon scale resolution. Measurements of the shape and size of the shadows cast by the black holes on the surrounding emission can test the cosmic censorship conjecture and the no-hair theorem and may find evidence for classical effects of the quantum structure of black holes. Observations of coherent structures in the accretion flows may lead to accurate measurements of the spins of the black holes and of other properties of their spacetimes. For Sgr A*, the black hole in the center of the Milky Way, measurements of the precession of stellar orbits and timing monitoring of orbiting pulsars offer complementary avenues to the gravitational tests with the Event Horizon Telescope.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

Scattering in the ionized interstellar medium is commonly observed to be anisotropic, with theories of magnetohydrodynamic (MHD) turbulence explaining the anisotropy through a preferred magnetic field direction throughout the scattering regions. In particular, the line of sight to the Galactic Center supermassive black hole, Sgr A*, exhibits strong and anisotropic scattering, which dominates its observed size at wavelengths of a few millimeters and longer. Therefore, inferences of the intrinsic structure of \sgra\ at these wavelengths are sensitive to the assumed scattering model. In addition, extrapolations of the scattering model from long wavelengths, at which its parameters are usually estimated, to 1.3 mm, where the Event Horizon Telescope (EHT) seeks to image Sgr A* on Schwarzschild-radius scales, are also sensitive to the assumed scattering model. Past studies of Sgr A* have relied on simple Gaussian models for the scattering kernel that effectively presume an inner scale of turbulence far greater than the diffractive scale; this assumption is likely violated for Sgr A* at 1.3 mm. We develop a physically motivated model for anisotropic scattering, using a simplified model for MHD turbulence with a finite inner scale and a wandering transverse magnetic field direction. We explore several explicit analytic models for this wandering and derive the expected observational properties --- scatter broadening and refractive scintillation --- for each. For expected values of the inner scale, the scattering kernel for all models is markedly non-Gaussian at 1.3 mm but is straightforward to calculate and depends only weakly on the assumed model for the wandering of the magnetic field direction. On the other hand, in all models, the refractive substructure depends strongly on the wandering model and may be an important consideration in imaging Sgr A* with the EHT.

We explore the use of principal component analysis (PCA) to characterize high-fidelity simulations and interferometric observations of the millimeter emission that originates near the horizons of accreting black holes. We show mathematically that the Fourier transforms of eigenimages derived from PCA applied to an ensemble of images in the spatial-domain are identical to the eigenvectors of PCA applied to the ensemble of the Fourier transforms of the images, which suggests that this approach may be applied to modeling the sparse interferometric Fourier-visibilities produced by an array such as the Event Horizon Telescope (EHT). We also show that the simulations in the spatial domain themselves can be compactly represented with a PCA-derived basis of eigenimages allowing for detailed comparisons between variable observations and time-dependent models, as well as for detection of outliers or rare events within a time series of images. Furthermore, we demonstrate that the spectrum of PCA eigenvalues is a diagnostic of the power spectrum of the structure and, hence, of the underlying physical processes in the simulated and observed images.

Gravitational waves from the coalescence of two neutron stars were recently detected for the first time by the LIGO-Virgo collaboration, in event GW170817. This detection placed an upper limit on the effective tidal deformability of the two neutron stars and tightly constrained the chirp mass of the system. We report here on a new simplification that arises in the effective tidal deformability of the binary, when the chirp mass is specified. We find that, in this case, the effective tidal deformability of the binary is surprisingly independent of the component masses of the individual neutron stars, and instead depends primarily on the ratio of the chirp mass to the neutron star radius. Thus, a measurement of the effective tidal deformability can be used to directly measure the neutron star radius. We find that the upper limit on the effective tidal deformability from GW170817 implies that the radius cannot be larger than ~13km, at the 90% level, independent of the assumed masses for the component stars. The result can be applied generally, to probe the stellar radii in any neutron star-neutron star merger with a measured chirp mass. The approximate mass-independence disappears for neutron star-black hole mergers. Finally, we discuss a Bayesian inference of the equation of state that uses the measured chirp mass and tidal deformability from GW170817 combined with nuclear and astrophysical priors and discuss possible statistical biases in this inference.

The need for a consistent quantum evolution for black holes has led to proposals that their semiclassical description is modified not just near the singularity, but at horizon or larger scales. If such modifications extend beyond the horizon, they influence regions accessible to distant observeration. Natural candidates for these modifications behave like metric fluctuations, with characteristic length and time scales set by the horizon radius. We investigate the possibility of using the Event Horizon Telescope to observe these effects, if they have a strength sufficient to make quantum evolution consistent with unitarity. We find that such quantum fluctuations can introduce a strong time dependence for the shape and size of the shadow that a black hole casts on its surrounding emission. For the black hole in the center of the Milky Way, detecting the rapid time variability of its shadow will require non-imaging timing techniques. However, for the much larger black hole in the center of the M87 galaxy, a variable black-hole shadow, if present with these parameters, would be readily observable in the individual snapshots that will be obtained by the Event Horizon Telescope.

We use VLBA+VLA observations to measure the sizes of the scatter-broadened images of 6 of the most heavily scattered known pulsars: 3 within the Galactic Centre (GC) and 3 elsewhere in the inner Galactic plane. By combining the measured sizes with temporal pulse broadening data from the literature and using the thin-screen approximation, we locate the scattering medium along the line of sight to these 6 pulsars. At least two scattering screens are needed to explain the observations of the GC sample. We show that the screen inferred by previous observations of SGR J1745-2900 and Sgr A*, which must be located far from the GC, falls off in strength on scales < 0.2 degree. A second scattering component closer to (< 2 kpc) or even (tentatively) within (< 700 pc) the GC produces most or all of the temporal broadening observed in the other GC pulsars. Outside the GC, the scattering locations for all three pulsars are ~2 kpc from Earth, consistent with the distance of the Carina-Sagittarius or Scutum spiral arm. For each object the 3D scattering origin coincides with a known HII region (and in one case also a supernova remnant), suggesting that such objects preferentially cause the intense interstellar scattering seen towards the Galactic plane. We show that the HII regions should contribute > 25% of the total dispersion measure (DM) towards these pulsars, and calculate reduced DM distances. Those distances for other pulsars lying behind HII regions may be similarly overestimated.

Fast and accurate integration of geodesics in Kerr spacetimes is an important tool in modeling the orbits of stars and the transport of radiation in the vicinities of black holes. Most existing integration algorithms employ Boyer-Lindquist coordinates, which have coordinate singularities at the event horizon and along the poles. Handling the singularities requires special numerical treatment in these regions, often slows down the calculations, and may lead to inaccurate geodesics. We present here a new general-purpose geodesic integrator, GRay2, that overcomes these issues by employing the Cartesian form of Kerr-Schild coordinates. By performing particular mathematical manipulations of the geodesic equations and several optimizations, we develop an implementation of the Cartesian Kerr-Schild coordinates that outperforms calculations that use the seemingly simpler equations in Boyer-Lindquist coordinates. We also employ the OpenCL framework, which allows GRay2 to run on multi-core CPUs as well as on a wide range of GPU hardware accelerators, making the algorithm more versatile. We report numerous convergence tests and benchmark results for GRay2 for both time-like (particle) and null (photon) geodesics.

Non-ideal MHD effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative GRMHD simulations of black hole accretion flows to identify and explore the structures and properties of current sheets as potential sites of magnetic reconnection. In the case of standard and normal (SANE) disks, we find that, in the reconnection sites, the plasma beta ranges from $0.1$ to $1000$, the magnetization ranges from $10^{-4}$ to $1$, and the guide fields are weak compared to the reconnecting fields. In magnetically arrested (MAD) disks, we find typical values for plasma beta from $10^{-2}$ to $10^3$, magnetizations from $10^{-3}$ to $10$, and typically stronger guide fields, with strengths comparable to or greater than the reconnecting fields. These are critical parameters that govern the electron energy distribution resulting from magnetic reconnection and can be used in the context of plasma simulations to provide microphysics inputs to global simulations. We also find that ample magnetic energy is available in the reconnection regions to power the fluence of bright X-ray flares observed from the black hole in the center of the Milky Way.

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. We present here a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parametrized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, rho_sat, can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ~30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass-radius relation can lead to biases of nearly 1 km in radius, towards larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias.

The Event Horizon Telescope (EHT) is a millimeter-wavelength, very-long-baseline interferometry (VLBI) experiment that is capable of observing black holes with horizon-scale resolution. Early observations have revealed variable horizon-scale emission in the Galactic Center black hole, Sagittarius A* (Sgr A*). Comparing such observations to time-dependent general relativistic magnetohydrodynamic (GRMHD) simulations requires statistical tools that explicitly consider the variability in both the data and the models. We develop here a Bayesian method to compare time-resolved simulation images to variable VLBI data, in order to infer model parameters and perform model comparisons. We use mock EHT data based on GRMHD simulations to explore the robustness of this Bayesian method and contrast it to approaches that do not consider the effects of variability. We find that time-independent models lead to offset values of the inferred parameters with artificially reduced uncertainties. Moreover, neglecting the variability in the data and the models often leads to erroneous model selections. We finally apply our method to the early EHT data on Sgr A*.

The increasing number and precision of measurements of neutron star masses, radii, and, in the near future, moments of inertia offer the possibility of precisely determining the neutron star equation of state. One way to facilitate the mapping of observables to the equation of state is through a parametrization of the latter. We present here a generic method for optimizing the parametrization of any physically allowed EoS. We use mock equations of state that incorporate physically diverse and extreme behavior to test how well our parametrization reproduces the global properties of the stars, by minimizing the errors in the observables mass, radius, and the moment of inertia. We find that using piecewise polytropes and sampling the EoS with five fiducial densities between ~1-8 times the nuclear saturation density results in optimal errors for the smallest number of parameters. Specifically, it recreates the radii of the assumed EoS to within less than 0.5 km for the extreme mock equations of state and to within less than 0.12 km for 95% of a sample of 42 proposed, physically-motivated equations of state. Such a parametrization is also able to reproduce the maximum mass to within 0.04 M_sun and the moment of inertia of a 1.338 M_sun neutron star to within less than 10% for 95% of the proposed sample of equations of state.

The observable quantities that carry the most information regarding the structures of the images of black holes in the interferometric observations with the Event Horizon Telescope are the closure phases along different baseline triangles. We use long time span, high cadence, GRMHD+radiative transfer models of Sgr A$^*$ to investigate the expected variability of closure phases in such observations. We find that, in general, closure phases along small baseline triangles show little variability, except in the cases when one of the triangle vertices crosses one of a small regions of low visibility amplitude. The closure phase variability increases with the size of the baseline triangle, as larger baselines probe the small-scale structures of the images, which are highly variable. On average, the jet-dominated MAD models show less closure phase variability than the disk-dominated SANE models, even in the large baseline triangles, because the images from the latter are more sensitive to the turbulence in the accretion flow. Our results suggest that image reconstruction techniques need to explicitly take into account the closure phase variability, especially if the quality and quantity of data allow for a detailed characterization of the nature of variability. This also implies that, if image reconstruction techniques that rely on the assumption of a static image are utilized, regions of the $u-v$ space that show a high level of variability will need to be identified and excised.

The Event Horizon Telescope is a global very-long baseline interferometer capable of probing potential deviations from the Kerr metric, which is believed to provide the unique description of astrophysical black holes. Here we report an updated constraint on the quadrupolar deviation of Sagittarius A* within the context of a radiatively inefficient accretion flow model in a quasi-Kerr background. We also simulate near-future constraints obtainable by the forthcoming eight-station array and show that in this model already a one-day observation can measure the spin magnitude to within $0.005$, the inclination to within $0.09^\circ$, the position angle to within $0.04^\circ$, and the quadrupolar deviation to within $0.005$ at $3\sigma$ confidence. Thus, we are entering an era of high-precision strong gravity measurements.

Significant X-ray variability and flaring has been observed from Sgr A* but is poorly understood from a theoretical standpoint. We perform GRMHD simulations that take into account a population of non-thermal electrons with energy distributions and injection rates that are motivated by PIC simulations of magnetic reconnection. We explore the effects of including these non-thermal electrons on the predicted broadband variability of Sgr A* and find that X-ray variability is a generic result of localizing non-thermal electrons to highly magnetized regions, where particles are likely to be accelerated via magnetic reconnection. The proximity of these high-field regions to the event horizon forms a natural connection between IR and X-ray variability and accounts for the rapid timescales associated with the X-ray flares. The qualitative nature of this variability is consistent with observations, producing X-ray flares that are always coincident with IR flares, but not vice versa, i.e., there are a number of IR flares without X-ray counterparts.

Many techniques for measuring neutron star radii rely on absolute flux measurements in the X-rays. As a result, one of the fundamental uncertainties in these spectroscopic measurements arises from the absolute flux calibrations of the detectors being used. Using the stable X-ray burster, GS 1826-238, and its simultaneous observations by Chandra HETG/ACIS-S and RXTE/PCA as well as by XMM-Newton EPIC-pn and RXTE/PCA, we quantify the degree of uncertainty in the flux calibration by assessing the differences between the measured fluxes during bursts. We find that the RXTE/PCA and the Chandra gratings measurements agree with each other within their formal uncertainties, increasing our confidence in these flux measurements. In contrast, XMM-Newton EPIC-pn measures 14.0$\pm$0.3 % less flux than the RXTE/PCA. This is consistent with the previously reported discrepancy with the flux measurements of EPIC-pn, compared to EPIC-MOS1, MOS2 and ACIS-S detectors. We also show that any intrinsic time dependent systematic uncertainty that may exist in the calibration of the satellites has already been implicity taken into account in the neutron star radius measurements.

A precise moment of inertia measurement for PSR J0737-3039A in the double pulsar system is expected within the next five years. We present here a new method of mapping the anticipated measurement of the moment of inertia directly into the neutron star structure. We determine the maximum and minimum values possible for the moment of inertia of a neutron star of a given radius based on physical stability arguments, assuming knowledge of the equation of state only at densities below the nuclear saturation density. If the equation of state is trusted up to the nuclear saturation density, we find that a measurement of the moment of inertia will place absolute bounds on the radius of PSR J0737-3039A to within $\pm$1 km. The resulting combination of moment of inertia, mass, and radius measurements for a single source will allow for new, stringent constraints on the dense-matter equation of state.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

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 present a comprehensive study of spectroscopic radius measurements of twelve neutron stars obtained during thermonuclear bursts or in quiescence. We incorporate, for the first time, a large number of systematic uncertainties in the measurement of the apparent angular sizes, Eddington fluxes, and distances, in the composition of the interstellar medium, and in the flux calibration of X-ray detectors. We also take into account the results of recent theoretical calculations of rotational effects on neutron star radii, of atmospheric effects on surface spectra, and of relativistic corrections to the Eddington critical flux. We employ Bayesian statistical frameworks to obtain neutron star radii from the spectroscopic measurements as well as to infer the equation of state from the radius measurements. Combining these with the results of experiments in the vicinity of nuclear saturation density and the observations of ~2 Msun neutron stars, we place strong and quantitative constraints on the properties of the equation of state between approximately 2-8 times the nuclear saturation density. We find that around M=1.5 Msun, the preferred equation of state predicts radii between 10.1 - 11.1 km. When interpreting the pressure constraints in the context of high density equations of state based on interacting nucleons, our results suggest a relatively weak contribution of the three-body interaction potential.

The Neutron-star Interior Composition Explorer (NICER) is an X-ray astrophysics payload that will be placed on the International Space Station. Its primary science goal is to measure with high accuracy the pulse profiles that arise from the non-uniform thermal surface emission of rotation-powered pulsars. Modeling general relativistic effects on the profiles will lead to measuring the radii of these neutron stars and to constraining their equation of state. Achieving this goal will depend, among other things, on accurate knowledge of the source, sky, and instrument backgrounds. We use here simple analytic estimates to quantify the level at which these backgrounds need to be known in order for the upcoming measurements to provide significant constraints on the properties of neutron stars. We show that, even in the minimal-information scenario, knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 allows for a measurement of the neutron star compactness to better than 10% uncertainty for most of the parameter space. These constraints improve further when more realistic assumptions are made about the neutron star emission and spin, and when additional information about the source itself, such as its mass or distance, are incorporated.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

The half opening angle of a Kerr black-hole shadow is always equal to (5+-0.2)GM/Dc^2, where M is the mass of the black hole and D is its distance from the Earth. Therefore, measuring the size of a shadow and verifying whether it is within this 4% range constitutes a null hypothesis test of General Relativity. We show that the black hole in the center of the Milky Way, Sgr A*, is the optimal target for performing this test with upcoming observations using the Event Horizon Telescope. We use the results of optical/IR monitoring of stellar orbits to show that the mass-to-distance ratio for Sgr A* is already known to an accuracy of +-4%. We investigate our prior knowledge of the properties of the scattering screen between Sgr A and the Earth, the effects of which will need to be corrected for in order for the black-hole shadow to appear sharp against the background emission. Finally, we explore an edge detection scheme for interferometric data and a pattern matching algorithm based on the Hough/Radon transform and demonstrate that the shadow of the black hole at 1.3 mm can be localized, in principle, to within ~9%. All these results suggest that our prior knowledge of the properties of the black hole, of scattering broadening, and of the accretion flow can only limit this General Relativistic null hypothesis test with Event Horizon Telescope observations of Sgr A* to 10%.

One avenue for testing the no-hair theorem is obtained through timing a pulsar orbiting close to a black hole and fitting for quadrupolar effects on the time-of-arrival of pulses. If deviations from the Kerr quadrupole are measured, then the no-hair theorem is invalidated. To this end, we derive an expression for the light travel time delay for a pulsar orbiting in a black-hole spacetime described by the Butterworth-Ipser metric, which has an arbitrary spin and quadrupole moment. We consider terms up to the quadrupole order in the black-hole metric and derive the time-delay expression in a closed analytic form. This allows for fast computations that are useful in fitting time-of-arrival observations of pulsars orbiting close to astrophysical black holes.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

We describe HEROIC, an upgraded version of the relativistic radiative post-processor code HERO described in a previous paper, but which now Includes Comptonization. HEROIC models Comptonization via the Kompaneets equation, using a quadratic approximation for the source function in the short characteristics radiation solver. It employs a simple form of accelerated lambda iteration to handle regions of high scattering opacity. In addition to solving for the radiation field, HEROIC also solves for the gas temperature by applying the condition of radiative equilibrium. We present benchmarks and tests of the Comptonization module in HEROIC with simple 1D and 3D scattering problems. We also test the ability of the code to handle various relativistic effects using model atmospheres and accretion flows in a black hole space-time. We present two applications of HEROIC to general relativistic MHD simulations of accretion discs. One application is to a thin accretion disc around a black hole. We find that the gas below the photosphere in the multi-dimensional HEROIC solution is nearly isothermal, quite different from previous solutions based on 1D plane parallel atmospheres. The second application is to a geometrically thick radiation-dominated accretion disc accreting at 11 times the Eddington rate. The multi-dimensional HEROIC solution shows that, for observers who are on axis and look down the polar funnel, the isotropic equivalent luminosity could be more than ten times the Eddington limit, even though the spectrum might still look thermal and show no signs of relativistic beaming.

The black hole in the center of the Milky Way, Sgr A*, has the largest mass-to-distance ratio among all known black holes in the Universe. This property makes Sgr A* the optimal target for testing the gravitational no-hair theorem. In the near future, major developments in instrumentation will provide the tools for high-precision studies of its spacetime via observations of relativistic effects in stellar orbits, in the timing of pulsars, and in horizon-scale images of its accretion flow. We explore here the prospect of measuring the properties of the black-hole spacetime using all these three types of observations. We show that the correlated uncertainties in the measurements of the black-hole spin and quadrupole moment using the orbits of stars and pulsars are nearly orthogonal to those obtained from measuring the shape and size of the shadow the black hole casts on the surrounding emission. Combining these three types of observations will, therefore, allow us to assess and quantify systematic biases and uncertainties in each measurement and lead to a highly accurate, quantitative test of the gravitational no-hair theorem.

We calculate the effects of spot size on pulse profiles of moderately rotating neutron stars. Specifically, we quantify the bias introduced in radius measurements from the common assumption that spots are infinitesimally small. We find that this assumption is reasonable for spots smaller than 10$^\circ$-18$^\circ$ and leads to errors that are $\le$10% in the radius measurement, depending on the location of the spot and the inclination of the observer. We consider the implications of our results for neutron star radius measurements with the upcoming and planned X-ray missions NICER and LOFT. We calculate the expected spot size for different classes of sources and investigate the circumstances under which the assumption of a small spot is justified.

Data selection and the determination of systematic uncertainties in the spectroscopic measurements of neutron star radii from thermonuclear X-ray bursts have been the subject of numerous recent studies. In one approach, the uncertainties and outliers were determined by a data-driven Bayesian mixture model, whereas in a second approach, data selection was performed by requiring that the observations follow theoretical expectations. We show here that, due to inherent limitations in the data, the theoretically expected trends are not discernible in the majority of X-ray bursts even if they are present. Therefore, the proposed theoretical selection criteria are not practical with the current data for distinguishing clean data sets from outliers. Furthermore, when the data limitations are not taken into account, the theoretically motivated approach selects a small subset of bursts with properties that are in fact inconsistent with the underlying assumptions of the method. We conclude that the data-driven selection methods do not suffer from the limitations of this theoretically motivated one.

Measuring neutron star radii with spectroscopic and timing techniques relies on the combination of multiple observables to break the degeneracies between the mass and radius introduced by general relativistic effects. Here, we explore a previously used frequentist and a newly proposed Bayesian framework to obtain the most likely value and the uncertainty in such a measurement. We find that, for the expected range of masses and radii and for realistic measurement errors, the frequentist approach suffers from biases that are larger than the accuracy in the radius measurement required to distinguish between the different equations of state. In contrast, in the Bayesian framework, the inferred uncertainties are larger, but the most likely values do not suffer from such biases. We also investigate ways of quantifying the degree of consistency between different spectroscopic measurements from a single source. We show that a careful assessment of the systematic uncertainties in the measurements eliminates the need for introducing ad hoc biases, which lead to artificially large inferred radii.

The current effort to test General Relativity employs multiple disparate formalisms for different observables, obscuring the relations between laboratory, astrophysical and cosmological constraints. To remedy this situation, we develop a parameter space for comparing tests of gravity on all scales in the universe. In particular, we present new methods for linking cosmological large-scale structure, the Cosmic Microwave Background and gravitational waves with classic PPN tests of gravity. Diagrams of this gravitational parameter space reveal a noticeable untested regime. The untested window, which separates small-scale systems from the troubled cosmological regime, could potentially hide the onset of corrections to General Relativity.

HERO (Hybrid Evaluator for Radiative Objects) is a 3D general relativistic radiative transfer code which has been tailored to the problem of analyzing radiation from simulations of relativistic accretion discs around black holes. HERO is designed to be used as a postprocessor. Given some fixed fluid structure for the disc (i.e. density and velocity as a function of position from a hydrodynamics or magnetohydrodynamics simulation), the code obtains a self-consistent solution for the radiation field and for the gas temperatures using the condition of radiative equilibrium. The novel aspect of HERO is that it combines two techniques: 1) a short characteristics (SC) solver that quickly converges to a self consistent disc temperature and radiation field, with 2) a long characteristics (LC) solver that provides a more accurate solution for the radiation near the photosphere and in the optically thin regions. By combining these two techniques, we gain both the computational speed of SC and the high accuracy of LC. We present tests of HERO on a range of 1D, 2D and 3D problems in flat space and show that the results agree well with both analytical and benchmark solutions. We also test the ability of the code to handle relativistic problems in curved space. Finally, we discuss the important topic of ray-defects, a major limitation of the SC method, and describe our strategy for minimizing the induced error.

We explore the variability properties of long, high cadence GRMHD simulations across the electromagnetic spectrum using an efficient, GPU-based radiative transfer algorithm. We focus on both disk- and jet-dominated simulations with parameters that successfully reproduce the time-averaged spectral properties of Sgr A* and the size of its image at 1.3mm. We find that the disk-dominated models produce short timescale variability with amplitudes and power spectra that closely resemble those inferred observationally. In contrast, jet-dominated models generate only slow variability, at lower flux levels. Neither set of models show any X-ray flares, which most likely indicate that additional physics, such as particle acceleration mechanisms, need to be incorporated into the GRMHD simulations to account for them. The disk-dominated models show strong, short-lived mm/IR flares, with short (<~ 1hr) time lags between the mm and IR wavelengths, that arise from strong-field gravitational lensing of magnetic flux tubes near the horizon. Such events provide a natural explanation for the observed IR flares with no X-ray counterparts.

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.

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black-holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast GPU-based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all twelve models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black-holes and their thermodynamic properties.

We present a new numerical algorithm for the calculation of pulse profiles from spinning neutron stars in the Hartle-Thorne approximation. Our approach allows us to formally take into account the effects of Doppler shifts and aberration, of frame dragging, as well as of the oblateness of the stellar surface and of its quadrupole moment. We confirm an earlier result that neglecting the oblateness of the neutron-star surface leads to ~5-30% errors in the calculated profiles and further show that neglecting the quadrupole moment of its spacetime leads to ~1-5% errors at a spin frequency of 600 Hz. We discuss the implications of our results for the measurements of neutron-star masses and radii with upcoming X-ray missions, such as NASA's NICER and ESA's LOFT.

Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

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 calculate the rotational broadening in the observed thermal spectra of neutron stars spinning at moderate rates in the Hartle-Thorne approximation. These calculations accurately account for the effects of the second-order Doppler boosts as well as for the oblate shapes and the quadrupole moments of the neutron stars. We find that fitting the spectra and inferring the bolometric fluxes under the assumption that a star is not rotating causes an underestimate of the inferred fluxes and, thus, radii. The correction depends on the stellar spin, radius, and observer's inclination. For a 10 km neutron star spinning at 600 Hz, the rotational correction to the flux is ~1-4%, while for a 15 km neutron star with the same spin period, the correction ranges from 2% for pole-on sources to 12% for edge-on sources. We calculate the inclination-averaged corrections to inferred radii as a function of the neutron-star radius and mass and provide an empirical formula for the corrections. For realistic neutron star parameters (1.4 M$_\odot$, 12 km, 600 Hz), the stellar radius is on the order of 4% larger than the radius inferred under the assumption that the star is not spinning.

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.

Black hole binaries exhibit a wide range of variability phenomena, from large-scale state changes to broadband noise and quasi-periodic oscillations, but the physical nature of much of this variability is poorly understood. We examine the variability properties of three GRMHD simulations of thin accretion disks around black holes of varying spin, producing light curves and power spectra as would be seen by observers. We find that the simulated power spectra show a broad feature at high frequency, which increases in amplitude with the inclination of the observer. We show that this high-frequency feature is a product of the Doppler effect and that its location is a function of the mass and spin of the black hole. This Doppler feature demonstrates that power spectral properties of the accretion disk can be tied to, and potentially used to determine, physical properties of the black hole.

The advent of the Event Horizon Telescope (EHT), a millimeter-wave very-long baseline interferometric array, has enabled spatially-resolved studies of the sub-horizon-scale structure for a handful of supermassive black holes. Among these, the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), presents the largest angular cross section. Thus far, these studies have focused upon measurements of the black hole spin and the validation of low-luminosity accretion models. However, a critical input into the analysis of EHT data is the structure of the black hole spacetime, and thus these observations provide the novel opportunity to test the applicability of the Kerr metric to astrophysical black holes. Here we present the first simulated images of a radiatively inefficient accretion flow (RIAF) around Sgr A* employing a quasi-Kerr metric that contains an independent quadrupole moment in addition to the mass and spin that fully characterize a black hole in general relativity. We show that these images can be significantly different from the images of a RIAF around a Kerr black hole with the same spin and demonstrate the feasibility of testing the no-hair theorem by constraining the quadrupolar deviation from the Kerr metric with existing EHT data. Equally important, we find that the disk inclination and spin orientation angles are robust to the inclusion of additional parameters, providing confidence in previous estimations assuming the Kerr metric based upon EHT observations. However, at present the limits upon potential modifications of the Kerr metric remain weak.

Using stellar structure calculations in the Hartle-Thorne approximation, we derive analytic expressions connecting the ellipticity of the stellar surface to the compactness, the spin angular momentum, and the quadrupole moment of the spacetime. We also obtain empirical relations between the compactness, the spin angular momentum, and the spacetime quadrupole. Our formulae reproduce the results of numerical calculations to within a few percent and help reduce the number of parameters necessary to model the observational appearance of moderately spinning neutron stars. This is sufficient for comparing theoretical spectroscopic and timing models to observations that aim to measure the masses and radii of neutron stars and to determine the equation of state prevailing in their interiors.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

According to the no-hair theorem, astrophysical black holes are fully characterized by their masses and spins and are described by the Kerr metric. This theorem can be tested observationally by measuring (at least) three different multipole moments of the spacetimes of black holes. In this paper, we calculate the profiles of fluorescent iron lines emitted from the accretion flows around black hole candidates within a framework that allows us to perform the calculation as a function of its mass and spin as well as of a free parameter that measures potential deviations from the Kerr metric. We show that such deviations lead to line profiles that are significantly altered and may exhibit a modified flux ratio of the two peaks in their characteristic double-peaked shape. We also show that the disk inclination can be measured independently of the spin and the deviation parameter at low to intermediate inclination angles as in the case of Kerr black holes. We estimate the precision that near-future X-ray missions such as Astro-H and ATHENA+ are required to achieve in order to resolve deviations from the Kerr metric in iron line profiles and show that constraints on such deviations will be strongest for rapidly spinning black holes. More generally, we show that measuring the line profile with a precision of ~5% at disk inclinations of 30 or 60 degrees constraints the deviation parameter to order unity for values of the spin a>0.5M.

White dwarfs, neutron stars and stellar mass black holes are key laboratories to study matter in most extreme conditions of gravity and magnetic field. The unprecedented effective area of Athena+ will allow us to advance our understanding of emission mechanisms and accretion physics over a wide range of mass accretion rates, starting from lower and sub-luminous quiescent X-ray binaries up to super-Eddington ultra-luminous sources. Athena+ will measure stellar black hole spins in a much higher number of binaries than achievable now, opening the possibility to study how spin varies with black hole history. The high throughput and energy resolution of the X-IFU will be instrumental in establishing how disc wind properties depend on accretion state, in determining wind launching mechanism and in quantifying the impact of the wind induced mass loss on binary evolution and environment. Triggers and high quality optical and radio data originating from large wide field contemporaneous instruments will provide essential complementary information on jet launching mechanisms and on the physics of rotation powered pulsars, for instance. In addition, Athena+ will furnish multiple, independent measurements of the neutron star mass/radius relation in a wide range of environments and conditions so as to constrain the debated equation of state.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

We introduce GRay, a massively parallel integrator designed to trace the trajectories of billions of photons in a curved spacetime. This GPU-based integrator employs the stream processing paradigm, is implemented in CUDA C/C++, and runs on nVidia graphics cards. The peak performance of GRay using single precision floating-point arithmetic on a single GPU exceeds 300 GFLOP (or 1 nanosecond per photon per time step). For a realistic problem, where the peak performance cannot be reached, GRay is two orders of magnitude faster than existing CPU-based ray tracing codes. This performance enhancement allows more effective searches of large parameter spaces when comparing theoretical predictions of images, spectra, and lightcurves from the vicinities of compact objects to observations. GRay can also perform on-the-fly ray tracing within general relativistic magnetohydrodynamic algorithms that simulate accretion flows around compact objects. Making use of this algorithm, we calculate the properties of the shadows of Kerr black holes and the photon rings that surround them. We also provide accurate fitting formulae of their dependencies on black hole spin and observer inclination, which can be used to interpret upcoming observations of the black holes at the center of the Milky Way, as well as M87, with the Event Horizon Telescope.

Neutron stars spinning at moderate rates (~300-600Hz) become oblate in shape and acquire a nonzero quadrupole moment. In this paper, we calculate profiles of atomic features from such neutron stars using a ray-tracing algorithm in the Hartle-Thorne approximation. We show that line profiles acquire cores that are much narrower than the widths expected from pure Doppler effects for a large range of observer inclinations. As a result, the effects of both the oblateness and the quadrupole moments of neutron stars need to be taken into account when aiming to measure neutron star radii from rotationally broadened lines. Moreover, the presence of these narrow cores substantially increases the likelihood of detecting atomic lines from rapidly spinning neutron stars.

Monitoring the orbits of stars around Sgr A* offers the possibility of detecting the precession of their orbital planes due to frame dragging, of measuring the spin and quadrupole moment of the black hole, and of testing the no-hair theorem. Here we investigate whether the deviations of stellar orbits from test-particle trajectories due to wind mass loss and tidal dissipation of the orbital energy compromise such measurements. We find that the effects of stellar winds are, in general, negligible. On the other hand, for the most eccentric orbits (e>0.96) for which an optical interferometer, such as GRAVITY, will detect orbital plane precession due to frame dragging, the tidal dissipation of orbital energy occurs at timescales comparable to the timescale of precession due to the quadrupole moment of the black hole. As a result, this non-conservative effect is a potential source of systematic uncertainty in testing the no-hair theorem with stellar orbits.

We investigate the distribution of neutron star masses in different populations of binaries, employing Bayesian statistical techniques. In particular, we explore the differences in neutron star masses between sources that have experienced distinct evolutionary paths and accretion episodes. We find that the distribution of neutron star masses in non-recycled eclipsing high-mass binaries as well as of slow pulsars, which are all believed to be near their birth masses, has a mean of 1.28 M_solar and a dispersion of 0.24 M_solar. These values are consistent with expectations for neutron star formation in core-collapse supernovae. On the other hand, double neutron stars, which are also believed to be near their birth masses, have a much narrower mass distribution, peaking at 1.33 M_solar but with a dispersion of only 0.05 M_solar. Such a small dispersion cannot easily be understood and perhaps points to a particular and rare formation channel. The mass distribution of neutron stars that have been recycled has a mean of 1.48 M_solar and a dispersion of 0.2 M_solar, consistent with the expectation that they have experienced extended mass accretion episodes. The fact that only a very small fraction of recycled neutron stars in the inferred distribution have masses that exceed ~2 M_solar suggests that only a few of these neutron stars cross the mass threshold to form low mass black holes.

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.

Dynamical mass measurements to date have allowed determinations of the mass M and the distance D of a number of nearby supermassive black holes. In the case of Sgr A*, these measurements are limited by a strong correlation between the mass and distance scaling roughly as M ~ D^2. Future very-long baseline interferometric (VLBI) observations will image a bright and narrow ring surrounding the shadow of a supermassive black hole, if its accretion flow is optically thin. In this paper, we explore the prospects of reducing the correlation between mass and distance with the combination of dynamical measurements and VLBI imaging of the ring of Sgr A*. We estimate the signal to noise ratio of near-future VLBI arrays that consist of five to six stations, and we simulate measurements of the mass and distance of Sgr A* using the expected size of the ring image and existing stellar ephemerides. We demonstrate that, in this best-case scenario, VLBI observations at 1 mm can improve the error on the mass by a factor of about two compared to the results from the monitoring of stellar orbits alone. We identify the additional sources of uncertainty that such imaging observations have to take into account. In addition, we calculate the angular diameters of the bright rings of other nearby supermassive black holes and identify the optimal targets besides Sgr A* that could be imaged by a ground-based VLBI array or future space-VLBI missions allowing for refined mass measurements.

A moderately spinning neutron star acquires an oblate shape and a spacetime with a significant quadrupole moment. These two properties affect its apparent surface area for an observer at infinity, as well as the lightcurve arising from a hot spot on its surface. In this paper, we develop a ray-tracing algorithm to calculate the apparent surface areas of moderately spinning neutron stars making use of the Hartle-Thorne metric. This analytic metric allows us to calculate various observables of the neutron star in a way that depends only on its macroscopic properties and not on the details of its equation of state. We use this algorithm to calculate the changes in the apparent surface area, which could play a role in measurements of neutron star radii and, therefore, in constraining their equation of state. We show that whether the spinning neutron star appears larger or smaller than its non-rotating counterpart depends primarily on its equatorial radius. For neutron stars with radii ~10 km, the corrections to the Schwarzschild spacetime cause the apparent surface area to increase with spin frequency. In contrast, for neutron stars with radii ~15 km, the oblateness of the star dominates the spacetime corrections and causes the apparent surface area to decrease with increasing spin frequency. In all cases, the change in the apparent geometric surface area for the range of observed spin frequencies is < 5% and hence only a small source of error in the measurement of neutron star radii.

We develop a new numerical algorithm to model collisional cascades in debris disks. Because of the large dynamical range in particle masses, we solve the integro-differential equations describing erosive and catastrophic collisions in a particle-in-a-box approach, while treating the orbital dynamics of the particles in an approximate fashion. We employ a new scheme for describing erosive (cratering) collisions that yields a continuous set of outcomes as a function of colliding masses. We demonstrate the stability and convergence characteristics of our algorithm and compare it with other treatments. We show that incorporating the effects of erosive collisions results in a decay of the particle distribution that is significantly faster than with purely catastrophic collisions.

We explore the evolution of the mass distribution of dust in collision-dominated debris disks, using the collisional code introduced in our previous paper. We analyze the equilibrium distribution and its dependence on model parameters by evolving over 100 models to 10 Gyr. With our numerical models, we confirm that systems reach collisional equilibrium with a mass distribution that is steeper than the traditional solution by Dohnanyi (1969). Our model yields a quasi steady-state slope of n(m) ~ m^{-1.88} [n(a) ~ a^{-3.65}] as a robust solution for a wide range of possible model parameters. We also show that a simple power-law function can be an appropriate approximation for the mass distribution of particles in certain regimes. The steeper solution has observable effects in the submillimeter and millimeter wavelength regimes of the electromagnetic spectrum. We assemble data for nine debris disks that have been observed at these wavelengths and, using a simplified absorption efficiency model, show that the predicted slope of the particle mass distribution generates SEDs that are in agreement with the observed ones.

The observed flat rotation curves of galaxies require either the presence of dark matter in Newtonian gravitational potentials or a significant modification to the theory of gravity at galactic scales. Detecting relativistic Doppler shifts and gravitational effects in the rotation curves offers a tool for distinguishing between predictions of gravity theories that modify the inertia of particles and those that modify the field equations. These higher-order effects also allow us in principle, to test whether dark matter particles obey the equivalence principle. We calculate here the magnitudes of the relativistic Doppler and gravitational shifts expected in realistic models of galaxies in a general metric theory of gravity. We identify a number of observable quantities that measure independently the special- and general-relativistic effects in each galaxy and suggest that both effects might be detected in a statistical sense by combining appropriately the rotation curves of a large number of galaxies.

Time resolved X-ray spectroscopy of thermonuclear bursts observed from low mass X-ray binaries offer a unique tool to measure neutron star masses and radii. In this paper, we continue our systematic analysis of all the X-ray bursts observed with RXTE from X-ray binaries. We determine the events which show clear evidence for photospheric radius expansion and measure the Eddington limits for these accreting neutron stars using the bolometric fluxes attained at the touchdown moments of each X-ray burst. We employ a Bayesian technique to investigate the degree to which the Eddington limit for each source remains constant between bursts. We find that for sources with a large number of radius expansion bursts, systematic uncertainties are at a 5-10% level. Moreover, in six sources with only pairs of Eddington-limited bursts, the distribution of fluxes is consistent with a ~10% fractional dispersion. This indicates that the spectroscopic measurements of neutron star masses and radii using thermonuclear X-ray bursts can reach the level of accuracy required to distinguish between different neutron star equations of state, provided that uncertainties related to the overall flux calibration of X-ray detectors are of comparable magnitude.

The masses and radii of low-magnetic field neutron stars can be measured by combining different observable quantities obtained from their X-ray spectra during thermonuclear X-ray bursts. One of these quantities is the apparent radius of each neutron star as inferred from the X-ray flux and spectral temperature measured during the cooling tails of bursts, when the thermonuclear flash is believed to have engulfed the entire star. In this paper, we analyze 13,095 X-ray spectra of 446 X-ray bursts observed from 12 sources in order to assess possible systematic effects in the measurements of the apparent radii of neutron stars. We first show that the vast majority of the observed X-ray spectra are consistent with blackbody functions to within a few percent. We find that most X-ray bursts follow a very well determined correlation between X-ray flux and temperature, which is consistent with the whole neutron-star surface emitting uniformly during the cooling tails. We develop a Bayesian Gaussian mixture algorithm to measure the apparent radii of the neutron stars in these sources, while detecting and excluding a small number of X-ray bursts that show irregular cooling behavior. This algorithm also provides us with a quantitative measure of the systematic uncertainties in the measurements. We find that those errors in the spectroscopic determination of neutron-star radii that are introduced by systematic effects in the cooling tails of X-ray bursts are in the range $\simeq 3-8$%. Such small errors are adequate to distinguish between different equations of state provided that uncertainties in the distance to each source and the absolute calibration of X-ray detectors do not dominate the error budget.

Observations of stars and pulsars orbiting the black hole in the center of the Milky Way offer the potential of measuring not only the mass of the black hole but also its spin and quadrupole moment, thereby providing observational verification of the no-hair theorem. The relativistic effects that will allow us to measure these higher moments of the gravitational field, however, are very small and may be masked by drag forces that stars and pulsars experience orbiting within the hot, tenuous plasma that surrounds the black hole. The properties of this plasma at large distances from the central object have been measured using observations of the extended X-ray emission that surrounds the point source. At distances comparable to the black-hole event horizon, the properties of the accretion flow have been constrained using observations of its long-wavelength emission and polarization, as well as of the size of the emitting region at 1.3 mm. I use models of the plasma density and temperature at various distances from the black hole to investigate the effect of hydrodynamic drag forces on future measurements of the higher moments of its gravitational field. I find that hydrodynamic drag does not preclude measurements of the black hole spin and quadrupole moment using high-resolution observations of stars and pulsars that orbit within a few thousand gravitational radii from its horizon.

According to the no-hair theorem, astrophysical black holes are uniquely characterized by their masses and spins and are described by the Kerr metric. Several parametric deviations from the Kerr metric have been suggested to study observational signatures in both the electromagnetic and gravitational-wave spectra that differ from the expected Kerr signals. Due to the no-hair theorem, however, such spacetimes cannot be regular everywhere outside the event horizons, if they are solutions to the Einstein field equations; they are often characterized by naked singularities or closed time-like loops in the regions of the spacetime that are accessible to an external observer. For observational tests of the no-hair theorem that involve phenomena in the vicinity of the circular photon orbit or the innermost stable circular orbit around a black hole, these pathologies limit the applicability of the metrics only to compact objects that do not spin rapidly. In this paper, we construct a Kerr-like metric which depends on a set of free parameters in addition to its mass and spin and which is regular everywhere outside of the event horizon. We derive expressions for the energy and angular momentum of a particle on a circular equatorial orbit around the black hole and compute the locations of the innermost stable circular orbit and the circular photon orbit. We demonstrate that these orbits change significantly for even moderate deviations from the Kerr metric. The properties of our metric make it an ideally suited spacetime to carry out strong-field tests of the no-hair theorem in the electromagnetic spectrum using the properties of accretion flows around astrophysical black holes of arbitrary spin.

The black hole in the center of the Milky Way has been observed and modeled intensely during the last decades. It is also the prime target of a number of new experiments that aim to zoom into the vicinity of its horizon and reveal the inner working of its spacetime. In this review we discuss our current understanding of the gravitational field of Sgr A* and the prospects of testing the Kerr nature of its spacetime via imaging, astrometric, and timing observations.

We describe a new numerical algorithm for ray tracing in the external spacetimes of spinning compact objects characterized by arbitrary quadrupole moments. Such spacetimes describe non-Kerr vacuum solutions that can be used to test the no-hair theorem in conjunction with observations of accreting black holes. They are also appropriate for neutron stars with spin frequencies in the 300-600 Hz range, which are typical of the bursting sources in low-mass X-ray binaries. We use our algorithm to show that allowing for the quadrupole moment of the spacetime to take arbitrary values leads to observable effects in the profiles of relativistic broadened fluorescent iron lines from geometrically thin accretion disks.

We use dynamical mass measurements of 16 black holes in transient low-mass X-ray binaries to infer the stellar black hole mass distribution in the parent population. We find that the observations are best described by a narrow mass distribution at 7.8 +/- 1.2 Msolar. We identify a selection effect related to the choice of targets for optical follow-ups that results in a flux-limited sample. We demonstrate, however, that this selection effect does not introduce a bias in the observed distribution and cannot explain the absence of black holes in the 2-5 solar mass range. On the high mass end, we argue that the rapid decline in the inferred distribution may be the result of the particular evolutionary channel followed by low-mass X-ray binaries. This is consistent with the presence of high-mass black holes in the persistent, high-mass X-ray binary sources. If the paucity of low-mass black holes is caused by a sudden decrease of the supernova explosion energy with increasing progenitor mass, this would have observable implications for ongoing transient surveys that target core-collapse supernovae. Our results also have significant implications for the calculation of event rates from the coalescence of black hole binaries for gravitational wave detectors.

According to the no-hair theorem, astrophysical black holes are uniquely described by their masses and spins. An observational test of the no-hair theorem can be performed by measuring at least three different multipole moments of the spacetime of a black hole and verifying whether their values are consistent with the unique combinations of the Kerr solution. In this paper, we study quasi-periodic variability observed in the emission from black holes across the electromagnetic spectrum as a test of the no-hair theorem. We derive expressions for the Keplerian and epicyclic frequencies in a quasi-Kerr spacetime, in which the quadrupole moment is a free parameter in addition to mass and spin. We show that, for moderate spins, the Keplerian frequency is practically independent of small deviations of the quadrupole moment from the Kerr value, while the epicyclic frequencies exhibit significant variations. We apply this framework to quasi-periodic oscillations in black-hole X-ray binaries in two different scenarios. In the case that a pair of quasi-periodic oscillations can be identified as the fundamental g- and c-modes in the accretion disk, we show that the no-hair theorem can be tested in conjunction with an independent mass measurement. If, on the other hand, the pairs of oscillations are identified with non-parametric resonance of dynamical frequencies in the accretion disk, then testing the no-hair theorem also requires an independent measurement of the black-hole spin. In addition, we argue that VLBI observations of Sgr A* may test the no-hair theorem through a combination of imaging observations and the detection of quasi-periodic variability.

The recent measurement of the Shapiro delay in the radio pulsar PSR J1614-2230 yielded a mass of 1.97 +/- 0.04 M_sun, making it the most massive pulsar known to date. Its mass is high enough that, even without an accompanying measurement of the stellar radius, it has a strong impact on our understanding of nuclear matter, gamma-ray bursts, and the generation of gravitational waves from coalescing neutron stars. This single high mass value indicates that a transition to quark matter in neutron-star cores can occur at densities comparable to the nuclear saturation density only if the quarks are strongly interacting and are color superconducting. We further show that a high maximum neutron-star mass is required if short duration gamma-ray bursts are powered by coalescing neutron stars and, therefore, this mechanism becomes viable in the light of the recent measurement. Finally, we argue that the low-frequency (<= 500 Hz) gravitational waves emitted during the final stages of neutron-star coalescence encode the properties of the equation of state because neutron stars consistent with this measurement cannot be centrally condensed. This will facilitate the measurement of the neutron star equation of state with Advanced LIGO/Virgo.

According to the no-hair theorem, astrophysical black holes are uniquely described by their mass and spin. In this paper, we review a new framework for testing the no-hair hypothesis with observations in the electromagnetic spectrum. The approach is formulated in terms of a Kerr-like spacetime containing a quadrupole moment that is independent of both mass and spin. If the no-hair theorem is correct, then any deviation from the Kerr metric quadrupole has to be zero. We show how upcoming VLBI imaging observations of Sgr A* as well as spectroscopic observations of iron lines from accreting black holes with IXO may lead to the first astrophysical test of the no-hair theorem.

We study the structure of neutron stars in f(R) gravity theories with perturbative constraints. We derive the modified Tolman-Oppenheimer-Volkov equations and solve them for a polytropic equation of state. We investigate the resulting modifications to the masses and radii of neutron stars and show that observations of surface phenomena alone cannot break the degeneracy between altering the theory of gravity versus choosing a different equation of state of neutron-star matter. On the other hand, observations of neutron-star cooling, which depends on the density of matter at the stellar interior, can place significant constraints on the parameters of the theory.

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 use the recent measurement of the velocity dispersion of star-forming, outer-disk knots by Herbert-Fort et al. in the nearly face-on galaxy NGC 628, in combination with other data from the literature, to execute a straightforward test of gravity at low accelerations. Specifically, the rotation curve at large radius sets the degree of non-standard acceleration and then the predicted scaleheight of the knots at that radius provides the test of the scenario. For our demonstration, we presume that the H alpha knots, which are young (age < 10 Myr), are distributed like the gas from which they have recently formed and find a marginal (> 97% confidence) discrepancy with a modified gravity scenario given the current data. More interestingly, we demonstrate that there is no inherent limitation that prevents such a test from reaching possible discrimination at the > 4 sigma level with a reasonable investment of observational resources.

According to the no-hair theorem, all astrophysical black holes are fully described by their masses and spins. This theorem can be tested observationally by measuring (at least) three different multipole moments of the spacetimes of black holes. In this paper, we analyze images of black holes within a framework that allows us to calculate observables in the electromagnetic spectrum as a function of the mass, spin, and, independently, the quadrupole moment of a black hole. We show that a deviation of the quadrupole moment from the expected Kerr value leads to images of black holes that are either prolate or oblate depending on the sign and magnitude of the deviation. In addition, there is a ring-like structure around the black-hole shadow with a diameter of about 10 black-hole masses that is substantially brighter than the image of the underlying accretion flow and that is independent of the astrophysical details of accretion flow models. We show that the shape of this ring depends directly on the mass, spin, and quadrupole moment of the black hole and can be used for an independent measurement of all three parameters. In particular, we demonstrate that this ring is highly circular for a Kerr black hole with a spin a<0.9M, independent of the observer's inclination, but becomes elliptical and asymmetric if the no-hair theorem is violated. Near-future very-long baseline interferometric observations of Sgr A* will image this ring and may allow for an observational test of the no-hair theorem.

According to the no-hair theorem, an astrophysical black hole is uniquely described by only two quantities, the mass and the spin. In this series of papers, we investigate a framework for testing the no-hair theorem with observations of black holes in the electromagnetic spectrum. We formulate our approach in terms of a parametric spacetime which contains a quadrupole moment that is independent of both mass and spin. If the no-hair theorem is correct, then any deviation of the black-hole quadrupole moment from its Kerr value has to be zero. We analyze in detail the properties of this quasi-Kerr spacetime that are critical to interpreting observations of black holes and demonstrate their dependence on the spin and quadrupole moment. In particular, we show that the location of the innermost stable circular orbit and the gravitational lensing experienced by photons are affected significantly at even modest deviations of the quadrupole moment from the value predicted by the no-hair theorem. We argue that observations of black-hole images, of relativistically broadened iron lines, as well as of thermal X-ray spectra from accreting black holes will lead in the near future to an experimental test of the no-hair theorem.

Neutron stars are sensitive laboratories for testing general relativity, especially when considering deviations where velocities are relativistic and gravitational fields are strong. One such deviation is described by dynamical, Chern-Simons modified gravity, where the Einstein-Hilbert action is modified through the addition of the gravitational parity-violating Pontryagin density coupled to a field. This four-dimensional effective theory arises naturally both in perturbative and non-perturbative string theory, loop quantum gravity, and generic effective field theory expansions. We calculate here Chern-Simons modifications to the properties and gravitational fields of slowly spinning neutron stars. We find that the Chern-Simons correction affects only the gravitomagnetic sector of the metric to leading order, thus introducing modifications to the moment of inertia but not to the mass-radius relation. We show that an observational determination of the moment of inertia to an accuracy of 10%, as is expected from near-future observations of the double pulsar, will place a constraint on the Chern-Simons coupling constant of \xi^{1/4} < 5 km, which is at least three-orders of magnitude stronger than the previous strongest bound.

The properties of matter at ultra-high densities, low temperatures, and with a significant asymmetry between protons and neutrons can be studied exclusively through astrophysical observations of neutron stars. We show that measurements of the masses and radii of neutron stars can lead to tight constraints on the pressure of matter at three fiducial densities, from 1.85 to 7.4 times the density of nuclear saturation, in a manner that is largely model-independent and that captures the key characteristics of the equation of state. We demonstrate that observations with 10% uncertainties of at least three neutron stars can lead to measurements of the pressure at these fiducial densities with an accuracy of 0.11 dex or ~ 30%. Observations of three neutron stars with 5% uncertainties are sufficient to distinguish at a better than 3-sigma confidence level between currently proposed equations of state. In the electromagnetic spectrum, such accurate measurements will become possible for weakly-magnetic neutron stars during thermonuclear flashes and in quiescence with future missions such as the International X-ray Observatory (IXO).

The Kerr spacetime of spinning black holes is one of the most intriguing predictions of Einstein's theory of general relativity. The special role this spacetime plays in the theory of gravity is encapsulated in the no-hair theorem, which states that the Kerr metric is the only realistic black-hole solution of the vacuum field equations. Recent and anticipated advances in the observations of black holes throughout the electromagnetic spectrum have secured our understanding of their basic properties while opening up new opportunities for devising tests of the Kerr metric. In this paper, we argue that imaging and spectroscopic observations of accreting black-holes with current and future instruments can lead to the first direct test of the no-hair theorem.

Observations of compact objects in the electromagnetic spectrum and the detection of gravitational waves from them can lead to quantitative tests of the theory of general relativity in the strong-field regime following two very different approaches. In the first approach, the general relativistic field equations are modified at a fundamental level and the magnitudes of the potential deviations are constrained by comparison with observations. In the second approach, the exterior spacetimes of compact objects are parametrized in a phenomenological way, the various parameters are measured observationally, and the results are finally compared against the general relativistic predictions. In this article, I discuss the current status of both approaches, focusing on the lessons learned from a large number of recent investigations.

Large extra dimensions have been proposed as a possible solution to the hierarchy problem in physics. One of the suggested models, the RS2 braneworld model, makes a prediction that black holes evaporate by Hawking radiation on a short timescale that depends on the black hole mass and on the asymptotic radius of curvature of the extra dimensions. Thus the size of the extra dimensions can be constrained by astrophysical observations. Here we point out that the black hole, recently discovered in a globular cluster in galaxy NGC 4472, places the strongest constraint on the maximum size of the extra dimensions, L < 0.003 mm. This black hole has the virtues of old age and relatively small mass. The derived upper limit is within an order of magnitude of the absolute limit afforded by astrophysical observations of black holes.

Major observational efforts in the coming decade are designed to probe the equation of state of dark energy. Measuring a deviation of the equation-of-state parameter w from -1 would indicate a dark energy that cannot be represented solely by a cosmological constant. While it is commonly assumed that any implied modification to the LambdaCDM model amounts to the addition of new dynamical fields, we propose here a framework for investigating whether or not such new fields are required when cosmological observations are combined with a set of minimal assumptions about the nature of gravitational physics. In our approach, we treat the additional degrees of freedom as perturbatively constrained and calculate a number of observable quantities, such as the Hubble expansion rate and the cosmic acceleration, for a homogeneous Universe. We show that current observations place our Universe within the perturbative validity of our framework and allow for the presence of non-dynamical gravitational degrees of freedom at cosmological scales.

Modifications to gravity that add additional functions of the Ricci curvature to the Einstein-Hilbert action -- collectively known as $f(R)$ theories -- have been studied in great detail. When considered as complete theories of gravity they can generate non-perturbative deviations from the general relativistic predictions in the solar system, and the simplest models show instabilites on cosmological scales. Here we show that it is possible to treat $f(R)=R\pm\mu^4/R$ gravity in a perturbative fashion such that it shows no instabilities on cosmological scales and, in the solar system, is consistent with measurements of the PPN parameters. We show that such a theory produces a spatially flat, accelerating universe, even in the absence of dark energy and when the matter density is too small to close the universe in the general relativistic case.

One of the plausible unification schemes in physics considers the observable universe to be a 4-dimensional surface (the "brane") embedded in a higher-dimensional curved spacetime (the "bulk"). In such braneworld gravity models with infinitely large extra dimensions, black holes evaporate fast through the emission of the additional gravitational degrees of freedom, resulting in lifetimes of stellar-mass black holes that are significantly smaller than the Hubble time. We show that the predicted evaporation rate leads to a change in the orbital period of X-ray binaries harboring black holes that is observable with current instruments. We obtain an upper limit on the rate of change of the orbital period of the binary A0620-00 and use it to constrain the asymptotic curvature radius of the extra dimension to a value comparable to the one obtained by table-top experiments. Furthermore we argue that any measurement of a period increase for low-mass X-ray binaries with a high mass ratio is evidence for new physics beyond general relativity and the standard model.

Bursting X-ray binaries in globular clusters are ideal sources for measuring neutron star masses and radii, and hence, for determining the equation of state of cold, ultradense matter. We use time-resolved spectroscopic data from EXO 1745-248 during thermonuclear bursts that show strong evidence for photospheric radius expansion to measure the Eddington flux and the apparent surface area of the neutron star. We combine this with the recent measurement of the distance to the globular cluster Terzan 5, where this source resides, to measure the neutron star mass and radius. We find tightly constrained pairs of values for the mass and radius, which are centered around M=1.4 M_sun and R=11 km or around M=1.7 M_sun and R=9 km. These values favor nucleonic equations of state with symmetry energy that is relatively low and has a weak dependence on density.

The redshift evolution of the Tully-Fisher Relation probes gravitational dynamics that must be consistent with any modified gravity theory seeking to explain the galactic rotation curves without the need for dark matter. Within the context of non-relativistic Modified Newtonian Dynamics (MOND), the characteristic acceleration scale of the theory appears to be related to the current value of either the Hubble constant, i.e., alpha ~ cH_0, or the dark energy density, i.e., alpha (8 pi G rho_lambda/3)^{1/2}. If these relations are the manifestation of a fundamental coupling of a_0 to either of the two cosmological parameters, the cosmological evolution would then dictate a particular dependence of the MOND acceleration scale with redshift that can be tested with Tully-Fisher relations of high-redshift galaxies. We compare this prediction to two sets of Tully-Fisher data with redshifts up to z=1.2. We find that both couplings are excluded within the formal uncertainties. However, when we take into account the potential systematic uncertainties in the data, we find that they marginally favor the coupling of the MOND acceleration scale to the density of dark energy.

The observational characteristics of quasi-periodic oscillations (QPOs) from accreting neutron stars strongly indicate the oscillatory modes in the innermost regions of accretion disks as a likely source of the QPOs. The inner regions of accretion disks around neutron stars can harbor very high frequency modes related to the radial epicyclic frequency $\kappa $. The degeneracy of $\kappa $ with the orbital frequency $\Omega $ is removed in a non-Keplerian boundary or transition zone near the magnetopause between the disk and the compact object. We show, by analyzing the global hydrodynamic modes of long wavelength in the boundary layers of viscous accretion disks, that the fastest growing mode frequencies are associated with frequency bands around $\kappa $ and $\kappa \pm \Omega $. The maximum growth rates are achieved near the radius where the orbital frequency $\Omega $ is maximum. The global hydrodynamic parameters such as the surface density profile and the radial drift velocity determine which modes of free oscillations will grow at a given particular radius in the boundary layer. In accordance with the peak separation between kHz QPOs observed in neutron-star sources, the difference frequency between two consecutive bands of the fastest growing modes is always related to the spin frequency of the neutron star. This is a natural outcome of the boundary condition imposed by the rotating magnetosphere on the boundary region of the inner disk.

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.

Neutron stars and black holes are the astrophysical systems with the strongest gravitational fields in the universe. In this article, I review the prospect of probing with observations of such compact objects some of the most intriguing General Relativistic predictions in the strong-field regime: the absence of stable circular orbits near a compact object and the presence of event horizons around black-hole singularities. I discuss the need for a theoretical framework within which future experiments will provide detailed, quantitative tests of gravity theories. Finally, I summarize the constraints imposed by current observations of neutron stars on potential deviations from General Relativity.

The detection of fast quasi-periodic variability from accreting black holes and neutron stars has been used to constrain their masses, radii, and spins. If the observed oscillations are linear modes in the accretion disks, then bounds can be placed on the properties of the central objects by assuming that these modes are locally sub-Keplerian. If, on the other hand, the observed oscillations correspond to non-linear resonances between disk modes, then the properties of the central objects can be measured by assuming that the resonant modes are excited at the same radial annulus in the disk. In this paper, we use numerical simulations of vertically integrated, axisymmetric hydrodynamic accretion disks to provide examples of situations in which the assumptions implicit in both methods are not satisfied. We then discuss our results for the robustness of the mass and spin measurements of compact objects from variability studies.

Eddington-limited X-ray bursts from neutron stars can be used in conjunction with other spectroscopic observations to measure neutron star masses, radii, and distances. In order to quantify some of the uncertainties in the determination of the Eddington limit, we analysed a large sample of photospheric radius-expansion thermonuclear bursts observed with the Rossi X-ray Timing Explorer. We identified the instant at which the expanded photosphere "touches down" back onto the surface of the neutron star and compared the corresponding touchdown flux to the peak flux of each burst. We found that for the majority of sources, the ratio of these fluxes is smaller than 1.6, which is the maximum value expected from the changing gravitational redshift during the radius expansion episodes (for a 2M_sun neutron star). The only sources for which this ratio is larger than 1.6 are high inclination sources that include dippers and Cyg X-2. We discuss two possible geometric interpretations of this effect and show that the inferred masses and radii of neutron stars are not affected by this bias. On the other hand, systematic uncertainties as large as ~50% may be introduced to the distance determination.

Numerical simulations of turbulent, magnetized, differentially rotating flows driven by the magnetorotational instability are often used to calculate the effective values of alpha viscosity that is invoked in analytical models of accretion discs. In this paper we use various dynamical models of turbulent magnetohydrodynamic stresses, as well as numerical simulations of shearing boxes, to show that angular momentum transport in MRI-driven accretion discs cannot be described by the standard model for shear viscosity. In particular, we demonstrate that turbulent magnetorotational stresses are not linearly proportional to the local shear and vanish identically for angular velocity profiles that increase outwards.

I show that several observable properties of bursting neutron stars in metric theories of gravity can be calculated using only conservation laws, Killing symmetries, and the Einstein equivalence principle, without requiring the validity of the general relativistic field equations. I calculate, in particular, the gravitational redshift of a surface atomic line, the touchdown luminosity of a radius-expansion burst, which is believed to be equal to the Eddington critical luminosity, and the apparent surface area of a neutron star as measured during the cooling tails of bursts. I show that, for a general metric theory of gravity, the apparent surface area of a neutron star depends on the coordinate radius of the stellar surface and on its gravitational redshift in the exact same way as in general relativity. On the other hand, the Eddington critical luminosity depends also on an additional parameter that measures the degree to which the general relativistic field equations are satisfied. These results can be used in conjunction with current and future high-energy observations of bursting neutron stars to test general relativity in the strong-field regime.

Considerable attention has recently focused on gravity theories obtained by extending general relativity with additional scalar, vector, or tensor degrees of freedom. In this paper, we show that the black-hole solutions of these theories are essentially indistinguishable from those of general relativity. Thus, we conclude that a potential observational verification of the Kerr metric around an astrophysical black hole cannot, in and of itself, be used to distinguish between these theories. On the other hand, it remains true that detection of deviations from the Kerr metric will signify the need for a major change in our understanding of gravitational physics.

We present a scaling law that predicts the values of the stresses obtained in numerical simulations of saturated MRI-driven turbulence in non-stratified shearing boxes. It relates the turbulent stresses to the strength of the vertical magnetic field, the sound speed, the vertical size of the box, and the numerical resolution and predicts accurately the results of 35 numerical simulations performed for a wide variety of physical conditions. We use our result to show that the saturated stresses in simulations with zero net magnetic flux depend linearly on the numerical resolution and would become negligible if the resolution were set equal to the natural dissipation scale in astrophysical disks. We conclude that, in order for MRI-driven turbulent angular momentum transport to be able to account for the large value of the effective alpha viscosity inferred observationally, the disk must be threaded by a significant vertical magnetic field and the turbulent magnetic energy must be in near equipartition with the thermal energy. This result has important implications for the spectra of accretion disks and their stability.

The low-luminosity black hole Sgr A* provides a testbed for models of Radiatively Inefficient Accretion Flows (RIAFs). Recent sub-millimeter linear polarization measurements of Sgr A* have provided evidence that the electrons in the accretion flow are relativistic over a large range of radii. Here, we show that these high temperatures result in elliptical plasma normal modes. Thus, polarized millimeter and sub-millimeter radiation emitted within RIAFs will undergo generalized Faraday rotation, a cyclic conversion between linear and circular polarization. This effect will not depolarize the radiation even if the rotation measure is extremely high. Rather, the beam will take on the linear and circular polarization properties of the plasma normal modes. As a result, polarization measurements of Sgr A* in this frequency regime will constrain the temperature, density and magnetic profiles of RIAF models.

Accretion disks are three-dimensional, turbulent, often self-gravitating, magnetohydrodynamic flows, which can be modeled in detail with numerical simulations. In this paper, we present a new algorithm that is based on a spectral decomposition method to simulate such flows. Because of the high order of the method, we can solve the induction equation in terms of the magnetic potential and, therefore, ensure trivially that the magnetic fields in the numerical solution are divergence free. The spectral method also suffers minimally from numerical dissipation and allows for an easy implementation of models for sub-grid physics. Both properties make our method ideal for studying MHD turbulent flows such as those found in accretion disks around compact objects. We verify our algorithm with a series of standard tests and use it to show the development of MHD turbulnce in a simulation of an accretion disk. Finally, we study the evolution and saturation of the power spectrum of MHD turbulence driven by the magnetorotational instability.

In braneworld gravity models with a finite AdS curvature in the extra dimension, the AdS/CFT correspondence leads to a prediction for the lifetime of astrophysical black holes that is significantly smaller than the Hubble time, for asymptotic curvatures that are consistent with current experiments. Using the recent measurements of the position, three-dimensional spatial velocity, and mass of the black hole XTE J1118+480, I calculate a lower limit on its kinematic age of 11 Myr (95% confidence). This translates into an upper limit for the asymptotic AdS curvature in the extra dimensions of 0.08 mm, which significantly improves the limit obtained by table-top experiments of sub-mm gravity.

High resolution observations of Sgr A* have revealed a wide variety of phenomena, ranging from intense rapid flares to quasi-periodic oscillations, making this object an ideal system to study the properties of low luminosity accreting black holes. In this paper, we use a pseudo-spectral algorithm to construct and evolve a three-dimensional magnetohydrodynamic model of the accretion disk in Sgr A*. Assuming a hybrid thermal-nonthermal emission scheme, we show that the MHD turbulence can by itself only produce factor of two fluctuations in luminosity. These amplitudes in variation cannot explain the magnitude of flares observed in this system. However, we also demonstrate that density perturbations in the disk do produce outbursts qualitatively similar to those observed by XMM-Newton in X-rays and ground-based facilities in the near infrared. Quasi-periodic oscillations emerge naturally in the simulated lightcurves. We attribute these to non-axisymmetric density perturbations that emerge as the disk evolves back toward its quiescent state.

We develop a local model for the exponential growth and saturation of the Reynolds and Maxwell stresses in turbulent flows driven by the magnetorotational instability. We first derive equations that describe the effects of the instability on the growth and pumping of the stresses. We highlight the relevance of a new type of correlations that couples the dynamical evolution of the Reynolds and Maxwell stresses and plays a key role in developing and sustaining the magnetorotational turbulence. We then supplement these equations with a phenomenological description of the triple correlations that lead to a saturated turbulent state. We show that the steady-state limit of the model describes successfully the correlations among stresses found in numerical simulations of shearing boxes.

The magnetorotational instability is thought to be responsible for the generation of magnetohydrodynamic turbulence that leads to enhanced outward angular momentum transport in accretion discs. Here, we present the first formal analytical proof showing that, during the exponential growth of the instability, the mean (averaged over the disc scale-height) Reynolds stress is always positive, the mean Maxwell stress is always negative, and hence the mean total stress is positive and leads to a net outward flux of angular momentum. More importantly, we show that the ratio of the Maxwell to the Reynolds stresses during the late times of the exponential growth of the instability is determined only by the local shear and does not depend on the initial spectrum of perturbations or the strength of the seed magnetic. Even though we derived these properties of the stress tensors for the exponential growth of the instability in incompressible flows, numerical simulations of shearing boxes show that this characteristic is qualitatively preserved under more general conditions, even during the saturated turbulent state generated by the instability.

Spectral methods are well suited for solving hydrodynamic problems in which the self-gravity of the flow needs to be considered. Because Poisson's equation is linear, the numerical solution for the gravitational potential for each individual mode of the density can be pre-computed, thus reducing substantially the computational cost of the method. In this second paper, we describe two different approaches to computing the gravitational field of a two-dimensional flow with pseudo-spectral methods. For situations in which the density profile is independent of the third coordinate (i.e., an infinite cylinder), we use a standard Poisson solver in spectral space. On the other hand, for situations in which the density profile is a delta function along the third coordinate (i.e., an infinitesimally thin disk), or any other function known a priori, we perform a direct integration of Poisson's equation using a Green's functions approach. We devise a number of test problems to verify the implementations of these two methods. Finally, we use our method to study the stability of polytropic, self-gravitating disks. We find that, when the polytropic index Gamma is <= 4/3, Toomre's criterion correctly describes the stability of the disk. However, when Gamma > 4/3 and for large values of the polytropic constant K, the numerical solutions are always stable, even when the linear criterion predicts the contrary. We show that, in the latter case, the minimum wavelength of the unstable modes is larger than the extent of the unstable region and hence the local linear analysis is inapplicable.

In the inner regions of accretion disks around compact objects, the orbital frequency of the gas deviates from the local Keplerian value. For long-wavelength modes in this region, the radial epicyclic frequency kappa is higher than the azimuthal frequency Omega. This has significant implications for models of the twin kHz QPOs observed in many neutron-star sources that traditionally identify the frequencies of the two kHz QPOs with dynamical frequencies in the accretion disk. The recognition that the highest frequency in the transition or boundary region of the disk is actually the epicyclic frequency also modifies significantly the constraints imposed by the observation of high-frequency QPOs on the mass and radius of the compact objects.

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)

During the last decade it has become evident that the magnetorotational instability is at the heart of the enhanced angular momentum transport in weakly magnetized accretion disks around neutron stars and black holes. In this paper, we investigate the local linear stability of differentially rotating, magnetized flows and the evolution of the magnetorotational instability beyond the weak-field limit. We show that, when superthermal toroidal fields are considered, the effects of both compressibility and magnetic tension forces, which are related to the curvature of toroidal field lines, should be taken fully into account. We demonstrate that the presence of a strong toroidal component in the magnetic field plays a non-trivial role. When strong fields are considered, the strength of the toroidal magnetic field not only modifies the growth rates of the unstable modes but also determines which modes are subject to instabilities. We find that, for rotating configurations with Keplerian laws, the magnetorotational instability is stabilized at low wavenumbers for toroidal Alfven speeds exceeding the geometric mean of the sound speed and the rotational speed. We discuss the significance of our findings for the stability of cold, magnetically dominated, rotating fluids and argue that, for these systems, the curvature of toroidal field lines cannot be neglected even when short wavelength perturbations are considered. We also comment on the implications of our results for the validity of shearing box simulations in which superthermal toroidal fields are generated.

We present a numerical method for studying the normal modes of accretion flows around black holes. In this first paper, we focus on two-dimensional, viscous, hydrodynamic disks, for which the linear modes have been calculated analytically in previous investigations. We use pseudo-spectral methods and low storage Runge-Kutta methods to solve the continuity equation, the Navier-Stokes equation, and the energy equation. We devise a number of test problems to verify the implementation. These tests demonstrate the ability of spectral methods to handle accurately advection problems and to reproduce correctly the stability criteria for differentially rotating hydrodynamic flows. They also show that our implementation is able to handle sound wave correctly with non-reflective boundary conditions, to recover the standard solution for a viscous spreading ring, and produce correctly the Shakura-Sunyaev steady disk solution. Finally, we have applied our algorithm to the problem of a non-axisymmetric viscous spreading ring and verify that such configuration is unstable to non-antisymmetric perturbations.

The five accreting, millisecond X-ray pulsars in ultracompact binaries that were recently discovered with the Rossi X-ray Timing Explorer provide excellent candidates for constraining the deviations from general relativity described by the Brans-Dicke parameter omega_BD. I calculate the expected rate of change of the orbital periods of these binaries and discuss the prospects of constraining omega_BD to values that are an order of magnitude larger than current constraints. Finally, I show how measurements of the orbital period derivative in ultracompact binaries can be used to place lower bounds on their orbital inclination.

This is a review of recent observational and theoretical developments related to accreting neutron stars and black holes in binary systems.

we propose a jet model for the low/hard state of black-hole X-ray sources which explains a) the X-ray spectra, b) the timelag spectra, c) the increase in the amplitude (QPO and high frequency) with increasing photon energy, and d) the narrowing of the autocorrelation function with increasing photon energy. The model (in its simplest form) assumes that i) there is a uniform magnetic field along the axis of the jet, ii) the electron density in the jet is inversely proportional to distance and iii) the jet is "hotter" near its center than at its periphery. We have performed Monte Carlo simulations of Compton upscattering of soft photons from the accretion disk and have found power-law high-energy spectra with photon number index in the range 1.5-2, power-law timelags versus Fourier frequency with index ~0.8, and an increase of the rms amplitude of the variability and a narrowing of the autocorrelation function with photon energy as they have been observed in Cygnus X-1.

The frequencies of quasi-periodic oscillations around neutron stars are believed to be related to characteristic frequencies in the gravitational fields of the compact objects. In different variability models, these include the Keplerian, epicyclic, and Lense-Thirring frequencies, which depend mostly on the properties of the stellar spacetimes. We argue that quasi-periodic oscillations in the X-ray flux of neutron stars can be used to map the external spacetimes of the compact objects and, therefore, lead to direct tests of general relativity in the strong-field regime. In particular, we show that particular extensions of General Relativity, in which the gravitational force felt by matter is mediated by both a rank-two tensor and a scalar field, can be constrained by current observations.

Fast variability studies of accreting black holes in the Galaxy offer us a unique opportunity to measure the spins of black holes and test the strong-field behavior of general relativity. In this review, I summarize the arguments often used in attempts of measuring the spins of black holes, concentrating on their theoretical foundations. I also argue that X-ray studies of accreting black holes will be able to provide in the future strong constraints on deviations from general relativity in the strong-field regime.

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.

In contrast to gravity in the weak-field regime, which has been subject to numerous experimental tests, gravity in the strong-field regime is largely unconstrained by observations. We show that gravity theories that cannot be rejected by solar system tests but that diverge from general relativity in the strong-field regime predict neutron stars with significantly different properties than their general relativistic counterparts. In particular, the range of redshifts of surface atomic lines predicted by such gravity theories is significantly larger than the uncertainty introduced by our lack of knowledge of the equation of state of ultra-dense matter. Measurements of such redshifted lines with current X-ray observatories such as Chandra and XMM-Newton can thus provide interesting new constraints on strong-field gravity.

We investigated the relation between the centroid frequency of the quasi-periodic oscillation observed in the power density spectra of a sample of galactic black-hole candidates with the power-law photon index obtained from spectral fits. Our aim is to avoid inner accretion disk radius determination directly from spectral fits, given the uncertainties of the absolute values obtained in that way, but to base our analysis on the likely association of QPO frequency to a characteristic radius. We used archival RXTE data of GRS 1915+105 and published parameters for GRO 1655-40, XTE J1550-564, XTE J1748-288 and 4U 1630-47. While for low values of the QPO frequency, the two parameters are clearly correlated for each source, there is evidence for a turnoff in the correlation above a characteristic frequency, different for different sources. We discuss the possible nature of this turnoff.

The line profiles from rotating neutron stars are affected by a number of relativistic processes such as Doppler boosts, strong self-lensing, frame-dragging, and the differential gravitational redshift arising from the stellar oblateness. In this Letter, we calculate line profiles taking into account the first two effects, which is accurate for rotation rates less than the breakup frequency. We show that the line profiles are not only broadened and weakened but are also significantly asymmetric, and allow for an independent measurement of both the mass and the radius of the neutron star. Furthermore, we investigate the case when a fraction of the neutron star surface contributes to the emission and find that the line profiles are typically doubly peaked. We discuss the implications of our results for searches for line features in the spectra of isolated neutron stars and X-ray bursters. We finally assess the systematic uncertainties introduced by the line asymmetry in inferring the compactness of neutron stars from the detection of redshifted lines.

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)

The illumination of an accretion disk around a black hole or neutron star by the central compact object or the disk itself often determines its spectrum, stability, and dynamics. The transport of radiation within the disk is in general a multi-dimensional, non-axisymmetric problem, which is challenging to solve. Here, I present a method of decomposing the radiative transfer equation that describes absorption, emission, and Compton scattering in an obliquely illuminated disk into a set of four one-dimensional transfer equations. I show that the exact calculation of the ionization balance and radiation heating of the accretion disk requires the solution of only one of the one-dimensional equations, which can be solved using existing numerical methods. I present a variant of the Feautrier method for solving the full set of equations, which accounts for the fact that the scattering kernels in the individual transfer equations are not forward-backward symmetric. I then apply this method in calculating the albedo of a cold, geometrically thin accretion disk.

We study an empirical model for a unified description of the power spectra of accreting neutron stars and black holes. This description is based on a superposition of multiple Lorentzians and offers the advantage that all QPO and noise components are dealt with in the same way, without the need of deciding in advance the nature of each component. This approach also allows us to compare frequencies of features with high and low coherences in a consistent manner and greatly facilitates comparison of power spectra across a wide range of source types and states. We apply the model to six sources, the low-luminosity X-ray bursters 1E 1724-3045, SLX 1735-269 and GS 1826-24, the high-latitude transient XTE J1118+480, the bright system Cir X-1, and the Z source GX 17+2. We find that it provides a good description of the observed spectra, without the need for a scale-free (1/f) component. We update previously reported correlations between characteristic frequencies of timing features in the light of this new approach and discuss similarities between different types of systems which may point towards similar underlying physics.

The spin periods of the ten observed Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters lie in the very narrow range 6-12 s. We use a point likelihood technique to assess the constraints this clustering imposes on the birth period and on the final period of such systems. We consider a general law for their spin evolution described by a constant braking index. We find that, for positive values of the braking index, the observed clustering requires an upper cut-off period that is very close to the maximum observed period of about 12 s. We also show that the constraint on the birth period depends very strongly on the assumed value of the braking index n, ranging from a few milliseconds for n>2 to a few seconds for n<2. We discuss the possible ways of tightening these constraints based on similarities with the population of radio pulsars and with future observations of such sources with current X-ray observatories.

We present the results of a multi-Lorentzian fit to the power spectra of two kilohertz QPO sources; 4U 0614+09 and 4U 1728-34. This work was triggered by recent results of a similar fit to the black-hole candidates (BHCs) GX 339-4 and Cyg X-1 by Nowak in 2000. We find that one to six Lorentzians are needed to fit the power spectra of our two sources. The use of exactly the same fit function reveals that the timing behaviour of 4U 0614+09 and 4U 1728-34 is almost identical at luminosities which are about a factor 5 different. As the characteristic frequency of the Lorentzians we use the frequency, nu_max, at which each component contributes most of its variance per log frequency as proposed by Belloni, Psaltis & van der Klis in 2001. When using nu_max instead of the centroid frequency of the Lorentzian, the recently discovered hectohertz Lorentzian is practically constant in frequency. We use our results to test the suggestions by, respectively, Psaltis Belloni and van der Klis in 1999 and Nowak in 2000 that the two Lorentzians describing the high-frequency end of the broad-band noise in BHCs in the low state can be identified with the kilohertz QPOs in the neutron star low mass X-ray binaries. We find, that when the two kilohertz QPOs are clearly present, the low-frequency part of the power spectrum is too complicated to draw immediate conclusions from the nature of the components detected in any one power spectrum. However, the relations we observe between the characteristic frequencies of the kilohertz QPOs and the band-limited noise, when compared to the corresponding relations in BHCs, hint towards the identification of the second-highest frequency Lorentzian in the BHCs with the lower kilohertz QPO. They do not confirm the identification of the highest-frequency Lorentzian with the upper kilohertz QPO.

Thermal emission from the surface of an ultramagnetic neutron star is believed to contribute significantly to the soft X-ray flux of the Anomalous X-ray Pulsars. We compare the detailed predictions of models of the surface emission from a magnetar to the observed properties of AXPs. In particular, we focus on the combination of their luminosities and energy-dependent pulsed fractions. We use the results of recent calculations for strongly magnetized atmospheres to obtain the angle- and energy-dependence of the surface emission. We include in our calculations the effects of general relativistic photon transport and interstellar extinction. We find that the combination of the large pulsed fractions and the high luminosities of AXPs cannot be accounted for by surface emission from a magnetar with two antipodal hot regions or a temperature distribution characteristic of a magnetic dipole. This result is robust for reasonable neutron star radii, for the range of magnetic field strengths inferred from the observed spin down rates, and for surface temperatures consistent with the spectral properties of AXPs. Models with a single hot emitting region can reproduce the observations, provided that the distance to one of the sources is ~30% less than the current best estimate, and allowing for systematic uncertainties in the spectral fit of a second source. Finally, the thermal emission models with antipodal emission geometry predict a characteristic strong increase of the pulsed fraction with photon energy, which is apparently inconsistent with the current data. The energy-dependence of the pulsed fraction in the models with one hot region shows a wider range of behavior and can be consistent with the existing data. Upcoming high-resolution observations with Chandra and XMM-Newton will provide a conclusive test.

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.

I study the formation of Comptonization spectra in spherically symmetric, fast moving media in a flat spacetime. I analyze the mathematical character of the moments of the transfer equation in the system-frame and describe a numerical method that provides fast solutions of the time-independent radiative transfer problem that are accurate in both the diffusion and free-streaming regimes. I show that even if the flows are mildly relativistic (V~0.1, where V is the electron bulk velocity in units of the speed of light), terms that are second-order in V alter the emerging spectrum both quantitatively and qualitatively. In particular, terms that are second-order in V produce power-law spectral tails, which are the dominant feature at high energies, and therefore cannot be neglected. I further show that photons from a static source are upscattered by the bulk motion of the medium even if the velocity field does not converge. Finally, I discuss these results in the context of radial accretion onto and outflows from compact objects.

Black holes and weakly-magnetic accreting neutron stars show strong quasi-periodic variability over timescales that span a very wide range. This property offer us a unique opportunity to improve our understanding of basic physical processes in strong gravitational fields as it reveals, for the first time, phenomena that occur on the most fundamental timescales near accreting compact objects. In this article, I review our current understanding of the variability properties of accretion flows onto neutron stars and black holes. I discuss a number of challenges faced by theoretical models, in order to identify open questions in accretion physics that need to be addressed. Finally, I discuss the relation to and implications for variability models of the detection of nearly coherent oscillations during type I X-ray bursts in several neutron-star sources.

We study the spectra generated as a result of bulk Comptonization by relativistic electrons in radial flows onto compact objects. We solve numerically the general-relativistic radiative transfer equation in the Schwarzschild spacetime, in steady-state and under minimal assumptions. We show that power-law spectra result from multiple scatterings, in a way similar to thermal Comptonization. We also find that photon-electron interactions taking place near the black-hole event horizon affect very little the emerging spectra. We, therefore, argue that bulk Comptonization spectra do not carry distinguishing signatures of the compact object around which they are produced. We examine the dependence of the spectra on simplifications often employed regarding the spacetime geometry, the distribution of photon sources, and the boundary conditions. We show that the existence of trapped characteristics around a black hole reduces the efficiency of Comptonization and that general relativistic effects identically cancel bulk Comptonization effects for a free-falling flow and in the limit of infinitesimal mean-free path. As a result, we find that neglecting the curved spacetime geometry leads to overestimating the high-energy flux by up to an order of magnitude. Finally, we demonstrate that the spectrum from accretion onto a neutron star depends sensitively on the imposed boundary conditions, while that from a black hole is immune to such choices.

The variability properties of accreting compact objects predicted by dynamical models are characterized by a number of distinct frequencies that are specific to each model. Because the accretion disks around neutron stars and black holes are hydrodynamic flows, these modulation frequencies cannot be strictly equal to test particle frequencies. I discuss several implications of hydrodynamic corrections to the modulation frequencies predicted by dynamical models. Finally, I show that the recent detection of yet a third kHz QPO in three neutron-star systems favors a previously developed model that attributes the various QPO frequencies to fundamental general relativistic frequencies in the accretion flow.

Anomalous X-ray pulsars, compact non-pulsing X-ray sources in supernova remnants, and X-ray bursters are three distinct types of sources for which there are viable models that attribute their X-ray emission to thermal emission from the surface of a neutron star. Inferring the surface area of the emitting regions in such systems is crucial in assessing the viability of different models and in providing bounds on the radii of neutron stars. We show that the inferred areas of the emitting regions may be over- or under-estimated by a factor of <=2, because of the geometry of the system and general relativistic light deflection, combined with the effects of phase averaging. Such effects make the determination of neutron-star radii uncertain, especially when compared to the ~5% level required for constraining the equation of state of neutron-star matter. We also note that, for a given spectral shape, the inferred source luminosities and pulse fractions are anticorrelated because they depend on the same properties of the emitting regions, namely their sizes and orientations, i.e., brighter sources have on average weaker pulsation amplitudes than fainter sources. We argue that this property can be used as a diagnostic tool in distinguishing between different spectral models. As an example, we show that the high inferred pulse fraction and brightness of the pulsar RXS J1708-40 are inconsistent with isotropic thermal emission from a neutron-star surface. Finally, we discuss the implication of our results for surveys in the soft X-rays for young, cooling neutron stars in supernova remnants and show that the absence of detectable pulsations from the compact source at the center of Cas A (at a level of >=30%) is not a strong argument againts its identification with a spinning neutron star.

Most models of anomalous X-ray pulsars (AXPs) account for the observed X-ray spectra and pulsations by means of radiation processes that occur on the surfaces of neutron stars. For any such model, general relativistic deflection of light severely suppresses the amplitude of the observed pulsations. We calculate the expected pulsation amplitudes of AXPs according to various models and compare the results with observations. We show that the high (<= 70%) pulse amplitudes observed in some AXPs can be accounted for only if the surface emission is localized (spot radius <40 degrees) and strongly beamed (cos^n[theta'] with n>2, where theta' is the angle to the normal). These constraints are incompatible with those cooling and magnetar models in which the observed X-rays originate as thermal emission from the neutron-star surface. Accretion models, on the other hand, are compatible with observations for a wide range of parameters. Finally, definitive conclusions cannot be reached on magnetospheric models, since their localization and beaming properties are not well understood.

We investigate the effect of a hybrid electron population, consisting of both thermal and non-thermal particles, on the synchrotron spectrum, image size, and image shape of a hot accretion flow onto a supermassive black hole. We find two universal features in the emitted synchrotron spectrum: (i) a prominent shoulder at low (< 10^11 Hz) frequencies that is weakly dependent on the shape of the electron energy distribution, and (ii) an extended tail of emission at high (> 10^13 Hz) frequencies whose spectral slope depends on the slope of the power-law energy distribution of the electrons. In the low-frequency shoulder, the luminosity can be up to two orders of magnitude greater than with a purely thermal plasma even if only a small fraction (< 1%) of the steady-state electron energy is in the non-thermal electrons. We apply the hybrid model to the Galactic center source, Sgr A*. The observed radio and IR spectra imply that at most 1% of the steady-state electron energy is present in a power-law tail in this source. This corresponds to no more than 10% of the electron energy injected into the non-thermal electrons and hence 90% into the thermal electrons. We show that such a hybrid distribution can be sustained in the flow because thermalization via Coulomb collisions and synchrotron self-absorption are both inefficient. The presence of non-thermal electrons enlarges the size of the radio image at low frequencies and alters the frequency dependence of the brightness temperature. A purely thermal electron distributions produces a sharp-edged image while a hybrid distribution causes strong limb brightening. These effects can be seen up to frequencies ~10^11 Hz and are accessible to radio interferometers.

Accreting neutron stars and black holes share a number of rapid variability characteristics, such as quasi-periodic oscillations and broad-band noise. The frequencies of these features were recently shown to be comparable to fundamental general relativistic frequencies of perturbed test-particle orbits near the compact objects. In this paper we propose a dynamical model for producing modulations in the inner disk properties at these characteristic frequencies. We postulate that a transition radius exists in the accretion disk close to the compact object and investigate its response to a broad spectrum of density fluctuations. We find that the transition radius acts as a low band-pass filter with a broadband response that is constant at frequencies lower than the inverse radial drift timescale and decreases above it as the inverse of the frequency. Moreover, the response shows strong resonances at frequencies nearly equal to the orbital, periastron-precession, and nodal-precession frequencies, thereby accounting for the observed power-density spectra of accreting compact objects and the correlations between their characteristic frequencies.

We calculate upper bounds on the inner radii of geometrically thin accretion disks in galactic black-hole systems by relating their rapid variability properties to those of neutron stars. We infer that the inner disk radii do not exhibit large excursions between different spectral states, in contrast with the concept that the disk retreats significantly during the soft-to-hard state transition. We find that, in the hard state, the accretion disks extend down to radii less than 6-25 GM/c^2 and discuss the implications of our results for models of black-hole X-ray spectra.

Recent X-ray variability observations of accreting neutron stars may provide the first evidence for frame dragging effects around spinning relativistic objects. Motivated by this possibility and its implications for neutron-star structural properties, we calculate new optimal bounds on the masses, radii, and moments of inertia of slowly rotating neutron stars that show kilohertz quasi-periodic oscillations (QPOs). These bounds are derived under minimal assumptions about the properties of matter at high densities and therefore are largely independent of the unknown equation of state. We further derive a semi-analytical upper bound on the neutron-star moment of inertia without making any assumptions about the equation of state of matter at any density. We use this upper bound to show that the maximum possible nodal precession frequency of an inclined circular orbit around a slowly spinning neutron star is nu_{NP}~45.2 (nu_s/300Hz) Hz, where nu_s is the spin frequency of the neutron star. We conclude that the nodal-precession interpretation of low-frequency QPOs in accreting neutron stars is inconsistent with their inferred spin frequencies or the identification of the highest-frequency QPO with that of a circular Keplerian orbit in the accretion disk.

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.

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 study systematically the ~0.1-1200 Hz quasi-periodic oscillations (QPOs) and broad noise components observed in the power spectra of non-pulsing neutron-star and black-hole low-mass X-ray binaries. We show that among these components we can identify two, occurring over a wide range of source types and luminosities, whose frequencies follow a tight correlation. The variability components involved in this correlation include neutron-star kilohertz QPOs and horizontal-branch oscillations, as well as black-hole QPOs and noise components. Our results suggest that the same types of variability may occur in both neutron-star and black-hole systems over three orders of magnitude in frequency and with coherences that vary widely but systematically. Confirmation of this hypothesis will strongly constrain theoretical models of these phenomena and provide additional clues to understanding their nature.

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.

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.

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.

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.

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.

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.