Differential Near-Field Scanning Optical Microscopy with THz quantum cascade laser sources R. Degl’Innocenti*1, , M. Montinaro1 , J. Xu1 , V. Piazza1 , P. Pingue1 , A. Tredicucci1 and F. Beltram1 H. E. Beere2 , D. A. Ritchie2 1 NEST,
CNR-INFM and Scuola Normale Superiore, Piazza San Silvestro 12, I-56127 Pisa, Italy 2 Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom ∗
[email protected]
Abstract: We have realized a differential Near-field Scanning Optical Microscope (NSOM) working with subwavelength resolution in the THz spectral region. The system employs a quantum cascade laser emitting at λ ∼ 105 μ m as source, and the method, differently from conventional NSOM, involves diffracting apertures with size comparable to the wavelength. This concept ensures a higher signal-to-noise level at the expense of an additional computational step. In the implementation here reported λ /10 resolution has been achieved; present limiting factors are investigated through finite difference time domain simulations. © 2009 Optical Society of America OCIS codes: (110.6795) Terahertz imaging; (180.4243) Near-field microscopy; (140.5965) Semiconductor laser, quantum cascade.
References and links 1. S. Barbieri, J. Alton, C. Baker, T. Lo, H. E. Beere, D. Ritchie, “Imaging with THz quantum cascade lasers using a Schottky diode mixer,” Opt. Express 13, 6497-6503 (2005). 2. J. Darmo, V. Tamosiunas, G. Fasching, J. Kr¨oll, K. Unterrainer, M. Beck, M. Giovannini, and J. Faist, C. Kremser, P. Debbage,“Imaging with a Terahertz quantum cascade laser,” Opt. Express 12, 1879-1884 (2004). 3. B.B. Hu, M.C. Nuss, “Imaging with terahertz wave,” Opt. Lett. 20, 1716-1718 (1995). 4. Y. C. Shen, T. Lo, P. F. Taday, B. E. Cole, W. R. Tribe, and M. C. Kemp, “Detection and identification of explosives using THz pulsed spectroscopic imaging,” Appl. Phys. Lett. 86, 241116-241118 (2005). 5. N. C. J. van der Valk and P. C. M. Planken, “Electro-optic detection of subwavelength terahertz spot sizes in the near field of a metal tip,” Appl. Phys. Lett. 81 1558-1561 (2002). 6. F. Buersgens, R. Kersting, H-T. Chen, “Terahertz microscopy of charge carriers in semiconductors,” Appl. Phys. Lett. 88, 112115-112118 (2006). 7. A. J. Huber, F. Keilmann, J. Wittborn, J. Aizpurua and R. Hillenbrand, “Terahertz Near-Field Nanoscopy of Mobile Carriers in Single Semiconductor Nanodevices,” Nano Lett. 8, 3766-3770 (2008). 8. J. R. Knab, A. J. L. Adam, M. Nagel, E. Shaner, M. A. Seo, D. S. Kim, and P. C. M. Planken, “Terahertz Near-Field Vectorial Imaging of Subwavelength Apertures and Aperture Arrays,” Opt. Express 17, 15072-15086 (2009). 9. A. J. L. Adam, J. M. Brok, M. A. Seo, K, J. Ahn, D. S. Kim, J. H. Kang, Q. H. Park, M. Nagel, P. C. M. Planken, “Advanced terahertz electric near-field measurements at sub-wavelength diameter metallic apertures,” Opt. Express 16, 7407-7417 (2008). 10. O. Mitrofanov, M. Lee, J. W. P. Hsu, L. N. Pfeiffer, K. W. West, and J. D. Wynn, J. F. Federici, “Terahertz pulse propagation through small apertures, ” Appl. Phys. Lett. 79, 907-909 (2001). 11. O. Mitrofanov, M. Lee, J. W. P. Hsu, I. Brener, R. Harel, J. F. Federici, J. D. Wynn, L. N. Pfeiffer, and K. W. West, “Collection-mode Near-Field Imaging With 0.5-THz Pulses,” IEEE J. Sel. Top. Quantum Electron. 7, 600-607 (2001).
#118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23785
12. R. K¨ohler, A. Tredicucci, F. Beltram, H. E. Beere, E. H. Linfield, A. G. Davies, D. A. Ritchie, R. Iotti, and F. Rossi, “THz semiconductor-heterostructure laser,” Nature 417, 156-159 (2002). 13. B. S. Williams, “Terahertz quantum-cascade lasers,” Nature Photon. 1, 517-525 (2007). 14. A. Ozcan, E. Cubucku, A. Bilenca, K. B. Crozier, B. E. Bouma, F. Capasso, and G. J. Tearney, “Differential Near-Field Scanning Optical Microscopy,” NanoLett. 6, 2609-2616 (2006). 15. A. Ozcan, E. Cubukcu, A. Bilenca, B. E. Bouma, F. Capasso, G. J. Tearney, “Differential Near-Field Scanning optical Microscopy Using Sensor Arrays,” IEEE J. Sel. Top. Quantum Electron. 13, 1721-1729 (2007). 16. H. A. Bethe, “Theory of diffraction by Small Holes,” Phys. Rev. 66, 163-182 (1944). 17. L. Novotny and B. Hecht,“Principle of Nano-Optics” Cambridge University Press, U. K. (2007). 18. T. Losco, J. H. Xu, R. P. Green, A. Tredicucci, H. E. Beere, and D. A. Ritchie, “THz quantum cascade designs for optimized injection,” Physica E 40, 2207-2209 (2008). 19. Kawano Y., Ishibashi K.,“An on chip near-field terahertz probe and detector,” Nature Photon. 2, 618-621, 2008.
1.
Introduction
THz technology has attracted a great interest in many different applications, ranging from biomedical imaging [1,2], to semiconductor device inspection [3] and security screening [4]. In particular, spectroscopy in the THz region can represent a non-destructive method and a valuable tool for molecular sensing or, alternatively, for retrieving information about the dielectric properties of the material. Near-Field Scanning Optical Microscopy (NSOM) in the THz range (corresponding to a wavelength λ between 30 μ m and 300 μ m) constitues a relevant technique to achieve imaging with sub-wavelength resolution, which is crucial for most of these applications, demanding a resolution of few micrometers or better. Although there are well known limits connected with conventional aperture-less NSOM, (collection of the total scattered field from the probe-sample system, the size of the scanning probe, etc . . .), near field THz imaging has already been successfully demonstrated with submicrometer resolution [5,6] using THz time-domain spectroscopy set-ups. Recently, by combining the strong concentration of the THz radiation on a metal tip with interferometric detection, imaging of the carriers in a semiconductor device has been achieved using a continuous-wave CH3 OH gas laser with the impressive resolution of λ /3000 at a wavelength λ =118 μ m [7]. On the other hand, THz NSOM has also been realized with sub-wavelength circular [8, 9] or square metallic aperture probes [10, 11], and a spatial resolution of 7 μ m was obtained over the broad spectral range of λ = 120-1500 μ m. Quantum cascade (QC) lasers have now become a very promising source of THz radiation, thanks to their high power output, spectral purity, compactness and reliability [12, 13]. A subwavelength imaging instrument based on QC lasers would permit to take advantage of the higher stability of these devices compared for example to gas lasers, or to the fs-pulse laser systems commonly used to generate THz radiation, reduce the power consumption of the whole system, and finally allow to perform the analysis with minimal changes at any chosen frequency in the 1-5 THz spectral range. Lately, a new microscopy approach has been proposed named differential NSOM [14, 15]. This concept differs from conventional NSOM because the acquisition is carried out using square apertures with dimensions comparable to the wavelength, which ensure a much higher signal/noise ratio. Here we adopted this new technique and developed a subwavelength resolution ( λ /10) differential NSOM with a QC laser source operating at a wavelength of 105 μ m. The paper is structured as follows. In the next section, we describe the experimental set-up and the main features of our system. In the third section we discuss briefly the basical theoretical background and present the numerical simulations carried out to validate our approach and to analyze polarization effects. In the last section we then present measurements realized with a CO2 gas laser emitting at a frequency 30 THz to study the dependence on sample-aperture distance, and the measurements performed at a frequency of 3 THz with the QC laser, comparing them #118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23786
with theoretical predictions. Finally, in the conclusions, we comment on the results and discuss possible future improvements. 2.
Experimental set-up
In standard aperture NSOM the light is diffracted by an aperture with size much smaller than the wavelength, which limits drastically the signal-to-noise ratio achievable, since in first approximation the electric field amplitude, in the limit of apertures much smaller than the wavelength, scales as the 3rd power of the aperture size [9, 10, 16]. In our procedure, we follow the theoretical approach described in [14, 15], which involves diffracting square apertures with size comparable with the irradiating wavelength. The aperture, illuminated by a laser source, is scanned in the near field of the object of interest, which is mounted on a xy translational stage, perpendicular to the incident radiation, driven with piezo-actuators, as shown in Fig. 1. The total transmitted power is then recorded in the far-field for each scanning position by a suitable detector. Differently from NSOM, recovery of the sample image is accomplished by performing a numerical second order derivative ∂ 2 /∂ x∂ y of the image acquired by the detector. The final result consists of four independent replica of the original object, two positive and two negative, originating from the diffraction at the aperture’s corners. In our measurements we have employed a THz QC laser emitting around 3 THz, operating at liquid Helium temperature. We fabricated, with standard photolithography, square apertures of different sizes varying from 70×70 μ m2 to 150×150 μ m2 in a Cr/Au metalization on a GaAs substrate, with typical thickness of 10/500 nm. The objects realized to test the system correspond to triangular or square metal shapes evaporated on quartz, with typical lateral size of 10→30 μ m, and approximately the same thickness of the aperture, sufficient to screen completely the laser wavelength. The objects are placed carefully in the near field of the aperture using a charge-coupled device (CCD) camera, which allows to determine the relative distance with an indetermination of Δz 3 μ m, and a z-piezo translation stage. In the inset of Fig. 1 we show a top-view visible picture of the sample, a metal triangle evaporated on a quartz substrate, placed in the near-field of the square aperture (100 μ m side). The THz signal passing through the aperture/sample system and collected by an off-axis parabolic mirror with numerical aperture 0.7 is focused on a Pyrodetector 1000 from Infrared System and then recorded as a function of the sample-aperture relative position. In order to ensure the maximum signal/noise ratio, we have implemented a lock-in technique for the far-field detection, and we have placed the whole set-up in a closed environment filled with Nitrogen gas, to reduce water absorption. 3.
Numerical simulations
The precision of these measurements in principle depends only on the sharpness of the corners, on the minimum step size, and on the relative distance Δz between object and aperture. In fact, the nonpropagating evanescent waves with high spatial frequency in the xy plane, originating at the corners of the aperture and carrying the subwavelength information, have an exponential decay in the z-direction [17]. The interaction of these high spatial frequencies with the object shifts the evanescent waves into oscillating waves that can be collected by the optical system in the far-field. In order to investigate the limits of our approach, the resolution, and the effect of the polarized QC laser emission, we have performed accurate simulations with the Finite Difference Time Domain (FDTD) method using the commercial software Remcom XFdtd, version 6.5 to calculate the intensiy profile emitted from the aperture, while we have performed the convolution between aperture and samples, and the bi-dimensional derivative process, with a custom C ANSI code. In our simulations, we assumed a linearly polarized plane wave incident on a 200 μ m thick GaAs substrate, and a square aperture with typical size of 100×100 μ m2 opened #118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23787
Fig. 1. Scheme of the experimental apparatus. The radiation of a cryogenically cooled QCL laser emitting at λ ∼ 105 μ m, is collected by a system of parabolic mirrors (not shown here) and then focused through a Winston cone on a square aperture in a metal evaporated on a GaAs substrate. The sample corresponds to a metallic layer with regular geometrical shape, a triangle in this case, on a quartz substrate, which is transparent in the THz range. The sample, mounted on a xyz piezoelectric traslational stage, is placed in the near-field of the aperture with the help of a CCD camera.
in a 400 nm thick Au layer evaporated on top of the GaAs substrate. The radiation frequency was kept fixed at 100 μ m in all the simulations. In Fig. 2(a) we present the total calculated intensity transmitted throught the aperture when the polarization is along a side of the aperture and when it is tilted to 45◦ with respect to the aperture side. At the interface between a metal and a dielectric material, boundary conditions are strongly dependent on the polarization, and surface plasmons are also preferentially excited accordingly. These effects, that critically affect the intensity distribution in the aperture, as showed in Fig. 2(a) left, are partly compensated by tilting the polarization to 45◦ , so that all four sides of the square aperture are now almost equivalent, as presented in Fig. 2(a) right. Clearly, also the convolution of the intensity transmitted through the aperture with the sample (a 30 μ m side totally opaque triangle similar to the one displayed in the inset of Fig. 1, placed at a distance Δz = 4 μ m from the aperture) presents a similar asymmetric response, as shown in Fig. 2(b). In turn, also the final recovery of the sample is affected, and the four images appear strongly altered in the parallel configuration, as seen in Fig. 2(c). 4.
Operation and measurements
In light of the simulation results described above, all our measurements have been carried out with the laser beam polarization positioned at 45◦ with respect to the aperture side. The minimum relative distance Δz = z0 attainable between sample and aperture has been estimated to be 3-4 μ m. Once the total collected signal passing through the aperture/sample system has been focused on the detector, a scan of the sample is performed in the xy plane, with an extension in the xy directions equal to twice the aperture size to avoid spatial aliasing [14], and the signal is recorded as a function of the aperture position. First measurements have been carried out with a CO2 gas laser emitting around λ 10.6 μ m in order to test the z-dependence of the system, which is obviously more sensitive at short wavelengths. We have recorded the same scan at different Δz between the sample, a 30 μ m side triangular metal shape on quartz shown in Fig. 1, #118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23788
Fig. 2. Numerical simulations obtained with the FDTD method. The aperture has 100×100 μ m2 size, and is defined in a 400 nm thick metal layer on top of a GaAs substrate illuminated with linearly polarized radiation with a wavelength of 100 μ m. The object is a metallic triangular shape completely opaque to the THz radiation, and the relative distance Δz has been kept fixed to 4 μ m. If the polarization is along one axis of the square aperture, the intensity profile of the aperture transmission, shown in a), is not uniform, due to the boundary conditions at the metal/dielectric interface related to the polarization and to the excitation of surface plasmons. This, in turns, affects the convolution with the sample, as it is presented in b), where the integrated transmitted intensity in the presence of the sample is displayed as function of the aperture-sample displacement, and the final recovery of the original object reported in c).
and a 50×50 μ m2 aperture in a 5/400 nm thick Cr/Au layer evaporated on a GaAs substrate. The total signal measured is presented in Fig. 3(a) for each Δz, and the corresponding recovery of the original image arising from the upper corner of the aperture, is reported in Fig. 3(b). While the total collected signal registered exhibits minimal changes at different Δz, the image obtained after the derivative process is more sensitive to this parameter. In fact, the original object can be distinguished only for the minimum relative distance Δz = z0 . All the measurements in the THz frequency range have been carried out using square apertures as large as 100×100 μ m2 , and a QC laser source emitting at λ 105 μ m [18]. As an example, we consider in Fig. 4 a sample composed of a square metallic layer with 28 μ m long sides, placed in the near #118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23789
Fig. 3. a) Total signal acquired using a CO2 laser source emitting at λ 10.6 μ m wavelength for different relative distance Δz between sample and aperture; b) Recovered image of the original sample (detail originating from one corner). The aperture, with size 50×50 μ m2 , is defined in a 5/400 nm Cr/Au evaporation on a GaAs substrate, and the sample corresponds to the one shown in Fig. 1. The minimum relative distance z0 between sample and aperture is estimated to be 3-4 μ m.
Fig. 4. The visible image of a square metallic evaporation with 28×28 μ m2 size on a quartz substrate, placed in the near field of the aperture, is shown in a). The result of the theoretical simulations, carried out assuming a relative distance object/aperture Δ z = 4 μ m, an incident wavelength of 100 μ m, a square opaque object of 30 μ m and an aperture of 100 μ m, are presented in b). The final recovery of the original sample, obtained from the measured THz signal after the second derivative process, is presented in c).
field of the aperture. Numerical simulations were carried out for comparison keeping the relative distance Δz fixed to 4 μ m between the 100×100 μ m2 aperture and a completely opaque #118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23790
square with 30 μ m side, representing the sample. According to theory, after having performed a numerical second order derivative ∂ 2 /∂ x∂ y of the total collected signal, the four replicas of the original object, two positive and two negative, appear, as reported in Fig. 4(b). The data obtained from the measurements are shown in Fig. 4(c) and show a good agreement with the calculated images. The final recovery in 4(c), has been obtained averaging over 3-4 pixels (corresponding to 8-10 μ m), in order to increase the signal/noise ratio, and using sharpening edge filters.
Fig. 5. a) The THz image of a 30-μ m square as recorded without further treatment from a single corner replica; b) the same for the object already presented in Fig. 4(a), but summing up the four replica originating from each corner of the aperture. For comparison a single simulated image is presented in c) selected among those shown in Fig. 4(b). The dashed lines indicate the chosen profiles extracted in the pictures below, yielding a resolution, respectively, of 16 ± 2 μ m in a), of 11 ± 2 μ m in b), and of 9 ± 1 μ m in c) (10 %- 90 % criterion).
In order to derive an estimate of the system resolution, we present instead in Fig. 5 untreated individual images. Panel a) displays the THz image of a 30-μ m side square as recorded from a single aperture corner. Alternatively, we report in panel b) the image of a 28-μ m side square obtained by summing the four independent replicas, appropriately cropped and considered with the right sign, originating from each aperture corner. The second procedure increases the signal/noise without altering the resolution of the image, provided that the aperture plane is parallel to the sample and that the four corners are identical. Figure 5(c) represents one of the replica selected from the simulation presented in Fig. 4. The intensity profiles extracted along the dashed lines are plotted below each corresponding image. The measured profiles are consistent with the simulated one, although there is no evident flat maximum plateau, which may again be the result of the coupling to the surface plasmons on the metal surfaces of aperture and object. This of course makes it difficult to define a proper resolution value; in any case, assuming conservatively the conventional 10%-90% criterion, we obtained a resolution of 16 ± 2 μ m for image a) and 11 ± 2 μ m for image b), with the theoretical case of panel c) only slightly better at 9 ± 1 μ m. The experimental values tend to vary over a few μ m range as the #118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23791
resolution depends quite sensitively on alignment, distance, illumination uniformity, etc. The reasonable agreement between the values of panel b) and c) is a good indication that the estimate of a distance Δz between aperture and sample of 4 μ m is indeed fairly correct. The resolution of 10 μ m, comparable to what was achieved by other aperture-based techniques [10,19], does not represent the final limit of this approach. According to our simulations, a relative distance sample-aperture of 1 μ m, still feasible without a closed servo-loop feedback system, would yield a resolution of 2 μ m. In our present experiment, the resolution is limited not only by the distance between sample and aperture, but also from the indetermination in the direction of the polarization with respect to the aperture position, from the aperture/sample materials, from the intensity distribution impinging on the aperture, and from the signal/noise ratio. An optimization of these factors, improving the aperture/sample alignment and using for instance a cryogenically cooled bolometer instead of the pyro-electric detector, is necessary for further enhancements. Due to the versatility of our approach, the system can operate in a wide frequency range, just apporting minimal changes to the set-up, thus profiting of the whole operating frequency range of QC lasers. Furthermore, the achieved resolution is already enough for investigating biological samples, or DNA chips, allowing a non-destructive spectroscopic analysis. 5.
Conclusions
In conclusion, we have experimentally realized a Differential Near-Field Scanning Optical Microscope, working in the THz spectral region with a resolution of ≈ λ /10, and relying on QC lasers as light sources. The results, which have been obtained using a laser emitting at a wavelength of λ 105 μ m, can be easily extended to a wide frequency range and have been validated through FDTD numerical simulations. Future developments will aim at improving the system resolution, and the results should further stimulate the investigation of biological samples. Acknowledgements This work was supported in part by the European Commission through the integrated project ”Teranova”. We also acknowledge support from the Italian Ministry of Research through the project ”National Laboratory for Nanotechnology applied to Genomics and Post-Genomics”.
#118476 - $15.00 USD Received 16 Oct 2009; revised 21 Nov 2009; accepted 29 Nov 2009; published 11 Dec 2009
(C) 2009 OSA
21 December 2009 / Vol. 17, No. 26 / OPTICS EXPRESS 23792