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.

Marginally outer trapped surfaces (MOTSs, or marginal surfaces in short) are routinely used in numerical simulations of black hole spacetimes. They are an invaluable tool for locating and characterizing black holes quasi-locally in real time while the simulation is ongoing. It is often believed that a MOTS can behave unpredictably under time evolution; an existing MOTS can disappear, and a new one can appear without any apparent reason. In this paper we show that in fact the behavior of a MOTS is perfectly predictable and its behavior is dictated by a single real parameter, the stability parameter, which can be monitored during the course of a numerical simulation. We demonstrate the utility of the stability parameter to fully understand the variety of marginal surfaces that can be present in binary black hole initial data. We also develop a new horizon finder capable of locating very highly distorted marginal surfaces and we show that even in these cases, the stability parameter perfectly predicts the existence and stability of marginal surfaces.

We find strong numerical evidence for a new phenomenon in a binary black hole spacetime, namely the merger of marginally outer trapped surfaces (MOTSs). By simulating the head-on collision of two nonspinning unequal mass black holes, we observe that the MOTS associated with the final black hole merges with the two initially disjoint surfaces corresponding to the two initial black holes. This yields a connected sequence of MOTSs interpolating between the initial and final state all the way through the non-linear binary black hole merger process. In addition, we show the existence of a MOTS with self-intersections formed immediately after the merger. This scenario now allows us to track physical quantities (such as mass, angular momentum, higher multipoles, and fluxes) across the merger, which can be potentially compared with the gravitational wave signal in the wave-zone, and with observations by gravitational wave detectors. This also suggests a possibility of proving the Penrose inequality mathematically for generic astrophysical binary back hole configurations.