• ### Shape and size dependence of dipolar plasmonic resonance of nanoparticles(1804.11092)

April 30, 2018 math.AP, math-ph, math.MP, math.SP
The aim of this paper is to present a new approach, based on singular volume integral equations, in order to compute the size dependency of plasmonic resonances. The paper also provides rigorous derivations of the extinction and absorption cross sections for elliptical particles.
• ### Wave-Field Shaping in Cavities by Tunable Metasurfaces(1804.10912)

April 29, 2018 math.AP
Cavities, because they trap waves for long times due to their reflecting walls, are used in a vast number of scientific domains. Indeed, in these closed media and due to interferences, the free space continuum of solutions becomes a discrete set of stationary eigenmodes. These enhanced stationary fields are commonly used in fundamental physics to increase wave-matter interactions. The eigenmodes and associated eigenfrequencies of a cavity are imposed by its geometrical properties through the boundary conditions. In this paper, we show that one can control the wave fields created by point sources inside cavities by tailoring only the boundaries of the cavities. This is achieved through the use of a tunable reflecting metasurface, which is part of the frontiers of the cavity, and can switch its boundary conditions from Dirichlet to Neumann. Based on the use of arrays of subwavelength resonators, a mathematical modeling of the physical mechanism underlying the concept of tunable metasurfaces is provided.
• ### Reconstructing fine details of small objects by using plasmonic spectroscopic data. Part II: The strong interaction regime(1712.02639)

April 16, 2018 math.AP, math-ph, math.MP, math.SP
This paper is concerned with the inverse problem of reconstructing a small object from far field measurements by using the field interaction with a plasmonic particle which can be viewed as a passive sensor. It is a follow-up of the work [H. Ammari et al., Reconstructing fine details of small objects by using plasmonic spectroscopic data, SIAM J. Imag. Sci., to appear], where the intermediate interaction regime was considered. In that regime, it was shown that the presence of the target object induces small shifts to the resonant frequencies of the plasmonic particle. These shifts, which can be determined from the far field data, encodes the contracted generalized polarization tensors of the target object, from which one can perform reconstruction beyond the usual resolution limit. The main argument is based on perturbation theory. However, the same argument is no longer applicable in the strong interaction regime as considered in this paper due to the large shift induced by strong field interaction between the particles. We develop a novel technique based on conformal mapping theory to overcome this difficulty. The key is to design a conformal mapping which transforms the two particle system into a shell-core structure, in which the inner dielectric core corresponds to the target object. We show that a perturbation argument can be used to analyze the shift in the resonant frequencies due to the presence of the inner dielectric core. This shift also encodes information of the contracted polarization tensors of the core, from which one can reconstruct its shape, and hence the target object. Our theoretical findings are supplemented by a variety of numerical results based on an efficient optimal control algorithm. The results of this paper make the mathematical foundation for plasmonic sensing complete.
• ### Subwavelength localized modes for acoustic waves in bubbly crystals with a defect(1804.00756)

April 2, 2018 math.AP
The ability to control wave propagation is of fundamental interest in many areas of physics. Photonic and phononic crystals proved very useful for this purpose but, because they are based on Bragg interferences, 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 bubble 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 increase 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.
• ### Dynamic Spike Super-resolution and Applications to Ultrafast Ultrasound Imaging(1803.03251)

March 8, 2018 math.NA, cs.IT, math.IT
We consider the dynamical super-resolution problem consisting in the recovery of positions and velocities of moving particles from low-frequency static measurements taken over multiple time steps. The standard approach to this issue is a two-step process: first, at each time step some static reconstruction method is applied to locate the positions of the particles with super-resolution and, second, some tracking technique is applied to obtain the velocities. In this paper we propose a fully dynamical method based on a phase-space lifting of the positions and the velocities of the particles, which are simultaneously reconstructed with super-resolution. We provide a rigorous mathematical analysis of the recovery problem, both for the noiseless case and in presence of noise. Several numerical simulations illustrate and validate our method, which shows some advantage over existing techniques. We then discuss the application of this approach to the dynamical super-resolution problem in ultrafast ultrasound imaging: blood vessels' locations and blood flow velocities are recovered with super-resolution.
• ### Double-negative acoustic metamaterials(1709.08177)

March 8, 2018 math.AP, math-ph, math.MP
The aim of this paper is to provide a mathematical theory for understanding the mechanism behind the double-negative refractive index phenomenon in bubbly fluids. The design of double-negative metamaterials generally requires the use of two different kinds of subwavelength resonators, which may limits the applicability of double-negative metamaterials. Herein we rely on media that consists of only a single type of resonant element, and show how to turn the acoustic metamaterial with a single negative effective property obtained in [H. Ammari and H. Zhang, Effective medium theory for acoustic waves in bubbly fluids near Minnaert resonant frequency. SIAM J. Math. Anal., 49 (2017), 3252--3276.] into a negative refractive index metamaterial, which refracts waves negatively, hence acting as a superlens. Using bubble dimers made of two identical bubbles, it is proved that both the effective mass density and the bulk modulus of the bubbly fluid can be negative near the anti-resonance of the two hybridized Minnaert resonances for a single constituent bubble dimer. A rigorous justification of the Minnaert resonance hybridization, in the case of a bubble dimer in a homogeneous medium, is established. The acoustic properties of a single bubble dimer are analyzed. Asymptotic formulas for the two hybridized Minnaert resonances are derived. Moreover, it is proved that the bubble dimer can be approximated by a point scatterer with monopole and dipole modes. For an appropriate volume fraction of bubble dimers with certain conditions on their configuration, a double-negative effective medium when the frequency is near the anti-resonance of the hybridized Minnaert resonances can be obtained.
• ### A mathematical and numerical framework for near-field optics(1803.01155)

March 3, 2018 math.AP
This paper is concerned with the inverse problem of reconstructing small and local perturbations of a planar surface using the field interaction between a known plasmonic particle and the planar surface. The aim is to perform a super-resolved reconstruction of these perturbations from shifts in the plasmonic frequencies of the particle-surface system. In order to analyze the interaction between the plasmonic particle and the planar surface, a well chosen conformal mapping, which transforms the particle-surface system into a coated structure, is used. Then the even Fourier coefficients of the transformed domain are related to the shifts in the plasmonic resonances of the particle-surface system. A direct reconstruction of the perturbations of the planar surface is proposed. Its viability and limitations are documented by numerical examples.
• ### Bloch waves in bubbly crystal near the first band gap: a high-frequency homogenization approach(1708.07955)

Jan. 25, 2018 math.AP, math-ph, math.MP
This paper is concerned with the high-frequency homogenization of bubbly phononic crystals. It is a follow-up of the works [H. Ammari et al., Sub-wavelength phononic bandgap opening in bubbly media, J. Diff. Eq., 263 (2017), 5610--5629] which shows the existence of a sub-wavelength band gap. This phenomena can be explained by the periodic inference of cell resonance which is due to the high contrast in both the density and bulk modulus between the bubbles and the surrounding medium. In this paper, we prove that the first Bloch eigenvalue achieves its maximum at the corner of the Brillouin zone. Moreover, by computing the asymptotic of the Bloch eigenfunctions in the periodic structure near that critical frequency, we demonstrate that these eigenfunctions can be decomposed into two parts: one part is slowly varying and satisfies a homogenized equation, while the other is periodic across each elementary crystal cell and is varying. They rigorously justify, in the nondilute case, the observed super-focusing of acoustic waves in bubbly crystals near and below the maximum of the first Bloch eigenvalue and confirm the band gap opening near and above this critical frequency.
• ### The point-interaction approximation for the fields generated by contrasted bubbles at arbitrary fixed frequencies(1801.07921)

Jan. 24, 2018 math.AP
We deal with the linearized model of the acoustic wave propagation generated by small bubbles in the harmonic regime. We estimate the waves generated by a cluster of $M$ small bubbles, distributed in a bounded domain $\Omega$, with relative densities having contrasts of the order $a^{\beta}, \beta>0,$ where $a$ models their relative maximum radius, $a\ll 1$. We provide useful and natural conditions on the number $M$, the minimum distance and the contrasts parameter $\beta$ of the small bubbles under which the point interaction approximation (called also the Foldy-Lax approximation) is valid. With the regimes allowed by our conditions, we can deal with a general class of such materials. Applications of these expansions in material sciences and imaging are immediate. For instance, they are enough to derive and justify the effective media of the cluster of the bubbles for a class of gases with densities having contrasts of the order $a^{\beta}$, $\beta \in (\frac{3}{2}, 2)$ and in this case we can handle any fixed frequency. In the particular and important case $\beta=2$, we can handle any fixed frequency far or close (but distinct) from the corresponding Minnaert resonance.
• ### Field expansions for systems of strongly coupled plasmonic nanoparticles(1709.08469)

Sept. 25, 2017 math.AP, math.NA
This paper is concerned with efficient representations and approximations of the solution to the scattering problem by a system of strongly coupled plasmonic particles. Three schemes are developed: the first is the resonant expansion which uses the resonant modes of the system of particles computed by a conformal transformation, the second is the hybridized resonant expansion which uses linear combinations of the resonant modes for each of the particles in the system as a basis to represent the solution, and the last one is the multipole expansion with respect to the origin. By considering a system formed by two plasmonic particles of circular shape, we demonstrate the relations between these expansion schemes and their advantages and disadvantages both analytically and numerically. In particular, we emphasize the efficiency of the resonant expansion scheme in approximating the near field of the system of particles. The difference between these plasmonic particle systems and the nonresonant dielectric particle system is also highlighted. The paper provides a guidance on the challenges for numerical simulations of strongly coupled plasmonic systems.
• ### Sub-wavelength focusing of acoustic waves in bubbly media(1706.08113)

June 25, 2017 math.AP
The purpose of this paper is to investigate acoustic wave scattering by a large number of bubbles in a liquid at frequencies near the Minnaert resonance frequency. This bubbly media has been exploited in practice to obtain super-focusing of acoustic waves. Using layer potential techniques we derive the scattering function for a single spherical bubble excited by an incident wave in the low frequency regime. We then derive the point scatter approximation for multiple scattering by N bubbles. We describe several numerical experiments based on the point scatterer approximation that demonstrate the possibility of achieving super-focusing using bubbly media.
• ### Plasmonic interaction between nanospheres(1604.02624)

When metallic (or plasmonic) nanospheres are nearly touching, strong concentration of light can occur in the narrow gap regions. This phenomenon has a potential application in nanophotonics, biosensing and spectroscopy. The understanding of the strong interaction between the plasmonic spheres turns out to be quite challenging. Indeed, an extremely high computational cost is required to compute the electromagnetic field. Also, the classical method of image charges, which is effective for dielectric spheres system, is not valid for plasmonic spheres because of their negative permittivities. Here we develop new analytical and numerical methods for the plasmonic spheres system by clarifying the connection between transformation optics and the method of image charges. We derive fully analytic solutions valid for two plasmonic spheres. We then develop a hybrid numerical scheme for computing the field distribution produced by an arbitrary number of spheres. Our method is highly efficient and accurate even in the nearly touching case and is valid for plasmonic spheres.
• ### Robust Cardiac Motion Estimation using Ultrafast Ultrasound Data: A Low-Rank-Topology-Preserving Approach(1611.02730)

April 25, 2017 cs.CV
Cardiac motion estimation is an important diagnostic tool to detect heart diseases and it has been explored with modalities such as MRI and conventional ultrasound (US) sequences. US cardiac motion estimation still presents challenges because of the complex motion patterns and the presence of noise. In this work, we propose a novel approach to estimate the cardiac motion using ultrafast ultrasound data. -- Our solution is based on a variational formulation characterized by the L2-regularized class. The displacement is represented by a lattice of b-splines and we ensure robustness by applying a maximum likelihood type estimator. While this is an important part of our solution, the main highlight of this paper is to combine a low-rank data representation with topology preservation. Low-rank data representation (achieved by finding the k-dominant singular values of a Casorati Matrix arranged from the data sequence) speeds up the global solution and achieves noise reduction. On the other hand, topology preservation (achieved by monitoring the Jacobian determinant) allows to radically rule out distortions while carefully controlling the size of allowed expansions and contractions. Our variational approach is carried out on a realistic dataset as well as on a simulated one. We demonstrate how our proposed variational solution deals with complex deformations through careful numerical experiments. While maintaining the accuracy of the solution, the low-rank preprocessing is shown to speed up the convergence of the variational problem. Beyond cardiac motion estimation, our approach is promising for the analysis of other organs that experience motion.
• ### Reconstructing fine details of small objects by using plasmonic spectroscopic data(1704.04870)

April 17, 2017 math.AP, math-ph, math.MP, math.SP
This paper is concerned with the inverse problem of reconstructing a small object from far field measurements. The inverse problem is severally ill-posed because of the diffraction limit and low signal to noise ratio. We propose a novel methodology to solve this type of inverse problems based on an idea from plasmonic sensing. By using the field interaction with a known plasmonic particle, the fine detail information of the small object can be encoded into the shift of the resonant frequencies of the two particle system in the far field. In the intermediate interaction regime, we show that this information is exactly the generalized polarization tensors associated with the small object, from which one can perform the reconstruction. Our theoretical findings are supplemented by a variety of numerical results. The results in the paper also provide a mathematical foundation for plasmonic sensing.
• ### Numerical determination of anomalies in multifrequency electrical impedance tomography(1704.04878)

April 17, 2017 math.NA
The multifrequency electrical impedance tomography consists in retrieving the conductivity distribution of a sample by injecting a finite number of currents with multiple frequencies. In this paper we consider the case where the conductivity distribution is piecewise constant, takes a constant value outside a single smooth anomaly, and a frequency dependent function inside the anomaly itself. Using an original spectral decomposition of the solution of the forward conductivity problem in terms of Poincar\'e variational eigenelements, we retrieve the Cauchy data corresponding to the extreme case of a perfect conductor, and the conductivity profile. We then reconstruct the anomaly from the Cauchy data. The numerical experiments are conducted using gradient descent optimization algorithms.
• ### Heat generation with plasmonic nanoparticles(1703.00422)

March 1, 2017 math.AP
In this paper we use layer potentials and asymptotic analysis techniques to analyze the heat generation due to nanoparticles when illuminated at their plasmonic resonance. We consider arbitrary-shaped particles and both single and multiple particles. For close-to-touching nanoparticles, we show that the temperature field deviates significantly from the one generated by a single nanoparticle. The results of this paper open a door for solving the challenging problems of detecting plasmonic nanoparticles in biological media and monitoring temperature elevation in tissue generated by nanoparticle heating.
• ### Subwavelength phononic bandgap opening in bubbly media(1702.05317)

Feb. 17, 2017 math.AP
The aim of this paper is to show both analytically and numerically the existence of a subwavelength phononic bandgap in bubble phononic crystals. The key is an original formula for the quasi-periodic Minnaert resonance frequencies of an arbitrarily shaped bubble. The main findings in this paper are illustrated with a variety of numerical experiments.
• ### Determining anisotropic conductivity using diffusion tensor imaging data in magneto-acoustic tomography with magnetic induction(1702.05187)

Feb. 16, 2017 math.NA
In this paper we present a mathematical and numerical framework for a procedure of imaging anisotropic electrical conductivity tensor by integrating magneto-acoutic tomography with data acquired from diffusion tensor imaging. Magneto-acoustic Tomography with Magnetic Induction (MAT-MI) is a hybrid, non-invasive medical imaging technique to produce conductivity images with improved spatial resolution and accuracy. Diffusion Tensor Imaging (DTI) is also a non- invasive technique for characterizing the diffusion properties of water molecules in tissues. We propose a model for anisotropic conductivity in which the conductivity is proportional to the diffusion tensor. Under this assumption, we propose an optimal control approach for reconstructing the anisotropic electrical conductivity tensor. We prove convergence and Lipschitz type stability of the algorithm and present numerical examples to illustrate its accuracy and feasibility.
• ### Two-Dimensional Elastic Scattering Coefficients and Enhancement of Nearly Elastic Cloaking(1606.08070)

Jan. 15, 2017 math-ph, math.MP
The concept of scattering coefficients has played a pivotal role in a broad range of inverse scattering and imaging problems in acoustic, and electromagnetic media. In view of their promising applications in inverse problems related to mathematical imaging and elastic cloaking, the notion of elastic scattering coefficients of an inclusion is presented from the perspective of boundary layer potentials and a few properties are discussed. A reconstruction algorithm is developed and analyzed for extracting the elastic scattering coefficients from multi-static response measurements of the scattered field in order to cater to inverse scattering problems. The decay rate, stability and error analyses, and the estimate of maximal resolving order in terms of the signal-to-noise ratio are discussed. Moreover, scattering-coefficients-vanishing structures are designed and their utility for enhancement of nearly elastic cloaking is elucidated.
• ### Mathematical modelling of the electric sense of fish: the role of multi-frequency measurements and movement(1612.03168)

Dec. 9, 2016 q-bio.TO
Understanding active electrolocation in weakly electric fish remains a challenging issue. In this article we propose a mathematical formulation of this problem, in terms of partial differential equations. This allows us to detail two algorithms: one for localizing a target using the multi-frequency aspect of the signal, and antoher one for identifying the shape of this target. Shape recognition is designed in a machine learning point of view, and take advantage of both the multi-frequency setup and the movement of the fish around its prey. Numerical simulations are shown for the computation of the electric field emitted and sensed by the fish; they are then used as an input for the two algorithms.
• ### Mathematical Analysis of Ultrafast Ultrasound Imaging(1604.04604)

Oct. 1, 2016 math.NA, math.PR
This paper provides a mathematical analysis of ultrafast ultrasound imaging. This newly emerging modality for biomedical imaging uses plane waves instead of focused waves in order to achieve very high frame rates. We derive the point spread function of the system in the Born approximation for wave propagation and study its properties. We consider dynamic data for blood flow imaging, and introduce a suitable random model for blood cells. We show that a singular value decomposition method can successfully remove the clutter signal by using the different spatial coherence of tissue and blood signals, thereby providing high-resolution images of blood vessels, even in cases when the clutter and blood speeds are comparable in magnitude. Several numerical simulations are presented to illustrate and validate the approach.
• ### Identification of an inclusion in Multifrequency Electric Impedance Tomography(1609.07091)

Sept. 21, 2016 math.AP
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.
• ### A signal separation technique for sub-cellular imaging using dynamic optical coherence tomography(1608.04382)

Aug. 15, 2016 math.NA, physics.med-ph
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.
• ### A mathematical and numerical framework for bubble meta-screens(1608.02733)

Aug. 9, 2016 math.AP
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].
• ### Mathematical framework for abdominal electrical impedance tomography to assess fatness(1607.05893)

July 20, 2016 math.NA
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.