BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Dimitrios Psaltis

Dimitrios Psaltis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.IMastro-ph.GAnucl-thastro-ph.COastro-ph.SRastro-phastro-ph.EPhep-phhep-th

Catalog footprint

What is connected

75works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

75 published item(s)

preprint2026arXiv

Limits to Extracting Neutron-Star Physics Constraints from NICER Pulse Profiles

Modeling energy-dependent X-ray pulse profiles from rotation-powered millisecond pulsars observed with NICER has emerged as a promising avenue for measuring neutron star radii and probing the equation of state of cold, ultra-dense matter. However, pulse profile models have often required an unwieldy number of parameters to account for complex surface emission geometries, introducing the risk of overfitting and degeneracies. To explore the number of model parameters that can be inferred uniquely, we perform a quantitative assessment of the information content in X-ray pulse profiles by applying Fourier methods. We determine the number of independent observables that can be reliably extracted from the pulse shapes, as well as from complementary X-ray spectral data obtained with XMM-Newton, for key NICER targets. Our analysis provides a framework for evaluating the match between model complexity and data constraints. It also demonstrates the importance of incorporating in the model the pulsed components of the magnetospheric non-thermal emission, which may often contribute significantly to the observed spectra. Our results highlight limitations in previous inferences of neutron-star radii from NICER observations, which may have incorporated model complexity not supported by the data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

MeqSilhouette v2: Spectrally-resolved polarimetric synthetic data generation for the Event Horizon Telescope

We present MeqSilhouette v2.0 (MeqSv2), a fully polarimetric, time-and frequency-resolved synthetic data generation software for simulating millimetre (mm) wavelength very long baseline interferometry (VLBI) observations with heterogeneous arrays. Synthetic data are a critical component in understanding real observations, testing calibration and imaging algorithms, and predicting performance metrics of existing or proposed sites. MeqSv2 applies physics-based instrumental and atmospheric signal corruptions constrained by empirically-derived site and station parameters to the data. The new version is capable of applying instrumental polarization effects and various other spectrally-resolved effects using the Radio Interferometry Measurement Equation (RIME) formalism and produces synthetic data compatible with calibration pipelines designed to process real data. We demonstrate the various corruption capabilities of MeqSv2 using different arrays, with a focus on the effect of complex bandpass gains on closure quantities for the EHT at 230 GHz. We validate the frequency-dependent polarization leakage implementation by performing polarization self-calibration of synthetic EHT data using PolSolve. We also note the potential applications for cm-wavelength VLBI array analysis and design and future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Brightness Asymmetry of Black Hole Images as a Probe of Observer Inclination

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimized Statistical Approach for Comparing Multi-Messenger Neutron Star Data

The neutron star equation of state 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 often mapped to a common domain for comparison with one another or are used to constrain the predictions of theoretical equations of state. We explore here the statistical biases that can arise when such multi-messenger data are compared or combined across different domains. We find that placing Bayesian priors individually in each domain of measurement can lead to biased constraints. We present a new prescription for defining Bayesian priors consistently across different experiments, which will allow for robust cross-domain comparisons. Using the first two binary neutron star mergers as an example, we show that a uniform prior in the tidal deformability can produce inflated evidence for large radii, while a uniform prior in the radius points towards smaller radii. Finally, using this new prescription, we provide a status update on multi-messenger constraints on the neutron star radius.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Parametric model for the shapes of black-hole shadows in non-Kerr spacetimes

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Continuum-Fitting the X-ray Spectra of Tidal Disruption Events

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Discretization and Filtering Effects on Black Hole Images Obtained with the Event Horizon Telescope

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interferometric Closure Phase Uncertainties in the Low Signal-to-Noise Ratio Regime

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting Theoretical Models to Event Horizon Telescope Observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Testing General Relativity with the Event Horizon Telescope

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bayesian techniques for comparing time-dependent GRMHD simulations to variable Event Horizon Telescope observations

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From Neutron Star Observables to the Equation of State. I. An Optimal Parametrization

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Measuring the neutron star equation of state using X-ray timing

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Model-Independent Inference of Neutron Star Radii from Moment of Inertia Measurements

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Particle Acceleration and the Origin of X-ray Flares in GRMHD simulations of Sgr A*

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Persistent Asymmetric Structure of Sagittarius A* on Event Horizon Scales

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Systematic Uncertainties in the Spectroscopic Measurements of Neutron-Star Masses and Radii from Thermonuclear X-ray Bursts. III. Absolute Flux Calibration

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing General Relativity with Accretion-Flow Imaging of Sgr A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A General Relativistic Null Hypothesis Test with Event Horizon Telescope Observations of the black-hole shadow in Sgr A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Quantitative Test of the No-Hair Theorem with Sgr A* using stars, pulsars, and the Event Horizon Telescope

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Data Selection Criteria for Spectroscopic Measurements of Neutron Star Radii with X-ray Bursts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of Spot Size on Neutron-Star Radius Measurements from Pulse Profiles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Variability and mm/IR flares in GRMHD Models of Sgr A* from Strong-Field Gravitational Lensing

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HERO: A 3D General Relativistic Radiative Postprocessor for Accretion Discs around Black Holes

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HEROIC: 3D General Relativistic Radiative Postprocessor with Comptonization for Black Hole Accretion Discs

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Linking Tests of Gravity On All Scales: from the Strong-Field Regime to Cosmology

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring Neutron Star Radii via Pulse Profile Modeling with NICER

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resolved Magnetic-Field Structure and Variability Near the Event Horizon of Sagittarius A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shapiro Delays at the Quadrupole Order for Tests of the No-Hair Theorem Using Pulsars around Spinning Black Holes

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistics of Measuring Neutron Star Radii: Assessing A Frequentist and A Bayesian Approach

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Dense Matter Equation of State from Neutron Star Radius and Mass Measurements

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Hard X-Ray Power-Law Spectral Cutoff in Centaurus X-4

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A High-Frequency Doppler Feature in the Power Spectra of Simulated GRMHD Black Hole Accretion Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Event-Horizon-Telescope Evidence for Alignment of the Black Hole in the Center of the Milky Way with the Inner Stellar Disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Imaging an Event Horizon: Mitigation of Scattering Toward Sagittarius A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Prospects for Measuring Neutron-Star Masses and Radii with X-Ray Pulse Profile Modeling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pulse Profiles from Spinning Neutron Stars in the Hartle-Thorne Approximation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rotational Corrections to Neutron-Star Radius Measurements from Thermal Spectra

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Power of Imaging: Constraining the Plasma Properties of GRMHD Simulations using EHT Observations of Sgr A*

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

GRay: a Massively Parallel GPU-Based Code for Ray Tracing in Relativistic Spacetimes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-Angular-Resolution and High-Sensitivity Science Enabled by Beamformed ALMA

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Narrow Atomic Features from Rapidly Spinning Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relations Between Neutron-Star Parameters in the Hartle-Thorne Approximation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing the No-Hair Theorem with Event Horizon Telescope Observations of Sagittarius A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum. IV. Relativistically Broadened Iron Lines

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hot and Energetic Universe: End points of stellar evolution

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. II. Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Deviation of Stellar Orbits from Test Particle Trajectories Around Sgr A* Due to Tides and Winds

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Masses of Nearby Supermassive Black Holes with Very-Long Baseline Interferometry

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling Collisional Cascades In Debris Disks: Steep Dust-Size Distributions

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling Collisional Cascades In Debris Disks: The Numerical Method

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Mass Distribution and Birth Masses of Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Special and General Relativistic Effects in Galactic Rotation Curves

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Systematic Uncertainties in the Spectroscopic Measurements of Neutron-Star Masses and Radii from Thermonuclear X-ray Bursts. I. Apparent Radii

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Systematic Uncertainties in the Spectroscopic Measurements of Neutron-Star Masses and Radii from Thermonuclear X-ray Bursts. II. Eddington Limit

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Metric for Rapidly Spinning Black Holes Suitable for Strong-Field Tests of the No-Hair Theorem

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Influence of Gas Dynamics on Measuring the Properties of the Black Hole in the Center of the Milky Way with Stellar Orbits and Pulsars

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Ray-Tracing Algorithm for Spinning Compact Object Spacetimes with Arbitrary Quadrupole Moments. I. Quasi-Kerr Black Holes

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining Parity Violation in Gravity with Measurements of Neutron-Star Moments of Inertia

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 ξ^{1/4} < 5 km, which is at least three-orders of magnitude stronger than the previous strongest bound.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Neutron Stars in f(R) Gravity with Perturbative Constraints

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reconstructing the Neutron-Star Equation of State from Astrophysical Measurements

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sgr A*: The Optimal Testbed of Strong-Field Gravity

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum: I. Properties of a Quasi-Kerr Spacetime

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum: II. Black-Hole Images

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing the No-Hair Theorem with Observations in the Electromagnetic Spectrum. III. Quasi-Periodic Variability

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing the No-Hair Theorem with Observations of Black Holes in the Electromagnetic Spectrum

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Black Hole Mass Distribution in the Galaxy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Incompatibility of Rapid Rotation with Narrow Photospheric X-ray Lines in EXO 0748-676

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Massive Pulsar PSR J1614-2230: Linking Quantum Chromodynamics, Gamma-ray Bursts, and Gravitational Wave Astronomy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Two approaches to testing general relativity in the strong-field regime

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Probes and Tests of Strong-Field Gravity with Observations in the Electromagnetic Spectrum

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Mass and Radius of the Neutron Star in EXO 1745-248

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.

0 citations0 reviewsNo code linkedOpen access