## Abstract

Photonic band properties are presented for a two-dimensional rectangular-groove grating of metal into air. The properties of the surface modes are shown and discussed with a perfect electric conductor, and compared to those of surface plasmons with real metal. The same structure is also studied with real metal in the near infrared. The results are obtained with a 3-D finite element numerical code.

©2006 Optical Society of America

## 1. Introduction

Nanoscale patterning of metal surfaces is a very powerful technique, opening new perspectives in the field of photonics, such as ultra-compact photonic integrated circuits, new kinds of optical sources and sub-diffraction limited imaging [1]. In particular, a complete photonic band gap in the dispersion of surface plasmons (SPs) [2] is obtained through periodic patterning of metallic surfaces. These structures are the surface equivalent of photonic crystals. Instead of being confined inside a dielectric, the optical mode is in vacuum, although tightly bounded to the surface. These properties imply that the local density of photonic states is modified and so the emission from nanoscale optical sources (the Purcell effect). This effect was exploited to improve the emission properties of nanocrystals [3]. By combining double periodicity with a “defect” a surface concentration effect was predicted, which can be exploited in detection or emission enhancement from a point source[4].

A very intuitive picture for explaining the behavior of these periodically patterned surfaces could be used: a bound mode, the surface plasmon (SP), exists due to the plasmon resonance into the metal for a flat surface ; therefore, the periodic modification of the properties of the surface (e.g. by etching the metal layer, or modifying its profile) induces the Bragg scattering pretty much as holes do in dielectric photonic crystals. This perturbative approach is not adequate for understanding the more general properties of patterned surfaces. First of all, periodic patterning not only generates Bragg scattering, but also modifies the degree of confinement of the mode at the surface. Indeed, the mode of the low frequency band is more strongly confined (stronger extinction) than the SPs mode on a flat surface, while the high frequency band mode is less [12, 13]. This effect becomes stronger as the depth of the patterning increases ; moreover the SP effect is reduced when moving from the visible spectral range to the infrared. Thus, it is expected that the high-frequency band disappears below a cut-off frequency, which depends on the materials (metal and dielectric) and on the depth of the patterning. In any case, surface wave exists in any spectral range, even in the limit of perfect electric conductors (PEC), at least in 1D periodic structures defined on the x-y plane which are translation-invariant along z.

The case of 2D periodic surfaces in the 3D space is more complex. Very recently, Pendry made the connection between surface plasmons and surface waves on PEC surfaces with holes[6]. Indeed, the cut-off frequency in metallic holes acts as a plasmon resonance and evanescent penetration of the field in the holes mimics the skin effect. The experimental demonstration of such effect followed [10].

The recent literature on patterned PEC surfaces focuses on a specific geometry, that is plane PEC surfaces with sub-wavelength holes [7, 8, 9]. In this Paper we consider a complementary case, that is a plane PEC surface with bumps or other kind of protuberances forming a regular array. This is the case investigated experimentally [2] and theoretically [15] in the visible spectral range. We focus on the long wavelength limit and we will show that this structure supports surface waves and, depending on the depth of the patterning, more than one band under the air light line. We discuss the nature and the properties of such modes in the case of a rectangular lattice and compare the case of PEC structures with real metal in the near infrared. We also show how the properties of this band structure depend on geometric parameters, which can be controlled by surface designers.

## 2. Discussion

In our structure, the metal/dielectric interface consists in a two-dimensional rectangular-groove grating of metal into air, as shown in Fig. 1. The structure is symmetric, with the same period, Λ, in both orthogonal directions x and y. The bumps have a square cross-section α×α, with height *h*.

We compute the photonic band structure by solving a suitable eigenvalue problem according to the Floquet-Bloch theorem. Thus, only one unit cell, with size Λ×Λ, is simulated, limited by two lateral planes orthogonal to x direction and two others orthogonal to y direction. Boundary conditions imposed on two parallel lateral surfaces respect the field’s property to be the same except for the phase, with a phase difference function of the wave vector. The eigenvalue problem is solved by HFSS (Ansoft), a 3-D finite-element commercial code. The height of the computational domain in the air is chosen to be large enough (equal to Λ) to ensure the existence of the surface modes, and to avoid artefact modes. We define Perfectly Matched Layers (PML) conditions in the z direction outside of the air domain. When non-perfect conductor is used, the metal is included in the computational domain and is meshed, and Perfect E boundary conditions are defined outside of the metal, to limit its thickness to $\frac{\Lambda}{15}$. The precision for eigenmode solutions in HFSS is given by the difference in the resonant frequencies from one adaptative solution to the next in the convergence process. The mesh size we used allowed a precision of 0.1% for PEC and 1% with real metal. To ensure that HFSS captures the physics associated with surface plasmons, we have computed the frequency of the SP constant on a metallic plane surface for metal with ε* _{m}*=-50. The theory predicts a propagation constant for a surface plasmon ${k}_{\mathit{sp}}=\frac{\omega}{c}\sqrt{{\epsilon}_{m}\u2044\left({\epsilon}_{m}+1\right)}$, corresponding to

*f*.

*π*/(

*k*.

_{sp}*c*)=0.495. The value calculated with HFSS is 0.497, which is reasonably close.

The calculated band structures are presented for two values of *α*, *α*=Λ/2 and α=Λ/√2, and for two values of the depth *h*, *h*=0.2×Λ (figure 2 a) and h=0.4×Λ (figure 2 b). The air light line, (thin dark line), defines the domain between two regions : the modes above the light line (shaded region) are leaky as they can couple with radiation. Thus, only modes with frequency below the light line are really confined at the surface and are of real interest for us.

The existence of surface modes with a corrugated PEC surface is clearly established by the diagrams, since for each structure there is at least one mode below the light line. Band flattening is remarkable between *h*=0.2×Λ and *h*=0.4×Λ.

Fig. 3 represents the Electric Field distribution of the surface Bloch modes at the high symmetry points. We represent the *E _{z}* component as it is the dominant polarisation. We also indicate the polarization of the field with arrows. Both figures correspond to a case with PEC, and

*h*=0.2×Λ and α=Λ/√2. The field is confined in the z direction at the interface PEC/air.

Let us first consider the low frequency mode at the *X _{X}* and

*X*points, $\overrightarrow{k}=(\frac{\pi}{\Lambda},0)$ and $\overrightarrow{k}=(0,\frac{\pi}{\Lambda})$ respectively, which are doubly degenerated in frequency (see the Brillouin’s zone sketched in Fig. 2). The

_{Y}*E*field at the

_{z}*X*point is periodic with period Λ along the

_{X}*x*direction and it is almost invariant along

*y*with anti-nodes centered on the bumps (Fig. 3(a)). Note the changes of sign every Λ in the

*x*direction. The mode at the

*M*point, $\overrightarrow{k}=(\frac{\pi}{\Lambda},\frac{\pi}{\Lambda})$, is represented in figure 3 b, the field changes sign every Λ in both x and y directions and the anti-nodes are centered on the bumps.

A high frequency band appears below the light line near the *M* point (see. Fig. 2(a)). The upper band mode at the *M* point is doubly-degenerated in frequency: the field distribution is consistent with the symmetries of the *M* point, but differs from that of the lower frequency band in the position of the anti-nodes. More specifically, the anti-nodes are centered *between two bumps*, which leaves two possibilities (Fig. 3(c),(d))

We have shown that a PEC surface with an arrays of bumps supports surface waves and exhibit a complete photonic band gap. In this respect, a patterned PEC surface keeps the properties of SPs on corrugated surfaces [2, 15], although the underlying physical mechanism is not related to the plasmon resonance. First of all, we notice that patterning of a PEC surface results in general into a penetration of the electric field with respect to an “average” flat surface. However, there is an important difference between an array of holes, which is considered in [6, 7, 8, 9] and an array of bumps. The first case corresponds to simply connected domains into the metal, forming waveguides with a cut-off frequency, which is the key to understand the properties of the surface waves on PEC surfaces. The case of array of bumps is different, as the grooves between the bumps does not form simply connected domains, therefore no cut-off exists.

Several papers on the dispersion properties of sinusoidal, mono-periodical gratings of real metal [12, 13, 2, 14] explain the onset of a photonic bandgap when the effective wavelength of the SP mode is twice the period of the grating, that is *k _{x}*=

*π*/Λ at the X point and

*k*=

_{x}*k*=

_{y}*π*/Λ at the M point. The lower and the upper bands correspond to different field distributions. The field of the low frequency mode is concentrated at the peaks of the grating, while the field of the high frequency mode is concentrated at the troughs; the field lines are more distorted in the latter case and consistently the greater energy stored in the fields by this mode corresponds to a higher frequency [12]. From this point of view, the modes of 2D PEC gratings have the same properties as in 1D gratings, no matter if PEC or not. We notice that, in our 2D case, even the high frequency band exists, with exactly the same kind of field distribution as it would be expected with metal with finite negative permittivity.

To support that, we compare the near infrared with the PEC case ; the real permittivity of the metal is -50 (complex permittivity, thus losses are not taken into account). Only points at the X and M points are calculated, and they are shown with solid dots in Fig. 2. The low and high frequency bands are lower with real metal in the near infrared. The Electric Field distribution of these modes is presented (Fig. 4) for α=Λ/√2 and *h*=0.2Λ, and is very similar to the case with PEC (Fig. 3).

It is also interesting to notice that a new mode appears at the X point for a deeper grating (*h*=0.4Λ, Fig. 2(b)). The Electric Field distribution of this mode is shown (Fig. 5) for *E _{x}*,

*E*and

_{y}*E*, since, oppositely to the previous cases,

_{z}*E*is not prevailing on the other fields components ; the antinodes are centered between two bumps. The periodicity Λ/2 for

_{z}*E*shows that this mode is of higher order than the mode of Fig. 4(a). We think that it could be the mode situated at (π/Λ, 2π/Λ) in the reciprocital lattice, corresponding to the

_{z}*X*point in the Brillouin zone (Fig. 5(d)). This mode being further from the G point, it has a higher frequency than the mode of Fig. 4(a), and thus appears for deeper gratings.

_{X}Finally, we show in Fig. 6 the positions of the modes at the X and M points for several values of $\frac{h}{\Lambda}$, for a PEC surface and for α=Λ/√2. The normalized frequency $\frac{f.\Lambda}{c}$ of each set of modes is decreasing when $\frac{h}{\Lambda}$ increases, because the modes become less well bound to the surface, and their energy and thus their frequency is decreased [12]. For nearly shallow gratings $\left(\frac{h}{\Lambda}=0.1\right)$, only the low frequency modes exist at the X and M points, and they are very close to the light line. For deeper grooves, the normalized frequency of these modes decreases and a new mode appears at the M point ($\frac{h}{\Lambda}=0.2$ and 0.4) then at the X point ($\frac{h}{\Lambda}=0.6$ and 0.75). This X mode has the same field distribution as for real metal (Fig. 5). The existence of this mode labeled *X*
_{1} is clearly due to the geometry of the grating, while it only appears for deep grooves. It is interesting to mention here that the causal role for the existence of this mode with real metal (Fig. 5) is played by the geometry of the structure and not by surface plasmons, although SPs may enhance the phenomena.

## 3. Conclusions

In conclusion, we have investigated the photonic bands of surface modes supported by a PEC structure consisting in an array of bumps. This structure is complementary with the array of holes considered very recently in the litterature. We find that very similar properties, typically associated with surface plasmons, still hold in the long wavelength limit. We conclude that surface gratings, allowing the penetration of the field in the structure, give birth to plasmon-like waves. We think that a better understanding of the properties of surface waves on PEC will be very important for infrared and THz photonics.

## References and links

**1. **E. Ozbay“Plasmonics: Merging Photonics and Electronics at Nanoscale Dimensions,” Science **311**, pages 189–193 (2006). [CrossRef] [PubMed]

**2. **S. C. Kitson, W. L. Barnes, and J. R. Sambles, “Full photonic band gap for surface modes in the visible,” Phys. Rev. Lett. **77**, 2670–2673 (1996). [CrossRef] [PubMed]

**3. **J. Zhang, Y.-H. Ye, X. Wang, P. Rochon, and M. Xiao, “Coupling between semiconductor quantum dots and two-dimensional surface plasmons,” Phys. Rev. B **72**, 201306 (2005). [CrossRef]

**4. **M. Carras and A. De Rossi, “Field concentration by exciting surface defect modes,” Opt. Lett. **31**, pages 47–49 (2006). [CrossRef] [PubMed]

**5. **W.L. Barnes, A. Dereux, and T.W. Ebbesen, “Surface plasmon subwavelength optics,” Nature **424**, 824–830 (2003). [CrossRef] [PubMed]

**6. **J. B. Pendry, L. Martin-Moreno, and F.J. Garcia-Vidal, “Mimicking Surface Plasmons with Structured Surfaces,” Science **305**, 847–848 (2004) [CrossRef] [PubMed]

**7. **F.J. Garcia-Vidal, L. Martin-Moreno, and J. B. Pendry, “Surfaces with holes in them : new plasmonic metamaterials,” J. Opt. A:Pure Appl. Opt. **7**, S97–S101 (2004) [CrossRef]

**8. **F.J. Garcia de Abajo and J.J. Saenz, “Electromagnetic surface modes in structured perfect-conductor surfaces,” Phys. Rev. Lett. **95**, 233901-1-4 (2005). [CrossRef]

**9. **Min Qiu, “Photonic band structures for surface waves on structured metal surfaces,” Opt. Express **13**, 7583–7588 (2005). [CrossRef] [PubMed]

**10. **A.P. Hibbins, B.R. Evans, and J.R. Sambles, “Experimental Verification of Designer Surface Plasmons,” Science **308**, 670–672 (2005). [CrossRef] [PubMed]

**11. **W. Barnes and R. Sambles, “Only Skin Deep,” Science **305**, 785–786 (2004). [CrossRef] [PubMed]

**12. **W. L. Barnes, T. W. Preist, S. C. Kitson, J. R. Sambles, N.P.K. Cotter, and D.J. Nash “Photonic gaps in the dispersion of surface plasmons on gratings,” Phys. Rev. B **51**, 11 164–11 168 (1995). [CrossRef]

**13. **W. L. Barnes, T. W. Preist, S. C. Kitson, and J. R. Sambles, “Physical origin of photonic energy gaps in the propagation of surface plasmons on gratings,” Phys. Rev. B **54**, 6227–6244 (1996). [CrossRef]

**14. **A. Giannattasio and W. L. Barnes, “Direct observation of surface plasmon-polariton dispersion,” Opt. Express **13**, 428–434 (2005). [CrossRef] [PubMed]

**15. **M. Kretschmann “Phase diagrams of surface plasmon polaritonic crystals” Phys. Rev. B **68**, 125419 (2003). [CrossRef]