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.

Comments: ApJ, in press. Revised version includes more expanded analyses and discussions of comparison to previous results

Similar Publications

The observed evolution of the broad-band spectral energy distribution (SED) in NS X-ray Nova Aql X-1 during the rise phase of a bright FRED-type outburst in 2013 can be understood in the framework of thermal emission from unstationary accretion disc with temperature radial distribution transforming from a single-temperature blackbody emitting ring into the multi-colour irradiated accretion disc. SED evolution during the hard to soft X-ray state transition looks curious, as it can not be repro- duced by the standard disc irradiation model with a single irradiation parameter for NUV, Optical and NIR spectral bands. NIR (NUV) band is correlated with soft (hard) X-ray flux changes during the state transition interval, respectively. Read More

We discuss a new method for unveiling the possible blazar AGN nature among the numerous population of Unassociated Gamma-ray sources (UGS) in the Fermi catalogues. Our tool relies on positional correspondence of the Fermi object with X-ray sources (mostly from Swift-XRT), correlated with other radio, IR and optical data in the field. We built a set of Spectral Energy Distributions (SED) templates representative of the various blazar classes, and we quantitatively compared them to the observed multi-wavelength flux density data for all Swift-XRT sources found within the Fermi error-box, by taking advantage of some well-recognised regularities in the broad-band spectral properties of the objects. Read More

I propose a scenario where the majority of the progenitors of type IIb supernovae (SNe IIb) lose most of their hydrogen-rich envelope during a grazing envelope evolution (GEE). In the GEE the orbital radius of the binary system is about equal to the radius of the giant star, and the more compact companion accretes mass through an accretion disk. The accretion disk is assumed to launch two opposite jets that efficiently remove gas from the envelope along the orbit of the companion. Read More

The compact and, with about 4.3+-0.3 million solar masses, very massive object located at the center of the Milky Way is currently the very best candidate for a supermassive black hole (SMBH) in our immediate vicinity. Read More

Gamma-ray bursts (GRBs) are powerful probes of early stars and galaxies, during and potentially even before the era of reionization. Although the number of GRBs identified at z>6 remains small, they provide a unique window on typical star-forming galaxies at that time, and thus are complementary to deep field observations. We report the identification of the optical drop-out afterglow of Swift GRB 120923A in near-infrared Gemini-North imaging, and derive a redshift of z=7. Read More

We develop a general approach to analytically calculate the shift $\Delta\mathcal{O}$ that observable quantities $\mathcal{O}$ pertaining binary systems experience because of the orbital perturbations affecting the motion around their center of mass. In particular, we deal with the disturbance $\Delta\delta\tau_\textrm{p}$ suffered by the otherwise periodic change $\delta\tau_\textrm{p}$ of the ratio of the projection $z_\textrm{p}$ of the barycentric orbital motion of the binary's visible member p along the line of sight $\boldsymbol{\hat{e}}_3$ to the speed of light $c$ because of some Newtonian and post-Newtonian (pN) non-central accelerations (mass quadrupole $Q_2$, 1pN static and stationary field). We apply our results to the double pulsar system PSR J0737-3039A/B and to the hypothetic scenario involving a pulsar orbiting the Supermassive Black Hole in in the Galactic Center at Sgr A$^\ast$. Read More

The first detected gravitational wave signal, GW150914, was produced by the coalescence of a stellar-mass binary black hole. Along with the subsequent detection of GW151226 and the candidate event LVT151012, this gives us evidence for a population of black hole binaries with component masses in the tens of solar masses. As detector sensitivity improves, this type of source is expected to make a large contribution to the overall number of detections, but this type of source has received little attention compared to binary neutron star systems in studies of projected network performance. Read More

We study rapidly rotating hybrid stars with the Dyson-Schwinger model for quark matter and the Brueckner-Hartree-Fock many-body theory with realistic two-body and three-body forces for nuclear matter. We determine the maximum gravitational mass, equatorial radius, and rotation frequency of stable stellar configurations by considering the constraints of the Keplerian limit and the secular axisymmetric instability, and compare with observational data. We also discuss the rotational evolution for constant baryonic mass, and find a spinup phenomenon for supramassive stars before they collapse to black holes. Read More

Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. Read More

Authors: Elena Pian1
Affiliations: 1INAF, IASF Bologna and Scuola Normale Superiore di Pisa, Italy

The INTEGRAL mission has played a major role in blazar science, thanks to its sensitive coverage of a spectral region (3-100 keV) that is critical for this type of sources, to its flexibility of scheduling and to the large field of view of its cameras. A number of flat-spectrum radio quasars (up to z ~ 3) and BL Lac objects were observed by INTEGRAL together with facilities at all wavelengths. These results have advanced our knowledge of blazars from a physical and cosmological point of view. Read More