Open Access
Add to BibSonomyAbstract
We design and fabricate ultra-low-loss tapered optical fibers (TOFs) with minimal lengths. We first optimize variations of the torch scan length using the flame-brush method for fabricating TOFs with taper angles that satisfy the adiabaticity criteria. We accordingly fabricate TOFs with optimal shapes and compare their transmission to TOFs with a constant taper angle and TOFs with an exponential shape. The highest transmission measured for TOFs with an optimal shape is in excess of 99.7% with a total TOF length of only 23mm, whereas TOFs with a constant taper angle of 2mrad reach 99.6% transmission for a 63mm TOF length.
© 2014 Optical Society of America
1. Introduction
Tapered optical fibers with a sub-wavelength waist radius have been utilized for diverse applications in photonics and quantum optics [1, 2]. Strong confinement of the guided mode to a sub-wavelength radius leads to efficient nonlinear effects such as supercontinuum generation [3–5] and third-harmonic generation [6–8]. Large evanescent fields of the guided mode are beneficial for optical sensing and spectroscopy [9–14], for realizing strong light-matter interaction [15–22], and for coupling to whispering-gallery-mode microresonators [23,24].
For these applications, minimizing TOF losses is important and sometimes crucial. In addition to extrinsic losses caused by surface roughness or contaminants, there are more fundamental TOF losses, caused by the power coupling from the fundamental mode to higher-order modes because of the non-zero taper angles. These power coupling losses depend on the TOF shape; TOFs with smaller taper angles, i.e., more adiabatic shapes, have lower power coupling losses. Transmissions of 99.4% and 99.95% have been achieved using TOFs with an exponential shape with a decay constant of 10mm [25] and using TOFs with a constant taper angle of 2mrad [26, 27], respectively. However, TOFs with small taper angles tend to be long (the lengths of the exponential TOF in Ref. 25 and the 2mrad TOF in Ref. 27 are 111mm and 84mm, respectively); there is a trade-off between the short length and high transmission of TOFs.
For some applications, e.g., fiber coupling to microresonators, shorter TOFs are desirable because of their higher mechanical stability. Moreover, there is sometimes a constraint on the length of the TOFs when they are installed in small spaces, such as cryostats [28] or vacuum chambers [29]. Therefore, it is desirable to restrict the length of the TOFs as much as possible while maintaining high transmission. Useful guidelines for designing such TOFs are the adiabaticity criteria [30–32], which give the upper limit for the local taper angle as a function of the local taper radius. There have been studies reporting the fabrication of TOFs containing three tapered sections where each section has a different taper angle such that the taper angle is smaller where the adiabaticity criterion gives a smaller upper limit [26, 27, 33]. Although transmissions higher than 99% have not been achieved with such trilinear TOFs, they can in principle have shorter TOF lengths than linear TOFs without sacrificing the transmission qualities. However, there have been no reports on the design and fabrication of TOFs with the shapes optimized in systematic ways for achieving minimal TOF lengths.
In this paper, we design and fabricate TOFs that optimally satisfy the adiabaticity criterion. We adopt the flame-brush method [34] using variable torch scan lengths, and we optimize the scan length such that the local taper angle satisfies the adiabaticity criterion by taking a systematic approach. We accordingly fabricate TOFs with optimal shapes and compare their transmission with those having a constant taper angle of 2mrad and those having an exponential shape with a decay constant of 3mm. To investigate the reproducibility, five samples are fabricated for each profile. The highest transmission measured for TOFs with an optimal shape reaches higher than 99.7% with a total TOF length of only 23mm, whereas TOFs with a constant taper angle of 2mrad reach 99.6% transmission for a 63mm TOF length.
2. Theory
2.1. Adiabaticity criteria
A TOF is adiabatic if the taper angle is small enough such that there is negligible power coupling from the fundamental mode to higher-order modes. Adiabaticity criteria have been proposed to delineate such conditions [30–32]. The delineation angle Ω(r) as a function of the local taper radius r is given by the following length-scale criterion [31]:
where β1(r) and β2(r) are the propagation constants of the fundamental (HE11) and the first excited (HE12) modes, respectively. If the local taper angle θ(r) is much smaller than Ω(r), the mode coupling is negligible and the fundamental mode propagates adiabatically. We introduce an adiabaticity factor F(> 0) to express this condition as follows: Smaller F values lead to reduced mode coupling. Figure 1 shows the delineation angle of a step-index fiber calculated with Mathematica (Wolfram Research) using the three-layer model [35] with the wavelength λ = 852nm and with the refractive indices and the radiuses of the core (cladding) of ncore = 1.4574 (nclad = 1.4525) and rcore = 2.4μm (rclad = 62.5μm), respectively.
Fig. 1 The calculated delineation angle Ω as a function of the fiber radius r of a step-index fiber for the wavelength λ = 852nm, where the refractive indices and the radii of the core (cladding) are ncore = 1.4574 (nclad = 1.4525) and rcore = 2.4μm (rclad = 62.5μm), respectively.
2.2. Model for the TOF shape
We model the shape of a TOF fabricated using the flame-brush method [34] with a variable torch scan length. In this method, a small flame heats a small region of a fiber with an initial radius r0 to the softening point. While the fiber is pulled from both ends with a constant velocity v1, the flame scans back and forth along the fiber axis with a constant velocity v2. The length of the k-th scan is Lk. Here we count k by the number of one-way scans, not by the number of roundtrip scans. This scanning process elongates the heated region of the fiber, and its radius is reduced to rk after the k-th scan of the torch.
Figure 2 shows the schematic of the TOF after the k-th scan. We assume that the flame width is small enough that a negligibly short length of the fiber is heated at any given time. As shown in Fig. 2, the cylinder with radius rk−1 and length between planes labeled A and B is heated during the k-th scan and become a cylinder with radius rk and length between planes labeled A′ and B′ at the end of the k-th scan. The conservation of the volume before and after the k-th scan,
leads to the a recurrence relation for the fiber radii: The local taper angle θk after the k-th scan can be derived from Fig. 2 as Note that this model is invalid if the torch passes the end of the heated region for the (k − 2)-th scan, . This requirement poses the limitation, The total TOF length Ltotal after the n-th scan can be written asFig. 2 Schematic of the TOF elongation after the k-th scan of the torch with the flame-brush method. The origin of the z-axis is taken as the position of the torch at the beginning of the k-th scan (i.e., at the end of the (k − 1)-th scan). The torch scans in the positive direction over a length of Lk, and the radius is reduced from rk−1 to rk after the k-th scan.
2.3. Optimization of the heat length profile
To design an adiabatic TOF with minimal length, Ltotal must be minimized under the condition θk < FΩ(r = rk) for any k. For given values of r0 and v1/v2, this minimization is an optimization of Lk (k = 1, 2, ··· ,n). Note that a smaller value of the ratio v1/v2 leads to smaller staircase-step changes in the fiber radius (see Eq. (4)), thus leading to a smoother taper at the cost of larger number of total scans n. In our typical experimental conditions, there are approximately 300 total scans, n ∼ 300, which has too many degrees of freedom for a straightforward optimization. To reduce the degrees of freedom, we assume that Lk is a smooth function of k, which should be a reasonable assumption, considering the smooth shape of Ω(r) (see Fig. 1). Specifically, we consider a Lagrange interpolating polynomial, which passes through nine points that are fixed with even intervals on the scan-number axis k and are variable on the L axis, as Lk.
In addition to the conditions given by Eqs. (2) and (6), we pose the following conditions. First, we introduce the upper limit for Lk of Lupper = 40000μm to take into account of the limitations of the experimental setup. Second, we introduce the lower limit for Lk of Llower = 400μm to take into account of the finite width of the flame. Furthermore, because the delineation angle Ω(r) abruptly increases as r becomes less than ∼ 4μm, we fix Lk as Lk = Llower for rk ≲ 4μm. This restriction gives better results than optimizing Lk for the whole indices of k.
We optimize Lk with the downhill simplex method [36], starting from ten sets of random initial values. The numerical calculation is performed with custom scripts using Matlab (Math-Works). Figure 3 shows the optimized profiles of Lk with the adiabaticity factors of F = 0.2, 0.4, 0.6, and 0.8. The TOFs with larger F values require longer scan lengths, which leads to longer TOF lengths.
3. Experiment
We experimentally fabricate TOFs according to the optimized Lk profile obtained in the previous section. A commercial step-index single-mode fiber (Thorlabs, SM800) is fixed with two fiber clamps (Thorlabs, T711/M-250) on stepper-motor-driven linear translation stages (SURUGA SEIKI, KX1250), which pull the fiber on both ends at a velocity of v1 = 15μm/s. The fiber is heated with a flame from a hydrogen torch whose stainless-steel nozzle has a single hole with an inner diameter of 135μm. We use pure hydrogen gas (no premix with oxygen), and the gas flow is adjusted to 10mL/min by a flow meter (KOFLOC, RK1250). The torch is also mounted on a stepper-motor-driven stage and scans along the fiber with a velocity of v2 = 750μm/s. Light from an 852nm laser diode is coupled to a 99/1-fiber beam splitter, and the 99% output port is fusion spliced to the fiber being pulled. We stabilize the laser output power to within ±0.1% by monitoring the power from the 1% port and feeding it back to the diode current. The TOF transmission is measured throughout the pull.
We fabricate TOFs optimized for the adiabaticity factors of F = 0.2, 0.4, 0.6, and 0.8. For comparison, we also fabricate TOFs with two other shapes. One TOF has a constant taper angle of 2mrad for r > 4μm and an exponential shape with a waist length of 400μm for r < 4μm. Note that all of the optimized TOFs also have an exponential shape with a waist length of 400μm for r ≲ 4μm because we fix Lk as Lk = 400μm for rk ≲ 4μm. The other TOF has an exponential shape with a waist length of 3mm. Five samples are fabricated for each profile. All samples are tapered to a waist radius that is thin enough to be single-mode guided (r ≲ 300nm).
Figure 4 shows the change of the transmission during the pull for a TOF with F = 0.4. The amplitude of the oscillation due to the coupling to higher-order modes is small, and the final transmission is in excess of 99.7%. The disappearance of the oscillation is clearly observed at the end of the pull, which is also confirmed with the short-time Fourier transform (spectrogram) of the transmission trace [37]. The disappearance of the oscillation indicates the single-mode guidance.
Fig. 4 The transmission change during the pull of the tapered fiber optimized with F = 0.4. The inset shows the magnified plot for the last ∼ 0.6mm of the pull.
Figure 5 summarizes the transmissions of all the fabricated TOFs. There is a tendency that the profiles optimized with larger F values tend to have larger losses. Note that there are no noticeable differences in transmissions among the profiles of F = 0.2, F = 0.4, and the 2mrad, in spite of the different taper angle profiles. This suggests that for these TOFs the losses due to mode coupling are sufficiently suppressed, and other extrinsic losses are dominant. Indeed, Ravets et al. reported the transmissions of 99.95% for 2mrad TOFs [26,27]; the highest transmission measured in our experiment for 2mrad TOFs is only 99.6%. We suspect that the extrinsic losses for our TOFs are primarily caused by contaminants on the TOF surface because we use a homemade clean hood of class 100000, compared to the class 100 cleanroom used in [26,27]. Furthermore, perturbations such as air, mechanical vibrations of the setup, and fluctuation of the gas flow may cause deviation of the TOF shape from the theoretical prediction, leading to less adiabaticity and larger losses.
Fig. 5 The transmissions of fabricated TOFs with the adiabaticity factors of F = 0.2, 0.4, 0.6, and 0.8, a constant taper angle of 2mrad, and an exponential shape with a decay constant of 3mm. Five TOFs are fabricated for each profile. The inset shows the transmissions for F = 0.2, 0.4, and 2mrad with a magnified scale for the vertical axis. The relative error in the transmission is ±0.1%, which is limited by the laser output power stability.
We measure the shape of the fabricated TOFs using a scanning electron microscope (SEM). Figure 6(a) and (b) show the measured shapes for the F = 0.4 and 2mrad TOFs. In practice, the TOFs have to be pulled over a longer distance than that predicted by the above model to achieve single-mode guidance at the waist. Longer pulling is required because of the effect of the finite flame width that is not negligible compared to the relatively small scan length Llower. Thus, we extend the scan length profiles Lk to a larger scan number until single-mode guidance is reached, which is confirmed by the observed disappearance of the oscillation in the transmission (see Fig. 4). The predicted shapes in Fig. 6 are calculated with experimentally executed scan numbers. The TOF radii are measured as the half-width at the half-maximum of the each fiber edge slope in the intensity profile of the SEM images, as shown in Fig. 6(c)–(f). The upper (lower) error bars are determined as the half-width at 10% (90%) of the maximum intensity of each slope. The SEM images at around the waist of TOFs (r ≲ 4μm) are blurry due to the mechanical vibrations of the TOF, and the measured radii inevitably have large errors and may be overestimated (Fig. 6(c) and (e)), compared to those at r ≳ 4μm (Fig. 6(d) and (f)). However, because the delineation angle Ω(r) abruptly increases at r ≲ 4μm, precise knowledge of the TOF radius at these regions is not crucial to our study. For the regions with r ≳ 4μm, the shapes of the fabricated TOFs agree well with those predicted by the model with a normalized residual of ≲ 10%. Moreover, in spite of the higher transmission, the F = 0.4 TOF has almost three times shorter length (23mm) than that for the 2mrad TOF (63mm). Note that the error bars around the waist of the F = 0.4 TOF are smaller than that for the 2mrad TOF because F = 0.4 TOF has smaller mechanical vibrations due to its shorter TOF length.
Fig. 6 (a) Comparison of the predicted and measured TOF shapes. Red and blue colored plots represent the shapes optimized for the F = 0.4 and the 2mrad profiles, respectively. The solid lines are the theoretically predicted curves. (b) The magnified plot for the region around the waist (r ≲ 4μm) of (a). Examples of SEM images of the 2mrad TOF (c) around the waist and (d) around r ≈ 20μm, respectively. (e) and (f) show the intensity profiles of (c) and (d), respectively.
Figure 7 shows the taper angles as a function of TOF radius for the fabricated F = 0.4 TOF and 2mrad TOF, along with the theoretically predicted taper angles for the optimized shape with F = 0.8, 0.6, 0.4, 0.2, 2mrad, and exponential TOFs. The measured taper angles of the F = 0.4 TOF exceed the theoretical angles at around 30 ≲ r ≲ 40μm. However, because the theoretical angles for the optimized profile coincidentally have much lower values than the targeted angles (0.4Ω) in this region, the measured angles do not considerably exceed the targeted angles.
Fig. 7 Comparison of the predicted and measured taper angles. The cyan, green, red, magenta, blue, and black solid lines are the predicted taper angles for F = 0.8, 0.6, 0.4, 0.2, 2mrad, and 3 mm exponential TOFs, respectively. Red circles and blue squares are measured taper angles for the fabricated F = 0.4 and 2mrad TOFs, respectively. The black dashed lines are FΩ with F = 1.0, 0.8, 0.6, 0.4, 0.2.
4. Conclusion
In conclusion, we have designed and fabricated TOFs with high transmissions and minimal lengths. Transmission in excess of 99.7% has been observed for a TOF with a total length of 23mm. We find that the losses due to the power coupling to higher-order modes are sufficiently suppressed, and the losses are dominated by other extrinsic losses, such as contaminants on the TOF surface. Further improvement of the optimization may be possible by increasing the degrees of freedom for the approximation of Lk at the cost of longer computation time. This work will enable the production of compact, stable setups for coupling whispering-gallery-mode microresonators to TOFs without sacrificing the transmission quality. We plan to build such a setup for the study of cavity quantum electrodynamics with laser-cooled atoms and microtoroidal resonators [29].
Acknowledgments
This work is partially supported by MATSUO FOUNDATION, JSPS KAKENHI Grant Number 26707022, Waseda University Grant for Special Research Projects (Project number: 2013A-501), and a research grant from The Murata Science Foundation.
References and links
1. G. Brambilla, “Optical fibre nanowires and microwires: a review,” J. Opt. A 12, 043001 (2010). [CrossRef]
2. M. J. Morrissey, K. Deasy, M. Frawley, R. Kumar, E. Prel, L. Russell, V. G. Truong, and S. Nic Chormaic, “Spectroscopy, manipulation and trapping of neutral atoms, molecules, and other particles using optical nanofibers: A review,” Sensors 13, 10449–10481 (2013). [CrossRef] [PubMed]
3. T. A. Birks, W. J. Wadsworth, and P. St. J. Russel, “Supercontinuum generation in tapered fibers,” Opt. Lett. 25, 1415–1417 (2000). [CrossRef]
4. S. Leon-Saval, T. Birks, W. Wadsworth, P. St. J. Russell, and M. Mason, “Supercontinuum generation in submicron fibre waveguides,” Opt. Express 12, 2864–2869 (2004). [CrossRef] [PubMed]
5. R. R. Gattass, G. T. Svacha, L. Tong, and E. Mazur, “Supercontinuum generation in submicrometer diameter silica fibers,” Opt. Express 14, 9408–9414 (2006). [CrossRef] [PubMed]
6. D. A. Akimov, A. A. Ivanov, O. A. Kolevatova, M. V. Alfimov, T. A. Birks, W. J. Wadsworth, P. St. J. Russel, A. A. Podshivalov, and A. M. Zheltikov, “Generation of a spectrally asymmetric third harmonic with unamplified 30-fs Cr:forsterite laser pulses in a tapered fiber,” Appl. Phys. B 76, 515–519 (2003). [CrossRef]
7. V. Grubsky and A. Savchenko, “Glass micro-fibers for efficient third harmonic generation,” Opt. Express 13, 6798–6806 (2005). [CrossRef] [PubMed]
8. T. Lee, Y. Jung, C. A. Codemard, M. Ding, N. G. R. Broderick, and G. Brambilla, “Broadband third harmonic generation in tapered silica fibres,” Opt. Express 20, 8503–8511 (2012). [CrossRef] [PubMed]
9. J. Villatoro and D. Monzón-Hernández, “Fast detection of hydrogen with nano fiber tapers coated with ultra thin palladium layers,” Opt. Express 13, 5087–5092 (2005). [CrossRef] [PubMed]
10. P. Polynkin, A. Polynkin, N. Peyghambarian, and M. Mansuripur, “Evanescent field-based optical fiber sensing device for measuring the refractive index of liquids in microfluidic channels,” Opt. Lett. 30, 1273–1275 (2005). [CrossRef] [PubMed]
11. L. Zhang, F. Gu, J. Lou, X. Yin, and L. Tong, “Fast detection of humidity with a subwavelength-diameter fiber taper coated with gelatin film,” Opt. Express 16, 13349–13353 (2008). [CrossRef] [PubMed]
12. F. Warken, E. Vetsch, D. Meschede, M. Sokolowski, and A. Rauschenbeutel, “Ultra-sensitive surface absorption spectroscopy using sub-wavelength diameter optical fibers,” Opt. Express 15, 11952–11958 (2007). [CrossRef] [PubMed]
13. A. Stiebeiner, O. Rehband, R. Garcia-Fernandez, and A. Rauschenbeutel, “Ultra-sensitive fluorescence spectroscopy of isolated surface-adsorbed molecules using an optical nanofiber,” Opt. Express 17, 21704–21711 (2009). [CrossRef] [PubMed]
14. R. Garcia-Fernandez, W. Alt, F. Bruse, C. Dan, K. Karapetyan, O. Rehband, A. Stiebeiner, U. Wiedemann, D. Meschede, and A. Rauschenbeutel, “Optical nanofibers and spectroscopy,” Appl. Phys. B 105, 3–15 (2011). [CrossRef]
15. V. V. Klimov and M. Ducloy, “Spontaneous emission rate of an excited atom placed near a nanofiber,” Phys. Rev. A 69, 013812 (2004). [CrossRef]
16. K. P. Nayak, P. N. Melentiev, M. Morinaga, Fam Le Kien, V. I. Balykin, and K. Hakuta, “Optical nanofiber as an efficient tool for manipulating and probing atomic fluorescence,” Opt. Express 15, 5431–5438 (2007). [CrossRef] [PubMed]
17. E. Vetsch, D. Reitz, G. Sagué, R. Schmidt, S. T. Dawkins, and A. Rauschenbeutel, “Optical interface created by laser-cooled atoms trapped in the evanescent field surrounding an optical nanofiber,” Phys. Rev. Lett. 104, 203603 (2010). [CrossRef] [PubMed]
18. A. Goban, K. S. Choi, D. J. Alton, D. Ding, C. Lacroûte, M. Pototschnig, T. Thiele, N. P. Stern, and H. J. Kimble, “Demonstration of a state-insensitive, compensated nanofiber trap,” Phys. Rev. Lett. 109, 033603 (2012). [CrossRef] [PubMed]
19. M. Fujiwara, K. Toubaru, T. Noda, H. -Q. Zhao, and S. Takeuchi, “Highly efficient coupling of photons from nanoemitters into single-mode optical fibers,” Nano Lett. 11, 4362–4365 (2011). [CrossRef] [PubMed]
20. R. Yalla, Fam Le Kien, M. Morinaga, and K. Hakuta, “Efficient channeling of fluorescence photons from single quantum dots into guided modes of optical nanofiber,” Phys. Rev. Lett. 109, 063602 (2012). [CrossRef] [PubMed]
21. S. Kato, S. Chonan, and T. Aoki, “High-numerical-aperture microlensed tip on an air-clad optical fiber,” Opt. Lett. 39, 773–776 (2014). [CrossRef] [PubMed]
22. S. Chonan, S. Kato, and T. Aoki, “Efficient single-mode photon-coupling device utilizing a nanofiber tip,” Sci. Rep. 4, 4785 (2014). [CrossRef] [PubMed]
23. J. C. Knight, G. Cheung, F. Jacques, and T. A. Birks, “Phase-matched excitation of whispering-gallery-mode resonances by a fiber taper,” Opt. Lett. 22, 1129–1131 (1997). [CrossRef] [PubMed]
24. M. Cai, O. Painter, and K. J. Vahala, “Observation of critical coupling in a fiber taper to a silica-microsphere whispering-gallery mode system,” Phys. Rev. Lett. 85, 74–77 (2000). [CrossRef] [PubMed]
25. T. Aoki, “Fabrication of ultralow-loss tapered optical fibers and microtoroidal resonators,” Jpn. J. Appl. Phys. 49, 118001 (2010). [CrossRef]
26. S. Ravets, J. E. Hoffman, P. R. Kordell, J. D. Wong-Campos, S. L. Rolston, and L. A. Orozco, “Intermodal energy transfer in a tapered optical fiber: optimizing transmission,” J. Opt. Soc. Am. A 30, 2361–2371 (2013). [CrossRef]
27. J. E. Hoffman, S. Ravets, J. A. Grover, P. Solano, P. R. Kordell, J. D. Wong-Campos, L. A. Orozco, and S. L. Rolston, “Ultrahigh transmission optical nanofibers,” AIP Adv. 4, 067124 (2014). [CrossRef]
28. R. Rivière, O. Arcizet, A. Schliesser, and T. J. Kippenberg, “Evanescent straight tapered-fiber coupling of ultra-high Q optomechanical micro-resonators in a low-vibration helium-4 exchange-gas cryostat,” Rev. Sci. Instrum. 84, 043108 (2013). [CrossRef] [PubMed]
29. T. Aoki, B. Dayan, E. Wilcut, W. P. Bowen, A. S. Parkins, T. J. Kippenberg, K. J. Vahala, and H. J. Kimble, “Observation of strong coupling between one atom and a monolithic microresonator,” Nature 443, 671–674 (2006). [CrossRef] [PubMed]
30. J. D. Love and W. M. Henry, “Quantifying loss minimisation in single-mode fibre tapers,” Electron. Lett. 22, 912–914 (1986). [CrossRef]
31. J. D. Love, W. M. Henry, W. J. Stewart, R. J. Black, S. Lacroix, and F. Gonthier, “Tapered single-mode fibres and devices. Part 1: Adiabaticity criteria,” IEE Proc. 138, 343–354 (1991).
32. T. A. Birks and Y. W. Li, “The shape of fiber tapers,” J. Lightwave Technol. 10, 432–438 (1992). [CrossRef]
33. A. Stiebeiner, R. Garcia-Fernandez, and A. Rauschenbeutel, “Design and optimization of broadband tapered optical fibers with a nanofiber waist,” Opt. Express 18, 22677–22685 (2010). [CrossRef] [PubMed]
34. F. Bilodeau, K. O. Hill, S. Faucher, and D. C. Johnson, “Low-loss highly overcoupled fused couplers: fabrication and sensitivity to external pressure,” J. Lightwave Technol. 6, 1476–1482 (1988). [CrossRef]
35. K. Karapetyan, “Single optical microfibre-based modal interferometer,” Ph.D. thesis, Rheinische Friedrich-Wilhelms-Universität Bonn (2012).
36. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C, 3rd ed. (Cambridge University, 2007).
37. F. Orucevic, V. Lefèvre-Seguin, and J. Hare, “Transmittance and near-field characterization of sub-wavelength tapered optical fibers,” Opt. Express 15, 13624–13629 (2007). [CrossRef] [PubMed]
