Tidal Deformability Of Neutron Stars Research Articles

Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.

Tidal deformability of neutron stars with realistic nuclear energy density functionals

We investigate the constraints on the mass and radius of neutron stars by considering the tidal deformability in the merge of neutron star binaries. In order to extract the most reliable range of uncertainty from theory, we employ models based upon the Skyrme force and density functional theory and select models that are consistent with empirical data of finite nuclei, measured properties of nuclear matter around the saturation density, and observation of the maximum mass of neutron stars. From the selected models, we calculate the Love number $k_2$, dimensionless tidal deformability $\Lambda$, and mass-weighted deformability $\tilde{\Lambda}$ in the binary system. We find that all the models considered in this work give $\tilde{\Lambda}$ less than 800 which is the upper limit obtained from the measurement of GW170817. However, the model dependence of tidal deformability is manifest such that our results on the tidal deformability exhibit critical sensitivity to the size of neutron stars.

Constraints on the Moment of Inertia of PSR J0737-3039A from GW170817

Abstract Continued observation of PSR J0737-3039, the double pulsar, is expected to yield a precise determination of its primary component’s moment of inertia in the next few years. Since the moment of inertia depends sensitively on the neutron star’s internal structure, such a measurement will constrain the equation of state of ultra-dense matter, which is believed to be universal. Independent equation-of-state constraints have already been established by the gravitational-wave measurement of neutron star tidal deformability in GW170817. Here, using well-known universal relations among neutron star observables, we translate the reported 90% credible bounds on tidal deformability into a direct constraint, , on the moment of inertia of PSR J0737-3039A. Should a future astrophysical measurement of I ⋆ disagree with this prediction, it could indicate a breakdown in the universality of the neutron star equation of state.

Astrophysical constraints on a parametric equation of state for neutron-rich nucleonic matter

Extracting the equation of state (EOS) and symmetry energy of dense neutron-rich matter from astrophysical observations is a long-standing goal of nuclear astrophysics. To facilitate the realization of this goal, the feasibility of using an explicitly isospin-dependent parametric EOS for neutron star matter was investigated recently in [1–3]. In this contribution, in addition to outlining the model framework and summarizing the most important findings from [1–3], we report a few new results regarding constraining parameters characterizing the high-density behavior of nuclear symmetry energy. In particular, the constraints on the pressure of neutron star matter extracted from combining the X-ray observations of the neutron star radius, the minimum–maximum mass $$M=2.01$$ $$\hbox {M}_\odot $$ , and causality condition agree very well with those extracted from analyzing the tidal deformability data by the LIGO + Virgo Collaborations. The limitations of using the radius and/or tidal deformability of neutron stars to constrain the high-density nuclear symmetry energy are discussed.

GW170817: Constraining the nuclear matter equation of state from the neutron star tidal deformability

Constraints set on key parameters of the nuclear matter equation of state (EoS) by the values of the tidal deformability, inferred from GW170817, are examined by using a diverse set of relativistic and non-relativistic mean field models. These models are consistent with bulk properties of finite nuclei as well as with the observed lower bound on the maximum mass of neutron star $\sim 2 ~ {\rm M}_\odot$. The tidal deformability shows a strong correlation with specific linear combinations of the isoscalar and isovector nuclear matter parameters associated with the EoS. Such correlations suggest that a precise value of the tidal deformability can put tight bounds on several EoS parameters, in particular, on the slope of the incompressibility and the curvature of the symmetry energy. The tidal deformability obtained from the GW170817 and its UV/optical/infrared counterpart sets the radius of a canonical $1.4~ {\rm M}_{\odot}$ neutron star to be $11.82\leqslant R_{1.4}\leqslant13.72$ km.

Comparing two models for measuring the neutron star equation of state from gravitational-wave signals

Observations of gravitational-wave signals from binary neutron-star mergers, like GW170817, can be used to constrain the neutron-star equation of state (EoS). One method involves modeling the EoS and measuring the model parameters through Bayesian inference. A previous study has demonstrated the effectiveness of using a phenomenologically parameterized piecewise polytrope to extract constraining information from a simulated population of binary neutron-star mergers. Despite its advantages compared to more traditional methods of measuring the tidal deformability of neutron stars, notable deficiencies arise when using this EoS model. In this work, we describe in detail the implementation of a model built from a spectral decomposition of the adiabatic index that was used by the LIGO-Virgo Collaboration in to constrain the neutron star EoS from GW170817. We demonstrate its overall consistency in recovering the neutron star EoS from a simulated signal to the piecewise polytropic implementation used in and explain any differences that arise. We find that both models recover consistent tidal information from the simulate signals with tightest constraints on the EoS around twice nuclear saturation density. As expected, the statistical error that arises in the piecewise polytropic representation near the fixed joining densities is greatly reduced by using the spectral model. In addition, we find that our choice of prior can have a dominant effect on EoS constraints

Neutron Star Tidal Deformabilities Constrained by Nuclear Theory and Experiment

We confront observational data from gravitational wave event GW170817 with microscopic modeling of the cold neutron star equation of state. We develop and employ a Bayesian statistical framework that enables us to implement constraints on the equation of state from laboratory measurements of nuclei and state-of-the-art chiral effective field theory methods. The energy density functionals constructed from the posterior probability distributions are then used to compute consistently the neutron star equation of state from the outer crust to the inner core, assuming a composition consisting of protons, neutrons, electrons, and muons. In contrast to previous studies, we find that the 95% credibility range of predicted neutron star tidal deformabilities (136<Λ<519) for a 1.4 solar-mass neutron star is already consistent with the upper bound deduced from observations of the GW170817 event. However, we find that lower bounds on the neutron star tidal deformability will very strongly constrain microscopic models of the dense matter equation of state. We also demonstrate a strong correlation between the neutron star tidal deformability and the pressure of beta-equilibrated matter at twice saturation density.

Gravitational Waves and Tidal Deformability of Neutron Stars

Gravitational Waves and Tidal Deformability of Neutron Stars

Measuring the neutron star tidal deformability with equation-of-state-independent relations and gravitational waves

Gravitational wave measurements of binary neutron star coalescences offer information about the properties of the extreme matter that comprises the stars. Despite our expectation that all neutron stars in the Universe obey the same equation of state, i.e. the properties of the matter that forms them are universal, current tidal inference analyses treat the two bodies as independent. We present a method to measure the effect of tidal interactions in the gravitational wave signal -- and hence constrain the equation of state -- that assumes that the two binary components obey the same equation of state. Our method makes use of a relation between the tidal deformabilities of the two stars given the ratio of their masses, a relation that has been shown to only have a weak dependance on the equation of state. We use this relation to link the tidal deformabilities of the two stars in a realistic parameter inference study while simultaneously marginalizing over the error in the relation. This approach incorporates more physical information into our analysis, thus leading to a better measurement of tidal effects in gravitational wave signals. Through simulated signals we estimate that uncertainties in the measured tidal parameters are reduced by a factor of at least 2 -- and in some cases up to 10 -- depending on the equation of state and mass ratio of the system.

On the ambiguity in relativistic tidal deformability

The LIGO collaboration recently reported the first gravitational-wave constraints on the tidal deformability of neutron stars. I discuss an inherent ambiguity in the notion of relativistic tidal deformability that, while too small to affect the present measurement, may become important in the future. I propose a new way to understand the ambiguity and discuss future prospects for reliably linking observed gravitational waveforms to compact object microphysics.

Upper limit set by causality on the tidal deformability of a neutron star

A principal goal of gravitational-wave astronomy is to constrain the neutron star equation of state (EOS) by measuring the tidal deformability of neutron stars. The tidally induced departure of the waveform from that of point-particle (or spinless binary black hole (BBH)) increases with the stiffness of the EOS. We show that causality (the requirement that the speed of sound is less than the speed of light for a perfect fluid satisfying a one-parameter equation of state) places an upper bound on tidal deformability as a function of mass. Like the upper mass limit, the limit on deformabity is obtained by using an EOS with $v_{sound} = c$ for high densities and matching to a low density (candidate) EOS at a matching density of order nuclear saturation density. We use these results and those of [B.D. Lackey et al., Phys. Rev. D 89, 043009 (2014)] to estimate the resulting upper limit on the gravitational-wave phase shift of a black hole-neutron star (BHNS) binary relative to a BBH. Even for assumptions weak enough to allow an maximum mass of $4\ M_\odot$ (a match at nuclear saturation density to an unusually stiff low-density candidate EOS), the upper limit on dimensionless tidal deformability is stringent. It leads to a still more stringent estimated upper limit on the maximum tidally induced phase shift prior to merger. We comment in an appendix on the relation between causality, the condition $v_{sound} < c$, and the condition $dp/d\epsilon < 1$ for the effective EOS governing the equilibrium star.

Measuring neutron star tidal deformability with Advanced LIGO: A Bayesian analysis of neutron star-black hole binary observations

The pioneering discovery of gravitational waves (GWs) by Advanced LIGO has ushered us into an era of observational GW astrophysics. Compact binaries remain the primary target sources for GW observation, of which neutron star-black hole (NSBH) binaries form an important subset. GWs from NSBH sources carry signatures of (a) the tidal distortion of the neutron star by its companion black hole during inspiral, and (b) its potential tidal disruption near merger. In this paper, we present a Bayesian study of the measurability of neutron star tidal deformability ${\mathrm{\ensuremath{\Lambda}}}_{\mathrm{NS}}\ensuremath{\propto}(R/M{)}_{\mathrm{NS}}^{5}$ using observation(s) of inspiral-merger GW signals from disruptive NSBH coalescences, taking into account the crucial effect of black hole spins. First, we find that if nontidal templates are used to estimate source parameters for an NSBH signal, the bias introduced in the estimation of nontidal physical parameters will only be significant for loud signals with signal-to-noise ratios greater than $\ensuremath{\simeq}30$. For similarly loud signals, we also find that we can begin to put interesting constraints on ${\mathrm{\ensuremath{\Lambda}}}_{\mathrm{NS}}$ (factor of 1--2) with individual observations. Next, we study how a population of realistic NSBH detections will improve our measurement of neutron star tidal deformability. For an astrophysically likely population of disruptive NSBH coalescences, we find that 20--35 events are sufficient to constrain ${\mathrm{\ensuremath{\Lambda}}}_{\mathrm{NS}}$ within $\ifmmode\pm\else\textpm\fi{}25%--50%$, depending on the neutron star equation of state. For these calculations we assume that LIGO will detect black holes with masses within the astrophysical mass gap. In case the mass gap remains preserved in NSBHs detected by LIGO, we estimate that approximately 25% additional detections will furnish comparable ${\mathrm{\ensuremath{\Lambda}}}_{\mathrm{NS}}$ measurement accuracy. In both cases, we find that it is the loudest 5--10 events that provide most of the tidal information, and not the combination of tens of low-SNR events, thereby facilitating targeted numerical-GR follow-ups of NSBHs. We find these results encouraging, and recommend that an effort to measure ${\mathrm{\ensuremath{\Lambda}}}_{\mathrm{NS}}$ be planned for upcoming NSBH observations with the LIGO-Virgo instruments.

Dynamical tides in general relativity: Effective action and effective-one-body Hamiltonian

Tidal effects have an important impact on the late inspiral of compact binary systems containing neutron stars. Most current models of tidal deformations of neutron stars assume that the tidal bulge is directly related to the tidal field generated by the companion, with a constant response coefficient. However, if the orbital motion approaches a resonance with one of the internal modes of the neutron star, this adiabatic description of tidal effects starts to break down, and the tides become dynamical. In this paper, we consider dynamical tides in general relativity due to the quadrupolar fundamental oscillation mode of a neutron star. We devise a description of the effects of the neutron star's finite size on the orbital dynamics based on an effective point-particle action augmented by dynamical quadrupolar degrees of freedom. We analyze the post-Newtonian and test-particle approximations of this model and incorporate the results into an effective-one-body Hamiltonian. This enables us to extend the description of dynamical tides over the entire inspiral. We demonstrate that dynamical tides give a significant enhancement of matter effects compared to adiabatic tides, at least for neutron stars with large radii and for low mass-ratio systems, and should therefore be included in accurate models for gravitational-wave data analysis.

Binary Love relations

When in a tight binary, the mutual tidal deformations of neutron stars get imprinted onto observables, encoding information about their internal structure at supranuclear densities and gravity in the extreme-gravity regime. Gravitational wave (GW) observations of their late binary inspiral may serve as a tool to extract the individual tidal deformabilities, but this is made difficult by degeneracies between them in the GW model. We here resolve this problem by discovering approximately equation-of-state (EoS)-insensitive relations between dimensionless combinations of the individual tidal deformabilities. We show that these relations break degeneracies in the GW model, allowing for the accurate extraction of both deformabilities. Such measurements can be used to better differentiate between EoS models, and improve tests of general relativity and cosmology.

On the validity of the adiabatic approximation in compact binary inspirals

Using a semi-analytical approach recently developed to model the tidal deformations of neutron stars in inspiralling compact binaries, we study the dynamical evolution of the tidal tensor, which we explicitly derive at second post-Newtonian order, and of the quadrupole tensor. Since we do not assume a priori that the quadrupole tensor is proportional to the tidal tensor, i.e. the so called "adiabatic approximation", our approach enables us to establish to which extent such approximation is reliable. We find that the ratio between the quadrupole and tidal tensors (i.e., the Love number) increases as the inspiral progresses, but this phenomenon only marginally affects the emitted gravitational waveform. We estimate the frequency range in which the tidal component of the gravitational signal is well described using the stationary phase approximation at next-to-leading post-Newtonian order, comparing different contributions to the tidal phase. We also derive a semi-analytical expression for the Love number, which reproduces within a few percentage points the results obtained so far by numerical integrations of the relativistic equations of stellar perturbations.

Tidal deformations of neutron stars: The role of stratification and elasticity

We discuss the response of neutron stars to the tidal interaction in a compact binary system, as encoded in the Love number associated with the induced deformation. This problem is of interest for gravitational-wave astronomy as there may be a detectable imprint on the signal from the late stages of binary coalescence. Previous work has focussed on simple barotropic neutron star models, providing an understanding of the role of the stellar compactness and overall density profile. We add realism to the discussion by developing the framework required to model stars with varying composition and an elastic crust. These effects are not expected to be significant for the next generation of detectors but it is nevertheless useful to be able to quantify them. Our results show that (perhaps surprisingly) internal stratification has no impact whatsoever on the Love number. We also show that crust elasticity provides a (predictably) small correction to existing models.

Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral

The early part of the gravitational wave signal of binary neutron star inspirals can potentially yield robust information on the nuclear equation of state. The influence of a star's internal structure on the waveform is characterized by a single parameter: the tidal deformability lambda, which measures the star's quadrupole deformation in response to the companion's perturbing tidal field. We calculate lambda for a wide range of equations of state and find that the value of lambda spans an order of magnitude for the range of equation of state models considered. An analysis of the feasibility of discriminating between neutron star equations of state with gravitational wave observations of the early part of the inspiral reveals that the measurement error in lambda increases steeply with the total mass of the binary. Comparing the errors with the expected range of lambda, we find that Advanced LIGO observations of binaries at a distance of 100 Mpc will probe only unusually stiff equations of state, while the proposed Einstein Telescope is likely to see a clean tidal signature.