On Quasi-static Cloaking Due to Anomalous Localized Resonance in \mathbb{R}^{3}

On Quasi-static Cloaking Due to Anomalous Localized Resonance in

Hongjie Li School of Mathematics and Statistics, Beijing Institute of Technology, Beijing, 100081, P. R. China.    Jingzhi Li Department of Mathematics, South University of Science and Technology of China, Shenzhen 518055, P. R. China. Email: li.jz@sustc.edu.cn    Hongyu Liu Department of Mathematics, Hong Kong Baptist University, Kowloon Tong, Hong Kong SAR, and HKBU Institute of Research and Continuing Education, Virtual University Park, Shenzhen, P. R. China. Email: hongyu.liuip@gmail.com
Abstract

This work concerns the cloaking due to anomalous localized resonance (CALR) in the quasi-static regime. We extend the related two-dimensional studies in [2, 11] to the three-dimensional setting. CALR is shown not to take place for the plasmonic configuration considered in [2, 11] in the three-dimensional case. We give two different constructions which ensure the occurrence of CALR. There may be no core or an arbitrary shape core for the cloaking device. If there is a core, then the dielectric distribution inside it could be arbitrary.

Key words. anomalous localized resonance, plasmonic material, invisibility cloaking

AMS subject classifications. 35R30, 35B30

1 Introduction

This work concerns the invisibility cloaking due to anomalous localized resonance (CALR) in the quasi-static regime, which has gained growing interest in the literature; see [1, 2, 3, 5, 6, 7, 8, 10, 11, 12, 13] and the references therein. Let and be bounded domains in , , such that . , and signify, respectively, the core, shell and matrix of a cloaking device, which hosts a dielectric object as follows

\hb@xt@.01(1.1)

In most of the existing studies, one takes and . is negatively valued, which denotes the plasmonic material parameter. Let denote the loss parameter and consider a material distribution given as

\hb@xt@.01(1.2)

where is given in (LABEL:eq:struc), and denotes the characteristic function of the domain , with or . For time-harmonic wave propagation in the quasi-static regime, the wave pressure satisfies the following equation

\hb@xt@.01(1.3)

where denotes a source term that is compactly supported in and satisfies

\hb@xt@.01(1.4)

Define

\hb@xt@.01(1.5)

where is the solution to (LABEL:eq:m1). denotes the rate at which the energy of the wave field is dissipated into heat. Then anomalous localized resonance (ALR) is said to occur if there holds

\hb@xt@.01(1.6)

In what follows, we sometimes simply refer to ALR as resonance. If in addition to (LABEL:eq:cond1), one further has that

\hb@xt@.01(1.7)

where is such that the central ball contains , then it is said that CALR occurs. Here and also in what follows, with denotes a central ball of radius in , . By (LABEL:eq:cond2), it is readily seen that if CALR occurs, then both the source and the cloaking device are invisible to the wave observation made from the outside of . If (LABEL:eq:cond1) is replaced by

\hb@xt@.01(1.8)

then it is said that weak CALR occurs. We refer to [2], [4] and [11] for more discussions on the anomalous localized resonance and its connection to invisibility cloaking.

The anomalous localized resonance phenomenon was first observed in [13] and connected to invisibility cloaking in [12]. Recently, a mathematical theory was developed in [2] by Ammari et al to rigorously explain the CALR observed in [13] and [12]. In their study, and is given by (LABEL:eq:d1) with

\hb@xt@.01(1.9)

Moreover, they let and , with . For the above plasmonic configuration, the solution to (LABEL:eq:m1) in [2] was shown to have a spectral representation associated with a Neumann-Poincaré-type operator. Using the spectral representation, it is shown that there exists a critical radius

\hb@xt@.01(1.10)

such that when a generic source term lies within , then CALR occurs; and when lies outside , then CALR does not occur. Here and also in what follows, when we say that CALR does not occur, it actually means that weak CALR does not occur. Later on, the CALR was considered from a variational perspective in [11] by taking and in (LABEL:eq:d1) with

\hb@xt@.01(1.11)

By using the primal and dual variational principles, it is shown in [11] that for a large class of sources, if , then resonance always occurs; whereas if and , then there exists the critical radius in (LABEL:eq:critical) for the occurrence and nonoccurrence of resonance.

The aim of this work is to extend the related results in [2, 11] to the three dimensional setting. Indeed, the three-dimensional CALR was considered in the literature and the situation becomes much more complicated. In [3], it is shown that if one takes a similar structure as that in (LABEL:eq:ammari) but with , then CALR does not occur. The same conclusion was draw in [10] without the quasi-static approximation. In [1], the anomalous localized resonance is shown to take place by using a folded geometry where the plasmonic material is spatially variable. In this paper, we show that CALR does not occur for the configuration (LABEL:eq:kohn) in . Then, we show that by properly choosing the plasmonic parameters, CALR can still happen, at least approximately. We follow both the spectral and variational arguments developed, respectively, in [2] and [11].

The rest of the paper is organized as follows. In Sections 2, by using the variational argument, we show the nonoccurrence and occurrence of ALR by taking the loss parameter to be given over the whole space . In Section 3, by following the relevant study in [2], we consider the occurrence and nonoccurrence of CALR by taking the loss parameter to be given only in the plasmonic layer.

2 Variational perspective on ALR in three dimensions

Henceforth, we let and be a constant. Throughout the present section, we let the source term be a real-valued distributional functional such that it is supported at a distance from the origin, and has a zero mean :

\hb@xt@.01(2.1)

where denotes the two-dimensional Hausdorff measure restricted to the set . Moreover, in this section, we let

be given by (LABEL:eq:d1) with and , \hb@xt@.01(2.2)

and without loss of generality, it is assumed that and . Indeed, the subsequent results derived in this section can be easily extended to the general case by a direct scaling argument. For the solution to (LABEL:eq:m1) with given in (LABEL:eq:d2), we set

\hb@xt@.01(2.3)

It is straightforward to verify that

\hb@xt@.01(2.4)
\hb@xt@.01(2.5)

Accordingly, the energy can be represented with and as

\hb@xt@.01(2.6)

The following variational principles were proved in [11] when , and can be extended to the three-dimensional case for the present study by straightforward modifications. Define

\hb@xt@.01(2.7)

Consider the optimization problem

\hb@xt@.01(2.8)

where , are assumed to be real-valued. (LABEL:eq:pv) is referred to as the primal variational problem, and the minimizing pair is attainable at such that is a solution to (LABEL:eq:m1). Similarly, we define

\hb@xt@.01(2.9)

and consider the following optimization problem

\hb@xt@.01(2.10)

where are assumed to be real-valued. (LABEL:eq:dv) is referred to as the dual variational problem, and the maximizing pair is attainable at such that is a solution to (LABEL:eq:m1).

We shall make use the variational principles introduced above to prove the resonance and non-resonance results. In doing so, the spherical harmonic functions for , and will be needed and they form an orthonormal basis to ; see [9]. In the rest of the current section, for ease of notations, we write instead of to signify the spherical harmonic functions of order . Set

Hence, the source in (LABEL:eq:source) can be written as

\hb@xt@.01(2.11)

where and also in what follows, for and , denotes the spherical coordinates. Moreover, in the subsequent arguments, we let and denote two generic positive constants that may change from one inequality to another, but should be clear from the context.

The following proposition will be needed and can be proved by direct verifications.

Proposition 2.1

Consider the PDE for : ,

\hb@xt@.01(2.12)

where

\hb@xt@.01(2.13)

with any . Then there exists a non-trivial solution which achieves its maximum value at a point with , given by

Moreover, one has that

\hb@xt@.01(2.14)

where denotes the jump of the normal flux of the function across .

2.1 Non-resonance result

In this subsection, we consider the non-resonance for the standard plasmonic configuration in (LABEL:eq:d2) with , and , . This is exactly the one considered in [11] for the two-dimensional case. We have

Theorem 2.1

Let be given in (LABEL:eq:d2) with , and , . Let be given in (LABEL:eq:source2). Then ALR does not occur.

Proof. We make use of the primal variational principle (LABEL:eq:pv) to prove the theorem. To that end, we first construct the test functions and that satisfy the PDE constraint in (LABEL:eq:pv) for . For , we set

\hb@xt@.01(2.15)

It is straightforward to verify that is continuous over , and satisfies

However, across , has a jump in its normal flux as follows,

\hb@xt@.01(2.16)

where denotes the exterior unit normal vector to . Therefore, if we let

\hb@xt@.01(2.17)

then one readily verifies that

\hb@xt@.01(2.18)

Hence, by setting

\hb@xt@.01(2.19)

with and , respectively, given in (LABEL:eq:decomp2) and (LABEL:eq:decomp1), then by (LABEL:eq:decomp4), one sees that there holds

\hb@xt@.01(2.20)

Finally, we let , then by virtue of (LABEL:eq:decomp5), clearly satisfies the PDE constraint in (LABEL:eq:pv). Hence, by the primal variational principle, together with straightforward calculations (though a bit tedious), one has that

That is, the resonance does not occur.

The proof is complete.     

2.2 ALR with no core

Theorem LABEL:thm:nr1 indicates that the standard plasmonic structure does not induce the ALR. In order for ALR to take place, one has to devise different plasmonic structures. Next, we first consider a construction without a core, which ensures that ALR can always occur.

Theorem 2.2

Consider the configuration described in (LABEL:eq:d2) with . Let be given in (LABEL:eq:source2) and assume that for some . Set . Then ALR occurs.

Proof. Since , we first assume that . Then we choose

\hb@xt@.01(2.21)

where satisfies and will be further chosen below. Clearly, by Proposition LABEL:prop:1, the pair satisfies the constraint . By the dual variational principle, we have

where the two positive constants and depend only on and . Choosing with , we obtain for .

Next, if , then by choosing

and using a similar argument as the previous case, one can show that resonance occurs.

The proof is complete.     

2.3 Approximate ALR with an arbitrary shape core

In this subsection, we assume that is a simply connected domain with a -smooth boundary. Let us consider a configuration given by (LABEL:eq:d2). We assume that , , is a symmetric-positive-definite matrix-valued function satisfying

\hb@xt@.01(2.22)

where and is the identity matrix. Fix and let . Let be the smallest integer satisfying

\hb@xt@.01(2.23)

and let , be given by

\hb@xt@.01(2.24)

Consider the source given in (LABEL:eq:source2). We shall prove that

Theorem 2.3

Consider the configuration described in (LABEL:eq:d2) with given in (LABEL:eq:assump1) and in (LABEL:eq:assump2). Let the source be given by (LABEL:eq:source2). Let and assume that the source satisfies

\hb@xt@.01(2.25)

Then for any , there exists a sufficiently small such that one has . That is, approximate ALR occurs.

Remark 2.4

Noting that and , we see that the condition (LABEL:eq:ca8) indicates that as long as the Fourier coefficient of the source in (LABEL:eq:source2) does not decay very quickly as , then approximate ALR occurs.

Proof. [Proof of Theorem LABEL:thm:r2] We make use of the dual variational principle to construct a sequence satisfying

First, we set

\hb@xt@.01(2.26)

where is to be chosen below. Let be the solution to By the standard elliptic estimate, one has

\hb@xt@.01(2.27)

where and are two positive constants depending on and in (LABEL:eq:assump1). Next, by straightforward calculations, we have

\hb@xt@.01(2.28)

By (LABEL:eq:assump2), we see that and , and hence the last two terms in the last inequality in (LABEL:eq:ca3) are of comparable order. Therefore, one further has from (LABEL:eq:ca3) that

\hb@xt@.01(2.29)

Using a completely similar argument by taking

\hb@xt@.01(2.30)

one can show that

\hb@xt@.01(2.31)

We choose to be

\hb@xt@.01(2.32)

where can be replaced by . Then one has

\hb@xt@.01(2.33)

Clearly, the estimate (LABEL:eq:ca7) also holds with replaced by . The proof can be immediately concluded by noting (LABEL:eq:ca8) and (LABEL:eq:ca7).     

2.4 Sensitivity and critical radius

By Theorem LABEL:thm:r2, we see that for any given , one can determine a sufficiently small and a sufficiently large according to (LABEL:eq:assump2) and (LABEL:eq:ca7), such that the configuration with given in (LABEL:eq:assump3) is “almost” resonant in the sense that . By (LABEL:eq:ca8) and (LABEL:eq:assump2), one has that as , and . Clearly, both and depend on , and hence of the plasmonic configuration depends on as well. It is natural to ask what would happen if one fixes the integer in (LABEL:eq:assump3). That is, in (LABEL:eq:assump3) is replaced by an integer , which can be as large as possible, but fixed. Next we show that resonance does not occur in such a case, and this indicates that the resonance is very sensitive to the plasmonic parameter.

Theorem 2.5

Let be any fixed positive integer and let . Let be given in (LABEL:eq:d2) with and . Suppose that the source is given in (LABEL:eq:source2). Then ALR does not occur.

Before giving the proof of Theorem LABEL:thm:sensitivity, we present another theorem whose proof would be more general than the one needed for Theorem LABEL:thm:sensitivity.

Theorem 2.6

Let be given in (LABEL:eq:d2) with and . Suppose that is given in (LABEL:eq:source2) and . Let be chosen according to (LABEL:eq:assump2). Then if with given in (LABEL:eq:critical), then ALR does not occur.

By Theorem LABEL:thm:r2, we see that ALR occurs if is chosen according to (LABEL:eq:assump3), namely, it is variable depending on the asymptotic parameter , and the source is located within the critical radius and satisfies the generic condition (LABEL:eq:ca8). However, by Theorem LABEL:thm:sensitivity, it is pointed out that the resonance phenomenon is very sensitive with respect to the plasmonic parameter , and if it is independent of the asymptotic parameter , then resonance does not occur. Theorem LABEL:thm:crc further shows that for the case with the variable plasmon parameter in Theorem LABEL:thm:r2, the resonance phenomenon is localized.

Proof. [Proof of Theorem LABEL:thm:crc] We make use of the primal variational principle to prove the theorem. To that end, we first construct test functions and that satisfy the PDE constraint in (LABEL:eq:pv).

Let be such that . Let be of the following form

\hb@xt@.01(2.34)

where , , satisfies

\hb@xt@.01(2.35)

and satisfies

\hb@xt@.01(2.36)

Define to be

\hb@xt@.01(2.37)

Using in (LABEL:eq:crc4), we then set

\hb@xt@.01(2.38)

where

\hb@xt@.01(2.39)

By straightforward calculations, though a bit tedious, one can verify that defined in (LABEL:eq:crc6) satisfies (LABEL:eq:crc2). Next, we define

and set

\hb@xt@.01(2.40)

It is directly verified that satisfies (LABEL:eq:crc3). Finally, we set

\hb@xt@.01(2.41)

Clearly, and satisfy the PDE constraint in (LABEL:eq:pv), and hence by the primal variational principle,

\hb@xt@.01(2.42)

where is defined in (LABEL:eq:v1).

We proceed to calculate the energy in (LABEL:eq:crc9) and show that it is bounded as , which readily implies that ALR does not occur. First, by (LABEL:eq:crc6), one can verify that

\hb@xt@.01(2.43)

Hence we have the following estimate

\hb@xt@.01(2.44)

Since , we have , which together with (LABEL:eq:crc10) implies that

\hb@xt@.01(2.45)

By (LABEL:eq:crc7), one can also calculate that

\hb@xt@.01(2.46)

Next we estimate the energy due to , and by (LABEL:eq:crc8) one has