arxiv:submit/ [math.ap] 18 Nov 2017

Transcription

1 ELECTROMAGNETIC FIELD ENHANCEMENT IN A SUBWAVELENGTH RECTANGULAR OPEN CAVITY YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN arxiv:submit/7549 [math.ap] 8 Nov 7 Abstract. Consider the transverse magnetic polarization of the electromagnetic scattering of a plane wave by a perfectly conducting plane surface, which contains a two-dimensional subwavelength rectangular cavity. The enhancement is investigated fully for the electric and magnetic fields arising in such an interaction. The cavity wall is assumed to be a perfect electric conductor, while the cavity bottom is allowed to be either a perfect electric conductor or a perfect magnetic conductor. We show that the significant field enhancement may be achieved in both nonresonant and resonant regimes. The proofs are based on variational approaches, layer potential techniques, boundary integral equations, and asymptotic analysis. Numerical experiments are also presented to confirm the theoretical findings.. Introduction The electromagnetic scattering by open cavities has significant applications in many scientific areas. For instance, the radar cross section RCS is a quantity which measures the detectability of a target by a radar system. Clearly, it is of high importance in military and civil use for a deliberate control in the form of enhancement or reduction of the RCS of a target. Since the cavity RCS can dominate the total RCS, it is critical to have a thorough understanding of the electromagnetic scattering by a target, particularly a cavity, for successful implementation of any desired control of its RCS. The underlying scattering problems have received much attention by many researchers in both the engineering and the mathematical communities. There has been a rapid development of the mathematical theories and computational methods in this area [3, 8, 9]. We refer to [6] for a survey of open cavity scattering problems. In optics, rough optical device surfaces containing subwavelength cavities have created a lot of interest in electromagnetic field enhancement and extraordinary optical transmission effects in recent years, due to their important potential applications, such as in biological and chemical sensing, spectroscopy, Terahertz semiconductor devices [6,,, 7, 8, 4, 3 33, 35]. The amplification of the confined fields in the cavity can be strikingly high. These amazing features of light localization and enhancement can be quite useful in imaging, microscopy, spectroscopy, and communication [7, 3]. However, there are still controversies on the mechanisms which contribute to anomalous field enhancement [8]. The complication arises from the multiscale nature of the structures and various enhancement behaviors that it induces. For instance, the enhancement can be attributed to surface plasmonic resonance [7, 8], non-plasmonic resonances [33, 35], or even without the resonant effect [8]. It is very helpful to provide a rigorous mathematical analysis for these phenomena, which turns out to be extremely challenging. The particular light patterns observed in the far-field or in the near-field may be the results of complex multi-interactions of surface waves, cavity resonances, resonant tunneling of plasmon waves, skin depth effects, etc. It requires the resolution of the full Mathematics Subject Classification. 45A5, 35C, 35Q6, 35C5. Key words and phrases. Cavity scattering problem, electromagnetic field enhancement, scattering resonances, Helmholtz equation, variational formulation, boundary integral equation, asymptotic analysis. The research of YG was supported in part by NSFC grant 5765, 677, National Research program of China Grant 3CB834 and FRFCU47FZ5. The research of PL was supported in part by the NSF grant DMS-538.

2 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Table. Electromagnetic field enhancement. : width of the cavity, d: depth of the cavity, λ: wavelength of the incident field. PEC-PMC cavity PEC-PEC cavity enhancement nonresonant resonant nonresonant resonant electric field Oλ/d O/ Oλ /d O/ magnetic field no O/ no O/ three-dimensional Maxwell equations in non-smooth geometries even for a qualitative description of the diffractive properties of electromagnetic waves in subwavelength structures. For most resonance phenomena, the localization or enhancement of light is very sensitive to geometrical parameters and the frequency. This motivates us to begin with simple geometries, where the analytical approaches are made possible to obtain an accurate description of the electromagnetic field. In [7, ], the field enhancement was considered for a single rectangle cavity and double rectangle cavities. Using asymptotic expansions of Green s function, they showed that the limiting Green function is a perfectly conducting plane with a dipole in place of the cavity when the width of the cavity shrinks to zero. In [3,,], an asymptotic expansion method was proposed to study the solution of the Helmholtz equation in a domain including a single thin slot. A quantitative analysis was made in [8, 9] for the field enhancement of the Helmholtz equation for a single narrow open slit. We refer [4,5] for the study of closely related Helmholtz resonators and resonances in bubbly media. In this paper, we consider the electromagnetic scattering of a plane wave by a perfectly electrically conducting plane surface, which contains a subwavelength rectangular cavity. More specifically, the cavity is assumed to be invariant in the x 3 -axis and only the transverse magnetic polarization TM is studied. Hence the three-dimensional Maxwell equations can be reduced into the two-dimensional Helmholtz equation. The cavity wall is assumed to be a perfect electric conductor PEC, and the cavity bottom is assumed to be either a perfect magnetic conductor PMC or a perfect electric conductor PEC. A cavity is called the PEC-PMC type when it has a PEC wall and a PMC bottom, while a cavity is called the PEC-PEC type when it has a PEC wall and a PEC bottom. We investigate the electromagnetic field enhancement for both types of cavities. Using a combination of variational approaches, layer potential techniques, boundary integral equations, and asymptotic analysis, we demonstrate the field enhancement mechanisms for the underlying scattering problems different boundary conditions. Denote by λ,, d the wavelength of the incident field, the width of the cavity, and the depth of the cavity, respectively. For the PEC-PMC cavity, we prove in the nonresonant regime that the electric field enhancement has an order Oλ/d and the magnetic field has no enhancement in the cavity provided that the length scales satisfy d λ. Moreover, if λ, we show that the Fabry Perot type resonance occurs. In the resonant regime, we derive an asymptotic expansion for the resonant frequencies and deduce an order O/ enhancement for both the electric and magnetic fields. For the PEC-PEC cavity, we carry the same analysis in both the nonresonant and resonant regimes. The electric field is shown to have a weak enhancement of order Oλ /d in the nonresonant regime with length scales satisfying d λ. In the resonant regime, a similar quantity analysis is done for the resonant frequencies and the same order O/ is achieved for the electric and magnetic field enhancement. These field enhancement results are summarized in Table. Numerical experiments are also done and used to confirm our theoretical findings. The rest of the paper is organized as follows. The PEC-PMC cavity is examined in sections 3 and the PEC-PEC cavity is studied in sections 4 5. Sections and 4 focus on the field enhancement for the non-resonant case. In each section, the model problems are introduced, the approximate problems are proposed, the error estimates are derived between the solutions of the exact and approximate problems, and the electric and magnetic field enhancement are discussed. Sections 3 and 5 address the field enhancement for the resonant case. In each section, the boundary integral

3 ELECTROMAGNETIC FIELD ENHANCEMENT 3 Γ D Γ Figure. The problem geometry of a rectangular open cavity. equations are given, the asymptotic analysis is carried for the Fabry Perot resonances, and the electric and magnetic field enhancement are presented under the resonant frequencies. The paper is concluded with some general remarks and directions for future work in section 6.. PEC-PMC cavity without resonance In this section, we discuss the electromagnetic field enhancement for the PEC-PMC cavity problem when the wavenumber is sufficiently small, i.e., no resonance occurs... Problem formulation. As is shown in Figure, the problem geometry is assumed to be invariant in the x 3 -axis. Let x = x, x R. Consider a two-dimensional rectangular open cavity D = {x R : < x <, d < x < }, where > and d > are constants. Denote by Γ = {x R : < x <, x = } and Γ = {x R : < x <, x = d} the open aperture and the bottom side of the cavity, respectively. Let R = {x R : x > } and Γ = {x R : x = }. Denote Ω = R D. We assume that the cavity and the upper half-space are filled with the same homogeneous medium. The electromagnetic wave propagation is governed by the time-harmonic Maxwell equations time dependence e iωt : E iωµh =, H iωεe = in Ω R,. where ω > is the angular frequency, ε > is the electric permittivity, µ > is the magnetic permeability, E and H are the electric field and the magnetic field, respectively. Furthermore, we assume that the boundary Ω \ Γ is a perfect electrical conductor PEC while the boundary Γ is a perfect magnetic conductor PMC. Hence E, H satisfy the following boundary conditions: ν E = on Ω \ Γ R, ν H = on Γ R,. where ν = ν, is the unit outward normal vector on Ω R, i.e., ν is the unit outward normal vector on Ω. We consider a polarized time-harmonic electromagnetic wave E inc, H inc that impinges on the cavity from above. In the transverse magnetic polarization TM, the magnetic field is parallel to the x 3 -axis, which implies that the incident magnetic wave H inc =,, u inc. Here u inc x = e iκx d is a plane wave propagating in the direction d = sin θ, cos θ, where π < θ < π is the incident angle and κ = ω εµ is the wavenumber. The corresponding incident electric field E inc is determined by Ampere s law E inc = iωε H inc = µ/εcos θ, sin θ, e iκx d. Denote by λ = π/κ the wavelength of the incident wave. In this section, we assume that the length scale of the underlying geometry satisfies d λ. Let H =,, u be the total magnetic field. It can be verified from.. that u satisfies u κ u = in Ω, ν u = on Ω \ Γ,.3 u = on Γ.

4 4 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Denote by u r the reflected field due to the interaction between the incident field u inc and the PEC plane Γ. It can be shown that u ref x = e iκx d, where d = sin θ, cos θ. In R, the total field u consists of the incident wave u inc, the reflected wave u ref, and the scattered field u sc. The scattered fields u sc is required to satisfy the Sommerfeld radiation condition: r r u sc lim r iκu sc =, r = x..4.. An approximate model. We introduce an approximate model for the problem.3.4 in order to estimate the field enhancement inside the cavity. Let φ x =, φ n x = cos nπx, n be an orthonormal basis on the interval,. It follows from the Neumann boundary condition in.3 that u can be expanded as the series of waveguide modes: u x = α n e iβnx α n e iβnx d φ n x, x D,.5 where the coefficients β n are defined as { κ, n =, β n = i nπ/ κ, n. If is small enough in fact it only needs < λ/ here, it follows that β n is are pure imaginary numbers for all n. It is easy to note that for each n, if nπ/ κ, the series consists of two propagating wave modes traveling upward and downward respectively; if nπ/ > κ the series consists of two evanescent wave modes decaying exponentially away from the bottom side and the open aperture of the cavity, respectively. Using the series.5, we may reformulate.3.4 into the following coupled problem: u κ u = in R, u = α n e iβnx αn e iβnx d φ n x in D, ν u = on Ω \ Γ,.7 u = on Γ, u x, = u x,, x u x, = x u x, on Γ, r r u s iκu s = in R, lim r where u sc = u u inc u ref in R. The continuity conditions are imposed along the open aperture Γ, where and indicate the limits taken from above and below Γ, respectively. For simplicity of notation, we define u x = u x, and u x = u x, d. Let u ±,n = u ±, φ n Γ ± be the Fourier coefficients for u and u, respectively, where the inner product, Γ ± is defined by u ±, φ n Γ ± := u ± x φ n x dx. For any s R, denote by H s Γ the trace functional space on Γ with the norm given by u H s Γ = nπ s u,n..8.6

5 ELECTROMAGNETIC FIELD ENHANCEMENT 5 Matching the series.5 on Γ, Γ and using the continuity conditions in.7, we get { α n αn e iβnd = u,n, α n e iβnd αn = u,n. Solving the above equations by Cramer s rule yields that α n = eiβnd u,n u,n e iβnd Therefore, in the cavity D, the total field e iβ nd u,n u,n u x = e iβnd, αn = eiβnd u,n u,n e iβnd..9 e iβnx eiβnd u,n u,n e iβnd e iβnx d φ n x.. Since u satisfies the homogeneous Dirichlet boundary condition on Γ, we have e iβ nd u,n u,n x u x, = iβ n e iβnd e iβnd eiβnd u,n u,n e iβnd φ n x on Γ,. e iβ nd u,n u,n u x, d = e iβnd e iβnd eiβnd u,n u,n e iβnd φ n x = on Γ.. It follows from. that Substituting.3 into. yields x u x, = u,n =..3 iβ n e iβnd e iβnd u,nφ n x on Γ..4 Given a function u H R and let u = ux,, we define a Dirichlet-to-Neumann DtN operator: B PMC e iβnd [u] = iβ n e iβnd u n φ n x on Γ,.5 where the Fourier coefficients u n = u, φ n Γ. Lemma.. The DtN operator B PMC is bounded from H / Γ H / Γ, i.e., B PMC [u] H / Γ Cλ, d u H / Γ u H / Γ, { } κ where the constant Cλ, d = max 4, sin κd. Proof. For any u, w H / Γ, they have the Fourier series expansions u = u n φ n x, w = w n φ n x, where u n = u, φ n Γ and w n = w, φ n Γ. It follows from.5 that B PMC e iβnd [u], w Γ = iβ n e iβnd u n w n. Note that the coefficients β n are given by.6. For n =, a straightforward calculation yields that iκeiκd e iκd u w κ cos κd = sin κd u w κ sin κd u w,.6

6 6 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN where we use the fact that d λ, which implies κd = πd λ n, using the fact that λ, we get iβ e iβnd n e iβnd u n w n nπ 4nπ u n w n 4 where we use the estimates nπ iβ n = κ nπ, tends to zero and cosκd /. For / nπ u n e iβnd e iβnd = e nπ/ κ d e nπ/ κ d 4. / w n,.7 Combining.6.7 and using.8, we obtain from the Cauchy Schwarz inequality that B PMC [u], w Γ Cλ, d u H / Γ w H / Γ, which completes the proof. Using.4.5 and the continuity of x u on Γ, we obtain the transparent boundary condition TBC: lim r x u = B PMC [u ] on Γ. Hence the problem.3.4 can be reduced to the following boundary value problem: u κ u = in R, ν u = on Γ \ Γ, x u = B PMC [u ] on Γ, r r u sc iκu sc = in R..8 By calculating the Fourier coefficients u,n, which give the coefficients u,n by.3, we may obtain the solution in the cavity D from the formula.. To find an approximate model to.8, we examine the series.5 more closely. Note that if λ, the wave modes e iβnx φ n x and e iβnxd φ n x decay exponentially in the cavity, with a decaying rate of Oe n/ for all n. Only the leading wave modes e iκx φ x and e iκxd φ x propagate in the cavity. The observation motivates us to approximate the DtN map.5 by dropping the high order modes and define an approximate DtN map by B PMC [v] = iκ eiκd e iκd v φ x..9 Now we arrive at an approximate model problem: v κ v = in R, ν v = on Γ \ Γ,. where v = on Γ single mode: and v sc x v = B PMC [v ] on Γ, r r v sc iκv sc = in R, v x = lim r = v u inc u ref in R. Accordingly, we may approximate u by one α e iκx α eiκx d φ x, x D,. where the coefficients α and α are given by α = eiκd v, v, e iκd, α = eiκd v, v, e iκd..

7 ELECTROMAGNETIC FIELD ENHANCEMENT 7 It follows from v = on Γ that α = α eiκd, which yields v, = Enhancement of the approximated field. This section introduces the estimates for the solution of the approximate model problem. Theorem.. Let v be the solution of the approximate model problem. and be given by. in the cavity, then there exist positive constants C, C independent of λ,, d such that C /d v L D C /d. Proof. In R, since the incident field u inc and the reflected field u ref satisfy the Helmholtz equation in., the scattered field v sc also satisfies v sc Noting that x u inc x u ref = on Γ, we have κ v sc = in R. x v s = on Γ \ Γ. Let G be the half-space Green function of the Helmholtz equation with Neumann boundary condition, i.e., it satisfies { Gx, y κ Gx, y = δx, y, x, y R, Gx,y ν y =, on Γ. It is easy to note that Gx, y = i 4 H κ x y H κ x y,.4 where y = y, y, H is the Hankel function of the first kind with order, and x is the reflection of the point x with respect to Γ, i.e., x = x, x. It follows from Green s identity that which gives v sc x = Γ v x = u inc x u r x Γ Gx, y vsc y ν y ds y, x R, Gx, y vsc y ν y ds y, x R. By the fact that x u inc x u ref = on Γ, especially on Γ, we have v sc ν y = v uinc u ref = v on Γ. ν y ν y ν y It follows from the continuity of single layer potentials [5, 6] and the above equality that v x = u inc x u ref x i H κ x y v y ds y, x Γ..5 ν y Using..3 yields that x v = iκ Γ α α eiκd φ x on Γ. Substituting the above equality into.5 and using the fact that φ x =, we get v x, = u inc x, u ref x, κ α α eiκd h x, x,,

8 8 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN where h x = H κ x y dy..6 Therefore, the Fourier coefficients v, may be expressed as v, = uinc, φ Γ u ref, φ Γ κ α α eiκd h, φ Γ. It follows from the fact u inc x, = u ref x, and. that α α eiκd = u inc, φ Γ κ h, φ Γ α α eiκd. By α = α eiκd and the above equation, we obtain α = u inc, φ Γ c c e iκd, α = eiκd u inc, φ Γ,.7 c c eiκd where c = κ h, φ Γ. Note that x, y Γ for small x y, asymptotically, it holds that cf. [34] H κ x y = i π ln x y i π ln κ γ O x y ln x y, where γ is the Euler constant defined by γ = lim c = κ h, φ Γ = κ It follows from the assumption d λ that τ τ m= m ln τ. A direct calculation yields H κ x y dy dx = iκ ln O..8 π c c e iκd e iκd = 4κd,.9 c c e iκd eiκd κd..3 Using the definition of the incident wave gives u inc, φ Γ = e iκ sin θx dx = O 3..3 Combining.7.3, we obtain α C κd, where the positive constant C is independent of λ,, d. It follows from. that v L D Substituting.7 into. yields v x = d u inc, φ Γ c c e iκd e iκx α C κd, κ α α φ x dx dx C /d. e iκd u inc, φ Γ c c e iκd eiκx d φ x..3

9 ELECTROMAGNETIC FIELD ENHANCEMENT 9 A simple calculation yields v = iκe iκd u inc, φ Γ x c c e iκd φ x e iκxd e iκx d = iκe iκd u inc, φ Γ c c e iκd cos κx dφ x..33 Noting that is small enough and κx d, we have from.9.3 and.33 that there exists a positive constant C such that v x C d, which yields that which completes the proof. v L D d v x dx dx C /d, Theorem.3. Let v be the solution of the approximate model problem. and be given by., then there exist two positive constant C 3 and C 4 independent of, d, λ such that Proof. It follows from.3 that d C 3 d v L D C 4 d. u inc, φ v x, x Γ e iκd dx dx = c c iκd Noting that d λ, we have and u inc, φ Γ = c c iκd d d d 4 sin κx ddx 4 4 sin κx ddx = d d 4 sin κx ddx. It follows from.9.3 that u inc, φ c κ d Γ e iκd c c iκd c κ d d e iκx d e iκx d dx κx d dx 4 3 κ d 3.34 sin κd κ = 4 3 κ d 3 Oκ 4 d 5 3 κ d for some positive constants c and c. Combining the above estimates, we obtain that there exist constants C 3 and C 4 independent on λ,, d such that C 3d v L D C 4d, which completes the proof.

10 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN B R B R Γ Γ R Γ R D Γ Figure. The problem geometry with a transparent boundary..4. Accuracy of the approximate model. To show the accuracy of the approximate model, we introduce a TBC in R and reformulate the problems.8 and. in a bounded domain. Let B R = {x R : x x c < R, x > } be the upper half disc with radius R centered at x c = /,, where R is a sufficiently large positive real number. Denote by B R = {x R : x x c = R, x > } be the upper half circle. Let Γ R = { x R : x / < R, x = } be the line segment on Γ with length R. The problem geometry is shown in Figure. In the exterior domain R \ B R, by noting the PEC condition on Γ \ Γ R and the radiation condition.4, the scattered field can be expressed as the Fourier series: u sc r, η = H n κr H n κr usc,n cosnη, < η < π, where u sc,n = π π usc R, η cosnηdη, H n is the Hankel function of the first kind with order n, and the polar coordinate r = x / x. Taking the normal derivative of usc on B R yields u s R, η r = κh n κr H n κr For any v H / B R with the Fourier expansion v = v n cosnη, v n = π we define the DtN operator on B R : T vη = κh n π u sc,n cosnη..36 vη cosnηdη, κr v n cosnη, < η < π,.37 H n κr It follows from.36 and.37 that the normal derivative of total field on B R where g = r u inc u ref T u inc u ref. can be written as r u = r u sc r u inc u ref = T u g,.38 Lemma.4. The DtN operator T : H / B R H / B R is bounded. Proof. Let A n = H n cf. [], it yields that κr H n κr. By the relation of the Hankel function and the modified Bessel function A n = K niκr K n iκr,

11 ELECTROMAGNETIC FIELD ENHANCEMENT where K n is the modified Bessel function of the second kind. Using the formula zk nz = nk n z zk n z, we obtain A n = K niκr K n iκr = n iκr K n iκr K n iκr n κr c n, where c is a positive constant. For any v, w H / B R, we have from the definition of.37 that Thus we have T v, w B R = κr T v H / B R = which completes the proof. Lemma.5. We have cκπr cκπr π A n v n cosnη wdη = n v n w n n v n = CR, λ v H / B R w H / B R. κπr A n v n w n n w n T v, w sup B R CR, λ v w H / B R w H / H / B R B, R Re T v, v B, Im T v, v R B >, v H / B R R. Proof. It follows from the definition.37 that T v, v B = κπr R A n v n..39 Note that H n = J n iy n, where J n and Y n are the Bessel functions of the first and second kind, respectively, and the modulus of each Hankel function is decreasing function [], then ReA n = J nkrj nkr Y n kry nkr J nkr Y n kr It follows from the Wronskian formula [], ImA n = J nκry nκr Y n κrj nκr JnκR Yn κr = W J nκr, Y n κr JnκR Yn = κr κπr JnκR Yn κr = Jn kr Yn kr JnkR Yn..4 kr >,.4 where W, denotes the Wronskian determinant of J n and Y n. The proof is completed by combining By the TBC operator.38, the model problem.8 can be reformulated as follows: u κ u = in B R, ν u = on Γ R \ Γ, x u = B PMC [u ] on Γ, r u = T [u ] g on B R..4

12 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN The variational formulation of.4 is to find u H B R such that au, w = g, w B R where the sesquilinear form au, w = u w κ u w dx B R B R w H B R,.43 T [u ] wdγ Γ Theorem.6. The variational problem.43 has a unique solution u H B R. Proof. By the definition.5, we have for sufficiently small that B PMC e iβnd [u ], u Γ = iβ n e iβnd u,n = κ cos κd sin κd u, n= It follows from Lemma.5 and the above inequality that B PMC [u ] wdγ..44 nπ/ κ e nπ/ κ e nπ/ κ u,n >. Re {au, u } = u L B R κ u L B R Re T [u ], u B R B PMC [u ], u Γ u L B R κ u L B R, which shows that the sesquilinear form a, satisfies a Gårding type inequality. Next is show the uniqueness. It suffices to show that u = if g =. A simple calculation yields Imau, u = Im T [u ], u B = κπr ImA n u,n =. R It follows from Lemma.5 that u,n = for all n. Therefore u = on B R. By the definition.37, we have T u =, which implies r u = T u = on B R by.38. Hence we have u = in R \ B R. It follows from the unique continuation [] that u = in B R. The proof is completed by applying the Fredholm alternative theorem. Correspondingly, the approximate model problem. can be written as follows: v κ v = in B R, ν v = on Γ R \ Γ, x v = B PMC [v ] on Γ, r v = T [v ] g on B R..45 The variational formulation of.45 is to find v H B R such that a v, w = g, w B R where the sesquilinear form a v, w = v κ v w dx B R B R w H B R,.46 T [v ] wdγ Γ B PMC [v ] wdγ..47 The well-posedness of the variational problem.46 can be shown similarly to Theorem.6. The proof is omitted for brevity. In the cavity D, the wave fields u and v are given by. and.., respectively. Their Fourier coefficients can be calculated from the solutions obtained in.4 and.45. The following lemma gives the accuracy of approximate model.45 and...

13 ELECTROMAGNETIC FIELD ENHANCEMENT 3 Lemma.7. Let u be the expansion of. inside the cavity with the Fourier coefficients given by the model.4, and let v be the expansion of.. with the Fourier coefficients given by the approximate model.45, then the following estimates hold: u v L D Cλ, d, R, where Cλ, d, R is a positive constant independent of. u v L D Cλ, d, R ln, Proof. The proof consists of five steps in order to show the estimates between u and v. Step : we first give the estimate of u u H B R and v u H B R, where u is the solution of.4 when there is no cavity, i.e., u satisfies the variational problem where the sesquilinear form bu, w = bu, w = g, w B R B R w H B R,.48 u w κ u w dx T [u ] wdγ. B R It follows from.43 and.48 that au u, w = g, w B au, w = bu, w au, w = R Using.5 and.9, we have from Lemma. that bu, w au, w B PMC [u ] wdγ where u, = u, φ Γ and w = w, φ Γ. Using the estimates [8, Lemma A., Lemma A.] we obtain from.5 that bu, w au, w Cλ, d, R Γ Γ Γ B PMC [u ] wdγ..49 B PMC B PMC [u ] wdγ κ sin κd u, w C u H / Γ w H / Γ,.5 u, CR u H B R if u H B R, w CR ln w H B R if w H B R, u H / Γ CR u H 3 B R if u H 3 B R, w H / Γ CR w H B R if w H B R, ln u H B R w H B R u H 3 B R w H B R. Since u = u inc u ref = e iκx d e iκx d, we have from straightforward calculations that bu, w au, w Cλ, d, R ln w H B R..5 Denote by L be the induced operator for the sesquilinear form.44 such that a u, w = L u, w, where, is the inner product on L B R. It follows from [8, Lemma A.3] that there exists a positive constant such that L LH B R Cλ, d, R < <,.5 where Cλ, d, R is a constant independent of. Combing.49.5 yields u u H B R Cλ, d, R ln.

14 4 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Since v satisfies the variational problem.46, following the same arguments as above, we may show for sufficiently small that v u H B R Cλ, d, R ln..53 Step : we estimate v u H 3/ B R. Let ξ = v u. It is easy to verify that ξ satisfies ξ κ ξ = in B R, ν ξ = on Γ R \ Γ, x ξ = B PMC [v ] on Γ, r ξ = T ξ on B R. It follows from the standard regularity estimate for the Neumann problem that ξ H 3/ B R Cλ, d, R BPMC [v ] L Γ Cλ, d, R u, v w, u, φ L, u Cλ, d, R H B R ln v u H B R Cλ, d, R,.54 where the last inequality is obtained by using.53. Step 3: This step is to estimate u v H B R. From.43 and.46, we have au v, w = a v, w av, w..55 Comparing.44 and.47 yields a v, w av, w B PMC B PMC Γ [v ] wdγ B PMC B PMC [u ] wdγ B PMC B PMC [v u ] wdγ..56 Γ Using the estimates.5.5 gives B PMC B PMC [u ] wdγ Cλ, d, R w H B. R Γ For the second term in.56, we have from the estimate of v u H 3/ Ω R that Γ B PMC B PMC [v u ] wdγ Cλ, d, R v u H / Γ w H / Γ Γ Cλ, d, R v u H 3/ B R w H B R Cλ, d, R w H B R. For the variational problem.55, a combination of the above estimates and the boundedness of inverse of the linear operator L in.5 lead to the following estimate u v H B R Cλ, d, R..57

15 ELECTROMAGNETIC FIELD ENHANCEMENT 5 Step 4: we estimate u v L D. We have from the expansions of u in.7 and v in. that u v L D = κ α α e iκx α α eiκx d dx d nπ α d n e iβnx αn e d iβnx n= iβ n α n e iβnx αn e iβnx d dx κ d α α α α nπ/ α nπ/ n αn, where the Parseval identity is used in the last inequality. For the first term n = in the above inequality, following from.9 and. and noting that d λ, we get n= α α = u, v, e iκd κ d u, v,, α α = eiκd u, v, e iκd κ d u, v,..58 For the second term n, by.9 and iβ n d = d nπ/ κ, we have Consequently, we obtain α n = u,n e iβnd u,n, α n = eiβnd u,n u v L D Cλ, d, R u, v, It follows from.54 and.57 that e iβnd u,n..59 n= nπ Cλ, d, R u, v, u H / Γ. / u,n u H / Γ u H / Γ v u H / Γ u v H / Γ Cλ, d, R u H 3 B R v u H 3/ B R u v H B R Cλ, d, R..6 On the other hand, we have from.57 that u, v, Cλ, d, R ln u v H B R Cλ, d, R ln..6 Combining the above estimates yields u v L D Cλ, d, R.

16 6 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Step 5: the last step is to estimate u v L D. Using the similar method as that in Step 4, we can obtain from the Parseval identity and that u v L D d α α α α α nπ/ n αn Cλ, d, R Cλ, d, R u, v, n= u, v, Cλ, d, R Cλ, d, R ln 3 n= nπ n= u,n nπ u, v, u H / Γ u,n Cλ, d, R ln, where the estimates.6.6 are used..5. Field enhancement without resonance. The electromagnetic field enhancement factor is defined as the ratio of the energy between the total field and the incident field in the cavity. Hence the electric and magnetic field enhancement factors are given by Q E = E L D E inc L D and Q H = H L D H inc. L D Theorem.8. There exists a constant depending on λ and d such that for all < < the following estimates hold: C λ/d Q E C λ/d, C Q H C, where C and C are some positive constants independent of λ,,, d. Remark.9. The electric field enhancement has an order O λ/d, which is enormous due to the length scale d λ, while the magnetic field has no significant enhancement in this situation. Proof. Recalling the scale assumption d λ, we have from Theorem. and Lemma.7 that u L D v L D u v L D C d Cd, λ, R 3C d, where we use the fact that Cλ, d, R d C for sufficiently small. On the other hand, A straightforward calculation yields u L D v L D u v L D C d Cλ, d, R C d. u inc L D = Combining the above estimates gives d C λ 4π d u L D u inc L D u inc x dx = κ d. 3C 4π λ d.

17 ELECTROMAGNETIC FIELD ENHANCEMENT 7 Noting that H =,, u = u and H = iωεe, we obtain C λ 4π d u L D u inc L D = E L D E inc L D = Q E 3C 4π For the magnetic field, it follows from Theorem.3 and Lemma.7 that λ d. and u L D v L D u v L D 3 C d Cλ, d, R ln C d u L D v L D u v L D C d Cλ, d, R ln C d, where we use Cλ, d, R ln /d min {C /, C /} for sufficiently small. Clearly, we have Therefore, we obtain which completes the proof. u inc L D = C u L D u inc L D d u inc x dx = d. = H L D H inc L D = Q H 3 C, 3. PEC-PMC cavity with resonances This section is to discuss the Fabry Perot type resonance and derive the asymptotic expansions for those resonance; analyze quantitatively the field enhancement at the resonance frequencies in a narrow cavity, i.e., λ. Recall the model problem.3.4: u κ u = in Ω, ν u = on Ω \ Γ, u = on Γ 3., lim r r r u s iκu s = in R. It is known that the scattering problem 3. has a unique solution for a complex wavenumber k with Imk cf., [3]. By analytic continuation, the solution has a meromorphic extension to the whole complex plane, except for a countable number of points, which are poles of the resolvent associated with the scattering problem 3.. These poles are called the resonances of the scattering problem, and the associated nontrivial solutions are called resonance modes. If the frequency of the incident wave is taken to be close to the real part of the resonance, an enhancement of the field is expected if the imaginary part of the resonance is small. 3.. Boundary integral equation. Define by g PMC x, y the Green function in the PEC-PMC cavity D, i.e., it satisfies g PMC x, y κ g PMC x, y = δx y, x, y D, g PMC x,y ν y = on D \ Γ, g PMC x, y = on Γ.

18 8 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Since the width of the cavity tends to zero, we may assume that κ is not the eigenvalue of the Laplacian with the above boundary conditions in the cavity. Thus the Green function of the Helmholtz operator in D exists and can be expressed as [4] g PMC x, y = m, c m,n φ m,n xφ m,n y, where c m,n = κ mπ/ n /π/d, αm,n φ m,n x = d cos mπx n /π sin x d, 3. d {, m =, n, α m,n = 4, m, n. Lemma 3.. The scattering problem 3. is equivalent to the following boundary integral equation: i H u y κ x y g x, y ds y u inc x u ref x =, x Γ. 3.3 ν y Γ Proof. Noting that ν u inc ν u ref = on Γ, especially on Γ, we have ν u sc = ν u ν u inc u ref = ν u on Γ. It follows from Green s identity in R that u x = u inc x u ref x u sc x = u inc x u ref x Gx, y usc y ds y Γ ν y = u inc x u ref x Gx, y u y ds y, x R ν, 3.4 y Γ where the half-space Green function G is given in.4. Using Green s identity inside the cavity yields u x = g Γ PMC x, y u y ds y, x D. 3.5 ν y By the continuity of the single layer potential, we have from 3.4 and 3.5 that u x = u inc x u ref x i H κ x y u y ds y, x Γ ν y and u x = Γ Γ g PMC x, y u y ds y, x Γ. ν y The proof is finished by imposing the continuity of the solution on the open aperture Γ. It is clear to note that ν u Γ = x u x,, u inc u ref = e iκ sin θx, x,.

19 ELECTROMAGNETIC FIELD ENHANCEMENT 9 Introduce new variables by rescaling with respect to : X = x /, Y = y /. Define three functions: ϕy = y u Y,, G e X, Y = i H k X Y, G i X, Y = g PMC X, ; Y, = and the boundary integral operators: T e ϕx = T i ϕx = m, c m,n α m,n d G e X, Y ϕy dy, X,, cosmπx cosmπy 3.6 G i X, Y ϕy dy, X,. 3.7 Hence the boundary integral equation 3.3 is equivalent to the following operator equation: where fx := u inc u ref X, = e iκ sin θx. T e T i ϕ = f/, Asymptotics of the integral operators. In this subsection, we study the asymptotic properties of the integral operators in 3.7. For clarity, we first introduce several notation below. Γ κ, = π ln κ γ ln, π Γ tan κd κ, = ln κ π, RX, Y = κ 4π X Y ln, RX, lnκ X Y Y = γ κ X Y oκ X Y, R 3X, Y = κ π π 4π m= X Y m3 4 ox Y X Y, lnπx Y X Y 4 lnπx Y where γ = γ ln πi, γ = i 4π i 8 4π ln γ, and γ is the Euler constant. For the fixed depth of the cavity d and the wave number κ, if small enough, we have the following asymptotic expansions for the kernels G e and G i. Lemma 3.. If λ, then the following estimates hold: where G e = Γ κ, π ln X Y R X, Y RX, Y, G i = Γ κ, ln πx Y π sin ln sin πx Y R3X, Y, R X, Y C ln, R X, Y C, R 3X, Y C 3. Here C j, j =,, 3 are positive constants independent of.

20 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Proof. It follows from the asymptotic of the Hankel function near zero cf. [] that G e X, Y = i H κ X Y = π γ ln κx Y 4π kx Y ln κx Y γ κ X Y oκ X Y = Γ κ, π ln X Y R X, Y R X, Y. The integral kernel G i X, Y can be expressed as G i X, Y = c m,n α m,n cosmπx cosmπy. 3.9 d m, Recalling the definitions of c m,n and α m,n in 3., we set α m,n C m, κ = c m,n α m,n = κ mπ n Using the fact cf. [9] that we may check for m = that Using the identity we have for m that m= tan πx = 4x π C κ = C m, κ = π d n x, x n, n Z, κ tanh πx = 4x π = n π d = d tanκd. κ n x, 4 κ mπ n π d mπ mπ d tanh κ = d mπ 3 κ d m 3 π 3 O d 5 κ 4 d 4!m 5 π 5. κ Using the following facts cf. [3, 5]: cos mπx m = ln ln sin πx, < x <, m= cos mπx m 3 = m 3 πx lnπx Ox, < x < m=.

21 ELECTROMAGNETIC FIELD ENHANCEMENT and substituting C κ, C m, κm into 3.9, we obtain G i X, Y = d C κ cosmπx cosmπy d mπ = tanκd κ m= 3 κ d m 3 π 3 cosmπx cosmπy κ 4 d O 5 4!m 5 π 5 m= ln ln πx Y π sin ln sin πx Y κ X Y X Y π π lnπx Y lnπx Y m3 4 4 m= OX Y X Y, which complete the proof. Let Γ = Γ κ, Γ κ,, k X, Y = RX, Y, k = RX, Y R3X, Y, kx, Y = ln X Y ln πx Y π sin ln sin πx Y. 3. Denote by K, K, K the integral operators corresponding to the Schwarz kernels kx, Y, k X, Y and k X, Y, respectively. Let I be a bounded open interval in R and define H s I := {u = U I, U H s R}, where s R. Then H s I is a Hilbert space with the norm { } u H s I = inf U H s R U Hs R and U I = u. Define H s I := { } u = U I U H s R and suppu Ī. It can be shown that the space H s I is the dual of H s I and the norm for H s I can be defined via the duality cf. []. We also define the operator P : H /, H /, by P φx = φ,, L,,, where the duality between H /, and H /, is defined by, L,,, is the function defined on the interval, and is equal to one. It is easy to note that, H /,. Lemma 3.3. The following conclusions hold: i The operator T e T i admits the decomposition: T e T i = ΓP K K K. Moreover, K and K are bounded from H /, to H /, with the operator norm satisfying K ln, K uniformly for wavenumber κ. ii The operator K is invertible from H /, to H /,. Moreover, the constant q := K,,, L,.

22 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Proof. The proof of i follows directly from the definitions of the operators T e, T i in 3.7 and the asymptotic expansions of their kernels given by Lemma 3.. The proof of ii follows directly from [, Theorem 4., Lemma 4.] Asymptotics of the resonances. The scattering resonance of 3. is defined as a complex wavenumber κ with negative imaginary part such that there is a nontrivial solution to 3. when the incident field is zero. This is the characteristic value of the operator ΓP K K K with respect to the variable κ. For simplicity, we write ΓP K K K := P L, where P = ΓP, L = K K K. By Lemma 3.3, it is easy to see that L is invertible for sufficiently small. Now assume that there exists ϕ such that Then P L ϕ =. L P I ϕ =, where I is the identity operator. It is straightforward to check that the eigenvalues of operator L P I are λκ, = Γκ, L,,, L,. Therefore, the characteristic values of the operator P L are the roots of the analytic functions λκ,, and the associated characteristic function is given by ϕ = Γκ, L,. Lemma 3.4. The resonance of the scattering problem 3. are the roots of the analytic functions λκ, =. Moreover, L, = K, O ln O, L,,, L, = q O ln O. 3. Proof. For given roots of λκ,, it is easy to check that they are the characteristic values of the operator P L with corresponding characteristic function defined above. It follows from the asymptotic expansions in Lemma 3.3 and the Neumann series that L = K K K = K j K K K j which gives 3.. = K O ln O, Next, we present the asymptotic expansion of resonances for the scattering problem 3.. Theorem 3.5. The scattering problem 3. has a set of resonances {k n }, which satisfy k n = nπ d nπ d π ln q π ln ln nπ d γ O ln, n =,,. 3. Proof. By Lemma 3.4, we consider the root of j= λκ, = Γ κ, Γ κ, L,,, L, =. Recall that Γ κ, = π ln κ γ π ln and Γ tan κd κ, = κ ln π. The above equation can be written as tan κd κ π ln ln κ γ π ln L,,, L, =.

23 ELECTROMAGNETIC FIELD ENHANCEMENT 3 Using Lemma 3.4, we obtain pκ, := λκ, = tan κd ρκ q κ π ln O ln O, 3.3 where ρκ = π ln ln κ γ. Note that the axis {z C : argz = ± π d } is the branch cut for pκ,, hence we choose a small number θ > and consider the domain {z C : π d θ < argz < π d θ or π d θ < argz < 3π d θ }. On the other hand, we are only interested in those resonances which are not in the high frequency regime. Therefore, we only need to find all the roots of pκ, in the domain { Ω θ,m = z C : z M, π d θ < argz < π d θ or π d θ < argz < 3π } d θ for some fixed number M >. Since pκ, blows up as κ or κ j /π/d, j Z. As a result, there exists δ > such that all the roots of p in Ω θ,m actually lies in the smaller domain Ω δ,θ,m := {z C : z δ or z j /π/d δ, j Z} Ω θ,m. tan κd It is clear that κ is analytic in Ω δ,θ,m and its roots are given by k n, = nπ/d, n =,,. Note that each root is simple. Denote by k n the roots of λκ,. From Rouche s theorem, we deduce that k n are also simple and are close to k n, if is sufficiently small. We now derive the leading order asymptotic terms for these roots. Define tan κd p κ, = ρκ κ π ln q. 3.4 Expanding p κ, at k n, yields tan κd p κ, = κ k n ρk n ρ k n κ k n Oκ k n κ κ=k n π ln d = nπ κ k n ln ln nπ π d γ d nπ κ k n Oκ k n π ln q. We can conclude that p κ, has simple roots in Ω δ,θ,m which are close to k n s. Moreover, these roots are analytic with respect to the variables and ln. Denoting the roots of p by k n and expanding them in term of and ln, we obtain k n =k n k n d π ln q π ln ln nπ d γ O ln. 3.5 Using Rouche s theorem, we claim that k n gives the leading order for the asymptotic expansion of the roots k n. More precisely, we have k n = k n O ln. 3.6 In the following, we give a proof of the claim. Note that tan κd pκ, p κ, = κ π ln ρκ O ln and tan κd p κ, = q q κ π ln ρκq. Hence there exists a constant C n such that pκ, p κ, < p κ,, κ { κ k n = C n ln }. which implies the claim 3.6 by Rouche s theorem. q

24 4 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN 3.4. The field enhancement with resonance. First we give an asymptotic expansion of pκ, at the resonant frequencies. Lemma 3.6. At the resonant frequencies κ = Rek n where k n is given in 3., we have pκ, = iq O ln. Proof. Assume that κ Rek n ln. It follows from the definition of p in 3.4 that pκ, = p κ, O ln = p k n κ k n Oκ k n O ln d = O k n k n q κ k n O ln k n πk n d = q κ Rek n iimk n O ln. k n πk n Recalling γ = γ ln πi, we have from the expansion of k n in 3.5 that Imk n = Imk n O ln = k n dπ Imγ = k n d. When κ = Rek n, we have pκ, = iq d O ln = iq O ln, which completes the proof. Lemma 3.7. The solution ϕ of equation 3.8 has the following asymptotic expansion in H /, : Moreover, ϕ = K, sin θ Oκ p ϕ,, L, = q sin θ Oκ p O ln. O ln. Proof. It follows from Lemma 3.3 that equation 3.8 is equivalent to L P I ϕ = L f/. Recall that the operator L P I has the eigenvalue λκ,. Thus ϕ,, L, = λκ, L f/,, L,, and ϕ = L f/ L ΓP ϕ = L f/ L Γ ϕ,,, = L f/ L, Γ λκ, L f/,, L,. Recall that λκ, = Γ L,,, L,. A simple calculation yields ϕ = L f/ L λκ, / L,,, L,, L f/, λκ,, L,. Using the definition of f and expanding it with respect to, we get f = eiκx sin θ = iκx sin θ Oκ,

25 ELECTROMAGNETIC FIELD ENHANCEMENT 5 which gives L f/ = sin θ Oκ K, O ln. Combining the above equation with 3., we obtain ϕ = sin θ Oκ K, O ln Thus λ/q O ln q q sin θ Oκ K λ, O ln sin θ Oκ = K λ, O ln. which completes the proof. ϕ = K sin θ Oκ, O ln, p ϕ,, L, = q sin θ Oκ O ln, p Combining Lemmas 3.6 and Lemma 3.7, we obtain at the resonant frequencies κ = Rek n that where f = K, ϕ = f O ln, ϕ,, L, = c O ln, 3.7 sin θ Oκ i q O ln H /,, c = q sin θ Oκ i q O ln Enhancement in the far field. We first investigate the scattered field in the domain R \ B R. Recall that u sc x = Gx, y usc y ds y = Gx, y u y ds y, x R Γ ν y Γ ν. y Here we use the fact ν u inc ν u ref = on Γ, especially on Γ. Using u y ν y = u y, y = ϕ, Γ y we have Since u sc x = Γ we get from Lemma 3.7 that Hence Gx, y, ϕ y dy = Gx, Y, ϕy dy. Gx, Y, = Gx,, O, x R \ B R, ϕy dy = ϕ,, L, = q sin θ Oκ p u sc x = Gx,, O ϕy dy = Gx,, q sin θ Oκ p O ln. O ln.

26 6 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN In the resonance case when κ = Rek n, the enhancement comes from the term p. It follows from Lemma 3.6 that Correspondingly, we have p = i q O ln. u sc x = igx,, sin θ Oκ O ln, which shows that the enhancement of the scattered magnetic field in the far field region has an order O/ compared to the non-resonance case when κ Rek n. It follows from the Ampere law H = iωεe and H =,, u that the enhancement of the scattered electric field also has an order O/ in the far field region when κ = Rek n Enhancement in the cavity. Now we study the field enhancement in the cavity. It follows from.5 that the total field u may be expanded as the sum of waveguide modes: mπx α m cos e iβmx u x = α e iκx α eiκx d m where β m = i mπ/ κ. m α m cos mπx Lemma 3.8. The coefficients in 3.8 have the following expansions: and α = q sin θ Oκe iκd O ln, iκ cosκdp α m C m, e iβmx d, 3.8 α = q sin θ Oκ O ln, 3.9 iκ cosκdp α m where the positive constant C independents of, κ and m. Proof. Recall that u = on Γ then equation 3.8 becomes u x = α e iκxd e iκx d m C m, α m e iβmxd e iβmx d cos mπx Taking the derivative of 3. with respect to x yields u x = α iκ e iκxd e iκx d x α miβ m e iβmxd e iβmx d cos m. 3. mπx, which gives u x, = α iκ cosκd x α miβ m cosβ m d cos m mπx. 3.

27 ELECTROMAGNETIC FIELD ENHANCEMENT 7 Multiplying 3. with φ x and integrating it over Γ, we get α iκ cosκd = u x, dx = x ϕ dx x = Γ Γ ϕxdx = ϕ,, L, = q sin θ Oκ p O ln, which gives the formulas for α and α in 3.9. For m, it follows from multiplying3. with cos mπx and integrating over Γ that α miβ m cosβ m d = u mπx x, cos dx x = Γ ϕx cosmπxdx. Note that iβ m = O m for m, and by 3.7 that ϕ H /,, cosmπx H /, m. Thus, we get which completes the proof. α m C m, In the following, we give the field enhancement inside the cavity. Theorem 3.9. Denote by D = { x D d x } the interior of the cavity. If κ = Rek n where k n is given in 3., then we have for d that i sin u x = sin θ O κx d Oln O ln, x κ cosκd D. Moreover, the electric and magnetic field enhancements are of an order O/. Proof. By the expansion 3.8 and Lemma 3.8, it is clear that in the region D u x = q sin θ Oκ e iκxd e iκx d O ln Oe / iκ cosκdp = q sin θ Oκ sin κx d κ cos κd p O ln Oe /. Noting when κ = Rek n, we have which yields that u x = p = i q O ln, i sin sin θ O κx d Oln O ln. κ cosκd Therefore, the magnetic field enhancement due to the resonance is of an order O/ in D. We conclude from Ampere s law that the electric field enhancement is also of an order O/.

28 8 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN Field enhancement on the open aperture. Now we present the enhancement of the scattered field on the open aperture of the cavity. Recall 3.4 for the scattered field on Γ : u sc x, = Gx,, y u y ds y. ν y Since x = X, y = Y, ϕ = ν u Γ, we have Γ u sc x, = G e X, Y ϕy dy. It follows from the asymptotic expansion of G e in Lemma 3. and Lemma 3.7 that u sc x, = Γ π ln X Y O ln ϕy dy sin θ Oκ = Γ ϕ,, L, p π h X O ln = π ln κ γ π ln q sin θ Oκ O ln p sin θ Oκ p π h X O ln = q π p ln π p q ln κ γ h X π q ln κ γ h X q sin θ O ln O ln. π p sin θ O p where h X := K, Y ln X Y dy. Note that at the resonant frequencies κ = Rek n, p = i q O ln. Thus we have u sc x, = 4i 4i ln ln κ γ h x / O ln 3. π π q It is clear to note that the leading order of the resonant mode is a constant on the open aperture with an order of Oln and the magnetic field enhancement is of an order O/. On the other hand, since x u sc sin θ Ok x, = p π h x / O ln, which implies that at the resonant frequencies κ = Rek n, on the top aperture Γ, the scatter field satisfies x u sc x, = 4i O ln h x / πq i sin θ O h πq x / O ln. By Ampere s law, we know the electric field enhancement is also of an order O/ on the top aperture Γ when κ = Rek n Numerical experiments. In this section, we present some numerical experiments to verify the theoretical studies of the field enhancement for the PEC-PMC cavity. As a representative example, we take the cavity width =.5, the cavity depth d =, and the angle of incidence θ = π/3 in the numerics. Figure 3 left and right show the plot of the electric field enhancement Q E the the magnetic field enhancement factor Q H against the wavenumber κ, respectively. We observe that both the electric and magnetic enhancement factors do reach peaks at the resonant frequencies

29 ELECTROMAGNETIC FIELD ENHANCEMENT 9 Enhancement factor Q E Enhancement factor Q H Wavenumber Wavenumber Figure 3. The field enhancement factor is plotted against the wavenumber κ for the PEC-PMC cavity: left the electric field enhancement factor Q E ; righ the magnetic field enhancement factor Q H. 3., which are close to κ = nπ, n =,,.... Besides the peaks at the resonant frequencies, the electric field enhancement factor Q E also displays a peak when the wavenumber κ is sufficiently small, i.e., when the wavelength of the incident field λ is sufficiently large. In this situation, the scales satisfy d λ. The electric field should have an enormous enhancement according to our theoretical result in Theorem.8. In contrast, the magnetic field enhancement factor Q H does not show a peak when the wavenumber approaches zero, which is also confirmed in Theorem.8 that the magnetic field has no enhancement when d λ. 4. PEC-PEC cavity without resonance In this section, assuming d λ, we show that there is a weak enhancement of the electric field and there is no enhancement of the magnetic field in the PEC-PEC cavity. 4.. Problem formulation. Consider the following model problem of the electromagnetic scattering by an open cavity: u κ u = in Ω, ν u = on Ω, 4. lim r r r u s iκu s = in R, where u sc = u u inc u ref in R. Here the incident field uinc = e iκx d is a plane wave propagating in the direction d = sin θ, cos θ and the reflected field u r = e iκx d with the propagating direction d = cos θ, sin θ. Let φ x = and φ n x = cos nπx n be an orthonormal basis on the interval,. Since the total field satisfies the Neumann boundary condition on the bottom and lateral sides of the cavity, we have from 4. that u can be expressed as the sum of waveguide modes: u x = α n e iβnx αn e iβnx d φ n x, x D, 4.

30 3 YIXIAN GAO, PEIJUN LI, AND XIAOKAI YUAN where the coefficients β n are defined as { κ, n =, β n = i nπ/ κ, n. If < λ/, then β n are pure imaginary numbers for all n. For each n, it can be seen if nπ/ κ, the expansion consists of two propagating wave modes traveling upward and downward respectively; if nπ/ > κ, the expansion consists of two evanescent wave modes decaying exponentially away from the bottom side and open aperture of the cavity, respectively. Using the expansion.5, we may reformulate the model 4. into the following coupled problem: u κ u = in R, u = α n e iβnx αn e iβnx d φ n x in D, ν u = on Ω \ Γ, 4.3 u x, = u x,, x u x, = x u x, on Γ, lim r r u sc iκu sc = in R. r By matching the expansion series 4. and using the continuity conditions in 4.3, we may obtain Therefore, the wave field u x = e iβ nd u,n u,n e iβnd A straightforward calculation yields x u x, = x u x, d = α n = eiβnd u,n u,n e iβnd, αn = eiβnd u,n u,n e iβnd. e iβnx eiβnd u,n u,n e iβnd e iβnx d e iβ nd u,n u,n iβ n e iβnd e iβnd eiβnd u,n u,n e iβnd It follows from 4.6 that Substituting 4.7 into 4.5 yields iβ n e iβ nd u,n u,n e iβnd x u x, = u,n = φ n x, x D. 4.4 φ n x on Γ, 4.5 eiβnd u,n u,n e iβnd e iβnd φ n x = on Γ. 4.6 u,n e iβnd. 4.7 eiβnd iβ n e iβnd e iβnd e iβnd e iβnd u,nφ n x on Γ. Given a function u H / Γ, we let u = ux, and define a DtN map on Γ : B PEC [u] = iβ n e iβnd e iβnd e iβnd e iβnd u n φ n x, 4.8 where the Fourier coefficients u n = u, φ n Γ. Following a similar proof to Lemma., we may show that the DtN map B PEC is also bounded. The proof is omitted for brevity.

31 ELECTROMAGNETIC FIELD ENHANCEMENT 3 Lemma 4.. The DtN map B PEC defined by.5 is bounded from H / Γ H / Γ, i.e., where the constant C = max By 4.8, we derive the TBC: B PEC [u] H / Γ C u H / Γ u H / Γ, { } κ sinκd cosκd, e nπ/ κ e nπ/ κ lim r x u = B PEC [u ] on Γ, < is independent of λ, d. which helps to reduce 4. into to a boundary value problem in R : u κ u = in R, ν u = on Γ \ Γ, x u = B PEC [u ] on Γ, r r u sc iκu sc = in R. The solution u in the cavity D may be obtained from 4.4 once the Fourier coefficients u,n are available note that u,n can be computed from from 4.7. To find an approximate the model problem to 4.9, we similarly drop all the nonzero wave modes and define lim r 4.9 B PEC [v] = iκ eiκd e iκd e iκd e iκd v φ x. 4. Hence we may consider the following problem: v κ v = in R, ν v = on Γ \ Γ, x v = B PEC [v ] on Γ, r r v sc iκv sc = in R, 4. where v s = v u inc u ref in R. Inside the cavity, we may approximate u with one single mode: v x = α e iκx α eiκx d φ x, 4. where the coefficients α and α are given by It follows from ν v = on Γ α = eiκd v, v, e iκd that α = α eiκd, which yields, α = eiκd v, v, e iκd. 4.3 v, = v, e iκd. 4.4 e iκd 4.. Enhancement of the approximated field. This section presents the estimates of the solution for the approximate model problem 4.. Theorem 4.. Let v be the solution of 4. and be given by 4. inside the cavity, then there exist positive constants C, C independent on λ,, d such that C κ d 3/ v L D C κ d 3/.