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\'{}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 =|m|=2, to the harmonics (,|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.