arxiv: v4 [math.sp] 11 Jan 2015

Transcription

1 Analysis of plasmon resonance on smooth domains using spectral properties of the Neumann-Poincaré operator Kazunori Ando Hyeonbae Kang arxiv: v4 [math.sp] Jan 205 January 24, 208 Abstract We investigate in a quantitative way the plasmon resonance at eigenvalues and the essential spectrum (the accumulation point of eigenvalues) of the Neumann-Poincaré operator on smooth domains. We first extend the symmetrization principle so that the single layer potential becomes a unitary operator from H /2 onto H /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 core-shell structure considered in [20]. It is shown that resonance does not occur at the essential spectrum on three dimensional balls. AMS subject classifications. 35J47 (primary), 35P05 (secondary) Key words. Neumann-Poincaré operator, symmetrization, plasmonic resonance, anomalous localized resonance, eigenvalues, essential spectrum Introduction The Neumann-Poincaré (NP) operator is a boundary integral operator which arises naturally when solving classical Dirichlet and Neumann boundary value problems using layer potentials. This operator can be realized as a self-adjoint operator using Plemelj s symmetrization principle (see the next section). If the boundary of the domain is smooth, the operator is compact (see [9, 6]) and its spectrum consists of the point spectrum (eigenvalues) and the essential spectrum which is an accumulation point of the eigenvalues. The purpose of this paper is to investigate resonance at eigenvalues and at the essential spectrum, and compare them in a quantitatively precise way. The resonance at the eigenvalues of the NP operator is the plasmon resonance [0]. We show that the resonance at the essential spectrum (on ellipses) is the anomalous localized resonance which was first discovered on a concentric core-shell structure in [22]. To be more precise, suppose that a bounded simply connected domain Ω in R d (d = 2, 3) is occupied with a plasmonic material of negative dielectric constant. In general the material property of the domain Ω is represented by ǫ c +iδ where ǫ c < 0 is the dielectric This work is supported by the Korean Ministry of Education, Sciences and Technology through NRF grants No Department of Mathematics, Inha University, Incheon , S. Korea (ando, hbkang@inha.ac.kr).

2 constant and δ > 0 indicates the dissipation. Let ǫ m > 0 be the dielectric constant of the matrix R d \Ω. So the distribution of the dielectric constant of the structure is given by { ǫ c +iδ in Ω, ǫ = in R d (.) \Ω. ǫ m We quantify the resonance through the following equation: { ǫ u = a δz in R d, u(x) = O( x d ) as x, (.2) where a is a constant vector and δ z is the Dirac mass at z R d \Ω. If u δ is the solution to (.2), the resonance is characterized by the blow-up of u δ L 2 (Ω): u δ L 2 (Ω) as δ 0. (.3) We are particularly interested in the blow-up rate of u δ L 2 (Ω) in terms of δ when the resonance takes place. We emphasize that as δ 0 ǫ tends to ǫ c < 0 in Ω, and so the problem is not an elliptic one. Resonance is connected to the spectrum of the NP operator because of the transmission conditions (continuity of the potential and the flux) along the interface Ω: u δ u δ + = u δ, ǫ m = (ǫ c +iδ) u δ on Ω, (.4) ν + ν where the subscripts + and denote limits to Ω from outside and inside Ω, respectively. We will make this connection clear in section 3. Thefindingsof thispapershowthat thegeneric rateof theresonanceat theeigenvalues is δ while that at the essential spectrum is weaker than δ. We then show that resonance occurs at the essential spectrum (0) on ellipses and exact rate of resonance is provided. It turns out that it is anomalous localized resonance (see subsection 5. for precise statements). We also show that resonance does not occur at the essential spectrum on the three dimensional balls. For the purpose of analysis on resonance we extend the symmetrization principle and show as a result that the single layer potential is a unitary operator from H /2 ( Ω) onto H /2 ( Ω) (H s is a Sobolev space). We also derive an expansion formula for the fundamental solution of the Laplacian in terms of eigenfunctions of the NP operator which seems of independent interest. It is worth mentioning that it is recently found in [2] that on a disk a complete resonance occurs at 0. This happens because 0 is the only eigenvalue of the NP operator on a disk. Recently there is rapidly growing interest in the spectral theory of the NP operator in relation to plasmonics. In [7] the Poincaré s variational program was revisited with modern prospective and symmetrization of the NP operator was proved. Eigenvalues on disks, ellipses, and balls were computed [4, 4]. Asymptotics of eigenvalues of the NP operator associated with closely located two-dimensional convex domain has been obtained in [5, 6]. There has been progress on spectral theory of the NP operator on non-smooth domains. A bound on the essential spectrum on the two-dimensional curvilinear domains has been obtained [23]. Quite recently, the complete spectral resolution of the NP operator on intersecting disks has been derived which in particular shows that only the absolutely continuous spectrum exists [3]. Interestingly, the results in [3] shows that the spectral bound obtained in [23] is optimal in the case of intersecting disks. The spectral theory of the NP operator has been applied to analysis of cloaking by anomalous localized resonance on the plasmonic structure [] and to the study of uniformity of elliptic estimates [2]. 2

3 This paper is organized as follows. In section 2 we review and extend the symmetrization of the NP operator. In section 3 we derive an expansion formula for the fundamental solution of the Laplacian in terms of the eigenfunctions of the NP operator and obtain a representation of the solution to (.2). Resonance at eigenvalues and at essential spectrum is studied in section 4 and section 5, respectively. In Appendix we derive explicit expansion formula for the fundamental solution in terms of spherical harmonics and elliptic harmonics. 2 Neumann-Poincaré operator and symmetrization Let us first fix some notation. We denotebyh /2 ( Ω) thedualspaceof H /2 ( Ω), andby, theduality pairing of H /2 and H /2, and s denotes the H s norm on Ω. Let H /2 0 ( Ω) be the space of ψ H /2 ( Ω) satisfying ψ, = 0. The notation A B means that A CB for some constant C, and A B means that both A B and B A hold. Let f(δ) and g(δ) be positive quantities depending on δ. We write f(δ) g(δ) as δ 0 if there are constants C and C 2 such that C < f(δ) g(δ) < C 2. (2.) Let Γ(x) be the fundamental solution to the Laplacian on R d (d = 2,3), i.e., ln x, d = 2, Γ(x) = 2π 4π x, d = 3. The Neumann-Poincaré (NP) operator on Ω, denoted by K Ω, is defined by K Ω [ϕ](x) := Γ(x y)ϕ(y)dσ(y), x Ω. (2.2) ν x Ω where ν x indicates the outward normal derivative in x-variable. The adjoint operator K Ω on L2 ( Ω) will be called the adjoint NP operator. Importance of the adjoint operator in dealing with interface problems lies in its relation with the single layer potential S Ω [ϕ] defined by S Ω [ϕ](x) := Γ(x y)ϕ(y)dσ(y), x R d. (2.3) The following jump relation holds: Ω ν S Ω[ϕ] ( ± (x) = ± ) 2 I +K Ω [ϕ](x), x Ω, (2.4) where the subscripts ± indicate the limits (to Ω) from outside and inside of Ω, respectively. 3

4 The NP operator is not self-adjoint on the usual L 2 -space, unless the domain Ω is a disc or a ball (see [9]). However, it is found in [7] (see also []) that the adjoint NP operator K Ω can be symmetrized using Plemelj s symmetrization principle (also known as Calderón s identity): S Ω K Ω = K ΩS Ω. (2.5) Define, for ϕ,ψ H /2 ( Ω), ϕ,ψ H := ϕ,s Ω [ψ]. (2.6) Since S Ω maps H /2 ( Ω) into H /2 ( Ω), the right hand side of (2.6) is well-defined. It is known that, H is an inner product on H /2 0 ( Ω), and the norm H induced by this inner product is equivalent to the H /2 ( Ω) norm, namely, ϕ H ϕ /2 (2.7) for all ϕ H /2 0 ( Ω) (see [2]). Let H0 be the space H /2 0 ( Ω) equipped with the inner product, H. Then the symmetrization principle (2.5) shows that K Ω is self-adjoint on H0. Letusnowconsidersymmetrization ofk Ω. TheNPoperatorK Ω canbesymmetrized using (2.5) expressed in a different form: K Ω S Ω = S Ω K Ω, (2.8) provided that S Ω exists. However, it was proved in [24] that S Ω exists only in three dimensions (or higher), and there are domains Ω in two dimensions where S Ω does not exist. Here we present a simple way to overcome this difficulty. To symmetrize K Ω, we use, as a replacement of S Ω in (2.5), the operator A : H /2 ( Ω) C H /2 ( Ω) C defined by A(ψ,a) := (S Ω [ψ]+a, ψ, ). (2.9) It is proved in [3, Theorem 2.3] and [2] that A is invertible for d 2. For f H /2 ( Ω) let (ψ f,a f ) = A (f,0). (2.0) Then ψ f H /2 0 ( Ω) and it holds that Moreover, we have S Ω [ψ f ]+a f = f (2.) ψ f /2 + a f f /2. (2.2) This shows in particular that the mapping f a is a bounded linear functional on H /2 ( Ω). So, there is a unique ϕ 0 H /2 ( Ω) such that a f = ϕ 0,f (2.3) for all f H /2 ( Ω). Since such a (ψ f,a f ) is unique, we see that a f = if f =, in other words, ϕ 0, =. (2.4) 4

5 Since K Ω [] = /2, by applying K Ω to both sides of (2.) we see that K Ω S Ω [ψ f ]+ 2 a f = K Ω [f]. It then follows from (2.5) that S Ω K Ω [ψ f]+ 2 a f = K Ω [f], which implies that a K Ω [f] = 2 a f, in other words, ϕ 0,K Ω [f] = ϕ 0,f /2. So, we have K Ω [ϕ 0] = 2 ϕ 0, (2.5) namely, ϕ 0 is an eigenfunction of K Ω (on H /2 ( Ω) corresponding to the eigenvalue /2 normalized by (2.4). We can infer from the jump formula (2.4) and (2.5) that S Ω [ϕ 0 ] is constant in Ω. We emphasize that the function ϕ 0 already appeared in literatures. It is proved in [24] that S Ω [ϕ 0 ] = 0 in Ω for some domain Ω in two dimensions, which is why S Ω is not invertible. We now define a variant of the single layer potential on H by { S Ω [ϕ] if ϕ, = 0, S Ω [ϕ] = (2.6) if ϕ = ϕ 0. Then we have an extension of (2.5): S Ω K Ω = K Ω S Ω. (2.7) Using this we can extend the inner product (2.6) (defined on H /2 0 ( Ω)) to H /2 ( Ω). For ϕ,ψ H /2 ( Ω), define We also have a new inner product on H /2 ( Ω): ϕ,ψ H := ϕ, S Ω [ψ]. (2.8) ϕ,ψ H := ϕ, S Ω [ψ]. (2.9) If we define H to be H /2 ( Ω) with this inner product, then K Ω is self-adjoint on H. We emphasize that norm H is equivalent to /2. Moreover, S Ω : H H is a unitary operator. Observe that if {ψ j } is an orthonormal basis of H 0, then {ψ j} {ϕ 0} is an orthonormal basis of H and {S Ω [ψ j ]} {} is an orthonormal basis of H. 3 Representation of the solution If ΩisC,α -smoothforsomeα > 0, thenitisknownthatk Ω isacompact operatoronh (see [6]). SinceK Ω is self-adjoint on H 0, its eigenvalues {λ j} are real andaccumulates to 0. We emphasize that λ j < /2 (see [7, 24]). Let {ψ j } with ψ j H = be the corresponding (real valued) eigenfunctions counting the multiplicities. Then we have shown in the previous section that {S Ω [ψ j ]} {} is an orthonormal basis of H. 5

6 Fix z R d \ Ω. Then Γ( z) belongs to H /2 ( Ω), and so admits the following decomposition: Γ(x z) = c j (z)s Ω [ψ j ](x)+c 0 (z), x Ω, (3.) for some constants c j (z) (depending on z) satisfying c j (z) 2 <. (3.2) Since S Ω [ψ j ],ψ i = δ ij where δ ij is the Kronecker s delta, we see that c j (z) = Γ(x z)ψ j (x)dσ(x) = S Ω [ψ j ](z), j =,2,3,... Ω We also see from (2.4) that So, we obtain the following formula: Observe that Γ(x z) = c 0 (z) = S Ω [ϕ 0 ](z), j =,2,3,... S Ω [ψ j ](z)s Ω [ψ j ](x)+s Ω [ϕ 0 ](z), x Ω, (3.3) S Ω [ψ j ](z)s Ω [ψ j ] 2 H = S Ω [ψ j ](z) 2 <. (3.4) Since H is equivalent to /2 as shown in the previous section, we find from the trace theorem that S Ω[ψ j ](z)s Ω [ψ j ] converges in H (Ω) and is harmonic in Ω. So, we obtain the following theorem on an expansion of the fundamental solution to the Laplacian. Theorem 3. It holds that Γ(x z) = S Ω [ψ j ](z)s Ω [ψ j ](x)+s Ω [ϕ 0 ](z), x Ω, z R d \Ω. (3.5) If Ω is a ball, then S Ω [ψ j ] is a spherical harmonics (see subsection 5.2). Therefore (3.5) is the expansion of Γ(x z) in terms of the spherical harmonics, which is well-known (see, for example, [6]). We will show explicit formula for (3.5) when Ω is a unit ball in R 3 or an ellipse in R 2 in the Appendix. The formula (3.5) shows that the fundamental solution can be expanded on any smooth domain in terms of the localized plasmon (the single layer potential of an eigenfunction is called a localized plasmon). The estimate (3.4) shows that S Ω [ψ j ](z) 0 as j for all z / Ω. So, if j is large then S Ω [ψ j ] is localized near Ω. It explains why S Ω [ψ j ] is called a localized plasmon. We now derive a representation of the solution to (.2). Let F z (x) := a x Γ(x z), x z. (3.6) 6

7 Then F z (x) = a δ z (x), and hence we see from transmission conditions (.4) on Ω and (2.4) that the solution u δ to (.2) can be represented as u δ (x) = F z (x)+s Ω [ϕ δ ](x), x R d, (3.7) where the potential ϕ δ H0 ( Ω) is the solution to the integral equation ( ν F z denotes the normal derivative of F z ). Here (λi K Ω )[ϕ δ] = ν F z on Ω (3.8) λ := ǫ c +ǫ m +iδ 2(ǫ c ǫ m )+2iδ. (3.9) ǫ c+ǫ m 2(ǫ c ǫ m) Note that λ ǫc+ǫm 2(ǫ as δ 0. The number ǫ c ǫ m) c/ǫ m such that is an eigenvalue of K Ω is called a plasmonic eigenvalue [0]. We emphasize that ǫ c is negative if and only ǫ if c+ǫ m 2(ǫ lies in ( /2,/2) where the spectrum of c ǫ m) K Ω lies. Since K Ω admits the spectral decomposition K Ω = λ j ψ j ψ j, (3.0) the solution ϕ δ to (3.8) can be expressed as where We can see from (3.6) that α j (z) = a It can seen from (2.4) and (3.5) that ν x Γ(x z) = ϕ δ = α j (z) λ λ j ψ j, (3.) α j (z) := ν F z,ψ j H. (3.2) Ω ν x Γ(x z)s Ω [ψ j ](x)dσ(x). S Ω [ψ j ](z) ν S Ω[ψ j ](x) = ( 2 λ j) S Ω [ψ j ](z)ψ j (x). It then follows that α j (z) = ( λ j ) a S Ω [ψ j ](z). (3.3) 2 4 Resonance at eigenvalues We investigate the behavior of the solution u δ when λ approaches to one of eigenvalues λ l 0 as δ 0, namely, when ǫ c +ǫ m 2(ǫ c ǫ m ) = λ l. (4.) We show that u δ L 2 (Ω) δ as δ 0, (4.2) 7

8 as one may expect. We first show that S Ω [ϕ] 2 L 2 (Ω) ϕ 2 H (4.3) for all ϕ H0. In fact, we see from (2.4) and (3.0) that S Ω [ϕ] 2 L 2 (Ω) = S Ω [ϕ] Ω ν S Ω[ϕ] dσ = ϕ,( 2 I +K Ω )[ϕ] = H ( 2 λ j) ϕ,ψ j H 2. Since λ j < /2 and accumulates to 0, we have (4.3). We now see from (3.) that (u δ F z ) 2 L 2 (Ω) ϕ δ 2 H = α j (z) 2 λ λ l 2 + α j (z) 2 λ λ j 2. λ j =λ l λ j λ l If λ j λ l, then λ j λ l > C for some positive constant C since λ l 0, and so the second term on the righthand side above is bounded as δ 0. Since λ λ l δ as δ 0, we obtain (u δ F z ) L 2 (Ω) δ. Since F z L 2 (Ω) is bounded, (4.2) follows. 5 Resonance at the essential spectrum In this section we assume that 0 is not an eigenvalue of K Ω. Since eigenvalues of K Ω converges to 0, {0} is the essential spectrum of K Ω. It is worth mentioning that we are not aware of any domain other than diskson which thenp operator has 0 as an eigenvalue. If Ω is a disk, then K Ω 0 on H 0. We consider the resonance when λ 0 as δ 0, in other words, when ǫ c +ǫ m = 0. (5.) In this case, we assume that λ = iδ for simplicity of presentation. We first obtain the following theorem which shows that the resonance, if it occurs, is never at the rate of δ unlike the resonance at eigenvalues. Theorem 5. It holds that Proof. We obtain from (3.) and (4.3) that lim δ u δ L δ 0 2 (Ω) = 0. (5.2) (u δ F z ) 2 L 2 (Ω) α j (z) 2 δ 2 +λ 2. (5.3) j 8

9 We then decompose the summation into two parts as α j (z) 2 δ 2 +λ 2 j = λ j δ α j (z) 2 δ 2 +λ 2 j + λ j >δ α j (z) 2 δ 2 +λ 2 j =: S +S 2. (5.4) Since α j (z) 2 <, we have δ 2 S λ j δ To show that δ 2 S 2 0, we express S 2 as k=0 S 2 = Then we see that δ 2 S k α j (z) 2 0 as δ 0. k=02 k δ< λ j 2 k+ δ 2 k δ< λ j 2 k+ δ α j (z) 2 α j (z) 2 δ 2 +λ 2. j k=0 +2 2k λ j 2 k+ δ α j (z) 2. For each fixed k, λ j 2 k+ δ α j(z) 2 0 as δ 0. So we infer that S 2 0 as δ 0 by Lebesgue dominated convergence theorem. The proof is complete. In order to derive estimates for the resonance, we need to investigate the asymptotic behavior of the quantity on the righthand side of (5.3) as δ 0. This seems a quite difficult task since it depends on the behavior of λ j and α j (z) as j. So, we deal with two specific domains, ellipses and three dimensional balls, where eigenvalues and eigenfunctions are known. It is worth mentioning that the disks are out of consideration since the NP operator on disks are Anomalous localized resonance on ellipses We first consider the resonance when Ω is an ellipse in R 2. The elliptic coordinates x = (x,x 2 ) = (x (ρ,ω),x 2 (ρ,ω)) is given by x (ρ,ω) = Rcosωcoshρ, x 2 (ρ,ω) = Rsinωsinhρ, ρ > 0, 0 ω < 2π. When the above holds, we denote ρ = ρ x and ω x. Then we can represent Ω as Ω = {x R 2 ;ρ x = ρ 0 } (5.5) for some ρ 0 > 0. The number ρ 0 is called the elliptic radius of Ω. The position z of the source is denote by ρ z and ω z in elliptic coordinates. We obtain the following theorem whose proof will be given in the next section. Theorem 5.2 Suppose that Ω is an ellipse given by (5.5). Then we have δ 3+ρz/ρ 0 logδ if ρ 0 < ρ z < 3ρ 0, u δ 2 L 2 (Ω) logδ 2 if ρ z = 3ρ 0, if ρ z > 3ρ 0, (5.6) as δ 0. 9

10 The quantity E δ := δ u δ 2 L 2 (Ω) is of particular interest since it represents the imaginary part of the energy, namely, δ u δ 2 L 2 (Ω) = I R d ǫ u δ 2 dx. Physically it represents the electro-magnetic energy dissipated into heat. Estimates (5.6) shows that E δ if ρ 0 < ρ z 2ρ 0 while it tends to 0 as δ 0 if ρ z > 2ρ 0. So the critical (elliptic) radius is 2ρ 0. This phenomenon is reminiscent of the anomalous localized resonance (ALR) discovered in [22, 20] (see also [, 8, 8]). There it is shown that a disk (a core) is coated by a concentric disk (a shell) of plasmonic material, then ALR occurs if the source is located within a critical radius, and does not occur for sources outside the radius. This is exactly what (5.6) shows. It also in accordance with the result in [8] where the core-shell structure of confocal ellipses of radii ρ e > ρ i is considered. There it is shown that the critical (elliptic) radius is given by { (3ρ e ρ i )/2, if ρ e 3ρ i, ρ = 2(ρ e ρ i ), if ρ e > 3ρ i. If the core shrinks, in other words ρ i tends to 0, then ρ tends to 2ρ e which is the critical radius found in this paper. ALR also requires the solution to be bounded outside a bounded set. This requirement is satisfied on ellipse as the following theorem shows. Theorem 5.2 and Theorem 5.3 show that ALR may occur not only on the coated structures but also on simply connected structure. So ALR may be regarded as resonance at the accumulation points of eigenvalues. We emphasize that the NP operator on the coated disk has 0 as its essential spectrum as proved in []. Theorem 5.3 Let Ω be an ellipse given by (5.5). It holds for all x satisfying ρ x +ρ z 4ρ 0 > 0 that u δ (x) F z (x) e n(ρx+ρz 4ρ0). (5.7) In particular, let ρ > 0 be such that ρ > 4ρ 0 ρ z, then there exists some C = C ρ > 0 such that sup u δ (x) F z (x) < C. (5.8) ρ x ρ To prove Theorem 5.2 and Theorem 5.3 let us recall some facts on the NP operator on ellipses. It is proved in [4] and [8] that eigenvalues of K Ω are and corresponding eigenfunctions are λ n = 2e 2nρ, n =,2,, (5.9) 0 φ c n(ω) := Ξ(ρ 0,ω) cosnω, φ s n(ω) := Ξ(ρ 0,ω) sinnω, n =,2,, (5.0) where Ξ = Ξ(ρ 0,ω) := R sinh 2 ρ 0 +sin 2 ω. (5.) 0

11 It is also proved in the same papers that coshnρ e nρ 0 cosnω, ρ < ρ S Ω [φ c 0, n](x) = n coshnρ 0 e nρ cosnω, ρ ρ 0, n and sinhnρ e nρ 0 sinnω, ρ < ρ S Ω [φ s 0, n](x) = n sinhnρ 0 e nρ sinnω, ρ ρ 0. n (5.2) (5.3) Here and afterwards (ρ,ω) is the elliptic coordinates of x. Since the length element dσ on Ω is given by dσ = Ξdω, (5.4) one can see that and So, the normalized eigenfunctions are φ c n 2 H = φc n,s Ω[φ c n ] = πcoshnρ 0 ne nρ, 0 φ s n 2 H = πsinhnρ 0. ne nρ 0 ψ c n := ne nρ 0 πcoshnρ 0 φ c n, ψs n := ne nρ 0 πsinhnρ 0 φ s n. (5.5) We see from (5.2) e nρ 0 coshnρcosnω, ρ < ρ 0, S Ω [ψn](x) c = nπcoshnρ 0 e nρ 0 coshnρ 0 e nρ cosnω, ρ ρ 0, nπ and from (5.3) e nρ 0 sinhnρsinnω, ρ < ρ 0, S Ω [ψn](x) s = nπsinhnρ 0 e nρ 0 sinhnρ 0 e nρ sinnω, ρ ρ 0. nπ (5.6) (5.7) Proof of Theorem 5.2. We see from (3.3) and (5.3) that = = (u δ F z ) 2 L 2 (Ω) δ 2 +λ 2 n δ 2 +λ 2 n α j (z) 2 δ 2 +λ 2 j [ a S Ω [ψ c n ](z) 2 + a S Ω [ψ s n ](z) 2] [ e nρ 0 coshnρ 0 nπ ( a z e nρ z ) cosnω z 2 e nρ 0 sinhnρ 0 ( + a z e nρ z ) ] sinnω z 2. nπ

12 Since coshnρ 0 sinhnρ 0 e nρ 0, we see that α j (z) 2 δ 2 +λ 2 j δ 2 +λ 2 n e2nρ 0[ ( a z e nρ z ) cosnω z 2 ( + a z e nρ z ) ] sinnω z 2 n as δ 0. Let U(ω) be the rotation by the angle ω, namely, [ ] cosω sinω U(ω) =. sinω cosω Using the change of variables formula and x = x 2 = R Ξ(ρ,ω) 2 ( cosωsinhρ ρ sinωcoshρ ω ( R Ξ(ρ,ω) 2 sinωcoshρ ρ +cosωsinhρ ) ω ) (5.8) where Ξ(ρ, ω) is given by (5.), we see through tedious but straightforward computations (we omit the computations) that and where a ( e nρ cosnω ) = Rne nρ Ξ(ρ,ω) 2 a U(nω)b(ρ,ω) a ( e nρ sinnω ) = Rne nρ a U(nω π/2)b(ρ,ω), Ξ(ρ,ω) 2 b(ρ,ω) = [ ] cosωsinhρ. sinωcoshρ Let θ n be the angle between the vectors a and U(nω)b(ρ,ω). Then we have a U(nω)b(ρ,ω) 2 = a 2 b(ρ,ω) 2 cos 2 θ n and which implies a U(nω π/2)b(ρ,ω) 2 = a 2 b(ρ,ω) 2 sin 2 θ n, a ( e nρ cosnω ) 2 + a ( e nρ sinnω ) 2 = n2 e 2nρ a 2 b(ρ,ω) 2 R 2( sinh 2 ρ+sin 2 ω ) 2. (5.9) It then follows from (5.8) that α j (z) 2 δ 2 +λ 2 j ne 2nρ 0 e 2nρz δ 2 +λ 2 n (5.20) as δ 0. 2

13 We now investigate the asymptotic behavior of the righthand side of (5.20) as δ 0. Let [ N = ] log2δ, (5.2) 2ρ 0 which is the first number such that δ > 2 e 2Nρ 0. Then one can easily see that ne 2nρ 0 e 2nρz δ 2 +λ 2 = + ne 2nρ 0 e 2nρz n e 4nρ + 0 δ 2 ne 2n(ρz ρ0). n N n>n n N Observe that n N ne 2nρ 0 e 2nρz e 4nρ 0 On the other hand, we have So we infer that n>n n Nne if ρ z > 3ρ 0, 2n(3ρ 0 ρ z) logδ 2 if ρ z = 3ρ 0, logδ δ 3+ρz/ρ 0 if ρ 0 < ρ z < 3ρ 0. δ 2 ne 2n(ρz ρ0) logδ δ 3+ρz/ρ 0. n>n ne 2nρ 0 e 2nρz δ 2 +λ 2 n Since F z 2 L 2 (Ω) is bounded, we obtain (5.6). if ρ z > 3ρ 0, logδ 2 if ρ z = 3ρ 0, logδ δ 3+ρz/ρ 0 if ρ 0 < ρ z < 3ρ 0. Proof of Theorem 5.3. One can see from (3.7), (3.) and (3.3) that ( ) [ u δ (x) F z (x) = iδ λ n 2 λ n (a z S Ω[ψn c ](z))s Ω[ψn c ](x) +(a z S Ω [ψ s n ](z))s Ω[ψ s n ](x) ]. It then follows from (5.6) and (5.7) that ( ) u δ (x) F z (x) = iδ λ n 2 λ n [ e nρ 0 coshnρ 0 ( ( a z e nρ z )) cosnω z e nρ x cosnω x nπ + enρ 0 sinhnρ 0 ( ( a z e nρ z )) ] sinnω z e nρ x sinnω x nπ where (ρ z,ω z ) is the elliptic coordinates of z. Therefore we have e 4nρ 0[ ( u δ (x) F z (x) a z e nρ z ) ( cosnω z + a z e nρ z ) ] sinnω z e nρx n We then see from (5.9) that u δ (x) F z (x) e 4nρ 0 n ne nρz e nρx = e n(ρx+ρz 4ρ0). (5.8) is an immediate consequence of (5.7). This completes the proof. 3

14 5.2 Anomalous localized resonance on balls The following theorem shows that resonance at the essential spectrum does not occur if Ω is a three dimensional ball. Theorem 5.4 Suppose that Ω is a three dimensional ball and f(x) = a δ z (x) for some z R 3 \Ω. Then there is a constant C such that for all δ. u δ L 2 (Ω) C (5.22) Proof. Assume that Ω is the unit ball centered at 0 for convenience. Then eigenvalue of K Ω are λ n =, n =,2,..., (5.23) 2(2n+) and corresponding eigenfunctions are Yn m (ˆx), m = n, n+,,n,n, the spherical harmonics of degree n, where ˆx = x/ x (see [2]). It is also proven in the same paper that S Ω [Yn m ](x) = 2n+ rn Yn m (ˆx), for x = r <, (5.24) 2n+ r (n+) Yn m (ˆx), for x = r. So the eigenfunctions normalized in H are ψ m n (x) = According to (3.3) and (5.3), we have 2n+ Y m n (x), m = n, n+,,n,n. (5.25) α j (z) 2 δ 2 +λ 2 j δ 2 +λ 2 n n m= n a S Ω [ψ m n ](z) 2. (5.26) We seefrom (5.24) thats Ω [ψ m n ] is ahomogeneous harmonicpolynomial of degree (n+) in R 3 \Ω. So, we have x j S Ω [ψ m n ](x) = r (n+2) n+ m= n for some constants a j nm, j =,2,3. We then have R 3 \Ω We then infer from (4.3) that S Ω [ψ m n ](x) 2 dx = 2(n+2) 3 n+ m= n 2(n+2) a j nm Y n+ m (ˆx), x > (5.27) 3 n+ m= n a j nm 2 ψ m n H =. a j nm 2. 4

15 It then follows from (5.27) that S Ω [ψn m ](z) ( 3 2 = r 2(n+2) 0 (n+)r 2(n+2) 0 where r 0 = z. Using the Unsöld s theorem we obtain n m= n We now see from (5.26) that n+ m= n n+ m= n a j nm 2)( n+ Y m n (ẑ) 2, m= n Y m n (ẑ) 2 = 2n+ 4π, n = 0,,2,, Y m n (ẑ) 2) S Ω [ψ m n ](z) 2 n 2 r 2(n+2) 0. (5.28) (u δ F z ) 2 L 2 (Ω) n 3 r 2(n+2) 0 δ 2 +λ 2 n n 5 r 2(n+2) 0 for any δ. Since F z 2 L 2 (Ω) is finite, the proof is completed. Appendix: Expansions of the fundamental solution In this appendix we write down the formula (3.5) explicitly when Ω is a unit ball in R 3 or an ellipse in R 2 to demonstrate the connection of (3.5) with the known expansion formula of the fundamental solution for the Laplace equation. Let Ω be the unit ball in R 3. We see from (5.25) that the eigenfunctions in H 0 are { Yn m 2n+ ; n =,2,... and m = n, n+,...,n,n,}. One can easily see that ϕ 0 (x) = 4π and, for x = r <, S Ω [ϕ 0 ](x) = 4π (5.29) 4π r, for x = r. Then, using (5.24) and (5.29), we have Γ(x z) = r n r n z (2n+) 3 n m= n Y n,m (ˆx)Y n,m (ẑ), r < r z, (5.30) 4π r z where r = x, r z = z and ˆx = x/ x, ẑ = z/ z. As mentioned right after Theorem 3., (5.30) is the well-known expansion of the fundamental solution in R 3 in terms of spherical harmonics. 5

16 Now let Ω be an ellipse with the elliptic radius ρ 0 in R 2. Then eigenfunctions in H 0 are given by {ψ c n,ψ s n} defined by (5.5). We can see from (5.6) and (5.7) that Ξ(ρ 0,ω) is orthogonal to {ψ c n,ψ s n} in H. So we have from (2.4) that ϕ 0 (ω) = 2πΞ(ρ 0,ω). Moreover, one can see that S Ω [ϕ 0 ](x) = 2π (ρ 0 +lnr ln2), ρ < ρ 0, 2π (ρ+lnr ln2), ρ ρ 0, where (ρ,ω) be the elliptic coordinates of x R 2. In fact, as ρ, we have x = R cosh 2 ρcos 2 ω +sinh 2 ρsin 2 ω = R 2 eρ +O(), (5.3) which implies that ρ+lnr ln2 = ln x +O( x ) as x. Note that ρ is a harmonic function. So, the righthand side of (5.3) is constant in ρ < ρ 0, harmonic in ρ > ρ 0, and behaves like ln x as x. So the equality in (5.3) holds. Using (5.6), (5.7) and (5.3), we now have Γ(x z) = ( coshnρcosnρe nρ z ) cosnρ z +sinhnρsinnρe nρz sinnρ z nπ + ( ( )) R ρ z +ln, ρ ρ 0 < ρ z, 2π 2 (5.32) where (ρ,ω) is the elliptic coordinate of x Ω and (ρ z,ω z ) is that of z R 2 \ Ω. An expansion of a derivative of Γ(x z) is also obtained in [8, (4.6)]. References [] H. Ammari, G. Ciraolo, H. Kang, H. Lee and G.W. Milton, Spectral theory of a Neumann-Poincaré-type operator and analysis of cloaking due to anomalous localized resonance, Arch. Ration. Mech. An. 208 (203), [2] H. Ammari, G. Ciraolo, H. Kang, H. Lee and G.W. Milton, Spectral theory of a Neumann-Poincaré-type operator and analysis of anomalous localized resonance II, Contemporary Math. 65, 4, 204. [3] H. Ammari and H. Kang, Polarization and moment tensors, Applied Mathematical Sciences, 62, Springer, New York [4] H. Ammari, H. Kang, and H. Lee, A boundary integral method for computing elastic moment tensors for ellipses and ellipsoids, J. Comp. Math. 25 () (2007), 2-2. [5] E. Bonnetier and F. Triki, Pointwise bounds on the gradient and the spectrum of the Neumann-Poincaré operator: The case of 2 discs, Contemp. Math. 577 (202),

17 [6] E. Bonnetier and F. Triki, On the spectrum of Poincaré variational problem for two close-to-touching inclusions in 2D, Arch. Ration. Mech. An. 209 (203), [7] T. Chang and K. Lee, Spectral properties of the layer potentials on Lipschitz domains, Illinois J. Math. 52 (2008), [8] D. Chung, H. Kang, K. Kim and H. Lee, Cloaking due to anomalous localized resonance in plasmonic structures of confocal ellipses, SIAM J. Appl. Math. 74 (204), [9] G. B. Folland, Introduction to partial differential equations, 2nd Ed., Princeton Univ. Press, Princeton, 995. [0] D. Grieser, The plasmonic eigenvalue problem, Rev. Math. Phys. 26 (204), [] H. Kang, Layer potential approaches to interface problems, Inverse problems and imaging, Panoramas et Syntheses, Societe Mathematique de France, to appear. [2] H. Kang, K. Kim, H. Lee and J. Shin, Spectral properties of the Neumann-Poincaré operator and uniformity of estimates for the conductivity equation with complex coefficients, submitted, arxiv [3] H. Kang, M. Lim and S. Yu, Spectral resolution of the Neumann-Poincaré operator on intersecting disks and analysis of plasmonic resonance, preprint. [4] H. Kang and J.-K. Seo, The inverse conductivity problem with one measurement: uniqueness for balls, SIAM J. of Applied Math. Vol 59 (999), [5] H. Kang and J-K. Seo, Recent progress in the inverse conductivity problem with single measurement, in Inverse Problems and Related Fields, CRC Press, Boca Raton, FL, 2000, [6] O. D. Kellogg, Foundations of potential theory, Dover, New York, 953 (Reprint from the first edition of Die Grundlehren der Mathematischen Wissenschaften, Band 3, Springer-Verlag, Berlin-New York, 929). [7] D. Khavinson, M. Putinar, and H. S. Shapiro, Poincaré s variational problem in potential theory, Arch. Ration. Mech. An. 85 (2007), [8] R. V. Kohn, J. Lu, B. Schweizer and M. I. Weinstein, A variational perspective on cloaking by anomalous localized resonance, Comm. Math. Phys. 328 (204), 27. [9] M. Lim, Symmetry of a boundary integral operator and a characterization of a ball, Illinois J. Math. 45 (200), [20] G.W. Milton and N.-A.P. Nicorovici, On the cloaking effects associated with anomalous localized resonance, Proc. R. Soc. A 462 (2006), [2] H.M. Nguyen and L.H. Nguyen, Localized and complete resonance in plasmonic structres, arxiv: [22] N.-A.P. Nicorovici, R.C. McPhedran, and G.W. Milton, Optical and dielectric properties of partially resonant composites, Phys. Rev. B 49 (994),

18 [23] K. Perfekt and M. Putinar, Spectral bounds for the Neumann-Poincaré operator on planar domains with corners, arxiv: [24] G.C. Verchota, Layer potentials and boundary value problems for Laplace s equation in Lipschitz domains, J. Funct. Anal. 59 (984),