Electromagnetic scattering beyond the weak regime: Solving the problem of divergent Born perturbation series by Padé approximants

Electromagnetic scattering is the main phenomenon behind all optical measurement methods where one aims to retrieve the shape or physical properties of an unknown object by measuring how it scatters an incident optical field. Such an inverse problem is often approached by solving, several times, the corresponding direct scattering problem and trying to find the best estimate of the object which is compatible with a set of measurements. In the direct scattering problem, two regimes can be distinguished depending on the size of the object and the permittivity contrast: the weak-scattering regime and the strong-scattering regime. Generally, the presence of the scatterer alters the form of the incident field inside the scatterer. If that effect is neglected in the physical model, then one speaks of the so-called single-scattering regime or, more often, the Born approximation. The regime in which this approximation is valid is the weak-scattering regime. The corresponding inverse problem, that aims to retrieve the object from scattering data, becomes linear in this case. Linearizing the problem simplifies the method to solve it, but also introduces limitations to the maximum spatial resolution achievable in the reconstruction of the object. In the strong-scattering regime, multiple-scattering effects are not neglected and the inverse problem is treated in its full non-linear nature, which makes finding its solution a far more challenging task. Despite the existence of numerical methods, a powerful way to solve those direct problems would be to use a perturbation approach where the field is expressed as a series, known as the Born series. The advantage of a perturbation approach stems from the fact that each term of the series has a clear physical meaning and can unveil much more about the scattering process than a purely numerical approach can offer. Unfortunately, the series solution turns out to be strongly divergent in the strong-scattering regime, making it an unpractical approach for problems under these strong-scattering conditions. Thus, despite the fact that multiple scattering could, in principle, allow resolving sub-wavelength details of the unknown object, this possibility is in practice hampered by the divergent nature of the higher-order terms of the Born series. In this work, we show how to solve this problem by employing Padé approximants and how to treat electromagnetic problems well beyond the weak-scattering regime and provide an accurate eval