Gravitational Waveforms
The gravitational wavetrain from an NSNS merger may be separated into four qualitatively different phases: the inspiral, merger, HMNS spindown and BH ringdown. During the inspiral phase, which takes up most of the binary's lifetime, gravity wave emission gradually reduces the binary separation. The merger phase of the gravitational wavetrain is characterized by tidal disruption of the neutron star, followed by the formation of a hypermassive neutron star (HMNS). As the HMNS spins down, it pulses and emits gravitational waves. After the HMNS undergoes delayed collapse to a spining BH, ringdown radiation is emitted as the distorted BH settles down to Kerrlike equilibrium (Note: Only in the case of a vacuum spacetime does the spinning BH obey the exact Kerr solution. The BHs formed here are surrounded by gaseous disks with small, but nonnegligible, rest mass). The h_{×} polarization mode of the gravitational wave is shown below.
Waveform Analysis
Figure 21 shows the waveform of this case. The plot can be divided in three
regions. Region I is from $t/M = 0$ to $t/M \approx 580$, where the system is two neutron stars
orbiting each other. Region II is from $t/M\approx 580$ to $t/M \approx 2800$, where the system is a
hypermassive neutron star. Note that going from region I to region II,
the amplitude of the gravitational wave is decreased, and is slowly diminishing in time. In region II, a weaker quadrupole wave is produced by the nonaxisymmetric spindown and pulsation of the HMNS.
Region III is from $t/M \approx 2800$ until the end. Here, the
amplitude goes to zero as the spinning BH rings downs.
Fig 21: A plot of $h_{\times}$ vs. $t/M$ of the canonical interior NSNS case
Figure 22 shows the FFT of Figure 21. Notice the peak at around $\Omega M = 0.07$.
This is the frequency of the gravitational wave in region I, just prior to tidal break up and it is indeed twice the
orbital frequency of the neutron stars.
Fig 22: Fast Fourier transform of figure 21. The frequency is in units of
$\Omega M$.
In order to see the low frequency mode at the postmerger stage (Region II), we can FFT the waveform at times $t/M > 650$, and figure 23 is this plot.
Fig 23: Fast Fourier transform of figure 21 for $t/M > 650$
Table 21 summarizes the computed frequencies in the NSNS system. $(\Omega M)_{GW}$ is the
peaktopeak frequency of the gravitational wave, and $(\Omega M)_{orb}$ is the orbital angular
velocity of the NSNS binary, and $(\Omega M)_{GW2}$ is the frequency of the secondary peaks in region II of the figure 21. $(\Omega M)_{GW}$ is calculated in both region I, just prior to tidal distruption, as well as in region II, just prior to HMNS delayed collapse. $(\Omega M)_{spin}$ is the spin angular velocity of nonaxisymmetric HMNS just prior to delayed collapse. $(\Omega M)_{orb}$ and $(\Omega M)_{spin}$ are calculated by visually determing the frequency in 3D movie.
Note that in region I, $(\Omega M)_{GW} = 2(\Omega M)_{orb}$, and in region II, $(\Omega M)_{GW} = 2(\Omega M)_{spin}$ as expected for quadrupole waves
NSNS 
$(\Omega M)_{GW} = 0.06597$ 
$(\Omega M)_{orb} = 0.03290$ 

HMNS 
$(\Omega M)_{GW} = 0.16969$ 
$(\Omega M)_{GW2} = 0.05839$ 
$(\Omega M)_{spin} = 0.08224$ 

Table 21: Table of frequency values for the interior canonical case.