  Since 2015 the advanced Laser Interferometer Gravitational-Wave Observatory
(aLIGO) has detected a large number of gravitational wave events, originating
from both binary neutron stars and binary black hole (BBH) mergers. In light of
these detections, we simulate the dynamics of ambient test particles in the
gravitational potential well of a BBH system close to its inspiral phase with
the goal of simulating the associated electromagnetic radiation and resulting
spectral energy distribution of such a BBH system. This could shed light on
possible detection ranges of electromagnetic counterparts to BBH mergers. The
potentials are numerically calculated using finite difference methods, under
the assumption of non-rotating black holes with the post-Newtonian
Paczynski-Wiita potential approximation in tandem with retarded time concepts
analogous to electrodynamics. We find that the frequencies of potential
electromagnetic radiation produced by these systems (possibly reaching Earth),
range between a few $\text{kHz}$ to a few $100 \text{MHz}$.
