In EM scattering problems in inhomogeneous (layered) media one may encounter Sommerfeld integrals of the form:(adsbygoogle = window.adsbygoogle || []).push({});

[tex]\int_{0}^{\infty}J_{n}(k_{\rho},\rho)k_{\rho}^{n+1}G(k_{\rho})dk_{\rho}[/tex]

where J is a Bessel function and G is a spectral Green's function, [itex]\rho[/itex] is source-observer distance and in my case [itex]n=0[/itex].

Generally the Green's function results in a branch cut on the real axis and poles either on or slightly off the real axis.

As I see it the integrand contains the mode structure of the field and the integration over wavevectors (k) leads to something like a total 'self-energy' from which dispersion and photonic mode density can be obtained.

I am trying to obtain numerical results for an integral of this type in the context of the problem of decay of an excited molecule above a metal dielectric interface (allowing creation of surface plasmons, hence the pole).

If anyone could give me pointers on how to approach this specific integral or, more generally, how to deal with poles and branch cuts in numerical integration, I would greatly appreciate it.

Thanks,

AR

**Physics Forums - The Fusion of Science and Community**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Numerical evaluation of Sommerfeld Integral

Loading...

Similar Threads - Numerical evaluation Sommerfeld | Date |
---|---|

I Numerical Diagonalization | Mar 9, 2017 |

A Computing solutions to the radial Schroedinger equation? | Mar 9, 2017 |

A Efficiently Computing Eigenvalues of a Sparse Banded Matrix | Nov 17, 2016 |

A Numerical solution to SE - variational method, many electrons | May 7, 2016 |

I Evaluating Scattering Integral | Apr 1, 2016 |

**Physics Forums - The Fusion of Science and Community**