Mid-infrared plasmonic multispectral filters

A miniaturized mid-infrared spectral analyzer will find a wide range of applications as a portable device in non-invasive disease diagnosis, environmental monitoring, food safety and others. In this work, we report an integrated spectral analyzer that can be constructed by using Au subwavelength hole arrays as multispectral filters. The hole arrays were fabricated with CMOS compatible processes. The transmission peak of the subwavelength hole arrays is continuously tuned from 3 μm to 14 μm by linearly increasing the periodicity of the holes in each array. Fourier transform infrared (FTIR) microscopy was applied to spatially map out the transmission of the hole arrays. The results show that each hole array can selectively allow for transmission at a specific wavelength. We further constructed an IR spectral analyzer model based on the microhole multispectral filters to retrieve IR spectral information of two test samples. Our experimental results show that the spectra from the integrated spectral analyzer follow nearly the same pattern of the FTIR spectra of the test samples, proving the potential of the miniaturized spectral analyzer for chemical analysis.

Most chemicals have distinct absorption signatures in mid-infrared spectral range. Spectral analysis provides a highly sensitive and selective method for chemical detection, which may find a wide range of applications in gas sensing[1] , noninvasive disease diagnosis[2] ,[3] , security monitoring[4] and others[5] -[7] . For example, two-channel spectral analyzers (InfraTec) for gas sensing have been commercialized by using two separate traditional Fabry-Perot spectral filters. The filtering spectrum of Fabry-Perot filters is determined by the cavity thickness. Each cavity thickness requires one step of microfabrication. It is extremely challenging to integrate a large number of Fabry-Perot filters on chip for constructing a smart multispectral analyzer. More importantly, the “stop band” of Fabry-Perot filters is only a few micrometers wide[8] due to the small contrast in refractive index of the available materials that are transparent in mid-infrared. It is difficult to avoid the cross-talk of the spectral filters in the wide mid-infrared spectrum range (3 µm-30 µm).

Surface plasmon resonances (SPR) of noble and transition metals are well known for their capability to enhance light intensity in visible and near infrared range[9] ,[10] . Nanostructured metal films have been widely used to construct RGB color filters in visible spectrum[11] ,[12] , aiming to replace the classical organic color filters in digital cameras. Analogically, metallic microhole arrays (MHAs) were also presented for enhancing transmission of mid-infrared light[13] ,[14] , and even terahertz radiation[15] ,[16] . Transmission peaks at specific wavelengths could be achieved by properly designing metal MHAs, providing a new approach to realize multispectral infrared photodetectors[17] . Here we report plasmonic multispectral filters in mid-infrared that are constructed using a subwavelength hole array in evaporated Au thin film. Nanohole arrays in Au thin film were previously reported as color filters in visible spectrum[18] . The plasmonic resonance in the metal film can be shifted from visible to mid-infrared range (3 µm to 14 µm) simply by increasing the size of nanoholes to microscale. The experimental results show that the microhole arrays on germanium (Ge) substrate have a main spectral peak of 60% transmittance with a full width at half maximum (FWHM) of ~1.5 µm in addition to some side spectral lobes. The transmission spatial distribution of 30 microhole arrays was also mapped out by a microscopic FTIR spectrometer. Each microhole array shows distinct spectral transmission although some crosstalk is observed for microhole arrays with close physical dimensions.

In 1998, Ebbeson et al. first investigated a periodic subwavelength hole array in a silver film and observed an extraordinary optical transmission (EOT) phenomenon such as enhanced transmission of light through the holes and wavelength filtering due to the excitation of SPR [10]. In case of normal incidence, the peak positions, λmax, of the transmission spectrum of a subwavelength hole array in a hexagonal lattice are approximately given byλmax=P43(i2+ij+j2)εmεdεm+εdwhere P is the period of the array, i and j are the scattering orders of the array, εm and εd are the dielectric constants of the metal and the dielectric material in contact with the metal, respectively[19] .

Following the above equation, multicolor filters in visible spectrum based on nanohole arrays has been extensively investigated in the past decade[20] -[23] . To design a multispectral filter in mid-infrared range, finite difference time domain (FDTD) simulations were performed on an 80 nm thick Au film with periodic micro-sized apertures. The Au film is on a mid-infrared transparent Ge substrate. Since the plasmonic resonances exist at the top and bottom surface of the Au film, a thin Ge film is also covered on top of the Au film to ensure that the plasmonic resonances at both surfaces are at the same wavelength. A light source with wavelength from 2 μm to 15 μm is set polarized in x direction (see SI Fig. S1 for other polarizations) and launched perpendicularly from the thin film top. Periodic MHAs in the Au film are simulated by setting periodic boundary conditions for a unit cell.

According to Eq. (1), the resonance wavelength is mainly determined by periodicity. We first set diameter at a fixed value of 1 µm and shifted the period from 1.5 μm to 2.5 μm at a step of 0.25 µm. It was found that each transmission spectrum has a main peak along with several small side lobes at shorter wavelengths. The main peak is the first order of scattering with (i = 1, j = 0) or (i = 0, j = 1), and the small side lobes are higher orders of scattering. For example, for the period of 1.5 μm, the first order is at ~5.8 μm and the second order (i = 1, j = 1) is at ~3.3 μm (5.8/3). As the period increases, the main peak position shifts from ~5.8 μm to ~9.8 μm but the peak amplitude shows a trend of decline [Fig. 1(a)]. This is not surprising since the relative aperture area becomes smaller for a larger period (the aperture diameter is fixed). As a result, a smaller portion of light will pass through at a given incident light intensity. To investigate the effect of diameter on transmission spectrum, we performed another set of simulations by fixing the period at 2 μm and changing the diameter from 0.8 μm to 1.2 μm at a step of 100 nm. The results are shown in Fig. 1(b), from which we find that the main peak position remains nearly constant, except that the peak amplitude becomes stronger due to the fact that a larger opening at a fixed period will allow more light to transmit.(a) Simulations as periodicity changes from 1.5 μm to 2.5 μm at a fixed diameter of 1 μm. (b) Simulations as diameter changes from 0.8 μm to 1.2 μm at a fixed periodicity of 2 μm. (c) Simulation results of Au MHAs with a fixed ratio 0.5 for aperture diameter to period. (d) Top view of the simulated E-field distribution inside the holes at peak position.

As a multispectral filter for spectral analysis, it is desirable to have a larger transmittance and a narrower FWHM to achieve a higher signal-to-noise ratio (SNR) and spectral resolution. However, as shown in Fig. 1(a,b), a higher transmittance always comes with a wider FWHM. To balance the transmittance and FWHM, we chose a ratio of 0.5 for the aperture diameter to period (see more discussions in SI Section 2). Figure 1(c) shows the simulation results based on this ratio, at which the transmittance maximum is approximately 60% and the FWHM is about 1 µm. As the period increases (diameter also increases accordingly), the main position linearly shifts from 3.5 µm to 13.5 µm and the main peak amplitude slightly increases from 50% to ~65%. For the following experimental work, we fixed the ratio of aperture diameter to period at this value. Figure 1(d) shows that the electric field intensity of the light is mainly localized near the metal edge of the aperture as the light transmits through the aperture due to metal surface plasmon polariton. This is consistent with the previously observed EOT phenomenon at visible spectral range.

To fabricate the multispectral filters according to the simulations above, the electron beam resist ZEP520 was first spin-cast onto the sample and patterned by electron-beam lithography (EBL, Vistec EBPG-5200) at a dosage of 155 μC/cm2. As shown in Fig. 2(a), we then evaporated an 80-nm thick Au film by an electron beam evaporator (Denton) on Ge substrate (1 cm × 1 cm). After that, we immersed the sample in acetone for hours, followed by a short time ultrasonic cleaning to strip the resist together with the gold film on top of it. As the final step, a 160-nm thick Ge protection layer was sputtered on the metal structure using a magnetic control sputtering system (Denton). The created apertures in the metal film are arranged hexagonally to form an array of 50 μm × 50 μm, shown in the scanning electron microscopic (SEM) image in Fig. 2(c). In total, we fabricated 30 microhole arrays with the aperture period tuning from 0.5 μm to 3.4 μm at a step of 100 nm, which are partially shown in Fig. 2(b). The circular aperture diameter is changed accordingly to ensure a constant ratio (1:2) of diameter to period.(a) Fabrication process of Au microhole arrays. Step I: applying EBL to write patterns in ZEP520 on Ge. Step II: depositing an 80 nm-thick Au film via electron beam evaporation. Step III: using acetone to pattern Au film. Step IV: sputter a 160 nm-thick Ge layer on the top of metal. (b) Overview micrograph of Au microhole arrays. (c) A close-up SEM image of the fabricated Au microhole array on Ge wafer.

The transmission spectra of the fabricated Au MHAs on a Ge substrate (“sample spectra”) were measured with a microscopic FTIR spectrometer (Thermo Scientific Nicolet iN10). To remove the effects caused by atmospheric conditions and substrate, the transmission spectra of a Ge substrate (“substrate spectra”) was also recorded. The transmission spectra of the Au MHAs were obtained by dividing the sample spectra with the substrate spectra. We plot six of the total 30 spectra in Fig. 3(a). The peak amplitudes show a trend of decline as the period increases. It is probably because the small aperture size is made slightly larger than designed due to the proximity effect of electron beam exposure. The transmittance peak position as a function of the aperture period is shown in Fig. 3(b). As the period increases from 0.7 μm to 3.2 μm (in 500 nm incremental), the transmittance peak shifts from 3.5 to 13 μm, following almost exactly the same dependence with the simulations. To visualize the wavelength selectivity of the MHA filters, we image the spatial transmittance distribution of the Au film at specific wavelengths (Fig. 3(c,d,e)). The transmittance is coded in color and the image has a spatial resolution of 2 μm × 2 μm. At a wavelength of 5.3 μm, the microhole array with a period of 1.2 μm has the maximum transmittance of 60%, which is consistent with the spectrum (red dash line) in Fig. 3(b). All the other arrays show a transmittance lower than ~40% except for the adjacent two arrays. Similar patterns can be also found in the other two images in Fig. 3(d,e).(a) Experimentally measured optical transmission spectra of six (out of 30 in total) microhole arrays. (b) Transmittance peak position as a function of the aperture period in simulation and experimental results. Mid-infrared transmission 2D maps of the MHAs are plotted at three different wavelengths: 5.3 μm (c), 8.9 μm (d) and 12.9 μm (e).

If each microhole array is equipped underneath with a broadband mid-infrared photodetector, a chip-size spectral analyzer can be realized with appropriate signal processing technology, as shown in Fig. 4(a). In such a system, we can assume a blackbody IR light source that has a spectral radiance density governed by eq. (2).BT(λ)=2hcλ31exp(hc/κTλ)−1where h is the Plank constant, c is the light speed, κ is the Boltzmann constant, T is the absolute temperate and λ is the wavelength. T = 450 k is used here so that the maximum radiation intensity is at λmax = 6.4 μm (Fig. 4(a)). Each array of microholes has a different filtering spectrum as shown in Fig. 3(a). The spectral profile for the ith array is represented by Fi(λ). All the photodetectors have the same wavelength-dependent quantum efficiency η(λ). The transmission spectrum BT(λ)Fi(λ) after the light source is filtered by a microhole array (P = 1.7 μm) is shown as the solid black line in Fig. 4(b). The transmission spectrum can be converted into the photon number per unit area, shown as the red slash profile in Fig. 4(b). Let us suppose that the sample under test has a transmission spectrum X(λ) which is the known. The transmission spectrum BT(λ)Fi(λ)X(λ) after filtered by the same microhole array is shown as the solid black line in Fig. 4(c). Similarly, the transmission spectrum can be converted into the photon number per unit area, shown as the red slash profile in Fig. 4(c).(a) IR spectrometric analyzer model using microhole array multispectral filters. Transmission spectrum of a microhole array (P = 1.7 μm) (b) without and (c) with a test sample, respectively.

The photon number will be collected as photocurrent by photodetectors. For simplicity, we assume that the internal quantum efficiency of the photodetectors is 100%, meaning that all the photons at different wavelengths will be collected by the photodetector and converted to a photocurrent. For the case when the sample under test does not present in the optical path, the photocurrent Ii0 for the ith photodetector will be expressed as Eq. (3) in which hc/λ is the photon energy.Ii0=∫λ1λ2η(λ)BT(λ)Fi(λ)hc/λAcdλwhere Ac is the area of each microhole array filter. For the case when the sample is in the optical path, the photocurrent Ii for the ith photodetector will be written asIi=∫λ1λ2η(λ)BT(λ)Fi(λ)hc/λX(λ)Acdλ

Clearly, the photocurrent ratio Ri=IiIi0 for the ith photodetector will be equal to the transmittance X(λi) at the wavelength λi when the multispectral filter has an ideal performance, that is, Fi(λ) = δ(λi). However, as the spectral profiles of our microhole arrays are not delta functions (Fig. 3(a)), the spectral information of a sample cannot be accurately given by the photocurrent ratios.

For this reason, we employ the following algorithm to reconstruct the spectral information of the sample. Eq. (4) is the photocurrent for the ith photodetector after integrating all wavelengths from 3 um to 14 um (2000 points). For all photodetectors, the photocurrent vector can be written as the following equation when the test sample is not present in the optical path and the quantum efficiency is a unity at all wavelengths.I0=FB,I0∈ℝ30×2000where F∈ℝ30×2000 is the filtering matrix for all 30 filters at the 2000 wavelengths. F is experimentally measured. B∈ℝ2000×1 is the vector of the light source spectrum and the ith element in the vector is equal to BT(λi)∗Δλ/(hc/λi).

Similarly, when the test sample is present in the optical path, the photocurrent vector can be written asI=FXB,I∈ℝ30×1where X∈ℝ2000×2000 is a diagonal matrix and the (i, i) element on the diagonal is the sample spectral information at the i th wavelength. The filtering spectral vector FX is experimental measured when the sample is in the optical path.

However, the matrix X cannot be directly solved from Eq. (6) because F and B are not reversible. Instead of directly solving Eq. (6), we can reformulate Eq. (6) into an optimization problem by finding X to minimize the difference between the two sides of the equation, i.e. min{I−FXB} which is equivalent to min{I−FXB}2=min{(I−FXB)T(I−FXB)}. Minimizing the difference optimized the reconstructed spectrum of the sample under test. Figure 5(a,b) show the FTIR (solid) and reconstructed (dotted) spectrum for CaCO3 and polyethylene (PE), respectively. (We use these two chemicals as examples simply because it is easy to prepare samples since they are solid.) Clearly, our multispectral analyzer can capture the main spectral feature of both materials. Due to the relatively wide FWHM of the filtering spectrum, however, some small side spectral dips are not resolved by the spectral analyzer. These dips are buried in the main valley, resulting in an asymmetric shape of the main valley.(a) FTIR spectrum of CaCO3 (solid) and reconstructed spectrum (dotted); (b) FTIR spectrum of polyethylene (solid) and reconstructed spectrum (dotted).

We should admit that such an IR spectrometric analyzer based on MHAs cannot compete with the Fourier Transform Infrared Spectroscopy (FTIR) in terms of spectral resolution. The measured filtering spectra for the MHA filters has a FWHM of ~1.5 μm on average. This limits the spectral resolution to 0.75 μm (half of the FWHM), meaning that our spectral analyzer will not be able to resolve a spectral feature narrower than 0.75 μm in FWHM (similar to the light not able to resolve an object smaller than half of the wavelength), no matter how many MHA filters are used. But the number of MHA filters, although the more the better, shall not be smaller than a lower limit. Otherwise, the integrated spectral analyzer cannot reach its full capability. To find this lower limit, we can regard the MHA filters as probes to sample the target spectral information. According to the sampling theory, the sampling rate should be at least twice the maximum bandwidth (BW) of the target information, that is,Δk≤12BWwhere Δk is the sampling interval in wavenumber domain. The bandwidth (BW) of the target IR transmission spectrum can be found[24] asBW≈0.35tτwhere tτ is the change of wavenumber as the IR spectrum varies from 10% to 90% of the spectral feature (absorption peak or transmission valley). From Eqs (7) and (8), we findΔk≤tτ0.7

Let us take the IR spectrum of CaCO3 in Fig. 5(a) as an example where tτ≈100cm−1, yielding a maximal sampling interval of 142 cm−1. In the spectral range from 4000 to 400 cm−1, the minimal number of filters is 26. In our case, 30 MHA filters are used, which is large enough to allow our spectral analyzer to reach its spectral resolution, although we designed the MHA filters to sample the target spectrum uniformly in wavelength.

Indeed, our integrated spectral analyzer cannot replace the FTIR spectral analyzer for the general purpose of spectral analysis. But it may be quite useful for some specific applications where the miniaturization of the device is highly desired while the spectral features are not complicated.

In this work, Au MHAs were successfully fabricated on Ge substrate using electron beam lithography and evaporation. The transmission spectra of the microhole arrays were measured with a microscopic FTIR spectrometer. Spatial mappings at fixed wavelengths were applied to visualize the wavelength selectivity of the Au MHA filters. An IR spectral analyzer model based on these filters was constructed to demonstrate the feasibility of retrieving the IR spectra of test materials. The successful development of such an integrated spectral analyzer may find a wide range of applications in non-invasive disease diagnosis, environmental monitoring, food safety and others.

The data that support the findings of this study are available from the authors on reasonable request, see author contributions for specific data sets.

This work is supported by the Key Laboratory of Infrared Imaging Materials and Detectors, Shanghai Institute of Technical Physics and Chinese Academy of Sciences. The microfabrication was conducted at the Center for Advanced Electronic Materials and Devices (AEMD), Shanghai Jiao Tong University. The mid-infrared spectroscopy and microscopic mapping were performed at the Instrumental Analysis Center, Shanghai Jiao Tong University.

Y.D. conceived the idea and directed the research. A.W. performed the experiments. A.W. and Y.D. analyzed the data and wrote the manuscript. All authors reviewed the manuscript.

The authors declare no competing interests.

1 Werle P, Near-and mid-infrared laser-optical sensors for gas analysis, Optics and lasers in engineering, 2002, 37, 101, 114, 10.1016/S0143-8166(01)00092-6

• 2 Klonoff DC, Mid-infrared spectroscopy for noninvasive blood glucose monitoring, IEEE Lasers Electro-optics Soc Newsletters, 1998, 12, 13, 14
• 3 Seddon AB, Mid‐infrared (IR)-A hot topic: The potential for using mid‐IR light for non‐invasive early detection of skin cancer in vivo, physica status solidi (b), 2013, 250, 1020, 1027, 10.1002/pssb.201248524
• 4 Willer U, Saraji M, Khorsandi A, Geiser P, Schade W, Near-and mid-infrared laser monitoring of industrial processes, environment and security applications, Optics and Lasers in Engineering, 2006, 44, 699, 710, 10.1016/j.optlaseng.2005.04.015
• 5 Guillen MD, Cabo N, Infrared spectroscopy in the study of edible oils and fats, Journal of the Science of Food and Agriculture, 1997, 75, 1, 11, 10.1002/(SICI)1097-0010(199709)75:1<1::AID-JSFA842>3.0.CO;2-R
• 6 Eerdenbrugh B, Taylor LS, Application of mid-IR spectroscopy for the characterization of pharmaceutical systems, International journal of pharmaceutics, 2011, 417, 3, 16, 10.1016/j.ijpharm.2010.12.011
• 7 Baker MJ, Using Fourier transform IR spectroscopy to analyze biological materials, Nature protocols, 2014, 9, 1771, 1791, 10.1038/nprot.2014.1104480339
• 8 Piegari A, Bulir J, Variable narrowband transmission filters with a wide rejection band for spectrometry, Applied optics, 2006, 45, 3768, 3773, 10.1364/AO.45.003768
• 9 Mulvaney P, Surface plasmon spectroscopy of nanosized metal particles, Langmuir, 1996, 12, 788, 800, 10.1021/la9502711
• 10 Ebbesen TW, Lezec HJ, Ghaemi H, Thio T, Wolff P, Extraordinary optical transmission through sub-wavelength hole arrays, Nature, 1998, 391, 667, 669, 10.1038/35570
• 11 Yokogawa S, Burgos SP, Atwater HA, Plasmonic color filters for CMOS image sensor applications, Nano Letters, 2012, 12, 4349, 4354, 10.1021/nl302110z
• 12 Burgos SP, Yokogawa S, Atwater HA, Color imaging via nearest neighbor hole coupling in plasmonic color filters integrated onto a complementary metal-oxide semiconductor image sensor, Acs Nano, 2013, 7, 10038, 10047, 10.1021/nn403991d
• 13 Chang C-Y, Wavelength selective quantum dot infrared photodetector with periodic metal hole arrays, Applied Physics Letters, 2007, 91, 163107, 10.1063/1.2800378
• 14 Jang W-Y, Experimental Demonstration of Adaptive Infrared Multispectral Imaging using Plasmonic Filter Array, Scientific reports, 2016, 6, 10.1038/srep348765056518
• 15 Cao H, Nahata A, Resonantly enhanced transmission of terahertz radiation through a periodic array of subwavelength apertures, Optics express, 2004, 12, 1004, 1010, 10.1364/OPEX.12.001004
• 16 McCrindle IJH, Grant JP, Gouveia LCP, Cumming DRS, Infrared plasmonic filters integrated with an optical and terahertz multi‐spectral material, physica status solidi (a), 2015, 212, 1625, 1633, 10.1002/pssa.201431943
• 17 Chang C-C, A surface plasmon enhanced infrared photodetector based on InAs quantum dots, Nano letters, 2010, 10, 1704, 1709, 10.1021/nl100081j
• 18 Brolo AG, Gordon R, Leathem B, Kavanagh KL, Surface plasmon sensor based on the enhanced light transmission through arrays of nanoholes in gold films, Langmuir, 2004, 20, 4813, 4815, 10.1021/la0493621
• 19 Genet C, Ebbesen T, Light in tiny holes, Nature, 2007, 445, 39, 46, 10.1038/nature05350
• 20 Chen Q, CMOS photodetectors integrated with plasmonic color filters, IEEE Photonics Technology Letters, 2012, 24, 197, 199, 10.1109/LPT.2011.2176333
• 21 Si G, Zhao Y, Chew AB, Liu YJ, Plasmonic color filters, Journal of Molecular and Engineering Materials, 2014, 2, 1440009, 10.1142/S2251237314400097
• 22 Tan SJ, Plasmonic color palettes for photorealistic printing with aluminum nanostructures, Nano letters, 2014, 14, 4023, 4029, 10.1021/nl501460x
• 23 Li, Z., Clark, A. W. & Cooper, J. M. Dual color plasmonic pixels create a polarization controlled nano color palette (2016).
• 24 Orwiler, B., Oscilloscope vertical amplifiers. (Tetronix, Inc., 1969).

Supplementary information accompanies this paper at 10.1038/s41598-018-29177-0.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR)

PHOTO (COLOR): supplementary information