We propose an automatic algorithm for 3D inverse electromagnetic scattering based on the combination of topological derivatives and regularized Gauss-Newton iterations. The algorithm is adapted to decoding digital holograms. A hologram is a two-dimensional light interference pattern that encodes information about three-dimensional shapes and their optical properties. The formation of the hologram is modeled using Maxwell theory for light scattering by particles. We then seek shapes optimizing error functionals which measure the deviation from the recorded holograms. Their topological derivatives provide initial guesses of the objects. Next, we correct these predictions by regularized Gauss-Newton techniques. In contrast to standard Gauss-Newton methods, in our implementation the number of objects can be automatically updated during the iterative procedure by new topological derivative computations. We show that the combined use of topological derivative based optimization and iteratively regularized Gauss-Newton methods produces fast and accurate descriptions of the geometry of objects formed by multiple components with nanoscale resolution, even for a small number of detectors and non convex components aligned in the incidence direction. The method could be applied in general imaging set-ups involving other waves (microwave imaging, elastography...) provided closed-form expressions for the topological and Fréchet derivatives are determined.

Abstract. In this paper we study the shape differentiability properties of a class of boundary integral operators and of potentials with weakly singular pseudo-homogeneous kernels acting between classical Sobolev spaces, with respect to smooth deformations of the boundary. We prove that the boundary integral operators are infinitely differentiable without loss of regularity. The potential operators are infinitely shape differentiable away from the boundary, whereas their derivatives lose regularity near the boundary. We study the shape differentiability of surface differential operators. The shape differentiability properties of the usual strongly singular or hypersingular boundary integral operators of interest in acoustic, elastodynamic or electromagnetic potential theory can then be established by expressing them in terms of integral operators with weakly singular kernels and of surface differential operators.

Abstract. We construct and analyze a family of well-conditioned boundary integral equations for the Krylov iterative solution of three-dimensional elastic scattering problems by a bounded rigid obstacle. We develop a new potential theory using a rewriting of the Somigliana integral representation formula. From these results, we generalize to linear elasticity the well-known Brakhage-Werner and Combined Field Integral Equation formulations. We use a suitable approximation of the Dirichlet-to-Neumann (DtN) map as a regularizing operator in the proposed boundary integral equations. The construction of the approximate DtN map is inspired by the On-Surface Radiation Conditions method. We prove that the associated integral equations are uniquely solvable and possess very interesting spectral properties. Promising analytical and numerical investigations, in terms of spherical harmonics, with the elastic sphere are provided.

This paper is devoted to a modal formulation of the Linear Sampling Method in elastic 2D or 3D waveguide, that is we use the guided modes (called Lamb modes in 2D) as incident waves and the corresponding far fields in order to retrieve some obstacles. We provide the mathematical background to tackle the problem of identifiability and the justification of the Linear Sampling Method for such case. The elastic waveguide raises a specific issue: it concerns the projection of the scattered field on a transverse basis, which requires the introduction of new variables that mix displacement and stress components and satisfy the so-called Fraser's biorthogonality condition. Some numerical experiments in 2D show the feasibility of the reconstruction in the case of a finite number of incident waves formed by Lamb modes.