Self-intersecting marginally outer trapped surfaces

We have shown previously that a merger of marginally outer trapped surfaces (MOTSs) occurs in a binary black hole merger and that there is a continuous sequence of MOTSs which connects the initial two black holes to the final one. In this paper, we confirm this scenario numerically and we detail further improvements in the numerical methods for locating MOTSs. With these improvements, we confirm the merger scenario and demonstrate the existence of self-intersecting MOTSs formed in the immediate aftermath of the merger. These results will allow us to track physical quantities across the non-linear merger process and to potentially infer properties of the merger from gravitational wave observations

A trapped surface in the higher-dimensional self-similar Vaidya spacetime

We investigate a trapped surface and naked singularity in a $D$-dimensional Vaidya spacetime with a self-similar mass function. A trapped surface is defined as a closed spacelike $(D-2)$-surface which has negative both null expansions. There is no trapped surface in the Minkowski spacetime. However, in a four-dimensional self-similar Vaidya spacetime, Bengtsson and Senovilla considered non-spherical trapped surfaces and showed that a trapped surface can penetrate into a flat region, if and only if the mass function rises fast enough [I. Bengtsson and J. M. M. Senovilla, Phys. Rev. D \textbf{79}, 024027 (2009).]. We apply this result to a $D$-dimensional spacetime motivated by the context of large extra dimensions or TeV-scale gravity. In this paper, similarly to Bengtsson and Senovilla's study, we match four types of $(D-2)$-surfaces and show that a trapped surface extended into the flat region can be constructed in the $D$-dimensional Vaidya spacetime, if the increasing rate of the mass function is greater than 0.4628. Moreover, we show that the maximum radius of the trapped surface constructed here approaches the Schwarzschild-Tangherlini radius in the large $D$ limit. Also, we show that there is no naked singularity, if the spacetime has the trapped surface constructed here.Comment: 13 pages, 5 figure

The runaway instability in general relativistic accretion disks

When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability.Comment: 5 pages, 3 figures, minor corrections to match published version in MNRAS, +link to animatio

Accurate Evolutions of Orbiting Binary Black Holes

We present a detailed analysis of binary black hole evolutions in the last orbit and demonstrate consistent and convergent results for the trajectories of the individual bodies. The gauge choice can significantly affect the overall accuracy of the evolution. It is possible to reconcile certain gauge-dependent discrepancies by examining the convergence limit. We illustrate these results using an initial data set recently evolved by Brügmann et al. [Phys. Rev. Lett. 92, 211101 (2004)]. For our highest resolution and most accurate gauge, we estimate the duration of this data set's last orbit to be approximately 59MADM

Numerical relativity with characteristic evolution, using six angular patches

The characteristic approach to numerical relativity is a useful tool in evolving gravitational systems. In the past this has been implemented using two patches of stereographic angular coordinates. In other applications, a six-patch angular coordinate system has proved effective. Here we investigate the use of a six-patch system in characteristic numerical relativity, by comparing an existing two-patch implementation (using second-order finite differencing throughout) with a new six-patch implementation (using either second- or fourth-order finite differencing for the angular derivatives). We compare these different codes by monitoring the Einstein constraint equations, numerically evaluated independently from the evolution. We find that, compared to the (second-order) two-patch code at equivalent resolutions, the errors of the second-order six-patch code are smaller by a factor of about 2, and the errors of the fourth-order six-patch code are smaller by a factor of nearly 50.Comment: 12 pages, 5 figures, submitted to CQG (special NFNR issue

Dynamics and gravitational wave signature of collapsar formation

We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formatio

Computational Relativistic Astrophysics With Adaptive Mesh Refinement: Testbeds

We have carried out numerical simulations of strongly gravitating systems based on the Einstein equations coupled to the relativistic hydrodynamic equations using adaptive mesh refinement (AMR) techniques. We show AMR simulations of NS binary inspiral and coalescence carried out on a workstation having an accuracy equivalent to that of a $1025^3$ regular unigrid simulation, which is, to the best of our knowledge, larger than all previous simulations of similar NS systems on supercomputers. We believe the capability opens new possibilities in general relativistic simulations.Comment: 7 pages, 16 figure

Neutrino-driven Turbulent Convection and Standing Accretion Shock Instability in Three-Dimensional Core-Collapse Supernovae

We conduct a series of numerical experiments into the nature of three-dimensional (3D) hydrodynamics in the postbounce stalled-shock phase of core-collapse supernovae using 3D general-relativistic hydrodynamic simulations of a $27$-$M_\odot$ progenitor star with a neutrino leakage/heating scheme. We vary the strength of neutrino heating and find three cases of 3D dynamics: (1) neutrino-driven convection, (2) initially neutrino-driven convection and subsequent development of the standing accretion shock instability (SASI), (3) SASI dominated evolution. This confirms previous 3D results of Hanke et al. 2013, ApJ 770, 66 and Couch & Connor 2014, ApJ 785, 123. We carry out simulations with resolutions differing by up to a factor of $\sim$4 and demonstrate that low resolution is artificially favorable for explosion in the 3D convection-dominated case, since it decreases the efficiency of energy transport to small scales. Low resolution results in higher radial convective fluxes of energy and enthalpy, more fully buoyant mass, and stronger neutrino heating. In the SASI-dominated case, lower resolution damps SASI oscillations. In the convection-dominated case, a quasi-stationary angular kinetic energy spectrum $E(\ell)$ develops in the heating layer. Like other 3D studies, we find $E(\ell) \propto \ell^{-1}$ in the "inertial range," while theory and local simulations argue for $E(\ell) \propto \ell^{-5/3}$. We argue that current 3D simulations do not resolve the inertial range of turbulence and are affected by numerical viscosity up to the energy containing scale, creating a "bottleneck" that prevents an efficient turbulent cascade.Comment: 24 pages, 15 figures. Accepted for publication in The Astrophysical Journal. Added one figure and made minor modifications to text according to suggestions from the refere