Boosting local field enhancement by on-chip nanofocusing and impedance-matched plasmonic antennas
Strongly confined surface plasmon-polariton modes can be used for efficiently delivering the electromagnetic energy to nano-sized volumes by reducing the cross sections of propagating modes far beyond the diffraction limit, i.e., by nanofocusing. This process results in significant local-field enhancement that can advantageously be exploited in modern optical nanotechnologies, including signal processing, biochemical sensing, imaging and spectroscopy. Here, we propose, analyze, and experimentally demonstrate on-chip nanofocusing followed by impedance-matched nanowire antenna excitation in the end-fire geometry at telecom wavelengths. Numerical and experimental evidences of the efficient excitation of dipole and quadrupole (dark) antenna modes are provided, revealing underlying physical mechanisms and analogies with the operation of plane-wave Fabry-Pérot interferometers. The unique combination of efficient nanofocusing and nanoantenna resonant excitation realized in our experiments offers a major boost to the field intensity enhancement up to , with the enhanced field being evenly distributed over the gap volume of , and promises thereby a variety of useful on-chip functionalities within sensing, nonlinear spectroscopy and signal processing. [This document is the unedited Author’s version of a Submitted Work that was subsequently accepted for publication in Nano Letters, ©American Chemical Society after peer review. To access the final edited and published work see http://dx.doi.org/10.1021/acs.nanolett.5b03593.]
The major aspect of focusing of electromagnetic radiation is the possibility of concentrating the energy in a small volume. Because of diffraction, the focusing of free-propagating optical waves is limited in size to the half of the light wavelength in the medium â the diffraction limit of light. (1) One approach to overcome this limit is to use surface plasmon polaritons (SPPs) â surface electromagnetic modes bound to and propagating along metal-dielectric interfaces, with electromagnetic fields in a dielectric being coupled to collective free electron oscillations in a metal. (2) Spatial confinement of SPP modes in the cross section perpendicular to the propagation direction depends on the material composition and geometric configuration of a waveguiding structure. Notably, some SPP modes (supported, for example, by metal nanowires (3)) exhibit a unique scaling property in their spatial confinement: the mode is progressively better confined for smaller lateral waveguide dimensions, opening thereby the possibility for guiding extremely confined (i.e., on a deep subwavelength scale) SPP modes(4) as well as for designing SPP-based nanoantennas. (5) This feature can further be used for nanofocusing, which is the process of reducing the cross sections of propagating optical modes far beyond the diffraction limit, simply by gradually decreasing lateral waveguide dimensions. (6)
Nanofocusing using SPP modes was first theoretically described for planar tapered metal-insulator-metal (MIM) and insulator-metal-insulator (IMI) geometries(7) and for tapered nanorods, (8); (9) instigating experiments with various configurations demonstrating this fascinating phenomenon (see review by Gramotnev and Bozhevolnyi(6) and references therein). As far as on-chip nanofocusing is concerned, the experimental demonstrations using various SPP configurations and material platforms have been reported. (10); (11); (12); (13); (14); (15); (16); (17); (18) Arguably the most impressive results were obtained with a MIM waveguide followed by a three-dimensional linear taper that facilitated the nanofocusing of infrared radiation, producing the intensity enhancement of 400 within a subwavelength area. (12) It should however be noted that the considered structure requires the use of rather sophisticated fabrication techniques as well as poses considerable difficulties in accessing the enhanced field, which is mainly distributed within the insulator between the metal stripes. Other configurations (which are less demanding in fabrication but also designed for considerably longer wavelengths) include the nanofocusing using tapered transmission lines, where the moderate intensity enhancement of was demonstrated within a hot spot of nm diameter, which is equivalent to , where is the free space wavelength. (17); (18) Finally, the nanofocusing has been demonstrated with even simpler (IMI-based) configurations, such as SPP stripe waveguides, in which linearly tapered stripes were used to concentrate the electromagnetic energy carried by short-range SPP modes(5) far beyond the diffraction limit with a relatively large efficiency of .(14); (15); (16)
In this paper, we propose, analyze, and experimentally demonstrate efficient on-chip nanofocusing at telecom wavelengths using a tapered-stripe nanofocuser (TSNF) coupled (in the end-fire excitation geometry) to an impedance-matched nanowire antenna across a gap. The nanofocusing effect in the TSNF is thereby combined with and enhanced by the resonance in the nanoantenna, which can be considered as an SPP nanowire analogue of the Fabry-Pérot interferometer, allowing us to further significantly boost the local intensity enhancement up to about 12,000, with the enhanced field being evenly distributed over the gap volume (). The taper shape is optimized numerically by maximizing its intensity enhancement and transmittance. The latter is additionally measured experimentally using the amplitude- and phase-resolved scanning near-field optical microscopy (SNOM). Finally, the experimentally studied antenna-coupled TSNF (a-TSNF) provides the evidence of resonant impedance-matched excitation of the nanoantenna by featuring a considerable local-field enhancement (FE) and substantial resonant decrease of back reflection from the tip of the taper, which is similar to the reflection suppression from a Fabry-Pérot interferometer under the resonance conditions.
The design of our nanofocusing structures (Figure 1a,b) is based on a 30-nm-thick gold stripe waveguide. The sample was fabricated using a combination of the electron-beam lithography and lift-off technique on a silica substrate (see Supporting Information for details). The width of the widest part of the taper () was chosen to be close to the free-space wavelength of the illumination ( nm). The width of the narrowest part of the taper ( nm) is equal to the metal thickness in order to avoid fabrication difficulties. The metal thickness was chosen as a trade-off between the field localization and Ohmic losses. It should be noted that, for the chosen range of the metal stripe width, only one guided mode is supported, which is convenient for observation of a nanofocusing effect. The SPP stripe mode was excited by the laser beam normally incident on the grating (Figure 1a) designed to have the period equal to the wavelength of the guided plasmonic mode . In order to suppress the effects of any leaky or radiating modes, excited by diffraction at the grating, a 5--long stripe was placed between the grating and the taper entrance.
The structures were experimentally investigated using a scattering-type SNOM, based on an atomic force microscope (AFM). We used our SNOM in the transmission mode(19); (20) with the sample illuminated normally from below (Figure 1c). SNOM with a platinum-coated AFM Si probe was operating in a tapping mode at frequency kHz, while the sample was moving during the scan. The lower parabolic mirror was moving synchronously with the sample in the -plane in order to maintain the excitation alignment during the scan. The signal, scattered by the probe, was detected and demodulated at the third harmonic () to filter the near-field contribution from the background. Additionally, our SNOM uses interferometric pseudoheterodyne detection, (21) which allows imaging of both the amplitude and the phase of the near-field. Finally, it should be mentioned that because of the AFM probe shape, elongated along the out-of-plane -axis, the measured optical signal corresponds mostly to the component of the near-field. (20); (22); (23)
It is known that there are certain constraints on the taper shape for adiabatic nanofocusing. (6); (9) However, the Ohmic losses of plasmonic modes increase with mode compression. Therefore, in order to maximize the efficiency of the taper, one should balance between Ohmic and scattering losses due to non-adiabaticity. We performed numerical optimization of the taper shape, fixing the input and output widths at and nm, correspondingly. First, the mode analysis for the two stripe widths and was performed, using the finite element method (Figure 2). The mode profile in the wide stripe (Figure 2a) demonstrates that the electric field is mostly concentrated in the substrate below the stripe, while on the air side it is noticeable only near the edges. For the narrow stripe (called nanowire below), the mode profile in Figure 2b demonstrates nearly uniform distribution of the electric field around the waveguide with âradialâ polarization.(3) Thus the nanofocusing of plasmons along the tapered stripe is topologically the same as nanofocusing along three-dimensional metal cones. (16) In order to quantify the nanofocusing effect correctly without influence of the lightning-rod effect near the sharp edges, (24) we set the probe point (blue circle in Figure 2a) just below the wide stripe in the silica substrate. For the nanowire (Figure 2b), the probe point was placed just above the nanowire in the air, where it is accessible for applications. From the middle cross-sections (Figure 2c) one can see that, for the same total transmitted power, the electric-field amplitude in the nanowire mode (the dashed curve in Figure 2c) is 18.5 times larger than the field at the wide section of the stripe (the solid curve in Figure 2c), meaning the intensity enhancement of (under the assumption of zero losses). Interestingly, this intensity enhancement was achieved upon just the 50-fold reduction of the taper width from 1500 nm to 30 nm. The mismatch between the reduction of the taper width and intensity enhancement () can be explained by the two-dimensional mode compression in the waveguide tapered in one dimension (compare Figures 2a and 2b).
A numerical optimization of the taper shape by maximizing its intensity enhancement and transmittance (see Supporting Information, Figures S1-3) resulted in a combined two-section taper (Figure 1a,b), whose transmittance was . Here, transmittance is defined as the ratio of the energy flow along the -axis in the guided SPP mode at the taper width of 30 nm to that at the beginning of the taper (Figure 1a,b). Interestingly, the losses of were distributed almost equally between Ohmic losses in the metal () and scattering into the substrate (31%), while only of the energy was reflected back. Thus, the expected intensity enhancement of a lossy TSNF is . Note that the optimum taper shape (including the angles of each section) depends on the stripe thickness, free-space wavelength, etc. Therefore, thorough numerical simulation should be performed for each particular case to find its own optimum shape.
To experimentally characterize the TSNF transmittance, we fabricated structures with two TSNFs, connected with a nanowire (Figure 3). The length of the middle nanowire was 1.5 in one case (Figure 3a) and zero in the other (Figure 3b). As can be seen from the near-field amplitude images, transmission through the two TSNFs connected by the nanowires is relatively high. The complex near-field pattern of each section with the wide stripe is due to interference of the guided, leaky and radiation modes. Similarly to our previous work, (19) the interference pattern was carefully analyzed and decomposed into individual modes by numerical fitting procedure (see Supporting Information, Figures S13-15). At the same time, the domination of the main guided mode (direct and reflected) is demonstrated by the bright fringes near the edges of the wide stripe, which agrees well with its mode profile (Figure 2a). Assuming the in-coupling and out-coupling transmittances of the TSNF are the same (supported by numerical simulations â see Supporting Information, Figure S3), the transmittance for a single TSNF is , which is in good agreement with the numerical simulations.
Once an SPP stripe mode has been nanofocused into a strongly confined guided mode of a nanowire, a few additional options can be used to further enhance the local electric field. Simple back-reflection of the focused SPP mode can cause additional 2-fold FE ( 4-fold for the intensity enhancement) because of constructive interference. (12) Another option is to introduce a small gap in the metal waveguide. In this case, the longitudinal component of the electric field inside metal can be enhanced in the gap (filled with dielectric) because of the boundary conditions by a factor equal to the ratio of the dielectric constants of the metal and dielectric. Considering a gold waveguide ( at = 1425 nm(25)) with an air gap (), one would expect 80-fold increase of the longitudinal electric field component inside the gap. Since for a nanowire the longitudinal component of the electric field inside the metal is 3.8 times smaller than the total electric field in the air on the nanowire surface (Figure 2b), the gap-induced FE is (440 for intensity enhancement). However, our numerical simulations (see Supporting Information, Figure S4) show that, to achieve these effects, the gap should be sufficiently small â for example, a 10 nm gap results in 60% reflectance and FE of only 5.8.
The gap-induced FE can be combined with nanofocusing by introducing a gap at the tip of the taper, forming a gapped TSNF (g-TSNF, Fig. 4a). Its intensity enhancement (solid blue line) can be calculated as a product of TSNF intensity enhancement and a square of the gap-induced FE, which agrees well with the simulations of the whole structure (blue dots). Thus, a 10 nm gap results in g-TSNF intensity enhancement of . Similarly, the reflectance (red) can be estimated as a product of the gap-induced reflectance in nanowire and a square of TSNF transmittance.
Finally, even greater FE could be achieved when using a dipole nanoantenna coupled to the tip of a TSNF across a narrow nanogap that plays the role of a semi-transparent mirror in a conventional Fabry-Pérot interferometer. In this case of an antenna-coupled TSNF (a-TSNF), the nanofocusing effect is combined with the antenna resonance, producing particularly strong local fields in the gap between the tip of the TSNF and the nanoantenna. This is also the geometry of the end-fire excitation of a dipole nanoantenna by means of a nanofocusing structure in the on-chip configuration. Expectedly, the maximum of the field inside the gap approximately corresponds to the minimum of the back reflection (Figure 4b) â small relative extremum displacements can be explained by the Ohmic and radiative losses. (26); (27) Note, however, that the indicated analogy with the Fabry-Pérot interferometer is not complete. For example, one of the significant physical differences is the strong enhancement of the local field in the nanogap compared to the rest of the structure, which is not typical for a conventional Fabry-Pérot resonator with mirrors and bulk waves.
The resonant excitation of a 240-nm-long nanoantenna leads to antenna-induced FE of 9.3 (see Supporting Information, Figure S5), which corresponds to the overall a-TSNF intensity enhancement of . This is 10 times larger than for the terminated TSNF (i.e., at zero antenna length in Figure 4b) and 2.5 times larger than for the g-TSNF (i.e., with an infinitely long nanoantenna). In addition to the described particularly strong local field enhancement in the gap, this local field is highly uniform within the volume of the gap (see Supporting Information, Figure S6), which is the consequence of the quasi-static approximation, where the nanogap has similarity with a parallel plate capacitor with nearly uniform field inside. This could be regarded as a practical advantage of the considered a-TSNF over other designs relying on the lightning-rod effect at sharp corners.
We fabricated and experimentally investigated the TSNF (Figure 5a), the terminated TSNF (t TSNF, Figure 5b), and the a-TSNF (Figure 5c) at = 1425 nm. SNOM measurements were compared to numerical simulations by taking a profile of the out-of-plane field component at the altitude of = 45 nm above the silica surface (i.e., 15 nm above the top gold surface). The field component is chosen due to the polarization selectivity of the SNOM, dictated by the elongation and corresponding polarizability of the AFM tip along the out-of-plane -axis. The final non-zero altitude of the profile section is chosen in order to take into account the complicated transfer function of the SNOM (see Supporting Information). (28) As in the previous section, the measured complex near-field pattern of each section with the wide stripe was decomposed into individual guided and leaky modes by the numerical fitting procedure (see Supporting Information, Figure S12). Then the amplitude of the direct-propagating guided stripe mode was used for normalization of the near-field amplitudes. Their x-profiles at the beginning of the taper (the dotted curves in Figure 5d) show close similarity, which reflects the robustness of both the fabrication and fitting procedures. The profiles with the maximum near-field amplitudes for all three structures (the solid curves in Figure 5d), taken along the white dashed lines in Figure 5a-c, demonstrate substantial FE compared to the field of the wide stripe mode. At the same time, the taper-gap-antenna structure demonstrates not only the excitation of the nanoantenna (Figure 5c), but also the strongest observed field enhancement among the three structures (Figure 5d). This observed relative FE (e.g., 1.2 for the a-TSNF, compared to the t-TSNF) should not be confused to the total FE, expected from the simulations (e.g., for the a-TSNF, relative to the t-TSNF), since SNOM measurements do not provide the field inside the gap, where the FE is the strongest. The appropriate line profiles, extracted from the simulated field distributions (inset in Figure 5d), agree well with the experimental results. Further evidence of resonant antenna excitation was the reduction of back reflection. As can be observed, the t-TSNF produces the most pronounced fringes in the near-field amplitude map (Figure 5b), indicating the largest back-reflection. The numerical fitting procedure yielded a substantially small reflectance of 4% for a-TSNF compared to 25% for t-TSNF.
The expected strong FE inside the gap between the nanoantenna and the tip of the TSNF is not observed in the near-field profiles (Figure 5c) because of the polarization selectivity of the SNOM, dictated by the elongation and corresponding polarizability of the AFM tip along the out-of-plane -axis. In order to overcome this restriction, we set a polarizer in front of the detector (Figure 1c) in order to suppress the scattered component and access the in-plane near-field component (the detected signal is denoted as s-polarization). The detected signal with the orthogonal p-polarization, achieved by 90-degrees rotation of the polarizer, corresponds to the maximum contribution of near-field. The recorded near-field maps of the nanoantenna in a TSNF (Figure 5e) clearly show the sharp peak in the profile at , exactly corresponding to the gap, that provides further evidence of the FE inside the gap. The comparison of the measured and profiles with the simulated field distribution is given in Supporting Information, Figure S17.
In summary, we have demonstrated a simple and robust planar configuration, which is amenable for one-step lithography fabrication, for efficient on-chip nanofocusing of radiation carried by stripe SPP modes. The considered structure, consisting of a tapered metal stripe coupled across a gap to a stripe nanoantenna, combines the benefits of near-adiabatic SPP nanofocusing in an optimized tapered structure (characterized by weak scattering losses) with the efficient resonant end-fire excitation of a dipole antenna mode, resulting in further boosting the local FE in the gap. Because of the unique combination of the nanoantenna resonance and optimized SPP nanofocusing (with 41% energy transmission through the stripe taper and the corresponding intensity enhancement of 140), the total intensity enhancement in the gap was deduced to be 12000 with a nearly uniform field over the gap volume (of ). Compared to the excitation of the dimer antenna of similar size under the normal illumination (see Supporting Information, Figure S8), the a-TSNF shows 1.3 times larger intensity enhancement and about 3 times larger inside the gap, normalized to the same total incident power. Supplemented with compactness and background-free hotspot excitation due to the on-chip realization, the a-TSNF design looks superior to the common dimer antenna design.
We would like to emphasize that our concept of a tapered structure (enabling radiation nanofocusing) being coupled (across a gap) to a resonant nanoantenna for further boosting the FE is rather general, and can be applied to many different SPP-based waveguide configurations, such as MIM waveguides, (12) slot waveguides, (10) transmission lines, (17); (18) hybrid plasmonic-photonic waveguides, (11) etc. Additionally, the proposed end-fire coupling can be used for excitation of any nanoantenna, including those supporting dark modes. The latter is illustrated with our simulations shown in Figure 4b, where the second resonance ( nm) corresponds to the second-order antenna mode, which is a dark mode for the far-field excitation. Finally, it should be borne in mind that both nanofocusing and antenna-coupling approaches are not limited to the visible and infrared wavelength ranges, but can also be applied in the THz range, (29); (30); (31) implying the possibility of making use of their combination considered in this work for efficiently concentrating and strongly enhancing the THz radiation in subwavelength-sized hot spots.
The considered combination of a tapered SPP waveguide and a gap-coupled nanoantenna opens up new avenues in the development of integrated nano-optical devices and circuits with a variety of useful on-chip functionalities, in which nanofocusing tapers could serve as efficient inputs/outputs for SPP signals that are to be used for further processing or detection. For example, a nanowire antenna (or a chain of such antennae) coupled to the input and output stripe tapers could be used for spectrally selective transmission (see Supporting Information, Figure S7), thus giving a possibility of combining SPP nanofocusing with wavelength filtering functionalities. In the case of two antennae (three gaps) we observe two local spectral maxima of transmittance corresponding to two different electric field distributions over the three gaps (see Supporting Information). This offers a possibility of adjusting the spatial position of subwavelength-sized hot spots, which could be useful for the realization of efficient spatio-temporal control on the nanoscale. A combination of strong local FE that can be obtained with our configuration and nonlinearities of the dielectric in the gap(s) or intrinsic nonlinearities of metal (gold) could open further opportunities for the development of new nanoscale non-linear components. Overall, bearing in mind the relative simplicity in fabrication and aforementioned features of the proposed structures, we believe that the unique combination of efficient nanofocusing and nanoantenna resonant excitation realized in our experiments promises a variety of useful on-chip functionalities within sensing, nonlinear spectroscopy and signal processing.
V.A.Z., I.P.R., V.S.V., and S.I.B. acknowledge financial support from the Danish Council for Independent Research (the FTP project ANAP, Contract No. 09-072949) and from the European Research Council, Grant No. 341054 (PLAQNAP). A.A. acknowledges financial support from the Danish Council for Independent Research via the GraTer project (Contract No. 0602-02135B).
Supporting Information Available
Details on the fabrication procedure, numerical simulations and optimization, processing of the experimental data. This material is available free of charge via the Internet at http://pubs.acs.org.
- Born, M.; Wolf, E. Principles of Optics 7th edn; Cambridge University Press: Cambridge, 1999; ch. 8.
- Raether, H. Surface Plasmons on Smooth and Rough Surfaces and on Gratings; Springer: Berlin, 1988; ch. 2.
- Takahara, J.; Yamagishi, S.; Taki, H.; Morimoto, A.; Kobayashi T. Opt. Lett. 1997, 22, 475-477.
- Gramotnev, D. K.; Bozhevolnyi, S. I. Nat. Photonics 2010, 4, 83-91.
- Schuller, J. A.; Barnard, E. S.; Cai, W.; Jun, Y. C.; White, J. S.; Brongersma, M. L. Nat. Mater. 2010, 9, 193-204.
- Gramotnev, D. K.; Bozhevolnyi, S. I. Nat. Photonics 2014, 8, 13-22.
- Nerkararyan, K. V. Phys. Lett. A 1997, 237, 103-105.
- Babadjanyan, A. J.; Margaryan, N. L.; Nerkararyan, K. V. J. Appl. Phys. 2000, 87, 3785-3788.
- Stockman, M. I. Phys. Rev. Lett. 2004, 93, 137404.
- Desiatov, B.; Goykhman, I.; Levy, U. Opt. Express 2011, 19, 13150-13157.
- Luo, Y.; Chamanzar, M.; Apuzzo, A.; Salas-Montiel, R.; Nguyen, K. N.; Blaize, S.; Adibi, A. Nano Lett. 2015, 15, 849-856.
- Choo, H.; Kim, M.-K.; Staffaroni, M.; Seok, T. J.; Bokor, J.; Cabrini, S.; Schuck, P. J.; Wu, M. C.; Yablonovitch, E. Nat. Photonics 2012, 6, 838â844.
- Volkov, V. S.; Bozhevolnyi, S. I.; Rodrigo, S. G.; Martin-Moreno, L.; GarcÃa-Vidal, F. J.; Devaux, E.; Ebbesen, T. W. Nano Lett. 2009, 9, 1278â1282.
- Verhagen, E.; Kuipers, L.; Polman, A. Nano Lett. 2007, 7, 334-337.
- Verhagen, E.; Polman, A.; Kuipers, L. K. Opt. Express 2008, 16, 45-57.
- Verhagen, E.; SpasenoviÄ, M.; Polman, A.; Kuipers, L. K. Phys. Rev. Lett. 2009, 102, 203904.
- Schnell, M.; Alonso Gonzalez, P.; Arzubiaga, L.; Casanova, F.; Hueso, L. E.; Chuvilin, A.; Hillenbrand, R. Nat. Photonics 2011, 5, 283â287.
- Sarriugarte, P.; Schnell, M.; Alonso-Gonzalez, P.; Arzubiaga, L.; Golmar, F.; Casanova, F.; Hueso, L. E.; Hillenbrand, R. Opt. Commun. 2012, 285, 3378â3382.
- Andryieuski, A.; Zenin, V. A.; Malureanu, R.; Volkov, V. S.; Bozhevolnyi, S. I.; Lavrinenko, A. V. Nano Lett. 2014, 14, 3925-3929.
- Schnell, M.; Garcia-Etxarri, A.; Huber, A. J.; Crozier, K. B.; Borisov, A.; Aizpurua, J.; Hillenbrand, R. J. Phys. Chem. C 2010, 114, 7341â7345.
- Ocelic, N.; Huber, A.; Hillenbrand, R. Appl. Phys. Lett. 2006, 89, 101124.
- Schnell, M.; Garcia-Etxarri, A.; Alkorta, J.; Aizpurua, J.; Hillenbrand, R. Nano Lett. 2010, 10, 3524â3528.
- Alonso-GonzÃ¡lez, P.; Schnell, M.; Sarriugarte, P.; Sobhani, H.; Wu, C.; Arju, N.; Khanikaev, A.; Golmar, F.; Albella, P.; Arzubiaga, L.; Casanova, F.; Hueso, L. E., Nordlander, P.; Shvets, G.; and Hillenbrand, R. Nano Lett. 2011, 11, 3922-3926
- Liao, P. F.; Wokaun, A. J. Chem. Phys. 1982, 76, 751-752.
- Palik, E. D. Handbook of Optical Constants of Solids; Academic Press: New York, 1985.
- Alonso-GonzÃ¡lez, P.; Albella, P.; Neubrech, F.; Huck, C.; Chen J.; Golmar, F.; Casanova, F.; Hueso, L. E.; Pucci, A.; Aizpurua, J.; Hillenbrand, R. Phys. Rev. Lett. 2013, 110, 203902.
- Zuloaga, J.; Nordlander, P. Nano Lett. 2011, 11, 1280â1283.
- Radko, I. P.; Bozhevolnyi, S. I.; Gregersen, N. Appl. Opt. 2006, 45, 4054-4061.
- Rusina, A.; Durach, M.; Nelson, K. A.; Stockman, M. I. Opt. Express 2008, 16, 18576-18589.
- Zhan, H.; Mendis, R.; Mittleman, D. M. Opt. Express 2010, 18, 9643-9650.
- Iwaszczuk, K.; Andryieuski, A.; Lavrinenko, A.; Zhang, X.-C.; Jepsen, P. U. Opt. Express 2012, 20, 8344-8355.