In earlier work concerning the sparsity of the Hawking flux, we found it necessary to re-examine what is known regarding the greybody factors of black holes, with a view to extending and expanding on some old results from the 1970s. Focusing specifically on Schwarzschild black holes, we have re-calculated and re-assessed the greybody factors using a path-ordered-exponential approach, a technique which has the virtue of providing a pedagogically useful semi-explicit formula for the relevant Bogoliubov coefficients. These path-ordered-exponentials, being based on a variant of the “transfer matrix” formalism, are closely related to so-called “product integrals”, leading to quite straightforward and direct numerical evaluation, while side-stepping any need for numerically solving the relevant ordinary differential equations. Furthermore, while considerable analytic information is already available regarding both the high-frequency and low-frequency asymptotics of these greybody factors, numerical approaches seem better adapted to finding suitable “global models” for these greybody factors in the intermediate frequency regime, where most of the Hawking flux is actually concentrated. Working in a more general context, these path-ordered-exponential techniques are also likely to be of interest for generic barrier-penetration problems.
Numerical Eulerian method for linearized gas dynamics in the high frequency regime
