Case E

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 435 M x 435 M, the innermost refinement level for the BH is 1.47 M x 1.47 M and the innermost refinement level for the NS is 6.8 M x 6.8 M. Here M is the total (ADM) mass of the initial system. The innermost resolution is ΔXmin/M = 1/58.8. In this simulation, the initial binary coordinate separation is D0/M = 8.61, the initial angular momentum of the system is J/M2 = 0.938, the mass ratio is MBH:MNS = 1: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.

Because of the low mass ratio in Case E, tidal disruption occurs at a farther binary separation. After tidal disruption, the NS matter curls around the BH forming a hot, low-density spiral that winds around the AH and smashes into the tidal tail, generating a large amount of shock heating. A fraction of the heated NS matter in the tail loses angular momentum and falls into the BH. The rest (2.3% of the rest mass) deforms into an inhomogeneous disk before settling into a quasistationary state. At the end of the simulation, the NS matter in the disk settles into a high density, low temperature torus of matter surrounding the BH.

Fig. 2-1 Color code for density profile

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

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

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

Download MPEG-4 (7.3 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 = 250

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

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

Download MPEG-4 (10.8 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.01%. Total angular momentum conservation is obeyed to within 0.9%.

Fig. 4-1 h+ Profile

Fig. 4-2 hx Profile

Download MPEG-4 (9.9 MB)

Play Online

Download MPEG-4 (10.3 MB)

Play Online

Fig. 4-3 h+ in orbital plane

Download MPEG-4 (26.7 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/J7.2%
&Delta EGW/M0.35%
Recoil velocity17 km/s

Download MPEG-4 (4.2 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