A generalized heat-brush approach for precise control of the waist profile in fiber tapers Chams Baker∗ and Martin Rochette McGill University, Department of Electrical and Computer Engineering, Montréal (PQ), H3A 2A7, Canada ∗ [email protected] Abstract: We present a generalized heat-brush tapering approach in which the ratio of the feed and draw velocities changes within each tapering sweep. This approach allows for controlled and precise shaping of tapers with an arbitrary waist profile and dissimilar transition regions as demonstrated experimentally. A quantitative analysis of the mismatch error after each tapering sweep is also provided. © 2011 Optical Society of America OCIS codes: (230.4000) Microstructure fabrication; (220.4610) Optical fabrication. References 1. P. Dumais, F. Gonthier, S. Lacroix, J. Bures, A. Villeneuve, P. G. J. Wigley, and G. I. Stegeman, “Enhanced self-phase modulation in tapered fibers,” Opt. Lett. 18, 1996–1998 (1993). 2. C. Baker and M. Rochette, “Highly nonlinear hybrid AsSe-PMMA microtapers,” Opt. Express 18, 12391–12398 (2010). 3. A. C. Boucouvalas and G. Georgiou, “Biconical taper coaxial optical fibre coupler,” Electron. Lett. 21, 864–865 (1985). 4. R. G. Lamont, D. C. Johnson, and K. O. Hill, “Power transfer in fused biconical-taper single-mode fiber couplers: dependence on external refractive index,” Appl. Opt. 24, 327–332 (1985). 5. J. V. Wright, “Wavelength dependence of fused couplers,” Electron. Lett. 22, 320–321 (1986). 6. S. Jensen, “The nonlinear coherent coupler,” J. Quantum Electron. 18, 1580–1583 (1982). 7. J. D. Love and W. M. Henry, “Quantifying loss minimisation in single-mode fibre tapers,” Electron. Lett. 22, 912–914 (1986). 8. J. D. Love, W. M. Henry, W. J. Stewart, R. J. Black, S. Lacroix, and F. Gonthier, “Tapered single-mode fibres and devices. I. adiabaticity criteria,” IEE Proc.-J: Optoelectron. 138, 343–354 (1991). 9. T. A. Birks and Y. W. Li, “The shape of fiber tapers,” J. Lightwave Technol. 10, 432–438 (1992). 10. F. Bilodeau, K. Hill, S. Faucher, and D. Johnson, “Low-loss highly overcoupled fused couplers: Fabrication and sensitivity to external pressure,” in Optical Fiber Sensors, (Optical Society of America, 1988), p. ThCC10. 11. R. P. Kenny, T. A. Birks, and K. P. Oakley, “Control of optical fibre taper shape,” Electron. Lett. 27, 1654–1656 (1991). 12. Y. Takeuchi and J. Noda, “Novel fiber coupler tapering process using a microheater,” Photon. Technol. Lett. 4, 465–467 (1992). 13. H. Yokota, E. Sugai, Y. Kashima, and Y. Sasaki, “Optical irradiation method for fiber coupler fabrications,” in Optical Fiber Sensors, (OSA, 1996), p. Th319. 14. A. J. C. Grellier, N. K. Zayer, and C. N. Pannell, “Heat transfer modelling in co laser processing of optical fibres,” Opt. Commun. 152, 324–328 (1998). 15. L. C. Ozcan, V. Treanton, F. Guay, and R. Kashyap, “Highly symmetric optical fiber tapers fabricated with a CO2 laser,” Photon. Technol. Lett. 19, 656–658 (2007). 16. S. Xue, M. van Eijkelenborg, G. W. Barton, and P. Hambley, “Theoretical, numerical, and experimental analysis of optical fiber tapering,” J. Lightwave Technol. 25, 1169–1176 (2007). 17. J. Dewynne, J. R. Ockendon, and P. Wilmott, “On a mathematical model for fiber tapering,” SIAM J. Appl. Math. 49, 983–990 (1989). 18. S. Pricking and H. Giessen, “Tapering fibers with complex shape,” Opt. Express 18, 3426–3437 (2010). #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1065 19. W. Sun, M. Yuan, X. Zeng, and T. Birks, “Theoretical shape analysis of tapered fibers using a movable large-zone furnace,” Optoelectron. Lett. 7, 154–157 (2011). 20. F. Geyling, “Basic fluid dynamic consideration in the drawing of optical fibers,” Bell Sys. Tech. J. 55, 1011–1056 (1976). 21. S. Leon-Saval, T. Birks, W. Wadsworth, P. S. J. Russell, and M. Mason, “Supercontinuum generation in submicron fibre waveguides,” Opt. Express 12, 2864–2869 (2004). 22. A. Kudlinski, A. K. George, J. C. Knight, J. C. Travers, A. B. Rulkov, S. V. Popov, and J. R. Taylor, “Zerodispersion wavelength decreasing photonic crystal fibers for ultraviolet-extended supercontinuum generation,” Opt. Express 14, 5715–5722 (2006). 23. G. Qin, X. Yan, C. Kito, M. Liao, T. Suzuki, A. Mori, and Y. Ohishi, “Zero-dispersion-wavelength-decreasing tellurite microstructured fiber for wide and flattened supercontinuum generation,” Opt. Lett. 35, 136–138 (2010). 24. A. C. Judge, O. Bang, B. J. Eggleton, B. T. Kuhlmey, E. C. Mägi, R. Pant, and C. M. de Sterke, “Optimization of the soliton self-frequency shift in a tapered photonic crystal fiber,” J. Opt. Soc. Am. B 26, 2064–2071 (2009). 25. A. Alkadery and M. Rochette, “Widely tunable soliton shifting for mid-infrared applications,” in IEEE Photonics Conference 2011 (IPC2011) (2011). 26. N. Vukovic, N. Broderick, M. Petrovich, and G. Brambilla, “Novel method for the fabrication of long optical fiber tapers,” Photon. Technol. Lett. 20, 1264–1266 (2008). 27. S. Chapra and R. Canale, Numerical Methods for Engineers, 5th ed. (McGraw-Hill, 2005). 28. A. W. Snyder, “Coupled-mode theory for optical fibers,” J. Opt. Soc. Am. 62, 1267–1277 (1972). 29. T. H. Cormen, C. E. Leiserson, R. L. Rivest, and C. Stein, Introduction to Algorithms, 2nd ed. (MIT Press, 2001). 1. Introduction Tapered optical fibers, illustrated in Fig. 1, made by a heat-and-draw approach are of interest in a wide range of applications. They have been used for enhancing nonlinear effects [1, 2], coaxial mode coupling [3], power splitting/combining [4], filtering optical spectra [5], and switching [6]. In all cases, a fine control of the taper shape is required to ensure an adiabatic transformation of the propagating mode [7, 8]. Fig. 1. Schematic of a fiber taper with a uniform waist and similar transition regions. A tapering model presented by Birks [9] provides an approach for shaping a fiber taper by changing the hot-zone length as the fiber is symmetrically stretched under tensile force at both ends. Birks’ model can be implemented using a stationary heater with a variable-length hot-zone, or using a heat-brush approach [10], where a heater travels back and forth within a variable-length brushing-zone [9]. The heat-brush implementation of Birks’ model provides better precision in shaping fiber tapers than the stationary heater implementation [11]. The heater in the heat-brush implementation can be a flame [10], a resistive heater [12], or a CO2 Laser [13–15]. The stationary heater implementation of Birks’ model has been analyzed theoretically and numerically [16] using a viscus fluid flow model [17]. There has also been a few heuristic theoretical and numerical analyses of taper shape evolution in the heat-brush implementation of Birks’ model [18, 19]. Precise taper shape evolution after each heater sweep can be analyzed using a viscus fluid flow model generally used for the study of fiber drawing [17, 20]. The effects of different parameters, such as the number of tapering sweeps and the hot-zone length, on the taper shape can thus be quantified. #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1066 Moreover, in the heat-brush implementation of Birks’ model, the tapering function s = v f /vd , each tapering where v f is the feed velocity and vd is the draw velocity, is constant throughout √ sweep. A constant s limits the lowest inverse tapering ratio ρ = φ j /φ j−1 = s, where φ j is the waist diameter after sweep j, that can be used in each sweep [21]. If ρ is less than 0.97 [11], the taper diameter in the transition region does not change smoothly, but rather it changes in steps. A generalized heat-brush method allows s to change as the heater sweeps along the brushingzone, and hence, the taper shape is carved within each sweep rather than having a sudden change in diameter. Just as in the heat-brush approach, the generalized heat-brush approach allows for precise shaping of the transition regions [11], a uniform waist profile [21], and a large contrast ratio between the initial and the final taper diameters [10]. In addition, the generalized approach allows for a smaller ρ in each sweep as well as controlled fabrication of tapers with an arbitrary waist profile and dissimilar transition regions [18]. A smaller ρ in each sweep reduces the number of sweeps required in the tapering process, and hence, reduces the taper fabrication duration. Replacing the uniform waist profile by one that follows an arbitrary function provides additional freedom in taper design and widens the range of taper applications. For example, a nonuniform waist profile in tapered fibers shifts the zero-dispersion wavelength along the microtaper waist for extended and flat supercontinuum generation [22, 23] and enhanced soliton self-frequency shift [24, 25]. Dissimilar transition regions also provide additional freedom in taper design. For example, in the case of soliton self-frequency shifting due to the Raman effect, the spectrum of a soliton slides towards longer wavelengths as it propagates from the input end to the output end of a taper waist. A design that minimizes the length of the taper has dissimilar adiabatic transition regions [8]. In this paper, we develop and demonstrate both by simulation and experiment a generalized heat-brush tapering method, and use it for the fabrication of tapers with a nonuniform waist profile and dissimilar transition regions. First, single-sweep tapering, the main constituent of the generalized heat-brush approach, is presented and simulated using a viscous fluid flow model to quantify the mismatch error between the targeted and the resulting taper profiles. Then, the generalized heat-brush approach is implemented by tapering a fiber over multiple sweeps, and the simulation results from the single-sweep tapering analysis are used to quantify the accumulated mismatch error after each tapernig sweep. Finally, we use of the generalized heat-brush approach to fabricate an As2 Se3 chalcogenide taper with a linearly decreasing waist profile and dissimilar transition regions. 2. Single-sweep tapering In this section, we present the single-sweep tapering method, an instance of the well-known fiber-drawing approach [17, 20, 26]. In the process of fiber drawing, mass conservation leads to φ (t) = φ0 s (t) where φ (t) is the taper diameter, φ0 is the initial fiber diameter, and s (t) = v f (t) /vd (t) is the tapering function. To draw a taper with a predefined profile φ (z), the tapering function s (t) must be determined accordingly. The replacement of the time variable t by the ´t drawing length ld (t) = 0 vd (τ ) d τ simplifies the implementation of the single-sweep tapering method because it can be readily used as a feedback parameter to control the draw velocity vd (ld ) = v f (ld ) /s (ld ). In this case, the tapering function s (ld ) is calculated from the taper profile φ (z) using φ 2 (z) s (ld ) = . φ2 0 z=ld Figure 2 provides an arbitrary taper profile φ (z) and its corresponding tapering function s (ld ). #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1067 1 φ μ 150 100 0.5 50 (b) (a) 0 0 10 20 30 z [mm] 0 0 10 20 30 ld [mm] Fig. 2. (a) A taper profile and (b) the resulting tapering function. 2.1. Single-sweep tapering simulations A general model of the viscous flow in the heat-softened region, or hot-zone, due to unidirectional stretching has been reported in [20]. A simplified model has been derived for the case when the fiber diameter is much smaller than the hot-zone length (Lhz ) [17]. In this model, the deformation of the hot-zone due to stretching is governed by ∂ ∂u 3μ A = 0, (1) ∂z ∂z ∂A ∂ + (uA) = 0, (2) ∂t ∂z where μ (z,t) is the viscosity distribution, u (z,t) is the axial velocity distribution, and A (z,t) is the cross-sectional area in the hot-zone [17]. For a Newtonian fluid, μ is independent of u, ¯ ∂ z × ∂ F/∂ z + F × ∂ 2 u/ ¯ ∂ z2 = 0, where u¯ = u/vd is the norand hence, Eq. (1) leads to ∂ u/ malized axial velocity and F = μ A. Using the centered differentiation formulas [27] ∂ F/∂ z = ¯ ∂ z = (u¯i+1 − u¯i−1 ) /2Δz, and ∂ 2 u/ ¯ ∂ z2 = (u¯i+1 − 2u¯i + u¯i−1 ) /Δz2 leads (Fi+1 − Fi−1 ) /2Δz, ∂ u/ to the finite difference form of Eq. (1) [Fi − 0.25 (Fi+1 − Fi−1 )] u¯i−1 − 2Fi u¯i + [Fi + 0.25 (Fi+1 − Fi−1 )] u¯i+1 = 0 (3) where Fi = F (ld , zi ), u¯i = u¯ (ld , zi ), and Δz is the separation between any two consecutive zi . Changing the variable t to ld in Eq. (2) leads to the equation vd ∂ A/∂ ld + ∂ (uA) /∂ z = 0, which ¯ ∂ z + u¯∂ A/∂ z = 0. Using the cenis expanded and divided by vd to obtain ∂ A/∂ ld + A∂ u/ ¯ ∂ z = (u¯i+1 − u¯i−1 ) /2Δz, ∂ A/∂ z = (Ai+1 − Ai−1 ) /2Δz and tered differentiation formulas ∂ u/ − Ai ] /Δld , the finite difference form the forward differentiation formula [27] ∂ A/∂ ld = [Anew i of Eq. (2) corresponding to the extension of the fiber by a distance Δld = 2Δz is given by = Ai − [Ai (u¯i+1 − u¯i−1 ) + u¯i (Ai+1 − Ai−1 )] Anew i (4) = A (ld + Δld , zi ). It is clear from Eq. (3) and Eq. (4) that, for a where Ai = A (ld , zi ), Anew i Newtonian fluid, the deformation of the hot-zone is independent of the actual drawing velocity. The flow-chart in Fig. 3 describes the program used to simulate the single-sweep experimental setup presented in Section 2.3. In this program, the taper profile is represented by an array of diameter values φk taken at points zk with any two consecutive points separated by Δz. The hot-zone is a subarray of the taper array and the starting point of the hot-zone subarray can change to simulate a moving heater as illustrated in Fig. 4(a). The cross-section area in the hot-zone is given by Ai where i = 1, 2, ..., N and the cross-section area of the extended hot-zone #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1068 that results from drawing the hot-zone, as illustrated in Fig. 4(b), is calculated as follows: first, Eq. (3) is used with the boundary conditions u¯i=0 = −1/2 and u¯i=N+1 = 1/2 to calculate the normalized axial velocity distribution u¯i in the hot-zone, and then, Eq. (4) is used to calculate the extended hot-zone profile. In the simulations that follow, the hot-zone is assumed to have a uniform viscosity distribution. Fig. 3. Flow-chart of the simulation program for the single-sweep tapering setup presented in Section 2.3. In this flow-chart, x is the displacement of both translation stages extending the fiber, y is the displacement of the heater translation stage, x previous and y previous are state variables, δ is a differential feed step, s is the tapering function, ld is the drawing length, and Δz is the longitudinal separation between any two consecutive diameter sampling points. Fig. 4. Single-sweep simulation schematics of (a) shifting the hot-zone by Δz, and (b) extension of the fiber by 2Δz. We simulate the fabrication of a step-taper where the diameter changes abruptly from the initial to the final taper diameter. Typical simulation results of step-taper fabrication show a transient response in the resulting taper with an overshoot and oscillations in the waist before the diameter settles to a final value, as shown in Fig. 5. The mismatch between the resulting and the targeted taper profiles is quantified by the percent error along the taper defined as ε (z) = [φr (z) − φt (z)] /φt (z) × 100% where φr is the resulting taper diameter and φt is the targeted taper diameter. The transient response is quantified by the percent overshoot εos = #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1069 (φt − φos ) /φt ×100% where φos is the overshoot diameter, and by the settling distance zs defined as the distance between the beginning of the waist and the point where the envelope of the absolute percent error is less than εs = 2%. Design Simulation Fig. 5. Simulation of step-taper fabrication using the single-sweep tapering method. The transient response parameters εos and zs represent the closeness of the of the resulting taper shape to the taper design, and the overall mismatch is reduced by reducing εos and zs . Steptaper simulation results in Fig. 6 show εos and zs as a function of Lhz and the inverse tapering ratio ρ = φmin /φ0 , where φmin is the minimum taper diameter. As expected, εos and zs decrease with increasing ρ (≤ 1) and shortening Lhz . With respect to optical propagation in the taper, the overshoot in the waist diameter acts as a perturbation that may lead to coupling between the fundamental mode and higher order modes, radiation modes, or reflection modes [28]. The values of εos and zs also represents the strength and the length of the perturbation region; therefore, a lower εos and a shorter zs reduces the perturbation impact. 30 10 Lhz = 2.0 mm (b) Lhz = 3.0 mm 30 Lhz = 3.0 mm Lhz = 4.0 mm 20 Lhz = 4.0 mm s εos 20 0 0.4 40 Lhz = 2.0 mm (a) 10 0.6 ρ 0.8 1 0 0.4 0.6 ρ 0.8 1 Fig. 6. (a) Percent overshoot, and (b) settling distance dependence on the inverse tapering ratio at different hot-zone lengths for the step-taper. 2.2. Single-sweep tapering optimization Simulation results in subsection 2.1 showed that εos and zs decrease when ρ → 1 and Lhz → 0 mm. However, applications such as the enhancement of the waveguide nonlinearity or the sensitivity require microtapers with a waist diameter on the order of 1 µm drawn from fibers with a diameter on the order of 100 µm leading to ρ ∼ 0.01. Also, Lhz is on the order of 1 mm and is limited by the temperature distribution in the fiber and the heater dimensions. Moreover, it turns out that εos and zs decrease when the taper slope decreases. As an example, Fig. 7 shows that as the slope decreases from 0.0105 to 0.0035, εos decreases from 8.8% to 3.8% #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1070 and zs decreases from 13.5 mm to 11.65 mm. In most cases, however, it is desirable to use the largest slope allowed by the adiabaticity criteria because using a small taper slope to reduce εos and zs leads to a long transition region and consequently increases the sensitivity of the taper to environmental variations [9] as well as increasing the device length. Section 3 shows that εos and zs are reduced by tapering a fiber over multiple sweeps leading to an implementation of the generalized heat-brush approach. φ φmax 0.8 0.6 0.8 0.6 (a) 0.4 Design Simulation |dφ / dz| = 0.0035 1 φ φmax Design Simulation |dφ / dz| = 0.0105 1 (b) -20 0 z [mm] 20 0.4 -20 0 20 z [mm] Fig. 7. Simulated fabrication results of taper profiles with linear transition regions at different slopes using the single-sweep tapering method. 2.3. Experimental setup Figure 8 illustrates the experimental implementation of the single-sweep tapering method where a translation stage moves the heater at a velocity vy and two other translation stages pull the fiber from opposite directions at equal velocities vw and vx . Using vd = vy + vw and v f = vy − vx = α , where α is a constant, the velocities of the heater and the translation stages pulling on the fiber at a drawing length ld = y + w are vd (ld ) + v f (ld ) α 1 = +1 , vy (ld ) = 2 2 s (ld ) vd (ld ) − v f (ld ) α 1 = −1 . vx (ld ) = vw (ld ) = 2 2 s (ld ) Fig. 8. Schematic of the experimental implementation of the single-sweep tapering method. #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1071 2.4. Single-sweep tapering experimental results Figure 9(a) shows the experimental results of a step-taper fabricated from an As2 Se3 fiber ◦ with an initial diameter of 170 µm using a 5 mm long resistive heater at 210 C with v f = max 0.72 mm/min and vd = max v f /s = 4.5 mm/min. The fabricated taper is removed from the tapering setup and placed straight on a flat plate, and then, an imaging system composed of a 20× lens and a CCD camera mounted on a motorized translation stage is used to measure the taper profile with a measurement taken every 1.0 mm. The measured step-taper profile clearly shows an overshoot in the fiber diameter arising from the finite length of the hot-zone. An effective hot-zone length of 2.7 mm is retrieved by simulating the step-taper fabrication and fitting the simulation results with the measured profile. The measured effective length is used to simulate the fabrication of the taper in Fig. 9(b) and the simulation results show good agreement with the experimental results within the measurement error of 1 µm. Design Simulation Measurements 0.8 0.6 1.0 φ φmax φ φmax 1 (a) 0.4 Design Simulation Measurements 0.8 0.6 (b) 0.4 -20 -10 0 10 20 -20 z [mm] 0 20 z [mm] Fig. 9. Experimentally measured profiles of (a) a step taper, and (b) an arbitrary taper fabricated using the single-sweep tapering method. 3. Multi-sweep tapering Multi-sweep tapering performed by systematic repetition of the single-sweep method as illustrated in Fig. 10 represents an implementation of the generalized heat-brush method. To taper a fiber over n sweeps, the taper profile is divided into subsections as shown in Fig. 11, where φn is the minimum taper diameter, and φ1 to φn−1 are the waist diameters for all intermediate tapering sweeps and are calculated using φ j = rφ j−1 with r = ρ 1/n and ρ = φn /φ0 . For every sweep j < n, the stage tapering function s( j) (l p ) is calculated from the stage taper profile ft φ ( j) (z) composed of a left transition region extracted from φ (z) between zlej−1 and zlej f t , a right transition region extracted from φ (z) between zright and zright j j−1 , and a uniform waist with a length ´ zright j le f t Lj = zj φ 2 (z) dz φ j2 , where L j makes the mass volume of the waist at stage j equal to the mass volume required to . The stage taper profile of the final sweep φ (n) (z) draw the taper section between zlej f t and zright j ft right is extracted from φ (z) between zle n−1 and zn−1 , and is used to calculate the final stage tapering (n) function s (l p ). Finally, for each stage j, a single tapering sweep is performed using thecalcu- right + L j. lated stage tapering function and then the heater is moved back a distance zright j−1 − z j #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1072 Fig. 10. Schematic of taper profile evolution using the multi-sweep tapering method. Fig. 11. Dividing the taper into sections for the determination of the tapering function of each tapering stage. 3.1. Quantitative analysis of multi-sweep tapering Based on the divide-and-conquer paradigm [29], tapering a fiber over multiple sweeps reduces the percent overshoot. For a step-taper, the worst-case overshoot diameter at sweep j is estimated using the recurrence relation ( j) ( j−1) φos = [1 − εos (ρ j ) /100%] × ρ j × φos (1) φos = [1 − εos (ρ1 ) /100%] × ρ1 × φ0 , #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1073 where εos (ρ j ) is provided in Fig. 6(a). By setting the inverse tapering ratio for all sweeps to r, the worst-case overshoot diameter becomes ( j) φos = [1 − εos (r) /100%] j × r j × φ0 , and the maximum percent overshoot at the end of tapering is (n) εos,max = [1 − (1 − εos (r) /100%)n ] × 100%, (5) (n) which is simplified to εos,max ≈ nεos (r) when εos (r) ≤ 1%. It is clear from Fig. 6(a) that (n) (n) εos,max < εos (ρ ) and that εos,max decreases as n increases. For, example, the fabrication of a steptaper with ρ = 0.5 over a single sweep using a 4 mm long hot-zone leads to εos (0.5) = 17%. However, when tapering is performed over 6 sweeps with r = 0.89 and εos (0.89) = 0.5%, the (6) maximum percent overshoot is εos,max = 3%. The use of a large number of sweeps increases the tapering duration. For the case of a stepmax is the maximum taper, the minimum time duration for stage j is T j = L j−1 /vmax f , where v f practical feed velocity, and the total tapering duration after n sweeps is T= 1 − ρ −2 L0 × , vmax 1 − ρ −2/n f and reducing n. In general, to keep the tapering duration which is reduced by increasing vmax f at a minimum, n is selected to be the minimum number of sweeps required to keep εos bellow a certain prescribed value. 3.2. Reduced mismatch in the transition regions using multi-sweep tapering The diameter decreases in steps in the heat-brush implementation of Birks’ model, limiting the minimum attainable mismatch between the resulting taper and the design. At any diameter φ , the diameter step is Δφ = (1 − ρ ) φ and the taper slope is approximated by ∂ φ /∂ z ≈ Δφ /Δz leading to Δz ≈ (1 − ρ ) φ / (∂ φ /∂ z). Setting Lhz |Δz| does not decrease the mismatch because the diameter steps in the transition region become more prominent; in fact, setting Lhz |Δz| is practical to keep the transition region smooth. For example, if the length of the brushing-zone is a constant L0 , then the taper profile is given by φ (z) = φ0 exp (−z/L0 ) [9] and |Δz| ≈ (1 − ρ ) L0 . Using typical values of ρ = 0.97 and L0 = 2.0 cm leads to |Δz| ≈ 0.6 mm, which requires Lhz 0.6 mm. In contrast, the diameter steps are eliminated in the multi-sweep tapering method because the transition region is carved within each tapering sweep; therefore, shortening Lhz always reduces the mismatch between the resulting taper and the design. 3.3. Multi-sweep tapering simulation Multi-sweep tapering simulation is performed by repeated application of the single sweep tapering program. Simulation results in Fig. 12 performed using Lhz = 3 mm for a step-taper with (n) ρ = 0.4 show that the percent overshoot εos decreases as n increases. Also shown in Fig. 12 is (n) (n) the worst-case percent overshoot, εos,max , calculated using Eq. (5). It is observed that εos does (n) (n) (n) not exceed εos,max , which is expected as εos,max estimates the upper limit of εos . #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1074 15 ρ = 0.4 εos(n) Lhz = 3.0 mm (n) εos,max εos 10 5 2 3 4 n 5 6 7 Fig. 12. Percent overshoot and maximum percent overshoot versus the number of tapering sweeps n for a step-taper with ρ = 0.4 using Lhz = 3 mm. Although increasing n reduces εos , Lhz must also be shortened to ensure that |ε (z)| is less than a prescribed value εtarget . Shortening Lhz is critical when the taper profile incorporates fine details such as a large ∂ φ /∂ z, a large change in ∂ φ /∂ z, or a short waist. For example, if the taper waist length is of the same order as Lhz , then the details of the waist can not be precisely shaped. The value of Lhz that ensures |ε (z)| < εtarget for a given taper profile can be determined through simulations. 3.4. Multi-sweep tapering experimental results Figure 13 shows the experimental results for the fabrication of an As2 Se3 taper with an initial fiber diameter of 170 µm, dissimilar left and right transition regions, and a nonuniform waist with a diameter decreasing linearly from 15 µm to 10 µm over a waist length of 2.0 cm. The taper is experimentally fabricated over 24 sweeps using the same resistive heater in the single-sweep experiment in Subsection 2.4 at 210◦ C with v f = 3.56 mm/min = 4.50 mm/min. The measurement error is 1 µm and the resulting taper matches the and vmax d design within the measurement error. Design Simulation Measurements φ φmax 1 0.5 0 -50 -30 -10 10 z [mm] 30 50 Fig. 13. Experimental results showing the profile of an As2 Se3 taper fabricated using the multi-sweep tapering method with n = 24. #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1075 4. Conclusion The multi-sweep tapering method has been used to implement the generalized heat-brush approach, which allows the ratio of the feed and draw velocities to change within each tapering sweep. A quantitative analysis showed that the mismatch error decreases by increasing the number of tapering sweeps and shortening the length of the hot-zone formed by the heater. An As2 Se3 chalcogenide taper with dissimilar transition regions and a waist diameter decreasing linearly from 15 µm to 10 µm over 2.0 cm was fabricated using the multi-sweep tapering method showing good agreement between the targeted and the measured taper profiles. #151542 - $15.00 USD (C) 2011 OSA Received 22 Jul 2011; revised 18 Aug 2011; accepted 29 Aug 2011; published 2 Sep 2011 1 October 2011 / Vol. 1, No. 6 / OPTICAL MATERIALS EXPRESS 1076
© Copyright 2024