Gravitational waves from the coalescence of two black holes carry the signature of the strong field dynamics of binary black holes. In this work we have used numerical relativity simulations and post-Newtonian theory to investigate this dynamics. Post-Newtonian theory is a low-velocity expansion that assumes the companion bodies to be point-particles, while numerical relativity treats black holes as extended objects with horizons and fully captures their dynamics. There is a priori no reason for the waveforms computed using these disparate methods to agree with each other, especially at late times when the black holes move close to the speed of light. We find, remarkably, that the leading order amplitudes in post-Newtonian theory agree well with the full general relativity solution for a large set of spherical harmonic modes, even in the most dynamical part of the binary evolution, with only some modes showing distinctly different behavior than that found by numerical relativity simulations. In particular, modes with spherical harmonic indices l = m as well as l = 2, m = 1 are least modified from their dominant post-Newtonian behavior. Understanding the nature of these modes in terms of the post-Newtonian description will aid in formulating better models of the emitted waveforms in the strong field regime of the dynamics.

The planned sensitivity upgrades to the LIGO and Virgo facilities could uniquely identify host galaxies of dark sirens-compact binary coalescences without any electromagnetic counterparts-within a redshift of z = 0.1. This is aided by the higher order spherical harmonic modes present in the gravitational-wave signal, which also improve distance estimation. In conjunction, sensitivity upgrades and higher modes will facilitate an accurate, independent measurement of the host galaxy's redshift in addition to the luminosity distance from the gravitational wave observation to infer the Hubble-Lema\^itre constant H0 to better than a few percent in five years. A possible Voyager upgrade or third generation facilities would further solidify the role of dark sirens for precision cosmology in the future.