Simultaneous retrieval of fluidic refractive index and ... - OSA Publishing

Simultaneous retrieval of fluidic refractive index and surface adsorbed molecular film thickness using silicon wire waveguide biosensors Yuki Atsumi,1,2 Dan-Xia Xu,1,* André Delâge,1 Jens H. Schmid,1,2 Martin Vachon,1 Pavel Cheben,1 Siegfried Janz,1 Nobuhiko Nishiyama,2 and Shigehisa Arai2,3 1 National Research Council Canada,1200 Montreal Rd, Ottawa, ON, K1A 0R6, Canada Department of Electrical and Electronic Engineering, Tokyo Institute of Technology, 2-12-1-S3-12 O-okayama, Meguro-ku, Tokyo 152-8550, Japan 3 Quantum Nanoelectronics Research Center, Tokyo Institute of Technology, 2-12-1-S9-5 O-okayama, Meguro-ku, Tokyo 152-8552, Japan * [email protected] 2

Abstract: For silicon wire based ring resonator biosensors, we investigate the simultaneous retrieval of changes in the fluidic refractive index ∆nc and surface adsorbed molecular film thickness ∆dF. This can be achieved by monitoring the resonance shifts of the sensors operating in the TE and TM polarizations at the same time. Although this procedure is straightforward in principle, significant retrieval errors can be introduced due to deviations in the sensor waveguide cross-sections from their nominal values in the range commonly encountered for silicon photonic wire devices. We propose a method of determining the fabricated waveguide size using the group indices derived from measured free spectral range (FSR) of the resonators. We further demonstrate that using experimentally measured group index values, the waveguide size can be determined to accuracies of ± 2 nm in width and ± 1 nm in height. By using this procedure, ∆nc and ∆dF can be obtained to a precision of within 10% of the true values using optically measurable parameters, improving the retrieval accuracy by more than 3 times. ©2012 Optical Society of America OCIS codes: (130.3120) Integrated optics devices; (130.6010) Sensors; (230.5750) Resonators.

References and links 1. 2. 3. 4. 5. 6. 7. 8.

X. Fan, I. M. White, S. I. Shopova, H. Zhu, J. D. Suter, and Y. Sun, “Sensitive optical biosensors for unlabeled targets: a review,” Anal. Chim. Acta 620(1-2), 8–26 (2008). K. Matsubara, S. Kawata, and S. Minami, “Optical chemical sensor based on surface plasmon measurement,” Appl. Opt. 27(6), 1160–1163 (1988). W. Lukosz, “Priciples and sensitivities of integrated optical and surface plasmon sensors for direct affinity sensing and immunosensing,” Biosens. Bioelectron. 6(3), 215–225 (1991). J. Escorihucla, M. J. Banuls, J. G. Castello, V. Toccafondo, J. G. Ruperez, R. Puchades, and A. Maquicira, “Chemical silicon surface modification and bioreceptor attachment to develop competitive integrated photonic biosensors,” in Analytical and Bioanalytical Chemistry, L. M. Lecuga ed. (Springer, 2012) S. Janz, A. Densmore, D.-X. Xu, P. Waldron, J. Lapointe, J. H. Schmid, T. Mischki, G. Lopinski, A. Delâge, and R. Mckinnon, “Silicon photonic wire waveguide sensors,” in Advanced Photonic Structures for Biological and Chemical Detection, F. Xudong, ed. (Springer, 2009). M. Iqbal, M. A. Gleeson, B. Spaugh, F. Tybor, W. G. Gunn, M. Hochberg, T. Baehr-Jones, R. C. Bailey, and L. C. Gunn, “Label-free biosensor arrays based on silicon ring resonators and high-speed optical scanning instrumentation,” IEEE J. Sel. Top. Quantum Electron. 16(3), 654–661 (2010). W. Bogaerts, P. De Heyn, T. Van Vaerenbergh, K. De Vos, S. Kumar Selvaraja, T. Claes, P. Dumon, P. Bienstman, D. Van Thourhout, and R. Baets, “Silicon microring resonators,” Laser Photonics Rev. 6(1), 47–73 (2012). D.-X. Xu, A. Densmore, A. Delâge, P. Waldron, R. McKinnon, S. Janz, J. Lapointe, G. Lopinski, T. Mischki, E. Post, P. Cheben, and J. H. Schmid, “Folded cavity SOI microring sensors for high sensitivity and real time measurement of biomolecular binding,” Opt. Express 16(19), 15137–15148 (2008).

#176366 - $15.00 USD

(C) 2012 OSA

Received 18 Sep 2012; revised 6 Nov 2012; accepted 6 Nov 2012; published 15 Nov 2012

19 November 2012 / Vol. 20, No. 24 / OPTICS EXPRESS 26969

9. 10. 11.

12. 13. 14. 15. 16. 17.

A. Densmore, M. Vachon, D.-X. Xu, S. Janz, R. Ma, Y.-H. Li, G. Lopinski, A. Delâge, J. Lapointe, C. C. Luebbert, Q. Y. Liu, P. Cheben, and J. H. Schmid, “Silicon photonic wire biosensor array for multiplexed realtime and label-free molecular detection,” Opt. Lett. 34(23), 3598–3600 (2009). D.-X. Xu, M. Vachon, A. Densmore, R. Ma, S. Janz, A. Delâge, J. Lapointe, P. Cheben, J. H. Schmid, E. Post, S. Messaoudène, and J.-M. Fédéli, “Real-time cancellation of temperature induced resonance shifts in SOI wire waveguide ring resonator label-free biosensor arrays,” Opt. Express 18(22), 22867–22879 (2010). A. Densmore, D.-X. Xu, N. A. Sabourin, H. Mcintosh, P. Cheben, J. H. Schmid, R. Ma, M. Vachon, A. Delâge, W. Sinclair, J. Lapointe, Y. Li, G. Lopinski, B. Lamontagne, and S. Janz, “A fully integrated silicon photonic wire sensor array chip and reader instrument,” in Proceedings of IEEE Conference on Group IV Photonics (The Royal Society, London, 2011), pp. 350–352. K. Tiefenthaler and W. Lukosz, “Sensitivity of grating couplers as integrated-optical chemical sensors,” J. Opt. Soc. Am. B 6(2), 209–220 (1989). D. R. Shankaran, K. V. Gobi, and N. Miura, “Recent advancements in surface plasmon resonance immunosensors for detection of small molecules of biomedical, food and environmental interest,” Sens. Actuators B Chem. 121(1), 158–177 (2007). S. Selvaraja, W. Bogaerts, P. Dumon, D. Van Thourhout, and R. Baets, “Sub-nanometer linewidth uniformity in silicon nano-photonic waveguide devices using CMOS fabrication technology,” IEEE J. Sel. Top. Quantum Electron. 16(1), 316–324 (2010). A. V. Krishnamoorthy, X. Zheng, G. Li, J. Yao, T. Pinguet, A. Mekis, H. Thacker, I. Shubin, Y. Luo, K. Raj, and J. E. Cunningham, “Exploiting CMOS manufacturing to reduce tuning requirements for resonant optical devices,” IEEE Photonics J. 3(3), 567–579 (2011). D.-X. Xu, A. Delâge, J. H. Schmid, R. Ma, S. Wang, J. Lapointe, M. Vachon, P. Cheben, and S. Janz, “Selecting the polarization in silicon photonic wire components,” Proc. SPIE 8266, 82660G, 82660G-9 (2012). A. Delâge, D.-X. Xu, R. W. McKinnon, E. Post, P. Waldron, J. Lapointe, C. Storey, A. Densmore, S. Janz, B. Lamontagne, P. Cheben, and J. H. Schmid, “Wavelength-dependent model of a ring resonator sensor excited by a directional coupler,” J. Lightwave Technol. 27(9), 1172–1180 (2009).

1. Introduction Precise detection of trace amounts of biomolecules is required for many applications including diagnostics, genomics, proteomics, and drug screening. Today optical biosensors are powerful analysis tools for such applications by means of fluorescence-based or label-free detection. In contrast to fluorescence-based techniques, label-free sensors allow for quantitative and kinetic measurement of molecular interaction [1]. While surface plasmon resonance (SPR) is commonly used in commercial sensing systems [2–4], recently silicon photonic wire waveguide evanescent field (PWEF) sensors have seen rapid development [5– 7]. Due to the high refractive index of silicon, the optical modal electric field is strongly localized near the waveguide surface, giving rise to high sensitivity to surface adsorption. Long cavity resonator sensors can be made with small footprint, making them easily scalable to large arrays [8]. Silicon PWEF sensors also benefit from established silicon photonics technologies developed in the past decade, allowing low cost fabrication. We have developed highly sensitive silicon-on-insulator (SOI) photonic wire sensors with folded path in ring resonator and Mach-Zehnder interferometer configurations, operating in the TM polarization [8–11]. With these structures, we have demonstrated the real-time detection of surface-adsorbed DNA and protein molecules with a sensitivity of ~0.2 pg/mm2, an order of magnitude better than the SPR technique, and the sensors are scaled up to × 128 arrays. The temperature-induced resonance shift is cancelled out using reference resonators which is protected from the sensing medium by a SU8 cladding layer [10]. In these label-free sensors, however, the changes in the refractive index ∆nc of the aqueous solution and the adsorbed molecular layer thickness ∆dF on the waveguide surface cannot be distinguished directly, similar to the case for SPR. In our previous experiments, the refractive index of the buffer solution was usually kept constant, and the measured sensing responses were attributed to the surface adsorption only. In general, several solutions of different refractive indices may be required for the immobilization and rinsing steps in sensing procedures, so the response of the PWEF sensors comes from both ∆nc and ∆dF [12]. In commercial SPR instruments, the fluidic index changes are monitored using a separate fluidic channel, requiring more complex instrumentation or regent preparation [13]. It is also difficult to ensure that the response of the reference channel solely comes from the fluidic index change. It is desirable to differentiate these responses optically for the same fluid flow.

#176366 - $15.00 USD

(C) 2012 OSA



In this work, we demonstrate the feasibility of retrieving these two changes (∆nc and ∆dF) simultaneously using experimentally measurable resonance-wavelength shifts of two sensors through a single fluidic channel, operating in the quasi transverse electric (TE) and transverse magnetic (TM) fundamental modes, respectively. The main challenge in this procedure comes from the waveguide size fluctuations common in silicon photonic wire fabrication, which can cause large errors in the retrieved values. We propose and numerically demonstrate that the actual waveguide size can be determined with high precision using measured group indices. Using this information, ∆nc and ∆dF can be obtained with good accuracy. 2. Theoretical background Water

W

dF dF

Adsorbed organic layer SiO2

Si

H Water (nc=1.325) Organic (n=1.5) Si (n=3.476) SiO2(n=1.444)

TM

TE

Fig. 1. (a) Cross-sectional view of a silicon photonic wire waveguide sensor. (b) Schematic layout of our spiral cavity resonator sensors for TE and TM polarizations.

The silicon wire waveguide structure in our resonator sensors is shown in Fig. 1(a). For comparison with our experimental work, the nominal waveguide size is chosen as 450 nm wide (W) and 260 nm high (H). The upper cladding is assumed to be an aqueous solution unless otherwise specified. Due to the different optical field distributions of the TE and TM modes, the sensing response to the variations in the cladding fluidic refractive index ∆nc and adsorbed molecular-film thickness ∆dF differ between the two polarizations. The correspondingly induced effective index changes for TE and TM modes are:  ∂N eff(TE) ∂N eff(TE)    ∂d F   Δnc   ΔN eff(TE)   ∂nc ⋅ (1)  = .    ΔN eff(TM)   ∂N eff(TM) ∂N eff(TM)   Δd F   ∂n ∂d F  c  In experiments, we can obtain the values of the resonance wavelength shift ∆λ and the free-spectral-range (FSR) for TE and TM modes in real time. The effective index changes ∆Neff can then be calculated from ∆λ using the relation ∆λ = λ(∆Neff / Ng), where the group index Ng = λ2/(FSR⋅L), and λ and L are the resonance wavelength and the resonator cavity length, respectively. The values of ∂Neff(TE,TM)/∂nc and ∂Neff(TE,TM)/∂dF can be calculated through simulations, provided that the waveguide size is known. Here the sensors operating in the TE and TM polarizations are assumed to have identical waveguide dimensions. Since these sensors are generally placed within a few hundred micrometers of each other, this is a reasonable assumption based on experimental observations. For a given waveguide crosssection, simulations show that ∂Neff(TE,TM)/∂nc is a constant for a wide range of fluidic refractive index. The value of ∂Neff(TE,TM)/∂dF can also be closely approximated by a constant for ∆dF < 10 nm, which is larger than the size of most proteins and DNA molecules encountered in biosensing and therefore adequately cover the range of interest. Using Eq. (1), it is straightforward to find the solutions for ∆nc and ∆dF as expressed using an inverse sensitivity matrix S −1:

#176366 - $15.00 USD

(C) 2012 OSA



∂N eff(TE)   ∂N eff(TM) −   ∂d F   ΔN eff(TE)   Δnc   ∂d F = ⋅ ⋅ S (1/ det )     ∂N eff(TM) ∂N eff(TE)   ΔN eff(TM)   Δd F   −  ∂nc ∂nc   ∂N eff(TE) ∂N eff(TM) ∂N eff(TE) ∂N eff(TM) ⋅ − ⋅ ≠0 det S = ∂nc ∂d F ∂d F ∂nc

(a)

(2)

(b)

(c)

(d)

Fig. 2. Elements of the inverse sensitivity matrix versus (a)(c) waveguide width at a height of 260 nm and (b)(d) waveguide height at a width of 450 nm.

3. Simulation results and discussion It is well known that the optical properties of silicon wire waveguides strongly depend on their cross-sectional dimensions. Figures 2(a) and 2(b) show the dependence of the inverse sensitivity matrix elements on the waveguide width for a height of 260 nm and the waveguide height for a width of 450 nm, calculated using the film mode matching (FMM) method at a wavelength of 1.55 µm. As expected, the matrix elements for the TE polarization vary with W more strongly than those for TM which have a stronger dependence on the height H. That is due to the high index contrast which leads to the large surface electric field magnitude and discontinuity at the waveguide interface in the surface normal direction. It is also evident that the sensing response of the TM mode is higher than that of the TE mode for W > 400 nm. The maximum response for the TM mode to ∆nc and ∆dF occurs at H of 220 nm and 245 nm, respectively. In order to obtain valid solutions from Eq. (2), it is required that the matrix determinant det S ≠0. In other words, when the combined responses of the TE and TM modes to ∆nc and ∆dF are the same, the differentiation between these two contributions becomes impossible. The dependence of 1/det S on the waveguide dimensions is shown in Figs. 2(c) and 2(d). The values of 1/det S should be as small as possible for best numerical accuracies. The waveguide

#176366 - $15.00 USD

(C) 2012 OSA



dimensions (W ~300 nm or H ~200 nm) for minimum 1/det S, however, leads to other tradeoffs such as large waveguide bend radius. Figures 3(a) and 3(b) show the overall coefficients for variables ∆Neff(TE) and ∆Neff(TM) in Eq. (2) as a function of W and H. For a given waveguide size, ∆nc and ∆dF can easily be calculated using these coefficients. These results also illustrate the high sensitivity to waveguide dimensions of the retrieved ∆nc and ∆dF using measured ∆Neff(TE) and ∆Neff(TM). Good structural tolerance is found for waveguides with e.g. W ~350 nm for H = 260 nm, or H ~200 nm for W = 450 nm, where the optical modes are rather delocalized for at least one polarization. On the other hand, the sensor waveguide choice also depends on requirements including the single mode operation, high sensing response and small bend radius to ensure compact footprint. Our nominal design of 450 nm × 260 nm represents a compromise between all these considerations.

(a)

(b) Fig. 3. Overall elements of the inverse sensitivity matrix S −1 versus (a) waveguide width at a height of 260 nm and (b) waveguide height at a width of 450 nm.

#176366 - $15.00 USD

(C) 2012 OSA



Using current fabrication technologies, the silicon thickness (H) uniformity on a wafer is typically ± 2 nm, while the variation between wafers is typically ± 5 nm. For state-of-the-art patterning techniques such as e-beam and 193 nm DUV lithography, the deviation of the waveguide width W is typically ± 10 nm, although a systematic bias can often be larger [14– 16]. The precise sensor waveguide dimensions are generally not available for each chip. In the following, we evaluate the retrieval accuracies using the nominal waveguide dimension of 450 nm × 260 nm in the presence of size fluctuations.

(a)

(b) Fig. 4. Evaluation of the retrieval accuracy as a function of waveguide dimension changes in (a) waveguide width and (b) waveguide height. Two cases were estimated: ∆Neff(∆nc): ∆Neff(∆nc) ~“10: 1” and “1: 1”.

As shown in Fig. 4, the vertical axis is defined as the normalized difference between the original (∆nc and ∆dF) and the retrieved values (∆nc’ and ∆dF’), with corresponding deviations of (∆nc − ∆nc’)/∆nc or (∆dF − ∆dF’)/∆dF. Two cases are considered. In the first case the maximum value for ∆nc is 0.003 which is comparable to the index difference between 2 × PBS (phosphate buffered saline) solution and water, and ∆dF of 1.2 nm is comparable to 1/3

#176366 - $15.00 USD

(C) 2012 OSA



of monolayer of streptavidin protein. The corresponding changes in the effective index for the TM mode are ∆Neff(∆nc) ~0.0012 and ∆Neff(∆dF) ~0.001, respectively. In other words, the contribution to the sensing signal from the fluidic index change and surface adsorption are comparable: ∆Neff(∆nc):∆Neff(∆dF) ~1:1. A more stringent case is also evaluated where ∆Neff(∆nc):∆Neff(∆dF) ~10:1. First, the response in ∆Neff(TE) and ∆Neff(TM) is calculated from Eq. (1) to emulate experimentally measured values using the waveguide dimensions as indicated in the Fig. 3. Then the values of ∆nc and ∆dF are retrieved using Eqs. (2) and (3) assuming the nominal waveguide size of 450 nm × 260 nm. As the waveguide size deviates from the nominal dimensions, significant errors in the retrieved fluidic index change or the surface adsorption are introduced, particularly for the component with the smaller sensing response. For the case of ∆Neff(∆nc):∆Neff(∆dF) ~10:1, the maximum errors are as large as about + 150% and −100%, respectively, for ∆W~ ± 10 nm and ∆H ~ ± 5 nm. Clearly, this level of errors is unacceptable. It is imperative that the actual waveguide dimensions are accurately known. In the following, we demonstrate that the fabricated waveguide dimensions can be determined with high accuracy using the group index Ng for the TE and TM polarizations that are deduced from the optically measurable FSR data. The measurement accuracy of FSR depends on the magnitude of FSR, the resonator quality factor and the extinction ratio. Therefore here we choose air as the upper cladding to avoid the degradation in the quality factor caused by the attenuation from water. It should be noted that air upper cladding is generally available during the setup procedure of sensing operations. Figure 5 shows the group indices for TE and TM modes as a function of W and H, respectively. As can be seen, the changes in Ng(TE) almost depends solely on W, while the changes in Ng(TM) depends mainly on H, allowing W and H being determined separately. The color-contours are separated by a group index interval of 0.01 which is experimentally achieved, as shown below.

Fig. 5. Group indices as the functions of waveguide width and height for (a) TE mode and (b) TM mode. The upper cladding is water and the contour interval was set to be 0.01 in both figures.

The transmission spectrum of a spiral cavity resonator (see Fig. 1(b)) in TM mode is shown in Fig. 6(a), measured with a wavelength scan step of 1 pm. Here the cavity length L = 700 µm which result in a FSR of ~0.7 nm. The quality factor is of a moderate value of Q ~20,000. The deduced Ng as well as the deviations from the polynomial fit is shown in Fig. 6(b). Clearly the group index can be determined to within an accuracy of ± 0.005. The measurement accuracy in Ng improves with decreasing L, but we found that similar precision can be achieved for resonators with cavity length L < 1000 µm.

#176366 - $15.00 USD

(C) 2012 OSA



Fig. 6. (a) Measured spectrum for spiral cavity resonator with the cavity length of 700 µm, and (b) calculated group index from the measured spectrum, and it’s deviation from the polynomial fit.

Within the range of waveguide dimensions considered, Ng is a monotonic function of both W and H that can be approximated by a 2nd order polynomial fit:  ∂N g(TE,TM)   ∂N g(TE,TM)  N g(TE,TM) = N g0(TE,TM) + ( H − H 0 )   + (W − W0 )   +  ∂H   ∂W  2 2   ∂ 2 N g(TE,TM)  2 ∂ N 2  ∂ N g(TE,TM)  W W H H W W + − + − − ( H − H 0 )  g(TE,TM) ( ) ( )( )     0 0 0  2 2  2∂H   2∂W   ∂H ∂W  (3)

where Ng0 is the group index for W0 × H0 = 440 nm × 255 nm as a reference point. The excellent agreement between the fitted values and the original data (standard deviation of 6.16 × 10−6 for TE and 1.11 × 10−5 for TM) suggests that H and W could also be retrieved from a second degree expansion in terms of TE and TM group indices as:  ∂W , H   ∂W , H  W , H = [W0 , H 0 ] + ( N g(TE) − N g0(TE) )   + ( N g(TM) − N g0(TM) )   +  ∂N g(TE)   ∂N g(TM)  2 2  2  2  ,H  ( N g(TE) − Ng0(TE) )  2∂∂WN 2, H  + ( Ng(TM) − Ng0(TM) )  2∂∂W  +  (4) 2   N g(TM)  g(TE)   2   ( Ng(TE) − N g0(TE) ) ⋅ ( N g(TM) − N g0(TM) )  ∂N ∂ W∂,NH  .  g(TE) g(TM) 

An optimization program, based on the Genetic Algorithm [17] is used to find the values of the 12 parameters of Eq. (4) represented in square brackets. The optimized coefficients are given in Table. 1. Table 1. Optimized fitting parameters for Eq. (4). Type

Ng0

TE4.415167 mode TM4.808605 mode * Standard deviation

∂Ng/∂H

∂Ng/∂W

½∂2Ng/∂H2

−7.16E-05

−2.86E-03

−8.37E-06

½∂2Ng/∂W2 7.19E-06

2.37E-02

2.00E-03

−5.828E-04

−1.230E-05

∂2Ng/∂H∂W

Error*

−1.03E-06

6.16E-06

−1.163E-04

1.11E-05

Figures 7(a) and 7(b) show the W and H matrix as a function of group indices for TE and TM mode which are calculated from Eq. (4) with the use of the fitting parameters listed in Table 1. Here the color contours are separated by a waveguide width ∆W of ± 2 nm in Fig.

#176366 - $15.00 USD

(C) 2012 OSA



7(a) and by a waveguide height ∆H was of ± 1 nm from Fig. 7(b). We can determine the size of the waveguide with these numerical expressions. With the obtained level of accuracy in measuring Ng, we can determine the waveguide size to an accuracy of ± 2 nm for W and ± 1nm for H.

Fig. 7. Distribution of (a) waveguide width and (b) waveguide height as the functions of group indices for TE and TM mode with the color interval of 2 nm and 1 nm, respectively. The squares with white dashed line show the range of measurement accuracy for Ng.

Finally, in order to evaluate the differentiation method, we evaluated the retrieval accuracy for fluidic index and film thickness changes again using the determined actual waveguide dimensions, as shown in the shaded areas of Fig. 4. In the case of ratio of 1:1, the maximum retrieval errors are reduced from 30% to 7%. Even when the effective index variation due to liquid index change exceeds the variation due to film thickness change by a factor of 10, the maximum errors for the smaller component are reduced from 150% to less than 40%, i.e. the retrieval accuracy is improved by more than 3 times by using our procedure. 4. Conclusion We have presented a robust method of simultaneous monitoring of the fluidic index change ∆nc and surface adsorption ∆dF using silicon photonic wire resonator sensors. The data retrieval is achieved by measuring the resonance wavelength shifts for the TE and TM modes at the same time, and correlating them to nc and ∆dF using calculated sensitivity coefficients. We also showed how to circumvent a critical challenge in ∆nc and ∆dF discrimination accuracy, namely the large dependence of the sensitivity coefficients on the cross-section dimensions of the sensor waveguide core. Waveguide size fluctuations commonly encountered in silicon photonic wire fabrication can introduce unacceptable retrieval errors. We propose a new method of accurately determining the size of fabricated waveguides using TE and TM mode group indices obtained from measured resonator FSR. Using this approach, the waveguide dimensions can be determined to a precision of ∆W ~ ± 2 nm and ∆H ~ ± 1 nm. By applying the actual waveguide size, the retrieval accuracies for ∆nc and ∆dF are improved by more than 3 times. The proposed method of determining waveguide dimensions indeed is not restricted to evanescent field waveguide sensors and it can also be used in other types of planar waveguide circuits. Acknowledgments This research was supported in part by the Genome and Health Initiative (GHI) program of the National Research Council Canada (NRC), by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) and the Japan Society for the Promotion of Science (JSPS) under Grants-in-Aid for Scientific Research (#19002009, #22360138, #21226010, and #11J08863) and by MEXT Global-COE program “Photonics integration-Core Electronics”.

#176366 - $15.00 USD

(C) 2012 OSA

