Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGW170817: Measurements of Neutron Star Radii and Equation of State
93
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function $p(rho)$ of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as $R_1=10.8^{+2.0}_{-1.7}$ km for the heavier star and $R_2= 10.7^{+2.1}_{-1.5}$ km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than $1.97 ,M_odot$ as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain $R_1= 11.9^{+1.4}_{-1.4}$ km and $R_2= 11.9^{+1.4}_{-1.4}$ km at the 90% credible level. Finally, we obtain constraints on $p(rho)$ at supranuclear densities, with pressure at twice nuclear saturation density measured at $3.5^{+2.7}_{-1.7}times 10^{34} ,mathrm{dyn}/mathrm{cm}^{2}$ at the 90% level.
rate research
Read More
We explore in a parameterized manner a very large range of physically plausible equations of state (EOSs) for compact stars for matter that is either purely hadronic or that exhibits a phase transition. In particular, we produce two classes of EOSs with and without phase transitions, each containing one million EOSs. We then impose constraints on the maximum mass, ($M < 2.16 M_{odot}$), and on the dimensionless tidal deformability ($tilde{Lambda} <800$) deduced from GW170817, together with recent suggestions of lower limits on $tilde{Lambda}$. Exploiting more than $10^9$ equilibrium models for each class of EOSs, we produce distribution functions of all the stellar properties and determine, among other quantities, the radius that is statistically most probable for any value of the stellar mass. In this way, we deduce that the radius of a purely hadronic neutron star with a representative mass of $1.4,M_{odot}$ is constrained to be $12.00!<!R_{1.4}/{rm km}!<!13.45$ at a $2$-$sigma$ confidence level, with a most likely value of $bar{R}_{1.4}=12.39,{rm km}$; similarly, the smallest dimensionless tidal deformability is $tilde{Lambda}_{1.4}!>!375$, again at a $2$-$sigma$ level. On the other hand, because EOSs with a phase transition allow for very compact stars on the so-called `twin-star branch, small radii are possible with such EOSs although not probable, i.e. $8.53!<!R_{1.4}/{rm km}!<!13.74$ and $bar{R}_{1.4}=13.06,{rm km}$ at a $2$-$sigma$ level, with $tilde{Lambda}_{1.4}!>!35.5$ at a $3$-$sigma$ level. Finally, since these EOSs exhibit upper limits on $tilde{Lambda}$, the detection of a binary with total mass of $3.4,M_{odot}$ and $tilde{Lambda}_{1.7}!>!461$ can rule out twin-star solutions.
Observations of the properties of multiple coalescing neutron stars will simultaneously provide insight into neutron star mass and spin distribution, the neutron star merger rate, and the nuclear equation of state. Not all merging binaries containing neutron stars are expected to be identical. Plausible sources of diversity in these coalescing binaries can arise from a broad or multi-peaked NS mass distribution; the effect of different and extreme NS natal spins; the possibility of NS-BH mergers; or even the possibility of phase transitions, allowing for NS with similar mass but strongly divergent radius. In this work, we provide a concrete algorithm to combine all information obtained from GW measurements into a joint constraint on the NS merger rate, the distribution of NS properties, and the nuclear equation of state. Using a concrete example, we show how biased mass distribution inferences can significantly impact the recovered equation of state, even in the small-$N$ limit. With the same concrete example, we show how small-$N$ observations could identify a bimodal mass and spin distribution for merging NS simultaneously with the EOS. Our concordance approach can be immediately generalized to incorporate other observational constraints.
We present results from three-dimensional general relativistic simulations of binary neutron star coalescences and mergers using public codes. We considered equal mass models where the baryon mass of the two Neutron Stars (NS) is $1.4M_{odot}$, described by four different equations of state (EOS) for the cold nuclear matter (APR4, SLy, H4, and MS1; all parametrized as piecewise polytropes). We started the simulations from four different initial interbinary distances ($40, 44.3, 50$, and $60$ km), including up to the last 16 orbits before merger. That allows to show the effects on the gravitational wave phase evolution, radiated energy and angular momentum due to: the use of different EOSs, the orbital eccentricity present in the initial data and the initial separation (in the simulation) between the two stars. Our results show that eccentricity has a major role in the discrepancy between numerical and analytical waveforms until the very last few orbits, where tidal effects and missing high-order post-Newtonian coefficients also play a significant role. We test different methods for extrapolating the gravitational wave signal extracted at finite radii to null infinity. We show that an effective procedure for integrating the Newman-Penrose $psi_4$ signal to obtain the gravitational wave strain $h$ is to apply a simple high-pass digital filter to $h$ after a time domain integration, where only the two physical motivated integration constants are introduced. That should be preferred to the more common procedures of introducing additional integration constants, integrating in the frequency domain or filtering $psi_4$ before integration.
Recently exploratory studies were performed on the possibility of constraining the neutron star equation of state (EOS) using signals from coalescing binary neutron stars, or neutron star-black hole systems, as they will be seen in upcoming advanced gravitational wave detectors such as Advanced LIGO and Advanced Virgo. In particular, it was estimated to what extent the combined information from multiple detections would enable one to distinguish between different equations of state through hypothesis ranking or parameter estimation. Under the assumption of zero neutron star spins both in signals and in template waveforms and considering tidal effects to 1 post-Newtonian (1PN) order, it was found that O(20) sources would suffice to distinguish between a hard, moderate, and soft equation of state. Here we revisit these results, this time including neutron star tidal effects to the highest order currently known, termination of gravitational waveforms at the contact frequency, neutron star spins, and the resulting quadrupole-monopole interaction. We also take the masses of neutron stars in simulated sources to be distributed according to a relatively strongly peaked Gaussian, as hinted at by observations, but without assuming that the data analyst will necessarily have accurate knowledge of this distribution for use as a mass prior. We find that especially the effect of the latter is dramatic, necessitating many more detections to distinguish between different EOS and causing systematic biases in parameter estimation, on top of biases due to imperfect understanding of the signal model pointed out in earlier work. This would get mitigated if reliable prior information about the mass distribution could be folded into the analyses.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
