Especially if one wishes to keep track of all the variables, like flux distribution vs. spatial, energy, wavevector, polarization, etc. for various geometries.
Nanophotonic scattering structures close to the emission zone make it possible to tap the “near field” of these “organic modes”. In the case of scattering the fields are the solution of an integral equation defined on the volume of the scatterers. Its kernel is given by the Green’s tensor. The numerical solution poses a number of implementation problems which can be overcome.
Finally had time to look into this question more seriously. I found the following resources quite valuable on addressing the issues of scattering events: