4 December 2012 Negative refraction via radiative surface phonon polaritons in a silicon carbide-based multilayer structure
Author Affiliations +
Optical Engineering, 51(12), 120502 (2012). doi:10.1117/1.OE.51.12.120502
Abstract
In this letter, we predict that radiative surface phonon polariton (RSPhP) modes, in a multilayer structure consisting of thin-film silicon carbide (SiC) bounded by silicon (Si) and diamond (Di), can be used to achieve negative refraction of mid-infrared light at approximately 11 μm. Dispersion relations, calculated for the Si/SiC/Di structure, show that the RSPhP mode exhibits negative dispersion and couples with incident light. Poynting vector calculations show how the energy flux may be refracted, negatively, in the SiC layer.
Adewuyi and Hammonds: Negative refraction via radiative surface phonon polaritons in a silicon carbide–based multilayer structure

Currently, there is great interest in materials that refract light negatively. Negative refraction can be characterized by a light beam whose wavefront is refracted normally (in a positive direction), while its energy flux is refracted negatively. Materials that negatively refract light can enable applications such as near-field superlensing and optical cloaking.1 Approaches to developing negative refraction materials include the design of photonic crystals with geometries that result in desired dispersion properties. The photonic crystal approach has produced some notable successes in the development of microwave to infrared negative refraction materials.2,3 Another approach is to exploit material resonance modes to achieve negative refraction.4,5 Recent examples include the use of infrared plasma resonance in free-carrier doped InGaAs films to suport negative refraction.6 Magnetic resonance, in gold, has also been applied to the development of mid-IR negative refraction materials.7 In this letter, we propose that coupled photon-phonon resonance modes, called radiative surface phonon polaritons (RSPhP), can be exploited in the development of mid-IR negative refraction materials. To demonstrate this idea, we present theoretical results for a prototype structure consisting of a thin silicon carbide (SiC) film, bounded by very thick layers of silicon (Si) and diamond (Di). Approaches to achieving mid-IR negative refraction materials are necessary for the development of next-generation thermal management devices such as thermal rectifiers and specialized thermal emitters.8,9

Physically, surface phonon polaritons are atomic lattice-waves coupled to transverse magnetic (p-polarized) electromagnetic fields that may be supported at the interface between a polar material bounded by a dielectric. For example, surface phonon polaritons can be excited at the interface between air (dielectric constant, ϵ01) and a semi-infinite layer of SiC [a polar material with a frequency varying dielectric function ϵ1(ω)]. The SiC/Air system is illustrated in Fig. 1(a), and the well-known dispersion relation for surface phonon polariton modes at an interface formed between a semi-infinite polar material and a dielectric is given by the expression

(1)

kx=ωcϵ0ϵ1(ω)ϵ0+ϵ1(ω),
where the dependency on frequency ω is made explicit and where c is the speed of light. The surface phonon polariton dispersion relation for the SiC/Air case is plotted as the dotted line in Fig. 2. The vertical axis of Fig. 2 is the frequency at which the surface mode exists and the horizontal axis is (the real part of) kx, which is the wavevector for the surface polariton that is propagating parallel to the SiC interface. The “light line” is the dashed diagonal line through the plot and gives the dispersion relation for light propagation in free space, or k=ω/c. The upper and lower limits of the frequency axis, in Fig. 2, approximately correspond to the spectral locations of the SiC transverse optical frequency ωTO and its longitudinal optical frequency ωLO, which gives the SiC/Air polariton dispersion cutoff. Between these frequency limits, the dielectric function of SiC is less than zero, i.e. ϵ1(ω)<0, and a negative dielectric indicates that the material may support electromagnetic resonance modes. More specifically, the SiC dielectric function can be described by the Lorentz model ϵSiC(ω)=ϵ(ωLO2ω2iγω)/(ωTO2ω2iγω) where ϵ=6.7, ωLO=18.3×1013s1, ωTO=14.9×1013s1, and the damping constant is given as γ=0.89×1012s1. In this letter, the loss term γ is included in all calculations.

Fig. 1

(a) Schematic of semi-infinite SiC/Air system. (b) Schematic of Si/SiC/Di structure studied in this work, where d is the SiC thickness, and θ is the angle of the incident light.

OE_51_12_120502_f001.png

Fig. 2

Frequency versus the real part of kx. Thin Dashed Line: Free space light line. Dotted Line: SPhP dispersion relation for semi-infinite SiC bounded by air [Fig. 1(a)]. Solid Thick Lines: SPhP dispersion relation for p-polarized light in multilayer structure shown in Fig. 1(b). Thick Dashed Line: Dispersion relation for s-polarized light in the multilayer structure. Thin Solid Line: Dispersion relation for semi-infinite SiC bounded by diamond.

OE_51_12_120502_f002.png

The wavevector representing the component of light propagating perpendicular to i’th interface is defined by kzi=ϵik2kx2. Therefore, surface modes with kx>ϵik cannot couple to propagating light, since kzi would be purely imaginary. Figure 2 shows that the surface phonon polariton mode for the SiC/Air structure is to the right of the light line and, thus, kx is too large for a beam of incident light to interact with the surface phonon polariton.

So, while it is possible for semi-infinite SiC to support mid-IR surface phonon polaritons, these modes are forbidden from direct coupling to propagating light. However, the situation completely changes for the case of a thin, 400 nm layer of SiC, sandwiched between semi-infinite thick layers of Si and diamond, as shown in Fig. 1(b). These materials were chosen because they are transparent in the infrared region and have high refractive indexes with nSi=3.5 and nDi=2.4. The dispersion relation, for surface phonon polaritons in the Si/SiC/Di structure, is determined by numerically finding values of kx(ω) that result in zeros for the denominator of the three-layer, p-polarized, reflectance equation,10

(2)

R210p=|r21p+r10pe2iα11+r21pr10pe2iα1|2,
where

(3)

rijp=(kziϵikzjϵj)/(kziϵi+kzjϵj),
α1=kz1d, and the subscripts refer to the layers i,j=0,1, and 2, where diamond is layer 0, SiC is layer 1, and Si is layer 3. The numerically determined surface phonon polariton dispersion relation for the Si/SiC/Di structure is given by the solid thick lines of Fig. 2, and shows three branches. Two branches are completely to the right of the light line and represent modes unaccessible by propagating light. As shown in the figure, the two modes are the antisymmetric and symmetric modes which are well-known to represent the surface polaritons at each interface the thin film.10 The plot also shows the dispersion relation for semi-infinite SiC bounded by diamond. While the horizontal axis range is limited in Fig. 2 for brevity, previous theoretical results suggest that for very large kx, degeneracy is lifted and the symmetric and anti-symmetric modes become equivalent to the semi-infinite SiC/Diamond case. The SPhP modes are known to only exist for the p-polarized case. For comparison, we also plot the dispersion relation for the s-polarized case in Fig. 2 to show that the effects discussed in this letter only occur for p-polarized light.

The third branch of the Si/SiC/Di dispersion relation extends to the left of the light line, crossing it at approximately ω=1.65×1014rads1, or 11.4 μm. Since this mode crosses the light line and is, thus, accessible by propagating light with a wavelength of about 11.4 μm, we refer to it as a radiative surface phonon polariton (RSPhP). As will be shown below, calculation of the magnetic field shows that the mode is symmetric since the transverse component of the electric field (Ez) is not zero inside the film [c.f. Fig. 3 and Eq. (4)]. Symmetric (and anti-symmetric) radiative surface polaritons have been predicted, before, for the case of multilayer structures. However, our calculations show that the radiative branch of the dispersion relation presented here has a negative slope. Since the group velocity is given by vg=ω/k, the negative slope implies a negative group velocity. Also, since the group velocity characterizes the direction of energy flow, which in this case is roughly opposite of the wave propagation direction (kx), negative refraction may be possible with the structure shown in Fig. 1(a). To validate that the radiative mode is a physical solution, we calculated the imaginary part of kx for the RSPhP mode as positive. Thus, the surface mode decays as it propagates along the film, radiating light into the silicon along the way.

Fig. 3

Magnetic field distribution in three layer structure at 11.4 μm. Inset: Details of magnetic field inside SiC layer.

OE_51_12_120502_f003.png

In the remainder of this letter, we will use theoretical evidence to support the idea that incident light coupling to radiative surface phonon polariton modes can be accomplished in the Si/SiC/Di structure. We also use theoretical results to show that negative refraction is occurring in the SiC layer of the Si/SiC/Di structure. We can characterize this light coupling by examining the magnetic field in the Si/SiC/Di structure. Referring to the system shown in Fig. 1(b), we start by describing the magnetic field in each layer as

(4)

Hy2=ei(kxxωt)ekz2zzd,

(5)

Hy1=ei(kxxωt)(Aekz1z+Bekz1z)0zd,

(6)

Hy0=Cei(kxxωt)ekz0zz0.
The boundary conditions require that the magnetic fields and their derivatives be continuous at each interface, e.g., Hy0=Hy1, and (1/ϵ1)(Hy1/z)=(1/ϵ2)(Hy2/z) at z=0. By manipulating Eqs. (4) through (6) and applying the attendant boundary conditions, we find

(7)

A=[ϵ1kz0ϵ0kz1+1],B=[ϵ1kz0ϵ0kz11],andC=A+B.
The magnetic field is plotted for the case of 11.4 μm light in Fig. 3, where the gray region represents the SiC layer.

Figure 3 shows the infrared light propagating through the diamond layer as an oscillating field below the gray region. The gray region is the SiC layer. Figure 3 also shows the field exponentially decaying within the SiC layer, which is characteristic of a surface phonon polariton mode. The inset of Fig. 3 shows the details of the exponential decay in the SiC layer. Finally, another oscillating field is shown in the region above the SiC region, which demonstrates that the RSPhP mode is leaking into the silicon layer. Furthermore, the results show that the RSPhP mode is not self-sustaining, but depends on infrared light incident from the diamond side.11

We now, theoretically, show that the 11.4 μm light energy flux is negatively refracted in the Si/SiC/Di structure. The Poynting vector describes the light energy flux as

(8)

S=E×H=(HyEz)i+(HyEx)k,
and thus the direction of the energy flux in the material is given by

(9)

tanθ=ExHyHyEz,
where the magnetic fields are given by Eqs. (4) through (6), and Maxwell’s equations give the associated electric fields as

(10)

Exi=icωϵiHyizandEzi=kxcωϵiHyi,
which when used with Eq. (9) for each layer, we obtain

(11)

tanθ=ikz2kxz>d,

(12)

tanθ=ikz1(Bekz1zCekz1z)kx(Bekz1z+Cekz1z)0<z<d,

(13)

tanθ=ikz0kxz<d.
The energy flux directions for p-polarized 11.4 μm light in Si/SiC/Di is plotted in Fig. 4.

Fig. 4

Energy flux direction of 11.4 μm mode in each layer of the Si/SiC/Di structure.

OE_51_12_120502_f004.png

Figure 4 shows that the energy flux has a positive incidence angle at the SiC/Di interface which we calculate from Eq. (11) to be θ0=66deg. The energy flux then, dramatically, changes its angular direction in the SiC such that it is negatively refracted as it propagates from the diamond to the Si/SiC interface. The flux angle in the SiC layer, θ1, is calculated using Eq. (12). The real part of θ1 varys from 80deg to 73deg, and is plotted in Fig. 4. The curvature of the energy flux in the SiC layer is consistent with Fermat’s principle, since there is significant variation in refractive index from the diamond to Si layers. The energy flux emerges from the SiC/Si interface again with a positive angle calculated, with Eq. (13), to be θ2=74deg.

In conclusion, this letter, theoretically, predicts that a submicron thin film of SiC, bounded by silicon and diamond, may be used as a negative refraction material. The key mechanism is the radiative surface phonon polariton mode coupled to incident light at the SiC interfaces.

Acknowledgments

This work was supported by AFRL through a subcontract with Universal Technology Corporation, subcontract No. 11-S567-0015-02-C6.

References

1. 

S. KawataY. InouyeP. Verma, “Plasmonics for near-field nano-imaging and superlensing,” Nat. Photonics 3(7), 388–394 (2009).1749-4885http://dx.doi.org/10.1038/nphoton.2009.111Google Scholar

2. 

A. Berrieret al., “Negative refraction at infrared wavelengths in a two-dimensional photonic crystal,” Phys. Rev. Lett. 93(7), 73902 (2004).PRLTAO0031-9007http://dx.doi.org/10.1103/PhysRevLett.93.073902Google Scholar

3. 

E. Schonbrunet al., “Wave front evolution of negatively refracted waves in a photonic crystal,” Appl. Phys. Lett. 90(4), 041113 (2007).APPLAB0003-6951http://dx.doi.org/10.1063/1.2435344Google Scholar

4. 

J. Yaoet al., “Optical negative refraction in bulk metamaterials of nanowires,” Science 321(5891), 930–930 (2008).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1157566Google Scholar

5. 

V. Shalaev, “Optical negative-index metamaterials,” Nat. photonics 1(1), 41–48 (2007).1749-4885http://dx.doi.org/10.1038/nphoton.2006.49Google Scholar

6. 

A. Hoffmanet al., “Negative refraction in semiconductor metamaterials,” Nat. Mater. 6(12), 946–950 (2007).NMAACR1476-1122http://dx.doi.org/10.1038/nmat2033Google Scholar

7. 

J. Bossardet al., “Synthesizing low loss negative index metamaterial stacks for the mid-infrared using genetic algorithms,” Opt. Express 17(17), 14771–14779 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.014771Google Scholar

8. 

S. BasuM. Francoeur, “Near-field radiative transfer based thermal rectification using doped silicon,” Appl. Phys. Lett. 98(11) 113106 (2011).APPLAB0003-6951http://dx.doi.org/10.1063/1.3567026Google Scholar

9. 

C. FuZ. Zhang, “Further investigation of coherent thermal emission from single negative materials,” Nanoscale Microsc. Thermophys. Eng. 12(1), 83–97 (2008).NMTECW1556-7265http://dx.doi.org/10.1080/15567260701866744Google Scholar

10. 

H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings, Springer, Berlin (1988).Google Scholar

11. 

J. BurkeG. StegemanT. Tamir, “Surface-polariton-like waves guided by thin, lossy metal films,” Phys. Rev. B 33(8), 5186 (1986).PRBMDO0163-1829http://dx.doi.org/10.1103/PhysRevB.33.5186Google Scholar

Biographies and photograph of the authors are not available.

Olalekan S. Adewuyi, James S. Hammonds, "Negative refraction via radiative surface phonon polaritons in a silicon carbide-based multilayer structure," Optical Engineering 51(12), 120502 (4 December 2012). http://dx.doi.org/10.1117/1.OE.51.12.120502
Submission: Received ; Accepted
JOURNAL ARTICLE
3 PAGES


SHARE
KEYWORDS
Silicon carbide

Polaritons

Phonons

Negative refraction

Silicon

Diamond

Interfaces

Back to Top