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.

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.

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.

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.

The nuclear symmetry energy plays a role in determining both the nuclear properties of terrestrial matter as well as the astrophysical properties of neutron stars. The first measurement of the neutron star tidal deformability, from gravitational wave event GW170817, provides a new way of probing the symmetry energy. In this work, we report on new constraints on the symmetry energy from GW170817. We focus in particular on the low-order coefficients: namely, the value of the symmetry energy at the nuclear saturation density, S_0, and the slope of the symmetry energy, L_0. We find that the gravitational wave data are relatively insensitive to S_0, but that they depend strongly on L_0 and point to lower values of L_0 than have previously been reported, with a peak likelihood near L_0 ~ 20 MeV. Finally, we use the inferred posteriors on L_0 to derive new analytic constraints on higher-order nuclear terms.

In this white paper we explore the capabilities required to identify and study supermassive black holes formed from heavy seeds ($\mathrm{M_{\bullet}} \sim 10^4 - 10^6 \, \mathrm{M_{\odot}}$) in the early Universe. To obtain an unequivocal detection of heavy seeds we need to probe mass scales of $\sim 10^{5-6} \, \mathrm{M_{\odot}}$ at redshift $z \gtrsim 10$. From this theoretical perspective, we review the observational requirements and how they compare with planned/proposed instruments, in the infrared, X-ray and gravitational waves realms. In conclusion, detecting heavy black hole seeds at $z \gtrsim 10$ in the next decade will be challenging but, according to current theoretical models, feasible with upcoming/proposed facilities. Their detection will be fundamental to understand the early history of the Universe, as well as its evolution until now. Shedding light on the dawn of black holes will certainly be one of the key tasks that the astronomical community will focus on in the next decade.

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.

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

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

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.

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

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.

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

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.

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.

The merger of a neutron star binary may result in the formation of a rapidly-spinning magnetar. The magnetar can potentially survive for seconds or longer as a supramassive neutron star before collapsing to a black hole if, indeed, it collapses at all. During this process, a fraction of the magnetar's rotational energy of ~10^53 erg is transferred via magnetic spin-down to the surrounding ejecta. The resulting interaction between the ejecta and the surrounding circumburst medium powers a >year-long synchrotron radio transient. We present a search for radio emission with the Very Large Array following nine short-duration gamma-ray bursts (GRBs) at rest-frame times of ~1.3-7.6 years after the bursts, focusing on those events which exhibit early-time excess X-ray emission that may signify the presence of magnetars. We place upper limits of <18-32 microJy on the 6.0 GHz radio emission, corresponding to spectral luminosities of <(0.05-8.3)x10^39 erg s^-1. Comparing these limits to the predicted radio emission from a long-lived remnant and incorporating measurements of the circumburst densities from broad-band modeling of short GRB afterglows, we rule out a stable magnetar with an energy of 10^53 erg for half of the events in our sample. A supramassive remnant that injects a lower rotational energy of 10^52 erg is ruled out for a single event, GRB 050724A. This study represents the deepest and most extensive search for long-term radio emission following short GRBs to date, and thus the most stringent limits placed on the physical properties of magnetars associated with short GRBs from radio observations.

We present a comprehensive study of interstellar X-ray extinction using the extensive Chandra supernova remnant archive and use our results to refine the empirical relation between the hydrogen column density and optical extinction. In our analysis, we make use of the large, uniform data sample to assess various systematic uncertainties in the measurement of the interstellar X-ray absorption. Specifically, we address systematic uncertainties that originate from (i) the emission models used to fit supernova remnant spectra, (ii) the spatial variations within individual remnants, (iii) the physical conditions of the remnant such as composition, temperature, and non-equilibrium regions, and (iv) the model used for the absorption of X-rays in the interstellar medium. Using a Bayesian framework to quantify these systematic uncertainties, and combining the resulting hydrogen column density measurements with the measurements of optical extinction toward the same remnants, we find the empirical relation NH = (2.87+/-0.12) x 10^21 AV cm^(-2), which is significantly higher than the previous measurements.

The mass function of neutron stars (NSs) contains information about the late evolution of massive stars, the supernova explosion mechanism, and the equation-of-state of cold, nuclear matter beyond the nuclear saturation density. A number of recent NS mass measurements in binary millisecond pulsar (MSP) systems increase the fraction of massive NSs (with $M > 1.8$ M$_{\odot}$) to $\sim 20\% $ of the observed population. In light of these results, we employ a Bayesian framework to revisit the MSP mass distribution. We find that a single Gaussian model does not sufficiently describe the observed population. We test alternative empirical models and infer that the MSP mass distribution is strongly asymmetric. The diversity in spin and orbital properties of high-mass NSs suggests that this is most likely not a result of the recycling process, but rather reflects differences in the NS birth masses. The asymmetry is best accounted for by a bimodal distribution with a low mass component centred at $1.393_{-0.029}^{+0.031}$ M$_{\odot}$ and dispersed by $0.064_{-0.025}^{+0.064}$ M$_{\odot}$, and a high-mass component with a mean of $1.807_{-0.132}^{+0.081}$ and a dispersion of $0.177_{-0.072}^{+0.115}$ M$_{\odot}$. We also establish a lower limit of $M_{max} \ge 2.018$ M$_{\odot}$ at 98% C.L. for the maximum NS mass, from the absence of a high-mass truncation in the observed masses. Using our inferred model, we find that the measurement of 350 MSP masses, expected after the conclusion of pulsar surveys with the Square-Kilometre Array, can result in a precise localization of a maximum mass up to 2.15 M$_{\odot}$, with a 5% accuracy. Finally, we identify possible massive NSs within the known pulsar population and discuss birth masses of MSPs.

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.

We summarize our current knowledge of neutron star masses and radii. Recent instrumentation and computational advances have resulted in a rapid increase in the discovery rate and precise timing of radio pulsars in binaries in the last few years, leading to a large number of mass measurements. These discoveries show that the neutron star mass distribution is much wider than previously thought, with 3 known pulsars now firmly in the 1.9-2.0 Msun mass range. For radii, large, high quality datasets from X-ray satellites as well as significant progress in theoretical modeling led to considerable progress in the measurements, placing them in the 9.9-11.2 km range and shrinking their uncertainties due to a better understanding of the sources of systematic errors. The combination of the massive neutron star discoveries, the tighter radius measurements, and improved laboratory constraints of the properties of dense matter has already made a substantial impact on our understanding of the composition and bulk properties of cold nuclear matter at densities higher than that of the atomic nucleus, a major unsolved problem in modern physics.

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.

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

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.

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.

We model the pulse profiles and the phase resolved spectra of the anomalous X-ray pulsar 1E 1048.1-5937 obtained with XMM-Newton to map its surface temperature distribution during an active and a quiescent epoch. We develop and apply a model that takes into account the relevant physical and geometrical effects on the neutron star surface, magnetosphere, and spacetime. Using this model, we determine the observables at infinity as a function of pulse phase for different numbers and sizes of hot spots on the surface. We show that the pulse profiles extracted from both observations can be modeled with a single hot spot and an antipodal cool component. The size of the hot spot changes from $\approx 80^{\circ}$ in 2007, 3 months after the onset of a dramatic flux increase, to $\approx 30^{\circ}$ during the quiescent observation in 2011, when the pulsed fraction returned to the pre-outburst $\approx$ 65\% level. For the 2007 observation, we also find that a model consisting of a single 0.4 keV hot spot with a magnetic field strength of $1.8 \times 10^{14}$ G accounts for the spectra obtained at three different pulse phases but under predicts the flux at the pulse minimum, where the contribution to the emission from the cooler component is non-negligible. The inferred temperature of the spot stays approximately constant between different pulse phases, in agreement with a uniform temperature, single hot spot model. These results suggest that the emitting area grows significantly during outbursts but returns to its persistent and significantly smaller size within a few year timescale.

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.

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.

The past three decades have seen prodigious advances in astronomy and astrophysics. Beginning with the exploration of our solar system and continuing through the pioneering Explorers and Great Observatories of today, NASA missions have made essential contributions to these advances. This roadmap presents a science-driven 30-year vision for the future of NASA Astrophysics that builds on these achievements to address some of our most ancient and fundamental questions: Are we alone? How did we get here? How does the universe work? The search for the answers constitutes the Enduring Quests of this roadmap. Building on the priorities identified in New Worlds, New Horizons, we envision future science investigations laid out in three Eras, with each representing roughly ten years of mission development in a given field. The immediate Near-Term Era covers ongoing NASA-led activities and planned missions. This will be followed by the missions of the Formative Era, which will build on the preceding technological developments and scientific discoveries, with remarkable capabilities that will enable breakthroughs across the landscape of astrophysics. These will then lay the foundations for the Daring Visions of the Visionary Era: missions and explorations that will take us deep into unchartered scientific and technological terrain. The roadmap outlined herein will require the vision and wherewithal to undertake highly ambitious programs over the next 30 years. The discoveries that emerge will inspire generations of citizen scientists young and old, and inspire all of humanity for decades to come.

(Abridged:) We present the identification of optical counterparts to 23 Galactic Bulge Survey X-ray sources. We report their accurate coordinates and optical spectra acquired at the VLT and Magellan. All sources are classified as accreting binaries according to their emission line characteristics. To distinguish accreting binaries from chromospherically active objects we develop and explain criteria based on Halpha and HeI 5786,6678 emission line properties available in the literature. The spectroscopic properties and photometric variability of all the objects are discussed and a classification of the source is given where possible. Among the 23 systems, at least 9 of them show an accretion-dominated optical spectrum (CX28, CX63, CX70, CX128, CX142, CX207, CX522, CX794, CX1011) and another 6 show photospheric lines from a late-type donor star in addition to accretion disc emission (CX44, CX93, CX137, CX154, CX377 and CX1004) indicating that they are probably accreting binaries in quiescence or in a low accretion rate state. Two sources are confirmed to be eclipsing: CX207 and CX794. CX207 shows a broad asymmetric Halpha profile blue-shifted by >300 km/s. Such line profile characteristics are consistent with a magnetic (Polar) CV. CX794 is an eclipsing nova-like CV in the period gap. Time-resolved photometry and the large broadening of the Halpha emission lines in CX446 (2100 km/s FWHM) suggest that this is also an eclipsing or high-inclination accreting binary. Finally, the low-accretion rate source CX1004 shows a double-peaked Halpha profile with a FWHM of 2100 km/s. This supports a high inclination or even eclipsing system. Whether the compact object is a white dwarf in an eclipsing CV or a black hole primary in a high-inclination LMXB remains to be established.

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.

We perform detailed magnetohydrodynamic simulations of the gas cloud G2 interacting with the accretion flow around the Galactic Center black hole Sgr A*. We take as our initial conditions a steady-state, converged solution of the accretion flow obtained earlier using the general-relativistic magnetohydrodynamic code HARM. Using the observed parameters for the cloud's orbit, we compute the interaction of the cloud with the ambient gas and identify the shock structure that forms ahead of the cloud. We show that for many configurations, the cloud front crosses orbit pericenter 7 to 9 months earlier than the center-of-mass.

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.

We calculate radio light curves produced by the bow shock that is likely to form in front of the G2 cloud when it penetrates the accretion disk of Sgr A*. The shock acceleration of the radio-emitting electrons is captured self-consistently by means of first-principles particle-in-cell simulations. We show that the radio luminosity is expected to reach maximum in early 2013, roughly a month after the bow shock crosses the orbit pericenter. We estimate the peak radio flux at 1.4 GHz to be 1.4 - 22 Jy depending on the assumed orbit orientation and parameters. We show that the most promising frequencies for radio observations are in the 0.1<nu<1 GHz range, for which the bow shock emission will be much stronger than the intrinsic radio flux for all the models considered.

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.

We use time resolved spectroscopy of thermonuclear X-ray bursts observed from SAX J1748.9-2021 to infer the mass and the radius of the neutron star in the binary. Four X-ray bursts observed from the source with RXTE enable us to measure the angular size and the Eddington limit on the neutron star surface. Combined with a distance measurement to the globular cluster NGC 6440, in which SAX J1748.9-2021 resides, we obtain two solutions for the neutron star radius and mass, R = 8.18 +/- 1.62 km and M = 1.78 +/- 0.3 M_\sun or R = 10.93 +/- 2.09 km and M = 1.33 +/- 0.33 M_\sun.

Neutron stars are associated with diverse physical phenomena that take place in conditions characterized by ultrahigh densities as well as intense gravitational, magnetic, and radiation fields. Understanding the properties and interactions of matter in these regimes remains one of the challenges in compact object astrophysics. Photons emitted from the surfaces of neutron stars provide direct probes of their structure, composition, and magnetic fields. In this review, I discuss in detail the physics that governs the properties of emission from the surfaces of neutron stars and their various observational manifestations. I present the constraints on neutron star radii, core and crust composition, and magnetic field strength and topology obtained from studies of their broadband spectra, evolution of thermal luminosity, and the profiles of pulsations that originate on their surfaces.

A dense ionized cloud of gas has been recently discovered to be moving directly toward the supermassive black hole, Sgr A*, at the Galactic Center. In June 2013, at the pericenter of its highly eccentric orbit, the cloud will be approximately 3100 Schwarzschild radii from the black hole and will move supersonically through the ambient hot gas with a velocity of v_p ~ 5400 km/s. A bow shock is likely to form in front of the cloud and could accelerate electrons to relativistic energies. We estimate via particle-in-cell simulations the energy distribution of the accelerated electrons and show that the non-thermal synchrotron emission from these electrons might exceed the quiescent radio emission from Sgr A* by a factor of several. The enhanced radio emission should be detectable at GHz and higher frequencies around the time of pericentric passage and in the following months. The bow shock emission is expected to be displaced from the quiescent radio emission of Sgr A* by ~33 mas. Interferometric observations could resolve potential changes in the radio image of Sgr A* at wavelengths < 6 cm.

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 present a systematic analysis of all archival Chandra observations of the soft-gamma repeater SGR 0526-66. Our results show that the X-ray flux of SGR 0526-66 decayed by about 20% between 2000 and 2009. We employ physically motivated X-ray spectral models and determine the effective temperature and the strength of the magnetic field at the surface as kT = 0.354_{-0.024}^{+0.031} keV and B = (3.73^{+0.16}_{-0.08})x10^{14} G, respectively. We find that the effective temperature remains constant within the statistical uncertainties and attribute the decrease in the source flux to a decrease in the emitting radius. We also perform timing analysis to measure the evolution of the spin period and the period derivative over the nine year interval. We find a period derivative of \.P = (4.0 +/- 0.5)x10^{-11} ss^{-1}, which allows us to infer the dipole magnetic field strength and compare it with the one determined spectroscopically. Finally, we compare the effective temperature of SGR 0526-66 with the expected cooling trends from magnetized neutron stars and suggest an initial magnetic field strength of 10^{15-16} G for the source.

We present optical, X-ray and radio observations of the black hole transient (BHT) XTE J1752-223 towards and in quiescence. Optical photometry shows that the quiescent magnitude of XTE J1752-223 is fainter than 24.4 magnitudes in the i'-band. A comparison with measurements of the source during its 2009-2010 outburst shows that the outburst amplitude is more than 8 magnitudes in the i'-band. Known X-ray properties of the source combined with the faintness of the quiescence optical counterpart and the large outburst optical amplitude point towards a short orbital period system (Porb<~6.8 h) with an M type (or later) mass donor, at a distance of 3.5<~d<~8 kpc. Simultaneous X-ray and radio data were collected with Chandra and the EVLA, allowing constraints to be placed on the quiescent X-ray and radio flux of XTE J1752-223. Furthermore, using data covering the final stage of the outburst decay, we investigated the low luminosity end of the X-ray - radio correlation for this source and compared it with other BHTs. We found that XTE J1752-223 adds to the number of outliers with respect to the `standard' X-ray - radio luminosity relation. Furthermore, XTE J1752-223 is the second source, after the BHT H1743-322, that shows a transition from the region of the outliers towards the `standard' correlation at low luminosity. Finally, we report on a faint, variable X-ray source we discovered with Chandra at an angular distance of ~2.9" to XTE J1752-223 and at a position angle consistent with that of the radio jets previously observed from the BHT. We discuss the possibility that we detected X-ray emission associated with a jet from XTE J1752-223.

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.

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 observed upper bound on the spin down rate of the otherwise typical Soft Gamma-ray Repeater SGR 0418+5729 has challenged the interpretation of this source as a neutron star with ultrastrong magnetic fields. Current limits imply a dipole magnetic field strength of less than 7.5 x 10^{12} G (Rea et al. 2010), which is significantly smaller than that of a typical SGR. Independent of the properties inferred from X-ray timing, the X-ray spectra of neutron stars allow a measurement of their magnetic field strengths because they are distorted from pure blackbodies due to the presence of a magnetic field in a radiative equilibrium atmosphere. In this paper, we model high signal-to-noise XMM-Newton spectra of SGR 0418+5729 to place constraints on the strength of the magnetic field at the surface of the neutron star. Our analysis shows that neutron star atmosphere models with moderate magnetic field strengths (10^{12-13} G) cannot fit the X-ray spectra, whereas, models with stronger magnetic fields are able to account for the observations. We find that the strength of the magnetic field at the surface is 1.0 x 10^{14} G. This value, although lower than all of the other SGRs analyzed to date, is still high enough to generate the observed X-ray bursts from the source. In connection to the spindown limits, it also implies a significantly non-dipolar structure of the magnetic field. We discuss the results of our spectral modeling and compare them with other SGRs.

Measurements of neutron star masses and radii are instrumental for determining the equation of state of their interiors, understanding the dividing line between neutron stars and black holes, and for obtaining accurate statistics of source populations in the Galaxy. We report here on the measurement of the mass and radius of the neutron star in the low-mass X-ray binary KS 1731-260. The analysis of the spectroscopic data on multiple thermonuclear bursts yields well-constrained values for the apparent angular area and the Eddington flux of the source, both of which depend in a distinct way on the mass and radius of the neutron star. The binary KS 1731-260 is in the direction of the Galactic bulge, allowing a distance estimate based on the density of stars in that direction. Making use of the Han & Gould model, we determine the probability distribution over the distance to the source, which is peaked at 8 kpc. Combining these measurements, we place a strong upper bound on the radius of the neutron star, R <= 12 km, while confining its mass to M <= 1.8 M_sun.

We introduce the Galactic Bulge Survey (GBS) and we provide the Chandra source list for the region that has been observed to date. Among the goals of the GBS are constraining the neutron star equation of state and the black hole mass distribution via the identification of eclipsing neutron star and black hole low-mass X-ray binaries. The latter goal will, in addition, be obtained by significantly enlarging the number of black hole systems for which a black hole mass can be derived. Further goals include constraining X-ray binary formation scenarios, in particular the common envelope phase and the occurrence of kicks, via source-type number counts and an investigation of the spatial distribution of X-ray binaries, respectively. The GBS targets two strips of 6x1 degrees (12 square degrees in total), one above (1<b<2 degrees) and one below (-2<b<-1 degrees) the Galactic plane in the direction of the Galactic Center at both X-ray and optical wavelengths. By avoiding the Galactic plane (-1<b<1 degrees) we limit the influence of extinction on the X-ray and optical emission but still sample relatively large number densities of sources. The survey is designed such that a large fraction of the X-ray sources can be identified from their optical spectra. The X-ray survey, by design, covers a large area on the sky while the depth is shallow using 2 ks per Chandra pointing. In this way we maximize the predicted number ratio of (quiescent) low-mass X-ray binaries to Cataclysmic Variables. The survey is approximately homogeneous in depth to an 0.5-10 keV flux of 7.7x10^-14 erg cm-2 s-1. So far, we have covered about two-thirds (8.3 square degrees) of the projected survey area with Chandra providing over 1200 unique X-ray sources. We discuss the characteristics and the variability of the brightest of these sources.

The prolific magnetar SGR 1900+14 showed two outbursts in the last decade and has been closely monitored in the X-rays to track the changes in its radiative properties. We use archival Chandra and XMM-Newton observations of SGR 1900+14 to construct a history of its spectrum and persistent X-ray flux spanning a period of about seven years. We show that the decline of its X-ray flux in these two outburst episodes follows the same trend. The flux begins to decline promptly and rapidly subsequent to the flares, then decreases gradually for about 600 days, at which point it resumes a more rapid decline. Utilizing the high quality spectral data in each epoch, we also study the spectral coevolution of the source with its flux. We find that neither the magnetic field strength nor the magnetospheric properties change over the period spanned by the observations, while the surface temperature as well as the inferred emitting area both decline with time following both outbursts. We also show that the source reached the same minimum flux level in its decline from these two subsequent outbursts, suggesting that this flux level may be its steady quiescent flux.

We present results from observations of the magnetar 1E 1547.0-5408 (SGR J1550-5418) taken with the Chandra X-ray Observatory and the Rossi X-ray Timing Explorer (RXTE) following the source's outbursts in 2008 October and 2009 January. During the time span of the Chandra observations, which covers days 4 through 23 and days 2 through 16 after the 2008 and 2009 events, respectively, the source spectral shape remained stable, while the pulsar's spin-down rate in the same span in 2008 increased by a factor of 2.2 as measured by RXTE. The lack of spectral variation suggests decoupling between magnetar spin-down and radiative changes, hence between the spin-down-inferred magnetic field strength and that inferred spectrally. We also found a strong anti-correlation between the phase-averaged flux and the pulsed fraction in the 2008 and 2009 Chandra data, but not in the pre-2008 measurements. We discuss these results in the context of the magnetar model.

We present the first astrophysical measurement of the pressure of cold matter above nuclear saturation density, based on recently determined masses and radii of three neutron stars. The pressure at higher densities are below the predictions of equations of state that account only for nucleonic degrees of freedom, and thus present a challenge to the microscopic theory of neutron star matter.

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.

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.

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 report on the measurement of the mass and radius of the neutron star in the low-mass X-ray binary 4U 1820-30. The analysis of the spectroscopic data on multiple thermonuclear bursts yields well-constrained values for the apparent emitting area and the Eddington flux, both of which depend in a distinct way on the mass and radius of the neutron star. The distance to the source is that of the globular cluster NGC 6624, where the source resides. Combining these measurements, we uniquely determine the probability density over the stellar mass and radius. We find the mass to be M = 1.58 +/- 0.06 M_sun and the radius to be R = 9.11 +/- 0.40 km.

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.

Low mass X-ray binaries that show thermonuclear bursts are ideal sources for constraining the equation of state of neutron star matter. The lack of independent distance measurements for most of these sources, however, prevents a systematic exploration of the masses and radii of the neutron stars, hence limiting the equation of state studies. We present here a measurement of the distance to the low mass X-ray binary 4U 1608-52 that is based on the study of the interstellar extinction towards the source. We first model the individual absorption edges of the elements Ne and Mg in the high resolution X-ray spectrum obtained with XMM-Newton. We then combine this information with a measurement of the run of reddening with distance using red clump stars and determine a minimum distance to the source of 3.9 kpc, with a most probable value of 5.8 kpc. Finally, we analyze time-resolved X-ray spectra of Type-I X-ray bursts observed from this source to measure the mass and the radius of the neutron star. We find a mass of M = 1.74 +/- 0.14 M_sun and a radius of R = 9.3 +/- 1.0 km, respectively. This mass and radius can be achieved by several multi-nucleon equations of state.

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

A linear relation between the hydrogen column density (N_H) and optical extinction (A_V) in the Galaxy has long been observed. A number of studies found differing results in the slope of this relation. Here, we utilize the data on 22 supernova remnants that have been observed with the latest generation X-ray observatories and for which optical extinction and/or reddening measurements have been performed and find N_H (cm^-2) = (2.21 \pm0.09) x10^{21} A_V (mag). We compare our result with the previous studies and assess any systematic uncertainties that may affect these results.

We study the topology of reionization using accurate three-dimensional radiative transfer calculations post-processed on outputs from cosmological hydrodynamic simulations. In our simulations, reionization begins in overdense regions and then "leaks" directly into voids, with filaments reionizing last owing to their combination of high recombination rate and low emissivity. This result depends on the uniquely-biased emissivity field predicted by our prescriptions for star formation and feedback, which have previously been shown to account for a wide array of measurements of the post-reionization Universe. It is qualitatively robust to our choice of simulation volume, ionizing escape fraction, and spatial resolution (in fact it grows stronger at higher spatial resolution) even though the exact overlap redshift is sensitive to each of these. However, it weakens slightly as the escape fraction is increased owing to the reduced density contrast at higher redshift. We also explore whether our results are sensitive to commonly-employed approximations such as using optically-thin Eddington tensors or substantially altering the speed of light. Such approximations do not qualitatively change the topology of reionization. However, they can systematically shift the overlap redshift by up to $\Delta z\sim 0.5$, indicating that accurate radiative transfer is essential for computing reionization. Our model cannot simultaneously reproduce the observed optical depth to Thomson scattering and ionization rate per hydrogen atom at $z=6$, which could owe to numerical effects and/or missing early sources of ionization.

The cores of neutron stars harbor the highest matter densities known to occur in nature, up to several times the densities in atomic nuclei. Similarly, magnetic field strengths can exceed the strongest fields generated in terrestrial laboratories by ten orders of magnitude. Hyperon-dominated matter, deconfined quark matter, superfluidity, even superconductivity are predicted in neutron stars. Similarly, quantum electrodynamics predicts that in strong magnetic fields the vacuum becomes birefringent. The properties of matter under such conditions is governed by Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED), and the close study of the properties of neutron stars offers the unique opportunity to test and explore the richness of QCD and QED in a regime that is utterly beyond the reach of terrestrial experiments. Experimentally, this is almost virgin territory.

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.

We measured the galactic hydrogen column densities to the neutron-star binaries GX 17+2, 4U 1705-44, and 4U 1728-34 by modeling the Mg and Si absorption edges found in high-resolution X-ray spectra obtained by the Chandra X-ray Observatory. We found for GX 17+2, N_H = (2.38 +/- 0.12) x 10^22 cm^-2, for 4U 1705-44, N_H = (2.44 +/- 0.09) x 10^22 cm^-2, and for 4U 1728-34, N_H = (2.49 +/- 0.14) x 10^22 cm^-2. These values are in reasonable agreement with the hydrogen column densities inferred earlier from modeling of the continuum spectra of the sources. Our results can be used to constrain the uncertainties of model parameters of the X-ray spectra of these sources that are correlated to the uncertainties of the hydrogen column density. In the case of continuum spectra obtained during thermonuclear X-ray bursts, they will significantly reduce the uncertainties in the spectroscopically measured masses and radii of the neutron stars.

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.

We investigate the dependence of pulse amplitudes of accreting millisecond pulsars on the masses of the neutron stars. Because the pulsation amplitudes are suppressed as the neutron stars become more massive, the probability of detection of pulsations decreases in systems that have been accreting for a long time. However, the probability of detectable pulsations is higher in transient systems where the mass accretion is sporadic and the neutron star is likely to have a low mass. We propose this mechanism as the explanation of the small number of millisecond X-ray pulsars found to date, as well as their emergence as fast pulsars mostly in transient, low-Mdot systems. This mechanism can also quantitatively explain the lack of pulsars in the majority of LMXBs.

We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation, closed by an Eddingtion tensor computed using a long characteristics method. We show that pure (i.e., not source-centered) short characteristics and the optically-thin approximation are inappropriate for computing Eddington factors for the problem of cosmological reionization. We evolve the non-equilibrium ionization field via an efficient and accurate (errors <1%) technique that switches between fully implicit or explicit finite-differencing depending on whether the local timescales are long or short compared to the timestep. We tailor our code for the problem of cosmological reionization. In tests, the code conserves photons, accurately treats cosmological effects, and reproduces analytic Stromgren sphere solutions. Its chief weakness is that the computation time for the long characteristics calculation scales relatively poorly compared to other techniques (t_{LC} \propto N_{cells}^1.5); however, we mitigate this by only recomputing the Eddington tensor when the radiation field changes substantially. Our technique makes almost no physical approximations, so it provides a way to benchmark faster but more approximate techniques. It can readily be extended to evolve multiple frequencies, though we do not do so here. Finally, we note that our method is generally applicable to any problem involving the transfer of continuum radiation through a periodic volume.

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.

Anomalous X-ray Pulsars (AXPs) belong to a class of neutron stars believed to harbor the strongest magnetic fields in the universe, as indicated by their energetic bursts and their rapid spindowns. We have developed a theoretical model that takes into account processes in the atmospheres and magnetospheres of ultramagnetic neutron stars, as well as the effects of their strong gravitational fields on the observable properties. Using this model, we have analyzed the X-ray spectra of a number of AXPs. We find that in all cases, the X-ray spectra are described very well with this emission model. The spectroscopically measured magnetic field strengths of these sources are in close agreement with the values inferred from their spindown properties and provide independent evidence for their magnetar nature. The analysis of spectral data using this physical model also sheds light on the long-term evolution of AXPs.

We analyze archival Chandra and XMM-Newton data of 4U 0142+61 within the context of the Surface Thermal Emission and Magnetospheric Scattering model. We show that 4U 0142+61 spectrum can be fit very well with this physical model that contains only four parameters. The system parameters can be tightly constrained from the fits, yielding a surface magnetic field strength of B=(4.75 \pm 0.02)x 10^{14} G, a surface temperature of kT = 0.309 \pm 0.001$ keV and a scattering optical depth of a few in the magnetosphere. These values do not vary between observations due to the stability of the source within the window of the observations. The detailed fits yield $\chi^2$ values that are statistically much better than the traditionally employed blackbody+power-law and two blackbody fits. The spectroscopically measured surface magnetic field strength is higher than, but within, the theoretical uncertainties of the value inferred from the dipole spindown formula.

Anomalous X-ray Pulsars (AXPs) belong to a class of neutron stars believed to harbor the strongest magnetic fields in the universe, as indicated by their energetic bursts and their rapid spindowns. However, an unambiguous measurement of their surface field strengths has not been made to date. It is also not known whether AXP outbursts result from changes in the neutron star magnetic field or crust properties. Here we report a spectroscopic measurement of the surface magnetic field strength of an AXP, XTE J1810$-$197, and solidify its magnetar nature. The field strength obtained from detailed spectral analysis and modeling, B = (2.72$\pm$0.03)$\times10^{14}$ G, is remarkably close to the value inferred from the rate of spindown of this source and remains nearly constant during numerous observations spanning over an order of magnitude in source flux. The surface temperature, on the other hand, declines steadily and dramatically following the 2003 outburst of this source. Our findings demonstrate that heating occurs in the upper neutron star crust during an outburst and sheds light on the transient behaviour of AXPs.

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.

We explore the hardness-intensity correlations observed in several AXPs and SGRs within the framework of a thermally emitting magnetar model. Using our detailed atmosphere models and taking into account reprocessing of the surface emission by the magnetosphere, we show that the hardness of the surface spectra increases with increasing temperature and hence the changes in the effective temperatures of the outer layers of the star alone can account for the observed correlations. We conclude that the slow release of the heat deposited in the deep crust during a magnetar burst naturally accounts for the spectral changes during the afterglow. The correlations are further enhanced by changes in the structures of the magnetic currents during or following a burst. However, the additional hardening produced by scattering of the surface photons off the magnetospheric charges saturates at moderate values of the scattering optical depth.

We present the first self-consistent theoretical models of magnetar spectra that take into account the combined effects of the stellar atmosphere and its magnetosphere. We find that the proton cyclotron lines that are already weakened by atmospheric effects become indistinguishable from the continuum for moderate scattering optical depths in the magnetosphere. Furthermore, the hard excess becomes more pronounced due to resonant scattering and the resulting spectra closely resemble the observed magnetar spectra. We argue that while the absence of proton cyclotron lines in the observed spectra are inconclusive about the surface field strengths of magnetars, the continuum carries nearly unique signatures of the field strength and can thus be used to infer this quantity. We fit our theoretical spectra with a phenomenological two-blackbody model and compare our findings to source spectra with existing two-blackbody fits. The field strengths that we infer spectroscopically are in remarkable agreement with the values inferred from the period derivative of the sources assuming a dipole spindown. These detailed predictions of line energies and equivalent widths may provide an optimal opportunity for measuring directly the surface field strengths of magnetars with future X-ray telescopes such as {\em Constellation}-X.

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.

The interiors of neutron stars contain matter at very high densities, in a state that differs greatly from those found in the early universe or achieved at terrestrial experiments. Matter in these conditions can only be probed through astrophysical observations that measure the mass and radius of neutron stars with sufficient precision. Here I report for the first time a unique determination of the mass and radius of the neutron star EXO 0748-676, which appears to rule out all the soft equations of state of neutron star matter. If this object is typical, then condensates and unconfined quarks do not exist in the centers of neutron stars.

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.

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.

We present precise localizations of AXP 1E1841-045 and SGR 1627-41 with Chandra. We obtained new infrared observations of SGR 1627-41 and reanalyzed archival observations of AXP 1E1841-045 in order to refine their positions and search for infrared counterparts. A faint source is detected inside the error circle of AXP 1E1841-045. In the case of SGR 1627-41, several sources are located within the error radius of the X-ray position and we discuss the likelihood of one of them being the counterpart. We compare the properties of our candidates to those of other known AXP and SGR counterparts. We find that the counterpart candidates for SGR 1627-41 and SGR 1806-20 would have to be intrinsically much brighter than AXPs to have detectable counterparts with the observational limits currently available for these sources. To confirm the reported counterpart of SGR 1806-20, we obtained new IR observations during the July 2003 burst activation of the source. No brightening of the suggested counterpart is detected, implying that the counterpart of SGR 1806-20 remains yet to be identified.

A number of Anomalous X-ray Pulsars (AXPs) have recently been detected in the optical/IR wavelengths. We use their inferred brightness to place general constraints on any model for this emission within the magnetar framework. We find that neutron-star surface emission cannot account for the observations and that the emission must be magnetospheric in origin. We propose a model for the optical/IR emission in which a distribution of energetic electrons in the neutron-star magnetosphere emits synchrotron radiation. This model can naturally reproduce the observed brightness and the rising spectra of AXPs as well as the observed pulsations at the stellar spin frequency and the correlation of the IR flux with their bursting activity.

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

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

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 consider the effects of vacuum polarization and proton cyclotron resonances on the propagation of radiation through a strongly magnetized plasma. We analyze the conditions under which the photons evolve adiabatically through the resonant density and find that the adibaticity condition is satisfied for most photon energies of interest, allowing for a normal-mode treatment of the photon propagation. We then construct radiative equilibrium atmosphere models of strongly magnetized neutron stars that includes these effects, employing a new numerical method that resolves accurately the sharp changes of the absorption and mode-coupling cross sections at the resonant densities. We show that the resulting spectra are modified by both resonances and are harder at all field strengths than a blackbody at the effective temperature. We also show that the narrow absorption features introduced by the proton cyclotron resonance have small equivalent widths. We discuss the implications of our results for properties of thermal emission from the surfaces of young neutron stars.

We study the pulse morphologies and pulse amplitudes of thermally emitting neutron stars with ultrastrong magnetic fields. The beaming of the radiation emerging from a magnetar was recently shown to be predominantly non-radial, with a small pencil and a broad fan component. We show that the combination of this radiation pattern with the effects of strong lensing in the gravitational field of the neutron star yields pulse profiles that show a qualitatively different behavior compared to that of the radially-peaked beaming patterns explored previously. Specifically, we find that: (i) the pulse profiles of magnetars with a single hot emission region on their surface exhibit 1-2 peaks, whereas those with an antipodal emission geometry have 1-4 peaks, depending on the neutron star compactness, the observer's viewing angle, and the size of the hot regions; (ii) the energy dependence of the beaming pattern may give rise to weakly or strongly energy-dependent pulse profiles and may introduce phase lags between different energy bands; (iii) the non-radial beaming pattern can give rise to high pulsed fractions even for very relativistic neutron stars; (iv) the pulsed fraction may not vary monotonically with neutron star compactness; (v) the pulsed fraction does not decrease monotonically with the size of the emitting region; (vi) the pulsed fraction from a neutron star with a single hot pole has, in general, a very weak energy dependence, in contrast to the case of an antipodal geometry. Comparison of these results to the observed properties of anomalous X-ray pulsars strongly suggests that they are neutron stars with a single hot region of ultrastrong magnetic field.

We construct radiative equilibrium models for strongly magnetized (B > 10^13 G) neutron-star atmospheres taking into account magnetic free-free absorption and scattering processes computed for two polarization modes. We include the effects of vacuum polarization in our calculations. We present temperature profiles and the angle-, photon energy-, and polarization-dependent emerging intensity for a range of magnetic field strengths and effective temperatures of the atmospheres. We find that for B < 10^14 G, the emerging spectra are bluer than the blackbody corresponding to the effective temperature, T_eff, with modified Planckian shapes due to the photon-energy dependence of the magnetic opacities. However, vacuum polarization significantly modifies the spectra for B~10^15 G, giving rise to power-law tails at high photon energies. The angle-dependence (beaming) of the emerging intensity has two maxima: a narrow (pencil) peak at small angles (<5 degrees) with respect to the normal and a broad maximum (fan beam) at intermediate angles (~20-60 degrees). The relative importance and the opening angle of the radial beam decreases strongly with increasing magnetic field strength and decreasing photon energy. We finally compute a T_eff-T_c relation for our models, where T_c is the local color temperature of the spectrum emerging from the neutron star surface, and find that T_c/T_eff ranges between 1.1-1.8. We discuss the implications of our results for various thermally emitting neutron star models.

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.

General relativistic bending of light dramatically alters the variability of X-ray emission originating from the surfaces of ultramagnetic neutron stars. We construct radiative equilibrium models of such strongly magnetic cooling neutron stars with light-element atmospheres to compute the angle- and energy-dependent intensity emerging from their surfaces and find that the beaming of surface emission is predominantly non-radial. The combination of this radiation pattern with the calculations of light bending yields pulse amplitudes that vary non-monotonically with the neutron star compactness and the size of the emitting region. The significant suppression of the pulse amplitude for large emitting areas provides very strong constraints on the mechanisms that can simultaneously produce high periodic variability and X-ray luminosity. We apply these results to the thermally-emitting magnetar models of anomalous X-ray pulsars (AXPs), which are bright slowly-rotating X-ray sources with large pulse amplitudes. We use the observed fluxes and pulse amplitudes for all known AXPs and show that thermal emission from two antipodal regions on their surfaces, as predicted by some magnetar models, is inconsistent with these observed properties.

We consider the X-ray emission due to bremsstrahlung processes from hot, low radiative-efficiency accretion flows around supermassive and galactic black holes. We calculate surface brightness profiles and Michelson visibility functions for a range of density profiles, rho ~ r^(-3/2+p), with 0 < p < 1, to allow for the presence of outflows. We find that although the 1 keV emitting region in these flows can always extend up to 10^6 Schwarzschild radii (R_S), their surface brightness profiles and visibility functions are strongly affected by the specific density profile. The advection-dominated solutions with no outflows (p=0) lead to centrally peaked profiles with characteristic sizes of only a few tens of R_S. Solutions with strong outflows (p~1) lead to flat intensity profiles with significantly larger characteristic sizes of up to 10^6 R_S. This implies that low luminosity galactic nuclei, such as M87, may appear as extended X-ray sources when observed with current X-ray imaging instruments. We show that X-ray brightness profiles and their associated visibility functions may be powerful probes for determining the relevant mode of accretion and, in turn, the properties of hot accretion flows. We discuss the implications of our results for observations with the Chandra X-ray Observatory and the planned X-ray interferometer MAXIM.

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.

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.