arxiv: v2 [physics.optics] 18 Dec 2018

Transcription

1 Extraordinar transmission through a narrow slit Jacob R. Holle and Or Schnitzer Department of Mathematics, Imperial College London, SW7 2AZ, UK arxiv: v2 [phsics.optics] 8 Dec 208 We revisit the problem of extraordinar transmission of acoustic (electromagnetic) waves through a slit in a rigid (perfectl conducting) wall. We use matched asmptotic expansions to stud the pertinent limit where the slit width is small compared to the wall thickness, the latter being commensurate with the wavelength. Our analsis focuses on near-resonance frequencies, furnishing elementar formulae for the field enhancement, transmission efficienc, and deviations of the resonances from the Fabr Pérot frequencies of the slit. We find that the apertures near fields pla a dominant role, in contrast with the prevalent approximate theor of Takakura [Phsical Review Letters, (200)]. Our theor agrees remarkabl well with numerical solutions and electromagnetic experiments [Suckling et al., Phsical Review Letters, (2004)], thus providing a paradigm for analzing a wide range of wave propagation problems involving small holes and slits. Introduction. The term extraordinar transmission originated in electromagnetism [] to describe enhanced transmission of wave energ through small apertures via excitation of surface plasmons, spoof plasmons and localized resonances, with analogous effects identified in acoustics [2] and water waves [3]. An elementar example, that has been widel studied experimentall and theoreticall, is transmission of TM-polarized electromagnetic waves through a single slit in a metallic wall whose thickness is comparable to the wavelength and large compared to the slit width [4 8]. When the metal is perfectl conducting then the problem is analogous to (lossless) acoustic transmission through a slit in a rigid wall [2, 9]. In these idealized settings enhanced transmission is attributed to excitation of standing waves in the slit, or Fabr Pérot resonances, which leak energ via diffraction at the slit ends; with decreasing slit width, the onresonance transmission efficienc and slit-field magnitude are enhanced while the transmission peaks approach the standing-wave frequencies. Takakura [4] was the first to put forward an approximate theor of ideal transmission through a narrow slit, based on an ad hoc truncation of an exact mode-matching scheme. Takakura s ke result, a simple closed-form approximation for the deviations of the transmission peaks from the Fabr Pérot resonances, was found to agree poorl with electromagnetic experiments [7]; the discrepanc could not be solel attributed to material losses or a skin effect as it persisted even for slits wide enough for the metal to be considered perfectl conducting (still exceedingl narrow relative to the wall thickness) [0]. More recent attempts to derive rigorous approximations for ideal transmission through a single slit [8, ] have not provided explicit formulae or phsical insight into the discrepanc. In this Letter we sstematicall develop an approximate theor of extraordinar transmission through a narrow slit. Our theor is based on an asmptotic analsis of the ideal transmission problem (described in the language of acoustics, for convenience) in the pertinent limit where the slit width is small compared to the wall thickness, the latter being comparable to the wavelength. We demonstrate that our results, which include asmptotic formulae for the transmission efficienc, field enhancement and frequenc shifts, are in excellent agreement with numerical and experimental data; moreover, we expose the logarithmic error rendering Takakura s formula inaccurate and point out the (ver common) flawed phsical assumption at the basis of that approximation. The notabl simple form ofour theorfollows from the manner in which we treat the disparate length scales and distinguished frequenc regimes in the problem. Namel, rather than attempting to find a single approximation for the wave field that is uniforml valid everwhere and for all frequencies, we use the method of matched asmptotic expansions [2] to sstematicall decompose the phsical domain into asmptoticall overlapping regions (representing the slit, the external regions and transition regions near the apertures) and separatel consider offand near-resonance frequencies; our analsis is thereb distinct from recent asmptotic analses of related transmission problems [3, 3 5]. Formulation. Consider acoustic transmission of a normall incident plane wave (angular frequenc ω, wave speedc)throughaslitofwidth2lhinarigidwallofthickness l. We suppress the time dependence e iωt in the usual wa and adopt a dimensionless formulation where lengthsarenormalisedbl andϕdenotesthevelocitpotential normalised b the amplitude of the incident wave. Figure shows a dimensionless schematic of the problem and defines the Cartesian coordinates (x, ). In the fluid domain ϕ satisfies the Helmholtz equation 2 ϕ+ω 2 ϕ = 0, () where we define the dimensionless frequenc Ω = ωl c. (2) On the walls, impermeabilit implies that the normal derivative vanishes, ϕ = 0. (3) n

2 2 Φ U(), where U() satisfies rigid/ PEC d 2 U d 2 +Ω2 U = 0 (6) normall incident plane wave 2h with solution (see figure 2(b)) 2/Ω x U() = Acos(Ω)+Bsin(Ω). (7) To determine the constants A and B we asmptoticall match the slit field to the outer domains; noting that (5) tends to a limiting value as /2 and that ϕ + is asmptoticall small, we find FIG.. Schematic of the dimensionless transmission problem. In addition, the scattered field ϕ ϕ (i), where ϕ (i) = e iω is the incident wave, must propagate outwards at large distances from the slit. (In the electromagnetic analog mentioned in the introduction, the same dimensionless problem holds with ϕ replaced b the scaled out-of-plane magnetic-field component.) In what follows we shall be interested in the resonant peaks of the transmission efficienc η, defined as the ratio of the acoustic (or electromagnetic) power transmitted through the slit to that transmitted through the same cross section in the absence of a wall: η 2hΩ h h Im { ϕ ϕ } dx, (4) where the asterisk denotes complex conjugation and the integrand is evaluated for fixed such that < /2. Analsis off-resonance. Let us first, naivel, consider the narrow-slit limit h holding Ω fixed. The slit openings reduce to points, as depicted in figure 2(a). Assuming that the slit potential is comparable in magnitude to the incident wave, as we shall later verif, diffraction from the slit is negligible in the domain left of the wall; the potential there, sa ϕ, is accordingl ϕ e iω +Re iω, (5) where the reflection coefficient R = e iω follows from (3) and denotes a relative asmptotic error of algebraic order, i.e., scaling with some positive power of h. Right of the wall, there is no incident nor reflected wave, onl weak diffraction from the slit; the potential there, ϕ +, is accordingl algebraicall small. The above approximations are outer expansions, namel with position held fixed. We next derive an inner approximation for the wave field inside the slit b considering the limit h 0 with the stretched transverse coordinate X = x/h fixed. Defining the slit field ϕ(x,) = Φ(X,), for < /2, it is straightforward to show from a regular perturbation of () and (3) that A = e iω/2 sec Ω 2, B = e iω/2 csc Ω 2. (8) Based on this solution the dimensionless flux densit in the direction, ϕ/, is O() in the slit; at the right end of the slit the corresponding net flux is 2hU (/2). This flux is approximatel conserved on the small scale of the aperture, impling that the outer potential ϕ + is given b a transmitted clindrical wave [6] ϕ + hqh 0 (Ωr + ), (9) where H 0 is the zeroth-order Hankel function, r ± = x2 +( /2) 2 and H 0 (t) 2i Q = 2ie iω/2 ΩcscΩ. (0) Here and later we use the asmptotic relation ( ln t 2 +γ E ) + as t 0, () wherein the algebraic error is in t and γ E is the Euler Mascheroni constant [7]. Our fixed-frequenc approximation gives no pronounced enhancement, viz. η = O(), except near the Fabr Pérot frequencies Ω = m, m =,2,..., (2) where (8) and (0) are singular; as Ω Ω the slit field U() diverges in amplitude and approaches a standingwave solution vanishing at = ±/2. The singular nature of the fixed-frequenc limit is evident on phsical grounds: as diffraction from the slit ends is neglected(and material losses have been discarded from the outset), the limiting standing waves constitute resonant modes. For Ω Ω = O(h), however, our present approximation suggests that the amplitude of the slit field is O(/h) and the diffracted waves are O(), comparable to the incident and reflected waves. Our asmptotic analsis therefore needs to be modified in this regime. Analsis near resonance. In light of the above, we now define Ω Ω Ω 2hΩ (3)

3 3 (a) (d) Y + E Y + = ( /2)/h Z x X = x/h q ζ A B X (b) /2 X = x/h /2 FIG. 2. Asmptotic outer regions (a) and inner slit region (b). Near resonances, the slit field is approximatel a Fabr Pérot mode enhanced b an O(h ) factor (c) with the transition regions (d) close to the apertures plaing a dominant role. and in what follows consider the near-resonance limit h 0 with Ω (and m) fixed. In this modified limit we anticipate the diffracted waves in the outer regions to be O(), i.e., and (c) ϕ e i Ω +e i Ω(+) +Q H 0 ( Ωr ) (4) ϕ + Q + H 0 ( Ωr + ), (5) where the coefficients Q ± remain to be determined. The slit potential is accordingl amplified, Φ(X,) h U(), (6) where U() satisfies equation (6) with Ω replacing Ω, whereas the relativel negligible magnitude of the outer potentials (4) and (5) implies the boundar conditions U(±/2) = 0. The resulting homogeneous problem for U(), b construction, has as non-trivial solutions U() = A { cos( Ω) sin( Ω) 2 }, (7) where A is a prefactor and henceforth the upper element of an arra corresponds to odd m (even standing wave) and the lower to even m (odd standing wave). The analsis of the slit region can be readil extended one further algebraic order b writing where V() satisfies Φ(X,) h U()+V(), (8) d 2 V d 2 + Ω 2 V = 4Ω Ω2 U. (9) A B C D E source p C D FIG. 3. Conformal mapping of the left half (X < 0) of the transition region into the upper half of an auxiliar ζ plane. (It is straightforward to verif that the slit potential deviates from a unidimensional profile onl at higher algebraic order.) Multipling equation (9) b U and subtracting the conjugate of equation (6) (with Ω replacing Ω) multiplied b V, followed b integration along the slit using U(±/2) = 0, gives ] /2 /2 [V du = 4Ω Ω2 U 2 d. (20) d /2 Substituting (7) we find A = im 2Ω Ω /2 { } i(v)/2 +i(v) /2. (2) (V) /2 (V) /2 As in the off-resonance analsis, the leading fluxes at the slit ends can be matched with the corresponding diffraction terms in the outer regions; we thereb find Q ± = i m ΩA { }. (22) i Unlike in the off-resonance analsis, however, the slit and outer potentials cannot be directl matched, since the leading-order outer potentials are singular as r ± 0. It is therefore necessar to consider transition regions near the apertures defined b intermediate limits with X = x/h and Y ± = ±( /2)/h fixed (we also define R ± = r ± /h); under these scalings, the slit appears infinite and the boundaries are Y ± = 0 for X > and X = ± for Y ± < 0 (see figure 2(d)). It is evident from the form of the slit and outer expansions that the potential is O() in the transitions regions; thus let ϕ(x,) T ± (X,Y ± ), where T ± satisfies Laplace s equation the Neumann boundar condition 2 T ± X T ± = 0, (23) Y ±2 T ± N = 0, (24)

4 FIG. 4. Numerical(black line) and near-resonance asmptotic (thick blue lines) transmission efficienc, for h = 0.0, plotted as a function of the normalized frequenc FIG. 5. Relative frequenc shifts from Fabr-Pérot values (numerics: smbols; asmptotics: solid lines) as a function of h for modes m = (circles) and m = 3 (squares). The dashed lines depict Takakura s approximation (35). where N denotes Y ± for the external walls and X for the inner slit walls, as well as asmptotic matching with the outer and slit regions. In the off-resonance limit the O() outerandslitpotentialswereregularhencet ± were constants; in the present near-resonance limit, however, matching with the outer expansions gives T ± iq ±2 [ lnr± +Q ± + 2i ( ln Ωh )] 2 +γ E +( )e i Ω/2 +o() as R ±, (25) while matching with the slit expansion gives T ± iq ± Y ± +(V) ±/2 +o() as Y ±, (26) where the leading term in (26) is written in terms of Q ± using (22). The Laplace problemfor T + (similarlt ) is uniquel solvable, up to an arbitrar additive constant, with just the leading terms in the far-field conditions (25) and (26) prescribed; that unique solution therefore determines the differences between the O() constant terms in those conditions. To treat both Laplace problems simultaneousl we note that T ± /iq ± solves a joint canonical problem where the difference ( ) T ± β lim Y ± iq ± Y± lim R ± ( T ± iq ± 2 lnr± ) (27) is a pure number, which depends onl on the geometr of the aperture. In the present geometr we obtain β b defining the complex variable Z = X +iy + and considering the conformal mapping [8] Z = 2i (+ζ)/2 + i log (+ζ)/2 (+ζ) /2 + (28) from the upper half-plane of an auxiliar complex variable ζ = p + iq to the left half (X < 0) of the XY + domain (see figure 3). The form of the solution is evident in terms of the auxiliar variable: T + iq + = ln ζ +const. (29) Using (29) and the limits of (28) as ζ 0 and ζ, definition (27) ields β = 2 (ln 4 ). (30) With β determined, the complex amplitude of the slit wave is found from (2), (22) and (25) (27) as { } / Ω i A = ( ). (3) Ω 2 ln 2 Ωh +γ E +i From this main result we find the transmission efficienc, η h Ω A 2, (32) and field enhancement in the slit, A /h, in terms of the amplitude A 2 = (m) 2 + [ Ω 2 (ln(2mh)+γ E ) ] 2, (33) a Lorenzian with a peak magnitude decreasing with m; the frequenc deviations f, of the resonant peaks from the Fabr Pérot frequencies f m = mc/(2l), are f f m 4h (lnh+ln(2m)+γ E ). (34) Discussion and concluding remarks. It is illuminating to compare our asmptotic theor with numerical calculations, experimental data and existing approximations. To begin with, we have validated our theor

5 (ignoring terms of quadratic order in h). Takakura s approximation is seen to be in poor agreement with the numerical data. Note that (34) and (35) disagree at O(h) which is practicall comparable to the leading O(h ln h) term[20]; our analsis suggests that the logarithmic error in Takakura s approximation is a result of his disregard of the details of the wave field close to the apertures, which are seen to be important specificall near resonance. Finall, in figure 6 we replicate the comparison carried out b Suckling et al. [7] between their electromagnetic measurements and Takakura s approximation (35), with our asmptotic prediction (34) overlaid. Our approximation is seen to be in excellent agreement with the data in an intermediate range of slit widths: narrow enough so that h and wide enough such that the metal can be modelled as a perfect conductor to a good approximation. In conclusion, our asmptotic analsis provides a simple, accurate and phsicall representative theor of ideal extraordinar transmission through a single narrow slit, which in particular improves agreement with experiments and highlights the importance of aperture effects close to resonance. As the slit width is reduced, holding the wall thickness and frequenc fixed, intrinsic (material) losses ultimatel become dominant over radiation damping. This is evident in the electromagnetic experimental data of Ref. [7] shown in the inset to figure 6 and more pronouncedl in recent acoustic experiments [9, 2]. It is therefore desirable to extend the present theor to incorporate losses, which are of an inherentl different nature in the electromagnetic and acoustic cases [7, 9, 22]. Our asmptotic approach ma offer simplification in both of these non-ideal regimes, which are fundamental to the practical design of acoustic and photonic metamaterials. OS is grateful for support from EPSRC through grant EP/R04458/. 5 FIG. 6. Comparison with experiments of Ref. [7] of electromagnetic transmission through a slit in a 9.58mm thick Aluminium plate (triangles using plot digitizer). Main plot shows shifts of resonance frequencies from the second to ninth Fabr-Pérot mode for a 500µm wide slit. Inset shows the resonance frequencies for mode m = 9 for varing slit width. In both we depict our asmptotic prediction b a solid line and Takakura s approximation (35) b a dashed line. against a mode-matching scheme solving the ideal transmission problem exactl [9]; in particular, figures 4 and 5 present excellent agreement between our on-resonance asmptotic predictions for η and f [cf. (32) (34)] and the corresponding numerical values. In figure 5 we have added a dashed line depicting Takakura s prediction for the frequenc shifts [4]: ( ) f = 4h ( lnh+ln(m) 3 ) (35) f m Ref. [4] 2 [] T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wolff, Nature 39, 667 (998). [2] J. Christensen, L. Martin-Moreno, and F. J. Garcia- Vidal, Phs. Rev. Lett. 0, 0430 (2008). [3] D. V. Evans, R. Porter, and J. R. Chaplin, in Proc 33rd International Workshop on Water Waves & Floating Bodies. [4] Y. Takakura, Phs. Rev. Lett. 86, 560 (200). [5] F. Yang and J. R. Sambles, Phs. Rev. Lett. 89, (2002). [6] F. J. García-Vidal, H. J. Lezec, T. W. Ebbesen, and L. Martin-Moreno, Phs. Rev. Lett. 90, 2390 (2003). [7] J. R. Suckling, A. P. Hibbins, M. J. Lockear, T. W. Preist, J. R. Sambles, and C. R. Lawrence, Phs. Rev. Lett. 92, 4740 (2004). [8] J. Lin and H. Zhang, SIAM J. Appl. Math. 77, 95 (207). [9] G. P. Ward, R. K. Lovelock, A. R. J. Murra, A. P. Hibbins, J. R. Sambles, and J. D. Smith, Phs. Rev. Lett. 5, (205). [0] Earlier, less accurate, experiments [5] seemed to show good agreement with Takakura s approximation. [] P. Jol and S. Tordeux, Multiscale Model. Simul. 5, 304 (2006). [2] E. J. Hinch, Perturbation methods (Cambridge universit press, 99). [3] D. V. Evans and R. Porter, Quart. J. Mech. Appl. Math. 70, 87 (207). [4] O. Schnitzer, Phs. Rev. B 96, (207). [5] R. Porter, Scattering in a waveguide with narrow side channels, Tech. Rep. [6] Higher-order fundamental solutions must be discarded as their algebraic singularit contradicts the order of magnitude of the slit potential. [7] M. Abramowitz and I. A. Stegun, Handbook of mathematical functions (Dover New York, 972). [8] R. V. Churchill and J. W. Brown, Complex variables and applications (McGraw-Hill, New York, 960). [9] Please contact authors for the supplementar information. [20] The leading O(hlnh) term was also found b Lin et al. [8]. Their expression for the O(h) term, however, includes a constant which was not determined. [2] G. P. Ward, A. P. Hibbins, J. R. Sambles, and J. D. Smith, Phs. Rev. B 94, (206). [22] M. Molerón, M. Serra-Garcia, and C. Daraio, New J. Phs. 8, (206).