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

Feryal Ozel

Feryal Ozel 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-qcnucl-thastro-ph.GAastro-ph.SRastro-phastro-ph.IMastro-ph.COhep-phhep-th

Catalog footprint

What is connected

52works
10topics
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

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

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

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
preprint2019arXiv

Measurement of the nuclear symmetry energy parameters from gravitational wave events

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

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

Probing X-ray Absorption and Optical Extinction in the Interstellar Medium Using Chandra Observations of Supernova Remnants

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radio Constraints on Long-Lived Magnetar Remnants in Short Gamma-Ray Bursts

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

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

The millisecond pulsar mass distribution: Evidence for bimodality and constraints on the maximum neutron star mass

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

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

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

Mapping the Surface of the Magnetar 1E 1048.1-5937 in Outburst and Quiescence Through Phase Resolved X-ray Spectroscopy

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

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

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

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

The mass and the radius of the neutron star in the transient low mass X-ray binary SAX J1748.9-2021

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

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

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: 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 Cooling Trend of SGR 0526-66

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

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

Radio Synchrotron Emission from a Bow Shock around the Gas Cloud G2 Heading toward the Galactic Center

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Surface Emission from Neutron Stars and Implications for the Physics of their Interiors

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

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 Magnetar Strength Surface Magnetic Field for the Slowly Spinning Down SGR 0418+5729

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Long Term Radiative Behavior of SGR 1900+14

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Mass and Radius of the Neutron Star in the Bulge Low-Mass X-ray Binary KS 1731-260

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Astrophysical Measurement of the Equation of State of Neutron Star Matter

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

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

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

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 Distance, Mass, and Radius of the Neutron Star in 4U 1608-52

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

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 Mass and Radius of the Neutron Star in 4U 1820-30

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

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

The Relation Between Optical Extinction and Hydrogen Column Density in the Galaxy

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

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
preprint2008arXiv

What Makes an Accretion-Powered Millisecond Pulsar?

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

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Hardness-Intensity Correlations in Magnetar Afterglows

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

0 citations0 reviewsNo code linkedOpen access