Multimode frequency-domain model for the gravitational wave signal from nonprecessing black-hole binaries

We present the imrphenomxhm frequency domain phenomenological waveform model for the inspiral, merger, and ringdown of quasicircular nonprecessing black hole binaries. The model extends the imrphenomxas waveform model [G. Pratten, S. Husa, C. Garc\'{\i}a-Quir\'os, M. Colleoni, A. Ramos-Buades, H. Estell\'es, and R. Jaume, preceding paper, Phys. Rev. D 102, 064001 (2020)], which describes the dominant quadrupole modes $\ensuremath{\ell}=|m|=2$, to the harmonics $(\ensuremath{\ell},|m|)=(2,1),(3,3),(3,2),(4,4)$, and includes mode mixing effects for the (3,2) spherical harmonic. imrphenomxhm is calibrated against hybrid waveforms, which match an inspiral phase described by the effective-one-body model and post-Newtonian amplitudes for the subdominant harmonics to numerical relativity waveforms and numerical solutions to the perturbative Teukolsky equation for large mass ratios up to 1000.