A General-relativistic Determination of the Threshold Mass to Prompt Collapse in Binary Neutron Star Mergers

TL;DR

This paper develops a universal relation between the threshold mass for prompt collapse in binary neutron star mergers and the maximum compactness, improving previous models and using GW170817 data to constrain neutron star radii.

Contribution

It introduces a novel, rigorous method to determine the threshold mass and establishes a nonlinear universal relation with maximum compactness, enhancing the accuracy over previous linear models.

Findings

Existence of a nonlinear universal relation between threshold mass and maximum compactness.

Improved accuracy of the threshold mass determination compared to previous linear models.

Set lower limits on neutron star radii using GW170817 data and the new universal relation.

Abstract

We study the lifetimes of the remnant produced by the merger of two neutron stars and revisit the determination of the threshold mass to prompt collapse, $M_{\rm th}$. Using a fully general-relativistic numerical approach and a novel method for a rigorous determination of $M_{\rm th}$, we show that a nonlinear universal relation exists between the threshold mass and the maximum compactness. For the temperature-dependent equations of state considered here, our results improve a similar linear relation found recently with methods that are less accurate but yield quantitatively similar results. Furthermore, exploiting the information from GW170817, we use the universal relation to set lower limits on the stellar radii for any mass.

Figures

Tables

Table 1. Table 1: Properties of the maximum-mass models for the EOSs considered here: BHB Λ ​ Φ Λ Φ \Lambda\Phi   (Banik et al., 2014 ) , (HS-) DD2  (Typel et al., 2010 ) , SFHo  (Steiner et al., 2013 ) , (HS-)TM1  (Hempel et al., 2012 ) , and Togashi + +   (Togashi et al., 2016 ) . Reported are: the maximum mass of a nonrotating star M TOV subscript 𝑀 TOV M_{{}_{\rm TOV}} , the corresponding radius R TOV subscript 𝑅 TOV R_{{}_{\rm TOV}} , the compactness C TOV subscript 𝐶 TOV C_{{}_{\rm TOV}} , and the free-fall-timescale τ TOV subscript 𝜏 TOV \tau_{{}_{\rm TOV}} . Also shown are the threshold masses M th subscript 𝑀 th M_{\rm th} and the corresponding errors.

EOS	$M_{{}_{\mathrm{TOV}}}$	$R_{{}_{\mathrm{TOV}}}$	$C{}_{\mathrm{TOV}}$	${\tau }_{{}_{\mathrm{TOV}}}$	$M_{\mathrm{th}}$	$\Delta M_{\mathrm{th}}$
	$\left[M_{\odot }\right]$	$\left[\mathrm{km}\right]$		$\left[\mu s\right]$	$\left[M_{{}_{\mathrm{TOV}}}\right]$	$\left[M_{{}_{\mathrm{TOV}}}\right]$
BHB$\Lambda \Phi$	2.10	11.64	0.26	83.31	1.503	0.005
DD2	2.42	11.94	0.30	80.60	1.364	0.020
SFHo	2.06	10.34	0.29	70.44	1.391	0.016
TM1	2.22	12.6	0.26	91.70	1.520	0.015
Togashi$+$	2.23	10.17	0.32	66.12	1.298	0.000
\botrule

Equations

$${\tau }_{\mathrm{ff}}(M,R):=\frac{\pi }{2}\sqrt{\frac{R^3}{2M}}.$$

$$t_{\mathrm{merg}}:$$

$$t_{{}_{\mathrm{BH}}}:$$

$$\frac{M_{\mathrm{th}}}{M_{{}_{\mathrm{TOV}}}}\approx 1.415,$$

$$\frac{M_{\mathrm{th}}}{M_{{}_{\mathrm{TOV}}}}=a-\frac{b}{1-c{\mathcal{C}}_{{}_{\mathrm{TOV}}}},$$

$$R_x=-0.88M^2+2.66M+8.91.$$

Full text

1 Introduction

The recent detection of gravitational waves from the merger of neutron-star binaries (Abbott et al., 2017) has heralded the new era of multi-messenger gravitational-wave astronomy. These observations offer new insight into the most extreme objects in universe, namely, neutron stars, and allow us to probe and constrain the properties of nuclear matter (Annala et al., 2018; Paschalidis et al., 2017; Bauswein et al., 2017; Radice et al., 2018; Most et al., 2018; Burgio et al., 2018; Montana et al., 2018).

When two neutron stars merge, they will produce an object that either collapses promptly to a black hole, or does not (Baiotti et al., 2008). In the latter case, the remnant may be a metastable object, e.g., a hypermassive neutron star (HMNS), eventually collapsing to a black hole on a secular timescale, or survive for much longer times, either as a rotating or a nonrotating star (see, e.g.,  Baiotti & Rezzolla, 2017, for a review). In the case of the first detection of merging neutron stars, GW170817 (Abbott et al., 2017), the precise fate of the merger remnant is presently unknown, although the formation of a black hole naturally matches the simultaneous observation of a short gamma-ray burst (Eichler et al., 1989; Rezzolla et al., 2011), and has been the working hypothesis to set new limits on the maximum mass of neutron stars (Margalit & Metzger, 2017; Shibata et al., 2017; Rezzolla et al., 2018; Ruiz et al., 2018).

Determining the time of collapse of the merger remnant is particularly challenging as there are a number of physical processes that either determine or undermine the stability of merger remnant. These include: the ejection of matter (Rosswog et al., 1999; Kyutoku et al., 2014; Radice et al., 2016; Lehner et al., 2016; Dietrich & Ujevic, 2017; Bovard et al., 2017), the angular-momentum transfer via magnetic fields (Kiuchi et al., 2015; Siegel et al., 2013; Kawamura et al., 2016), the evolution of the degree of differential rotation (Kastaun et al., 2016; Hanauske et al., 2017), and possible viscous effects mediated either by neutrinos or magnetic fields (Duez et al., 2004; Shibata & Kiuchi, 2017; Radice, 2017; Alford et al., 2018).

The determination of the critical (threshold) mass to a prompt collapse, $M_{\rm th}$, is much simpler, although it still poses numerical and conceptual challenges. Bauswein et al. (2013) have been the first to explore this problem by employing a smooth-particle approximation for the hydrodynamics and a conformally flat approximation to general relativity. In this way, they were able to find a linear universal relationship between $M_{\rm th}$ and the compactness of the maximum-mass model, $\mathcal{C}_{{}_{\rm TOV}}:=M_{{}_{\rm TOV}}/R_{{}_{\rm TOV}}$, where $M_{{}_{\rm TOV}}$ and $R_{{}_{\rm TOV}}$ are respectively the mass and radius of the maximum-mass nonrotating star. Here, we improve on this result by using a fully general-relativistic approach, a wider range of compactnesses, and a rigorous definition of the threshold mass. As a result, we find a nonlinear relation between $M_{\rm th}$ and $\mathcal{C}_{{}_{\rm TOV}}$, which offers a better match to the numerical-relativity results. Furthermore, exploiting the information from GW170817, we use the new relation to set more stringent lower bounds on the radii neutron stars.

2 Methods and setup

To describe the evolution of the merging system, we solve the coupled Einstein-hydrodynamics system (Rezzolla & Zanotti, 2013) using the Einstein Toolkit (Löffler et al., 2012). In particular, we evolve the spacetime with the McLachlan code (Brown et al., 2009), with the same gauges as in Hanauske et al. (2017). On the other hand, we evolve the matter with the high-order relativistic-hydrodynamics code WhiskyTHC (Radice & Rezzolla, 2012; Radice et al., 2014). The numerical grid uses the fixed-mesh refinement driver Carpet (Schnetter et al., 2004), with a total of six refinement levels having a highest resolution of $215\,{\rm m}$ covering the two stars and a total extent of $700\,{\rm km}$. For one equation of state (EOS), we also considered different resolutions of $215,287$ and $573\,{\rm m}$, obtaining threshold masses within a variance of $\Delta M_{\rm th}\lesssim 0.005\,M_{\odot}$.

A number of zero-temperature (“cold”) EOSs are available for numerical simulations and rely on nuclear-physics calculations following a variety of theoretical approaches. While these EOSs are suitable to describe the inspiral phase, they become obviously inadequate after the merger, when the temperatures reach values of several tens of MeV. To counter this, it is customary to model the post-merger dynamics by modifying these EOS and adding a “thermal” contribution via an ideal-fluid EOS (Rezzolla & Zanotti, 2013) so as to account for the shock heating (Janka et al., 1993). This approach, while not self-consistent, is rather robust and such “hybrid EOSs” have been employed extensively in the literature (Baiotti & Rezzolla, 2017).

However, since our goal here is that of determining as accurately as possible the threshold mass to prompt gravitational collapse, it is essential that the description of the thermal effects in the matter is as realistic and self-consistent as possible. In turn, this forces us to consider EOSs that have a physically consistent dependence on temperature. Unfortunately, the number of EOSs that can be employed for this scope and that do not violate some basic nuclear-physics requirement (as it is the case for the widely employed LS220 EOS, Tews et al., 2017), is much more restricted. Here, we have employed all of the five “hot” EOSs that can be used with confidence to determine the threshold mass. The corresponding properties, when expressed in terms of the masses and radii of the maximum-mass nonrotating configuration (hereafter TOV), are reported in Table 1. Similarly, Fig. 1, shows the masses and radii of the TOV equilibrium solutions, both stable (solid lines) and unstable (dashed lines). Solid squares mark the maximum-mass solutions, while open circles the models used.

Using these EOSs, we have modelled the initial data under the assumption of irrotational quasi-circular equilibrium (Gourgoulhon et al., 2001) and computed it via the LORENE library for a total of 15 equal-mass (i.e., $q=1$) binaries. The initial separation is of $45\,{\rm km}$, so that the binaries perform around five orbits before the merger. We note that since the threshold mass for equal-mass binaries is always larger than for unequal-mass binaries, i.e., $M_{\rm th}(q=1)>M_{\rm th}(q<1)$, the use of equal-mass binaries is not a restriction but optimises the search for $M_{\rm th}$ (see also Bauswein et al., 2017, for a discussion).

The measurement of the threshold mass inevitably imposes a clear definition, but also a straightforward procedure to extract this information from the numerical-relativity simulations. Quite generically, one expects that the lifetime of the HMNS, or collapse time $t_{\rm coll}$, will decrease as the mass of the binary is increased, so that the threshold mass will represent the shortest possible lifetime. This definition, is however inconvenient since different EOSs will yield different “shortest lifetimes” and comparing different EOSs may introduce clear biases. We resolve this problem by building our analysis around two important logical steps. First, we consider the collapse time as a dimensionless quantity by expressing it in terms of the free-fall time $\tau_{\rm ff}$ that, for an object of mass $M$ and radius $R$, is given by (Rezzolla & Zanotti, 2013)

[TABLE]

Since $R\sim 1/M$ for stable models, the smallest free-fall time will be achieved for the maximum-mass model, so that the shortest free-fall timescale is $\tau_{{}_{\rm TOV}}:=\tau_{\rm ff}(M_{{}_{\rm TOV}},R_{{}_{\rm TOV}})$. Second, we define the threshold mass $M_{\rm th}$ as the one for which the merger remnant will collapse over such a timescale, i.e., $M/M_{{}_{\rm TOV}}\to M_{\rm th}/M_{{}_{\rm TOV}}$ for $t_{\rm coll}/\tau_{{}_{\rm TOV}}\to 1$.

We next discuss the procedure followed for the measurement of the collapse time, $t_{\rm coll}$ as given by the difference between the (coordinate) time of merger, $t_{\rm merg}$, and that of collapse to a black hole, $t_{{}_{\rm BH}}$, i.e., $t_{\rm coll}:=t_{{}_{\rm BH}}-t_{\rm merg}$. Both of these times can be measured in a number of different ways. The first one involves the emission of gravitational waves, with $t_{\rm merg}$ being given by the time of the first maximum of the gravitational-wave strain amplitude $h$, or of the Weyl scalar $\Psi_{4}$. The time of collapse, on the other hand, can be estimated as the time when the ringdown signal starts. While overall robust, measurements involving gravitational waves are prone to errors as $t_{\rm merg}$ depends sensitively on the phase evolution of the binary in its most nonlinear stage111As the threshold mass is approached, the merger can take place a fraction of a radian earlier/later than expected, biasing the measurement.. Similarly, determining $t_{{}_{\rm BH}}$ is complicated by the fact that the beginning of the ringdown is somewhat arbitrary and ringdown itself can be modified by the matter infalling onto the black hole.

To counter these difficulties, $t_{\rm merg}$ could be measured via the proper separation between the two stellar cores, marking the merger as the time when such a separation is below a fraction of the initial diameters of the two stars. Similarly, $t_{{}_{\rm BH}}$ could also be measured in terms of the first appearance of an apparent horizon. Although these measurements show an overall consistency, they suffer in from the fact that the proper separation is very sensitive to the properties of the EOSs, and that the first appearance of the apparent horizon is ultimately set by the frequency at which it is searched during the simulations.

A third approach for measuring $t_{\rm coll}$ and $t_{{}_{\rm BH}}$ involves instead the monitoring of the minimum of the lapse function $\alpha$, which we evolve employing a singularity-avoiding “$1+\log$” slicing (Alcubierre, 2008). This quantity has been shown to be a very good proxy for the tracking and appearance of an apparent horizon (Alcubierre, 2008) and has the advantage of being extremely robust. As a result, we can mark the two times respectively when

[TABLE]

These values for $\min(\alpha)$ are less arbitrary than they may appear as $\alpha_{\rm merg}=0.35$ and $\alpha_{{}_{\rm BH}}=0.2$ systematically represent the first minimum and zero of second derivative of the function $\min(\alpha(t))$, respectively; also, at $\alpha_{{}_{\rm BH}}=0.2$ the apparent horizon is also normally first found. More importantly, the results are invariant under a change of the values in (2).

In summary, all methods to compute $t_{\rm coll}$ provide results that are consistent across different choices. However, when it comes to robustness and simplicity of implementation, the monitoring of the minimum of the lapse function represents the optimal choice and is the one employed in the results that will be presented next.

3 Results

Figure 2 reports the collapse times normalised to the free-fall timescale of the maximum-mass models, $t_{\rm coll}/\tau_{{}_{\rm TOV}}$, for all of the EOSs considered; these times are then shown as a function of the initial half-mass of the binary normalized to the maximum-mass, $M/M{{}_{\rm TOV}}$. The adoption of such set of dimensionless quantities has the goal of revealing a behaviour of the threshold mass that is universal, i.e., only weakly dependent on the EOS  (see, e.g.,  Breu & Rezzolla, 2016; Weih et al., 2018, for some examples).

Filled circles of different colours in Fig. 2 report the numerical data for the various EOSs and mark the values of the initial masses in the binaries. Note that as $t_{\rm{coll}}$ decreases, even small differences in the initial masses can lead to rather large differences in the survival time. Also note that only two values are reported for the Togashi$+$ EOS and these differ by only $3.7\%$ in mass (i.e., $M=1.440,\,1.435\,M_{\odot}$); a binary with a slightly smaller mass (e.g., $M=1.430\,M_{\odot}$) leads to a HMNS effectively stable over the timescales investigated here. Finally, since $t_{\rm{coll}}$ should diverge for vanishingly small values of $M$, we fit the numerical data with a simple exponentially decaying function of the type $M/M_{{}_{\rm{TOV}}}=\tilde{a}\exp[-\tilde{b}(t_{\rm coll}/\tau_{{}_{\rm TOV}})^{2}]$. A rapid inspection of Fig. 2 shows that the exponential fit is rather good and that a linear approximation would overestimate the threshold mass. Indeed, near the free-fall limit, the behaviour of $M/M_{{}_{\rm{TOV}}}$ should not be linear as in this limit infinitesimal changes in $M$ are sufficient to yield a prompt collapse. Such a behaviour, frequently encountered in critical-collapse calculations, requires the function $M/M_{{}_{\rm{TOV}}}$ to have vanishing derivative for $t_{\rm coll}/\tau_{{}_{\rm TOV}}\to 1$. Clearly, our nonlinear fitting reflects this behaviour while a linear one does not.

Figure 2 also reveals that the threshold mass is roughly given by

[TABLE]

with an uncertainty of $\Delta M_{\rm th}=0.05M_{\odot}$, i.e., a relative error of $\sim 4\%$. Hence, Eq. (3) provides the lowest-order approximation between the threshold mass and the corresponding maximum mass. The existence of a relation of this type has been suggested initially by Bauswein et al. (2013), who concluded that ${M_{\rm th}}/{M_{{}_{\rm TOV}}}=k$, with $k$ a linear function of the maximum compactness $\mathcal{C}_{{}_{\rm TOV}}$. More specifically, the universal ansatz proposed by Bauswein et al. (2013) is that $k=\hat{a}\,\mathcal{C}_{{}_{\rm TOV}}+\hat{b}$, with $\hat{a}=3.38,\,\hat{b}=2.43$, independent of the EOS.

Such a linear ansatz does represent a reasonable first approximation to the data but not the most general one. In particular, if it represents a universal behaviour of compact self-gravitating objects it should be valid for all possible compactnesses and provide the expected black-hole limit, for which ${M_{\rm th}}/{M_{{}_{\rm TOV}}}\to 0$ for $\mathcal{C}_{{}_{\rm TOV}}\to 1/2$. Hence, we correct the linear approximation via a nonlinear fit of the type

[TABLE]

where $a,b,c$ are to be determined from the data. However, imposing the fulfilment of the black-hole limit removes one free parameter and sets $a=2b/(2-c)$.

Figure 3 reports with a solid blue line the fit (4) with $b=1.01,\,c=1.34$, against the numerical-relativity data shown with stars of the same colorcode as in Fig. 2 (see Table 1 for the errors in the fit). Also shown with a red-dashed line is the linear approximation of Bauswein et al. (2017), which clearly suggests larger threshold masses, most probably because the conformally flat approximation used by Bauswein et al. (2013) underestimates the strong-curvature behaviour that characterises the threshold to black hole collapse. At the same time, the difference with the linear approximation of is not enormous and is of $8\%$ at most for the cases considered here.

Additionally, we show in the inset of Fig. 3 a comparison of the linear and nonlinear fittings for the threshold mass in the whole range of possible compactnesses, i.e., $\mathcal{C}_{{}_{\rm TOV}}\in[0,1/2]$. The light-green area refers to neutron stars, with the lower limit $\mathcal{C}_{{}_{\rm TOV}}\gtrsim 0.2$ deduced from a large sample of EOSs (Most et al., 2018) and the upper limit $\mathcal{C}_{{}_{\rm TOV}}\lesssim 0.35$ set by causality; instead, the upper limit of the dark-green area $\mathcal{C}_{{}_{\rm TOV}}<4/9$ is set by the Buchdahl limit for compact stars (Rezzolla & Zanotti, 2013).

Besides providing an improvement over the linear approximation of Bauswein et al. (2013), the nonlinear expression (4) can now be used to provide more stringent (i.e., larger) lower limits on the radii of possible stellar models in the light of the recent detection of GW170817 (Abbott et al., 2017). In particular, following Bauswein et al. (2017), we use (4) to plot the threshold mass for different potential choices of $R_{{}_{\rm TOV}}$; this is shown in the left panel of Fig. 4 with black solid lines and for $R_{{}_{\rm TOV}}=10,11,12\,{\rm km}$. Also reported in Fig. 4 with a grey-shaded area is the limit set by causality and that requires $M_{{}_{\rm TOV}}/R_{{}_{\rm TOV}}\lesssim 0.354$ (Koranda et al., 1997). As noted by Bauswein et al. (2017), given the merger of a neutron-star binary with total mass $M_{\rm tot}$, it is possible to set a lower limit on $M_{\rm th}$. This is shown in the left panel Fig. 4, where we report with a horizontal blue-dashed line the total gravitational mass estimated for GW170817, $M_{\rm tot}=2.74^{+0.04}_{-0.01}\,M_{\odot}$ (Abbott et al., 2017).

The corresponding uncertainty band (blue-shaded area) gives a lower constraint on $M_{\rm th}$, since GW170817 did not lead to a prompt collapse. The blue band thus constraints the red-shaded area from below, yielding a lower limit for the radius of the maximum-mass star, $R_{{}_{\rm TOV}}\geq 9.74^{+0.14}_{-0.04}\,{\rm km}$ (red solid line); this is to be contrasted with the value deduced by Bauswein et al. (2017), i.e., $R_{{}_{\rm TOV}}\geq 9.26^{+0.17}_{-0.03}\,{\rm km}$, on the basis of their linear approximation. Interestingly, to obtain a more stringent constraint similar to the one derived here, Bauswein et al. (2017) required a hypothetical detection of a binary with $M_{\rm tot}\simeq 2.9\,M_{\odot}$.

The logical approach and the mathematical procedure followed so far to derive the nonlinear fit (4) for $\mathcal{C}_{{}_{\rm TOV}}$ can be repeated for the compactness of a fixed mass $M_{x}$, i.e., $\mathcal{C}_{x}:=M_{x}/R_{x}$, thus allowing us to set constraints not only on $R_{{}_{\rm TOV}}$, but on any radius $R_{x}$ within a reasonable range (the fit becomes increasingly bad for large masses since the EOSs do not have solutions well above two solar masses). This is shown in the right panel of Fig. 4, where the values of $M_{x}$ and $R_{x}$ are indicated with blue crosses, while in black is the quadratic fit

[TABLE]

The importance of expression (5) is that it offers a handy expression for the lower limit of stellar models as deduced from GW170817. A similar procedure has been followed by Bauswein et al. (2017), but only for a fixed mass of $1.6\,M_{\odot}$, from which it was deduced that $R_{1.6}\geq 10.30\,{\rm km}$; this result should be contrasted with the value derived from (5), which is instead $R_{1.6}\geq 10.90\,{\rm km}$. Similarly, for a reference star of $1.4\,M_{\odot}$ we obtain $R_{1.4}\geq 10.92\,{\rm km}$, which is worryingly close to the estimate by Bauswein et al. (2017) for $1.6\,M_{\odot}$222Bauswein et al. (2017) provide an estimate only for $R_{1.6}$.. On the other hand, our $R_{1.4}$ estimate is in good agreement with that of Most et al. (2018), who have explored a large number of possible EOSs and built a set of one billion stellar models from which they deduced that $12.00<R_{1.4}/{\rm km}<13.45$.

4 Conclusions

The detection of the merger of a binary system of neutron stars has concretely initiated the process of extracting information on the most extreme state of matter from gravitational-wave signals. An important phase of this process lies in understanding the post-merger behaviour of the binary and the stability of the remnant. Indeed, knowing whether the system promptly forms a black hole is critical to understanding and interpreting the electromagnetic signals that may be observed.

Using a fully general-relativistic approach and a novel method for the determination of the threshold mass, we have carried out simulations making use of all of the realistic EOSs available to describe this process. In this way, we have found a nonlinear universal relation for the threshold mass as a function of the maximum compactness and which is potentially valid for all compactnesses. At least for the temperature-dependent EOSs considered here, this universal relation improves the linear relation found recently with methods that are less accurate, but that also yield quantitatively similar results. Furthermore, exploiting the detection of GW170817, we have used the universal relation to set lower limits on the stellar radii for any any mass.

These results can be improved in at least two ways. First, as new hot EOSs becomes available for numerical simulations it will be possible to extend the analysis carried here, reducing its uncertainty. Second, as new detections from binary neutron-star mergers will be revealed, the masses of these systems and their electromagnetic counterparts will be used to set ever more precise lower bounds on the radii of neutron stars.

Acknowledgements

Support comes from: the ERC synergy grant “BlackHoleCam: Imaging the Event Horizon of Black Holes” (Grant No. 610058), “PHAROS”, COST Action CA16214; the LOEWE-Program in HIC for FAIR; the European Union’s Horizon 2020 Research and Innovation Programme (Grant 671698) (call FETHPC-1-2014, project ExaHyPE). The simulations were performed on the clusters SuperMUC (LRZ, Garching), LOEWE (CSC, Frankfurt), and HazelHen (HLRS, Stuttgart). The EOSs employed can be found on stellarcollapse.org.