Issue 
EPJ Applied Metamaterials
Volume 1, 2014



Article Number  2  
Number of page(s)  16  
DOI  https://doi.org/10.1051/epjam/2014003  
Published online  13 June 2014 
https://doi.org/10.1051/epjam/2014003
Review
Deep subwavelength FabryPerot resonances
^{1}
Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, P.R. China
^{2}
Department of Applied Physics, Nanjing Tech University, Nanjing
210009, P.R. China
^{*} email: cphuang@njtech.edu.cn
Received:
10
April
2014
Accepted:
1
May
2014
Published online: 13 June 2014
Confinement of light by subwavelength objects facilitates the realization of compact photonic devices and the enhancement of lightmatter interactions. The FabryPerot (FP) cavity provides an efficient tool for confining light. However, the conventional FP cavity length is usually comparable to or larger than the light wavelength, making them inconvenient for many applications. By manipulating the reflection phase at the cavity boundaries, the FP cavity length could be made much smaller than the wavelength. In this review, we consider the subwavelength FP resonance in a plasmonic system composed of a slit grating backed with a ground plane, covering the spectral range from microwave to THz and infrared regime. For very narrow slit width and spacer thickness, a typical zeroorder and deep subwavelength FP resonance in the metallic slits can be strongly induced. Moreover, due to the subwavelength FP resonance, greatly enhanced electromagnetic pressure can also be induced in the system. The sign and magnitude of the electromagnetic pressure are dominated by the field penetration effect in the metal as well as the field enhancement in the FP cavities. The effect promises a variety of potential applications, such as detecting tiny motions and driving the mechanical oscillations.
Key words: Plasmonics / Subwavelength / FabryPerot cavity / Resonance
© C.P. Huang & C.T. Chan, Published by EDP Sciences, 2014
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Waves in the periodic media are subject to the Bragg scattering, giving rise to a band structure. Such effect presents in both the electronic and photonic or phononic crystals [1–4], where the lattice constants are of the order of the wavelength. For waves with a long wavelength, the manipulation of wave motion typically requires a bulky device. On the other hand, by utilizing the strong resonance effect, waves can also be manipulated by the subwavelength or even deepsubwavelength objects. One example is the piezo electric superlattice composing of periodic ferroelectric domains [5–7], where the coupling between the electromagnetic wave and superlattice vibration may lead to a polaritonic bandgap in the longwavelength region. Another example is the socalled “locally resonant sonic material” [8, 9], which can block acoustic waves with wavelength much larger than the characteristic size of the crystal. Recently, plasmonic materials have attracted much research interests [10, 11]. It is well known that the plasmonic resonance of nanostructures can be employed to confine light or guide the flow of light [12, 13]. Moreover, the electric and/or magnetic resonances in metamaterials with the subwavelength components have been demonstrated to enable exotic wave manipulation properties such as the negative refraction and superlensing [14–16]. The concept of resonance may provide a way for miniaturizing the devices and realizing novel optical effects with subwavelength components.
There are many different resonators which can be employed to realize optical resonances, and examples include metallic nanoparticles, metallic or dielectric microdisks, and photoniccrystal cavities [12, 17–19]. Here we would like to discuss two typical optical resonators: the inductancecapacitance (LC)like resonator and the the FabryPerot (FP)like resonator. A wellknown example of LC resonator is the metallic splitring resonator (SRR), which plays a crucial role in metamaterial design [14, 20–22]. The LC resonator has an effective inductance and capacitance, inducing the free electrons to oscillate like mechanical harmonic oscillators [23]. The resonance effect can be extended from the microwave band to THz and infrared range by reducing the resonator sizes [21, 22]. For example, the optical properties of nanometals or flat metal surfaces have been studied using the LC model [24–27]. However, the scaling of resonance frequency with size will break down in the optical regime [28], due to the dominance of internal inductance. The LC resonance can be excited by the electric or magnetic field of light or by the neighboring resonators via the nearfield coupling. The coupling between the LClike resonance and electric dipole mode can mimic the electromagneticallyinduced transparency [29].
In contrast, the FPlike resonator is an optical cavity based on the electromagnetic interference effect [30, 31]. A conventional FP resonator consists of an air cavity of length L sandwiched between two parallel reflecting mirrors. The constructive or destructive interference of forward and backward propagating waves results in the formation of a standing wave pattern in the cavity. The FP cavity has been widely used in applications such as the laser cavity, signal filtering, directive emission, and optomechanical coupling [31–33]. However, compared with the LClike resonator which is frequency selective and can be deep subwavelength in size, the FPlike resonator operates with the wavelength that scales with the cavity length. The resonance condition of a conventional FP cavity requires that the total phase shift of one round trip equals a multiple of 2π, i.e., 4πnL/λ + φ_{1} + φ_{2} = 2πm (here, n is the index of cavity filling medium, φ_{1} and φ_{2} are the reflection phases, and m is an integer). Generally speaking, the reflection phase at the interface between two media has two possible values 0 or π, depending on the ratio of medium refractive index. Thus, it can be inferred that the cavity length must be comparable to or larger than the resonance wavelength. This makes the device unsuitable for applications that require small formfactor components.
In recent years, the FPlike effect in plasmonic circumstances has received much attention [34–41]. It was found that the FPlike resonances in the subwavelength apertures of a metal film can increase the transmission of light [34, 35]. Moreover, a simple plasmonic FP cavity can be constructed by using the plasmonic mirrors and/or employing the plasmonic confinement [36, 38]. Although the quality factor is compromised due to the Ohmic loss, the presence of plasmonic effect can squeeze the light wavelength and reduce the mode area drastically. Consequently, the combination of FP resonance and surfaceplasmon resonance yields an efficient localization of light and a significant enhancement of electro magnetic fields in the cavity [36]. By inserting the metallic nanostructures into an FP cavity, one can also observe the strong coupling between the optical mode and localized plasmon mode [40, 41]. These effects may be used for increasing the optical nonlinearity or enhancing the cavitybased plasmonic sensing. Although a simple plasmonic FP cavity may provide a solution for the miniaturization of device, the fundamental limit on the cavity length cannot be circumvented in this way.
It was noted by Trentini that the FP resonators could be created with a partially reflecting reflector that induces a phase shift different from 0 or π [42]. This causes a slight variation of cavity length with respect to the standard value. Although the effect was not very obvious in the original structure, significant advances can be achieved now thanks to the modern metamaterial designs. In 1999, Sievenpiper et al. suggested that a mushroomlike structure, consisting of an array of metal plates connected to a back metal sheet, can act as a highimpedance surface [43, 44]. The reflection phase of such a surface may experience a transition with frequency from π to –π by crossing through zero at the LC resonance frequency. With this unusual reflection property, a subwavelength FP cavity with cavity length one order of magnitude smaller than the resonance wavelength has been designed and constructed operating in the microwave band [45]. By going from one to all three dimensions [46], a cubic cavity with each dimension about a quarter of the resonance wavelength was also realized. The FP cavity length can be further squeezed by using the double metamaterialbased reflectors. The reflection phases at both reflector surfaces decrease gradually with the frequency; in a narrow frequency range, the two phases cancel each other partially, allowing for an ultrathin cavity with the length as small as λ/60 [47]. With such cavities, an ultrasmall FP antenna, i.e., a subwavelength FP cavity coupled with an internal antenna, can provide a highly directive emission at the resonance [45, 47, 48].
The subwavelength FP cavities constructed with the multilayer metamaterial reflectors are of the complicated geometry and thus not easy to realize at high frequencies. In this review, we discuss the designs and properties of a special kind of subwavelength FP resonator, which is composed of a metallic slit grating backed with a ground plane. Due to the unique phase shift, the metallic slits within the structure may act as “deep subwavelength” FP cavities [49], where the resonance wavelength can be two orders of magnitude larger than the cavity length. Compared with the previous designs, our FP cavity is openended without the microstructured physical reflectors at the cavity openings. Thus the cavity can be reduced with the nanoscale sizes and can work at high frequencies. The deep subwavelength FPlike resonance is also accompanied by many interesting physical properties [50], such as enhanced absorption and enhanced electromagnetic pressure.
A systematic discussion will be presented in the following sections. In Section 2, several mechanisms for the nontrivial phase shift at the reflecting interface are analyzed. This could be helpful for understanding Section 3, where the reflection phase and resonance condition of the proposed subwavelength FP (SFP) resonator are presented. In Section 4, we discuss the SFP resonance in the microwave band. The enhanced electromagnetic pressure due to the SFP resonance was emphasized. The effect is extended in Section 5 to the higher frequency THz and infrared regime. The role of the field penetration effect in the metal is also suggested. Some additional comments are provided in Section 6 and the findings are summarized in Section 7.
2. Reflection phase at the boundaries
2.1 Phase shift due to resonant impedances
For materials involving both electric and magnetic response, the reflection coefficient of light at normal incidence is determined by the impedance of the interface medium,(1)
Here, is the impedance of a bulk medium.
For natural materials, the impedance and the phase shift are usually monotonous functions of frequency. However, by using the designed artificial materials or metamaterials with the subwavelength constituents, the electric and magnetic response can be engineered at will. For example, an array of metallic SRRs shows a magnetic resonance [20–22] near which the permeability can be enhanced, thus giving rise to higher impedance. Moreover, a negative permeability may also be induced in a frequency gap. On the other hand, an array of metallic nanoparticles has an electric response [23, 51, 52]. In the effectivemedium model, the impedance is then modulated depending on the permittivity which can be positive, negative, or zero in the specific frequency range. The negative (or complex) permeability or permittivity due to the resonance (or absorption) gives complex impedance and thus changes the phase shift [41].
Interestingly, an array of subwavelength holes (the complementary structure of the particles) in a semiinfinite thick and perfectelectric conducting (PEC) metal, as shown in Figure 1a, may act as an effective medium with high impedance. Following Pendry et al. [53, 54], the perforatedmetal medium has an effective permeability μ_{eff} = 8a ^{2}/π ^{2} d ^{2} and permittivity ε_{eff} = (π ^{2} d ^{2}/8a ^{2})(ε_{h} − π ^{2} c ^{2}/a ^{2} ω ^{2}) in the longwavelength approximation. Here, a is the width of the square hole, d is the lattice constant, c is the light speed, and ε_{h} (ε_{h} ≫ 1) is the permittivity of hole filling medium. The effective permittivity is of the Drudemodel form with a “plasma” frequency of the waveguide cutoff frequency. The effective impedance shows a resonance at the hole cutoff. Below the cutoff frequency, the impedance is imaginary and causes a phase shift. The effect can also be understood in terms of the evanescentwave effect, as indicated below. Figure 1d shows an example for this case.
Figure 1. Reflection of light at various interfaces: (a) A semiinfinite PEC metal milled with subwavelength square holes behaves like an effective medium (λ ≫ d) with a resonant impedance (the hole width is a, the period is d, and the permittivity of hole filling medium is set as ε_{h} = 25); (b) A semiinfinite PEC metal milled with periodic slits (the slit width is a and the period is d) can induce the diffracted evanescent waves at the interface (λ ~ d); (c) Light reflection by a dielectric slab (with a thickness h) sandwiched between two semiinfinite media. The reflection efficiency and phase shift for the case (a)–(c) are shown in (d)–(f) respectively. We set in (d) and (e) a = d/2 and in (f), n_{1} = 1, n_{2} = 2, and n_{3} = 1.5. 
The device is bulky using a bulk medium as the reflectors. In the longwavelength limit, a flat surface textured with periodic metallic features can be treated as an effective medium, which may present unusual effective surface impedance [43]. Here, the surface impedance relates the tangential electric field, including both the incoming and reflecting components, to the magnetic field at a boundary, thus determining the reflection of light directly. For example, the mushroomlike structure can be modeled as an equivalent resonant LC circuit, which has a surface impedance of Z_{s} = iωL/(ω ^{2} LC − 1) [43]. This impedance is inductive in the low frequencies, capacitive in the high frequencies, and exhibits a larger value at the LC resonance. Accordingly, the reflection phase of the surface decreases continuously with the frequency from π to –π. Such a structure may also be treated as an effective medium with a resonant permeability [55].
2.2 Phase shift due to evanescent waves
For transversemagnetically (TM) polarized incident light, the magnetic Fresnel’s coefficient of reflection between the dielectric medium 1 and 2 (the permittivity is ε_{1} and ε_{2}, respectively; the permeability is μ_{1} = μ_{2} = 1) reads [56]:(2)
Here, represents the wavevector component along the normal of interface, k_{//} is the parallel wavevector component along the boundary, and k_{0} is the wavevector in free space.
Generally, the reflecting wave presents a trivial phase shift 0 or π. When the normal wavevector is imaginary, however, the reflection phase will deviate from 0 or π. This case can be encountered in the total reflection when the incident angle is larger than a critical value. Alternatively, the normal wavevector is also imaginary when the metal with a negative permittivity constitutes one of the interface media. The imaginary wavevector means that the wave is evanescent and decays exponentially from the interface. Accordingly, a phase shift proportional to the decay length, , will be resulted. Hence, the evanescent wave will play a crucial role in tailoring the reflection phase at the boundaries although the evanescent wave cannot carry the energy away from the interface. When the transmitted wave is propagating or has a real momentum in the normal direction, the reflection phase will be 0 or π; when the transmitted wave is evanescent or has an imaginary normal momentum, the reflection phase is then changed.
The evanescent wave can also present in many other systems, such as the near field of a subwavelength nanostructure, the surface of a periodic structure, etc. [57–59]. For light incidence upon a metal surface perforated with the periodic subwavelength apertures, the propagating and evanescent waves can be excited simultaneously and the reflection phase may be changed depending on the weight of the evanescent components. We take a PEC and semiinfinite thick metal, which is milled with periodic subwavelength slits, as an example (see Figure 1b; the slit depth is also semiinfinite). By expanding the fields on the incident side as and that in the slits as , the reflection coefficient for the TM polarization can be obtained with the boundary conditions as:(3)
Here, G_{m} = 2πm/d, g_{m} = sin c(mπa/d), d is the lattice period, and a is the slit width. Note that the evanescent waves on the incident side are marked by the diffraction orders with G_{m} ≥ k_{0} and that the evanescent effect on the transmission side can be accounted for by using a complex fillingmedium index n_{s} = n′ + in″.
Equation (3) suggests that the reflection phase at the periodic interface is still dominated by the propagating or evanescent characteristics of the waves. If we neglect the evanescent components on both the incident and transmission sides, that is, n_{s} is real and G_{m} < k_{0}, the reflection coefficient will be real, indicating a common phase shift of 0 or π. However, the presence of evanescent effect on either side may contribute to a deviation of phase shift. For example, Figure 1e presents the calculation results for n_{s} = 1 and a/d = 0.5, where a nonzero phase shift (due to the evanescent effect on the incident side) can be seen clearly.
In the longwavelength region, the semiinfinite structure may be approximated by a homogeneous medium [54], which has an effective permittivity and a permeability μ_{z} = a/d. Consequently, the impedance of the effective medium is Z = a/n_{s}d and the reflection coefficient becomes . This result can be recovered with the equation (3) by neglecting the evanescent terms. Thus, the phase variation due to the diffracted evanescent effect will be excluded in the simplified model. This also implies that one should be careful to use the effectivemedium model when dealing with the reflection phase at the boundaries.
2.3 Phase shift due to resonant cavities
The reflection phase at the boundary of a slab can also be induced, which is correlated with the oscillation of waves in the slab. Consider a dielectric slab of the permittivity ε_{2} and thickness h sandwiched between two semiinfinite medium 1 (ε_{1}) and 3 (ε_{2}) (see Figure 1c). The TMpolarized light is incident with the incident angle θ upon the slab and transmits from the medium 1 to medium 3. The reflection coefficient of the slab can be expressed as follows [56],(4)
Here, r_{12} (or r_{23}), which is determined according to the form of equation (2), is the Fresnel’s reflection coefficient between the medium ε_{1} and ε_{2} (or ε_{2} and ε_{3}); is the normal wavevector component in the slab.
Equation (4) shows that the reflection phase at the boundary is modulated by the dielectric slabcavity effect. The phase shift is generally not 0 or π even though there is no evanescent mode in the structure. To see this point more clearly, we consider the case that n_{1} = 1, n_{2} = 2, n_{3} = 1.5, and θ = 0. The reflection efficiency and phase shift as a function of wavelength (normalized by the optical path of slab) are shown in Figure 1f by the open and solid circles, respectively. It can be seen that the phase shift is modified relying on the wavelength or the slab thickness h. Note that, when h is much smaller than the wavelength, the effect is obvious as well. However, at the slabcavity resonance where or λ/n_{2}h = 4/m, the reflection phase will go back to zero.
Similar cavitydependent phase shift at the boundaries can also be present in systems such as a metal film (of finite thickness) perforated with the subwavelength slits or holes [34, 35], a thick metal milled with the periodic grooves or holes of finite depth [54, 60, 61], etc. In the case of holes, the wavelength should be smaller than the cutoff wavelength. In these cases, the evanescent modes due to the periodic textures can also generate a contribution to the phase shift which is included in the Fresnel’s coefficient r_{ij}. This evanescentmode effect will mix with the slabcavity effect to induce a total reflection phase at the boundaries.
3 SFP resonators based on metallic slits
With the modulated reflection phase, the FP cavity length could be reduced. However, for a closeended FP cavity defined by two physical reflectors, the lateral size of cavity is usually much larger than the wavelength. It is interesting to find that single slits or slit array cut into a metal film can also behave like the common FP cavities [34, 35, 38, 62, 63], as long as the slit length is large enough. The effect originates from the fact that the slit mode couples with the surface charges on the metal walls and two effective reflection interfaces can be introduced because of the discontinuities at the slit openings. Thus, the metallic slits may act as the openended FP cavities which have no physical reflectors. Moreover, the width of slits can be squeezed to a deep subwavelength scale.
Compared with the FP resonance with perfect reflecting boundaries, the FPlike resonance in the openended slit cavities only shows a slight wavelength shift [34, 38]. To realize the deep subwavelength effect, a modified sandwiched plasmonic structure can be used (see Figure 2). The sandwiched structure is usually composed of a thin metallic pattern (e.g., a periodic array of wellseparated disks, patches, or stripes) and a back ground plane, separated by a dielectric spacer [64–73]. In such structures, horizontal or vertical resonances in the dielectric spacer can be excited, leading to enhanced electromagnetic fields and absorption. For the slit grating with a narrow slit width and a larger thickness, however, the vertical metallic slits may behave as the FPlike cavities [49, 74, 75].
Figure 2. Schematic view and physical parameters of the plasmonic system under study. The slit grating (with a narrow slit width) and ground plane are separated by a thin dielectric spacer. The vertical narrow slits support a deep subwavelength FPlike resonance. Such resonance can be excited by a TMpolarized incident light. 
According to the Section 2, the upper periodic interface provides the diffracted evanescent waves in this composite cavity structure, which modify the reflection phase at the top slitcavity openings. The dielectric spacer sandwiched between the metal grating and the ground plane acts as a thin slab cavity, which incorporates both the cavity and evanescent effects. Thus, the reflection phase at the bottom slitcavity openings can also be modulated.
To determine the reflection phase of slit mode at the slit openings, one can write the fundamental slit mode in the slits and determine the reflection coefficients with the boundary conditions. For very narrow slits, the fundamental slit mode, consisting of downward and upward components, can be expressed as , where q_{0} = k_{0}n_{h} is the propagation constant and is the effective index (ε_{M} is the permittivity of the metal) [50]. By expanding the fields on the upper side as and using the boundary conditions, we have(5)
Here, w = n_{h}a/d, , and γ_{m} = mλ/d. Consequently, the slitmode electricfield reflection coefficient at the upper slit opening (y = 0), r_{1} ≡ −A/B, can be deduced as:(6)
Equation (6) confirms that the phase shift at the upper slit openings is associated with the evanescent effect.
Similarly, by expanding the fields in the dielectric spacer and using the boundary conditions [50], one obtain(7)
where and . Accordingly, the slit mode electricfield reflection coefficient at the lower slit opening (y = h), , can be derived as(8)
where is related to the spacer thickness. Thus, besides the evanescent effect, the phase shift at the bottom slit openings also benefits from a cavity effect of the spacer. When the mthorder diffraction mode is propagating, v_{m} is real and thus μ_{m} is complex, giving rise to a phase contribution due to a pure cavity effect. In addition, when the diffraction mode is evanescent, a phase contribution due to the evanescent effect will emerge (here v_{m} is imaginary, ρ_{m} ≈ −1, and μ_{m} is real). In the latter case, the phase contribution can also be modulated by the spacer cavity (dependent on t).
The SFP resonance condition can thus be written as or φ_{1} + φ_{3} + 2q_{0}h = 2πm, where φ_{1} and φ_{3} is the phase shift of the slit mode at the top and bottom slit openings, respectively. Such a SFP resonance can be excited efficiently by a TMpolarized light incident from the upper side. At the SFP resonance, the electromagnetic fields will be trapped in the slit cavities and the reflection of light is greatly suppressed. Moreover, the SFP resonance can squeeze the electromagnetic energy into the spacer, giving rise to enhanced radiation pressure in the system. These effects may find interesting potential applications in future. In the following, the SFP resonance and the resulted physical effects in the microwave band, THz and infrared range will be discussed respectively. For simplicity, throughout the rest of article, the metal is assumed to be freestanding (ε_{1} = ε_{3} = 1), and the metal permittivity is modeled with the Drude dispersion , where ω_{p} = 1.37 × 10^{16} rad/s and γ = 4.05 × 10^{13} rad/s.
4 SFP resonance in the microwave band
To illustrate the SFP resonance in the microwave band, we use the following structural parameters: d = 20 mm, a = 0.3 mm, h = 5 mm, and t = 0.1 mm. Figure 3a presents the reflection phases of the slit mode at the cavity openings, which are dispersive in the studied spectral range: the reflection phase shift at the top (φ_{1}) is positive but small while that at the bottom (φ_{3}) goes from positive to negative values by crossing a zerophase point A. The negative phase shift enables a zeroorder FP resonance with a long wavelength λ = −4πh/(φ_{1} + φ_{3}). To confirm this point, the zeroorder reflection spectrum of the whole structure at normal incidence was calculated by the equation [50]:(9)
Figure 3. The SFP resonance in the microwave band for the structure with d = 20 mm, a = 0.3 mm, h = 5 mm, and t = 0.1 mm. (a) Phase shift of slit mode at the slit openings (the arrow indicates the SFP resonance); (b) Zeroorder reflection (the line) and the absolute value of the FPlike resonance factor (the circles). The resonance around 20 mm is dominated by the grating surface modes. 
The result is shown in Figure 3b by the solid line, which exhibits a significant reflection dip at the wavelength 55.4 mm (h/λ = 9%). This reflection dip coincides with the peak of the FPlike resonance factor (see the circles), where . Thus, a subwavelength FP resonance can be induced.
4.1 Phase shift vs. structural sizes
To employ the reflection phase to modulate the SFP resonance, the dependence of phase shift on the structural parameters a and t are investigated. In Figure 4a, two slit widths a = 1.0 mm and 0.3 mm are used and the other parameters are fixed. With the reduction of slit width, φ_{1} decreases slightly. However, an obvious change occurs for φ_{3}: the slope of phase curve around the zerophase point A becomes smaller (the point A shows a weak redshift). This could push the SFP resonance to the longer wavelength. To understand this effect, we go back to equations (6) and (8), which show that the reflection coefficients of the slit mode are strongly dependent on the parameter w or duty cycle a/d. The smaller is the duty cycle, the closer the r_{1} and r_{3} to unity. That means very narrow slits exhibit excellent confinement ability, as the slit mode can be reflected almost perfectly at the openings. Simultaneously, the reflection phases of slit mode also approach zero, giving rise to more flat phase curves.
Figure 4. Dependence of phase shift on the wavelength for different slit width a and spacer thickness t: (a) t = 0.1 mm, and a = 1.0 or 0.3 mm; (b) a = 1.0 mm, and t = 0.1 or 0.3 mm. Here d = 20 mm. The correspondence between the phase shift φ_{3} (the line) and (m = 0, 1, 2) are shown in (c), where d = 20 mm, a = 0.5 mm, and t = 0.2 mm. The inset of (c) shows the zerophase point A (the circles) as a function of lattice period d (where we fix a = 1.0 mm and t = 0.2 mm), which track the position of HFP mode in the spacer closely (the line). 
The reflection phases for two spacer thickness t = 0.3 mm and 0.1 mm are shown in Figure 4b. The result is similar to that of Figure 4a, but the larger the spacer thickness, the gentler the phase variation around the point A (the position of point A is almost unchanged). This suggests that the spacer thickness can also be used to manipulate the reflection phase and the SFP resonance. Equation (8) suggests that this phenomenon is related to the spacer cavity effect, which is governed by μ_{m}. For the small spacer thickness (t ≪ λ), μ_{m} can be simplified as . By neglecting the term which contributes little to the phase in the microwave region, we have and the summation term in equation (8) becomes(10)
Equation (10) shows that each diffraction term (which is imaginary) will make a contribution to the reflection phase. Since the sum is inversely proportional to the spacer thickness t, the larger the t, the smaller the phase shift and the gentler the variation of phase curve.
To further distinguish the role of spacerconfined diffraction modes in the reflection phase, φ_{3} and μ_{m} (here m = 0, 1, and 2) as a function of wavelength have been mapped in Figure 4c. One can see that evolves linearly with the wavelength whereas and present peaks at the phase jump around the wavelength 10 mm and 20 mm. This indicates that the phase jump of φ_{3} is related to some type of resonance mode. By checking the pole of , one obtain(11)
where is the skin depth of metal. This corresponds to the mthorder “SPP resonance” in the spacer region. In the microwave band, the skin depth (~μm) is very small compared with the slit width. Thus, the resonance wavelength, mλ = ±n_{3}d(1 + δ/2a), appears very close to the Rayleigh anomaly.
In the longwavelength region, however, the role of zeroorder diffraction mode becomes significant. We note from equation (10) that the zeroorder term and the highorder terms are just of the opposite role for the reflection phase. In a particular wavelength (A point), the role of all terms may be canceled completely, giving rise to a null phase shift. This particular case satisfies the following condition:(12)
We found that the wavelength thus determined is very close to , i.e., the wavelength of horizontal FP (HFP) resonance in the spacer (beneath the metal stripes). This can be seen from the inset of Figure 4c, where the grating period d is changed (a is fixed as 0.5 mm) and the theoretically obtained positions (the circles) of the point A approach that of the HFP resonance (the line). We thus attribute the zero phaseshift at the point A to the HFP resonance in the spacer. On the left or righthand side of point A, the role of zeroorder mode will be defeated by or surpass that of highorder modes, yielding the positive or negative reflection phase.
4.2 Reflection spectrum vs. structural sizes
The influence of structural parameters on the SFP resonance can be further demonstrated by studying the reflection spectrum. The reflection spectrum as a function of wavelength and slit cavity length is shown in Figure 5a. With the increase of cavity length, the reflection minimum caused by the SFP resonance (the redcolored part) redshifts obviously. For the cavity length h = 1 ~ 10 mm, the resonance wavelength varies from 42 mm to 64 mm with h/λ = 2.4% ~ 15.6%. The increases of spacer thickness t can also redshift the SFP resonance due to the variation of reflection phase. This point is verified in Figure 5b. When t is changed from 0.1 mm to 0.7 mm, for example, the resonance wavelength will shift from 50 mm to around 90 mm. However, the reflection dip becomes less significant. In addition, we can also manipulate the SFP resonance with the slit cavity width a (see Figure 5c). The smaller is the cavity width, the larger the resonance wavelength. For a = 0.1 mm with a duty cycle as small as 0.5%, the resonance wavelength is around 76 mm (h = 5 mm) and the absorption efficiency of microwave is up to 95%. When normalized to the duty cycle of grating, the enhancement of absorption attains 190.
Figure 5. Reflection spectrum as a function of wavelength and structure sizes: (a) a = 0.5 mm, t = 0.1 mm; (b) a = 0.5 mm, h = 5 mm; and (c) h = 5 mm, t = 0.1 mm. The lattice constant is fixed as d = 20 mm. 
Experimentally, we constructed two samples to demonstrate such an effect. The samples are composed of an aluminum slit grating and a planar ground plane, which are separated by a tunable air gap (the grating covers an area of 300 mm × 300 mm). For both samples, the lattice period is d = 20 mm, slit width a = 1 mm, and the slit cavity length is 5 mm and 10 mm, respectively (see Figure 6a). The TMpolarized microwave (the magnetic field is parallel to the slits) is incident with an angle of 5°. The reflected microwave is collected and measured with a microwave network analyzer (HP8150C). Figure 6b shows the measured (the circles) and calculated (the line) reflection spectra for the sample with h = 5 mm and t = 0.15 mm. A reflection dip corresponding to the SFP resonance at 6.3 GHz (or 47.6 mm) has been observed. By tuning the air gap thickness t (from 0.15 mm to 0.9 mm), a redshift of SFP from 6.3 GHz to 3.8 GHz can be seen (Figure 6c). The spectral position of experiments agrees well with the theory. However, the measured reflection efficiency at the resonance is lower than the theoretical prediction. This can be attributed to the finite area of the samples, where the diffraction modes in the air gap can escape from the four sample ridges, thus reducing the actual reflection of microwave. In addition, for h = 10 mm and t = 0.3 mm (see Figure 6d), two reflection dips around 4.3 GHz (70 mm) and 10.8 GHz (27.8 mm) can be found. Compared with h = 5 mm, the lowfrequency dip for h = 10 mm shows an obvious redshift. We also measured the microwave reflection for different spacer thickness and the spectral position as a function of t is plotted in Figure 6e. A good agreement between theory and experiment is seen.
Figure 6. Experimental results for the microwave band: (a) Photographs of the fabricated samples. The lattice constant of the two samples is d = 20 mm, the slit width is a = 1 mm, and the grating thickness is h = 5 mm and 10 mm, respectively. The measured (the circles) and calculated (the line) reflection spectra are presented in (b) for h = 5 mm, t = 0.15 mm and (d) for h = 10 mm, t = 0.3 mm, respectively. The measured (the circles) and calculated (the lines) resonance wavelength of the two samples as a function of spacer thickness t are shown in (c), h = 5 mm, and (e), h = 10 mm, respectively. 
4.3 Electromagnetic pressure vs. structural sizes
The optical force may provide a convenient tool for manipulating small objects [76, 77]. In the cavity opto mechanics, the optical force can provide a strong coupling between the optical mode and mechanical oscillation [78–81]. Optical forces can be enhanced employing the resonance modes in the plasmonic systems [82–88]. Here we investigate the electromagnetic pressure associated with the SFP resonance which supports an efficient confinement of electromagnetic energy in the slit cavities. The coupling between the slits and spacer also enhances the electromagnetic fields in the spacer at the SFP resonance (but not at the HFP resonance). Such field enhancement may lead to an amplified electromagnetic pressure.
The numerical simulations show that the antiparallel (horizontal) currents can be excited near the spacer, which generates a repulsive or positive force [50]. Accompanying the oscillating currents, positive and negative charges accumulate on the upper or lower side of the spacer, giving rise to an attractive or negative force. The photonic repulsive pressure, (proportional to the magnetic energy density in the spacer), may be balanced by the attractive pressure (proportional to the electric energy density in the spacer) [50] due to the equipartition of energy. However, the situation will be quite different considering the following two effects.
One is the field penetration effect in the metal [89], which suppresses the magnetic field/energy due to a sharing of energy by the kinetic energy of free electrons. The larger the share of the electronic kinetic energy, the smaller the magnetic energy and the weaker the repulsive force. This effect is important in the visible frequency range. In the microwave band, the electron kinetic energy u_{K} = L_{e}i ^{2}/2, where the kinetic inductance , (k_{p} = 2π/λ_{p}, λ_{p} is the plasma wavelength, and δ is the skin depth) is small compared with the magnetic energy u_{M} = L_{m}i ^{2}/2, where the Faraday inductance L_{m} ∝ μ_{0}t (t is the spacer thickness). Thus, the kinetic energy (leading to internal inductance) will can be neglected in the microwave regime . The other effect is the field leakage effect [90], where the “leakage” of electric or magnetic energy from the spacer to the surrounding environment (especially the slit FP cavities) plays a crucial role. Hence, the subwavelength FP cavities may modulate the electromagnetic pressure strongly.
Figure 7 presents the calculated (normalized) electromagnetic pressure as a function of wavelength and structural parameters. The pressure normalized by the common value induced on a flat PEC surface can be written as [50]:(13)
Figure 7. The SFPresonance induced electromagnetic pressure in the microwave band as a function of wavelength and structure sizes: (a) a = 0.5 mm, t = 0.1 mm; (b) a = 0.5 mm, h = 10 mm; and (c) h = 10 mm, t = 0.1 mm. The lattice constant is fixed as d = 20 mm. 
Here, T_{m} is the amplitude of the mthorder diffraction mode in the spacer:(14)
Note that T_{m} is proportional to the FPlike resonance factor f_{p}. Figure 7a shows that the pressure for the longwavelength (zeroorder FP) resonance is positive whereas that for the shorterwavelength resonance is negative. In the former case, the energy “leaked” to the FP cavities is mainly the electric energy, which suppresses the attractive force. In the latter case, the antiparallel currents can also be induced in the metallic slit walls, due to the larger phase retardation in the slits; thus the energy “leaked” to the FP cavities becomes the magnetic energy, which suppresses the repulsive force. The effect becomes more significant with the increase of FP cavity length. For h = 10 mm and a = 0.5 mm, the normalized pressure attains ~1000 at the wavelength 64 mm and −500 at 30 mm.
The positive and negative pressure will decrease with the spacer thickness t (see Figure 7b), as the increase of t may reduce the electromagnetic energy density in the spacer. This may be useful for the dynamic control of the force. In addition, a reduction of slit width a can also boost the electromagnetic pressure, especially that for the zeroorder FP resonance (see Figure 7c). For example, when a is reduced from 1 mm to 0.1 mm, the normalized pressure changes from less than 500 (at 58 mm) to around 2000 (at 100 mm). With the decrease of slit width, the coupling between the adjacent metal stripes becomes stronger. Thus the currents at the bottoms of metal stripes may induce the horizontal electric dipoles crossing the narrow slits. Accordingly, the original (vertical) electric dipoles developed crossing the spacer are suppressed, which leads to a decrease of attractive force and enhancement of the total electromagnetic pressure.
5. SFP resonances at high frequencies
The SFP resonance effect may be extended to higher frequencies such as the THz or infrared regime by scaling the structural sizes with the wavelength [49, 91]. The difference is that the permittivity of metal becomes significantly smaller than that in the microwave band. Consequently, the skin depth scaled by the wavelength, , will be increased. This will modify the SFP resonance effect and the electromagnetic pressure. In the following, the THz and infrared region will be considered using two examples.
5.1 Phase shift at the cavity openings
The reflection phases at the cavity openings are checked again, considering the variation of slit width a and spacer thickness t. The result for varying the slit width a is similar to that obtained in the microwave band (see Figure 4a), but the case will be somewhat different concerning the variation of t. Figure 8a shows the wavelength dependence of reflection phase at the THz band for two values of t = 0.1 μm and 0.5 μm (d = 20 μm and a = 0.5 μm). We see that, besides the slope of phase curve, the position of zerophase point A varies with t (in the microwave band, the point A is almost independent of t). To see the effect more clearly, the inset of Figure 8a plots the position of A as a function of spacer thickness t (the circles). For the smaller t (t < 0.3 μm), a redshift with the decrease of t is seen. However, for the larger t (t > 0.3 μm), an opposite result can be obtained. As a comparison, Figure 8b also shows the wavelength dependence of reflection phase at the infrared band for two values of t = 10 nm and 30 nm (d = 2000 nm and a = 50 nm). The dependence of point A on the spacer thickness t is more significant, as shown by the circles in inset of Figure 8b. Thus, both the position of A and the slope of phase curve are governed by the spacer thickness t, which may in turn modulate the reflection phase φ_{3}.
Figure 8. (a) Phase shift in the THz band for two spacer thickness t = 0.1 and 0.5 μm (d = 20 μm, a = 0.5 μm); (b) Phase shift in the infrared region for t = 10 and 30 nm (d = 2000 nm, a = 50 nm). The inset of (a) and (b) shows, respectively, the zerophase point A as a function of t for the two frequency band, where the circles represent the theoretical results and the line a simple HFP model. The correspondence (in the infrared region) between the phase shift φ_{3} (the line) and (m = 0, 1, 2) are shown in (c), where d = 2000 nm, a = 50 nm, and t = 20 nm. 
To highlight the spacercavity effect in the high frequencies, φ_{3} and μ_{m} (here m = 0, 1, and 2) as a function of wavelength are also calculated. Here we only take the infrared regime as an example and the results are shown in Figure 8c. For the infrared regime and the smaller t, , where δ_{0} = c/ω_{p} is the skin depth. Thus, μ_{0} is suppressed due to the field penetration effect; μ_{1} and μ_{2} exhibit zero at the Rayleigh anomaly (v_{m} = 0) and maxima at the SPP resonance (2πm/d = ±k_{0}n_{t}, where a phase jump occurs; n_{t} is the effective refractive index of air spacer). Moreover, due to the field penetration effect, the zerophase point A, corresponding to , becomes dependent on t significantly.
We have also fitted the position of point A with a simple HFP model with λ ~ 2n_{t}(d − a). In the THz regime, ; and in the infrared regime, . The predictions of the simple HFP model are indicated by the lines in inset of Figures 8a and 8b, respectively. For the smaller t, the circles and the line agree well with each other. For larger t, a slight deviation can be seen which is due to the fact that the simple HFP model neglects the field decaying length (or reflection phase) at the periodic spacer openings. This decaying length is negligible for the smaller t but becomes noticeable for the larger t. Moreover, the agreement in the infrared regime seems better than that in the THz band, as the fields can be well confined due to the strong plasmonic effect.
5.2 Reflection spectrum
The reflection spectra for the THz band as a function of wavelength and structure parameters h, t, and a are plotted in Figures 9a–9c, respectively. The qualitative results are very similar to that for the microwave band (see Figure 5). However, due to a relatively larger penetration length, the absorption of THz wave becomes more significant, yielding a wider width of absorption spectrum. In Figure 9b, the absorption increases with decreasing spacer thickness t; but for smaller t (t < 0.2 μm), a degeneration of absorption can be seen instead (similar effect, not shown in Figure 5b, also exists in microwave band for very small t). According to equation (14), we have , where f_{p} was found to reduce with a decrease of t. Thus, the field amplitude and absorption will be suppressed in the case of small as well as large t. For completeness, the corresponding results for the infrared band are shown in Figures 9d–9f. The results, except for Figure 9e, are also similar to that obtained in the THz band. In Figure 9e, with the decrease of spacer thickness t, the SFP resonance wavelength blueshifts when t > 40 nm and redshifts when t < 40 nm. The redshift of SFP resonance at the smaller t is mainly caused by the shift of the zerophase point. This effect is more significant in the infrared region. The dependence of resonance on the spacer thickness may be used to modulate the SFP effect or measure the tiny motion with a high accuracy.
Figure 9. Reflection spectra as a function of wavelength and structure sizes for the THz (a–c) and infrared (d–f) regime. In the THz band, we fix d = 20 μm and set: (a) a = 1 μm, t = 0.5 μm; (b) a = 1 μm, h = 5 μm; and (c) h = 5 μm, t = 0.5 μm. In the infrared regime, we fix d = 2000 nm and set: (d) a = 50 nm, t = 20 nm; (e) a = 50 nm, h = 100 nm; and (f) h = 100 nm, t = 20 nm. 
5.3 Electromagnetic pressure
The electromagnetic pressure at the THz band as a function of wavelength and structural parameters h, t, and a are shown in Figures 10a–10c, respectively. We found that the pressure at the THz band can also be enhanced at the SFP resonance. At the longwavelength FP resonance, the pressure is positive; and at the shortwavelength resonance, the pressure is negative. To understand this effect, we should note that the kinetic energy of free electrons will play a role. At the longwavelength SFP resonance, the electric energy is “leaked” to the slit FP cavities, thus reducing the attractive force. However, the kinetic energy of free electrons suppresses the magnetic energy and the repulsive force slightly. Thus, a total positive but smaller pressure can be induced, which increases with the cavity length. At the shortwavelength resonance, on the other hand, the magnetic energy is “leaked” due to the slitFP cavity effect and field penetration effect (electron kinetic energy). This leads to a negative pressure.
Figure 10. Enhanced electromagnetic pressure as a function of wavelength and structure sizes for the THz (a–c) and infrared (d–f) regime. The structure parameters are the same as those set in Figure 9. 
In Figure 10b, a decrease of pressure (at longwavelength resonance) with decreasing t can be observed. For the smaller t, although the electric and magnetic field energy density can be enhanced, the ratio between the kinetic energy of free electrons and magnetic energy increases, which suppresses the total force (here is the skin depth). At the shortwavelength SFP resonance, nonetheless, the enhancement of the role of free electrons gives rise to an increase of negative pressure. However, when t is large enough (not shown here), the total pressure will be reduced eventually because of the decrease of field amplitude and the attractive/repulsive force.
The electromagnetic pressure at the infrared regime is quite different, as shown in Figure 10d–10f. The pressure associated with the SFP resonance is negative and can be enhanced by two orders of magnitude. In this regime, the kinetic energy of free electrons may exceed the magnetic energy and thus play a dominant role (suppresses the repulsive force and yields a negative pressure). In Figure 10d, with the increase of slit length and “leakage” of electric (or magnetic) energy, the pressure at the long (or short) wavelength SFP resonance will be reduced (or enhanced). In Figure 10e, a drastic decrease of optical pressure with the spacer thickness t can be found. This can be attributed mainly to the decrease of ratio between the electron kinetic energy and magnetic field energy. When t is large enough, a positive pressure at the SFP resonance can even be achieved. In addition, the decrease of slit width a results in a slight decrease of negative pressure at the longwavelength resonance (see Figure 10f). This is contrary to the result in the THz or microwave band. However, the physical origin is common: the coupling between the adjacent metal stripes increases, leading to an enhanced “leakage” of electric energy to the slit FP cavities.
6 Some additional discussions
6.1 Incidentangle dependence
The discussion presented above has been focused on the normal incidence of incident light. What will happen then when the incident angle is increased? Taking the infrared regime as an example, we plot in Figure 11 the phase shift and reflection spectrum for normal (a) and offnormal (b) incidence. For offnormal incidence, γ_{m} in previous equations should be revised as . For normal incidence, the oddorder HFP resonance modes exist, which cause in the spectrum two zerophase points and allow for the formation of two SFP resonances. For offnormal incidence, the positions of the oddorder HFP resonance modes remain unchanged, but the SPP resonance mode corresponding to the phase jump shifts to the long wavelength. This increases the slope of phase curve and leads to a slight blueshift of the SFP resonance around 8 μm. Moreover, when the incident angle is nonzero, the evenorder HFP resonance modes (which are inactive in the normal incidence) can also be resulted due to the breaking of system symmetry. Consequently, new reflection minima associated with the evenorder HFP modes are produced. The increase of incident angle can also reduce the electromagnetic pressure [50].
Figure 11. Incidentangle dependence of phase shift and reflection spectrum: (a) θ = 0^{o} and (b) θ = 30°. Here, φ_{1}, φ_{3}, and r_{0} are labeled by the open circles, solid circles, and the line, respectively. As an example, the structure parameters are set as d = 2000 nm, a = 50 nm, h = 100 nm, and t = 20 nm. 
6.2 Transition from SFP to FP resonance
Our discussion has also been restricted to the small spacer thickness t (t ≪ h), which is sufficient for the SFP resonances. However, one may ask, what will happen when the spacer thickness is large? To answer the question, Figure 12 plot the reflection spectrum with the spacer thickness varying from 20 nm to 3 μm (h = 100 nm). When t is small, two SFP resonances can be seen in the spectrum clearly. When t > 1 μm, new resonant reflection minima appear, which redshift with the wavelength scaling with the spacer thickness. The resonances are just the common (vertical) FP resonance in the spacer defined by the two metal reflectors. This point can be seen with the factor μ_{m}, where, for the very larger t, μ_{m} ≈ i/tan(k_{0}v_{m}t). Thus, the vertical FP resonances can be induced when k_{0}v_{0}t = mπ or 2t = mλ. Such resonance can also generate an enhanced (positive) electromagnetic pressure, but the magnitude and spectral width are significantly smaller than that of SFP resonance.
Figure 12. Reflection spectrum as a function of wavelength and spacer thickness, where a transition from the SFP resonance (in the slits) to the common FP resonance (in the spacer) is present. Here, d = 2000 nm, a = 50 nm, and h = 100 nm. 
6.3 Resonances in visible/nearinfrared band
For the metal thickness or slit length much larger than the skin depth, the SFP resonances usually locate in the midinfrared or even long wavelength regime. To conclude, we relax the constraint in the “SFP” and explore the properties of the structure in the visible or nearinfrared region. Figure 13a presented the reflection spectrum of the structure, where d = 600 nm, a = 40 nm, h = 155 nm, and t = 20 nm. Above the Rayleigh anomaly, two reflection minima locating around the wavelength 660 nm and 1064 nm are induced. Simultaneously, strongly enhanced negative pressure can also be found at the resonance. Figure 13b plots the phase shift at the slit opening as a function of wavelength. One can see that the reflection dips appear very close to the SPP resonance modes above or below the grating. Considering that the dip position is slightly modified by the slit length h, the resonance effect may be treated as a mixture of surface mode and slit cavity mode. However, the surface mode plays a dominant role in this case. The resonance effect in the visible or nearinfrared regime may be boosted by using the laser light.
Figure 13. (a) Reflection spectrum (the line) and enhanced pressure (the circles) in the visible and nearinfrared region; (b) Phase shift as a function of wavelength. Here d = 600 nm, a = 40 nm, h = 155 nm, and t = 20 nm. The resonance is mainly due to the grating surface mode rather than the SFP effect. 
7 Conclusions
In summary, a special kind of deep subwavelength FP resonance in a plasmonic structure has been discussed. The plasmonic system consists of a metallic slit grating backed with a metallic ground plane. It is crucial that the structure can present a unique phase shift at the bottom slit openings, which varies gently from the positive to negative values by crossing zero. This phase shift is correlated with both the evanescent effect of diffraction components and the cavity effect of the spacer sandwiched between the two metals. The negative phase shift can compensate the positive phase accumulation in the slit cavity, thus forming, typically, a zeroorder and deep subwavelength FP resonance, where the cavity length may be two orders of magnitude smaller than the resonance wavelength. By scaling the structural sizes with the wavelength, the effect can operate from the microwave band to the THz and infrared regime.
The SFP resonance can also induce an enhanced electromagnetic pressure. In the plasmonic system, the positive and negative forces depend on the magnetic and electric energy stored in the spacer, and a competition between them determines the sign and amplitude of the total pressure. In addition to the kinetic energy of free electrons, the “leakage” of electromagnetic energy from the spacer to the SFP resonators can modulate the total pressure. In the microwave band, the electron kinetic energy can be neglected. The “leakage” of electric (or magnetic) energy at the SFP resonances results in strong positive (or negative) pressure. The effect is similar in the THz band, but the electron energy begins to play a role. In the infrared regime, however, the kinetic energy of free electrons may surpass the magnetic energy and plays a crucial role. Consequently, the pressure becomes negative.
The phase shift of slit mode and the SFP resonance may be manipulated with the improved structures or designs. One can introduce a quasiperiodic rather than periodic sequence of slits to the structure, which may generate multiple SFP resonances. In addition, due to the dependence of phase shift on the slit width, the spectral width of SFP resonance can be widened by using the slits with varying slit width. This can be achieved by varying the slit width gradually along the slit array direction (the slit width is homogeneous in single slits) or employing the fanshaped slits with the inhomogeneous width in each slit (or with the inhomogeneous lattice period). On the other hand, the SFP resonance can be manipulated dynamically by varying the incident angle or the spacer thickness. For example, by connecting the planar metal film with a piezoelectric crystal, the airspacer thickness and SFP resonance can be controlled with an external electric field. It should also be noticed that one may fill the spacer with the specific materials and control the SFP resonance using the thermal, optical, electric, or magnetic methods.
In addition to the reflectionphase modulation, ultrasmall FP resonators can also be constructed by using a cavityindex design. For example, a combination of positive and negative index materials sandwiched between two mirrors may function as an ultrasmall cavity [92, 93]. On the other hand, the ultrasmall FP resonators can be realized with the highindex materials. For example, the indefinite medium with a unique hyperbolic dispersion can provide a larger kvector or high refractive index [94]. In this case, the cavity size can be much smaller than the vacuum wavelength (but not the wavelength in the cavity) [95, 96]. Currently, the reflectionphase and cavityindex designs for the compact FP resonators are employed independently. It may be interesting and simultaneously a challenge to consider both designs in an FP cavity. This may represent a further breakthrough towards the deep subwavelength FP resonators. A possible and simple scheme is to insert the highindex metameterials into a cavity, which is composed of the metamaterial reflectors. This would be feasible in the low frequencies such as the THz or microwave band.
8 Implications and influences
The results provide a method for designing ultrasmall FP resonators. Since no microstructured reflectors are required at the cavity opening, the cavity can be squeezed to small dimensions and operate at high frequencies. The SFP resonance may find applications in constructing ultrasensitive sensors. For example, by measuring the spectral shift of SFP resonance, one can detect the tiny motion or oscillation of ground plane with high accuracy. The SFP resonance is also accompanied by strong absorption and fieldenhancement, which may boost the lightmatter interactions, such as enhancing the molecule fluorescence and optical nonlinearity. In addition, the enhanced electromagnetic force can be employed to enhance the optomechanic coupling effect. With the focused and pulsed incident wave, mechanical motions of the devices or elastic waves can be excited efficiently. The SFP resonance effect may also be extended to the acoustic regime.
Acknowledgments
The authors would like to thank Dr. X.G. Yin, J.M. Zhao, and Prof. Y.J. Feng for their help in the microwave measurements. This work was supported by the Hong Kong Research Grants Council (Grant No. AOE/P02/12), the National Basic Research Program of China (Grant No. 2012CB921502), HKUST internal (Grant No. SRFI11SC07), and the National Natural Science Foundation of China (Grant No. 11174146).
References
 E. Yablonovitch, Phys. Rev. Lett. 58 (1987) 2059. [CrossRef] [PubMed] [Google Scholar]
 K.M. Ho, C.T. Chan, C.M. Soukoulis, Phys. Rev. Lett. 65 (1990) 3152. [CrossRef] [PubMed] [Google Scholar]
 F.R. Montero de Espinosa, E. Jimenez, M. Torres, Phys. Rev. Lett. 80 (1998) 1208. [CrossRef] [Google Scholar]
 R. Sainidou, N. Stefanou, A. Modinos, Phys. Rev. B 66 (2002) 212301. [CrossRef] [Google Scholar]
 Y.Q. Lu, et al., Science 284 (1999) 1822. [CrossRef] [Google Scholar]
 Y.Y. Zhu, et al., Phys. Rev. Lett. 90 (2003) 053903. [CrossRef] [Google Scholar]
 C.P. Huang, Y.Y. Zhu, AIP Advances 2 (2012) 042117. [CrossRef] [Google Scholar]
 Z. Liu, et al., Science 289 (2000) 1734. [CrossRef] [PubMed] [Google Scholar]
 G. Wang, et al., Phys. Rev. Lett. 93 (2004) 154302. [CrossRef] [PubMed] [Google Scholar]
 W.A. Murray, W.L. Barnes, Adv. Mat. 19 (2007) 3771. [CrossRef] [Google Scholar]
 E. Ozbay, Science 311 (2006) 189. [CrossRef] [PubMed] [Google Scholar]
 J.A. Schuller, et al., Nat. Mater. 9 (2010) 193. [CrossRef] [Google Scholar]
 D.K. Gramotnev, S.I. Bozhevolnyi, Nat. Photon. 4 (2010) 83. [CrossRef] [Google Scholar]
 V.M. Shalaev, Nat. Photon. 1 (2007) 41. [Google Scholar]
 X. Zhang, Z. Liu, Nat. Mater. 7 (2008) 435. [Google Scholar]
 S. Xi, et al., Phys. Rev. Lett. 103 (2009) 194801. [CrossRef] [Google Scholar]
 B. Min, et al., Nature 457 (2009) 455. [CrossRef] [PubMed] [Google Scholar]
 K.J. Vahala, Nature 424 (2003) 839. [Google Scholar]
 Y. Akahane, et al., Nature 425 (2003) 944. [CrossRef] [PubMed] [Google Scholar]
 J.B. Pendry, et al., IEEE Trans. Microwave Theor. Tech. 47 (1999) 2075. [Google Scholar]
 T.J. Yen, et al., Science 303 (2004) 1494. [CrossRef] [PubMed] [Google Scholar]
 S. Linden, et al., Science 306 (2004) 1351. [CrossRef] [PubMed] [Google Scholar]
 C.P. Huang, et al., Phys. Rev. Lett. 104 (2010) 016402. [CrossRef] [Google Scholar]
 N. Engheta, A. Salandrino, A. Alù, Phys. Rev. Lett. 95 (2005) 095504. [CrossRef] [Google Scholar]
 V. Delgado, et al., Metamaterials 3 (2009) 57. [CrossRef] [Google Scholar]
 C.P. Huang, et al., Opt. Express 17 (2009) 6407. [CrossRef] [Google Scholar]
 M. Staffaroni, et al., Photon. Nanostruct. Fund. Appl. 10 (2012) 166. [CrossRef] [Google Scholar]
 J. Zhou, et al., Phys. Rev. Lett. 95 (2005) 223902. [CrossRef] [PubMed] [Google Scholar]
 S. Zhang, et al., Phys. Rev. Lett. 101 (2008) 047401. [CrossRef] [PubMed] [Google Scholar]
 C. Fabry, A. Perot, Ann. Chim. Phys. 16 (1899) 115. [Google Scholar]
 J.M. Vaughan, The FabryPerot interferometer: history, theory, practice, and applications, Adam Hilger, Bristol, England, 1989. [Google Scholar]
 A. Kumar, H.D. Hristov, Microwave cavity antennas, Artech House, Norwood, MA, 1989. [Google Scholar]
 T.J. Kippenberg, K.J. Vahala, Opt. Express 15 (2007) 17172. [CrossRef] [PubMed] [Google Scholar]
 B. Hou, et al., Phys. Rev. B 76 (2007) 054303. [CrossRef] [Google Scholar]
 J.A. Porto, F.J. GarcıaVidal, J.B. Pendry, Phys. Rev. Lett. 83 (1999) 2845. [CrossRef] [Google Scholar]
 V.J. Sorger, et al., Nano Lett. 9 (2009) 3489. [CrossRef] [Google Scholar]
 A. Artar, A.A. Yanik, H. Altug, Appl. Phys. Lett. 95 (2009) 051105. [CrossRef] [Google Scholar]
 A. Chandran, et al., Phys. Rev. B 85 (2012) 085416. [CrossRef] [Google Scholar]
 V.I. Belotelov, et al., J. Opt. Soc. Am. B 29 (2012) 294. [CrossRef] [Google Scholar]
 R. Ameling, H. Giessen, Nano Lett. 10 (2010) 4394. [CrossRef] [Google Scholar]
 R. Ameling, H. Giessen, Laser & Photon. Rev. 7 (2013) 141. [CrossRef] [Google Scholar]
 G.V. Trentini, IRE Trans. Antennas Propag. 53 (1956) 666. [CrossRef] [Google Scholar]
 D. Sievenpiper, et al., IEEE Trans. Microwave Theor. Tech. 47 (1999) 2059. [Google Scholar]
 M.J. Lockyear, A.P. Hibbins, J.R. Sambles, Phys. Rev. Lett. 102 (2009) 073901. [CrossRef] [Google Scholar]
 L. Zhou, et al., Appl. Phys. Lett. 86 (2005) 101101. [CrossRef] [Google Scholar]
 H.Q. Li, et al., Appl. Phys. Lett. 89 (2006) 104101. [CrossRef] [Google Scholar]
 A. Ourir, A. de Lustrac, J. Lourtioz, Appl. Phys. Lett. 88 (2006) 084103. [CrossRef] [Google Scholar]
 F. Costa, A. Monorchio, Progress Electromagn. Res. 111 (2011) 467. [CrossRef] [Google Scholar]
 C.P. Huang, et al., Phys. Rev. B 85 (2012) 235410. [CrossRef] [Google Scholar]
 C.P. Huang, et al., Phys. Rev. B 86 (2012) 085446. [CrossRef] [Google Scholar]
 J. Li, G. Sun, C.T. Chan, Phys. Rev. B 73 (2006) 075117. [CrossRef] [Google Scholar]
 B. Auguie, W.L. Barnes, Phys. Rev. Lett. 101 (2008) 143902. [CrossRef] [Google Scholar]
 J.B. Pendry, L. MartinMoreno, F.J. GarciaVidal, Science 305 (2004) 847. [CrossRef] [PubMed] [Google Scholar]
 F.J. GarciaVidal, L. MartinMoreno, J.B. Pendry, J. Opt. A: Pure Appl. Opt. 7 (2005) S97. [CrossRef] [Google Scholar]
 J.M. Hao, L. Zhou, C.T. Chan, Appl. Phys. A 87 (2007) 281. [CrossRef] [Google Scholar]
 H. Raether, Surface plasmons on smooth and rough surfaces and on gratings, SpringerVerlag, New York, 1988. [Google Scholar]
 C. Girard, A. Dereux, Rep. Prog. Phys. 59 (1996) 657. [CrossRef] [Google Scholar]
 L. Salomon, F. de Fornel, P.M. Adam, J. Opt. Soc. Am. A 16 (1999) 2695. [CrossRef] [Google Scholar]
 A.V. Zayats, I.I. Smolyaninov, A.A. Maradudin, Phys. Rep. 408 (2005) 131. [CrossRef] [Google Scholar]
 C.P. Huang, et al., Appl. Phys. Lett. 93 (2008) 081917. [CrossRef] [Google Scholar]
 T. LópezRios, et al., Phys. Rev. Lett. 81 (1998) 665. [CrossRef] [Google Scholar]
 Y. Takakura, Phys. Rev. Lett. 86 (2001) 5601. [CrossRef] [PubMed] [Google Scholar]
 A.P. Hibbins, M.J. Lockyear, J.R. Sambles, J. Appl. Phys. 99 (2006) 124903. [CrossRef] [Google Scholar]
 A.P. Hibbins, et al., Phys. Rev. Lett. 92 (2004) 143904. [CrossRef] [Google Scholar]
 J.R. Brown, et al., J. Appl. Phys. 104 (2008) 043105. [CrossRef] [Google Scholar]
 M.J. Lockyear, et al., Appl. Phys. Lett. 94 (2009) 041913. [CrossRef] [Google Scholar]
 Y. Todorov, et al., Phys. Rev. Lett. 102 (2009) 186402. [CrossRef] [PubMed] [Google Scholar]
 Y. Todorov, et al., Opt. Express 18 (2010) 13886. [CrossRef] [Google Scholar]
 P. Jouy, et al., Appl. Phys. Lett. 98 (2011) 021105. [CrossRef] [Google Scholar]
 Z. Wei, et al., Opt. Express 18 (2010) 12119. [CrossRef] [Google Scholar]
 Y.Q. Ye, Y. Jin, S.L. He, J. Opt. Soc. Am. B 27 (2010) 498. [CrossRef] [Google Scholar]
 J. Hao, L. Zhou, M. Qiu, Phys. Rev. B 83 (2011) 165107. [CrossRef] [Google Scholar]
 H.T. Chen, Opt. Express 20 (2012) 7165. [CrossRef] [PubMed] [Google Scholar]
 Y.G. Ma, C.K. Ong, X.S. Rao, J. Appl. Phys. 103 (2008) 123510. [CrossRef] [Google Scholar]
 N. Mattiucci, et al., Appl. Phys. Lett. 101 (2012) 141115. [CrossRef] [Google Scholar]
 J. Chen, et al., Nat. Photon. 5 (2011) 531. [CrossRef] [Google Scholar]
 A. Ashkin, Optical trapping and manipulation of neutral particles using lasers, World Scientific, Singapore, 2006. [Google Scholar]
 A. Schliesser, et al., Phys. Rev. Lett. 97 (2006) 243905. [CrossRef] [PubMed] [Google Scholar]
 T.J. Kippenberg, K.J. Vahala, Science 321 (2008) 1172. [CrossRef] [PubMed] [Google Scholar]
 M. Eichenfield, et al., Nature 459 (2009) 550. [CrossRef] [PubMed] [Google Scholar]
 D. Van Thourhout, J. Roels, Nat. Photon. 4 (2010) 211. [CrossRef] [Google Scholar]
 H. Xu, M. Kall, Phys. Rev. Lett. 89 (2002) 246802. [CrossRef] [PubMed] [Google Scholar]
 L. Huang, O.J.F. Martin, Opt. Lett. 33 (2008) 3001. [CrossRef] [Google Scholar]
 J. Kohoutek, et al., Nano Lett. 11 (2011) 3378. [CrossRef] [Google Scholar]
 R. Sainidou, F.J. Garcıa de Abajo, Phys. Rev. Lett. 101 (2008) 136802. [CrossRef] [Google Scholar]
 X. Yang, et al., Nano Lett. 11 (2011) 321. [CrossRef] [Google Scholar]
 D. Woolf, M. Loncar, F. Capasso, Opt. Express 17 (2009) 19996. [CrossRef] [Google Scholar]
 J. Zhang, K.F. MacDonald, N.I. Zheludev, Phys. Rev. B 85 (2012) 205123. [CrossRef] [Google Scholar]
 H. Liu, et al., Phys. Rev. Lett. 106 (2011) 087401. [CrossRef] [Google Scholar]
 S.B. Wang, et al., Phys. Rev. B 84 (2011) 075114. [CrossRef] [Google Scholar]
 Y. Zhang, M. Han, C.P. Huang, Europhys. Lett. 102 (2013) 34001. [CrossRef] [Google Scholar]
 N. Engheta, IEEE Antennas Wirel. Propag. Lett. 1 (2002) 10. [Google Scholar]
 L. Shen, S.L. He, S. Xiao, Phys. Rev. B 69 (2004) 115111. [CrossRef] [Google Scholar]
 D.R. Smith, D. Schurig, Phys. Rev. Lett. 90 (2003) 077405. [CrossRef] [PubMed] [Google Scholar]
 X. Yang, et al., Nat. Photon. 6 (2012) 450. [CrossRef] [Google Scholar]
 J. Yao, et al., PNAS 108 (2011) 11327. [CrossRef] [Google Scholar]
Cite this article as: Huang CP & Chan CT: Deep subwavelength FabryPerot resonances. EPJ Appl. Metamat. 2014, 1, 2.
All Figures
Figure 1. Reflection of light at various interfaces: (a) A semiinfinite PEC metal milled with subwavelength square holes behaves like an effective medium (λ ≫ d) with a resonant impedance (the hole width is a, the period is d, and the permittivity of hole filling medium is set as ε_{h} = 25); (b) A semiinfinite PEC metal milled with periodic slits (the slit width is a and the period is d) can induce the diffracted evanescent waves at the interface (λ ~ d); (c) Light reflection by a dielectric slab (with a thickness h) sandwiched between two semiinfinite media. The reflection efficiency and phase shift for the case (a)–(c) are shown in (d)–(f) respectively. We set in (d) and (e) a = d/2 and in (f), n_{1} = 1, n_{2} = 2, and n_{3} = 1.5. 

In the text 
Figure 2. Schematic view and physical parameters of the plasmonic system under study. The slit grating (with a narrow slit width) and ground plane are separated by a thin dielectric spacer. The vertical narrow slits support a deep subwavelength FPlike resonance. Such resonance can be excited by a TMpolarized incident light. 

In the text 
Figure 3. The SFP resonance in the microwave band for the structure with d = 20 mm, a = 0.3 mm, h = 5 mm, and t = 0.1 mm. (a) Phase shift of slit mode at the slit openings (the arrow indicates the SFP resonance); (b) Zeroorder reflection (the line) and the absolute value of the FPlike resonance factor (the circles). The resonance around 20 mm is dominated by the grating surface modes. 

In the text 
Figure 4. Dependence of phase shift on the wavelength for different slit width a and spacer thickness t: (a) t = 0.1 mm, and a = 1.0 or 0.3 mm; (b) a = 1.0 mm, and t = 0.1 or 0.3 mm. Here d = 20 mm. The correspondence between the phase shift φ_{3} (the line) and (m = 0, 1, 2) are shown in (c), where d = 20 mm, a = 0.5 mm, and t = 0.2 mm. The inset of (c) shows the zerophase point A (the circles) as a function of lattice period d (where we fix a = 1.0 mm and t = 0.2 mm), which track the position of HFP mode in the spacer closely (the line). 

In the text 
Figure 5. Reflection spectrum as a function of wavelength and structure sizes: (a) a = 0.5 mm, t = 0.1 mm; (b) a = 0.5 mm, h = 5 mm; and (c) h = 5 mm, t = 0.1 mm. The lattice constant is fixed as d = 20 mm. 

In the text 
Figure 6. Experimental results for the microwave band: (a) Photographs of the fabricated samples. The lattice constant of the two samples is d = 20 mm, the slit width is a = 1 mm, and the grating thickness is h = 5 mm and 10 mm, respectively. The measured (the circles) and calculated (the line) reflection spectra are presented in (b) for h = 5 mm, t = 0.15 mm and (d) for h = 10 mm, t = 0.3 mm, respectively. The measured (the circles) and calculated (the lines) resonance wavelength of the two samples as a function of spacer thickness t are shown in (c), h = 5 mm, and (e), h = 10 mm, respectively. 

In the text 
Figure 7. The SFPresonance induced electromagnetic pressure in the microwave band as a function of wavelength and structure sizes: (a) a = 0.5 mm, t = 0.1 mm; (b) a = 0.5 mm, h = 10 mm; and (c) h = 10 mm, t = 0.1 mm. The lattice constant is fixed as d = 20 mm. 

In the text 
Figure 8. (a) Phase shift in the THz band for two spacer thickness t = 0.1 and 0.5 μm (d = 20 μm, a = 0.5 μm); (b) Phase shift in the infrared region for t = 10 and 30 nm (d = 2000 nm, a = 50 nm). The inset of (a) and (b) shows, respectively, the zerophase point A as a function of t for the two frequency band, where the circles represent the theoretical results and the line a simple HFP model. The correspondence (in the infrared region) between the phase shift φ_{3} (the line) and (m = 0, 1, 2) are shown in (c), where d = 2000 nm, a = 50 nm, and t = 20 nm. 

In the text 
Figure 9. Reflection spectra as a function of wavelength and structure sizes for the THz (a–c) and infrared (d–f) regime. In the THz band, we fix d = 20 μm and set: (a) a = 1 μm, t = 0.5 μm; (b) a = 1 μm, h = 5 μm; and (c) h = 5 μm, t = 0.5 μm. In the infrared regime, we fix d = 2000 nm and set: (d) a = 50 nm, t = 20 nm; (e) a = 50 nm, h = 100 nm; and (f) h = 100 nm, t = 20 nm. 

In the text 
Figure 10. Enhanced electromagnetic pressure as a function of wavelength and structure sizes for the THz (a–c) and infrared (d–f) regime. The structure parameters are the same as those set in Figure 9. 

In the text 
Figure 11. Incidentangle dependence of phase shift and reflection spectrum: (a) θ = 0^{o} and (b) θ = 30°. Here, φ_{1}, φ_{3}, and r_{0} are labeled by the open circles, solid circles, and the line, respectively. As an example, the structure parameters are set as d = 2000 nm, a = 50 nm, h = 100 nm, and t = 20 nm. 

In the text 
Figure 12. Reflection spectrum as a function of wavelength and spacer thickness, where a transition from the SFP resonance (in the slits) to the common FP resonance (in the spacer) is present. Here, d = 2000 nm, a = 50 nm, and h = 100 nm. 

In the text 
Figure 13. (a) Reflection spectrum (the line) and enhanced pressure (the circles) in the visible and nearinfrared region; (b) Phase shift as a function of wavelength. Here d = 600 nm, a = 40 nm, h = 155 nm, and t = 20 nm. The resonance is mainly due to the grating surface mode rather than the SFP effect. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.