Gravitational Waveforms


The gravitational wavetrain from SMS collapse may be separated into two qualitatively different phases: the collapsing phase and the BH ringdown. During the early collapsing phase, the gravitational wave has smaller, growing intensity and lower frequency. As the matter source becomes denser, the wave intensity increases and right before the BH forms, the intensity reaches its peak. Once the BH forms, ringdown radiation is emitted as the distorted BH settles down to Kerr-like 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 $M_{\text{disk}} \sim 0.1 M_{\text{BH}} $. Since the collapse is nearly axisymmetric, the dominant amplitude is the h+ polarization, with nearly zero h× polarization. The plot for h+ is shown below.

h+ Polarization (Lower Hemisphere)

Shown below is h+ in the lower hemisphere. Superimposed is the star, which collapses to a BH immersed in a disk. Due to the scale of the movie and the dense disk, the BH shown here is hard to observe. The interior-only B-field case is depicted here, but all other cases are similar. We show only the waveform near the time that the BH forms. We observe at $r/M \geq 100$, and the scale indicates the amplitude at $r/M = 100$.
Fig. 1-1: t/M = 224
Fig. 1-1: t/M = 224
Fig. 1-2: t/M = 454
Fig. 1-2: t/M = 454
Fig. 1-3: t/M = 698
Fig. 1-3: t/M = 750
Fig. 1-4: t/M = 1224
Fig. 1-4: t/M = 1224
Fig. 1-5: t/M = 1553
Fig. 1-5: t/M = 2593
Fig. 1-6: t/M = 1823
Fig. 1-6: t/M = 2803
Play Online

Waveform Analysis

Figure 2-1 shows the real part of the quadrupole mode with $l = 2$, $m = 0$ of Weyl scalar $\Psi_{2,0}$ for the three cases in this project. No significant differences of waveform amplitudes and wavelength exist among threes cases. The lower panel displays the $h_{+}$ for the $l = 2$, $m =0 $ mode with observer sitting on the equator.

Fig. 1-6: t/M = 1801
Fig. 2-1: Top: Real part of Weyl scalar $\Psi_{2,0}$ for $3$ cases as a function of retarded time. Gravitational waves are extracted at r = $100.6$ M. The time has been shifted by the BH-formation time $t_{BH}$ for each case respectively.

For a SMBH with mass $\sim 10^6 M_{\odot}$, the frequency is located in the most sensitive band of eLISA for $z \lesssim 7$. To confirm detectability, we calculate the strain amplitude $|h(f)| = (h(f)_+^2 + h(f)_{\times}^2)^{1/2}$ and compare with the sensitivity curve of current GW detectors. Fig. 2-2 shows the Fourier spectrum of strain times freqeuncy $|h(f)|f$ versus frequency for three cases at redshift z = 1 and Case C at z = 2 and 3 alongside the characteristic strain amplitude, $S(f)f^{1/2}$, where $S(f)$ is the linear spectral density, from two different modes of LISA. Both modes have 4 laser links and arm length $L = 5 \times 10^6$ km with one having an optimistically functioning LISA Pathfinder(N2) and the other 10 times lower(N1). As expected from the Weyl scalar, waveforms of the three cases overlap closely at the same redshift. With $M = 10^{6} M_{\odot}$, the peak strain at z=1 is: $$ h \approx 9.2 \times 10^{-21} \left(\frac{M}{10^6 M_{\odot}}\right) \left( \frac{6.8 \text{Gpc}}{D_L}\right) $$ where $D_L = 6.8$ Gpc, corresponding to redshift z = 1 in a $\Lambda$CDM cosmology model with $H_0 = 67.6 \text{km s}^{-1} \text{Mpc}^{-1}$ and $\Omega_{M} = 0.311$. In addition, the plot shows that the peak strain is located near the most sensitive frequency range of LISA. For the N1 configuration, the SNR is about 5 at z =1. If the sensitivity of LISA is enhanced to that of LISA project(N2), the observation of GW signals for a similar object as far as z = 3 is possible with SNR $\gtrsim 10$. Due to the extremely low gas metalicity, a SMS is unlikely to form for $z \lesssim 2.5$. Hence, a GW detector as sensitive as LISA is required to detect such signals.

Fig. 1-6: t/M = 1801
Fig. 2-2: $2|h(f)|f$ vs. frequency diagram for collapses models. The upper (3 overlapped), middle and lower curves denote the signal strength from z = 1, 2, and 3, respectively. The bottom curve is the characteristic strain amplitudes of LISA.