Exploring properties of highdensity matter through remnants of neutronstar mergers
Abstract
Remnants of neutronstar mergers are essentially massive, hot, differentially rotating neutron stars, which are initially strongly oscillating. As such they represent a unique probe for highdensity matter because the oscillations are detectable via gravitationalwave measurements and are strongly dependent on the equation of state. The impact of the equation of state for instance is apparent in the frequency of the dominant oscillation mode of the remnant. For a fixed total binary mass a tight relation between the dominant postmerger oscillation frequency and the radii of nonrotating neutron stars exists. Inferring observationally the dominant postmerger frequency thus determines neutron star radii with high accuracy of the order of a few hundred meters. By considering symmetric and asymmetric binaries of the same chirp mass, we show that the knowledge of the binary mass ratio is not critical for this kind of radius measurements. We perform simulations which show that initial intrinsic neutron star rotation is unlikely to affect this method of constraining the highdensity equation of state. We also summarize different possibilities about how the postmerger gravitationalwave emission can be employed to deduce the maximum mass of nonrotating neutron stars. We clarify the nature of the three most prominent features of the postmerger gravitationalwave spectrum and argue that the merger remnant can be considered to be a single, isolated, selfgravitating object that can be described by concepts of asteroseismology. We sketch how the consideration of the strength of secondary gravitationalwave peaks leads to a classification scheme of the gravitationalwave emission and postmerger dynamics. The understanding of the different mechanisms shaping the gravitationalwave signal yields a physically motivated analytic model of the gravitationalwave emission, which may form the basis for templatebased gravitationalwave data analysis. We explore the observational consequences of a scenario of two families of compact stars including hadronic and quark matter stars. We find that this scenario leaves a distinctive imprint on the postmerger gravitationalwave signal. In particular, a strong discontinuity in the dominant postmerger frequency as function of the total mass will be a strong indication for two families of compact stars.
pacs:
26.60.KpEquations of state of neutronstar matter and 97.60.JdNeutron stars and 04.30.DbWave generation and sources and 95.85.SzGravitational radiation, magnetic fields, and other observations and 04.25.dkNumerical studies of other relativistic binaries and 21.65.QrQuark matter1 Introduction
By means of intense experimental efforts gravitational waves (GWs) are expected to be detected within the next years. To date, their existence is only indirectly proven by the observations of neutron star (NS) binaries (119); (114); (127), in particular the HulseTaylor binary pulsar PSR B1913+16 (71). The orbital period and orbital separation of this binary decrease precisely as determined by General Relativity, which predicts that the GW emission continuously reduces the system’s angular momentum and energy. As the orbital frequency increases with time the GW emission becomes stronger, and thus the decay of the orbit proceeds increasingly faster. For PSR B1913+16 the orbital period will decrease from currently 7.75 hours to a few milliseconds in the next years. The corresponding GW emission reflects the dynamics of the orbital motion and thus increases in frequency and amplitude resulting in a chirplike signal, which is determined by the binary masses (see e.g. the first milliseconds in Fig. 12 or Figs. 5 and 6 in (18)). Only during the last seconds of this socalled inspiral (named after the shape of the stellar trajectories) the signal will enter the sensitivity band of the existing and upcoming GW detectors between roughly 10 and 10000 Hz (3); (11); (120). Importantly, the total population of NS binaries within a distance of a few 100 Mpc to our Galaxy (i.e. the range that will be covered by GW instruments) is expected to be sufficiently high such that per year of the order of 40 binaries in the last phase of their inspiral will enter the sensitivity band of GW detectors (2).
The continuously decreasing orbital separation of NS binaries inevitably leads to the merging of the binary components. On a dynamical time scale of the order of a millisecond the two stars form a single, massive, differentially rotating object (see e.g. Fig. 9 and Fig. 2 in (22)). Mostly because of its rotation the merger remnant can be supported against the gravitational collapse even if the total binary mass significantly exceeds the maximum mass of nonrotating NSs. Only for very high total binary masses a prompt gravitational collapse occurs. The stellar object forming during a coalescence of two NSs is very interesting from the perspective of highdensity matter physics
Because of its high mass, the postmerger remnant is particularly interesting since the central densities are relatively high compared to those of many observed NSs in binary systems, which have masses in the range of 1.2 to 1.5 (see e.g. (81)). This implies that the GW emission of the postmerger phase probes a density regime that is hardly accessible by other observational methods. In particular, attempts to infer NS properties via finitesize effects in the GW signal during the inspiral phase preceeding the merger are restricted to the mass range of the individual components of the binary (see e.g. (67); (39); (40); (54); (105); (42); (126); (5); (80); (34)). In this sense the postmerger oscillations represent a complementary approach by providing insights to NS properties at higher masses and higher densities (e.g. (23); (24)). Moreover, the collapse behavior of NS merger remnants may reveal the properties of the very high density regime, which is decisive for determining the threshold of the gravitational collapse (15). The consideration of the collapse behavior may be particularly rewarding for establishing limits on the maximum mass of nonrotating NSs, which is challenging given the paucity of NS systems with masses close to the collapse threshold, while the mass range of merger remnants partially overlaps with the mass range where the gravitational collapse is expected to take place for various EoSs. Therefore, an upper limit on the maximum mass may be established through NS mergers, which may be hard to derive through other types of observations.
The basic impact of the EoS on the merger dynamics and GW signal can be described as follows. Stiff EoSs lead to NSs with large radii. Such stars can be deformed more easily under the influence of an external tidal field. Consequently, finite size effects during the inspiral set in at a larger orbital separation, i.e. a lower orbital frequency, and the stars finally merge at a relatively low orbital frequency. In contrast, soft EoSs yield more compact NSs, which in comparison to larger NSs behave more like point particles during the late inspiral phase. Stars described by soft EoSs are harder to deform and during the inspiral they reach higher orbital velocities before they merge. This also implies a higher linear velocity before merging and results in a merger with a higher impact velocity. The stiffness also affects the dynamics of the postmerger phase and, in particular, the frequencies of the excited oscillation modes. A stiff EoS leads to a relatively large merger remnant, whose quadrupolar fluid oscillation frequency is relatively low, since it scales approximately with the mean density. In the case of a soft EoS, the merger remnant is more compact and thus oscillates at higher frequencies. (For such EoSs the higher impact velocity during merging additionally leads to a stronger excitation of the quasiradial oscillation mode of the remnant.)
In this paper we summarize the current status of possibilities to deduce NS properties and EoS constraints from GW observations of the postmerger phase of NS mergers. We focus mostly on the most prominent oscillation mode in the GW spectrum of the postmerger phase, the fundamental quadrupolar fluid mode, since this is the feature which is most likely to be detected. We prepend a discussion of mass measurements through GW chirp signals of NS mergers, which is essential to NS radius measurements by the observation of GWs from the postmerger phase as well as for constraints on the maximum mass of nonrotating NSs. We supplement a review of previous findings by new results, which for instance clarify the role of intrinsic NS rotation of the binary components. Also, we provide further evidence that the merger remnant can be considered to be a single, selfgravitating object and can thus be investigated by concepts of asteroseismology (116). By considering certain oscillation modes and secondary features of the postmerger GW spectrum, we arrive at a unified classification scheme of the postmerger dynamics and GW emission. We also present an analytic model of the GW emission, which may represent a step towards the construction of GW templates that can be employed in GW data analysis. Along the lines of the current Topical Issue of “Exotic Matter in Neutron Stars” we also discuss a scenario of the existence of two families of compact stars (46) and identify certain observational features which can test these ideas. This discussion represents an example of how the postmerger phase can be employed to test highdensity matter properties at densities which are not accessible through NSs of masses below .
The results presented in this paper rely mostly on the calculations discussed in (17); (18); (15); (23); (22), where further information can be found. Additional new findings laid out here are obtained within the same physical and numerical model, for which further details are provided in (96); (95); (19); (21); (18). We only consider mergers from quasicircular orbits (for a discussion of tidal capture events see e.g. (82)). This paper is written in a modular way such that the individual sections are mostly selfcontained and the reader can skip certain topics. Introductory remarks on mass measurements via GWs are presented in Sect. 2. NS radius determinations via the postmerger phase are discussed in Sect. 3. We address the collapse behavior of NS merger remnants and outline possibilities to measure the maximum mass of nonrotating NSs in Sect. 4. The origin of the major and of secondary features of the GW spectrum are explained in Sect. 5 as well as the dependencies of the secondary GW peaks, which lead to a classification scheme of the postmerger dynamics and GW emission. In Sect. 6 we provide details on our analytic model of the postmerger GW emission. The impact of intrinsic NS rotation is considered in Sect. 7. The possible observational implications of a scenario of two distinct families of compact stars is explored in Sect. 8. We summarize and conclude in Sect. 9.
2 Mass measurements
The strongest part of the GW signal of a NS merger originates from the preceding inspiral phase, where the system continuously losses angular momentum and energy by the emission of GWs. The losses lead to inspiralling trajectories of the binary components, while the GW signal is determined by the orbital motion of the NSs. The inspiral accelerates, which results in a chirplike GW signal with an increasing amplitude and an increasing frequency. Finitesize effects become important only for the last orbits before merging and may yield information on the NS EoS, specifically on the tidal deformability, which for fixed mass approximately scales with the stellar radius (see e.g. (67); (39); (40); (54); (105); (42); (126); (5); (80); (34)).
The key parameter determining the GW signal during the inspiral phase is the socalled chirp mass apart from other parameters like the distance to the source, the inclination of the binary, the orientation of the instruments and so on. The chirp mass is given by
(1) 
with the NS masses and . Being the crucial quantity to determine the signal, the chirp mass will be measured with very high accuracy since it is one of the parameters used to construct the template bank for matched filtering searches that will be employed to actually detect GWs from binaries (55); (38); (73); (10); (123); (124); (59); (1); (107); (125); (53). In contrast, the mass ratio can only be measured with higher accuracy for nearby events (55); (38); (73); (10); (123); (124); (59); (1); (107); (125); (53).
It can be seen from Eq. (1) that a measurement of the chirp mass alone does not suffice to determine the individual binary masses, which can only be achieved by determining the mass ratio with sufficient precision. However, for NS binaries the chirp mass alone provides already a very good estimate of the total binary mass . Figure 1 shows the total binary mass as a function of the mass ratio for sequences of constant chirp mass (solid lines). For a constant chirp mass the total binary mass depends only weakly on the mass ratio. The dashed lines define the binary systems, where the less massive component has a gravitational mass of 1.1 and 1.2 , respectively. Under the reasonable assumption that NSs cannot be less massive than 1.2 (see discussion below), a measured chirp mass restricts the possible binary parameters to the right of the dashed curve for . This implies that for instance a measured chirp mass of 1.1752 restricts the total binary mass of this event to the range (without any information about the mass ratio). From observed NS binaries (see e.g. listing in (81)) and theoretical population synthesis studies (e.g. (45)) it is expected that most binaries have a chirp mass close to 1.1752 (corresponding to a total binary mass of roughly 2.7 ). This means that for the most likely merger event the total binary mass can be determined very accurately from the chirp mass measurement only. Even for very massive binaries (with a chirp mass of 1.3058 ) the total binary mass can be inferred with a precision of .
To support the arguments above we remark that recent corecollapse simulations for a large number of progenitor stars find a minium NS rest mass of about 1.3 (51). A NS rest mass of 1.3 corresponds to a gravitational mass of roughly 1.2 , somewhat dependent on the EoS. These findings justify to assume that NS binaries are unlikely to host NSs less massive than 1.2 . Also other formation channels like an accretioninduced collapse of a white dwarf are unlikely to form less massive NSs. See also the discussion of the minimum NS mass in (81). Finally, we note that a white dwarfNS binary might have a chirp mass in the range which is typical of NSNS binaries, and a more extreme mass ratio if the white dwarf is less massive than 1.2 . The merger of a white dwarfNS binary will, however, lead to significantly altered inspiral dynamics (103) and thus a different GW signal prior to merging such that a white dwarfNS merger can be easily distinguished from a NS binary merger.
The considerations above play a role for the following discussions, because methods to infer EoS properties from NS mergers rely on the ability to measure at least the total binary mass from the inspiral GW signal. Information on the mass ratio and thus the individual binary component masses is not critical, but improves the constraints on NS properties. As outlined above, the determination of the total binary mass via the chirp mass alone represents the absolute minimum of what is achievable with existing and upcoming GW detectors. However, it is important to stress that for distances of order 50 Mpc, which allow to deduce EoS properties, the individual binary masses are expected to be recovered with a precision of a few per cent (55); (38); (73); (10); (123); (124); (59); (1); (107); (125); (53). For instance, in (107) the individual NS masses have been recovered within about 10 per cent at a distance of 100 Mpc. Assuming that the error scales linearly with the distance, an accuracy of a few per cent could be reached in determining the individual masses of the merging NSs at a distance of a Mpc. Therefore, we will work in the following under the condition that the individual masses can be measured sufficiently accurately for our purposes and that, in particular, the total binary mass can be determined very well.
3 Radius measurements
The outcome of a NS merger essentially depends on the EoS of highdensity matter and the total binary mass. The binary mass ratio has a modest impact, while the initial NS rotation and magnetic fields have only a small influence.) Simulations show that for binary masses of about 2.7 the merging leads to the formation of a massive, hot, differentially rotating NS remnant for most EoSs, even for EoS models which yield a maximum mass of nonrotating NSs significantly below 2.7 , e.g. (130); (109); (14); (87); (111); (113); (112); (95); (94); (7); (84); (12); (76); (116); (58); (70); (110); (17); (18); (102); (108); (15); (69); (29); (23); (117); (75); (104); (118); (74); (22); (28); (101); (43); (44); (41); (56). The rapid rotation stabilizes the remnant against gravitational collapse. According to pulsar observations in NS binaries most systems are expected to have a total mass of roughly 2.7 (see e.g. the compilation of binary masses in (81)), which is also supported by theoretical studies of the binary population, e.g. (45). Given this mass range, the most likely outcome of a merger event should be expected to be the formation of a NS remnant (see (15) for the EoS dependence of collapse behavior).
A typical GW spectrum of a NS merger resulting in the formation of a NS remnant is shown in Fig. 2. The spectrum is computed from a simulation of a 1.351.35 merger with the DD2 EoS (63); (122). The dominant GW oscillation frequency of the postmerger phase is clearly visible as a pronounced peak in the kHz range, which is produced by the dominant quadrupolar remnant oscillation. Apart from the main peak one recognizes several additional peaks, whose nature will be discussed in Section 5. By computing the spectrum of the GW signal of the postmerger phase only, one can show that the different features are related to the postmerger phase (see thin lines in Fig. 2).
The dominant oscillation frequency depends in a specific way on the highdensity EoS (17); (18). This is understandable since the EoS affects the size of the remnant, which in turn determines its oscillation frequency (see Fig. 13 in (18)). The structure of the remnant is also influenced by its angular momentum. The available angular momentum, however, is given by the dynamics of the late inspiral/merging phase, which is fully determined by the stellar structure of the inspiralling stars and thus also depends on the EoS in a particular way. The strong EoS dependence of the peak frequency can be expressed as follows. An EoS which is used in a given simulation, can be conveniently characterized by the radii of nonrotating NSs, which are uniquely determined by this EoS through the stellar structure equations (TolmanOppenheimerVolkoff equations (121); (98)). Specifically, for a set of calculations
with
a fixed total binary mass but different EoSs
The empirical relation between and is very tight, which implies that a measurement of the peak frequency can be used to determine the unknown radius of a nonrotating NS with a fixed mass by simply inverting the empirical relation (17); (18). Thus, a future detection of the GW postmerger phase and extraction of the peak frequency (see (35); (36)) will yield strong constraints on the highdensity EoSs. In (17); (18) the largest deviation of the empirical data from a fit is only a few hundred meters. The accuracy of a radius determination by the postmerger GW signal is mostly affected by two sources of error. One error is the uncertainty of the measurement of the peak frequency. Apart from this, one should take into account deviations between the data and the fit to the data describing the empirical relation. A measurement of the peak frequency (of the true EoS) does not reveal in which way the measured frequency slightly deviates from the empirical relation. Hence, one conservatively has to assume that the true data point may deviate as much as the largest deviation found in the large sample of candidate EoSs.
The peak frequency has been shown to be measurable with very high precision by a coherent burst search analysis (35). In this study waveforms from numerical models were superimposed with the recorded data stream of previous GW detector science runs, which simulates the noise of the future instruments. The model waveforms were injected at random times and the noise was rescaled to the anticipated sensitivity of the secondgeneration GW detectors Advanced LIGO and Virgo. The existing GW data analysis pipeline was able to recover the injected signal and to determine the peak frequency with an accuracy of Hz, which is smaller than the spread in the empirical relation between and the NS radius. This implies that the radii of the inspiralling stars can be determined with a precision of a few hundred meters.
These considerations show that the larger contribution to the error of a radius measurement originates from the scatter in the empirical relation between and . In this context, the following observation is important. One has the freedom to choose any fiducial NS mass different from for characterizing a given EoS by the TOV radius . Empirically, it turns out that using a fiducial NS mass somewhat larger than leads to tighter relations between and . This is exemplified in Fig. 3. For (circles in Fig. 3) the maximum deviation between the data and a fit amounts to only meters if is chosen. This implies that the measurement of the dominant postmerger frequency for determines the radius of a nonrotating 1.6 NS with an accuracy of better than 200 meters.
It is natural that a fiducial mass of is somewhat more appropriate than for characterizing the postmerger oscillations of 1.351.35 mergers ( ). The maximum densities in the massive, rotating merger remnant are higher than in the initial NSs and they are comparable to the central densities of nonrotating, static NSs with a mass of roughly 1.6 (see e.g. Fig. 15 in (18)). For this reason, nonrotating NSs with better represent the density regime encountered in the merger remnant and thus provide a better description of the EoS.
Figure 3 also shows that similar empirical relations hold for other binary masses. The different symbols display the peak frequencies from calculations with total binary masses of 2.4 (plus symbol), 2.7 (circles), and 3.0 (crosses) for symmetric mergers. The solid lines are leastsquare fits to the data. One recognizes that also for other binary masses there are only small deviations from the empirical relations on the order of only a few 100 meters. Hence, these relations can be used for determining NS radii with similar precision after the total binary mass was deduced from the inspiral signal (see discussion in Sect. 2). The total binary mass determines which of the relations has to be employed to convert the measured peak frequency to a NS radius. Note that all data are plotted as a function of and hence the fiducial mass is not optimized to yield the tightest relations for total binary masses of 2.4 or 3.0 . Here we assume that the determination of the binary masses does not contribute a significant error, but that the binary masses are known with sufficient precision as argued in Sect. 2. We also point out that very similar relations can be constructed for unequalmass binaries with a fixed total mass and mass ratio.
We here do not further discuss asymmetric binaries in great detail, but refer for instance to (18). There it is shown that the peak frequencies of asymmetric binaries deviate only somewhat from the ones of the symmetric binaries of the same total mass. This can be also seen in Fig. 4, which shows the peak frequency as function of for symmetric and asymmetric mergers. The symmetric and asymmetric binary systems are chosen such that they have the same chirp mass of (see Sect. 2), which will be measured very accurately from the GW inspiral signal. This chirp mass corresponds to a symmetric setup with two stars of 1.34 (black symbols
The precise form of the fits of the data in Fig. 3 and the deviations from the fits of course depend slightly on the chosen sample of candidate EoSs, which should preferentially include all EoSs that are acceptable based on current knowledge. Moreover, uncertainties in the numerical or physical model may also lead to slight differences in the empirical relations, which, however, will diminish in the future by more detailed simulations. Also, it is possible to choose different functional forms which may to some small extent affect the accuracy of the empirical relations. It is worth noting that in contrast to (118), we do not see difficulties in fitting a function to their data for the dominant oscillation frequency. Using the listed frequencies
In this context we note that the peak frequencies from our numerical simulations agree well with those of (110); (69); (117); (118); (56). Hence, it would be surprising if the data of (117); (118) did not yield a close relation for a fixed binary mass while our data show a good correlation and the frequencies essentially agree.
In (18) it was argued that relations between and NS radii are not unexpected because the dominant emission is generated by the fundamental quadrupolar oscillation mode (116), whose frequency is known to scale tightly with for nonrotating NSs (see (9)). This scaling for nonrotating stars suggests that it may be possible to rescale the frequency by dividing by to describe the data of different total binary masses by a single relation. While this is indeed possible, empirically we find that results in an even tighter relation. Figure 5 shows the rescaled peak frequency as a function of for total binary masses of 2.4 , 2.7 , and 3.0 . The rescaled frequencies can be described by the quadradic leastsquares fit
(2) 
with in kHz, the binary mass in and in km. The data in Fig. 5 deviate by at most 500 meters from the fit. The reason why yields a better universality than (which results in deviations of up to km) may be that the radius is also mass dependent. A perfect scaling behavior with the mass cannot be expected considering the findings in (23), which show that different EoSs can lead to different dependencies on the total binary mass (see Fig. 1 in (23)). While such a universal relation of the rescaled is theoretically interesting, from a practical point of view it is preferable to consider different sets of massdependent relations as in Fig. 3 (or possible interpolations between them) because they yield tighter relations, and the total binary mass will always be known from the GW inspiral signal. Concerning massindependent relations for we also note that a relation between the dominant postmerger frequency and the tidal coupling constant, pointed out in (28), is interesting because it connects the inspiral and the postmerger signals. Such a relation is understandable because essentially both quantities are known to depend on the EoS and to scale with the NS radius. As in the case introduced above (Fig. 5) the practical use of such relations depends on their tightness and a spread of Hz may be large compared to a typical width of the main postmerger peak.
4 Gravitational collapse and estimates of the maximum mass
Apart from determining NS radii as discussed in the previous section, the postmerger phase also offers the possibility to constrain the maximum mass of nonrotating NSs. The maximum mass is a key parameter regarding the properties of highdensity matter because its value is determined by the very highdensity regime that is usually not encountered in most observed NSs, which have a lower mass. The maximum mass is thus particularly interesting for probing exotic phases of matter, which potentially occur at higher densities. Here we outline three different methods to determine or at least constrain the maximum mass of nonrotating NSs. An important quantity in the following discussion is the threshold binary mass to prompt black hole (BH) collapse. If one considers binary mergers with different total binary masses for a given EoS, one finds that for a binary mass above a certain threshold the merging leads directly to the formation of a BH and no (transiently stable) NS remnant is formed because the merged object is unstable against gravitational collapse. Note that the gravitational collapse in a merger remnant sets in at densities below the maximum density of nonrotating NSs.
The connection between this threshold mass and the maximum mass of nonrotating NSs has the following background. The maximum NS mass is the threshold for BH formation of static, nonrotating stars, while the threshold mass represents the threshold for BH formation of differentially rotating, hot NSs. Thus, a quantitive relation between those two quantities is expected; and a determination of is possible since the threshold mass can be determined from the binary masses of different merger events considering their respective outcome.
Method 1: An analysis of a large number of simulations with different EoSs and different total binary masses has revealed that the threshold binary mass depends in a particular way on the EoS (15). This dependence can be described by TOV properties, which are uniquely defined by the EoS. With very good accuracy the threshold mass is given by
(3) 
with the gravitational constant and the speed of light . Similarly, it can also be expressed as
(4) 
i.e. using the maximum mass and the compactness of the maximummass configuration. The coefficients in both equations are obtained by fits to the ratio , where is given by the results from simulations (15). Equations (3) and (4) reproduce the numerical results with a precision better than 0.1 (on average the deviations between the numerical and the estimated threshold mass are 0.024 for Eq. (3) and 0.031 for Eq. (4)). Bear in mind that the numerical value of can only be determined with a finite precision because apart from uncertainties of the numerical and physical model only a limited set of models with different has been analyzed. (74) tested the collapse behavior for two EoS models, which were also employed in (15), and found compatible results.
Using Eq. (3) Fig. 6 visualizes how depends on the maximum mass of nonrotating NSs for different chosen values of . If the radius has been determined observationally, for instance through the detection of the peak frequency of a 1.351.35 merger (see Sect. 3), then the threshold mass becomes an unambiguous function of . The threshold mass may be observationally determined by confirming or excluding the presence of postmerger GW emission of a NS remnant for several merger events with different total binary masses. (Note that the prompt formation of a BH leads to very weak postmerger GW emission, which can be distinguished from the formation of a NS remnant for nearby events and sufficient sensitivity.) The inversion of the relation between and for a fixed can then be employed for an estimate of the maximum mass of nonrotating NSs.
We remark that a distinction between the prompt collapse and the formation of a NS remnant may also be possible by the observation of electromagnetic counterparts of mergers. In particular, the thermal emission of the ejecta, which are heated by radioactive decays, may be observable (83); (79); (91); (93). The ejecta mass depends on the EoS, binary mass and mass ratio (see e.g. (16)). The prompt collapse to a BH leads to smaller amounts of unbound matter compared to mergers which form a NS remnant. This should lead to a noticeable difference in the light curves of these two different outcomes and may allow an observational identification of prompt collapse events. The detectability of electromagnetic counterparts to GW events has been addressed for instance in (90); (93); (89) and depends on (currently) uncertain details like the observer angle, the available instruments, the search strategy, the observing conditions, the ejecta mass, the outflow velocity, and the neutronrichness of the ejecta determining the opacity. The last three source properties are affected by the highdensity EoS and the binary parameters.
Method 2: We point out that already a single GW event may provide a relatively precise estimate of the threshold mass and a constraint on the maximum mass of nonrotating NSs. Figure 7 displays the threshold mass as a function of the peak frequency of equalmass mergers with for different EoSs. We use here the threshold masses estimated via Eq. (3), which are very close to the numerical values. The blue symbols show the results for additional simulations with representative EoSs from (62), which are derived within a chiral effective field theory at lower densities (up to roughly nuclear saturation density) and supplemented by extrapolations of the EoS at higher densities with physically motivated parameter variations (minimum NS mass of 2 and causality). For more information on this set of EoSs see (18). Since these EoSs describe only cold NS matter, the EoSs are supplemented with an approximate treatment of thermal effects for the merger simulations (see e.g. (19) for details and an assessment of this approximation). In these calculations the “idealgas” index , which regulates the thermal pressure support, was chosen to be 5/3. See e.g. (19) and (37) for a motivation of this value. One should keep in mind that the exact choice of introduces an ambiguity and thus an uncertainty of the numerical values of certain quantities extracted from simulations. This is the reason why we did not include these models in the discussion in Sect. 3. Note that from the original set of six EoSs in (18), two models lead to the prompt formation of a BH for 1.51.5 binaries. This is fully compatible with the theoretically estimated threshold mass from the TOV properties of these EoSs via Eqs. (3) and (4).
Figure 7 shows that the determination of of a 1.51.5 merger yields the threshold mass with a precision of about . Note that a total binary mass of 3.0 is only slightly larger than the masses of observed binaries for which precise mass measurements exist. Thus, a detection of a merger with may not be improbable.
At first glance it may seem tempting to use the peak frequency of a 1.51.5 merger to fix (Fig. 7), then to determine as discussed in Sect. 3 (Fig. 3), and then to estimate via Eq. (3). This, however, represents an attempt to determine two independent quantities ( and ) by only one observable. In fact, a thorough error analysis for only one peak frequency measurement reveals that depending on the EoS the maximum mass can be determined only with a large error bar. (Notice the flat slope for small radii in Fig. 6.) These considerations are summarized in Fig. 8, which shows directly as a function of the peak frequency of 1.51.5 mergers. Remarkably, immediately yields an upper limit for , which for soft EoSs implies even a good estimate of (recall the lower bound on by the measurement of NSs with about 2 ). In Fig. 8, the dashed line, which is given by
(5) 
with the frequency in kHz and the mass in , indicates an exclusion region and thus defines an upper limit on , which can be derived from a single detection of a merger with . We stress that some of the EoSs within our sample (especially some of the models by (62), blue symbols) have a speed of sound equal to the speed of light at higher densities. These models are thus maximally stiff, which is favorable for yielding high maximum masses. We thus conjecture that the dashed line represents a true upper limit and is not an artifact of the chosen sample of EoSs. We also refer to Fig. 16 in (18), which shows a plot similar to Fig. 8 but for simulations of 1.351.35 mergers.
Fig. 8 shows that a single merger event can only constrain the maximum mass. However, we emphasize that the estimate of the threshold mass, e.g. through Fig. 7, represents a viable goal on its own. The threshold mass is the crucial quantity to judge the outcome of a detected merger event. For instance, at larger distances only the chirp mass can be measured (which determines the total mass with some accuracy, see Sect. 2). Having an estimate of available from a previous nearby event, the total mass of a distant merger may be crucial information for interpreting a possibly detected electromagnetic counterpart, for instance from thermal emission by the ejecta, which is powered by the radioactive decay of nucleosynthesis products (83); (79); (91); (89). It is known that the prompt collapse leads to a reduced ejecta mass and thus to dimmer electromagnetic counterparts, which reach their peak emission on a shorter time scale e.g. (16). Also, for the coincident detection of a short gammaray burst (99); (50); (27) and a GW signal of a merger the information on the collapse behavior will provide valuable information and help to understand the conditions leading to this gammaray burst.
Method 3: While a single event can only yield an upper limit on the maximum mass , we have recently shown that several events with a measurement of the dominant postmerger GW frequency may suffice to determine the maximum mass (23). Specifically, mergers with different binary masses (but in the most likely range of binary masses) can be employed for an estimate. The procedure has the following background. Measuring the peak frequencies for two different total binary masses allows to estimate how the dominant GW frequency depends on . This can be used to extrapolate and thus to estimate the behavior at higher binary masses, which is particularly sensitive to the maximum mass and the radius of the maximummass configuration. In essence, the absolute values of and its slope are very sensitive to and and thus yield an estimate of the maximum mass of nonrotating NSs with an accuarcy of roughly . For details we refer to (23) and note that a slightly different approach is described in (24), which results in a similar precision. We also point out that this extrapolation procedure determines the radius of the maximummass configuration and the maximum density of nonrotating NSs, which are quantities that are highly characteristic for the highdensity regime of the EoS.
The advantage of this approach lies in the fact that it relies only on detections of mergers with binary masses in the most likely range, i.e. roughly between 2.4 and 3.0 (see e.g. (81)). This contrasts the procedure to estimate by a direct measurement of the threshold mass as sketched for Method 1. Estimating the threshold mass directly by determining the outcome of mergers with binary masses above and below requires several merger detections at higher binary masses, which are possibly very unlikely.
5 Characterization of GW peak frequencies
5.1 Origin of main and secondary GW peaks
We start with a brief summary of our current understanding of the mechanisms which produce certain peaks in the GW spectrum (22). For the interpretation of the main and secondary peaks it is important to realize that the merger remnant can be considered as an isolated selfgravitating object (116). This view is supported by Fig. 9 showing the lapse function in the equatorial plane for a 1.351.35 merger with the DD2 EoS (63); (122). The lapse function may be seen as the relativistic analog of the gravitational potential. In the upper left panel the system is still composed of two cores corresponding to the initial NSs. Shortly later, the two cores have merged into a single potential well. The time steps of the snapshots in Fig. 9 are the same as the ones in Fig. 2 of (22), which shows the evolution of the restmass density in the equatorial plane for the same model. It is remarkable that in the density evolution, a doublecore structure persists for several milliseconds after merging, whereas the lapse function exhibits already a single core with partially strong deformations. (In the density plots corresponding to the upper right and lower panels of Fig. 9 a clear doublecore structure is still visible.) These findings imply that the doublecore structure visible in the density field for many milliseconds should not be interpreted as two gravitationally interacting objects, but rather as local overdensities moving in a single gravitational potential. Hence, the double cores should be viewed as tracers of the dynamics of a single, isolated object, rather than two independent dynamical objects.
The time evolution of the remnant structure is also summarized in Fig. 10. It shows the time evolution of the central lapse function, which may be interpreted as a measure for the compactness of the stellar object. The two initial NSs first touch at about ms. The vertical line indicates the time when the two cores in the lapse function (upper left panel in Fig. 9) merge into a single core, which occurs already during the first compression phase, i.e. right during the final plunge. The subsequent oscillations in the central lapse function indicate quasiradial oscillations (bounces and compressions of the remnant). The figure implies that matter accumulates into a single, gravitational trough already at a very early time in the remnant’s evolution. Consequently, right from its formation the remnant’s evolution can be described by the (nonlinear) dynamics of a single, selfgravitating object. For instance, the dominant osillation can be associated with the fundamental quadrupolar fluid oscillation mode. This can be seen by adding a velocity perturbation to the remnant at late times, when it reaches a quasiequilibrium and the GW emission essentially has diminished. We add an instantaneous perturbation of
(6) 
to the component of the coordinate velocity as a function the polar angle , the azimuthal angle and the radial coordinate . (Geometrical units are adopted.) The radial coordinate of the surface, which only depends on , is defined by the coordinate at which the density drops below some threshold. We evolve the remnant with the added velocity perturbation for several milliseconds and extract the GW signal. The spectrum of the GW signal of the orginal simulation is compared to the spectrum of the perturbed model in Fig. 11. It is evident that the perturbation excites an oscillation mode with the same frequency as the dominant remnant oscillation, which strongly suggests that is the frequency of the fundamental mode
Also other (secondary) peaks in the GW spectrum can be explained by oscillation modes of the remnant. For instance, the peak at 1.5 kHz in Fig. 2 is produced by a quasilinear interaction between the fundamental quadrupolar mode and the quasiradial oscillation of the remnant (116). This can be shown by determining the frequency of the quasiradial mode, which itself does not occur in the GW spectrum. However, it is very pronounced in the time evolution of the central lapse function (Fig. 10) and in other characteristic properties of the remnant, such as the maximum density and the size of the remnant. Again, by adding a suitable velocity perturbation to the latetime remnant, one can predominantly excite the quasiradial mode and can extract its frequency from the time evolution of the central lapse function. This is visible in Fig. 10, where the green dashed curve shows the evolution of the radially perturbed model. The determination of the frequencies of the quadrupolar mode () and of the radial mode () reveals that a secondary peak is expected to occur at , which indeed is the case. The frequency coincidence confirms the nature of this secondary peak as being a coupling of two modes, which is why we refer to this feature as the peak. Another corresponding combination frequency can be recognized in Fig. 2 at approximately . The peak at is observationally less interesting because of its weakness and the smaller sensitivity of GW detectors at higher frequencies, but it substantiates the importance of mode couplings in NS merger remnants.
Finally, there is one more secondary peak visible in the GW spectrum displayed in Fig. 2. Recently, we provided evidence that this feature is generated by a spiral deformation which is created during merging. This deformation cannot follow the faster rotation of the inner remnant. The spiral deformation forms antipodal bulges, which orbit around the central part of the remnant for several milliseconds. Being a strong, nonaxisymmetric, orbiting deformation, the antipodal bulges generate a GW signal at a frequency which is twice the orbital frequency of the bulges. A deeper analysis of the simulation data confirms this origin of the GW peak at 2 kHz for teh particular model (e.g. by extracting the orbital motion of the antipodal bulges and estimating their mass, by comparing the presence of the bulges and the presence of the secondary peak in the GW spectrum, and by computing GW spectra for the inner and outer parts of the remnant separately to estimate the contribution of the different remnant components to the different GW features). See (22) for more details. Clearly, this socalled feature cannot be explained within a perturbative approach. One can show by an analytic model that a peak with an appropriate strength in the GW spectrum can be produced by point particles of a few 0.1 , which orbit for only a few milliseconds with an orbital separation of roughly the diameter of the inner remnant (see also Sect. 6).
5.2 Classification of postmerger GW emission
Considering models of NS mergers with varied binary masses and with different EoSs, one realizes that the three most prominent features in the GW spectrum can be explained by the three mechanisms detailed above: the fundamental quadrupolar mode, the coupling of the quadrupolar and the quasiradial mode, and the orbital motion of antipodal bulges. The peak of the fundamental mode is present in all models and it is always the strongest feature. The presence and the strength of the different secondary features, however, are sensitively affected by the total binary mass. Depending on the total binary mass relative to the threshold mass (see Sect. 4) the different secondary peaks are more or less pronounced. This leads to a classification scheme that relies on the presence and strength of the different secondary features. One can identify three different types of NS mergers and corresponding GW spectra.
For relatively high , i.e. close to but below , the quasiradial mode is strongly excited during merging, and consequently the feature is the most prominent secondary feature (Type I). In comparison, the peak is weaker and may even be hardly visible in the spectrum. For moderate binary masses both secondary features are clearly visible and can be clearly distinguished in the spectrum (Type II). The secondary peaks have a roughly comparable strength. In a third type of GW spectra, the feature is absent (or hardly visible) and the GW peak by the spiral deformation is the most prominent secondary feature. This Type III occurs for relatively low total binary masses, that means for much below . (An additional case is the prompt collapse to a BH as discussed in Sect. 4 for .)
Since the threshold mass depends on the EoS, the notation of a “relatively high” or “relatively low” binary mass is EoSdependent as well. As a consequence, for all three types of mergers can occur depending on the EoS (see Fig. 5 in (22)). For EoSs which lead to compact NSs, the initial stars merge with a higher impact velocity (see Fig. 3 in (16)), which is why the quasiradial mode is strongly excited resulting in a pronounced peak (Type I). For stiff EoSs, i.e. less compact stars, the merging proceeds more gently with a lower impact velocity, which suppresses the strong excitation of the quasiradial mode. Such models favor the formation of pronounced spiral deformations, and thus the peak becomes particularly pronounced. The same reasoning explains the occurence of the different merger types for a fixed EoS dependending on the total binary mass.
We note that the classification of different merger simulations in (22) is done based on the GW spectrum of the full signal (inspial and postmerger phase). Considering the postmerger spectrum only (i.e. windowing
Finally, we mention that windowing the GW data has an impact in particular on the precise frequency of the peak. Choosing the postmerger spectrum only or even excluding the very early postmerger phase, the peak shifts to higher frequencies and its strength decreases. The latter is understandable because one excludes a part of the signal, when the mechanism generating the feature is operating. The frequency shift may not be unexpected, given that the peak is produced by a highly dynamical feature. For instance, one may expect that the antipodal bulges orbit at slightly smaller radii and thus lead to slightly higher orbital frequencies.
5.3 Frequency dependence of secondary peaks (and EoS constraints)
The dependence of the dominant oscillation frequency on the EoS was extensively discussed in Sect. 3, and especially the potential for an accurate determination of NS radii (or, equivalently, NS compactnesses at fixed masses) was pointed out. The frequencies of the secondary peaks and follow a behavior very similar to that of , i.e. for a fixed total binary mass the frequencies are higher for EoSs which yield smaller NS radii. Relating the frequencies with the radii of nonrotating NSs (as in Sect. 3) one finds similarly tight correlations for and somewhat less tight for (see (22), see also (36) for a plot relating the secondary frequencies directly to the dominant frequency ). The frequencies follow the order .
In principle, secondary frequencies may be employed for EoS constraints (117); (118); (22), however, a few issues are worth being mentioned. The secondary peaks are weaker than the main peak even if one takes into account the better sensitivity of GW detectors at lower frequencies. Hence, secondary frequencies will be more difficult to measure. In addition, the secondary peaks are broader in comparison to and often do not stand out clearly from the background. This will further impede an accurate determination of the secondary frequencies. Also, a detection of a single secondary peak will require additional information, e.g. from a measurement of , to safely associate the secondary peak with either or . In comparison, a detection of the dominant peak frequency is more likely because of its strength, and, as discussed in Sect. 3, a single detection of is sufficient to yield an accurate determination of NS radii.
Given the similarity to we do not further discuss the dependencies of the secondary peak frequencies here, but we refer to (22); (36) for more details. Regarding the implications for EoS constraints we, however, note that our simulations in (22) do not confirm the existence of a massindependent universal relation for the strongest secondary peak as claimed in (117); (118)
6 Analytic model for postmerger GW emission
The understanding of the most prominent GW emission mechanisms as detailed in Sect. 5 motivate us to set up an analytic model for the postmerger GW signal. This model may form the basis for GW templates to be used in matched filtering GW searches (see (35) for the potential of templatebased searches in comparison to morphologyindependent burst search algorithms). We report the model here by specifying the xy component of the reduced quadrupole moment of the GW source, which can easily be interpreted as the cross polarization of the GW signal along the polar direction. The other components of the quadrupole moment can be deduced with the same set of parameters, from which the complete GW signal can be derived (i.e. both polarizations in all emission directions). For instance, the plus polarization can be obtained by adding an appropriate phase shift. In (22) we identify three main mechanisms that produce the dominant postmerger GW emission. As also detailed in Sect. 5 this includes the dominant oscillation by the fundamental quadrupole mode, the coupling of this mode with the quasiradial mode of the remnant and the transient emission from a spiral deformation. All three mechanisms can be modelled by individual sine functions with given initial amplitude, initial phases and exponential damping behavior. Hence, our analytic model reads
(7) 
for with a starting time . The particular advantage of our model (for instance in comparison to the fits to GW signals described by (69)) is that the parameters are motivated by the underlying physical mechanisms. Hence, it is not difficult to choose appropriate values for these parameters if one wants to reproduce a given GW signal, e.g. from a numerical simulation. The frequencies , and can be chosen as found in the spectra. The amplitudes and the damping time scales can be estimated from the time evolution of the GW signal and adjusted such that the correspoding peaks in the GW spectrum coincide with numerically obtained spectra. Generally, the amplitude and the damping timescale of the fundamental mode (“peak”) should be the largest. However, for the model discussed above (Fig. 2), may have a comparable strength but a shorter damping time scale than the contribution from the fundamental mode. As argued in (22), the antipodal bulges of the spiral deformation disappear after a few orbits, which suggests a damping time scale of the order of a few milliseconds.
The physical background of the different parameters is also advantageous because the individual parameters can be constrained to certain ranges. For instance, the damping time scales and are always smaller than . The restricted range of certain parameters is an important property of our model because it signficantly reduces the parameter space and thus the computational costs in future applications as templates in matched filtering searches. The range and number of parameters (for instance the amplitudes and the frequencies) may be constrained if information about the total binary mass (or at least the chirp mass) is available from the GW inspiral signal. In this context, we point out that for a fixed total binary mass and a fixed mass ratio the three frequencies are highly correlated (22); (36). For symmetric binaries with , and are tight functions of , which can be well approximated by
(8) 
and
(9) 
with the frequencies in kHz. Since these relations are accurate with deviations of typically only a few 10 Hz, and can be essentially eliminated from our model, which reduces the dimensionality of the parameter space. If this precision is not sufficient for template searches one can at least define a very narrow range of frequencies around the estimates of Eqs. (8) and (9). Further correlations and constraints of the different parameters will be explored in future work. For an application of this model in a templatebased GW search, which we leave to future work, possible refinements of the damping behavior or small frequency drifts may be considered as well. Also a certain rise time of the signal may be implemented or the model may be directly connected to an inspiral signal (see (28)). For strongly aymmetric binaries modifications may be necessary.
We note that our model naturally leads to a time varying instantaneous GW frequency if the signal is interpreted as being produced by a single instantaneous frequency (see (69); (118); (74)). The interpretation of a single instantaneous frequency, however, cannot be supported by the underlying physical mechanisms producing the GW signal (see Sect. 5 and (22); (36)). Instead, there is evidence that several mechanisms with different frequencies contribute simultaneously to the signal. (See also the timefrequency map in (36), which clearly shows distinct frequencies being simultaneously present and relatively stable in time.) This feature also distinguishes our analytic model from the model in (118), which describes only a single instantaneous orbital frequency which contributes to the GW signal and thus neglects the presence of the orbiting antipodal bulges generating the distinct peak.
To exemplify the potential of our analytic model we show in Fig. 12 the postmerger GW signal of a simulation with the DD2 EoS (black) together with the signal of our analytic model (dashed blue curve, Eq. (7)) computed for a chosen set of parameters. As described above, the frequencies are taken from the GW spectrum and appropriate values are chosen for the remaining parameters. We stress that for this example the parameters were not determined by a fitting procedure but simply by a visual trialanderror comparison of the numerical waveform and the model signal. The very good match between the waveforms in Fig. 12 is obtained without extensive finetuning of the parameters, but only by a crude, physically motivated choice of the parameters. An even better match can be achieved by an appropriate fitting procedure or more elaborate tuning of the parameters of the model (see (69)). This is shown in Fig. 13, where we employ a
publicly
available
algorithm
Finally, we note that for the feature in the model shown in Figs. 12 and 14 we have chosen an amplitude of and a damping time scale of ms. An amplitude of this magnitude corresponds to two point particles of about 0.2 to 0.3 orbiting at roughly the surface of the remnant. This further substantiates the finding in (22) that the peak in the GW spectrum is generated by antipodal bulges which form during merging and which orbit around the inner remnant for only a few milliseconds.
7 Impact of rotation
Most NS merger simulations have focussed on binaries with irrotational velocity profile, i.e. binaries with NSs that do not rotate intrinsically. This is justified by the conclusion that the viscosity is not sufficient to enforce corotation of the binary components by tidal interactions during the inspiral phase (30); (78). In addition, the intrinsic rotation of NSs in binaries is expected to be relatively slow since the stars cannot be spun up by accretion because of a missing donor star. In fact, the fastest known pulsar in a NS binary system has a rotation period of only 22 ms (see e.g. list in (85)), which is slow compared to the orbital period of the binary prior to merging, which is of the order of roughly 2 ms. An intrinsic rotation period of 22 ms is also slow in the sense that the stellar structure is practically unaffected by the rotation, which only becomes important for spin periods below ms. Finally, one should bear in mind that most NSs in binaries rotate even more slowly than 22 ms, and that in the particular case of a NS spin of 22 ms the rotation will be further reduced by pulsar spin down due to magnetic dipole radiation until the binary components merge on a time scale of roughly 100 Myrs.
Here we are mostly interested in the impact of intrinsic NS rotation on the GW signal, in particular on the dominant peak frequency of the GW spectrum, which can be employed for accurate NS radius constraints. With regard to the discussion in Sect. 3, the most crucial question is whether intrinsic NS rotation can substantially alter for expected rotation rates compared to models without intrinsic rotation. (The relations like the ones in Fig. 3 could still be used to determine NS radii but with some loss of precision if the intrinsic rotation cannot be measured from the GW inspiral phase (but see e.g. (59) for a discussion of spin measurements).) Comparing irrotational, corotating and counterrotating spins (with respect to the orbital motion) only very little influence on was observed by (94). More recent studies in full general relativity found a somewhat larger effect of the order of 100 Hz for intrinsic rotation rates roughly comparable to 22 ms (29). A comparable influence was reported by (74) but for roughly three times faster rotation rate, suggesting a somewhat weaker effect. However, these results and the observed impact of rotation may not be representative because the chosen binary systems are close to the threshold to prompt BH formation, and most of these simulated NS merger remnants experience the delayed collapse to a BH a few milliseconds after merging. For such binary setups the postmerger object undergoes strong changes in the stellar structure, and therefore the GW signal may be more sensitively affected by small changes, e.g. by a small difference in the angular momentum of the remnant. Moreover, the calculations of (29) employed an idealgas EoS, which may not represent a proper description of NS matter at all densities. (43) considered unequalmass binaries with fast intrinsic rotation ( exceeding the rotation rate of the fastest known pulsar in binary systems) and found a shift of by only a few 10 Hz (in addition one setup with a misalignment between orbital and intrinsic rotation was investigated).
In order to assess the impact of intrinsic rotation for a more representative setup we performed calculations of 1.351.35 binaries with the DD2 EoS and varied the initial intrinsic rotation rate using spin periods of 19 ms, 7 ms, 4 ms, and 2.2 ms. These simulations were then compared to the results from a calculation with irrotational velocity profile, which represents the usual setup used for deriving the results in Sects. 3 to 5. Our simulations were performed with the code described in (96); (95); (19); (18), which imposes the conformal flatness condition on the spatial metric for solving the Einstein equations (72); (128). For simulations of irrotational binaries this approximation showed a very good quantitative agreement with fully general relativistic calculations (comparisons were reported in e.g. (17); (18); (117); (56)).
Figure 15 demonstrates that the impact of initial NS rotation on the dominant postmerger GW frequency is practically negligible. Note that all tested models with NS spins rotate faster than the fastest known pulsar in a binary. For the model with the very fast rotation period of about 2.2 ms the peak frequency shifts by roughly 50 Hz relative to the model without initial NS spin, whereas the calculation with an initial NS spin period of 19 ms is practically identical to the computation without NS spin. One observes that models with faster initial spin lead to slightly lower peak frequencies and somewhat lower peak heights. This behavior is understandable from the fact that the remnants in such models have slightly more angular momentum. Hence, they are less compact and oscillate at lower frequencies. Also, the higher angular momentum slightly damps the excitation of the fmode oscillation because of the centrifugal barrier, which is why the height of the main postmerger peak is slightly reduced. This can also be seen by considering the minimum of the central lapse function in the first compression phase (see Fig. 10 for the model without intrinsic NS rotation), which reveals that mergers with initial spin are not as much compressed during the plunge and during the first compression phase as the corresponding model without spin.
Finally, the insensitivity of the peak frequency with respect to variations of the initial NS spin is understandable by considering the angular momentum of the remnant. Compared to the model without NS spin, the remnant’s angular momentum is only slightly higher in models with spin because most of the angular momentum originates from the orbital motion of the inspiralling NSs. The contribution from the NS spins is only at the level of a few per cent (only for the model with a spin period of 2.2 ms the additional angular momentum slightly exceeds 10 per cent of the orbital angular momentum). Hence, differences in the remnant’s structure are marginal, which explains why the peak frequencies are hardly affected by initial rotation even if the NSs spin relatively fast in comparison to measured spins in NS binaries. Therefore, the consideration of intrinsic NS spins can be safely neglected for the purpose of measuring the NS radii via the dominant postmerger oscillation frequency (Sect. 3).
The different models displayed in Fig. 15 show that the secondary peaks may be affected by very fast intrinsic NS rotation. While the model with a rotation period of 19 ms does not exhibit important differences compared to the calculation with an irrotational velocity profile, the binaries with faster spinning NSs lead to a stronger feature. This makes sense considering our explanations about the origin of this peak (Sect. 5). Clearly, initial NS rotation favors the formation of antipodal bulges and thus leads to higher peaks. Future work should also consider cases with different spin orientations and magnitudes of the individual binary components.
8 Twofamilies scenario
Recently, it has been proposed by (46) that two different families of compact stars, ordinary hadronic NSs and absolutely stable strange stars (SSs), might coexist under the assumption that the strange matter hypothesis is correct (31); (129) (see e.g. (26); (32); (25) for similar ideas and (77) for arguments against absolutely stable strange matter). Specifically, the authors argued for a soft hadronic EoS and a stiff quark matter EoS. In this scenario, lowmass compact stars are NSs and compact stars with masses above a certain threshold are absolutely stable SSs, i.e. objects entirely composed of quark matter except for a dynamically unimportant crust of nucleonic matter with densities below the neutron drip density. A conversion of NSs to SSs may be triggered by different processes, e.g. mass accretion onto a NS. The massradius diagram in Fig. 16 illustrates such a jump from the NS branch to the SS branch. During the conversion (e.g. (97); (86); (64); (100); (46); (48)) the baryon number of the compact star is conserved, and the change in the binding energy leads to a different gravitational mass. While clearly speculative, the scenario of two families of compact stars represents a possible solution to the hyperon puzzle, i.e. the scenario is compatible with a hadronic EoS which is strongly softened by the occurrence of hyperons, and a maximum mass of compact stars above . Here we do not perform actual simulations which treat the conversion process and the detailed postmerger evolution, but based on theoretical arguments we discuss the observational signature of such a scenario without further commenting on the likelihood of the scenario or considering possible counterarguments. For instance, it still remains to be explored to which extent the existence of SSs could lead to a cosmic flux of strangelets, which would prevent the existence of NSs because seed strangelets would convert a NS into a SS (6); (88); (33). For SS EoSs which are compatible with the lower observational bound of the maximum mass of nonrotating compact stars, we expect that significant amounts of strange quark matter become gravitationally unbound during a merger of SSs (20). Whether there exists a cosmic flux of strangelets, however, depends also on the fragmentation and survival probability of the unbound strange quark matter (see e. g. (68) for a discussion of a possible reconversion of quark matter nuggets to ordinary nucleonic matter). We focus on the GW emission of mergers of compact stars within the described scenario.
We consider only one exemplary model of the twofamilies scenario, i.e. one pair of nucleonic and quark matter EoSs (see Fig. 16). (Specifically, for the hadronic EoS we use the relativistic mean field model of (47), which includes hyperons and delta resonances and which is based on the model of (115) with densitydependent coupling constants. We consider the model where the ratios between the coupling constants of the mesons to the isobars and the coupling constants of the mesons to the nucleons are chosen to be . The EoS of absolutely stable strange quark matter adopts the model of (57) with the parameter fixed to 3.5. The EoS tables were provided by G. Pagliara.) Clearly, variations to this model are possible, but the case discussed here is sufficient for the sake of outlining the special observational signatures of such a scenario and the qualitative behavior. We expect no qualitative differences for other realizations of the twofamily scenario. We consider only the GW emission of the postmerger phase and assume that several GW observations of compactstar mergers are made with different total binary masses (see also (21) for an earlier discussion of the GW signature of SS mergers). We restrict the discussion to symmetric mergers.
For the particular model discussed here the transition of a NS to a SS occurs for NSs more massive than (gravitational mass), which corresponds to a central energy density of (see Fig. 16). This means that for total binary masses below 2.9 two NSs will merge, since the mass of each binary component is below 1.45 . This should be considered as the generic case given that many observed NS binaries have a total mass below . Prior to merging the maximum density in the initial stars starts to decrease because of the deformations in the very last inspiral phase. During merging the densities then increase again (see e.g. Fig. 3 in (21)). The maximum density in the merger remnant shows some oscillations and typically the maximum density in the remnant is larger than the one of the initial stars. However, in particular for lowmass mergers with below , the maximum densities do not significantly exceed the one of the inspiralling stars because of the strong rotational (and thermal) effects in the remnant. Hence, for lower total binary masses the densities in the remnants remain below and are thus not sufficient to trigger the conversion of the NS matter to absolutely stable quark matter (here we neglect a possible influence of thermal effects on the conversion threshold). For larger the densities in the remnanst may exceed the central density of a nonrotating NS with , and thus the conversion of NS matter to quark matter will set in. We note that a first density maximum is reached directly in the very first compression phase during merging. During the subsequent evolution of the remnant the densities are typically lower than this first maximum. Only after many ms, the redistribution of angular momentum in the remnant may lead to an increase of the maximum density above the first maximum. Hence, the conversion to absolutely stable strange quark matter is triggered right after merging or it occurs only on a longer time scale, after the remnant oscillations and thus the GW emission have ceased. In summary, this results in the following picture. Binaries below a certain binary mass threshold lead to a merger remnant made of NS matter. For total binary masses above the remnant encounters conditions that trigger the immediate conversion of NS matter to quark matter right during merging. Here we assume that the conversion takes place on a very short time scale that is shorter or at least comparable to the dynamical time scale of the remnant (64); (48). Hence, a quark matter remnant forms essentially immediately, which has the same total rest mass (but a somewhat reduced gravitational mass compared to the corresponding NS merger remnant).
By performing a simulation for a 1.21.2 binary with the hadronic EoS, we find that the maximum density after merging is only somewhat lower than the transition density , which triggers the conversion of a NS to a SS in a static configuration. From this we estimate that the binary mass which triggers a conversion is only slightly higher than 2.4 , and thus we adopt in the following . We stress that thermal effects may actually trigger the conversion (during merging) already at lower densities. Hence, may be somewhat lower, which we neglect for simplicity because it would require a detailed microphysical calculation of the temperaturedependent threshold.
Based on these remarks we can now estimate the postmerger GW emission of these binary mergers without actually performing simulations, which are technically nontrivial because the burning of NS matter to quark matter needs to be treated in an appropriate way including for instance an implementation of turbulent combustion. For the model discussed here we find km for the hadronic EoS and km for the SS EoSs by computing the TOV solutions (see Fig. 16). Using these properties we can estimate the peak frequency of the postmerger phase by employing the fits discussed in Sect. 3, which relate and . From the fits we can compute the expected peak frequency for total binary masses of 2.4 , 2.7 and 3.0 for both EoSs. With these data we can then interpolate linearly to predict the GW peak frequency as a function of , where we distinguish the cases and . In the former case we employ the estimates from the hadronic EoS, while for the latter case we use the peak frequency derived from km, i.e. the quark matter EoS. In Fig. 17 one can read off which GW peak frequency is expected for a given total binary mass, which is measured by the inspiral GW signal (see Sect. 2). It is clear that the formation of a SS remnant leads to dramatic changes of the postmerger GW signal compared to the GW emission from a NS remnant, which occurs for low total binary masses. Figure 17 demonstrates that two families of compact stars result in a jump to lower peak frequencies at . Such a scenario can be clearly distinguished from only one family of compact stars, because in this case one expects a continuous increase of with (see Fig 1 in (23)). Observationally, this requires the detection of several merger events with different total binary masses to probe whether or not a jump to lower peak frequencies occurs with increasing binary mass.
The vertical dashed line in Fig. 17 displays the threshold mass for prompt BH formation for the hadronic EoS estimated via Eq. (4). Thus, at least up to this binary mass we expect the remnant to be stable independent of how fast the onset of the conversion to strange matter can stabilize the remnant.
We note that also SS EoSs roughly follow the relation, which justifies the use of the relations discussed above (see (18)
The results in Fig. 17 represent the realization of only one possible model of the twofamily scenario. However, we expect a qualitatively similar effect for other EoS models within the twofamily scenario. In this context it is worth mentioning that for the particular model discussed here (i.e. the pair of hadronic and quark EoSs) the total binary mass for formation of a SS remnant is relatively low compared to the expected mass range of NS binaries. As argued in Sect. 2, binaries with total masses below 2.4 should be considered unlikely. Hence, the probability to detect a binary with may not be high. However, we point out that there may be a second way to discern the twofamily scenario from the ordinary model of only one family of compact stars, which may succeed with only one detection of a binary in the most likely mass range of about 2.7 . GW detectors may measure finitesize effects during the last cycles of the inspiral (67); (39); (40); (54); (105); (42); (126); (5); (80); (34). Therefore, if the EoS information revealed from the inspiral phase is in stark conflict with the EoS constraint from the postmerger phase, this could be a strong indication for the twofamily scenario, too. Note that for binaries in the most likely mass range, the inspiral probes only the hadronic EoS as long as for the model discussed here. Thus, the GW signal from the inspiral alone may not be sufficient to exclude the twofamily scenario if no binary merger with is measured. Only the consideration of both merger phases may be conclusive in this respect. Finally, we note that also the possibility of SSNS mergers for asymmetric binaries with individual masses above and below should be considered. We leave this for future work.
9 Conclusions
We summarize the main ideas to infer EoS and NS properties from the postmerger GW emission of NS binaries and in particular we review the new findings of this work.

The minimum requirement for accurate constraints of NS properties from the GW emission of NS merger remnants is the measurement of the total mass of the binary via the GW inspiral signal. The knowledge of the binary mass ratio is not critical, since the remnant’s dominant oscillation frequency depends only weakly on the initial binary mass ratio. We illustrate that in the context of NS mergers the total binary mass can be estimated very well from the chirp mass, which will be measured with very high accuracy.

For a fixed total mass the dominant oscillation frequency of the NS merger remnant scales tightly with the radii of nonrotating NSs of chosen fiducial mass for different EoSs. This in turn allows to determine the NS radius from a measurement of the dominant postmerger GW frequency, . The error is given by the maximum deviation from the empirical relation between the GW frequency and the NS radius found among all EoSs. Choosing a fiducial NS mass somewhat higher than the mass of one of the merging NSs minimizes the deviations in the empirical relation and allows NS radius measurements with an accuracy of about 100 to 200 meters. Consequently, already a single detection of only the dominant peak frequency is already sufficient for a tight constraint on the NS radius and thus the highdensity EoS through the frequencyradius relation.

We investigate asymmetric and symmetric binaries of the same chirp mass, which can be measured very accurately during the inspiral phase. The frequencyradius relation for the combined data from symmetric and asymmetric binaries is very tight (with maximum deviations below 300 meters). Therefore, even under the pessimistic assumption that only the chirp mass was measured accurately and no information on the binary mass ratio is available, the frequencyradius relation can be employed for accurate NS radius measurements.

For a broad range of binary masses (in the representative range between 2.4 and 3.0 ) a single relation between the dominant postmerger GW frequency and the NS radius can be constructed by rescaling the frequency with the total binary mass. Deviations from this relation are of the order of 500 meters. Thus, for practical purposes, relations for fixed total binary masses may be more useful because they are tighter and the total binary mass will be known.

From our representative simulations we conclude that intrinsic rotation of the initial NSs has no significant impact on the dominant oscillation frequency for initial rotation rates as expected in NS binaries. Thus, intrinsic NS rotation can be safely neglected for deriving NS properties and EoS constraints from GW signals of the postmerger phase.

We present three different possibilities to constrain the maximum mass of nonrotating NSs from the collapse behavior of NS mergers or from the measurement of the dominant GW frequency. One of the key observations is that the binary threshold mass for prompt BH formation is well described by pure TOV properties, specifically by the maximum mass of nonrotating NSs and by the compactness of the maximummass TOV configuration. Thus, the threshold mass, which can be observationally determined, depends in a specific way on pure EoS properties. Also other properties of the maximummass configuration of nonrotating NSs like its radius or the maximum central density can be obtained from measuring postmerger GW emission, thus probing the very high density regime of NSs.

Our analysis provides further evidence that the dominant oscillation mode of the merger remnant is the fundamental quadrupolar fluid mode. In addition, our calculations show that the merger remnant forms a single, selfgravitating object right from the time of its formation. Although there are overdensities visible in the evolution of the density in the equatorial plane, which appear as a rotating doublecore structure, these features should not be considered as independent dynamical features, but as tracers of the nonaxisymmetries of a single star.

Besides the dominant GW emission frequency (discussed above), we identify two distinct processes which produce prominent secondary peaks in the GW spectrum. One GW peak is prodcued by the orbital motion of an outer antipodal spiral pattern that forms during merging and persists for several milliseconds. Another GW peak is generated by the coupling of the fundamental quadrupolar mode with the quasiradial mode. Depending on the exact system parameters (EoS and total binary mass) one or the other secondary feature is more prominent or, for some cases, both GW peaks can have comparable strengths. The presence and the prominence of the different features follows a clear behavior and can thus be embedded in a classification scheme of the postmerger dynamics and GW emission.

We present details of an analytic model for describing the structure of the postmerger GW signal that is based on the physical mechanisms producing the GW emission. A given postmerger GW signal can be well described by our analytic model. The model parameters are physically motivated and thus specific, physically motivated bounds on the parameters can be imposed, which is advantageous for a templatebased GW search by reducing the computational effort. We point out the existence of certain parameter correlations which further reduce the dimensionality of the problem, and thus our description may form the basis for future matchedfiltering GW data analysis.

We finally explore a more speculative scenario of two families of compact stars. Within this scenario massive compact stars are made of absolutely stable strange matter, whereas lowmass compact stars are ordinary NSs. Observationally this scenario can be revealed by a strong discontinuity in the dependence of the dominant postmerger GW frequency on the total binary mass. For the falsification or verification of this scenario the postmerger GW signal may be crucial, since the inspiral phase may only test the hadronic regime for the expected range of binary masses.
Overall, these conclusions show that the postmerger evolution is a highly interesting phase for understanding NS and EoS properties. In particular, it seems that certain characteristics are only accessible through the postmerger phase. Clearly, apart from advances in the modelling of the different aspects of this object, further studies of the data analysis capabilities of the advanced GW detectors are required. The highly rewarding prospects sketched in this work strongly motivate the construction of even more sensitive GW detectors, such as upgrades of Advanced LIGO and Advanced Virgo or the Einstein Telescope, e.g. (66); (65); (4); (92); (36).
Acknowledgements.
We thank Matthias Hempel and Giuseppe Pagliara for providing EoS tables, and James Clark, Sebastiano Bernuzzi, Giuseppe Pagliara and Alessandro Drago for helpful discussions and suggestions. A.B. is a Marie Curie IntraEuropean Fellow within the 7th European Community Framework Programme (IEF 331873). A.B. acknowledges support by the Klaus Tschira Foundation. Partial support came from âNewCompStarâ, COST Action MP1304. The computations were performed at the Max Planck Computing and Data Facility (MPCDF), the Max Planck Institute for Astrophysics, and the Cyprus Institute under the LinkSCEEM/CyTera project.Footnotes
 offprints:
 Within this paper the term “high density” refers to densities roughly above nuclear saturation density, which is the regime most relevant for the stellar structure of NSs.
 Except for some models considered in Sect. 4, the EoSs discussed in this study are temperature dependent and include electrons, positrons and photons, while neutrino contributions are neglected. With regard to the resulting stellar properties these EoSs cover a representative range, which, for example, can be seen from the range of radii in Fig. 3 and maximum masses in Fig. 8.
 We did not simulate 1.341.34 mergers, but interpolated linearly the relation, which is given by the results from 1.21.2 binaries and 1.351.35 binaries. In this binary mass range a linear interpolation is a very good approximation for symmetric systems (see Fig. 1 in (23)).
 Note that the dominant oscillation frequency is called in (118) as well as in some other references.
 For this reason it was called in (116) in the context of mode identification.
 By windowing we mean the application of a window function, e.g. a Tukey window, to the GW data, which can be used to exclude certain parts of the whole signal.
 In (117); (118) no distinction between and was made, but their secondary frequency should correspond to the strongest secondary peak, which for most cases may be the feature.
 We use the Covariance Matrix Adaptation Evolution Strategy downloaded from https://www.lri.fr/hansen/cmaes_inmatlab.html, which is described in (60).
 The SS EoSs considered in (18) describe bare SSs. While the inclusion of a nuclear crust would not have an impact on the peak frequency, which is mostly determined by the highdensity regime of the EoS, the nuclear crust would lead to somewhat larger stellar radii (a few hundred meters) and thus would fulfill the relation even better than indicated in Fig. 11 of (18).
References
 Aasi, J., Abadie, J., Abbott, B. P., et al. 2013, Phys. Rev. D, 88, 062001
 Abadie, J., et al. 2010, Classical Quantum Gravity, 27, 173001
 Acernese, F., Agathos, M., Agatsuma, K., et al. 2015, Classical and Quantum Gravity, 32, 024001
 Adhikari, R. X. 2014, Reviews of Modern Physics, 86, 121
 Agathos, M., Meidam, J., Del Pozzo, W., et al. 2015, Phys. Rev. D, 92, 023012
 Alcock, C., Farhi, E., & Olinto, A. 1986, Astrophys. J., 310, 261
 Anderson, M., Hirschmann, E. W., Lehner, L., et al. 2008, Phys. Rev. D, 77, 024006
 Andersson, N., Ferrari, V., Jones, D. I., et al. 2011, General Relativity and Gravitation, 43, 409
 Andersson, N., & Kokkotas, K. D. 1998, Mon. Not. Roy. Astron. Soc., 299, 1059
 Arun, K. G., Iyer, B. R., Sathyaprakash, B. S., & Sundararajan, P. A. 2005, Phys. Rev. D, 71, 084008
 Aso, Y., Michimura, Y., Somiya, K., et al. 2013, Phys. Rev. D, 88, 043007
 Baiotti, L., Giacomazzo, B., & Rezzolla, L. 2008, Phys. Rev. D, 78, 084033
 Baumgarte, T. W., & Shapiro, S. L. 2010, Numerical relativity: Solving Einstein’s Equations on the Computer (Cambridge University Press, Cambridge)
 Baumgarte, T. W., Shapiro, S. L., & Shibata, M. 2000, Astrophys. J. Lett., 528, L29
 Bauswein, A., Baumgarte, T. W., & Janka, H.T. 2013a, Phys. Rev. Lett., 111, 131101
 Bauswein, A., Goriely, S., & Janka, H.T. 2013b, Astrophys. J., 773, 78
 Bauswein, A., & Janka, H.T. 2012, Phys. Rev. Lett., 108, 011101
 Bauswein, A., Janka, H.T., Hebeler, K., & Schwenk, A. 2012, Phys. Rev. D, 86, 063001
 Bauswein, A., Janka, H.T., & Oechslin, R. 2010a, Phys. Rev. D, 82, 084043
 Bauswein, A., Janka, H.T., Oechslin, R., et al. 2009, Physical Review Letters, 103, 011101
 Bauswein, A., Oechslin, R., & Janka, H.T. 2010b, Phys. Rev. D, 81, 024012
 Bauswein, A., & Stergioulas, N. 2015, Phys. Rev. D, 91, 124056
 Bauswein, A., Stergioulas, N., & Janka, H.T. 2014, Phys. Rev. D, 90, 023002
 Bauswein, A., Stergioulas, N., & Janka, H. T. 2015, Phys. Part. Nucl., 46, 835
 Benić, S., Blaschke, D., AlvarezCastillo, D. E., Fischer, T., & Typel, S. 2015, Astron. Astrophys., 577, A40
 Berezhiani, Z., Bombaci, I., Drago, A., Frontera, F., & Lavagno, A. 2003, Astrophys. J., 586, 1250
 Berger, E. 2014, Annual Review of Astronomy and Astrophysics, 52, 43
 Bernuzzi, S., Dietrich, T., & Nagar, A. 2015, Phys. Rev. Lett., 115, 091101
 Bernuzzi, S., Dietrich, T., Tichy, W., & Brügmann, B. 2014, Phys. Rev. D, 89, 104021
 Bildsten, L., & Cutler, C. 1992, Astrophys. J., 400, 175
 Bodmer, A. R. 1971, Phys. Rev. D, 4, 1601
 Bombaci, I., Parenti, I., & Vidaña, I. 2004, Astrophys. J., 614, 314
 Caldwell, R. R., & Friedman, J. L. 1991, Physics Letters B, 264, 143
 Chatziioannou, K., Yagi, K., Klein, A., Cornish, N., & Yunes, N. 2015, ArXiv eprints, arXiv:1508.02062
 Clark, J., Bauswein, A., Cadonati, L., et al. 2014, Phys. Rev. D, 90, 062004
 Clark, J. A., Bauswein, A., Stergioulas, N., & Shoemaker, D. 2015, ArXiv eprints, arXiv:1509.08522
 Constantinou, C., Muccioli, B., Prakash, M., & Lattimer, J. M. 2015, Phys. Rev. C, 92, 025801
 Cutler, C., & Flanagan, É. E. 1994, Phys. Rev. D, 49, 2658
 Damour, T., & Nagar, A. 2010, Phys. Rev. D, 81, 084016
 Damour, T., Nagar, A., & Villain, L. 2012, Phys. Rev. D, 85, 123007
 De Pietri, R., Feo, A., Maione, F., & Löffler, F. 2015, ArXiv eprints, arXiv:1509.08804
 Del Pozzo, W., Li, T. G. F., Agathos, M., Van Den Broeck, C., & Vitale, S. 2013, Phys. Rev. Lett., 111, 071101
 Dietrich, T., Moldenhauer, N., JohnsonMcDaniel, N. K., et al. 2015, ArXiv eprints, arXiv:1507.07100
 Dionysopoulou, K., Alic, D., & Rezzolla, L. 2015, Phys. Rev. D, 92, 084064
 Dominik, M., Belczynski, K., Fryer, C., et al. 2012, Astrophys. J., 759, 52
 Drago, A., Lavagno, A., & Pagliara, G. 2014a, Phys. Rev. D, 89, 043014
 Drago, A., Lavagno, A., Pagliara, G., & Pigato, D. 2014b, Phys. Rev. C, 90, 065809
 Drago, A., & Pagliara, G. 2015, Phys. Rev. C, 92, 045801
 Duez, M. D. 2010, Classical Quantum Gravity, 27, 114002
 Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 340, 126
 Ertl, T., Janka, H.T., Woosley, S. E., Sukhbold, T., & Ugliano, M. 2015, ArXiv eprints, arXiv:1503.07522
 Faber, J. A., & Rasio, F. A. 2012, Living Rev. Relativity, 15, 8
 Farr, B., Berry, C. P. L., Farr, W. M., et al. 2015, ArXiv eprints, arXiv:1508.05336
 Favata, M. 2014, Phys. Rev. Lett., 112, 101101
 Finn, L. S., & Chernoff, D. F. 1993, Phys. Rev. D, 47, 2198
 Foucart, F., Haas, R., Duez, M. D., et al. 2015, ArXiv eprints, arXiv:1510.06398
 Fraga, E. S., Kurkela, A., & Vuorinen, A. 2014, Astrophys. J. Lett., 781, L25
 Giacomazzo, B., Rezzolla, L., & Baiotti, L. 2011, Phys. Rev. D, 83, 044014
 Hannam, M., Brown, D. A., Fairhurst, S., Fryer, C. L., & Harry, I. W. 2013, Astrophys. J. Lett., 766, L14
 Hansen, N. 2006, in Towards a new evolutionary computation. Advances on estimation of distribution algorithms, ed. J. Lozano, P. Larranaga, I. Inza, & E. Bengoetxea (Springer), 75–102
 Harry, G. M., & LIGO Scientific Collaboration. 2010, Classical Quantum Gravity, 27, 084006
 Hebeler, K., Lattimer, J. M., Pethick, C. J., & Schwenk, A. 2010, Phys. Rev. Lett., 105, 161102
 Hempel, M., & SchaffnerBielich, J. 2010, Nucl. Phys. A, 837, 210
 Herzog, M., & Röpke, F. K. 2011, Phys. Rev. D, 84, 083002
 Hild, S. 2012, Class.Quant.Grav., 29, 124006
 Hild, S., Chelkowski, S., Freise, A., et al. 2010, Classical Quantum Gravity, 27, 015003
 Hinderer, T., Lackey, B. D., Lang, R. N., & Read, J. S. 2010, Phys. Rev. D, 81, 123016
 Horvath, J. E., & Paulucci Marinho, L. 2015, ArXiv eprints, arXiv:1504.03365
 Hotokezaka, K., Kiuchi, K., Kyutoku, K., et al. 2013, Phys. Rev. D, 88, 044026
 Hotokezaka, K., Kyutoku, K., Okawa, H., Shibata, M., & Kiuchi, K. 2011, Phys. Rev. D, 83, 124008
 Hulse, R. A., & Taylor, J. H. 1975, Astrophys. J. Lett., 195, L51
 Isenberg, J., & Nester, J. 1980, in General Relativity and Gravitation (Plenum Press, New York), 23
 Jaranowski, P., Kokkotas, K. D., Królak, A., & Tsegas, G. 1996, Classical and Quantum Gravity, 13, 1279
 Kastaun, W., & Galeazzi, F. 2015, Phys. Rev. D, 91, 064027
 Kiuchi, K., Kyutoku, K., Sekiguchi, Y., Shibata, M., & Wada, T. 2014, Phys. Rev. D, 90, 041502
 Kiuchi, K., Sekiguchi, Y., Shibata, M., & Taniguchi, K. 2009, Phys. Rev. D, 80, 064037
 Klähn, T., & Fischer, T. 2015, Astrophys. J., 810, 134
 Kochanek, C. S. 1992, Astrophys. J., 398, 234
 Kulkarni, S. R. 2005, ArXiv eprints, arXiv:astroph/0510256
 Lackey, B. D., & Wade, L. 2015, Phys. Rev. D, 91, 043002
 Lattimer, J. M. 2012, Annu. Rev. Nucl. Part. Sci., 62, 485
 Lee, W. H., RamirezRuiz, E., & van de Ven, G. 2010, Astrophys. J., 720, 953
 Li, L.X., & Paczyński, B. 1998, Astrophys. J. Lett., 507, L59
 Liu, Y. T., Shapiro, S. L., Etienne, Z. B., & Taniguchi, K. 2008, Phys. Rev. D, 78, 024012
 Lorimer, D. R. 2008, Living Reviews in Relativity, 11, 8
 Lugones, G., Benvenuto, O. G., & Vucetich, H. 1994, Phys. Rev. D, 50, 6100
 Lyford, N. D., Baumgarte, T. W., & Shapiro, S. L. 2003, Astrophys. J., 583, 410
 Madsen, J. 1988, Physical Review Letters, 61, 2909
 Metzger, B. D., Bauswein, A., Goriely, S., & Kasen, D. 2015, Mon. Not. R. Astron. Soc., 446, 1115
 Metzger, B. D., & Berger, E. 2012, Astrophys. J., 746, 48
 Metzger, B. D., MartínezPinedo, G., Darbha, S., et al. 2010, Mon. Not. R. Astron. Soc., 406, 2650
 Miller, J., Barsotti, L., Vitale, S., et al. 2015, Phys. Rev. D, 91, 062005
 Nissanke, S., Kasliwal, M., & Georgieva, A. 2013, Astrophys. J., 767, 124
 Oechslin, R., & Janka, H.T. 2007, Phys. Rev. Lett., 99, 121102
 Oechslin, R., Janka, H.T., & Marek, A. 2007, Astron. Astrophys., 467, 395
 Oechslin, R., Rosswog, S., & Thielemann, F.K. 2002, Phys. Rev. D, 65, 103005
 Olinto, A. V. 1987, Physics Letters B, 192, 71
 Oppenheimer, J. R., & Volkoff, G. M. 1939, Phys. Rev., 55, 374
 Paczynski, B. 1986, Astrophys. J. Lett., 308, L43
 Pagliara, G., Herzog, M., & Röpke, F. K. 2013, Phys. Rev. D, 87, 103007
 Palenzuela, C., Liebling, S. L., Neilsen, D., et al. 2015, Phys. Rev. D, 92, 044045
 Paschalidis, V., Etienne, Z. B., & Shapiro, S. L. 2012, Phys. Rev. D, 86, 064032
 Paschalidis, V., Liu, Y. T., Etienne, Z., & Shapiro, S. L. 2011, Phys. Rev. D, 84, 104032
 Radice, D., Rezzolla, L., & Galeazzi, F. 2014, Mon. Not. Roy. Astron. Soc., 437, L46
 Read, J. S., Baiotti, L., Creighton, J. D. E., et al. 2013, Phys. Rev. D, 88, 044042
 Rezzolla, L., & Zanotti, O. 2013, Relativistic Hydrodynamics (Oxford University Press, New York)
 Rodriguez, C. L., Farr, B., Raymond, V., et al. 2014, Astrophys. J., 784, 119
 Rosswog, S., Piran, T., & Nakar, E. 2013, Mon. Not. Roy. Astron. Soc., 430, 2585
 Ruffert, M., Janka, H.T., & Schaefer, G. 1996, Astron. Astrophys., 311, 532
 Sekiguchi, Y., Kiuchi, K., Kyutoku, K., & Shibata, M. 2011, Phys. Rev. Lett., 107, 051102
 Shibata, M. 2005, Phys. Rev. Lett., 94, 201101
 Shibata, M., & Taniguchi, K. 2006, Phys. Rev. D, 73, 064027
 Shibata, M., Taniguchi, K., & Uryū, K. 2005, Phys. Rev. D, 71, 084021
 Stairs, I. H. 2003, Living Reviews in Relativity, 6, 5
 Steiner, A. W., Hempel, M., & Fischer, T. 2013, Astrophys. J., 774, 17
 Stergioulas, N., Bauswein, A., Zagkouris, K., & Janka, H.T. 2011, Mon. Not. Roy. Astron. Soc., 418, 427
 Takami, K., Rezzolla, L., & Baiotti, L. 2014, Physical Review Letters, 113, 091104
 —. 2015, Phys. Rev. D, 91, 064001
 Taylor, J. H., & Weisberg, J. M. 1989, Astrophys. J., 345, 434
 The LIGO Scientific Collaboration, Aasi, J., Abbott, B. P., et al. 2015, Classical and Quantum Gravity, 32, 074001
 Tolman, R. C. 1939, Phys. Rev., 55, 364
 Typel, S., Röpke, G., Klähn, T., Blaschke, D., & Wolter, H. H. 2010, Phys. Rev. C, 81, 015803
 van der Sluys, M., Raymond, V., Mandel, I., et al. 2008, Classical Quantum Gravity, 25, 184011
 Veitch, J., Mandel, I., Aylott, B., et al. 2012, Phys. Rev. D, 85, 104045
 Veitch, J., Raymond, V., Farr, B., et al. 2015, Phys. Rev. D, 91, 042003
 Wade, L., Creighton, J. D. E., Ochsner, E., et al. 2014, Phys. Rev. D, 89, 103012
 Weisberg, J. M., Nice, D. J., & Taylor, J. H. 2010, Astrophys. J., 722, 1030
 Wilson, J. R., Mathews, G. J., & Marronetti, P. 1996, Phys. Rev. D, 54, 1317
 Witten, E. 1984, Phys. Rev. D, 30, 272
 Zhuge, X., Centrella, J. M., & McMillan, S. L. W. 1994, Phys. Rev. D, 50, 6247