Case D

Evolution of Density Profile
Evolution of Density Profile with Velocity Field
Evolution of Gravitational Radiation Profile
Final Black Hole Parameters


Fig. 1-1 Initial Configuration of Binary

Evolution is performed on an AMR (Adaptive Mesh Refinement) grid. The outermost grid is 393 M x 393 M, the innermost refinement level is 2.2 M x 2.2 M. Here M is the total (ADM) mass of the initial system. The innermost resolution is ΔXmin/M = 1/48.9. In this simulation, the initial binary coordinate separation is D0/M = 8.79, the initial angular momentum of the system is J/M2 = 0.518, the mass ratio is MBH:MNS = 5:1, and the black hole spin parameter is JBH/MBH2 = 0 .00. (Note: The open circle in the lower right-hand corner of the above figure is a clock.)

Evolution of the Density Profile

In the clip from the equatorial plane, the rest-mass density of the neutron star is plotted on a logarithmic scale normalized to the initial central density. The gravitational field is evolved via the BSSN scheme using "moving puncture" gauge conditions. The relativistic hydrodynamic equations are solved using a high-resolution shock-capturing (HRSC) method.

After 4.5 orbits the NS begins to shed its outer layers due to tidal disruption. The NS basically plunges into the BH, leaving < 1% of the NS rest mass as a disk.

Fig. 2-1 Color code for density profile

Fig. 2-2 Density Profile at t/M = 570

Fig. 2-3 Density Profile at t/M = 610

Fig. 2-4 Density Profile at t/M = 851

Download MPEG-4 (6.8 MB)

Play Online

Evolution of Density Profile with Velocity Field

Fig. 3-1 Density Profile at t = 0

Fig. 3-2 Density Profile at t = 570

Fig. 3-3 Density Profile t/M = 610

Fig. 3-4 Density Profile at t/M = 851

Download MPEG-4 (7.2 MB)

Play Online

Evolution of Gravitational Radiation Profile

The gravitational wavetrain from a compact binary system may be separated into three qualitatively different phases: the inspiral, merger, and 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. Finally, ringdown radiation is emitted as the distorted black hole settles down to Kerr-like equilibrium (Note: Only in the case of a vacuum spacetime does the spinning BH obey the Kerr solution. The BHs formed here are surrounded by gaseous disks with small, but nonnegligible, rest mass). Both polarization modes (h+ and hx) are shown. Total energy conservation is obeyed to within 0.02%. Total angular momentum conservation is obeyed to within 1.9%.

Fig. 4-1 h+ Profile

Fig. 4-2 hx Profile

Download MPEG-4 (10.8 MB)

Play Online

Download MPEG-4 (9.3 MB)

Play Online

Fig. 4-3 h+ in orbital plane

Download MPEG-4 (28.8 MB)

Play Online

Final Black Hole Parameters

Listed in the table below is the dimensionless spin of the Kerr-like black hole at the end of our simulation. Also shown are the radiated energy, angular momentum, and linear recoil velocity resulting from gravitational wave emission. Additional parameters for the disk are given in the film clip.

&Delta JGW/J19.2%
&Delta EGW/M0.98%
Recoil velocity73 km/s

Download MPEG-4 (3.8 MB)

Play Online

Back to Index

last updated 15 December 2014 by aakhan3

Center for Theoretical Astrophysics---University of Illinois at Urbana-Champaign

Home | Research | Activities | Faculty | Postdocs | Graduate | Undergraduate | Movies