Higher-order modes in range calculations
lal
provides access to waveform families such as IMRPhenomXHM, which can generate waveforms beyond the usual (2, 2) mode that we use for range estimation. In principle, higher-order modes will extend the heavy end of the detection band by a modest amount at low redshift (see image below reproduced from doi:10.3390/galaxies10040090, where dashed lines show this extension).
Generally the HOM amplitude varies with the inclination \iota
of the system, meaning that solving for the horizon now also requires optimizing over \iota
, and hence multiple calls to the waveform generation function. I am not sure if there is some trick to speed this up.