Catching supermassive black hole binaries without a net

The gravitational wave signals from coalescing Supermassive Black Hole Binaries are prime targets for the Laser Interferometer Space Antenna (LISA). With optimal data processing techniques, the LISA observatory should be able to detect black hole mergers anywhere in the Universe. The challenge is to find ways to dig the signals out of a combination of instrument noise and the large foreground from stellar mass binaries in our own galaxy. The standard procedure of matched filtering against a grid of templates can be computationally prohibitive, especially when the black holes are spinning or the mass ratio is large. Here we develop an alternative approach based on Metropolis-Hastings sampling and simulated annealing that is orders of magnitude cheaper than a grid search. For the first time, we show that it is possible to detect and characterize the signals from binary systems of Schwarzschild Black Holes that are embedded in instrument noise and a foreground containing millions of galactic binaries. Our technique is computationally efficient, robust, and applicable to both high and low signal-to-noise ratio systems.