
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.

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 wavematter 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.

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 followup
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 shellcore
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.

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), 56105629], 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.

We consider the dynamical superresolution problem consisting in the recovery
of positions and velocities of moving particles from lowfrequency static
measurements taken over multiple time steps. The standard approach to this
issue is a twostep process: first, at each time step some static
reconstruction method is applied to locate the positions of the particles with
superresolution and, second, some tracking technique is applied to obtain the
velocities. In this paper we propose a fully dynamical method based on a
phasespace lifting of the positions and the velocities of the particles, which
are simultaneously reconstructed with superresolution. 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
superresolution problem in ultrafast ultrasound imaging: blood vessels'
locations and blood flow velocities are recovered with superresolution.

The aim of this paper is to provide a mathematical theory for understanding
the mechanism behind the doublenegative refractive index phenomenon in bubbly
fluids. The design of doublenegative metamaterials generally requires the use
of two different kinds of subwavelength resonators, which may limits the
applicability of doublenegative 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),
32523276.] 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 antiresonance 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 doublenegative effective medium when the
frequency is near the antiresonance of the hybridized Minnaert resonances can
be obtained.

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
superresolved reconstruction of these perturbations from shifts in the
plasmonic frequencies of the particlesurface system. In order to analyze the
interaction between the plasmonic particle and the planar surface, a well
chosen conformal mapping, which transforms the particlesurface 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
particlesurface system. A direct reconstruction of the perturbations of the
planar surface is proposed. Its viability and limitations are documented by
numerical examples.

This paper is concerned with the highfrequency homogenization of bubbly
phononic crystals. It is a followup of the works [H. Ammari et al.,
Subwavelength phononic bandgap opening in bubbly media, J. Diff. Eq., 263
(2017), 56105629] which shows the existence of a subwavelength 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 superfocusing 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.

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 FoldyLax 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.

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.

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
superfocusing 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 superfocusing using bubbly media.

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.

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 L2regularized class. The displacement is
represented by a lattice of bsplines 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 lowrank data
representation with topology preservation. Lowrank data representation
(achieved by finding the kdominant 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 lowrank 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.

This paper is concerned with the inverse problem of reconstructing a small
object from far field measurements. The inverse problem is severally illposed
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.

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.

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 arbitraryshaped particles and both single and
multiple particles. For closetotouching 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.

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 quasiperiodic Minnaert resonance
frequencies of an arbitrarily shaped bubble. The main findings in this paper
are illustrated with a variety of numerical experiments.

In this paper we present a mathematical and numerical framework for a
procedure of imaging anisotropic electrical conductivity tensor by integrating
magnetoacoutic tomography with data acquired from diffusion tensor imaging.
Magnetoacoustic Tomography with Magnetic Induction (MATMI) is a hybrid,
noninvasive 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.

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 multistatic 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 signaltonoise ratio are discussed. Moreover,
scatteringcoefficientsvanishing structures are designed and their utility for
enhancement of nearly elastic cloaking is elucidated.

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 multifrequency
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 multifrequency 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.

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
highresolution 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.

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.

This paper aims at imaging the dynamics of metabolic activity of cells. Using
dynamic optical coherence tomography, we introduce a new multiparticle
dynamical model to simulate the movements of the collagen and the cell
metabolic activity and develop an efficient signal separation technique for
subcellular imaging. We perform a singularvalue 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.

The aim of this paper is to provide a mathematical and numerical framework
for the analysis and design of bubble metascreens. An acoustic metascreen 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
superabsorption behavior observed in [V. Leroy et al., Phys. Rev. B, 2015].

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
depthbased reconstruction method that uses a specially chosen current pattern
to construct referencelike 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 32channel EIT system and human
like domain.