Observing Gravitational Waves From The PostMerger Phase Of Binary Neutron Star Coalescence
Abstract
We present an effective, lowdimensionality frequencydomain template for the gravitational wave signal from the stellar remnants from binary neutron star coalescence. A principal component decomposition of a suite of numerical simulations of binary neutron star mergers is used to construct orthogonal basis functions for the amplitude and phase spectra of the waveforms for a variety of neutron star equations of state and binary mass configurations. We review the phenomenology of late merger / postmerger gravitational wave emission in binary neutron star coalescence and demonstrate how an understanding of the dynamics during and after the merger leads to the construction of a universal spectrum. We also provide a discussion of the prospects for detecting the postmerger signal in future gravitational wave detectors as a potential contribution to the science case for third generation instruments. The template derived in our analysis achieves match across a wide variety of merger waveforms and strain sensitivity spectra for current and potential gravitational wave detectors. A Fisher matrix analysis yields a preliminary estimate of the typical uncertainty in the determination of the dominant postmerger oscillation frequency as Hz. Using recently derived correlations between and the neutron star radii, this suggests potential constraints on the radius of a fiducial neutron star of m. Such measurements would only be possible for nearby ( Mpc) sources with advanced LIGO but become more feasible for planned upgrades to advanced LIGO and other future instruments, leading to constraints on the high density neutron star equation of state which are independent and complementary to those inferred from the premerger inspiral gravitational wave signal. We study the ability of a selection of future gravitational wave instruments to provide constraints on the neutron star equation of state via the postmerger phase of binary neutron star mergers.
pacs:
04.80.Nn, 07.05.Kf, 97.60.Jd, 04.25.dk1 Introduction
The second generation of gravitational wave (GW) observatories has now become operational with the first observations by advanced LIGO (aLIGO) underway [1]. Instruments such as advanced Virgo (advVirgo) [2] and KAGRA [3] will soon come online, eventually culminating in a worldwide network of GW observatories. The GW signal from the inspiral stage of binary neutron star (BNS) coalescence is amongst the most promising sources for this second generation of GW detectors. Observations of BNS GW inspiral signals from relatively nearby events (a few tens of Mpc) can lead to strong constraints on the supranuclear equation of state (EoS) via the impact on the phase evolution of the signal from tidal interactions during the latter stages of the merger [4, 5, 6, 7, 8, 9, 10] ^{1}^{1}1Reviews of the subject may also be found in [11, 12, 13, 14, 15].
For example, Read et al [5] showed that neutron star radii could be constrained with an uncertainty of for a single nearby (100 Mpc, assuming optimal orientation and skylocation) source, based on Fisher matrix estimates. More recently, a number of full Bayesian analyses have been carried out which have used astrophysicallymotivated simulated populations of BNS merger events to develop and probe EoS constraints in the low signaltonoise ratio (SNR) regime. In [16] it was found that only a few tens of inspiral events are required to measure neutron star (NS) tidal deformability to . Similar results are found in [17, 18], where the NS radius is determined to km and tidal deformability is determined to % after a few tens of GW detections.
The focus of this work, however, is on the independent and complementary constraints on the EoS which may be obtained from the postmerger signal. Depending on the mass configuration of the system and the EoS, a BNS merger may result in prompt collapse to a black hole (highmass, soft EoS) or the formation of a stable or quasistable NS remnant which again, may or may not collapse to a black hole depending on its mass and the EoS, while transient nonaxisymmetric deformations and quadrupolar oscillations in this remnant typically give rise to a richly structured, highfrequency (1–4 kHz) GW spectrum and a signal lasting ms, [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39]. Characterising the frequency content of GW signals from the postmerger system provides unique opportunities for GW asteroseismology: the dominant postmerger oscillation frequency exhibits a tight correlation with the radius of nonrotating neutron stars, with an overall uncertainty of a few hundred meters, depending on the total binary mass. For example, for a total binary mass of 2.7 the uncertainty in the radius of a cold, nonrotating NS of mass 1.6 (denoted as ) is about 220m [34, 40]. Similar relationships between the dominant spectral features and stellar parameters have been confirmed elsewhere [36, 38]. A deeper understanding of the features of postmerger GW spectra has been provided in [29, 41, 42], where it was shown that the spectrum is dominated by a linear feature (quadrupolar oscillations), a quasilinear feature (a coupling between quadrupolar and quasiradial oscillations) and a fully nonlinear feature (a transient spiral deformation), leading to a classification scheme of the postmerger GW emission depending on the EoS and binary mass. More recently, efforts have been made to find correlations between the pre and postmerger signals. In Ref. [43] the authors derive a relation between the tidal coupling constant that determines the tidal interactions before and during the merger and the peak frequency in the postmerger spectrum. Thus, measurements of the inspiral signal (which determine ) could be used to constrain by restricting its range of possible values and by combining measurements with those of the postmerger signal.
This connection between the tidal interactions and the postmerger oscillations highlights the complementarity of pre and postmerger GW observations. The constraints arising from inspiral observations may be subject to systematic biases induced by errors in the phase due to missing high PNorder terms or insufficiently accurate descriptions of spin or tidal effects. These systematic errors can be as large as the statistical uncertainty in characterising the inspiral signal [44, 18]. While inspiral waveform models will continue to improve and incorporate such effects, we note that analyses of the postmerger signal are subject to a completely independent source of systematic error (e.g., the precise relationship). Moreover, since the majority of premerger NSs are likely to have masses in the range , the premerger waveforms are limited to probing the structure of NSs in that massrange, while the postmerger signal allows us to probe the regime of higher masses (this is because, e.g. the central density of the remnant of a merger is close to the central density of a nonrotating star).
This high frequency component of the merger signal, however, will be somewhat challenging to observe in the upcoming generation of GW detectors. Typically, the most sensitive frequencies of groundbased GW instruments lie around 10–1000 Hz, with a rapidly diminishing sensitivity in the kHz regime. Additionally, the absence of a complete waveform model for the full pre and postmerger signal, or even for the postmerger signal alone, currently prohibits the use of matched filtering and one must turn to more robust, but ultimately less sensitive unmodelled burst searches. For example, the study in [45] revealed that a typical realistic burst analysis yielded an effective range approximately 30–40% of that which could be possible with an optimal matched filter.
Clearly then, there is great motivation and opportunity to develop more sensitive, more targeted analysis techniques and effective models which will bring us closer to the sensitivity offered by a matched filtering analysis. It is the goal of this work to explore a principal component analysis (PCA) based approach to constructing precisely such an effective model for the highfrequency component of the BNS merger signal. We construct a catalogue of 50 numerical waveforms from the merger and postmerger evolution of a variety of BNS systems with various EoSs and mass configurations. The magnitude and phase spectra of the waveforms in the catalogue are then decomposed into orthogonal bases using PCA. These basis functions can then be used to construct a frequencydomain waveform template which provides, on average, a 93% match for the waveforms in the catalogue for both aLIGO and a variety of potential upgrades and new GW instruments.
The structure of this paper is as follows: in section 2 we provide a detailed review of BNS merger and postmerger phenomenology, focussing on the resulting features in the GW spectrum and hence how one may constrain the NS EoS. Section 3 summarises the expected detectability of the highfrequency BNS waveforms used in this study, assuming a matchedfiltering approach and a variety of current and potential future GW instruments. In section 4 we describe and characterise our PCAbased frequencydomain waveform template in terms of waveform match and provide Fishermatrix estimates of the uncertainties in and based on this template as a preliminary guide to its potential. Finally, section 5 provides a summary and some concluding remarks relating to the planned applications of this model and the potential for similar approaches to enhance unmodelled burst analyses.
2 Properties of Postmerger GW Spectra and Constraining the Neutron Star EOS
2.1 Types of merger dynamics and GW spectra
For symmetric (i.e., equal component masses) and mildly asymmetric binaries the GW postmerger spectra of NS mergers (see e.g. right panel of Fig. 1 or Fig. 1 in Ref. [41]) show a generic behavior in the sense that certain features of the spectrum depend in a particular way on the total binary mass and the highdensity EoS [41]. Specifically, distinct peaks in the spectrum can be associated with distinct mechanisms generating those features, and the frequency and strength of the different GW peaks are determined by the total binary mass and EoS. The presence or absence of certain secondary peaks in the spectrum, together with their relative strengths is determined by the quasilinear coupling between the quasiradial and quadrupolar oscillation modes and by the orbital motion of antipodal bulges of a spiral deformation in the remnant. The characteristics of these distinct spectral features can be used to classify the postmerger dynamics of the system [41].
The most striking feature of the postmerger spectrum is a major peak generated by the dominant quadrupolar oscillation of the remnant, which is present in all models that form a NS merger remnant. The determination of the frequency of this peak in a GW measurement is the focus of this work because the peak frequency scales tightly with the radii of nonrotating NSs (see Fig. 4 below and discussion in [33, 46]) and thus provides strong constraints on the only incompletely known EoS of NS matter. Apart from this main peak, there can be up to two pronounced secondary peaks at frequencies below .
One of the secondary features is a peak generated by the quasilinear interaction between the dominant quadrupolar oscillation and the quasiradial mode of the remnant (the latter does not appear strongly in the GW spectrum on its own) [29]. The corresponding peak of the quasilinear mode coupling has an amplitude proportional to the product of the amplitudes of the quadrupolar mode and of the quasiradial mode, while its frequency, which we denote as , is equal to the difference of the frequencies of these two modes, i.e. , where is the frequency of the quasiradial mode. The feature is particularly pronounced for relatively high total binary masses and soft EoSs. Another secondary spectral peak is produced by the orbital motion of antipodal bulges, which form during the merging as a spiral deformation and then orbit around the inner remnant for a few milliseconds [41] (see Fig. 2). This dynamical feature is present in addition to the main emission at for the first few milliseconds after merging. Bulges moving with an orbital frequency result in a peak in the GW spectrum at . This finding receives further support by the timefrequency map of the GW signal shown in Fig. 1 for the 1.351.35 merger with the TM1 EoS [47, 48]. One can clearly recognize that in the early postmerger phase there are two distinct frequencies simultaneously contributing to the GW signal. The frequency of the dominant remnant oscillation is present for many milliseconds. The secondary peak at is generated within the first few milliseconds, when the antipodal bulges are pronounced (see Fig. 2). There is no evidence for a strong time variation of the frequencies, especially of the dominant frequency, which was suggested as an explanation for the structure of the GW spectrum in [49, 50].
The information in the timefrequency map of the GW signal can be related to the dynamical behavior of the remnant, which we illustrate by the evolution of the restmass density in the equatorial plane for the same simulation (see Fig. 2). The time step of the different snapshots are marked in the timefrequency map (Fig. 1) by vertical lines. Evidently, the presence of antipodal bulges at the outer remnant coincides with the presence of power at in the timefrequency map. It is apparent that the feature is initially particularly strong exceeding even the emission at ; the antipodal bulges are strongest during and immediately after merging and the spiral deformation forming the bulges initially comprises large parts of the remnant (see upper right panel in Fig. 2). In Fig. 2, the antipodal bulges complete approximately one orbit from the top right to the bottom left panel in about 1.2 ms. Thus, the orbital frequency kHz is expected to produce a peak at kHz, where a peak is found in the spectrum (see Fig. 1). For comparison we also show the timefrequency analysis for the SFHO EoS and component masses 1.351.35 in figure 3. Here the secondary peak at 2.2 kHz likely arises from the feature. An examination of the hydrodynamical data for this model reveals an of about 1.25 kHz (resulting in kHz), whereas the frequency of the quasiradial mode is kHz, and thus the peak is expected to occur at about 2.2 kHz.
The above findings on the timefrequency characteristics of the peak are consistent with the explanations of its origin presented in [41]. The feature is particularly strong for mergers with relatively low binary masses and stiff EoSs because less compact NSs favor the spiral deformation and the formation of the antipodal bulges during merging. In contrast, binaries with more massive components, i.e. very compact stars, merge with a higher impact velocity, which favours a strong excitation of the quasiradial mode of the remnant, leading to a strong feature, while the spiral deformation becomes less pronounced. For intermediate cases, i.e. moderately high binary masses, both secondary peaks are clearly present with comparable strength and distinguishable in frequency. Overall, this implies that for a given EoS the binary mass determines the presence and strength of the different secondary features. According to the classification scheme introduced in [41], one can identify three different types of spectra: Highmass/soft EOS binaries produce spectra where the dominant secondary peak is (Type I mergers). For intermediate binary masses and EOS stiffness, both the and features are present with roughly comparable amplitude (Type II mergers). Lowmass/stiff EOS binaries produce spectra with a strong peak and an absent feature (Type III mergers). See [41, 42] for further discussion.
2.2 Universal Postmerger Spectra & Measuring The Neutron Star Radius
For a fixed total binary mass the frequencies of the three different peaks depend in a particular way on the EoS, which can be characterized by the radius or compactness of nonrotating NSs [33, 46, 41] (see also [51, 50] for the dependence of the strongest secondary feature on compactness, without distinguishing the different nature of secondary peaks). The two secondary frequencies show a tight correlation with the dominant postmerger frequency . This is shown in Fig. 4 for 1.351.35 mergers.
The existence of generic spectral features with predictable behavior suggests that the construction of a universal spectrum should be feasible through the appropriate alignment of the main peaks from spectra for various EoSs. Figure 5 shows the GW spectra for equal mass binaries (1.351.35) with different EoSs. The waveforms have been normalised such that the rootsumsquared amplitude is unity,
(1) 
In the right panel, we rescale the frequency axis such that the dominant quadrupolar oscillation peak feature is located at a common reference value (2.6 kHz) for all models. Apart from small variations of the secondary features a remarkable universality of the spectra is found, which can be explained as follows: We choose a reference peak frequency of kHz. Thus, for a spectrum with the main peak at the factor for rescaling the frequency is . This factor is also applied to the frequencies of the secondary peaks. Therefore, a rescaled secondary peak (i.e. or ) is located at . Since the fraction is approximately constant and similar for both secondary features (see Fig. 4), a rescaled secondary feature occurs at approximately the same frequency for all models.
This universality of the scaled spectra suggests that it should be possible to produce a model from the mean spectrum, computed over a number of numerical simulations, plus some small deviations. In Sec. 4, we demonstrate that PCA provides an approach to solve exactly this problem by producing an orthornomal basis constructed from a superposition of the meancentered spectra. Furthermore, we find that the perturbations from the mean spectrum are generally well described by a small number of basis functions.
It is important to stress here that, for distances which allow for the detection of postmerger GW emission, the individual masses of the binary components can be determined with an accuracy of a few per cent [52, 53, 54]. Furthermore, current observations suggest that BNS mass configurations will not be dramatically asymmetric (see e.g. the compilation of NS masses in [55]). Peak frequencies which are recovered within this data analysis study, are converted to NS radii via
(2) 
describing the empirical relation between NS radii and the dominant postmerger oscillation frequency for symmetric mergers for total binary masses of 2.7 ^{2}^{2}2We note that, even in the abseence of a measurement of mass ratio, the relation is quite robust for a constant chirp mass, which is generally recovered to high precision. See e.g., [42]. Here we adopt the coefficients , and from previous work [40]. For this particular total binary mass the maximum uncertainty in the empirical relation is 175 m. In addition to this systematic error, the measurement of in noisy GW data will introduce a statistical error, whose determination is one of the major goals of this work and is quantified in Sect. 4.
2.3 Binary Neutron Star Merger Waveforms Used In This Study
The numerical waveforms used in this study rely mostly on the calculations discussed in [33, 46, 40, 45, 41], where further information can be found. Additional waveform models employed here are obtained within the same physical and numerical model, for which further details are provided in [56, 23, 57, 34]. The EoS models for the hydrodynamical simulations are chosen to cover a large variety including very stiff and very soft EoSs (see Table 1) and a variety of binary mass configurations are used. All EoSs are compatible with a maximum NS mass of [58, 59]. Figure 6 illustrates the selection of waveforms used in this study in terms of their EoSs and mass configuration.
EoS  Ref.  [km] 

NL3  [60, 61]  14.75 
LS375  [62]  13.65 
DD2  [63, 61]  13.21 
TM1  [47, 48]  14.49 
SFHX  [64]  11.98 
GS2  [65]  13.38 
SFHO  [64]  11.92 
LS220  [62]  12.73 
TMA  [66, 48]  13.86 
APR  [67]  11.33 
BHBLP  [68]  13.21 
Shen  [69]  14.64 
Heb6  [70]  13.33 
Heb5  [70]  12.38 
Heb4  [70]  12.51 
Heb3  [70]  12.03 
3 Detectability
We now discuss the expected detectability of the postmerger GW signal in current and planned GW instruments. A natural, preliminary, measure of detectability is the matchedfilter SNR one would obtain given a perfect model, or template, for the signal waveform in GW detector data . The matchedfilter SNR is defined as,
(3) 
The optimal SNR, where the template exactly matches detector output is then simply,
(4) 
where is the usual inner product [71]:
(5) 
and is the noise spectrum of a given GW detector and the asterix indicates complex conjugation. Note that we impose a lower bound on the frequency over which the inner product is evaluated in order to target the detectability of the highfrequency part of the signal. In this study we use kHz. The inner product is evaluated up to the Nyquist frequency of the spectrum, 8192,Hz in this study. We also characterise detectability in terms of horizon distance : the distance at which an optimally oriented source yields an SNR at least as large as some nominal threshold, . For GW searches in which the time of arrival of the signal and the source skylocation are unknown, it is typical to evaluate horizon distances with . In our application, however, we envisage a hierachical ‘triggered’ analysis, similar to that described in [45], wherein the earlier, lowerfrequency inspiral portion of the coalescence signal has already been detected at high confidence. It is likely then that the time of coalescence has been determined to an accuracy of a few or a few 10’s of milliseconds and we can significantly reduce the threshold used to define the horizon distance. Following [45], we choose . Finally, we can determine the rate with which we will obtain signals with from the expected number of BNS mergers which are accessible to a search with a given horizon distance [72]. For the purposes of this study, we assume the ‘realistic’ rate of BNS coalescence from [72]: MWEGMyr.
We now compute each of the figures of merit (the SNR for an optimally oriented source at 50 Mpc; the horizon distance assuming an SNR threshold and the expected detection rate ) for aLIGO [73, 1], as well as the following selection of proposed upgrades to aLIGO and new facilities. The following descriptions emphasise the expected increases in sensitivity relative to aLIGO only over 1–4 kHz; the band of interest for the postmerger signal. Comparisons with the increased range and sensitivity to the earlier inspiral part of the signal are left to future studies. Note also that we take aLIGO to be the most sensitive of the second generation GW detectors; instruments such as advVirgo and Kagra offer comparable or reduced sensitivity in the frequency regime of interest to this study. It should be noted, however, that a network of detectors with comparable sensitivity could improve the range of an search by a factor of up to with respect to the single detector expectation, assuming stationary Gaussian noise and an optimal analysis. We restrict our estimates to single detector ranges and rates in the interests of conservatism and simplicity.
 LIGO A+ [74, 75]

a set of upgrades to the existing LIGO facilities, including frequencydependent squeezed light, improved mirror coatings and potentially increased laser beam sizes. Noise amplitude spectral sensitivity would be improved by a factor of over 1–4 kHz. A+ could begin operation as early as 2017–18.
 LIGO Voyager (LV) [75]

a major upgrade to the existing LIGO facilities, including higher laser power, changes to materials used for suspensions and mirror substrates and, possibly, low temperature operation. LV would become operational around 2027–28 and offer noise amplitude spectral sensitivity improvements of over 1–4 kHz.
 LIGO Cosmic Explorer (CE) [75]

a new LIGO facility rather than an upgrade, with operation envisioned to commence after 2035, probably as part of a network with LIGO Voyager. In its simplest incarnation, Cosmic Explorer would be a straightforward extrapolation of A+ technology to a much longer arm length of 40 km, referred to as CE1 which would be more sensitive than aLIGO over 1–4 kHz. An alternative extrapolation is that of Voyager technology to the 40 km arm length, referred to as CE2. CE2 is only more sensitive than aLIGO for the frequency range of interest in this study. For simplicity, we consider only CE1.
 Einstein Telescope (ETD) [76, 77]

the European thirdgeneration GW detector. In this work, we consider the ETD configuration which is comprised of two individual inteferometers where one targets low frequency sensitivity and the other high frequency sensitivity. Both interferometers will be of 10 km arm length and housed in an underground facility. Furthermore, the full observatory will consist of three such detectors in a triangle arrangement. ETD is more sensitive than aLIGO over 1–4 kHz. Due to the network configuration (i.e., the alignment of the component instruments) the effective sensitivity of ETD is % higher than that for a single ETD detector.
Figure 7 shows the design sensitivity spectra for each of these instruments, again focussing on the 1–4 kHz range of interest for the postmerger BNS GW signal. For comparison, we also show the amplitude spectrum of a typical BNS waveform (the TM1 1.35+1.35 example discussed in section 2) for an optimally oriented source at 50 Mpc. Finally, the figures of merit describing the detectability of the postmerger signal for each instrument are summarised in table 2. Note that, we compute two measures of SNR: SNR, where we simply evaluate equation 5 over 1–4 kHz for the full merger waveform; as well as SNR, where the timedomain waveform has been windowed to suppress power prior to the merger (taken to occur at the peak strain amplitude), in order to yield an estimate of the contribution to the SNR from the postmerger oscillations. Since there are 5 instruments, 50 waveforms and 4 figures of merit, we choose to summarise the results for each instrument in terms of the 10, 50 and 90 percentiles, evaluated over the 50 waveforms used in the study.
Instrument  SNR  SNR  [Mpc]  [year] 

aLIGO  2.99  1.48  29.89  0.01 
A+  7.89  4.19  78.89  0.13 
LV  14.06  7.28  140.56  0.41 
ETD  26.65  12.16  266.52  2.81 
CE  41.50  20.52  414.62  10.59 
4 A Waveform Model Using Principal Component Analysis
The optimal data analysis method for the identification and characterisation of a GW signal in noisy data is matchedfiltering, wherein an exact analytic model, or template, for the waveform is convolved with data stream from a network of GW detectors. Unfortunately, the physical complexity of the merging binary neutron star system is such that detailed numerical simulations are required to produce even an approximate waveform. Furthermore, since the physical parameters of the system are essentially unknown, many such simulations would be required in order to build a template bank to maximise the likelihood of signal detection.
We are, therefore, confronted with a similar data analysis problem to that in the analysis of GWs from core collapse supernovae: the absence of an accurate analytic waveform template, a limited number of computationally expensive and approximate simulations and a requirement to significantly reduce the complexity of the modelling problem to faciliate the use of an approximate matchedfilter. Motivated by the work in [78, 79], we find that we can construct an effective waveform model from a basis constructed using PCA of a suite of merger simulations comprised of systems with different equations of state, masses and mass ratios.
Our goal is to reduce the complexity of the modelling problem from a highdimensional physical parameter space, where the waveforms are modelled directly through numerical simulation, to a lowerdimensional problem to model the dominant features of the waveform. PCA of a catalogue of simulated waveforms provides a solution to precisely this problem. Denoting the timedomain merger waveform as , its complex Fourier spectrum is given by,
(6) 
where and are the magnitude and phase spectra of signal , respectively, In a similar spirit to the approach described in [80] we construct orthnormal bases for the amplitude and phase spectra separately, using similar a PCA decomposition to that described in [78, 79]. Principal component analysis forms a basis from the eigenvectors of the covariance matrix of some set of data. The procedure is as follows:

Collate a representative sample of binary merger GW waveforms, sampled at 16384 Hz. This sample of waveforms is hereafter referred to as our training catalogue. Each waveform is normalised to unit rootsumsquared amplitude (equation 1) to reduce catalogue variance from different amplitude scales and emphasise morphological differences.

Compute the complex Fourier spectra of the timedomain waveforms in our catalogue. A Tukey window is initially applied to the timedomain signals to minimise spectral ringing and the the waveforms are zeropadded to a uniform samples. The complex spectra are computed using the fast Fourier transform. The amplitude and phase spectra are computed from the absolute values and arguments of the complex frequency series and the phase spectra are unwrapped to yield smooth functions, each of samples.

The unique feature to the analysis presented in this work is our choice of feature alignment, an absolutely key component to PCA. In [78], for example, the GW waveforms are aligned such that the peak amplitudes lie at a common reference time removing the need for the PCA to account for trivial variance in the catalogue. The analogous procedure in our application is to align features in the frequency domain. Each amplitude spectrum is rescaled such that the dominant postmerger peak, labelled , is aligned to a common reference value . This alignment is achieved by computing a set of frequencies , where are the angular frequencies of the original spectrum. We then interpolate the original spectrum to the new frequencies where the dominant spectral feature (the postmerger oscillation peak) is aligned. Although it is not perfect, this geometric scaling (as opposed to a simple linear shift) also helps to align the subdominant and features. Three examples of original and aligned amplitude spectra are shown in figure 8.
Figure 8: Left: Example magnitude spectra; Right: Example spectra after alignment to a common peak frequency . 
Next, we construct an matrix where each row correponds to the sample featurealigned amplitude spectrum of each waveform after subtracting the mean spectrum (averaged over the waveforms)^{3}^{3}3Note that this matrix is transposed relative to the descriptions in [78, 79] for more straightforward comparison with other PCA literature and use with software packages such as that offered in [81].. The mean amplitude spectrum, evaluated over our waveforms is shown in the left panel of figure 9.

Finally, we perform the PCA decomposition in which we compute the eigenvectors of the empirical covariance matrix . Following [79] and noting the change in row/column convention for the data matrix, the centered data matrix , of dimension , can be factorised using singular value decomposition,
(7) where and are orthonormal matrices with dimensions and , respectively; is a diagonal matrix of singular values of in descending order and . The columns of , , contain the eigenvectors of the covariance matrix , our principal components, and the singular values in are the square roots of the eigenvalues of the covariance matrix . Finally, the columns of contain the eigenvectors of . The principal components , comprise an orthonomal basis of the rows (i.e., the aligned and centered waveforms) in so that each of the aligned waveform amplitude (or phase) spectra can be represented as a linear sum of principal components and the mean. For example, the aligned amplitude spectrum of the first waveform can be constructed as:
(8) where is the mean amplitude spectrum over the aligned waveform catalogue and are weighting coefficients given by the projection of the centered onto the principal component basis,
(9) and are the elements of . The original waveform magnitude spectrum is, at last, obtained by applying the inverse of the alignment procedure in step (3). Figure 9 shows the mean aligned magnitude spectrum and the first principal component as computed for the 50 waveforms described in section 2.3. Note that we have chosen to align the dominant postmerger peak to a value of 2710 Hz; this choice is essentially arbitrary and simply corresponds to the mean of the peak frequencies in the catalogue. It is important to note here that the value of is a free parameter in the spectral model; in practice, its value must be inferred from GW observations.
Figure 9: Left: Mean magnitude spectrum; Right: First principal component for the magnitude spectrum
Ultimately, our goal is to construct a reduced basis from which any postmerger waveform can be reconstructed to some accuracy. It is, therefore, helpful to understand the relative importance of each principal component. A measure of the total variance in the centered catalogue is given by the trace of the covariance matrix . The variance explained by principal components is, then, the sum of the first eigenvalues:
(10) 
where are the singular values from equation 7. The postmerger waveforms can then be approximated by using a reduced basis with , with the choice of based on capturing a reasonable degree of variance in the catalogue, and equation 8. Figure 10 shows the cumulative explained variance for both the magnitude and phase spectra of the waveforms in our catalogue (i.e., equation 10) as a function of the number of principal components. One can immediately see that the variation between the waveforms is dominated by the rich and varied structure in the magnitude spectra; only of the total variance is explained by the first principal component of the magnitude spectra, while of the variance in the phase spectra is described by the first component.
4.1 PCA Templates: Characterisation & Expected Performance
Remembering that our goal is to build an approximate waveform template for matched filtering, a useful figure of merit to characterise PCA–based model is the waveform match , which describes the fraction of the optimal signaltonoise ratio for a given signal which is captured by the waveform template :
(11) 
where is the innerproduct, defined by equation 5, maximised over the start time and initial phase offset of the signal. The match is normalised such that for a perfect template and zero for an template which is orthogonal to the target signal. In the following examples, the match is computed assuming the aLIGO noise curve. Figure 11 shows an example of the reconstructed time series and magnitude spectrum for the TM1 1.35+1.35 system considered earlier in section 2. The time series is given by the inverse Fourier transform of the complex spectrum constructed from the separate amplitude and phase PCA. As expected, when we use the full PCA basis with and include the waveform in the data matrix , we obtain a complete basis which allows a perfect reconstruction such that .
It is unlikely that nature will provide us with a signal which exactly matches one of those contained in the set of training data . The right panel of figure 11 again shows the original and reconstructed TM1 1.35+1.35 waveform, except now this waveform has been excluded from the training data. In addition, we use only the first principal components in amplitude and phase. With this more realistic example and a much smaller parameter space, we are still able to reconstruct the target signal with a match ^{4}^{4}4Again, the value of is assumed known here; the match here represents the the best case scenario.
We now compute similar matches for all of the waveforms in our catalogue and for the different instruments described in § 3. To begin, we compute match using the aLIGO noise curve as a function of the number of principal components used and compare the results of including and excluding each waveform from the data used to compute the PCA. These results are summarised in figure 12 with the mean, minimum, maximum and the tenth and ninetieth percentiles over the matches computed for each of the fifty waveforms. The left panel shows the results when all of the waveforms are used while the right panel summarises the matches when each waveform is removed from the catalogue prior to computing the PCA. We see that, as before, perfect reconstruction fidelity is attained using the full basis when all waveforms are used. In contrast, the match remains approximately constant with respect to the number of principal components used when each waveform being matched is excluded from the training data. This is a reflection of the fact that the lowerorder principal components represent the most common generic features in the catalogue, while the higherorder components are essentially minor corrections to the mean which may not be present in the waveform which is excluded. Given the quite respectable matches obtained with just the first principal component and its apparent robustness, we propose modelling the highfrequency GW spectrum for binary neutron star mergers using equation 8 with .
We now repeat the match calculation for each of the instrument noise curves described in § 3, using just the first principal component. The 10, 50 (i.e., the median) and 90 percentiles, computed over the matches for different waveforms, are summarised in figure 13 and listed explicitly in table 3. We find that the PCA templates yield a match of across all of the instruments considered. Variations in the match arise from differences in the shapes of the noise curve, i.e., the denominator in equation 5; in the kHz regime, where sensitivity is limited by photon shotnoise, the noise curves mostly only differ in their overall amplitude scale and we do not expect significant variations in match quality.
4.2 Implications For Parameter Estimation
Given this approximate waveform template it is useful to determine its effectiveness in parameter estimation and, ultimately, the extraction of astrophysics. Recall from § 2 that the single most robust feature of the GW spectrum for these signals is the presence of a dominant spectral peak due to excitation of the postmerger remnant’s quadrupolar mode oscillation. Recall also that the peak frequency of this excitation in systems with total binary masses of 2.7 correlates strongly with the radius of a fiducial nonspinning NS across a wide variety of equations of state. Our goal then, is to determine how accurately we might expect to measure using our PCA–based waveform model. In this section, we review the waveform template and derive Fishermatrix estimates for the accuracy with which we may measure and hence given current and planned GW observatories.
The aligned magnitude spectrum at frequency of our waveform model is given by
(12) 
where the spectrum is aligned such that the dominant peak lies at frequency , is the mean magnitude spectrum of the catalogue, is the first principal component and is the coefficient an arbitrary waveform’s projection onto . The final magnitude spectrum is given by interpolating the aligned spectrum to a set of new frequencies :
(13) 
and an identical procedure is applied to the phase spectrum . In this prescription then, the peak frequency is a direct parameter of the model.
We estimate the expected accuracy of the estimation from the Fisher matrix and by considering a oneparameter family of waveform templates in which only varies. We assume for simplicity that other important parameters, such as the start time of the signal, have already been determined we hold the value of at its nominal value. While does indeed play a role in determining the detailed shape of the spectrum and, particularly, the degree of asymmetry in the main spectral peak, its will not have a strong correlation with the location of the maximum. A more detailed and realistic Bayesian analysis will be conducted in the near future to account for possible correlations.
The expected error in some parameter can be determined from the Fisher matrix [71]:
(14) 
Following the procedure in [5], we estimate the error in to first order from,
(15) 
where and are our waveform templates evaluated at peak frequencies and , which lie below and above the true , respectively and Hz. We find that this expansion is quite stable to the choice of , such that varies by only a few percent up to Hz.
Figure 14 and table 3 summarise the expected frequency errors obtained across waveforms and instruments using the Fisher matrix estimate. Errors are evaluated at SNR=5, corresponding to a source at the horizon distance. As with the match summary from earlier the results for each instrument are summarised with the 10, 50 and 90 percentiles computed over the different waveforms. Again, the expected frequency error is fairly consistent between the different instruments and we find Hz. We can propagate the expected error in the determination to that in the NS radius using equation 2:
(16)  
(17) 
The errors thus obtained represent the statistical uncertainty in the radius, arising from the measurement of a signal in noisy data. The fit given by equation 2 is also subject to a systematic error which, as described in section 2.2, we take to be the maximum deviation in the relationship across a variety of EoSs, m. To arrive at a total expected error in the determination of the radius then, we quote the quadrature sum of the statistical and systematic errors:
(18) 
The expected radius errors thus obtained are summarised in the right panel of figure 14 and in table 3 with the usual breakdown by instrument and percentile summary statistics.
Instrument  [Hz]  [m]  [m]  

aLIGO  
A+  
LV  
CE  
ETD 
5 Summary & Outlook
The wealth of information contained in the highfrequency spectrum of BNS mergers means that there is a strong motivation to develop effective models for the merger and postmerger phase of the coalescence GW signal. Through consideration of the general morphology of the postmerger spectrum and the phenomenology during and after the merger, we have determined that the highfrequency complex spectrum is remarkably well modelled by an orthogonal basis constructed from a catalogue of numerical simulations using a PCA decomposition.
Typically, the waveform templates thus constructed yield a match of , over the frequency range 1–4 kHz, with the majority of the waveforms used in this study. While the typical desideratum in most matchedfiltering analyses is , it is worth noting that the only other systematic and wellquantified estimate of GW search effectiveness to date has been the burst analysis reported in [45]. While the burst analysis is robust to uncertainties in the waveform it was found that its effective range was only that of an optimal matched filter analysis. The PCA model presented in this work therefore holds the potential to double or even triple the sensitivity offered by existing analyses. Furthermore, a preliminary Fisher matrix analysis reveals that the uncertainty in the determination of the peak postmerger oscillation frequency is Hz, implying a statistical uncertainty on the radius of a 1.6 NS of m for sources with sufficient power at 1–4 kHz or proximity to Earth to produce SNR=5. Assuming a conservative estimate of for the systematic error in the relation when the binary masses are known, we find the total error in the radius is m. For comparison, the analysis in [45] found m. Both cases assume an relationship appropriate for a symmetric mass configuration with total mass . Note that a) the estimate in [45] included only the statistical error and b) the burst analysis requires a relatively large SNR before sufficient GW signal power is acquired to generate a detection candidate; by this time, the peak frequency itself can be quite easily resolved. Furthermore, our estimates here are based on a singledetector analysis; those in [45] considered a threedetector network operating with comparable sensitivity in each instrument.
For aLIGO the horizon distance with an optimal template for SNR=5 is generally Mpc with a plausible signal rate of approximately 1 event per 100 years, comparable to the rate of Galactic supernovae (see Table 2). Our template, however, will lose % of this SNR, resulting in a proportional decrease in the horizon distance. In fact, thanks to the local overdensity of galaxies, the impact on signal rate from this mismatch is rather negligible (%) until we consider the ETD or CE sensitivities.
In addition to the obvious benefit of potentially yielding a greater detection horizon than an unmodelled search, it is worth highlighting a number of other advantages of using even a rather ad hoc waveform template such as ours. The strain sensitivity spectrum of GW detectors and the short duration of postmerger GW signals suggest that the premerger inspiral signal will always be observed at high SNR for any source which is sufficiently close to observe the postmerger signal. This will lead to a quite precise determination of the time of coalescence, potentially a constraint on the skylocation and constraints on the binary masses. If we also assume that the merger does indeed result in the formation of a stable or quasistable NS remnant then the analysis need only consist of inferring the parameters of our waveform model and does not necessarily require a signal to produce an SNR above some detection threshold. Instead, our estimate of the parameters (e.g., ) simply have greater uncertainty for low SNR signals. This is in contrast to typical burst analyses which require signals to be sufficiently loud that they produce statistically significant loud pixels in the timefrequency plane. We note, however, that a timefrequency PCA could very easily be used to ‘inform’ burst clustering algorithms to better target signals such as these where there is a rich timefrequency structure.
For the purposes of this study, we adopted an SNR threshold as a fiducial point of reference; in practice, however, it may be possible to determine at larger distances (although with correspondingly greater uncertainty) than suggested in this work. Since the Fisher matrix approximation is not valid in the low SNR regime, this point will be investigated via a full Bayesian analysis using our PCA templates in a future study.
Finally, it is worth mentioning that the construction of PCA templates leads to an intriguing and natural way to feed GW observations back to the numerical modelling community to produce a feedback loop for the estimation of the NS EoS and refinement of our waveform models. For now we simply sketch the basic algorithm as follows, leaving an example implementation to future work:

Construct a PCA template from a coarselysampled catalogue of merger waveforms, whose span the full frequency space permitted by allowed EoSs and masses.

Following a nearby BNS detection, determine the probable component masses from the inspiral signal and the best estimate of from the PCA template constructed in (i).

Produce a refined, more finelysampled catalogue of merger waveforms (with new simulations if necessary) which correspond only to those EoSs and mass configurations which are compatible with the observations in (ii). Construct a new PCA template from this catalogue.
This process could then be iterated until some desirable stopping criterion, such as reaching some critical value of the minimum match between the PCA templates and waveforms used, is reached. This approach may provide an avenue to go beyond simply determining and allow an accurate reconstruction of the full spectrum of the underlying signal. We see then that the application of PCA to construct robust and simple phenomenological templates for the characterisation of postmerger BNS signals holds great promise on its own, and may provide a useful tool to augment other approaches such as those in [45, 43, 42].
Acknowledgements
The authors thank Tjonnie Li for helpful input and careful reading of this manuscript. J. C. and D. S. gratefully acknowledge support from NSF grants PHY0955825, PHY1212433, PHY1333360, PHY1505824 and PHY1505524. A. B. is a Marie Curie IntraEuropean Fellow within the 7th European Community Framework Programme (IEF 331873). 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.
References
References
 [1] The LIGO Scientific Collaboration, J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Abernathy, K. Ackley, C. Adams, T. Adams, P. Addesso, and et al. Advanced LIGO. Classical and Quantum Gravity, 32(7):074001, April 2015.
 [2] F. Acernese et al. Advanced Virgo: a secondgeneration interferometric gravitational wave detector. Class. Quant. Grav., 32(2):024001, 2015.
 [3] Y. Aso, Y. Michimura, K. Somiya, M. Ando, O. Miyakawa, T. Sekiguchi, D. Tatsumi, and H. Yamamoto. Interferometer design of the KAGRA gravitational wave detector. Phys. Rev. D., 88(4):043007, August 2013.
 [4] É. É. Flanagan and T. Hinderer. Constraining neutronstar tidal Love numbers with gravitationalwave detectors. Phys. Rev. D., 77(2):021502, January 2008.
 [5] J. S. Read, C. Markakis, M. Shibata, K. Uryū, J. D. E. Creighton, and J. L. Friedman. Measuring the neutron star equation of state with gravitational wave observations. Phys. Rev. D., 79(12):124033, June 2009.
 [6] T. Hinderer, B. D. Lackey, R. N. Lang, and J. S. Read. Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral. Phys. Rev. D., 81(12):123016, June 2010.
 [7] L. Baiotti, T. Damour, B. Giacomazzo, A. Nagar, and L. Rezzolla. Analytic Modeling of Tidal Effects in the Relativistic Inspiral of Binary Neutron Stars. 105(26):261101, December 2010.
 [8] S. Bernuzzi, A. Nagar, M. Thierfelder, and B. Brügmann. Tidal effects in binary neutron star coalescence. Phys. Rev. D., 86(4):044030, August 2012.
 [9] T. Damour, A. Nagar, and L. Villain. Measurability of the tidal polarizability of neutron stars in lateinspiral gravitationalwave signals. Phys. Rev. D., 85(12):123007, June 2012.
 [10] J. S. Read, L. Baiotti, J. D. E. Creighton, J. L. Friedman, B. Giacomazzo, K. Kyutoku, C. Markakis, L. Rezzolla, M. Shibata, and K. Taniguchi. Matter effects on binary neutron star waveforms. Phys. Rev. D., 88(4):044042, August 2013.
 [11] M. D. Duez. Numerical relativity confronts compact neutron star binaries: a review and status report. Classical and Quantum Gravity, 27(11):114002, June 2010.
 [12] N. Andersson, V. Ferrari, D. I. Jones, K. D. Kokkotas, B. Krishnan, J. S. Read, L. Rezzolla, and B. Zink. Gravitational waves from neutron stars: promises and challenges. General Relativity and Gravitation, 43:409–436, February 2011.
 [13] T. W. Baumgarte and S. L. Shapiro. Numerical Relativity: Solving Einstein’s Equations on the Computer. June 2010.
 [14] Joshua A. Faber and Frederic A. Rasio. Binary neutron star mergers. Living Reviews in Relativity, 15(8), 2012.
 [15] L. Rezzolla and O. Zanotti. Relativistic Hydrodynamics. September 2013.
 [16] Walter Del Pozzo, Tjonnie G. F. Li, Michalis Agathos, Chris Van Den Broeck, and Salvatore Vitale. Demonstrating the feasibility of probing the neutronstar equation of state with secondgeneration gravitationalwave detectors. Phys. Rev. Lett., 111:071101, Aug 2013.
 [17] M. Agathos, J. Meidam, W. Del Pozzo, T. G. F. Li, M. Tompitak, J. Veitch, S. Vitale, and C. Van Den Broeck. Constraining the neutron star equation of state with gravitational wave signals from coalescing binary neutron stars. Phys. Rev. D., 92(2):023012, July 2015.
 [18] B. D. Lackey and L. Wade. Reconstructing the neutronstar equation of state with gravitationalwave detectors from a realistic population of inspiralling binary neutron stars. Phys. Rev. D., 91(4):043002, February 2015.
 [19] X. Zhuge, J. M. Centrella, and S. L. W. McMillan. Gravitational radiation from coalescing binary neutron stars. Phys. Rev. D., 50:6247–6261, November 1994.
 [20] M. Ruffert, H.T. Janka, and G. Schaefer. Coalescing neutron stars  a step towards physical models. I. Hydrodynamic evolution and gravitationalwave emission. Astron. Astrophys., 311:532–566, July 1996.
 [21] M. Shibata. Constraining Nuclear Equations of State Using Gravitational Waves from Hypermassive Neutron Stars. 94(20):201101, May 2005.
 [22] M. Shibata, K. Taniguchi, and K. Uryū. Merger of binary neutron stars with realistic equations of state in full general relativity. Phys. Rev. D., 71(8):084021, April 2005.
 [23] R. Oechslin, H.T. Janka, and A. Marek. Relativistic neutron star merger simulations with nonzero temperature equations of state. I. Variation of binary parameters and equation of state. Astron. Astrophys., 467:395–409, May 2007.
 [24] R. Oechslin and H.T. Janka. Gravitational Waves from Relativistic NeutronStar Mergers with Microphysical Equations of State. 99(12):121102, September 2007.
 [25] M. Anderson, E. W. Hirschmann, L. Lehner, S. L. Liebling, P. M. Motl, D. Neilsen, C. Palenzuela, and J. E. Tohline. Simulating binary neutron stars: Dynamics and gravitational waves. Phys. Rev. D., 77(2):024006, January 2008.
 [26] Y. T. Liu, S. L. Shapiro, Z. B. Etienne, and K. Taniguchi. General relativistic simulations of magnetized binary neutron star mergers. Phys. Rev. D., 78(2):024012, July 2008.
 [27] L. Baiotti, B. Giacomazzo, and L. Rezzolla. Accurate evolutions of inspiralling neutronstar binaries: Prompt and delayed collapse to a black hole. 78:084033, Oct 2008.
 [28] K. Kiuchi, Y. Sekiguchi, M. Shibata, and K. Taniguchi. Longterm general relativistic simulation of binary neutron stars collapsing to a black hole. Phys. Rev. D., 80(6):064037, September 2009.
 [29] N. Stergioulas, A. Bauswein, K. Zagkouris, and H.T. Janka. Gravitational waves and nonaxisymmetric oscillation modes in mergers of compact object binaries. Mon. Not. Roy. Astron. Soc., 418:427–436, November 2011.
 [30] B. Giacomazzo, L. Rezzolla, and L. Baiotti. Accurate evolutions of inspiralling and magnetized neutron stars: Equalmass binaries. Phys. Rev. D., 83(4):044014, February 2011.
 [31] K. Hotokezaka, K. Kyutoku, H. Okawa, M. Shibata, and K. Kiuchi. Binary neutron star mergers: Dependence on the nuclear equation of state. Phys. Rev. D., 83(12):124008, June 2011.
 [32] Y. Sekiguchi, K. Kiuchi, K. Kyutoku, and M. Shibata. Gravitational Waves and Neutrino Emission from the Merger of Binary Neutron Stars. prl, 107(5):051102, July 2011.
 [33] A. Bauswein and H.T. Janka. Measuring NeutronStar Properties via Gravitational Waves from NeutronStar Mergers. Physical Review Letters, 108(1):011101, January 2012.
 [34] A. Bauswein, H.T. Janka, K. Hebeler, and A. Schwenk. Equationofstate dependence of the gravitationalwave signal from the ringdown phase of neutronstar mergers. Phys. Rev. D., 86(6):063001, September 2012.
 [35] A. Bauswein, T. W. Baumgarte, and H.T. Janka. Prompt Merger Collapse and the Maximum Mass of Neutron Stars. Phys. Rev. Lett., 111(13):131101, September 2013.
 [36] K. Hotokezaka, K. Kiuchi, K. Kyutoku, T. Muranushi, Y. Sekiguchi, M. Shibata, and K. Taniguchi. Remnant massive neutron stars of binary neutron star mergers: Evolution process and gravitational waveform. Phys. Rev. D., 88(4):044026, August 2013.
 [37] S. Bernuzzi, T. Dietrich, W. Tichy, and B. Bruegmann. Mergers of binary neutron stars with realistic spin. ArXiv eprints, November 2013.
 [38] K. Takami, L. Rezzolla, and L. Baiotti. Constraining the equation of state of neutron stars from binary mergers. ArXiv eprints, March 2014.
 [39] T. W. Baumgarte, S. L. Shapiro, and M. Shibata. On the Maximum Mass of Differentially Rotating Neutron Stars. Astrophys. J. Lett., 528:L29–L32, January 2000.
 [40] A. Bauswein, N. Stergioulas, and H.T. Janka. Revealing the highdensity equation of state through binary neutron star mergers. Phys. Rev. D., 90(2):023002, July 2014.
 [41] A. Bauswein and N. Stergioulas. Unified picture of the postmerger dynamics and gravitational wave emission in neutron star mergers. Phys. Rev. D, 91:124056, Jun 2015.
 [42] A. Bauswein, N. Stergioulas, and H.T. Janka. submitted to the EPJA topical issue ”Exotic Matter in Neutron Stars”, ArXiv eprints, 2015.
 [43] S. Bernuzzi, T. Dietrich, and A. Nagar. Modeling the complete gravitational wave spectrum of neutron stars mergers. ArXiv eprints, April 2015.
 [44] M. Favata. Systematic Parameter Errors in Inspiraling Neutron Star Binaries. Phys. Rev. Lett., 112(10):101101, March 2014.
 [45] J. Clark, A. Bauswein, L. Cadonati, H.T. Janka, C. Pankow, and N. Stergioulas. Prospects for high frequency burst searches following binary neutron star coalescence with advanced gravitational wave detectors. Phys. Rev. D, 90(6):062004, September 2014.
 [46] A. Bauswein, H.T. Janka, K. Hebeler, and A. Schwenk. Equationofstate dependence of the gravitationalwave signal from the ringdown phase of neutronstar mergers. Phys. Rev. D., 86(6):063001, September 2012.
 [47] Y. Sugahara and H. Toki. Relativistic meanfield theory for unstable nuclei with nonlinear and terms. Nuclear Physics A, 579:557–572, October 1994.
 [48] M. Hempel, T. Fischer, J. SchaffnerBielich, and M. Liebendörfer. New Equations of State in Simulations of Corecollapse Supernovae. Astrophys. J., 748:70, March 2012.
 [49] W. Kastaun and F. Galeazzi. Properties of hypermassive neutron stars formed in mergers of spinning binaries. Phys. Rev. D., 91(6):064027, March 2015.
 [50] K. Takami, L. Rezzolla, and L. Baiotti. Spectral properties of the postmerger gravitationalwave signal from binary neutron stars. Phys. Rev. D., 91(6):064001, March 2015.
 [51] K. Takami, L. Rezzolla, and L. Baiotti. Constraining the Equation of State of Neutron Stars from Binary Mergers. Physical Review Letters, 113(9):091104, August 2014.
 [52] M. Hannam, D. A. Brown, S. Fairhurst, C. L. Fryer, and I. W. Harry. When can Gravitationalwave Observations Distinguish between Black Holes and Neutron Stars? Astrophys. J. Lett., 766:L14, March 2013.
 [53] C. L. Rodriguez, B. Farr, V. Raymond, W. M. Farr, T. B. Littenberg, D. Fazi, and V. Kalogera. Basic Parameter Estimation of Binary Neutron Star Systems by the Advanced LIGO/Virgo Network. Astrophys. J., 784:119, April 2014.
 [54] B. Farr, C. P. L. Berry, W. M. Farr, C.J. Haster, H. Middleton, K. Cannon, P. B. Graff, C. Hanna, I. Mandel, C. Pankow, L. R. Price, T. Sidery, L. P. Singer, A. L. Urban, A. Vecchio, J. Veitch, and S. Vitale. Parameter estimation on gravitational waves from neutronstar binaries with spinning components. ArXiv eprints, August 2015.
 [55] J. M. Lattimer. The Nuclear Equation of State and Neutron Star Masses. Annu. Rev. Nucl. Part. Sci., 62:485–515, November 2012.
 [56] R. Oechslin, S. Rosswog, and F.K. Thielemann. Conformally flat smoothed particle hydrodynamics application to neutron star mergers. Phys. Rev. D., 65(10):103005–+, May 2002.
 [57] A. Bauswein, H.T. Janka, and R. Oechslin. Testing approximations of thermal effects in neutron star merger simulations. Phys. Rev. D., 82(8):084043–+, October 2010.
 [58] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels. A twosolarmass neutron star measured using Shapiro delay. Nature, 467:1081–1083, October 2010.
 [59] John Antoniadis, Paulo C. C. Freire, Norbert Wex, Thomas M. Tauris, Ryan S. Lynch, Marten H. van Kerkwijk, Michael Kramer, Cees Bassa, Vik S. Dhillon, Thomas Driebe, Jason W. T. Hessels, Victoria M. Kaspi, Vladislav I. Kondratiev, Norbert Langer, Thomas R. Marsh, Maura A. McLaughlin, Timothy T. Pennucci, Scott M. Ransom, Ingrid H. Stairs, Joeri van Leeuwen, Joris P. W. Verbiest, and David G. Whelan. A massive pulsar in a compact relativistic binary. Science, 340(6131), 2013.
 [60] G. A. Lalazissis, J. König, and P. Ring. New parametrization for the Lagrangian density of relativistic mean field theory. Phys. Rev. C., 55:540–543, January 1997.
 [61] M. Hempel and J. SchaffnerBielich. A statistical model for a complete supernova equation of state. Nucl. Phys. A, 837:210–254, June 2010.
 [62] J. M. Lattimer and F. D. Swesty. A generalized equation of state for hot, dense matter. Nucl. Phys. A, 535:331–376, December 1991.
 [63] S. Typel, G. Röpke, T. Klähn, D. Blaschke, and H. H. Wolter. Composition and thermodynamics of nuclear matter with light clusters. Phys. Rev. C., 81(1):015803, January 2010.
 [64] A. W. Steiner, M. Hempel, and T. Fischer. Corecollapse Supernova Equations of State Based on Neutron Star Observations. Astrophys. J., 774:17, September 2013.
 [65] G. Shen, C. J. Horowitz, and E. O’Connor. Second relativistic mean field and virial equation of state for astrophysical simulations. Phys. Rev. C., 83(6):065808, June 2011.
 [66] H. Toki, D. Hirata, Y. Sugahara, K. Sumiyoshi, and I. Tanihata. Relativistic many body approach for unstable nuclei and supernova. Nuclear Physics A, 588:357–363, February 1995.
 [67] A. Akmal, V. R. Pandharipande, and D. G. Ravenhall. Equation of state of nucleon matter and neutron star structure. Phys. Rev. C., 58:1804–1828, September 1998.
 [68] S. Banik, M. Hempel, and D. Bandyopadhyay. New Hyperon Equations of State for Supernovae and Neutron Stars in Densitydependent Hadron Field Theory. Astrophys. J. Supp., 214:22, October 2014.
 [69] H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi. Relativistic equation of state of nuclear matter for supernova and neutron star. Nucl. Phys. A, 637:435–450, July 1998.
 [70] K. Hebeler and A. Schwenk. Chiral threenucleon forces and neutron matter. Phys. Rev. C., 82(1):014314, July 2010.
 [71] C. Cutler and É. E. Flanagan. Gravitational waves from merging compact binaries: How accurately can one extract the binary’s parameters from the inspiral waveform? Phys. Rev. D., 49:2658–2697, March 1994.
 [72] J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, P. Ajith, B. Allen, and et al. TOPICAL REVIEW: Predictions for the rates of compact binary coalescences observable by groundbased gravitationalwave detectors. Classical and Quantum Gravity, 27(17):173001, September 2010.
 [73] LIGO Scientific Collaboration. Advanced LIGO anticipated sensitivity curves (LIGOT0900288). https://dcc.ligo.org/LIGOT0900288, 2015.
 [74] J. Miller, L. Barsotti, S. Vitale, P. Fritschel, M. Evans, and D. Sigg. Prospects for doubling the range of Advanced LIGO. Phys. Rev. D., 91(6):062005, March 2015.
 [75] LIGO Scientific Collaboration. Instrument Science White Paper (LIGOT1500290). https://dcc.ligo.org/LIGOT1500290, 2015.
 [76] M. Punturo, M. Abernathy, F. Acernese, B. Allen, N. Andersson, K. Arun, F. Barone, B. Barr, M. Barsuglia, M. Beker, N. Beveridge, S. Birindelli, S. Bose, L. Bosi, S. Braccini, C. Bradaschia, T. Bulik, E. Calloni, G. Cella, E. Chassande Mottin, S. Chelkowski, A. Chincarini, J. Clark, E. Coccia, C. Colacino, J. Colas, A. Cumming, L. Cunningham, E. Cuoco, S. Danilishin, K. Danzmann, G. De Luca, R. De Salvo, T. Dent, R. Derosa, L. Di Fiore, A. Di Virgilio, M. Doets, V. Fafone, P. Falferi, R. Flaminio, J. Franc, F. Frasconi, A. Freise, P. Fulda, J. Gair, G. Gemme, A. Gennai, A. Giazotto, K. Glampedakis, M. Granata, H. Grote, G. Guidi, G. Hammond, M. Hannam, J. Harms, D. Heinert, M. Hendry, I. Heng, E. Hennes, S. Hild, J. Hough, S. Husa, S. Huttner, G. Jones, F. Khalili, K. Kokeyama, K. Kokkotas, B. Krishnan, M. Lorenzini, H. Lück, E. Majorana, I. Mandel, V. Mandic, I. Martin, C. Michel, Y. Minenkov, N. Morgado, S. Mosca, B. Mours, H. MüllerEbhardt, P. Murray, R. Nawrodt, J. Nelson, R. Oshaughnessy, C. D. Ott, C. Palomba, A. Paoli, G. Parguez, A. Pasqualetti, R. Passaquieti, D. Passuello, L. Pinard, R. Poggiani, P. Popolizio, M. Prato, P. Puppo, D. Rabeling, P. Rapagnani, J. Read, T. Regimbau, H. Rehbein, S. Reid, L. Rezzolla, F. Ricci, F. Richard, A. Rocchi, S. Rowan, A. Rüdiger, B. Sassolas, B. Sathyaprakash, R. Schnabel, C. Schwarz, P. Seidel, A. Sintes, K. Somiya, F. Speirits, K. Strain, S. Strigin, P. Sutton, S. Tarabrin, J. van den Brand, C. van Leewen, M. van Veggel, C. van den Broeck, A. Vecchio, J. Veitch, F. Vetrano, A. Vicere, S. Vyatchanin, B. Willke, G. Woan, P. Wolfango, and K. Yamamoto. The third generation of gravitational wave observatories and their science reach. Classical and Quantum Gravity, 27(8):084007, April 2010.
 [77] S. Hild et al. Sensitivity studies for thirdgeneration gravitational wave observatories. Classical and Quantum Gravity, 28(9):094013, May 2011.
 [78] I. S. Heng. Rotating stellar corecollapse waveform decomposition: a principal component analysis approach. Classical and Quantum Gravity, 26(10):105005, May 2009.
 [79] C. Röver, M.A. Bizouard, N. Christensen, H. Dimmelmeier, I. S. Heng, and R. Meyer. Bayesian reconstruction of gravitational wave burst signals from simulations of rotating stellar core collapse and bounce. Phys. Rev. D., 80(10):102004, November 2009.
 [80] M. Pürrer. Frequencydomain reduced order models for gravitational waves from alignedspin compact binaries. Classical and Quantum Gravity, 31(19):195010, October 2014.
 [81] F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cournapeau, M. Brucher, M. Perrot, and E. Duchesnay. Scikitlearn: Machine learning in Python. Journal of Machine Learning Research, 12:2825–2830, 2011.