Source author record

Habib Ammari

Habib Ammari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

80works
20topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

80 published item(s)

preprint2024arXiv

Scattering from time-modulated subwavelength resonators

We consider wave scattering from a system of highly contrasting resonators with time-modulated material parameters. In this setting, the wave equation reduces to a system of coupled Helmholtz equations that models the scattering problem. We consider the one-dimensional setting. In order to understand the energy of the system, we prove a novel higher-order discrete, capacitance matrix approximation of the subwavelength resonant quasifrequencies. Further, we perform numerical experiments to support and illustrate our analytical results and show how periodically time-dependent material parameters affect the scattered wave field.

preprint2023arXiv

Transmission properties of time-dependent one-dimensional metamaterials

We solve the wave equation with periodically time-modulated material parameters in a one-dimensional high-contrast resonator structure in the subwavelength regime exactly, for which we compute the subwavelength quasifrequencies numerically using Muller's method. We prove a formula in the form of an ODE using a capacitance matrix approximation. Comparison of the exact results with the approximations reveals that the method of capacitance matrix approximation is accurate and significantly more efficient. We prove various transmission properties in the aforementioned structure and illustrate them with numerical simulations. In particular, we investigate the effect of time-modulated material parameters on the formation of degenerate points, band gaps and k-gaps.

preprint2022arXiv

Dynamic super-resolution in particle tracking problems

Particle tracking in biological imaging is concerned with reconstructing the trajectories, locations, or velocities of the targeting particles. The standard approach of particle tracking consists of two steps: first reconstructing statically the source locations in each time step, and second applying tracking techniques to obtain the trajectories and velocities. In contrast, the dynamic reconstruction seeks to simultaneously recover the source locations and velocities from all frames, which enjoys certain advantages. In this paper, we provide a rigorous mathematical analysis for the resolution limit of reconstructing source number, locations, and velocities by general dynamical reconstruction in particle tracking problems, by which we demonstrate the possibility of achieving super-resolution for the dynamic reconstruction. We show that when the location-velocity pairs of the particles are separated beyond certain distances (the resolution limits), the number of particles and the location-velocity pair can be stably recovered. The resolution limits are related to the cut-off frequency of the imaging system, signal-to-noise ratio, and the sparsity of the source. By these estimates, we also derive a stability result for a sparsity-promoting dynamic reconstruction. In addition, we further show that the reconstruction of velocities has a better resolution limit which improves constantly as the particles moving. This result is derived by an observation that the inherent cut-off frequency for the velocity recovery can be viewed as the total observation time multiplies the cut-off frequency of the imaging system, which may lead to a better resolution limit as compared to the one for each diffraction-limited frame. It is anticipated that this observation can inspire new reconstruction algorithms that improve the resolution of particle tracking in practice.

preprint2022arXiv

Nearly optimal resolution estimate for the two-dimensional super-resolution and a new algorithm for direction of arrival estimation with uniform rectangular array

In this paper, we develop a new technique to obtain nearly optimal estimates of the computational resolution limits introduced in Appl. Comput. Harmon. Anal. 56 (2022) 402-446; IEEE Trans. Inf. Theory 67(7) (2021) 4812-4827; Inverse Probl. 37(10) (2021) 104001 for two-dimensional super-resolution problems. Our main contributions are fivefold: (i) Our work improves the resolution estimate for number detection and location recovery in two-dimensional super-resolution problems to nearly optimal; (ii) As a consequence, we derive a stability result for a sparsity-promoting algorithm in two-dimensional super-resolution problems (or Direction of Arrival problems (DOA)). The stability result exhibits the optimal performance of sparsity promoting in solving such problems; (iii) Our techniques pave the way for improving the estimate for resolution limits in higher-dimensional super-resolutions to nearly optimal; (iv) Inspired by these new techniques, we propose a new coordinate-combination-based model order detection algorithm for two-dimensional DOA estimation and theoretically demonstrate its optimal performance, and (v) we also propose a new coordinate-combination-based MUSIC algorithm for super-resolving sources in two-dimensional DOA estimation. It has excellent performance and enjoys many advantages compared to the conventional DOA algorithms. The coordinate-combination idea seems to be a promising way for multi-dimensional DOA estimation.

preprint2022arXiv

Topological phenomena in honeycomb Floquet metamaterials

We dedicate this paper to the topological analysis of subwavelength solutions in Floquet metamaterials. This work should be considered as a basis for further investigation on whether topological properties of the bulk materials are linked to the occurrence of edge modes. The subwavelength solutions being described by a periodically parameterized time-periodic linear ordinary differential equation $\left\{\frac{d}{dt}X = A_α(t)X\right\}_{α\in \mathbb{T}^d}$, we put ourselves in the general setting of periodically parameterized time-periodic linear ordinary differential equations and introduce a way to (topologically) classify a Floquet normal form $F,~P$ of the associated fundamental solution $\left\{X_α(t) = P(α,t)\exp(tF_α)\right\}_{α\in \mathbb{T}^d}$. This is achieved by analysing the topological properties of the eigenvalues and eigenvectors of the monodromy matrix $X_α(T)$ and the Lyapunov transformation $P(α,t)$. The corresponding topological invariants can then be applied to the setting of Floquet metamaterials. In this paper these general results are considered in the case of honeycomb Floquet metamaterials. We provide two interesting examples of topologically non-trivial time-modulated honeycomb structures.

preprint2022arXiv

Transmission properties of space-time modulated metamaterials

We prove the possibility of achieving exponentially growing wave propagation in space-time modulated media and give an asymptotic analysis of the quasifrequencies in terms of the amplitude of the time modulation at the degenerate points of the folded band structure. Our analysis provides the first proof of existence of k-gaps in the band structures of space-time modulated systems of subwavelength resonators.

preprint2022arXiv

Unidirectional edge modes in time-modulated metamaterials

We prove the possibility of achieving unidirectional edge modes in time-modulated supercell structures. Such finite structures consist of two trimers repeated periodically. Because of their symmetry, they admit degenerate edge eigenspaces. When the trimers are time-modulated with two opposite orientations, the degenerate eigenspace splits into two one-dimensional eigenspaces described by an analytical formula, each corresponds to a mode which is localized at one edge of the structure. Our results on the localization and stability of these edge modes with respect to fluctuations in the time-modulation amplitude are illustrated by several numerical simulations.

preprint2020arXiv

A mathematical and numerical framework for gradient meta-surfaces built upon periodically repeating arrays of Helmholtz resonators

In this paper a mathematical model is given for the scattering of an incident wave from a surface covered with microscopic small Helmholtz resonators, which are cavities with small openings. More precisely, the surface is built upon a finite number of Helmholtz resonators in a unit cell and that unit cell is repeated periodically. To solve the scattering problem, the mathematical framework elaborated in [Ammari et al., Asympt. Anal., 2019] is used. The main result is an approximate formula for the scattered wave in terms of the lengths of the openings. Our framework provides analytic expressions for the scattering wave vector and angle and the phase-shift. It justifies the apparent absorption. Moreover, it shows that at specific lengths for the openings and a specific frequency there is an abrupt shift of the phase of the scattered wave due to the subwavelength resonances of the Helmholtz resonators. A numerically fast implementation is given to identify a region of those specific values of the openings and the frequencies.

preprint2020arXiv

Close-to-touching acoustic subwavelength resonators: eigenfrequency separation and gradient blow-up

In this paper, we study the behaviour of the coupled subwavelength resonant modes when two high-contrast acoustic resonators are brought close together. We consider the case of spherical resonators and use bispherical coordinates to derive explicit representations for the capacitance coefficients which, we show, capture the system's resonant behaviour at leading order. We prove that the pair of resonators has two subwavelength resonant modes whose frequencies have different leading-order asymptotic behaviour. We, also, derive estimates for the rate at which the gradient of the scattered pressure wave blows up as the resonators are brought together.

preprint2020arXiv

Edge modes in active systems of subwavelength resonators

Wave scattering structures with amplification and dissipation can be modelled by non-Hermitian systems, opening new ways to control waves at small length scales. In this work, we study the phenomenon of topologically protected edge states in acoustic systems with gain and loss. We demonstrate that localized edge modes appear in a periodic structure of subwavelength resonators with a defect in the gain/loss distribution, and explicitly compute the corresponding frequency and decay length. Similarly to the Hermitian case, these edge modes can be attributed to the winding of the eigenmodes. In the non-Hermitian case, the topological invariants fail to be quantized, but can nevertheless predict the existence of localized edge modes.

preprint2020arXiv

High-order exceptional points and enhanced sensing in subwavelength resonator arrays

Systems exhibiting degeneracies known as exceptional points have remarkable properties with powerful applications, particularly in sensor design. These degeneracies are formed when eigenstates coincide, and the remarkable effects are exaggerated by increasing the order of the exceptional point (that is, the number of coinciding eigenstates). In this work, we use asymptotic techniques to study PT-symmetric arrays of many subwavelength resonators and search for high-order asymptotic exceptional points. This analysis reveals the range of different configurations that can give rise to high-order exceptional points and provides efficient techniques to compute them. We also show how systems exhibiting high-order exceptional points can be used for sensitivity enhancement.

preprint2020arXiv

Honeycomb-lattice Minnaert bubbles

The ability to manipulate the propagation of waves on subwavelength scales is important for many different physical applications. In this paper, we consider a honeycomb lattice of subwavelength resonators and prove, for the first time, the existence of a Dirac dispersion cone at subwavelength scales. As shown in [H. Ammari et al., A high-frequency homogenization approach near the Dirac points in bubbly honeycomb crystals, arXiv:1812.06178], near the Dirac points, honeycomb crystals of subwavelength resonators has a great potential to be used as near-zero materials. Here, we perform the analysis for the example of bubbly crystals, which is a classic example of subwavelength resonance, where the resonant frequency of a single bubble is known as the Minnaert resonant frequency. Our first result is an asymptotic formula for the quasi-periodic Minnaert resonant frequencies. We then prove the linear dispersion relation of a Dirac cone. Our findings in this paper are numerically illustrated in the case of circular bubbles, where the multipole expansion method provides an efficient technique for computing the band structure.

preprint2020arXiv

Mathematical modelling of plasmonic strain sensors

We provide a mathematical analysis for a metasurface constructed of plasmonic nanoparticles mounted periodically on the surface of a microcapsule. We derive an effective transmission condition, which exhibits resonances depending on the inter-particle distance. When the microcapsule is deformed, the resonances are shifted. We fully characterise the dependence of these resonances on the deformation of the microcapsule, enabling the detection of strains at the microscale level. We present numerical simulations to validate our results.

preprint2020arXiv

Quantitative estimates for stress concentration of the Stokes flow between adjacent circular cylinders

When two inclusions with high contrast material properties are located close to each other in a homogeneous medium, stress may become arbitrarily large in the narrow region between them. In this paper, we investigate such stress concentration in the two-dimensional Stokes flow when inclusions are the two-dimensional cross sections of circular cylinders of the same radii and the background velocity field is linear. We construct two vector-valued functions which completely capture the singular behavior of the stress and derive an asymptotic representation formula for the stress in terms of these functions as the distance between the two cylinders tends to zero. We then show, using the representation formula, that the stress always blows up by proving that either the pressure or the shear stress component of the stress tensor blows up. The blow-up rate is shown to be $δ^{-1/2}$, where $δ$ is the distance between the two cylinders. To our best knowledge, this work is the first to rigorously derive the asymptotic solution in the narrow region for the Stokes flow.

preprint2020arXiv

Quantum ergodicity and localization of plasmon resonances

We are concerned with the geometric properties of the surface plasmon resonance (SPR). SPR is a non-radiative electromagnetic surface wave that propagates in a direction parallel to the negative permittivity/dielectric material interface. It is known that the SPR oscillation is very sensitive to the material interface. However, we show that the SPR oscillation asymptotically localizes at places with high magnitude of curvature in a certain sense. Our work leverages the Heisenberg picture of quantization and quantum ergodicity first derived by Shnirelman, Zelditch, Colin de Verdière and Helffer-Martinez-Robert, as well as certain novel and more general ergodic properties of the Neumann-Poincaré operator to analyze the SPR field, which are of independent interest to the spectral theory and the potential theory.

preprint2020arXiv

Super-resolution in recovering embedded electromagnetic sources in high contrast media

The purpose of this work is to provide a rigorous mathematical analysis of the expected super-resolution phenomenon in the time-reversal imaging of electromagnetic (EM) radiating sources embedded in a high contrast medium. It is known that the resolution limit is essentially determined by the sharpness of the imaginary part of the EM Green's tensor for the associated background. We first establish the close connection between the resolution and the material parameters and the resolvent of the electric integral operator, via the Lippmann-Schwinger representation formula. We then present an insightful characterization of the spectral structure of the integral operator for a general bounded domain and derive the pole-pencil decomposition of its resolvent in the high contrast regime. For the special case of a spherical domain, we provide some quantitative asymptotic behavior of the eigenvalues and eigenfunctions. These mathematical findings shall enable us to provide a concise and rigorous illustration of the super-resolution in the EM source reconstruction in high contrast media. Some numerical examples are also presented to verify our main theoretical results.

preprint2020arXiv

Topologically protected edge modes in one-dimensional chains of subwavelength resonators

The goal of this paper is to advance the development of wave-guiding subwavelength crystals by developing designs whose properties are stable with respect to imperfections in their construction. In particular, we make use of a locally resonant subwavelength structure, composed of a chain of high-contrast resonators, to trap waves at deep subwavelength scales. We first study an infinite chain of subwavelength resonator dimers and define topological quantities that capture the structure's wave transmission properties. Using this for guidance, we design a finite crystal that is shown to have wave localization properties, at subwavelength scales, that are robust with respect to random imperfections.

preprint2019arXiv

Optimization of Steklov-Neumann eigenvalues

This paper examines the Laplace equation with mixed boundary conditions, the Neumann and Steklov boundary conditions. This models a container with holes in it, like a pond filled with water but partly covered by immovable pieces on the surface. The main objective is to determine the right extent of the covering pieces, so that any shock inside the container yields a resonance. To this end, an algorithm is developed which uses asymptotic formulas concerning perturbations of the partitioning of the boundary pieces. Proofs for these formulas are established. Furthermore, this paper displays some results concerning bounds and examples with regards to the governing problem.

preprint2019arXiv

Subwavelength localized modes for acoustic waves in bubbly crystals with a defect

The ability to control wave propagation is of fundamental interest in many areas of physics. Photonic and phononic crystals have proved very useful for this purpose but, because they are based on Bragg interference, these artificial media require structures with large dimensions. In [Ammari et al., Subwavelength phononic bandgap opening in bubbly media, J. Diff. Eq., 263 (2017), 5610--5629], it has been proved that a subwavelength bandgap opening occurs in bubbly phononic crystals. To demonstrate the opening of a subwavelength phononic bandgap, a periodic arrangement of bubbles is considered and their subwavelength Minnaert resonance is exploited. In this paper, this subwavelength bandgap is used to demonstrate cavities, very similar to those obtained in photonic and phononic crystals, albeit of deeply subwavelength dimensions. The key idea is to perturb the size of a single bubble inside the crystal, thus creating a defect. The goal is then to analytically and numerically show that this crystal has a localized eigenmode close to the defect bubble.

preprint2019arXiv

The equivalent media generated by bubbles of high contrasts: Volumetric metamaterials and metasurfaces

We deal with the point-interaction approximations for the acoustic wave fields generated by a cluster of highly contrasted bubbles for a wide range of densities and bulk moduli contrasts. We derive the equivalent fields when the cluster of bubbles is appropriately distributed (but not necessarily periodically) in a bounded domain $Ω$ of $\mathbb{R}^3$. We handle two situations. 1. In the first one, we distribute the bubbles to occupy a $3$ dimensional domain. For this case, we show that the equivalent speed of propagation changes sign when the medium is excited with frequencies smaller or larger than (but not necessarily close to) the Minnaert resonance. As a consequence, this medium behaves as a reflective or absorbing depending on whether the used frequency is smaller or larger than this resonance. In addition, if the used frequency is extremely close to this resonance, for a cluster of bubbles with density above a certain threshold, then the medium behaves as a 'wall', i.e. allowing no incident sound to penetrate. 2. In the second one, we distribute the bubbles to occupy a $2$ dimensional (open or closed) surface, not necessarily flat. For this case, we show that the equivalent medium is modeled by a Dirac potential supported on that surface. The sign of the surface potential changes for frequencies smaller or larger than the Minnaert resonance, i.e. it behaves as a smart metasurface reducing or amplifying the transmitted sound across it. As in the $3$D case, if the used frequency is extremely close to this resonance, for a cluster of bubbles with density above an appropriate threshold, then the surface allows no incident sound to be transmitted across the surface, i.e. it behaves as a white screen.

preprint2018arXiv

Subwavelength resonances of encapsulated bubbles

The aim of this paper is to derive a formula for the subwavelength resonance frequency of an encapsulated bubble with arbitrary shape in two dimensions. Using Gohberg-Sigal theory, we derive an asymptotic formula for this resonance frequency, as a perturbation away from the resonance of the uncoated bubble, in terms of the thickness of the coating. The formula is numerically verified in the case of circular bubbles, where the resonance can be efficiently computed using the multipole method.

preprint2016arXiv

A mathematical and numerical framework for bubble meta-screens

The aim of this paper is to provide a mathematical and numerical framework for the analysis and design of bubble meta-screens. An acoustic meta-screen is a thin sheet with patterned subwavelength structures, which nevertheless has a macroscopic effect on the acoustic wave propagation. In this paper, periodic subwavelength bubbles mounted on a reflective surface (with Dirichlet boundary condition) is considered. It is shown that the structure behaves as an equivalent surface with Neumann boundary condition at the Minnaert resonant frequency which corresponds to a wavelength much greater than the size of the bubbles. Analytical formula for this resonance is derived. Numerical simulations confirm its accuracy and show how it depends on the ratio between the periodicity of the lattice, the size of the bubble, and the distance from the reflective surface. The results of this paper formally explain the super-absorption behavior observed in [V. Leroy et al., Phys. Rev. B, 2015].

preprint2016arXiv

A signal separation technique for sub-cellular imaging using dynamic optical coherence tomography

This paper aims at imaging the dynamics of metabolic activity of cells. Using dynamic optical coherence tomography, we introduce a new multi-particle dynamical model to simulate the movements of the collagen and the cell metabolic activity and develop an efficient signal separation technique for sub-cellular imaging. We perform a singular-value decomposition of the dynamic optical images to isolate the intensity of the metabolic activity. We prove that the largest eigenvalue of the associated Casorati matrix corresponds to the collagen. We present several numerical simulations to illustrate and validate our approach.

preprint2016arXiv

Effective medium theory for acoustic waves in bubbly fluids near Minnaert resonant frequency

We derive an effective medium theory for acoustic wave propagation in bubbly fluid near Minnaert resonant frequency. We start with a multiple scattering formulation of the scattering problem of an incident wave by a large number of identical small bubbles in a homogeneous fluid. Under certain conditions on the configuration of the bubbles, we justify the point interaction approximation and establish an effective medium theory for the bubbly fluid as the number of bubbles tends to infinity. The convergence rate is also derived. As a consequence, we show that near and below the Minnaert resonant frequency, the obtained effective media can have a high refractive index, which is the reason for the super-focusing experiment observed in "Subwavelength focusing in bubbly media using broadband time reversal", Physical Review, B (2015), by M. Lanoy et al. Moreover, our results indicate that the obtained effective medium can be dissipative above the Minnaert resonant frequency, while at that frequency, effective medium theory does not hold. Our theory sheds light on the mechanism of the extraordinary wave properties of metamaterials, which include bubbly fluid as an example, near resonant frequencies.

preprint2016arXiv

Identification of an inclusion in Multifrequency Electric Impedance Tomography

The multifrequency electrical impedance tomography is considered in order to image a conductivity inclusion inside a homogeneous background medium by injecting one current. An original spectral decomposition of the solution of the forward conductivity problem is used to retrieve the Cauchy data corresponding to the extreme case of perfect conductor. Using results based on the unique continuation we then prove the uniqueness of multifrequency electrical impedance tomography and obtain rigorous stability estimates. Our results in this paper are quite surprising in inverse conductivity problem since in general infinitely many input currents are needed in order to obtain the uniqueness in the determination of the conductivity.

preprint2016arXiv

Mathematical framework for abdominal electrical impedance tomography to assess fatness

This paper presents a static electrical impedance tomography (EIT) technique that evaluates abdominal obesity by estimating the thickness of subcutaneous fat. EIT has a fundamental drawback for absolute admittivity imaging because of its lack of reference data for handling the forward modeling errors. To reduce the effect of boundary geometry errors in imaging abdominal fat, we develop a depth-based reconstruction method that uses a specially chosen current pattern to construct reference-like data, which are then used to identify the border between subcutaneous fat and muscle. The performance of the proposed method is demonstrated by numerical simulations using 32-channel EIT system and human like domain.

preprint2016arXiv

Minnaert resonances for acoustic waves in bubbly media

Through the application of layer potential techniques and Gohberg-Sigal theory we derive an original formula for the Minnaert resonance frequencies of arbitrarily shaped bubbles. We also provide a mathematical justification for the monopole approximation of scattering of acoustic waves by bubbles at their Minnaert resonant frequency. Our results are complemented by several numerical examples which serve to validate our formula in two dimensions.

preprint2016arXiv

Shape reconstruction of nanoparticles from their associated plasmonic resonances

We prove by means of a couple of examples that plasmonic resonances can be used on one hand to classify shapes of nanoparticles with real algebraic boundaries and on the other hand to reconstruct the separation distance between two nanoparticles from measurements of their first collective plasmonic resonances. To this end, we explicitly compute the spectral decompositions of the Neumann-Poincaré operators associated with a class of quadrature domains and two nearly touching disks. Numerical results are included in support of our main findings.

preprint2016arXiv

Theory of plasmonic metasurfaces

In this paper we derive an impedance boundary condition to approximate the optical scattering effect of an array of plasmonic nanoparticles mounted on a perfectly conducting plate. We show that at some resonant frequencies the impedance blows up, allowing for a significant reduction of the scattering from the plate. Using the spectral properties of a Neumann-Poincare type operator, we investigate the dependency of the impedance with respect to changes in the nanoparticle geometry and configuration.

preprint2016arXiv

Towards monitoring critical microscopic parameters for electropermeabilization

Electropermeabilization is a clinical technique in cancer treatment to locally stimulate the cell metabolism. It is based on electrical fields that change the properties of the cell membrane. With that, cancer treatment can reach the cell more easily. Electropermeabilization occurs only with accurate dosage of the electrical field. For applications, a monitoring for the amount of electropermeabilization is needed. It is a first step to image the macroscopic electrical field during the process. Nevertheless, this is not complete, because electropermeabilization depends on critical individual properties of the cells such as their curvature. From the macroscopic field, one cannot directly infer that microscopic state. In this article, we study effective parameters in a homogenization model as the next step to monitor the microscopic properties in clinical practice. We start from a physiological cell model for electropermeabilization and analyze its well-posedness. For a dynamical homogenization scheme, we prove convergence and then analyze the effective parameters, which can be found by macroscopic imaging methods. We demonstrate numerically the sensitivity of these effective parameters to critical microscopic parameters governing electropermeabilization. This opens the door to solve the inverse problem of rreconstructing these parameters.

preprint2015arXiv

A mathematical and numerical framework for magnetoacoustic tomography with magnetic induction

We provide a mathematical analysis and a numerical framework for magnetoacoustic tomography with magnetic induction. The imaging problem is to reconstruct the conductivity distribution of biological tissue from measurements of the Lorentz force induced tissue vibration. We begin with reconstructing from the acoustic measurements the divergence of the Lorentz force, which is acting as the source term in the acoustic wave equation. Then we recover the electric current density from the divergence of the Lorentz force. To solve the nonlinear inverse conductivity problem, we introduce an optimal control method for reconstructing the conductivity from the electric current density. We prove its convergence and stability. We also present a point fixed approach and prove its convergence to the true solution. A new direct reconstruction scheme involving a partial differential equation is then proposed based on viscosity-type regularization to a transport equation satisfied by the electric current density field. We prove that solving such an equation yields the true conductivity distribution as the regularization parameter approaches zero. Finally, we test the three schemes numerically in the presence of measurement noise, quantify their stability and resolution, and compare their performance.

preprint2015arXiv

Mathematical analysis of plasmonic nanoparticles: the scalar case

Localized surface plasmons are charge density oscillations confined to metallic nanoparticles. Excitation of localized surface plasmons by an electromagnetic field at an incident wavelength where resonance occurs results in a strong light scattering and an enhancement of the local electromagnetic fields. This paper is devoted to the mathematical modeling of plasmonic nanoparticles. Its aim is threefold: (i) to mathematically define the notion of plasmonic resonance and to analyze the shift and broadening of the plasmon resonance with changes in size and shape of the nanoparticles; (ii) to study the scattering and absorption enhancements by plasmon resonant nanoparticles and express them in terms of the polarization tensor of the nanoparticle. Optimal bounds on the enhancement factors are also derived; (iii) to show, by analyzing the imaginary part of the Green function, that one can achieve super-resolution and super-focusing using plasmonic nanoparticles. For simplicity, the Helmholtz equation is used to model electromagnetic wave propagation.

preprint2015arXiv

Mathematical analysis of plasmonic resonances for nanoparticles: the full Maxwell equations

In this paper we use the full Maxwell equations for light propagation in order to analyze plasmonic resonances for nanoparticles. We mathematically define the notion of plasmonic resonance and analyze its shift and broadening with respect to changes in size, shape, and arrangement of the nanoparticles, using the layer potential techniques associated with the full Maxwell equations. We present an effective medium theory for resonant plasmonic systems and derive a condition on the volume fraction under which the Maxwell-Garnett theory is valid at plasmonic resonances.

preprint2015arXiv

Mathematical framework for multi-frequency identification of thin insulating and small conductive inhomogeneities

We are aiming to identify the thin insulating inhomogeneities and small conductive inhomogeneities inside an electrically conducting medium by using multi-frequency electrical impedance tomography (mfEIT). The thin insulating inhomogeneities are considered in the form of tubular neighborhood of a curve and small conductive inhomogeneities are regarded as circular disks. Taking advantage of the frequency dependent behavior of insulating objects, we give a rigorous derivation of the potential along thin insulating objects at various frequencies. Asymptotic formula is given to analyze relationship between inhomogeneities and boundary potential at different frequencies. In numerical simulations, spectroscopic images are provided to visualize the reconstructed admittivity at various frequencies. For the view of both kinds of inhomogeneities, an integrated reconstructed image based on principle component analysis (PCA) is provided. Phantom experiments are performed by using Swisstom EIT-Pioneer Set.

preprint2015arXiv

Phased and phaseless domain reconstruction in inverse scattering problem via scattering coefficients

In this work we shall review the (phased) inverse scattering problem and then pursue the phaseless reconstruction from far-field data with the help of the concept of scattering coefficients. We perform sensitivity, resolution and stability analysis of both phased and phaseless problems and compare the degree of ill-posedness of the phased and phaseless reconstructions. The phaseless reconstruction is highly nonlinear and much more severely ill-posed. Algorithms are provided to solve both the phased and phaseless reconstructions in the linearized case. Stability is studied by estimating the condition number of the inversion process for both the phased and phaseless cases. An optimal strategy is suggested to attain the infimum of the condition numbers of the phaseless reconstruction, which may provide an important guidance for efficient phaseless measurements in practical applications. To the best of our knowledge, the stability analysis in terms of condition numbers are new for the phased and phaseless inverse scattering problems, and are very important to help us understand the degree of ill-posedness of these inverse problems. Numerical experiments are provided to illustrate the theoretical asymptotic behavior, as well as the effectiveness and robustness of the phaseless reconstruction algorithm.

preprint2015arXiv

Surface Plasmon Resonance of Nanoparticles and Applications in Imaging

In this paper we provide a mathematical framework for localized plasmon resonance of nanoparticles. Using layer potential techniques associated with the full Maxwell equations, we derive small-volume expansions for the electromagnetic fields, which are uniformly valid with respect to the nanoparticle's bulk electron relaxation rate. Then, we discuss the scattering and absorption enhancements by plasmon resonant nanoparticles. We study both the cases of a single and multiple nanoparticles. We present numerical simulations of the localized surface plasmonic resonances associated to multiple particles in terms of their separation distance.

preprint2015arXiv

Well-posedness of a Pulsed Electric Field Model in Biological Media and its Finite Element Approximation

This work aims at providing a mathematical and numerical framework for the analysis on the effects of pulsed electric fields on biological media. Biological tissues and cell suspensions are described as having a heteregeneous permittivity and a heteregeneous conductivity. Well-posedness of the model problem and the regularity of its solution are established. A fully discrete finite element scheme is proposed for the numerical approximation of the potential distribution as a function of time and space simultaneously for an arbitrary shaped pulse, and it is demonstrated to enjoy the optimal convergence order in both space and time. The proposed numerical scheme has potential applications in the fields of medicine, food sciences, and biotechnology.

preprint2014arXiv

A mathematical and numerical framework for ultrasonically-induced Lorentz force electrical impedance tomography

We provide a mathematical analysis and a numerical framework for Lorentz force electrical conductivity imaging. Ultrasonic vibration of a tissue in the presence of a static magnetic field induces an electrical current by the Lorentz force. This current can be detected by electrodes placed around the tissue; it is proportional to the velocity of the ultrasonic pulse, but depends nonlinearly on the conductivity distribution. The imaging problem is to reconstruct the conductivity distribution from measurements of the induced current. To solve this nonlinear inverse problem, we first make use of a virtual potential to relate explicitly the current measurements to the conductivity distribution and the velocity of the ultrasonic pulse. Then, by applying a Wiener filter to the measured data, we reduce the problem to imaging the conductivity from an internal electric current density. We first introduce an optimal control method for solving such a problem. A new direct reconstruction scheme involving a partial differential equation is then proposed based on viscosity-type regularization to a transport equation satisfied by the current density field. We prove that solving such an equation yields the true conductivity distribution as the regularization parameter approaches zero. We also test both schemes numerically in the presence of measurement noise, quantify their stability and resolution, and compare their performance.

preprint2014arXiv

A mathematical theory of super-resolution by using a system of sub-wavelength Helmholtz resonators

A rigorous mathematical theory is developed to explain the super-resolution phenomenon observed in the experiment by F.Lemoult, M.Fink and G.Lerosey (Acoustic resonators for far-field control of sound on a subwavelength scale, Phys. Rev. Lett., 107 (2011)). A key ingredient is the calculation of the resonances and the Green function in the half space with the presence of a system of Helmholtz resonators in the quasi-stationary regime. By using boundary integral equations and generalized Rouche's theorem, the existence and the leading asymptotic of the resonances are rigorously derived. The integral equation formulation also yields the leading order terms in the asymptotics of the Green function. The methodology developed in the paper provides an elegant and systematic way for calculating resonant frequencies for Helmholtz resonators in assorted space settings as well as in various frequency regimes. By using the asymptotics of the Green function, the analysis of the imaging functional of the time-reversal wave fields becomes possible, which clearly demonstrates the super-resolution property. The result provides the first mathematical theory of super-resolution in the context of a deterministic medium and sheds light to the mechanism of super-resolution and super-focusing for waves in deterministic complex media.

preprint2014arXiv

Admittivity imaging from multi-frequency micro-electrical impedance tomography

The aim of this paper is to propose an optimal control optimization algorithm for reconstructing admittivity distributions (i.e., both conductivity and permittivity) from multi-frequency micro-electrical impedance tomography. A convergent and stable optimization scheme is shown to be obtainable from multi-frequency data. The results of this paper have potential applicability in cancer imaging, cell culturing and differentiation, food sciences, and biotechnology.

preprint2014arXiv

Detection and classification from electromagnetic induction data

In this paper we introduce an efficient algorithm for identifying conductive objects using induction data derived from eddy currents. Our method consists of first extracting geometric features from the induction data and then matching them to precomputed data for known objects from a given dictionary. The matching step relies on fundamental properties of conductive polarization tensors and new invariants introduced in this paper. A new shape identification scheme is introduced and studied. We test it numerically in the presence of measurement noise. Stability and resolution capabilities of the proposed identification algorithm are quantified in numerical simulations.

preprint2014arXiv

Electrical impedance tomography-based pressure-sensing using conductive membrane

This paper presents a mathematical framework for a flexible pressure-sensor model using electrical impedance tomography (EIT). When pressure is applied to a conductive membrane patch with clamped boundary, the pressure-induced surface deformation results in a change in the conductivity distribution. This change can be detected in the current-voltage data ({\it i.e.,} EIT data) measured on the boundary of the membrane patch. Hence, the corresponding inverse problem is to reconstruct the pressure distribution from the data. The 2D apparent conductivity (in terms of EIT data) corresponding to the surface deformation is anisotropic. Thus, we consider a constrained inverse problem by restricting the coefficient tensor to the range of the map from pressure to 2D-apparent conductivity. This paper provides theoretical grounds for the mathematical model of the inverse problem. We develop a reconstruction algorithm based on a careful sensitivity analysis. We demonstrate the performance of the reconstruction algorithm through numerical simulations to validate its feasibility for future experimental studies.

preprint2014arXiv

Magnetic resonance-based reconstruction method of conductivity and permittivity distributions at the Larmor frequency

Magnetic resonance electrical property tomography is a recent medical imaging modality for visualizing the electrical tissue properties of the human body using radio-frequency magnetic fields. It uses the fact that in magnetic resonance imaging systems the eddy currents induced by the radio-frequency magnetic fields reflect the conductivity ($σ$) and permittivity ($ε$) distributions inside the tissues through Maxwell's equations. The corresponding inverse problem consists of reconstructing the admittivity distribution ($γ=σ+iωε$) at the Larmor frequency ($ω/2π=$128 MHz for a 3 tesla MRI machine) from the positive circularly polarized component of the magnetic field ${\bf H}=(H_x,H_y,H_z)$. Previous methods are usually based on an assumption of local homogeneity ($\nablaγ\approx 0$) which simplifies the governing equation. However, previous methods that include the assumption of homogeneity are prone to artifacts in the region where $γ$ varies. Hence, recent work has sought a reconstruction method that does not assume local-homogeneity. This paper presents a new magnetic resonance electrical property tomography reconstruction method which does not require any local homogeneity assumption on $γ$. We find that $γ$ is a solution of a semi-elliptic partial differential equation with its coefficients depending only on the measured data $H^+$, which enable us to compute a blurred version of $γ$. To improve the resolution of the reconstructed image, we developed a new optimization algorithm that minimizes the mismatch between the data and the model data as a highly nonlinear function of $γ$. Numerical simulations are presented to illustrate the potential of the proposed reconstruction method.

preprint2014arXiv

Mathematical modeling in full-field optical coherence elastography

We provide a mathematical analysis of and a numerical framework for full-field optical coherence elastography, which has unique features including micron-scale resolution, real-time processing, and non-invasive imaging. We develop a novel algorithm for transforming volumetric optical images before and after the mechanical solicitation of a sample with sub-cellular resolution into quantitative shear modulus distributions. This has the potential to improve sensitivities and specificities in the biological and clinical applications of optical coherence tomography.

preprint2014arXiv

Mathematical modeling of mechanical vibration assisted conductivity imaging

This paper aims at mathematically modeling a new multi-physics conductivity imaging system incorporating mechanical vibrations simultaneously applied to an imaging object together with current injections. We perturb the internal conductivity distribution by applying time-harmonic mechanical vibrations on the boundary. This enhances the effects of any conductivity discontinuity on the induced internal current density distribution. Unlike other conductivity contrast enhancing frameworks, it does not require a prior knowledge of a reference data. In this paper, we provide a mathematical framework for this novel imaging modality. As an application of the vibration-assisted impedance imaging framework, we propose a new breast image reconstruction method in electrical impedance tomography (EIT). As its another application, we investigate a conductivity anomaly detection problem and provide an efficient location search algorithm. We show both analytically and numerically that the applied mechanical vibration increases the data sensitivity to the conductivity contrast and enhances the quality of reconstructed images and anomaly detection results. For numerous applications in impedance imaging, the proposed multi-physics method opens a new difference imaging area called the vibration-difference imaging, which can augment the time-difference and also frequency-difference imaging methods for sensitivity improvements.

preprint2014arXiv

Reconstruction and stability in acousto-optic imaging for absorption maps with bounded variation

The aim of this paper is to propose for the first time a reconstruction scheme and a stability result for recovering from acoustic-optic data absorption distributions with bounded variation. The paper extends earlier results on smooth absorption distributions. It opens a door for a mathematical and numerical framework for imaging, from internal data, parameter distributions with high contrast in biological tissues.

preprint2014arXiv

Stability Analysis for Magnetic Resonance Elastography

We consider the inverse problem of finding unknown elastic parameters from internal measurements of displacement fields for tissues. The measurements are made on the entirety of a smooth domain. Since tissues can be modeled as quasi-incompressible fluids, we examine the Stokes system and consider only the recovery of shear modulus distributions. Our main result is to establish Lipschitz stable estimates on the shear modulus distributions from internal measurements of displacement fields. These estimates imply convergence of a numerical scheme known as the Landweber iteration scheme for reconstructing the shear modulus distributions.

preprint2014arXiv

Super-resolution in Imaging High Contrast Targets from the Perspective of Scattering Coefficients

In this paper we consider the inverse scattering problem for high-contrast targets. We mathematically analyze the experimentally-observed phenomenon of super-resolution in imaging the target shape. This is the first time that a mathematical theory of super-resolution has been established in the context of imaging high contrast inclusions. We illustrate our main findings with a variety of numerical examples. Our analysis is based on the novel concept of scattering coefficients. These findings may help in developing resonant structures for resolution enhancement.

preprint2014arXiv

Time-domain multiscale shape identification in electro-sensing

This paper presents premier and innovative time-domain multi-scale method for shape identification in electro-sensing using pulse-type signals. The method is based on transform-invariant shape descriptors computed from filtered polarization tensors at multi-scales. The proposed algorithm enjoys a remarkable noise robustness even with far-field measurements at very limited angle of view. It opens a door for pulsed imaging using echolocation and induction data.

preprint2013arXiv

Anomalous localized resonance using a folded geometry in three dimensions

If a body of dielectric material is coated by a plasmonic structure of negative dielectric material with nonzero loss parameter, then cloaking by anomalous localized resonance (CALR) may occur as the loss parameter tends to zero. It was proved in other papers by authors that if the coated structure is circular (2D) and dielectric constant of the shell is a negative constant (with loss parameter), then CALR occurs, and if the coated structure is spherical (3D), then CALR does not occur. The aim of this paper is to show that the CALR takes place if the spherical coated structure has a specially designed anisotropic dielectric tensor. The anisotropic dielectric tensor is designed by unfolding a folded geometry.

preprint2013arXiv

Backpropagation Imaging in Nonlinear Harmonic Holography in the Presence of Measurement and Medium Noises

In this paper, the detection of a small reflector in a randomly heterogenous medium using second-harmonic generation is investigated. The medium is illuminated by a time-harmonic plane wave at frequency omega. It is assumed that the reflector has a non-zero second-order nonlinear susceptibility, and thus emits a wave at frequency two omega in addition to the fundamental frequency linear scattering. It is shown how the fundamental frequency signal and the second-harmonic signal propagate in the medium. A statistical study of the images obtained by migrating the boundary data is performed. It is proved that the second-harmonic image is more stable with respect to medium noise than the one obtained with the fundamental signal. Moreover, the signal-to-noise ratio for the second-harmonic image does not depend neither on the second-order susceptibility tensor nor on the volume of the particle.

preprint2013arXiv

Mathematical modeling of fluorescence diffuse optical imaging of cell membrane potential changes

The aim of this paper is to provide a mathematical model for spatial distribution of membrane electrical potential changes by fluorescence diffuse optical tomography. We derive the resolving power of the imaging method in the presence of measurement noise. The proposed mathematical model can be used for cell membrane tracking with the resolution of the optical microscope.

preprint2013arXiv

Modeling active electrolocation in weakly electric fish

In this paper, we provide a mathematical model for the electrolocation in weakly electric fishes. We first investigate the forward complex conductivity problem and derive the approximate boundary conditions on the skin of the fish. Then we provide a dipole approximation for small targets away from the fish. Based on this approximation, we obtain a non-iterative location search algorithm using multi-frequency measurements. We present numerical experiments to illustrate the performance and the stability of the proposed multi-frequency location search algorithm. Finally, in the case of disk- and ellipse-shaped targets, we provide a method to reconstruct separately the conductivity, the permittivity, and the size of the targets from multi-frequency measurements.

preprint2013arXiv

Optimal Shape Design by Partial Spectral Data

In this paper, we are concerned with a shape design problem, in which our target is to design, up to rigid transformations and scaling, the shape of an object given either its polarization tensor at multiple contrasts or the partial eigenvalues of its Neumann-Poincaré operator, which are known as the Fredholm eigenvalues. We begin by proposing to recover the eigenvalues of the Neumann-Poincaré operator from the polarization tensor by means of the holomorphic functional calculus. Then we develop a regularized Gauss-Newton optimization method for the shape reconstruction process. We present numerical results to demonstrate the effectiveness of the proposed methods and to illustrate important properties of the Fredholm eigenvalues and their associated eigenfunctions. Our results are expected to have important applications in the design of plasmon resonances in nanoparticles as well as in the multifrequency or pulsed imaging of small anomalies.

preprint2013arXiv

Passive Array Correlation-Based Imaging in a Random Waveguide

We consider reflector imaging in a weakly random waveguide. We address the situation in which the source is farther from the reflector to be imaged than the energy equipartition distance, but the receiver array is closer to the reflector to be imaged than the energy equipartition distance. As a consequence, the reflector is illuminated by a partially coherent field and the signals recorded by the receiver array are noisy. This paper shows that migration of the recorded signals cannot give a good image, but an appropriate migration of the cross correlations of the recorded signals can give a very good image. The resolution and stability analysis of this original functional shows that the reflector can be localized with an accuracy of the order of the wavelength even when the receiver array has small aperture, and that broadband sources are necessary to ensure statistical stability, whatever the aperture of the array.

preprint2013arXiv

Radiative transfer and diffusion limits for wave field correlations in locally shifted random media

The aim of this paper is to develop a mathematical framework for opto-elastography. In opto-elastography, a mechanical perturbation of the medium produces a decorrelation of optical speckle patterns due to the displacements of optical scatterers. To model this, we consider two optically random media, with the second medium obtained by shifting the first medium in some local region. We derive the radiative transfer equation for the cross-correlation of the wave fields in the media. Then we derive its diffusion approximation. In both the radiative transfer and the diffusion regimes, we relate the correlation of speckle patterns to the solutions of the radiative transfer and the diffusion equations. We present numerical simulations based on our model which are in agreement with recent experimental measurements.

preprint2013arXiv

Reconstruction of Inhomogeneous Conductivities via the Concept of Generalized Polarization Tensors

This paper extends the concept of generalized polarization tensors (GPTs), which was previously defined for inclusions with homogeneous conductivities, to inhomogeneous conductivity inclusions. We begin by giving two slightly different but equivalent definitions of the GPTs for inhomogeneous inclusions. We then show that, as in the homogeneous case, the GPTs are the basic building blocks for the far-field expansion of the voltage in the presence of the conductivity inclusion. Relating the GPTs to the Neumann-to-Dirichlet (NtD) map, it follows that the full knowledge of the GPTs allows unique determination of the conductivity distribution. Furthermore, we show important properties of the the GPTs, such as symmetry and positivity, and derive bounds satisfied by their harmonic sums. We also compute the sensitivity of the GPTs with respect to changes in the conductivity distribution and propose an algorithm for reconstructing conductivity distributions from their GPTs. This provides a new strategy for solving the highly nonlinear and ill-posed inverse conductivity problem. We demonstrate the viability of the proposed algorithm by preforming a sensitivity analysis and giving some numerical examples.

preprint2013arXiv

Shape identification and classification in echolocation

The paper aims at proposing the first shape identification and classification algorithm in echolocation. The approach is based on first extracting geometric features from the reflected waves and then matching them with precomputed ones associated with a dictionary of targets. The construction of such frequency-dependent shape descriptors is based on some important properties of the scattering coefficients and new invariants. The stability and resolution of the proposed identification algorithm with respect to measurement noise and the limited-view aspect are analytically and numerically quantified.

preprint2013arXiv

Shape recognition and classification in electro-sensing

This paper aims at advancing the field of electro-sensing. It exhibits the physical mechanism underlying shape perception for weakly electric fish. These fish orient themselves at night in complete darkness by employing their active electrolocation system. They generate a stable, high-frequency, weak electric field and perceive the transdermal potential modulations caused by a nearby target with different admittivity than the surrounding water. In this paper, we explain how weakly electric fish might identify and classify a target, knowing by advance that the latter belongs to a certain collection of shapes. Our model of the weakly electric fish relies on differential imaging, i.e., by forming an image from the perturbations of the field due to targets, and physics-based classification. The electric fish would first locate the target using a specific location search algorithm. Then it could extract, from the perturbations of the electric field, generalized (or high-order) polarization tensors of the target. Computing, from the extracted features, invariants under rigid motions and scaling yields shape descriptors. The weakly electric fish might classify a target by comparing its invariants with those of a set of learned shapes. On the other hand, when measurements are taken at multiple frequencies, the fish might exploit the shifts and use the spectral content of the generalized polarization tensors to dramatically improve the stability with respect to measurement noise of the classification procedure in electro-sensing. Surprisingly, it turns out that the first-order polarization tensor at multiple frequencies could be enough for the purpose of classification. A procedure to eliminate the background field in the case where the permittivity of the surrounding medium can be neglected, and hence improve further the stability of the classification process, is also discussed.

preprint2013arXiv

Spectroscopic imaging of a dilute cell suspension

A rigorous homogenization theory is derived to describe the effective admittivity of cell suspensions. A new formula is reported for dilute cases that gives the frequency-dependent effective admittivity with respect to the membrane polarization. Different microstructures are shown to be distinguishable via spectroscopic measurements of the overall admittivity using the spectral properties of the membrane polarization. The Debye relaxation times associated with the membrane polarization tensor are shown to be able to give the microscopic structure of the medium. A natural measure of the admittivity anisotropy is introduced and its dependence on the frequency of applied current is derived. A Maxwell-Wagner-Fricke formula is given for concentric circular cells, and the results can be extended to the random cases. A randomly deformed periodic medium is also considered and a new formula is derived for the overall admittivity of a dilute suspension of randomly deformed cells.

preprint2013arXiv

Target Detection and Characterization from Electromagnetic Induction Data

The goal of this paper is to contribute to the field of nondestructive testing by eddy currents. We provide a mathematical analysis and a numerical framework for simulating the imaging of arbitrarily shaped small volume conductive inclusions from electromagnetic induction data. We derive, with proof, a small-volume expansion of the eddy current data measured away from the conductive inclusion. The formula involves two polarization tensors: one associated with the magnetic contrast and the second with the conductivity of the inclusion. Based on this new formula, we design a location search algorithm. We include in this paper a discussion on data sampling, noise reduction, and on probability of detection. We provide numerical examples that support our findings.

preprint2013arXiv

Target Identification Using Dictionary Matching of Generalized Polarization Tensors

The aim of this paper is to provide a fast and efficient procedure for (real-time) target identification in imaging based on matching on a dictionary of precomputed generalized polarization tensors (GPTs). The approach is based on some important properties of the GPTs and new invariants. A new shape representation is given and numerically tested in the presence of measurement noise. The stability and resolution of the proposed identification algorithm is numerically quantified.

preprint2013arXiv

The Concept of Heterogeneous Scattering Coefficients and Its Application in Inverse Medium Scattering

This work investigates the scattering coefficients for inverse medium scattering problems. It shows some fundamental properties of the coefficients such as symmetry and tensorial properties. The relationship between the scattering coefficients and the far-field pattern is also derived. Furthermore, the sensitivity of the scattering coefficients with respect to changes in the permittivity and permeability distributions is investigated. In the linearized case, explicit formulas for reconstructing permittivity and permeability distributions from the scattering coefficients is proposed. They relate the exponentially ill-posed character of the inverse medium scattering problem at a fixed frequency to the exponential decay of the scattering coefficients. Moreover, they show the stability of the reconstruction from multifrequency measurements. This provides a new direction for solving inverse medium scattering problems.

preprint2012arXiv

Boundary perturbations due to the presence of small linear cracks in an elastic body

In this paper, Neumann cracks in elastic bodies are considered. We establish a rigorous asymptotic expansion for the boundary perturbations of the displacement (and traction) vectors that are due to the presence of a small elastic linear crack. The formula reveals that the leading order term is ε^2 where εis the length of the crack, and the ε^3-term vanishes. We obtain an asymptotic expansion of the elastic potential energy as an immediate consequence of the boundary perturbation formula. The derivation is based on layer potential techniques. It is expected that the formula would lead to very effective direct approaches for locating a collection of small elastic cracks and estimating their sizes and orientations.

preprint2012arXiv

Enhancement of near-cloaking for the full Maxwell equations

In this paper, we consider near cloaking for the full Maxwell equations. We extend the recent results, where the quasi-static limit case and the Helmholtz equation are considered, to electromagnetic scattering problems. We construct very effective near cloaking structures for the electromagnetic scattering problem at a fixed frequency. These new structures are, before using the transformation optics, layered structures and are designed so that their first scattering coefficients vanish. Inside the cloaking region, any target has near-zero scattering cross section for a band of frequencies. We analytically show that our new construction significantly enhances the cloaking effect for the full Maxwell equations.

preprint2012arXiv

Localization, Stability, and Resolution of Topological Derivative Based Imaging Functionals in Elasticity

The focus of this work is on rigorous mathematical analysis of the topological derivative based detection algorithms for the localization of an elastic inclusion of vanishing characteristic size. A filtered quadratic misfit is considered and the performance of the topological derivative imaging functional resulting therefrom is analyzed. Our analysis reveals that the imaging functional may not attain its maximum at the location of the inclusion. Moreover, the resolution of the image is below the diffraction limit. Both phenomena are due to the coupling of pressure and shear waves propagating with different wave speeds and polarization directions. A novel imaging functional based on the weighted Helmholtz decomposition of the topological derivative is, therefore, introduced. It is thereby substantiated that the maximum of the imaging functional is attained at the location of the inclusion and the resolution is enhanced and it proves to be the diffraction limit. Finally, we investigate the stability of the proposed imaging functionals with respect to measurement and medium noises.

preprint2012arXiv

Quantitative thermo-acoustic imaging: An exact reconstruction formula

This paper aims to mathematically advance the field of quantitative thermo-acoustic imaging. Given several electromagnetic data sets, we establish for the first time an analytical formula for reconstructing the absorption coefficient from thermal energy measurements. Since the formula involves derivatives of the given data up to the third order, it is unstable in the sense that small measurement noises may cause large errors. However, in the presence of measurement noise, the obtained formula, together with a noise regularization technique, provides a good initial guess for the true absorption coefficient. We finally correct the errors by deriving a reconstruction formula based on the least square solution of an optimal control problem and prove that this optimization step reduces the errors occurring and enhances the resolution.

preprint2012arXiv

Reconstruction of a piecewise smooth absorption coefficient by an acousto-optic process

The aim of this paper is to tackle the nonlinear optical reconstruction problem. Given a set of acousto-optic measurements, we develop a mathematical framework for the reconstruction problem in the case where the optical absorption distribution is supposed to be a perturbation of a piecewise constant function. Analyzing the acousto-optic measurements, we establish a new equation in the sense of distributions for the optical absorption coefficient. For doing so, we introduce a weak Helmholtz decomposition and interpret in a weak sense the cross-correlation measurements using the spherical Radon transform. We next show how to find an initial guess for the unknown coefficient and finally construct the true coefficient by providing a Landweber type iteration and proving that the resulting sequence converges to the solution of the system constituted by the optical diffusion equation and the new equation mentioned above. Our results in this paper generalize the acousto-optic process proposed in [Ammari, Garnier, Nguyen and Seppecher, 2012] for piecewise smooth optical absorption distributions.

preprint2012arXiv

Spectral theory of a Neumann-Poincaré-type operator and analysis of anomalous localized resonance II

If a body of dielectric material is coated by a plasmonic structure of negative dielectric constant with nonzero loss parameter, then cloaking by anomalous localized resonance (CALR) may occur as the loss parameter tends to zero. The aim of this paper is to investigate this phenomenon in two and three dimensions when the coated structure is radial, and the core, shell and matrix are isotropic materials. In two dimensions, we show that if the real part of the permittivity of the shell is -1 (under the assumption that the permittivity of the background is 1), then CALR takes place. If it is different from -1, then CALR does not occur. In three dimensions, we show that CALR does not occur. The analysis of this paper reveals that occurrence of CALR is determined by the eigenvalue distribution of the Neumann-Poincaré-type operator associated with the structure.

preprint2012arXiv

Spectral theory of a Neumann-Poincare-type operator and analysis of cloaking due to anomalous localized resonance

The aim of this paper is to give a mathematical justification of cloaking due to anomalous localized resonance (CALR). We consider the dielectric problem with a source term in a structure with a layer of plasmonic material. Using layer potentials and symmetrization techniques, we give a necessary and sufficient condition on the fixed source term for electromagnetic power dissipation to blow up as the loss parameter of the plasmonic material goes to zero. This condition is written in terms of the Newtonian potential of the source term. In the case of concentric disks, we make the condition even more explicit. Using the condition, we are able to show that for any source supported outside a critical radius CALR does not take place, and for sources located inside the critical radius satisfying certain conditions CALR does take place as the loss parameter goes to zero.

preprint2012arXiv

Stability estimates for an inverse scattering problem at high frequencies

We consider an inverse scattering problem and its near-field approximation at high frequencies. We first prove, for both problems, Lipschitz stability results for determining the low-frequency component of the potential. Then we show that, in the case of a radial potential supported sufficiently near the boundary, infinite resolution can be achieved from measurements of the near-field operator in the monotone case.

preprint2012arXiv

Strong convergence of the solutions of the linear elasticity and uniformity of asymptotic expansions in the presence of small inclusions

We consider the Lamé system of linear elasticity when the inclusion has the extreme elastic constants. We show that the solutions to the Lamé system converge in appropriate $H^1$-norms when the shear modulus tends to infinity (the other modulus, the compressional modulus is fixed), and when the bulk modulus and the shear modulus tend to zero. Using this result, we show that the asymptotic expansion of the displacement vector in the presence of small inclusion is uniform with respect to Lamé parameters.

preprint2012arXiv

Tracking of a Mobile Target Using Generalized Polarization Tensors

In this paper we apply an extended Kalman filter to track both the location and the orientation of a mobile target from multistatic response measurements. We also analyze the effect of the limited-view aspect on the stability and the efficiency of our tracking approach. Our algorithm is based on the use of the generalized polarization tensors, which can be reconstructed from the multistatic response measurements by solving a linear system. The system has the remarkable property that low order generalized polarization tensors are not affected by the error caused by the instability of higher orders in the presence of measurement noise.

preprint2011arXiv

Enhancement of Near Cloaking Using Generalized Polarization Tensors Vanishing Structures. Part I: The Conductivity Problem

The aim of this paper is to provide an original method of constructing very effective near-cloaking structures for the conductivity problem. These new structures are such that their first Generalized Polarization Tensors vanish. We show that this in particular significantly enhances the cloaking effect. We then present some numerical examples of Generalized Polarization Tensors vanishing structures.

preprint2011arXiv

Enhancement of near-cloaking. Part II: the Helmholtz equation

The aim of this paper is to extend the method of improving cloaking structures in the conductivity to scattering problems. We construct very effective near-cloaking structures for the scattering problem at a fixed frequency. These new structures are, before using the transformation optics, layered structures and are designed so that their first scattering coefficients vanish. Inside the cloaking region, any target has near-zero scattering cross section for a band of frequencies. We analytically show that our new construction significantly enhances the cloaking effect for the Helmholtz equation.

preprint2011arXiv

Sharp estimates for the Neumann functions and applications to quantitative photo-acoustic imaging in inhomogeneous media

We obtain sharp $L^p$ and Hölder estimates for the Neumann function of the operator $\nabla \cdot γ\nabla - ik$ on a bounded domain. We also obtain quantitative description of its singularity. We then apply these estimates to quantitative photo-acoustic imaging in inhomogeneous media. The problem is to reconstruct the optical absorption coefficient of a diametrically small anomaly from the absorbed energy density.

preprint2010arXiv

Layer Potential Techniques for the Narrow Escape Problem

The narrow escape problem consists of deriving the asymptotic expansion of the solution of a drift-diffusion equation with the Dirichlet boundary condition on a small absorbing part of the boundary and the Neumann boundary condition on the remaining reflecting boundaries. Using layer potential techniques, we rigorously find high-order asymptotic expansions of such solutions. We explicitly show the nonlinear interaction of many small absorbing targets. Based on the asymptotic theory for eigenvalue problems developed in \cite{book}, we also construct high-order asymptotic formulas for eigenvalues of the Laplace and the drifted Laplace operators for mixed boundary conditions on large and small pieces of the boundary.