We present the first numerical-relativity simulation of a compact-object binary whose gravitational waveform is long enough to cover the entire frequency band of advanced gravitational-wave detectors, such as LIGO, Virgo, and KAGRA, for mass ratio 7 and total mass as low as 45.5M ⊙ . We find that effective-one-body models, either uncalibrated or calibrated against substantially shorter numericalrelativity waveforms at smaller mass ratios, reproduce our new waveform remarkably well, with a negligible loss in detection rate due to modeling error. In contrast, post-Newtonian inspiral waveforms and existing calibrated phenomenological inspiral-merger-ringdown waveforms display greater disagreement with our new simulation. The disagreement varies substantially depending on the specific post-Newtonian approximant used. [6]. The detection of GWs from compact-object binaries, as well as the determination of source properties from detected GW signals, relies on the accurate knowledge of the expected gravitational waveforms via matchedfiltering [7] and Bayesian methods [8].The need for accurate waveforms has motivated intense research. Early waveform models based on the postNewtonian (PN) formalism [9] were limited to the early inspiral. The effective-one-body (EOB) formalism [10,11] extended waveform models to the late inspiral, merger, and ringdown. Since 2005, research has greatly benefited from numerical-relativity (NR) simulations [12][13][14]. (Besides its importance for GW astronomy, NR has also deepened the understanding of general relativity in topics such as binary BH recoil [15,16], gravitational self-force [17], highenergy physics, and cosmology [18,19].) Current inspiral-merger-ringdown (IMR) waveform models [20][21][22][23]