We present a metamaterial that acts as a strongly resonant absorber at terahertz frequencies. Our design consists of a bilayer unit cell which allows for maximization of the absorption through independent tuning of the electrical permittivity and magnetic permeability. An experimental absorptivity of 70% at 1.3 terahertz is demonstrated. We utilize only a single unit cell in the propagation direction, thus achieving an absorption coefficient α=2000 cm-1. These metamaterials are promising candidates as absorbing elements for thermally based THz imaging, due to their relatively low volume, low density, and narrow band response.
©2008 Optical Society of America
The electromagnetic response of natural materials forms the basis for the construction of most modern optoelectronic devices. However, this EM response is not evenly distributed across the electromagnetic spectrum. At frequencies of a few hundred gigahertz and lower, electrons are the principle particles which serve as the workhorse of devices. On the other hand, at infrared through optical/UV wavelengths, the photon is the fundamental particle of choice. In-between these two fundamental response regimes there exists a region comparatively devoid of material response, commonly referred to as the “terahertz gap” (0.1-10 THz, λ=3mm-30µm) [1, 2]. Although enormous efforts have focused on the search for “terahertz” materials or alternative novel techniques to enable the construction of device components, much work remains. There is a wide range of natural phenomena that could be probed with terahertz (THz) devices. Specifically, a THz detector would be useful for imaging in areas such as biology [3, 4] and security [5, 6, 7, 8, 9].
Recently, there has been considerable effort to construct engineered electromagnetic materials for operation specifically within the void of natural material response described above [10, 11, 12, 13]. These artificial systems, called metamaterials (MMs), are composites whose EM properties originate from oscillating electrons in unit cells comprised of highly conductive and shaped metals such as gold or copper. The sub-wavelength unit cell is replicated to form a material, which allows for a designed resonant response of the metamaterial’s electrical and magnetic properties. Metamaterials can be regarded as effective media and characterized by a complex electric permittivity (ω)=ε 1(ω)+iε 2(ω) and complex magnetic permeability (ω)=µ 1(ω)+iµ 2(ω). Resonant structures that couple strongly to either the electric  or magnetic  fields have been demonstrated at terahertz. Significant growth in metamaterial research has been due to efforts to create negative refractive index (NRI) materials [14, 15, 16] and, more recently, invisibility cloaks [17, 18]. As such, the primary focus has been on the index of refraction defined as , where one desires n 1<0 for negative index or 0<n 1<1 for cloaks. To create such structures, it is important to minimize losses over the operating frequency range, which is associated with the imaginary portion of the index, and thus strive for n 2→0. Conversely, for many other applications it would be desirable to maximize the metamaterial loss which is an aspect of metamaterial research that, to date, has received very little attention. A recent example is the creation of a resonant high absorber which has been demonstrated at microwave frequencies . Such an absorber would be of particular importance at terahertz frequencies where it is difficult to find naturally occurring materials with strong absorption coefficients that are also compatible with standard microfabrication techniques. By fabricating bilayer metamaterial structures it becomes possible to simultaneously tune (ω) and (ω) such that a high absorptivity can be achieved. In principle, this tunability could lead to near unity absorptivity. In practice this is limited by achievable fabrication tolerances.
We present a first generation terahertz metamaterial absorber which achieves a resonant absorptivity of 70% at 1.3 THz. Given the 6 µm thickness of our metamaterial, this corresponds to a power absorption coefficient of α=2000 cm-1 which is significant at THz frequencies. The strong absorption coefficient makes this low volume structure a promising candidate for the realization of enhanced, spectrally selective, thermal detectors. A single unit of the absorber consists of two distinct metallic elements: an electrical ring resonator (ERR) Fig. 1 (a) and a split wire Fig. 1 (b). The electrical ring resonator (ERR) consists of two single split rings sitting back to back. The two inductive loops are of opposite handedness and thus couple strongly to a uniform electric field, and negligibly to magnetic fields [10, 20]. The magnetic component of light couples to both the center section of the electric resonator and the cut wire, thus generating antiparallel currents resulting in resonant µ(ω) response. The magnetic response can therefore be tuned independently of the electric resonator by changing the geometry of the cut wire and the distance between elements. By tuning each of the resonances it is possible to approximately match the impedance ( ) to free space, i.e. (ε=µ)⇒(Z=Z 0) and minimize the reflectance at a specific frequency. When the material is impedance-matched, the transmission is governed by the quantity n 2 kd, which can be simultaneously tuned with Z to obtain high absorption.
Computer simulations were performed using the commercial finite-difference time domain solver CST Microwave Studio™ 2006B and 2008. The metamaterials depicted in Fig. 1 were modeled as lossy gold with a conductivity of σ=1.0×107 S/m. The bottom substrate was modeled as gallium arsenide with a dielectric constant of 10.75. A 8 µm thick layer of dielectric, =3.5+i0.02 was used as the spacer between the two metallic metamaterial elements. We first investigated the S-parameters of transmission (S̃ 21) and reflection (S̃ 11) of a single unit cell with Perfect Electric (PE) and Perfect Magnetic (PM) boundary conditions along the x̂ and ŷ directions, respectively, (see Fig. 1). The absorptivity was calculated using the equation A=1-|S 21|2-|S 11|2. The electric and magnetic fields were examined at resonance to verify that we were coupling to the correct resonant mode of each metamaterial element. In Fig. 2, the resonant component of the electric field at resonance is plotted for the electric ring resonator (ERR) (a) and the split wire (b). The electric field is concentrated strongly in the gaps of the ring resonator and at the edges of the split wire in accord with previous results [10, 20]. Figures 2 (c) and (d) show a vector plot of the surface current density for the ERR and the split wire, respectively. Notice that at resonance currents are anti-parallel in the two metamaterial elements, which is the basis of the magnetic response and consistent with previous results [19, 22].
By changing the electric and magnetic resonances individually, we were able to create a condition such that the material was at an impedance near the free space value in a region of very low T(ω). The simulated transmission is relatively low across the entire range shown in Fig. 2 (e), whereas the reflectivity is relatively high except near the resonance at 1.12 THz where it drops to a value of 2%. Near unity absorbtion is theoretically possible and here we achieve a simulated value of 98% at 1.12 THz, as shown in Fig. 2(e). It should be noted that at THz frequencies the magnetic response of both natural materials and metamaterials is significantly weaker than the electric. Thus, matching the exact form of the ε(ω) and µ(ω) resonances becomes increasingly difficult at these and higher frequencies. Further, from a viewpoint of spectrally selective thermal imaging, it is desirable to have a narrow-band absorber. Thus we strive for two requirements of our metamaterials at our target frequency, i.e. Z=Z 0 so that R=0 and n 2 kd≫0 so that T→0. With these two goals we can achieve a significant and narrowband A(ω), but tolerance limits associated with microfabrication can reduce the absorptivity from unity.
We fabricated the metamaterial shown in Fig. 1 using a surface micromachining process, as shown in Fig. 3. A semi-insulating GaAs wafer was chosen because it is highly transmissive at THz frequencies. AZ5214e image reversal photoresist was spin-coated and patterned using standard photolithography. A 200 nm-thick Au/Ti film was E-beam evaporated to create the cut wire on the bottom layer. Lift off of the photoresist was achieved by rinsing in acetone for several minutes. The liquid polyimide, HD Microsystems™ PI-5878G, was spin-coated at 2,000 rpm on the GaAs wafer to form an insulating spacer with a thickness of 8 µm, and cured for five hours in an oven at 275°C in a nitrogen environment after the soft bake at 110°C for 6 minutes on a hot plate. AZ5214e image reversal photoresist was spin-coated, aligned, and patterned using standard photolithography. Another 200 nm-thick Au/Ti was E-beam evaporated as the material of the electric resonant ring on the top layer and then lifted off. Microscopic images of the as-fabricated samples are shown in Fig. 3 (right).
We experimentally verified the behavior of the absorber by measuring the transmission and reflectance of a large (1 cm×1 cm) planar array. We used an evacuated Fourier transform infrared (FTIR) spectrometer in the range from 300 GHz - 3 THz (10 - 100 cm-1) with 15 GHz (0.5 cm-1) spectral resolution. For transmission measurements the sample was mounted in the FTIR at normal incidence with the electric field perpendicular to the gap of ERR, as depicted in Fig. 1(a). Reflection was performed at an angle of 30° due to experimental limitations. The blue curves in Figs. 4(a) and (b) show, respectively, the measured reflectivity and transmission. Measured R(ω) and T(ω) differ significantly from that simulated as shown in Fig. 2(e). However, it should be noted that values used in simulation for the polyimide spacing layer were estimated based on published values at lower GHz frequencies. Further, the thickness of the polyimide layer was measured to be closer to 6µm rather than the 8µm used in simulations. Taking both of these factors into account, we are able to match measurements using an experimentally determined value for polyimide of ε=2.5+i0.2. The red curves are the corresponding simulations which are in good agreement with experiment. The simulated reflectance matches reasonably well near the resonance with slight deviations at lower and higher frequencies. The simulated transmission also agrees well with experiment, particularly in the vicinity of the resonance. Further, the simulated T(ω) reproduces the same qualitative features as experiment, including a distinct kink near 1.75 THz.
From experimental data presented in Fig. 4, the corresponding absorptivity is determined as shown in Fig. 5. These results demonstrate that the as-fabricated MMs have a strong resonance at around 1.3 THz and a high absorptivity of approximately 70%. Simulations show good agreement with experiment (red). The simulated absorptivity at resonance matches very well with the measured value. The off-resonance absorptivity is higher in experiment than in simulation due to differences in the experimental and simulated S-parameters. The response of absorber could be further improved through refinement and optimization of the fabrication process.
The absorber presented in this work absorbs strongly for light polarized along the x̂ direction, as shown in Fig. 1, but poor for ŷ polarized light, as shown in Fig. 6. In this polarization, the electric field is perpendicular to the center stalk of ERR, so an electric response cannot be driven. Similarly, there are no parallel wires for the magnetic field to develop a flux through and thus no net magnetic response. Such a polarization-sensitive device is desirable for both mm-wave and THz imaging as reflections from metallic objects often saturate the imager, thus significantly degrading it performance, in a problem known as “glint”. Additionally, polarization-sensitive detection has been shown to aid in discrimination of objects in a scene [23, 24]. However, by incorporating higher symmetry metamaterials, (similar to those presented here), a polarization-insensitive design could be acheived.
The performance of a THz radiation detector depends on the efficiency of converting radiation energy to an output signal. Therefore, maximizing the THz radiation absorption efficiency is integral to the development of a functional THz detector/imager. It is difficult to find strongly absorbing materials at THz frequencies that are compatible with standard photolithography. Thus, a potential application of these metamaterial structures is as the absorbing elements in thermal detectors. A strong absorption coefficient is also necessary to have a small thermal mass. This is important for optimizing the temporal response of thermal detectors. The metamaterial presented here has a 6 micron thick film and 70% absorptivity, which yields an absorption coefficient of 2000 cm-1. With better fabrication tolerances we could approach the simulated results (Fig. 2). This would yield an absorption coefficient three times greater than experimentally demonstrated here.
Many micro thermal detectors (typically based on bolometric detection and appropriately termed microbolometers) utilize several materials as the sensing element such as VOx , poly-Si-Ge , YBCO , or metal resistors such as, titanium  and niobium . Some of these materials are not fully compatible with microfabrication processing. For those that are, however, it is difficult and thus expensive to prepare or deposit a high quality film. Additionally, most of these materials show broad-band absorption. This limits potential applications, such as spectroscopic detection of explosive materials, which show unique responses at varied frequencies . The narrowband absorptivity of metamaterials presented here, on the other hand, enable spectrally selective detection. Furthermore, MMs are geometrically scalable and have been demonstrated over many decades of frequency. Thus, our results are not limited to terahertz frequencies and may be used over much of the electromagnetic spectrum. Another salient feature of the design presented here is that it may be combined with semiconducting materials or ferroelectrics to enable optically or electrically tunable frequency agile terahertz metamaterials. This would further permit a hyperspectral metamaterial focal plane array imager able to imaging over a relatively large band . Planar metamaterial absorbers consisting of different unit cells with distinct resonance frequencies may permit “multi-color” imaging.
In summary, we have demonstrated that the electromagnetic response of metamaterials can be tailored by manipulating the geometries of electric and magnetic resonators individually to create a highly selective absorber over a narrow band at THz frequencies. The successful demonstration of the high absorber holds great promise for future applications which includes metamaterial-based structures for creating a narrow-band, low thermal mass absorber as required for thermal sensing applications.
NIL, CMB, RDA and WJP acknowledge support from the Los Alamos National Laboratory LDRD program. This project has been supported in part by the DOD/Army Research Laboratory through grant W911NF-06-2-0040, DOE Los Alamos National Laboratory subcontract 50332-001-07 and DOE Los Alamos National Laboratory subcontract 50335-001-07. The authors would like to thank the Photonics Center at Boston University for technical support throughout the course of this research.
References and links
1. G. P. Williams, “Filling the THz gap - high power sources and applications,” Rep. Prog. Phy. 69, 301–326 (2006). [CrossRef]
2. M. Tonouchi, “Cutting-edge terahertz technology,” Nat. Photonics 1, 97–105 (2007). [CrossRef]
4. T. W. Crowe, T. Globus, D. L. Woolard, and J. L. Hesler, “Terahertz sources and detectors and their application to biological sensing,” Philosophical Transactions of the Royal Society of London A 362, 265–377 (2004).
5. F. Oliveira, R. Barat, B. Schulkin, F. Huang, J. Federici, and D. Gary, “Neural network analysis of terahertz spectra of explosives and bio-agents,” Proc. SPIE 5070, 60–70 (2003). [CrossRef]
6. D. Zimdars, “Fiber-pigtailed terahertz time-domain spectroscopy instrumentation for package inspection and security imaging,” Proc. SPIE 5070, 108–116 (2003). [CrossRef]
7. J. F. Federici, B. Schulkin, F. Huang, D. Gary, R. Barat, F. Oliveira, and D. Zimdars, “THz imaging and sensing for security applications - explosives, weapons, and drugs,” Semicond. Sci. Technol. 20, S266–S280 (2005). [CrossRef]
8. H.-B. Liu, Y. Chen, G. J. Bastiaans, and X.-C. Zhang, “Detection and identification of explosive RDX by THz diffuse reflection spectroscopy,” Opt. Express 11, 2549–2554 (2003).
9. J. Barber, D. E. Hooks, D. J. Funk, R. D. Averitt, A. J. Taylor, and D. Babikov, “Temperature-dependent far-infrared spectra of single crystals of high explosives using terahertz time-domain spectroscopy,” J. Phys. Chem. A 109, 3501–3505 (2005). [CrossRef]
10. W. J. Padilla, M. T. Aronsson, C. Highstrete, M. Lee, A. J. Taylor, and R. D. Averitt, “Electrically resonant terahertz metamaterials: theoretical and experimental investigations,” Phys. Rev. B 75, 041102R (2007). [CrossRef]
11. T. J. Yen, W. J. Padilla, N. Fang, D. C. Vier, D. R. Smith, J. B. Pendry, D. N. Basov, and X. Zhang, “Terahertz Magnetic Response from Artificial materials,” Science 303, 1494–1496 (2004). [CrossRef] [PubMed]
13. W. J. Padilla, A. J. Taylor, C. Highstrete, Mark Lee, and R. D. Averitt, “Dynamical electric and magnetic metamaterial response at terahertz frequencies,” Phys. Rev. Lett. 96, 107401 (2006). [CrossRef] [PubMed]
14. D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “A composite medium with simultaneously negative permeability and permitivity,” Phys. Rev. Lett. 84, 4184–4187 (2000). [CrossRef] [PubMed]
16. J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 86, 3996 (2000).
18. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial Electromagnetic Cloak at Microwave Frequencies,” Science 314, 977–980 (2006). [CrossRef] [PubMed]
19. N. I. Landy, S. Sajuyigbe, J. J. Mock, D. R. Smith, and W. J. Padilla, “A Perfect Metamaterial Absorber,” Submitted to Phys. Rev. Lett. [PubMed]
20. D. Schurig, J. J. Mock, and D. R. Smith, “Electric-field-coupled resonators for negative permittivity metamaterials,” Appl. Phys. Lett. 88, 041109 (2006). [CrossRef]
21. D. R. Smith, S. Schultz, P. Markos, and C. M. Soukoulis, “Determination of effective permittivity and permeability of metamaterials from reflection and transmission coefficients,” Phys. Rev. B 65, 195104 (2001). [CrossRef]
22. G. Dolling, M. Wegener, C. M. Soukoulis, and S. Linden, “Negative-index metamaterial at 780 nm wavelength,” Opt. Lett. 32, 53–55 (2007). [CrossRef]
23. F. Cremer, P. B.W. Schwering, W. Jong, K. Scutte, and A. N. Jong, “Infrared polarization measurements of targets and background marine environment,” Proc. SPIE 4370, 169–179 (2001). [CrossRef]
26. J. Wauters, “Doped silicon creates new bolometer material,” Laser Focus World 33, 145–149 (1997).
27. M. Almasri, D. P. Butler, and Z. Celik-Butler, “Self-supporting uncooled infrared bolometers with low thermal mass,” J. Microelectromechanical Syst. 10, 469–476 (2001). [CrossRef]
28. H. K. Lee, J. B. Yoon, E. Yoon, S. B. Ju, Y. J. Wong, W. Lee, and S. G. Kim, “A high fill-factor infrared bolometer using micromachined multilevel electrothermal structures,” IEEE Trans. Electron. Devices 46, 1489–1491 (1999). [CrossRef]
29. L. Baorino, E. Monticone, G. Amato, R. Steni, G. Benedetto, A. M. Rossi, V. Lacquaniti, R. Spagnolo, V. Lysenko, and A. Dittmar, “Design and fabrication of metal bolometers on high porosity silicon layers,” Microelectron. J. 30, 1149–1154 (1999). [CrossRef]
30. D. M. Mittleman, M. Gupta, R. Neelamani, R. G. Baraniuk, J. V. Rudd, and M. Koch, “Recent advances in terahertz imaging,” Appl. Phys. Lett. B 68, 1085–1094 (1999).
31. S. Nishizawa, K. Sakai, T. Hangyo, T. Nagashima, M. W. Takeda, K. Tominaga, A. Oka, K. Tanaka, and O. Morikawa, “Terahertz time-domain spectroscopy,” Terahertz Optoelectronics 97, 203–269 (2005). [CrossRef]
33. H.-T. Chen, J. F. O’Hara, A. K. Azad, A.J. Taylor, R. D. Averitt, D. B. Shrekenhamer, and W. J. Padilla, “Experimental Demonstration of Frequency Agile Terahertz Metamaterials,” Nat. Photonics, in press.