Phenomenological gravitational-wave model for precessing black-hole binaries with higher multipoles and asymmetries

In this work we introduce phenomxo4a, the first phenomenological, frequency-domain gravitational waveform model to incorporate multipole asymmetries and precession angles tuned to numerical relativity. We build upon the modeling work that produced the phenompnr model and incorporate our additions into the imrphenomx framework, retuning the coprecessing frame model and extending the tuned precession angles to higher signal multipoles. We also include, for the first time in frequency-domain models, a recent model for spin-precession-induced multipolar asymmetry in the coprecessing frame to the dominant gravitational-wave multipoles. The accuracy of the full model and its constituent components is assessed through comparison to numerical relativity and numerical relativity surrogate waveforms by computing mismatches and performing parameter estimation studies. We show that, for the dominant signal multipole, we retain the modeling improvements seen in the phenompnr model. We find that the relative accuracy of current full IMR models varies depending on location in parameter space and the comparison metric, and on average they are of comparable accuracy. However, we find that variations in the pointwise accuracy do not necessarily translate into large biases in the parameter estimation recoveries.