Generating Higher Order Modes from Binary Black Hole mergers with Machine Learning

We introduce a machine learning model designed to rapidly and accurately predict the time domain gravitational wave emission of non-precessing binary black hole coalescences, incorporating the effects of higher order modes of the multipole expansion of the waveform. Expanding on our prior work, we decompose each mode by amplitude and phase and reduce dimensionality using principal component analysis. An ensemble of artificial neural networks is trained to learn the relationship between orbital parameters and the low-dimensional representation of each mode. We train our model on 10⁵ signals with mass ratio q 1, 10 and dimensionless spins ᵢ -0. 9, 0. 9, generated with the state-of-the-art approximant SEOBNRv4HM. We find that it achieves a median faithfulness of 10^-4 averaged across the parameter space. We show that our model generates a single waveform two orders of magnitude faster than the training model, with the speed up increasing when waveforms are generated in batches. This framework is entirely general and can be applied to any other time domain approximant capable of generating waveforms from aligned spin circular binaries, possibly incorporating higher order modes.