
The NeumannPoincar\'e operator on the sphere has $\frac{1}{2(2k+1)}$,
$k=0,1,2,\ldots$, as its eigenvalues and the corresponding multiplicity is
$2k+1$. We consider the bifurcation of eigenvalues under deformation of
domains, and show that Frech\'et derivative of the sum of the bifurcations is
zero. We then discuss the connection of this result with some conjectures
regarding the NeumannPoincar\'e operator.

Upon insertion of an inclusion into a medium with the uniform field, if the
field is not perturbed at all outside the inclusion, then it is called a
neutral inclusion. It is called a weakly neutral inclusion if the field is
perturbed weakly. The inclusions neutral to multiple uniform fields are of
circular shape if the medium is isotropic, and any other shape cannot be
neutral. We consider in this paper the problem of constructing inclusions of
general shape which are weakly neutral to multiple fields in two dimensions. We
show that a simply connected domain satisfying a certain geometric condition
can be realized as a weakly neutral inclusion to multiple fields by introducing
an imperfect interface parameter on the boundary. The geometric condition on
the domain and the imperfect interface parameter are determined by the first
coefficient of the conformal mapping from the exterior of the unit disk onto
the exterior of the domain. We provide some numerical examples to compare field
perturbations by weakly neutral inclusions and perfectly bonding interfaces.
They clearly show that the field perturbation by weakly neutral inclusions is
much weaker.

In the region between closetotouching hard inclusions, the stress may be
arbitrarily large as the inclusions get closer. The stress is represented by
the gradient of a solution to the Lam\'e system of linear elasticity. We
consider the problem of characterizing the gradient blowup of the solution in
the narrow region between two inclusions and estimating its magnitude. We
introduce singular functions which are constructed in terms of nuclei of strain
and hence are solutions of the Lam\'{e} system, and then show that the singular
behavior of the gradient in the narrow region can be precisely captured by
singular functions. As a consequence of the characterization, we are able to
regain the existing upper bound on the blowup rate of the gradient, namely,
$\epsilon^{1/2}$ where $\epsilon$ is the distance between two inclusions. We
then show that it is in fact an optimal bound by showing that there are cases
where $\epsilon^{1/2}$ is also a lower bound. This work is the first to
completely reveal the singular nature of the gradient blowup in the context of
the Lam\'{e} system with hard inclusions. The singular functions introduced in
this paper play essential roles to overcome the difficulties in applying the
methods of previous works. Main tools of this paper are the layer potential
techniques and the variational principle. The variational principle can be
applied because the singular functions of this paper are solutions of the
Lam\'{e} system.

This paper deals with the field enhancement, that is, the gradient blowup,
due to presence of a bowtie structure of perfectly conducting inclusions in
two dimensions. The bowtie structure consists of two disjoint bounded domains
which have corners with possibly different aperture angles. The domains are
parts of cones near the vertices, and they are nearly touching to each other.
We characterize the field enhancement using explicit functions and, as
consequences, derive optimal estimates of the gradient in terms of the distance
between two inclusions and aperture angles of the corners. The estimates show
that the field is enhanced beyond the corner singularities due to the
interaction between two inclusions.

We prove in a mathematically rigorous way the asymptotic formula of Flaherty
and Keller on the effective property of densely packed periodic elastic
composites with hard inclusions. The proof is based on the primaldual
variational principle, where the upper bound is derived by using the
Kellertype test functions and the lower bound by singular functions made of
nuclei of strain. Singular functions are solutions of the Lam\'e system and
capture precisely singular behavior of the stress in the narrow region between
two adjacent hard inclusions.

We prove that the elastic NeumannPoincar\'e operator defined on the smooth
boundary of a bounded domain in three dimensions, which is known to be
noncompact, is in fact polynomially compact. As a consequence, we prove that
the spectrum of the elastic NeumannPoincar\'e operator consists of three
nonempty sequences of eigenvalues accumulating to certain numbers determined
by Lam\'e parameters. These results are proved using the surface Riesz
transform, calculus of pseudodifferential operators and the spectral mapping
theorem.

We investigate anomalous localized resonance on the circular coated structure
and cloaking related to it in the context of elastostatic systems. The
structure consists of the circular core with constant Lam\'e parameters and the
circular shell of negative Lam\'e parameters proportional to those of the core.
We show that the eigenvalues of the NeumannPoincar\'e operator corresponding
to the structure converges to certain nonzero numbers determined by Lam\'e
parameters and derive precise asymptotics of the convergence. We then show with
estimates that cloaking by anomalous localized resonance takes place if and
only if the dipole type source lies inside critical radii determined by the
radii of the core and the shell.

We show that the eigenvalues of the NeumannPoincar\'e operator on analytic
boundaries of simply connected bounded planar domains tend to zero
exponentially fast, and the exponential convergence rate is determined by the
maximal Grauert radius of the boundary. We present a few examples of boundaries
to show that the estimate is optimal.

We study spectral properties of the NeumannPoincar\'e operator on planar
domains with corners with particular emphasis on existence of continuous
spectrum and pure point spectrum. We show that the rate of resonance at
continuous spectrum is different from that at eigenvalues, and then derive a
method to distinguish continuous spectrum from eigenvalues. We perform
computational experiments using the method to see whether continuous spectrum
and pure point spectrum appear on domains with corners. For the computations we
use a modification of the Nystr\"om method which makes it possible to construct
highorder convergent discretizations of the NeumannPoincar\'e operator on
domains with corners. The results of experiments show that all three possible
spectra, absolutely continuous spectrum, singularly continuous spectrum, and
pure point spectrum, may appear depending on domains. We also prove rigorously
two properties of spectrum which are suggested by numerical experiments:
symmetry of spectrum (including continuous spectrum), and existence of
eigenvalues on rectangles of high aspect ratio.

The neutral inclusion problem in two dimensional isotropic elasticity is
considered. The neutral inclusion, when inserted in a matrix having a uniform
applied field, does not disturb the field outside the inclusion. The inclusion
consists of the core and shell of arbitrary shapes, and their elasticity
tensors are isotropic. We show that if the coated inclusion is neutral to a
uniform bulk field, then the core and shell must be concentric disks, provided
that the shear and bulk moduli satisfy certain conditions.

We consider wellposedness of the boundary value problem in presence of an
inclusion with complex conductivity $k$. We first consider the transmission
problem in $\mathbb{R}^d$ and characterize solvability of the problem in terms
of the spectrum of the NeumannPoincar\'e operator. We then deal with the
boundary value problem and show that the solution is bounded in its $H^1$norm
uniformly in $k$ as long as $k$ is at some distance from a closed interval in
the negative real axis. We then show with an estimate that the solution depends
on $k$ in its $H^1$norm Lipschitz continuously. We finally show that the
boundary perturbation formula in presence of a diametrically small inclusion is
valid uniformly in $k$ away from the closed interval mentioned before. The
results for the single inclusion case are extended to the case when there are
multiple inclusions with different complex conductivities: We first obtain a
complete characterization of solvability when inclusions consist of two
disjoint disks and then prove solvability and uniform estimates when imaginary
parts of conductivities have the same signs. The results are obtained using the
spectral property of the associated NeumannPoincar\'e operator and the
spectral resolution.

We first investigate spectral properties of the NeumannPoincar\'e (NP)
operator for the Lam\'e system of elastostatics. We show that the
elastostatic NP operator can be symmetrized in the same way as that for
Laplace operator. We then show that even if elastostatic NP operator is not
compact even on smooth domains, its spectrum consists of eigenvalues which
accumulates to two numbers determined by Lam\'e constants. We then derive
explicitly eigenvalues and eigenfunctions on disks and ellipses. We then
investigate resonance occurring at eigenvalues and anomalous localized
resonance at accumulation points of eigenvalues. We show on ellipses that
cloaking by anomalous localized resonance takes place at accumulation points of
eigenvalues.

We study resonance for the Helmholz equation with a finite frequency in a
plasmonic material of negative dielectric constant in two and three dimensions.
We show that the quasistatic approximation is valid for diametrically small
inclusions. In fact, we quantitatively prove that if the diameter of a
inclusion is small compared to the loss parameter, then resonance occurs
exactly at eigenvalues of the NeumannPoincar\'e operator associated with the
inclusion.

The purpose of this paper is to propose a noniterative method for the
inverse conductivity problem of recovering multiple small anomalies from the
boundary measurements. When small anomalies are buried in a conducting object,
the electric potential values inside the object can be expressed by integrals
of densities with a common sparse support on the location of anomalies. Based
on this integral expression, we formulate the reconstruction problem of small
anomalies as a joint sparse recovery and present an efficient noniterative
recovery algorithm of small anomalies. Furthermore, we also provide a slightly
modified algorithm to reconstruct an extended anomaly. We validate the
effectiveness of the proposed algorithm over the linearized method and the
MUSIC algorithm by numerical simulations.

We study an elliptic equation with measurable coefficients arising from
photoacoustic imaging in inhomogeneous media. We establish Holder continuity
of weak solutions and obtain pointwise bounds for Green's functions subject to
Dirichlet or Neumann condition.

We consider the neutral inclusion problem in three dimensions which is to
prove if a coated structure consisting of a core and a shell is neutral to all
uniform fields, then the core and the shell must be concentric balls if the
matrix is isotropic and confocal ellipsoids if the matrix is anisotropic. We
first derive an overdetermined boundary value problem in the shell of the
neutral inclusion, and then prove in the isotropic case that if the over
determined problem admits a solution, then the core and the shell must be
concentric balls. As a consequence it is proved that the structure is neutral
to all uniform fields if and only if it consists of concentric balls provided
that the coefficient of the core is larger than that of the shell.

The purpose of this paper is to investigate the spectral nature of the
NeumannPoincar\'e operator on the intersecting disks, which is a domain with
the Lipschitz boundary. The complete spectral resolution of the operator is
derived, which shows in particular that it admits only the absolutely
continuous spectrum, no singularly continuous spectrum and no pure point
spectrum. We then quantitatively analyze using the spectral resolution the
plasmon resonance at the absolutely continuous spectrum.

We investigate in a quantitative way the plasmon resonance at eigenvalues and
the essential spectrum (the accumulation point of eigenvalues) of the
NeumannPoincar\'e operator on smooth domains. We first extend the
symmetrization principle so that the single layer potential becomes a unitary
operator from $H^{1/2}$ onto $H^{1/2}$. We then show that the resonance at the
essential spectrum is weaker than that at eigenvalues. It is shown that
anomalous localized resonance occurs at the essential spectrum on ellipses, but
cloaking does not occur on ellipses unlike the coreshell structure considered
in [20]. It is shown that resonance does not occur at the essential spectrum on
three dimensional balls.

If stiff inclusions are closely located, then the stress, which is the
gradient of the solution, may become arbitrarily large as the distance between
two inclusions tends to zero. In this paper we investigate the asymptotic
behavior of the stress concentration factor, which is the normalized magnitude
of the stress concentration, as the distance between two inclusions tends to
zero. For that purpose we show that the gradient of the solution to the case
when two inclusions are touching decays exponentially fast near the touching
point. We also prove a similar result when two inclusions are closely located
and there is no potential difference on boundaries of two inclusions. We then
use these facts to show that the stress concentration factor converges to a
certain integral of the solution to the touching case as the distance between
two inclusions tends to zero. We then present an efficient way to compute this
integral.

We present a new systematic method to construct the conformal mapping from
outside the unit disc to outside of a simply connected domain using the
generalized polarization tensors. We also present some numerical results to
validate effectiveness of the method.

The size estimation problem in electrical impedance tomography is considered
when the conductivity is a complex number and the body is twodimensional.
Upper and lower bounds on the volume fraction of the unknown inclusion embedded
in the body are derived in terms of two pairs of voltage and current data
measured on the boundary of the body. These bounds are derived using the
translation method. We also provide numerical examples to show that these
bounds are quite tight and stable under measurement noise.

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 farfield expansion of the voltage in the presence of the
conductivity inclusion. Relating the GPTs to the NeumanntoDirichlet (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 illposed inverse conductivity
problem. We demonstrate the viability of the proposed algorithm by preforming a
sensitivity analysis and giving some numerical examples.

If a core of dielectric material is coated by a plasmonic structure of
negative dielectric material with nonzero loss parameter, then anomalous
localized resonance may occur as the loss parameter tends to zero and the
source outside the structure can be cloaked. It has been proved that the
cloaking due to anomalous localized resonance (CALR) takes place for structures
of concentric disks and the critical radius inside which the sources are
cloaked has been computed. In this paper, it is proved that CALR takes place
for structures of confocal ellipses and the critical elliptic radii are
computed. The method of this paper uses the spectral analysis of the
NeumannPoincar\'e type operator associated with two interfaces (the boundaries
of the core and the shell).

When two perfectly conducting inclusions are located closely to each other,
the electric field concentrates in a narrow region in between two inclusions,
and becomes arbitrarily large as the distance between two inclusions tends to
zero. The purpose of this paper is to derive an asymptotic formula of the
concentration which completely characterizes the singular behavior of the
electric field, when inclusions are balls of the same radii in three
dimensions.

We consider the problem of neutral inclusions for twodimensional
conductivity and antiplane elasticity. The neutral inclusion, when inserted in
a matrix having a uniform field, does not disturb the field outside the
inclusion. The inclusion consists of a core and a shell. We show that if the
inclusion is neutral to two linearly independent fields, then the core and the
shell are confocal ellipses.