Narrowband thermal emission from a uniform tungsten surface critically coupled with a photonic crystal guided resonance

We numerically demonstrate narrowband thermal emission with unity emissivity peak in the near-infrared range by critically coupling a flat tungsten surface with guided resonances of a dielectric photonic crystal slab. The tungsten surface is separated from the photonic crystal slab by a vacuum gap. The structure possesses significant tunability for both the center frequency and the linewidth of the thermal emission band. Moreover, the tungsten surface, being un-structured, should exhibit enhanced thermal stability at elevated temperature as compared to tungsten nanostructures. c © 2016 Optical Society of America OCIS codes: (290.6815) Thermal emission; (350.4238) Nanophotonics and photonic crystals. References and links 1. J.-J. Greffet, R. Carminati, K. Joulain, J.-P. Mulet, S. Mainguy, and Y. Chen, “Coherent emission of light by thermal sources,” Nature 416, 61–64 (2002). 2. G. Chen and A. Shakouri, “Heat transfer in nanostructures for solid-state energy conversion,” J. Heat Transfer 124, 242–252 (2001). 3. V. Rinnerbauer, S. Ndao, Y. X. Yeng, W. R. Chan, J. J. Senkevich, J. D. Joannopoulos, M. Soljačić, and I. Celanovic, “Recent developments in high-temperature photonic crystals for energy conversion,” Energy Environ. Sci. 5, 8815– 8823 (2012). 4. M. U. Pralle, N. Moelders, M. P. McNeal, I. Puscasu, A. C. Greenwald, J. T. Daly, E. A. Johnson, T. George, D. S. Choi, I. El-Kady, and R. Biswas, “Photonic crystal enhanced narrow-band infrared emitters,” Appl. Phys. Lett. 81, 4685–4687 (2002). 5. S. Y. Lin, J. Moreno, and J. G. Fleming, “Three-dimensional photonic-crystal emitter for thermal photovoltaic power generation,” Appl. Phys. Lett. 83, 380–382 (2003). 6. M.-W. Tsai, T.-H. Chuang, C.-Y. Meng, Y.-T. Chang, and S.-C. Lee, “High performance midinfrared narrow-band plasmonic thermal emitter,” Appl. Phys. Lett. 89, 173116 (2006). 7. I. Celanovic, N. Jovanovic, and J. Kassakian, “Two-dimensional tungsten photonic crystals as selective thermal emitters,” Appl. Phys. Lett. 92, 193101 (2008). 8. J. A. Schuller, T. Taubner, and M. L. Brongersma, “Optical antenna thermal emitters,” Nat. Photonics 3, 658–661 (2009). 9. X. Liu, T. Tyler, T. Starr, A. F. Starr, N. M. Jokerst, and W. J. Padilla, “Taming the Blackbody with Infrared Metamaterials as Selective Thermal Emitters,” Phys. Rev. Lett. 107, 045901 (2011). 10. T. Inoue, T. Asano, M. De Zoysa, A. Oskooi, and S. Noda, “Design of single-mode narrow-bandwidth thermal emitters for enhanced infrared light sources,” J. Opt. Soc. Am. B 30, 165 (2013). 11. B. Zhao, L. Wang, Y. Shuai, and Z. M. Zhang, “Thermophotovoltaic emitters based on a two-dimensional grating/thinfilm nanostructure,” Int. J. Heat Mass Transfer 67, 637–645 (2013). 12. S. Molesky, C. J. Dewalt, and Z. Jacob, “High temperature epsilon-near-zero and epsilon-near-pole metamaterial emitters for thermophotovoltaics,” Opt. Express 21, A96–A110 (2013). 13. T. Inoue, M. D. Zoysa, T. Asano, and S. Noda, “Single-peak narrow-bandwidth mid-infrared thermal emitters based on quantum wells and photonic crystals,” Appl. Phys. Lett. 102, 191110 (2013). 14. Y. Guo, S. Molesky, H. Hu, C. L. Cortes, and Z. Jacob, “Thermal excitation of plasmons for near-field thermophotovoltaics,” Appl. Phys. Lett. 105, 073903 (2014). 15. P. N. Dyachenko, J. J. do Rosário, E. W. Leib, A. Y. Petrov, M. Störmer, H. Weller, T. Vossmeyer, G. A. Schneider, and M. Eich, “Tungsten band edge absorber/emitter based on a monolayer of ceramic microspheres,” Opt. Express 23, A1236 (2015). 16. D. Costantini, A. Lefebvre, A.-L. Coutrot, I. Moldovan-Doyen, J.-P. Hugonin, S. Boutami, F. Marquier, H. Benisty, and J.-J. Greffet, “Plasmonic Metasurface for Directional and Frequency-Selective Thermal Emission,” Phys. Rev. Appl. 4, 014023 (2015). Vol. 24, No. 26 | 26 Dec 2016 | OPTICS EXPRESS 29896 #275275 Journal © 2016 http://dx.doi.org/10.1364/OE.24.029896 Received 13 Sep 2016; revised 31 Oct 2016; accepted 31 Oct 2016; published 15 Dec 2016 17. Z. Wang, T. S. Luk, Y. Tan, D. Ji, M. Zhou, Q. Gan, and Z. Yu, “Tunneling-enabled spectrally selective thermal emitter based on flat metallic films,” Appl. Phys. Lett. 106, 101104 (2015). 18. T. Inoue, M. De Zoysa, T. Asano, and S. Noda, “Electrical tuning of emissivity and linewidth of thermal emission spectra,” Phys. Rev. B 91, 235316 (2015). 19. C. Blanchard, Q. Lévesque, D. Costantini, C. Jamois, J.-L. Leclercq, A.-L. Coutrot, F. Marquier, L. Milord, C. Grillet, H. Benisty, P. Viktorovitch, X. Letartre, and J.-J. Greffet, “Directional and selective mid-infrared thermal emitters for sensing applications,” in Advanced Photonics 2015, OSA Technical Digest (OSA, 2015), paper SeW2B.2. 20. A. Narayanaswamy and G. Chen, “Surface modes for near field thermophotovoltaics,” Appl. Phys. Lett. 82, 3544– 3546 (2003). 21. S. Basu, Z. M. Zhang, and C. J. Fu, “Review of near-field thermal radiation and its application to energy conversion,” Int. J. Energy Res. 33, 1203–1232 (2009). 22. E. Rephaeli and S. Fan, “Absorber and emitter for solar thermo-photovoltaic systems to achieve efficiency exceeding the Shockley-Queisser limit,” Opt. Express 17, 15145 (2009). 23. W. R. Chan, P. Bermel, R. C. N. Pilawa-Podgurski, C. H. Marton, K. F. Jensen, J. J. Senkevich, J. D. Joannopoulos, M. Soljačić, and I. Celanovic, “Toward high-energy-density, high-efficiency, and moderate-temperature chip-scale thermophotovoltaics,” PNAS 110, 5309–5314 (2013). 24. A. Lenert, D. M. Bierman, Y. Nam, W. R. Chan, I. Celanović, M. Soljačić, and E. N. Wang, “A nanophotonic solar thermophotovoltaic device,” Nat. Nano 9, 126–130 (2014). 25. K. A. Arpin, M. D. Losego, A. N. Cloud, H. Ning, J. Mallek, N. P. Sergeant, L. Zhu, Z. Yu, B. Kalanyan, G. N. Parsons, G. S. Girolami, J. R. Abelson, S. Fan, and P. V. Braun, “Three-dimensional self-assembled photonic crystals with high temperature stability for thermal emission modification,” Nat. Commun. 4, 2630 (2013). 26. I. Celanovic, D. Perreault, and J. Kassakian, “Resonant-cavity enhanced thermal emission,” Phys. Rev. B 72, 075127 (2005). 27. M. Ghebrebrhan, P. Bermel, Y. X. Yeng, I. Celanovic, M. Soljačić, and J. D. Joannopoulos, “Tailoring thermal emission via Q matching of photonic crystal resonances,” Phys. Rev. A 83, 033810 (2011). 28. Y. X. Yeng, M. Ghebrebrhan, P. Bermel, W. R. Chan, J. D. Joannopoulos, M. Soljačić, and I. Celanovic, “Enabling high-temperature nanophotonics for energy applications,” Proc. Natl. Acad. Sci. U. S. A. 109, 2280–2285 (2012). 29. V. Liu and S. Fan, “S4 : A free electromagnetic solver for layered periodic structures,” Comput. Phys. Commun. 183, 2233–2244 (2012). 30. O. Ilic, P. Bermel, G. Chen, J. D. Joannopoulos, I. Celanovic, and M. Soljačić, “Tailoring high-temperature radiation and the resurrection of the incandescent source,” Nat. Nano 11, 320–324 (2016). 31. S. Fan, W. Suh, and J. D. Joannopoulos, “Temporal coupled-mode theory for the Fano resonance in optical resonators,” J. Opt. Soc. Am. A 20, 569 (2003). 32. J. R. Piper and S. Fan, “Total Absorption in a graphene Monolayer in the optical regime by critical coupling with a photonic crystal guided resonance,” ACS Photonics 1, 347–353 (2014). 33. I. Celanovic, F. O’Sullivan, M. Ilak, J. Kassakian, and D. Perreault, “Design and optimization of one-dimensional photonic crystals for thermophotovoltaic applications,” Opt. Lett. 29, 863 (2004). 34. S. Fan and J. D. Joannopoulos, “Analysis of guided resonances in photonic crystal slabs,” Phys. Rev. B 65, 235112 (2002). 35. W. G. Spitzer, D. Kleinman, and D. Walsh, “Infrared properties of hexagonal silicon carbide,” Phys. Rev. 113, 127–132 (1959). 36. Z. Yu, N. P. Sergeant, T. Skauli, G. Zhang, H. Wang, and S. Fan, “Enhancing far-field thermal emission with thermal extraction,” Nat. Commun. 4, 1730 (2013).


Introduction
Conventional thermal light emitting sources, such as the sun or the filament of an incandescent light bulb, typically generate broadband and incoherent electromagnetic radiations [1].Yet in various applications of light sources, including sensors, spectroscopy and energy conversion [2,3], it is essential to have a thermal light source with a narrowband thermal emission.Therefore, in recent years, there have been extensive works aiming to achieve narrowband thermal emission using a wide variety of nanophotonic structures [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19].In addition, it is desirable to control both the frequency and the linewidth of the narrowband emission while maintaining high peak emissivity.For example, in solar thermophotovoltaics [3,[20][21][22][23][24], the emission wavelength needs to matches the band gap of the photovoltaic cell.And moreover there is an optimal emission bandwidth in order to achieve maximum power output [22].
In addition to the need for thermal emission control and energy conversion applications, there is a challenge to design thermal emitters that maintain structural stability at elevated temperatures.From basic thermodynamic efficiency considerations, the energy conversion efficiency improves Thermal emission spectra from the structure shown in Fig. 1(a) with tunable linewidth.Here we vary the geometric parameters of the photonic crystal slab and then tune the vacuum gap size to achieve critical coupling.The structured parameters for each of the curves are as follows: yellow curve: d=0.06µm, g=0.12µm; red curve: d=0.04µm, g=0.16µm; blue curve: d=0.02µm, g=0.23µm.Structures with smaller gap size has higher intrinsic losses, leading to larger linewidths.(c) Thermal emission spectrum from a tungsten photonic crystal slab as shown in the inset.The structure has a period of 2µm and a grating depth of 3µm.(d) Thermal emission spectra from structures consisting of multiple layers of dielectric films placed on top of a flat tungsten surface separated by a vacuum gap.We use two sets of dielectrics to form the top mirror.The first set consists of high index layers with n H = 3.34, low index layer with n L = 1.45, with 3 layers in total (HLH), and the second set consists of high index layers with n H = 3.57, low index layers with n L = 2.90, with 5 layers in total (HLHLH).The emissivity spectra and hence the emission linewidth of the two structures are almost the same.
as the temperature of the thermal source increases.For solar thermophotovoltaics applications, the optimal temperature for the emitter may exceed 2000K [22].Since tungsten has a high melting point, there have been significant efforts in exploring tungsten nanostructures for the control of thermal emissivity [3,5,7,11,15,16,22,[24][25][26][27][28].Unlike bulk tungsten however, it has been experimentally observed that tungsten nanostructure is far less stable at elevated temperature [25].There is therefore a need to develop a mechanism to control the thermal emissivity of a flat tungsten surface without nanostructuring the tungsten.
In this paper, using the rigorous coupled wave analysis [29], we numerically demonstrate the control of thermal emissivity from a flat tungsten surface, by critical coupling to a dielectric photonic crystal guided resonance in the near field.We show that narrowband thermal emission with peak emissivity approaching unity can be achieved in these structures.Both the bandwidth and the center frequency of the thermal emission can be controlled by designing the structure.Importantly, in this design, the photonic crystal, which provides the capability for spectral control, can in principle be completely lossless and need not be maintained at high temperature.The tungsten region, which is the source of thermal emission and therefore needs to be at high temperature, is flat and therefore should not suffer from the thermal instability issues of tungsten nanostructures.While it has been previously shown that narrowband emission from a flat tungsten surface can be achieved by placing a one-dimensional crystal on top of it [26,30], our design here possesses a much wider tunability in its emission bandwidth.
The paper is organized as follows: In Section 2, we briefly summarize previous methods on achieving narrowband thermal emission and point out the advantages of our design.In Section 3, we discuss the simulation results for our design using a one-dimensional dielectric photonic crystal slab in the vicinity of the flat tungsten surface.We first briefly review the physics of guided resonances.We then show how to tune the resonant frequency and linewidth of thermal emission based on the properties of guided resonances.We present the results for the emissivity spectra of two-dimensional photonic crystal slab coupled to tungsten surface in Section 4 and conclude in Section 5.

The concept of critical coupling and its applications for achieving narrowband thermal emission
To achieve narrowband thermal emission, it is natural to consider structures that support resonances.From coupled mode theory [31,32], for a single-mode resonator coupled to a single input/output port, its spectral absorption coefficient A is where ω 0 is the resonant frequency, δ is the intrinsic loss rate of the resonator that arises from the material loss inside the resonator, and γ is the external leakage rate that arises from the coupling of the resonance to the output port.From this equation, we see that the peak absorptivity reaches unity at the resonance frequency ω 0 , if the resonance satisfies the condition that the external leakage and intrinsic loss rates are the same (γ = δ).This condition of γ = δ is the critical coupling condition.At critical coupling, the resonance has a total linewidth of 2δ.
Based on the concept of critical coupling, we now briefly review several existing approaches as well as our proposed approach in this paper for achieving narrowband thermal emission from tungsten structures.In all the structures shown in Fig. 1, one assume a tungsten region that is sufficiently thick, such that a resonance created near the top surface can only emit to the half space on top of the tungsten structure.And moreover, the structures are either uniform (Fig. 1(d)), or have a structural period that is sub-wavelength at the peak emission wavelength (Fig. 1(a) and  (c)), such that no diffraction occurs.In these structures, for a normally incident light the reflected wave can only occur in the normal direction, the resonance therefore couples only to a single input/output port.The emission property of all these structures to the normal direction therefore can all be understood in terms of a single resonance coupling to an input/output port.For these structures then, peak unity emissivity is achieved by satisfying the critical coupling condition.
Most existing approaches for designing tungsten thermal emitter can be broadly categorized into two categories.In the first category, one patterns either a one-dimensionally or twodimensionally periodic structure on a tungsten surface [7,23,27], forming a tungsten photonic crystal slab.The structure shown in Fig. 1(c) represents an example with a one-dimensional periodic structure introduced into tungsten.The tungsten photonic crystal slab has a period of 2µm and grating depth of 3µm.This structure is designed to achieve unity emissivity at 0.6eV.In general, in these structures, the peak frequency of the resonance can be controlled by changing the period.However, the achievable linewidth is relatively broad.The resonant modes in these structures are localized in the tungsten region as shown in Fig. 2(b), therefore the intrinsic loss rate is almost entirely controlled by the properties of tungsten with relatively little tunablity.Also, since the loss of tungsten is substantial, the linewidth of the resonance is also relatively broad.In addition, such patterned structures are not as stable as their bulk counterparts at high temperatures.It was experimentally shown that tungsten-coated colloidal crystal nanostructures were completely destroyed at 1500K [25].Thus protective coatings are often required to enhance the thermal stability.
In the second category, a dielectric multilayer stack is put on top of a flat tungsten surface with a vacuum gap in between [22,26,33].The dielectric multilayer consists of alternating layers of a high index material and a low index material.In our example a quarter-wave layer design is chosen for illustration purposes.This structure is capable of realizing narrowband thermal emission with tunable peak frequency.The linewidth is much narrower compared to the structure shown in Fig. 1(c).This can be understood by examining the corresponding field distributions.For the structure shown in Fig. 1(c), its electric field strongly penetrates into the tungsten photonic crystal slab as shown in Fig. 2 is a Fabry-Perot mode formed through reflection between the tungsten layer and the multi-layer dielectric mirror as shown in Fig. 2(c).The field is therefore located somewhat away from the tungsten region, resulting in a smaller intrinsic loss rate and a narrower linewidth.Moreover, for such a Fabry-Perot mode, the intrinsic linewidth is controlled by the absorption coefficient for light normally incident from vacuum onto the flat tungsten surface, which is a constant at a given resonant frequency.Therefore, at critical coupling, the linewidth of the lowest-order Fabry-Perot resonance at a given frequency is completely fixed, as explained by Eq. ( 1).As an illustration, in Fig. 1(d), we choose two sets of dielectrics to form the top dielectric multilayer.The first set consists of high index layers with n H = 3.34, low index layer with n L = 1.45, with 3 layers in total (HLH), and the second set consists of high index layers with n H = 3.57, low index layers with n L = 2.90, with 5 layers in total (HLHLH).Despite the differences in dielectric constants and number of layers, the emissivity spectra of the two structures are almost the same, and the linewidths are nearly identical, which is consistent with the argument above.
In this work we also use the concept of critical coupling to achieve narrowband thermal emission.The geometry is shown in Fig. 1(a), where we place a silicon carbide (SiC) photonic crystal slab on top of a tungsten slab separated by vacuum gap.A similar geometry of a photonic crystal slab separated from a metal surface by a vacuum gap has been recently used as a mid-infrared thermal emitter [19].The SiC photonic crystal slab supports guided resonances.Unlike the two previous cases above, the energy density of the resonances is located in the SiC region away from the tungsten surface.Since in the wavelength range of interest SiC is largely transparent, the intrinsic loss rate of the guided resonances arises from the overlap of the evanescent tail of these resonances with the tungsten region and thus is very sensitive to the size of the vacuum gap.At a given vacuum gap size, which largely fixes the intrinsic loss rate of the resonance, one can then achieve critical coupling by adjusting the geometric parameters of the photonic crystal slab which control the external loss rate.Alternatively, one can also vary the vacuum gap size to achieve critical coupling for given geometric parameters of the photonic crystal slab.As a result, unlike previous methods, here at critical coupling we can achieve a much wider range of resonance linewidth, including resonances with a much narrower linewidth, as compared to previous methods.Also, similar to the structure shown in Fig. 1(d), in our case the tungsten itself needs not to be patterned, and the structure that is patterned, i.e. the SiC region, needs not to be heated.As a result the structure should be more thermally stable as compared to the tungsten photonic crystal slab shown in Fig. 1(c).

Critical coupling with one-dimensional photonic crystal slabs
In this Section we provide an in-depth discussion of the thermal emission properties of the structure shown in Fig. 1(a).In Section 3.1, we first review the properties of guided resonance of the SiC photonic crystal slab.In Section 3.2, we then utilize the properties of such photonic crystal slab in the vicinity of a tungsten slab to enable critical coupling.Finally, in Sections 3.3 and 3.4, we show that the critical coupling can be maintained in this structure as we vary either the resonance linewidth or the resonant frequencies.

A brief review of guided resonances
We first examine the modal properties of a uniform dielectric slab, to which the spectral features of guided resonances of photonic crystal slabs are closely related [34].In Fig. 3(a), we plot the four lowest order modes, TE 0 , TM 0 , TE 1 , and TM 1 of a SiC slab with a thickness of 0.4µm.In the modal calculation, the dielectric constant of SiC is fixed to be 6.6, while in all the absorption calculations, we have considered both the dispersion and loss of SiC, with dielectric constant data taken from Ref. [35].Periodic patterning of the slab breaks the continuous translational symmetry, and as a result some of the guided modes can couple to free space, resulting in the creation of guided resonances.
For applications such as thermophotovoltaics, one typically would like to achieve narrowband emission at the band gap of the photovoltaic cell, while suppress emission at all frequencies.Therefore, in our design, we focus on the lowest order guided resonance.To design a photonic crystal slab that has the frequency of its lowest order guided resonance located at ω 0 , we examine the lowest order band of the underlying uniform structure and identify the wavevector k associated with the mode of frequency ω 0 in this band.We then choose the period of photonic crystal slab to be approximately 2π/k.This enables the phase matching of the guided mode with normally incident plane wave.Moreover this period, by construction, is smaller than the free space wavelength, ensuring that the guided mode is phase matched only to the zeroth diffraction order, which is important in order to achieve critical coupling with such guided resonances.
For the uniform slab considered here, the lowest order mode is the TE 0 mode (Fig. 3).To place the lowest order guided resonance at the target frequency of 0.6eV, we choose a period of 1µm following the procedure outlined above.Such a guided resonance couples to s-polarized wave, which has electric field along the y-direction in Fig. 1(a), and as a result the incident wave is chosen to be s-polarized for one-dimensional photonic crystal slabs.The presence of a guided resonance in a photonic crystal slab manifests in the transmission spectrum as a Fano line shape superimposed on an otherwise smooth background.In Fig. 3(b), we show the transmission spectrum of a photonic crystal slab with a period of 1µm and a grating depth of 0.04µm.We do observe a Fano resonance at the target frequency, and no resonances below it.

Critical coupling
We next place the SiC photonic crystal slab as discussed above, with a period of 1µm and a grating depth of 0.04µm, on top of a flat tungsten surface separated by a vacuum gap (Fig. 1(a)).We find that the critical coupling into the lowest order guided resonance is achieved with a gap size of 0.16µm.In Fig. 4(a), we show the absorption spectrum of the tungsten region in this structure at critical coupling.We find a distinct peak with near unity peak absorption at the target energy of 0.6eV, and no resonances existing at lower frequencies.In Fig. 4(a), we also compare the spectrum with the SiC photonic crystal slab to that of a flat tungsten surface.Away from the two main resonances near 0.6eV, the structure with photonic crystal slab has lower absorptivity as compared to flat tungsten surface, in the frequency range below 0.55eV, and almost in the entire frequency range from 0.8eV to 1.4eV which is above the main resonance.Such a suppression of emissivity away from the resonance is advantageous from efficiency considerations in thermophotovoltaic systems.In Fig. 2(a), we plot the electric field intensity when the structure is excited by an incident plane wave with an energy of 0.6eV.We observe the excitation of the guided resonance as evident from strong field concentration inside the SiC regions.
To achieve critical coupling, the key control parameter is the gap size.In Fig. 4(b), we show the absorption spectra for the structures that are the same as the one shown in Fig. 4(a), except that the gap sizes take the values of 0.10µm, 0.16µm and 0.22µm.As the gap sizes increases, the spectra show the transition from under coupling, through critical coupling, to the over coupling regime.
In Fig. 4(c), we plot the absorptivity for the s-polarization as a function of photon energy and angle of incidence.We see that the peak of absorptivity can shift significantly as the incident angle varies.In particular, the resonance at 0.6eV splits into two modes as the angle of incidence increases.Also, as angle of incidence increases, some of these modes will deviate from the critical coupling condition.In Fig. 4(d), we plot the angle averaged absorptivity spectrum as a function of photon energy.The average is defined as ᾱ = 2 π π/2 0 α(θ) cos θdθ, where θ is the polar angle.We see that, in spite of the strong angular dependence of the resonances, the angle averaged absorptivity still exhibits a spectral peak.Further improvement of angle selectivity can be accomplished, for example, by following the method as discussed in Ref. [16].

Control of linewidth
A key advantage of our design is the capability to achieve a wide variety of emission linewidth while maintaining critical coupling.In this system, the intrinsic loss rate of a guided resonance is strongly influenced by the gap size, whereas the external loss rate is largely controlled by the grating depth.Therefore, to vary the emission linewidth, we can vary the grating depth and then adjust the gap size to maintain the critical coupling condition.In this Section, when varying the geometric parameters, we fix the grating width (w) to be half of the period, and the effective thickness of the slab (b + d/2) to be 0.4µm.
As an illustration, in Fig. 1(b), we vary the grating depth d at 0.02µm, 0.04µm and 0.06µm, the gap sizes at critical coupling are 0.23µm, 0.16µm and 0.12µm, respectively.All three structures achieve critical couping at the photon energy of 0.6eV.The linewidth of the resonance gets larger as the grating depth increases.Alternatively, we can also choose a particular gap size, and then adjust the grating depth to achieve critical coupling.In Fig. 5(a), we show the linewidth at gap sizes of 0.05µm (blue curve), 0.10µm (red curve), 0.20µm (yellow curve), where the grating depths at critical coupling are 0.147µm, 0.073µm and 0.027µm, the periods are 1.048µm, 1.013µm and 0.991µm, respectively.The linewidth of the resonance gets smaller as the vacuum gap size increases.
The examples above demonstrate the ability to design a wide range of resonance linewidth while maintaining critical coupling in our system.In the two examples as illustrated in Fig. 1(b) and Fig. 5(a), we are able to straightforwardly vary the linewidth by more than a factor of five.
Here we comment on the range of emission linewidth that can be achieved in our structure.In principle there is no lower limit of the external loss rate because an unpatterned slab has an external loss rate of exactly zero.Also, assuming the photonic crystal slab itself is lossless, and the loss arises only from the tungsten, there is then no lower limit of the intrinsic loss rate, either, since we can always increase the vacuum gap between the photonic crystal slab and tungsten.Thus theoretically we can achieve resonances with arbitrarily small linewidth while maintaining critical coupling.On the other hand, there is an upper limit for the linewidth that can be achieved.Fundamentally, the upper limit arises from the absorption coefficient of tungsten at the target photon energy.And moreover, in our planar geometry, there is a large impedance mismatch between vacuum and tungsten which limits the field penetration into the tungsten region.In our numerical simulation, the largest full width at half maximum for the resonance is approximately 0.02eV.This is comparable to linewidth in the multilayer structure shown in Fig. 1 since both have a flat tungsten surface and is significantly smaller than the tungsten photonic crystal slab shown in Fig. 1(c), where the fields penetrate significantly into the tungten region as shown in Fig. 2(b).

Control of peak frequency
The resonant frequency can also be varied over a broad range by varying the period of the photonic crystal slab.As the period a increases, the wavevector of the lowest order guided mode that the normal incident wave can couple to, which is 2π/a, decreases.As a result, the resonant frequency decreases as well, as seen from the modal dispersion curves shown in Fig. 3(a).This redshift of resonant frequency affects the intrinsic loss rate, since both the material loss of tungsten, and the penetration of the evanescent waves into the tungsten region, vary as a function of frequency.Therefore, in order to maintain critical coupling, the external loss rate should be designed correspondingly to match the change of intrinsic loss rate.
We show an example for such design in Fig. 5(b), where we vary the period of the photonic crystal slab from 0.8µm (red curve), through 1.0µm (blue curve), to 1.5µm (yellow curve).For each period, we vary the grating depth as well as the vacuum gap size until the critical coupling is reached.As a result of such a tuning process, we are able to achieve critical coupling at the wavelengths of 0.71eV, 0.60eV, and 0.46eV, respectively.This examples indicates that our structure can be designed to exhibit narrowband emission with unity-emissivity peak located at a wide range of energies.Such a capability of wavelength tunability is potentially important for the design of narrowband thermal sources.

Critical coupling with two-dimensional photonic crystal slabs
In previous sections, we have illustrated the concepts by focusing on a class of thermal emitters where the photonic crystal structure has one-dimensional periodicity.Such a system has strong polarization dependency.On the other hand, one can apply the same concept to create thermal emitters that achieve narrowband emission with unity emissivity for both polarizations at normal emission direction.For this purpose, we instead consider a photonic crystal slab consisting of a square array of air holes introduced into the SiC slab, as illustrated in Fig. 6(a) where we show four unit cells.At normal incidence the two polarizations are degenerate by symmetry.The slab is again placed on top of a flat tungsten surface separated by a vacuum gap.
In general, such a photonic crystal slab supports both TE-like and TM-like guided resonances.These guided resonances correspond to TE and TM modes in the dielectric slab waveguides.The TE and TM-like resonances have mostly in-plane electric or magnetic field within the slabs, respectively [34].Both of these resonances can couple to normally incident plane wave.In addition, the frequencies of these resonances are typically close to each other, since the effective index for the TE and TM modes of the same order in the dielectric slab waveguides are relatively close to each other for modes that are strongly confined within the slab.Therefore, for the lowest energy emission peak one may expect to see two resonances adjacent to each other.On the other hand, the quality factors of these resonances are in general quite different, and therefore by controlling the intrinsic loss rates of the resonance one can design the system such one of the resonances reaches critical coupling and hence has unity emissivity at its peak, whereas the other resonance does not satisfy the critical coupling condition and hence has a suppressed emissivity.
As an illustration, we consider two examples shown in Fig. 6.In Fig. 6(b), the thickness of the SiC slab is 0.45µm, the period is 1.08µm, the radius of vacuum hole is 0.18µm, the vacuum gap size is 0.07µm.In this structure, the normal incident wave is critically coupled to the lowest-order TE-like resonance, which results in a unity emissivity peak at 0.6eV, whereas the lowest-order energy(eV) Narrowband thermal emission with unity peak emissivity by critical coupling to the lowest order TE-like guided resonance located at 0.6eV.The thickness of the SiC slab is 0.45µm, the period is 1.08µm, the radius of vacuum hole is 0.18µm, the vacuum gap size is 0.07µm.(c) Narrowband thermal emission with unity peak emissivity by critical coupling to the lowest order TM-like resonance located at 0.6eV.The thickness of the SiC slab is 0.58µm, the period is 1.03µm, the radius of vacuum hole is 0.10µm, the vacuum gap size is 0.13µm.TM-like resonance located at 0.637eV is under coupled and hence has a strongly suppressed thermal emissivity.In Fig. 6(c), the thickness of the SiC slab is 0.58µm, the period is 1.03µm, the radius of vacuum hole is 0.10µm, the vacuum gap size is 0.13µm.Here the lowest order TM-like resonance satisfies the critical coupling condition, resulting in an unity emissivity peak centered at 0.6eV, whereas the lowest-order TE-like resonance, located at 0.57eV, is over coupled with suppressed emissivity.These examples show that one can achieve critical coupling with two-dimensional photonic crystal slab system, and moreover generate a single emissivity peak that has its energy well separated from all other emissivity peaks in the system.
In general, a two-dimensional lattice has a more extensive set of reciprocal lattice vectors.And consequently there are more resonant peaks in the emissivity spectrum, as we indeed observe when we compare the plots shown in Fig. 6 with any of the spectra that we show for the one-dimensional case.

Conclusion
In summary, we show that narrowband thermal emission can be achieved, by placing a photonic crystal slab on top of a flat tungsten surface separated by a vacuum gap, and by designing the parameters of the structure such that the lowest order guided resonance reaches critical coupling.This method achieves strong modifications of the spectral emissivity without the need to create tungsten nanostructures, which is expected to significantly improve the structure stability of tungsten based thermal emitters.This class of structures also offer wide tunablity range for both the linewidth and the position of the emissivity peak while maintaining unity emissivity.Such a capability is of importance for thermal source design in general.While for illustration purposes we have focused on creating a narrowband tungsten emitter that emits strongly near an energy of 0.6eV, the same principle can be applied for other wavelength ranges and other material systems as well.Finally, the aspect of near-field coupling between a thermal emitter and a dielectric structure has been previously used to demonstrate a thermal exaction scheme [36].Therefore, an interesting future opportunity is to combine the concept here with the thermal exaction scheme to achieve narrowband thermal emission with power beyond the apparent blackbody limit.

Fig. 1 .
Fig. 1.(a) Schematic of the emitter structure consisting of a photonic crystal slab with one-dimensional periodicity separated from a flat tungsten surface by a vacuum gap.(b)Thermal emission spectra from the structure shown in Fig.1(a) with tunable linewidth.Here we vary the geometric parameters of the photonic crystal slab and then tune the vacuum gap size to achieve critical coupling.The structured parameters for each of the curves are as follows: yellow curve: d=0.06µm, g=0.12µm; red curve: d=0.04µm, g=0.16µm; blue curve: d=0.02µm, g=0.23µm.Structures with smaller gap size has higher intrinsic losses, leading to larger linewidths.(c) Thermal emission spectrum from a tungsten photonic crystal slab as shown in the inset.The structure has a period of 2µm and a grating depth of 3µm.(d) Thermal emission spectra from structures consisting of multiple layers of dielectric films placed on top of a flat tungsten surface separated by a vacuum gap.We use two sets of dielectrics to form the top mirror.The first set consists of high index layers with n H = 3.34, low index layer with n L = 1.45, with 3 layers in total (HLH), and the second set consists of high index layers with n H = 3.57, low index layers with n L = 2.90, with 5 layers in total (HLHLH).The emissivity spectra and hence the emission linewidth of the two structures are almost the same.

Fig. 2 .
Fig. 2. Distributions of the electric field intensity at critical coupling of (a) structure shown in Fig. 1(a) with a gap size of 0.16µm, (b) tungsten photonic crystal slab in Fig. 1(c), and (c) one-dimensional multilayer structure in Fig. 1(d).The dashed lines represent the outline of various elements of the structures.

Fig. 3 .
Fig. 3. (a) Optical modes of a uniform SiC slab with a thickness of 0.4µm.(b) Transmission spectrum of a SiC photonic crystal slab with a period of 1µm and a grating depth of 0.04µm.The spectrum has its lowest order resonance located around 0.6eV.

Fig. 4 .
Fig. 4. (a) Narrowband thermal emission peaking at 0.6eV.Parameters: period a = 1.00µm, grating depth d = 0.04µm, vacuum gap size g = 0.16µm.(b) The system goes from under coupling (yellow curve, g = 0.10µm), through critical coupling (blue curve, g = 0.16µm), to over coupling (red curve, g = 0.22µm) as the gap size increases.(c) Absorption spectrum as a function of photon energy and angle of incidence for s-polarized wave, where the structure has same parameters as subplot (a).(d) Angle averaged absorptivity as a function of photon energy for structure in subplot (c).

Fig
Fig. (a) Tuning of linewidth by varying vacuum gap size.Here we vary the vacuum gap size and then adjust the geometric parameters of the photonic crystal slab to achieve critical coupling.The structured parameters for each of the curves are as follows: yellow curve: g = 0.20µm, d = 0.027µm, a = 0.991µm; red curve: g = 0.10µm, d = 0.073µm, a = 1.013µm; blue curve: g = 0.05µm, d = 0.147µm, a = 1.048µm.(b) Tuning of peak position by varying the period.The structured parameters for each of the curves are as follows: red curve: a = 0.8µm, g = 0.18µm, d = 0.02µm; blue curve: a = 1.0µm, g = 0.16µm, d = 0.04µm; yellow curve: a = 1.5µm, g = 0.23µm, d = 0.24µm.

Fig. 6 .
Fig. 6.(a) Schematic of the emitter structure consisting of a photonic crystal slab with two-dimensional periodicity separated from a flat tungsten surface by a vacuum gap.(b)Narrowband thermal emission with unity peak emissivity by critical coupling to the lowest order TE-like guided resonance located at 0.6eV.The thickness of the SiC slab is 0.45µm, the period is 1.08µm, the radius of vacuum hole is 0.18µm, the vacuum gap size is 0.07µm.(c) Narrowband thermal emission with unity peak emissivity by critical coupling to the lowest order TM-like resonance located at 0.6eV.The thickness of the SiC slab is 0.58µm, the period is 1.03µm, the radius of vacuum hole is 0.10µm, the vacuum gap size is 0.13µm.