“…The first detection of gravitational waves (GWs), GW150914 from binary black holes (BBHs) [1,2,3,4] with the succeeding detection, GW151226 [5], GW170104 [6], GW170814 [7], and a candidate event, LVT151012 [8], recorded by Advanced LIGO detectors [9,10,11] opened a new window on physics and the Universe. To perform such GW astrophysics with very high precision in the context of ground-based GW detectors, including Advanced LIGO, Advanced Virgo [12] and KAGRA [13,14] as well as planned space-based GW detectors such as LISA [15] and (B-)DECIGO [16,17], it is now crucial to have extremely accurate predictions of GWs emitted from BBHs to maximize the extraction of physical information from noisy GW signals through the well-known technique of matched filtering; cross correlating the noisy detector output with the theoretical GW waveforms for the expected GW signal (see, e.g., [18] for the algorithm used by LIGO Scientific Collaboration).…”