Abstract
Propagation of whispering gallery modes in a microfiber coil having an n-fold helical symmetry is considered. The n-fold helically symmetric deformations of a coiled microfiber is similar to a long period grating introduced by periodic microfiber bending with the azimuthal angle period 2π/n. The considered modes are localized near a geodesic situated at the peripheral part of the microfiber surface. Using the perturbation theory, a simple condition for the stability of this geodesic is found. Violation of this condition happens in a small region localized near a point determined by the phase matching condition for the introduced long period grating. In the classical limit, the unstable region corresponds to the parametric resonance and stochastization of whispering gallery rays. Generally, this region corresponds to resonance intermode coupling, which may result either in the periodic transmission of radiation power between two modes or in the aperiodic wave parametric resonance, which causes simultaneous coupling and power transfer between numerous whispering gallery modes. The obtained results are important for engineering of miniature optical fiber coils and analysis of their propagation loss.
©2010 Optical Society of America
1. Introduction
There has been an increasing interest in development of theory, fabrication methods, and applications of optical fiber and microfiber coils [1–15]. Coiling an optical fiber results in significant reduction of its dimensions, which is important for applications. Not less importantly, a coiled waveguide exhibits remarkable physical properties. Initially, special interest in these devices was caused by a class of problems related to the polarization of light propagating along a fiber coil, e.g., to the geometric rotation of polarization [1–3,6]. Later, a new type of miniature optical resonator composed of tightly coiled optical microfiber with coupling turns was considered [4,8–12]. Recently, a miniature multimode fiber coil (MFC) with uncoupled turns illustrated in Fig. 1(a) was suggested to be used as a compact optical delay line having small losses, small dimensions, and, at the same time, broadband andsmooth transmission characteristics [14]. The performance of this delay line is based on the fundamental whispering gallery mode (WGM) propagating along the periphery of the coiled microfiber. Besides, a similar device can be used as a very perceptive evanescent optical sensor based on the single fundamental WGM propagation [13] and multi-WGM interference [15]. The advantage of this type of optical sensor is in its miniature dimensions, which are significantly smaller than the length of the microfiber waveguide.
In practice, the shape of an MFC is never perfectly uniform. The MFC nonuniformities are introduced in the process of microfiber and central rod fabrication as well as in the process of wounding of an MFC. Special MFC nonuniformities can be introduced on purpose in order to provide certain device functionalities. Understanding of the effect of MFC nonuniformities is critical both for the experimental realization of a low-loss MFC and for engineering of advanced nonuniform MFC devices.
Consider an MFC created of a microfiber having a diameter of a few tens of microns. This MFC supports a series of WGMs propagating near the peripheral part of the coiled microfiber. Ideal MFC, which consists of a uniform microfiber coiled on an axially symmetric and uniform rod, possesses the infinite helical symmetry and does not exhibit intermode transitions. Deviation from this symmetry causes coupling between modes. This paper addresses the problem of propagation of light along the deformed MFC having an n-fold helical symmetry, which can be fabricated by coiling a microfiber on a central rod having an n-fold axial symmetry as illustrated in Fig. 1(b) for . The effect of the introduced periodic bending of the microfiber is similar to the effect of a long period grating (LPG) in a straight optical waveguide [16–18]. Due to the special properties of the WGMs (e.g., the quasi-equal spacing between propagation constants) the introduced perturbation is responsible for a quite complex dynamics of WGMs investigated below. Beside the general interest in engineering of the MFCs, the obtained results are useful for understanding and evaluation of the MFC transmission loss.
In Section 2, the structure and analytical expressions of the WGMs propagating along the nonuniform MFC are described. In Section 3, the classical perturbation theory for the whispering gallery rays propagating along the peripheral part of the coiled microfiber is developed and the condition of stability of these rays is found. In Section 4, it is shown that breaking of the condition of stability coincides with the condition of the classical parametric resonance and, simultaneously, with the phase matching condition of the introduced LPG. In Section 5, a numerical example of the classical ray stochastization in the neighborhood of parametric resonance is given. In Section 6, the behavior of WGMs in the neighborhood of parametric resonance is demonstrated numerically using the beam propagation method and analyzed using the coupled wave equations. It is shown that, in certain cases, there exists a periodic power transfer between two modes, while, in other cases, the wave parametric resonance causes complex simultaneous and aperiodic multi-mode coupling. Finally, in Section 7, the results of the paper are discussed and summarized.
2. Whispering gallery modes (WGMs) in a microfiber coil (MFC)
Propagation of WGMs in a uniform MFC was considered in Ref [14]. Here it is assumed that the nonuniformities of the MFC are relatively smooth having a characteristic length much greater than the radiation wavelength. In addition, the mode behavior near a thin layer surrounding the MFC surface is approximated by the zero boundary condition at the surface. Then, the propagation of a WGM localized near a bent surface can be described through the parameters of classical geodesics on this surface [19]. The surface geodesics serve as “classical bones” near which a WGM “wave flesh” can be constructed [20]. The geodesic of our interest is the helical trajectory s shown in Fig. 1(c). Introduce the local orthogonal coordinate system near this geodesic where s is the coordinate along the geodesic and p and q are, respectively, the coordinates along the normal and binormal vector of this geodesic at point s. The asymptotic expression for the modes (or their linear combinations) localized near geodesic s has the form [19]:
Here, is the mode propagation constant, which is determined asIn Eq. (1) and (2), is the local radius of curvature of the geodesic, is the microfiber refractive index, λ is the radiation wavelength, is the Hermite polynomial, is the Airy function, is the root of Airy function (…) andwhere functions satisfy the linear differential equation:In Eq. (4), is the Gaussian curvature [21] of the MFC surface at the geodesic s and coefficients are free parameters connected by the equality:As a particular case, consider a uniform MCR with microfiber radius and central rod radius . The radius of curvature of the geodesic in this case is constant:and solutions of Eq. (4) are where the Gaussian curvatureGenerally, Eqs. (1) and (2) describe a linear combination of WGMs. Individual WGMs can be obtained by setting , so that, from Eq. (5), and, from Eq. (3), . Substitution of these values into Eqs. (1) and (2) yields the expression for the WGMs considered in Ref [14]:From Eq. (1), the characteristic widths of the WGMs along axes x and y are and , respectively. For the uniform MFC, these widths are and . WGMs which are localized near the geodesic for arbitrarily large s and experience only small changes for a weakly perturbed MFC are called the stable WGMs. In other words, we call a WGM stable if (a) its characteristic widths are bounded for any s and (b) if any small s-dependent perturbation of the MFC parameters does not cause unlimited growth of these widths. The condition of stability along the n-axis is simple: it holds if the radius of curvature is bounded. The condition of stability along the t-axis is less trivial: it coincides with the condition of stability of the geodesic s at the MFC surface, which is determined by the boundedness of the solution of Eq. (4). This differential equation is the Jacobi equation, which determines the geodesics at the MFC surface adjacent to the geodesic s [22].Equations (1)-(5) define the evolution of WGMs along the MFC with varying parameters. It allows to obtain simple analytical expressions for the WGMs and for the transmission amplitudes between WGMs provided that the shape of the geodesic s, its curvature, , and the Gaussian curvature of the MFC surface at this geodesic,, are expressed through the MFC parameters. A model of a perturbed MFC is considered in the next Section.
3. Stable geodesic and the Gaussian curvature
In this section, an MFC surface is modeled in the following simplified form:
where is the local radius of the central rod and is the local radius of the coiled microfiber. For simplicity, the local asymmetry of microfiber, i.e., the dependence of r on θ is ignored. For a uniform MFC, and are constants and the coil pitch is . Samples of MFCs depicted with Eq. (9) are shown in Fig. 1.The shape of an MFC is assumed to be close to a uniform MFC, so that
The WGMs of our interest are localized in the neighborhood of the geodesic s located at the peripheral part of the MFC surface as shown in Fig. 1(c). For a uniform MFC, s is a helix defined by the equations:For the relatively weak dependence of and on azimuthal angle φ determined by Eq. (10), the shape of the trajectory s is obtained from Eq. (9) by substitution where is a solution of the nonuniform Jacobi equation:Here s is the coordinate along the unperturbed helical trajectory, , so that and .Below, only MFCs with a relatively small pitch will be considered. It is seen from Eqs. (10) and (12) that, in the approximation linear in , , and , the shape of the trajectory s is not changed.
As it is explained in Section 2, the propagating WGMs are localized only near a geodesic, which is stable on the MFC surface. To check for stability, it is sufficient to determine how the adjacent geodesics spread away from s. The deviation of these geodesics is determined by the Jacobi Eq. (4) where, in the linear in, , and γ approximation, the Gaussian curvature is reduced to
4. Long period grating: mode coupling and the classical parametric resonance
Assume that the microfiber radius is uniform, , while the central rod is axially asymmetric having the n-fold rotational symmetry:
In this case, the surface of the coiled microfiber, which supports the geodesic s, experiences periodic bending with the period . This bending is a particular case of a long period grating explicitly investigated in fiber optics [16–18]. It is known that the long period grating, which period matches the difference between two propagation constants of two modes so that , causes resonant coupling between these modes. For a uniform MFC, the propagation constants are approximately uniformly spaced in quantum number m (see Eq. (8)) and the matching condition for these propagation constants causes multiple coupling between modes. From Eqs. (8) and (14), due to the symmetry of perturbation, only coupling between the adjacent modes with the same parity is not zero. For these modes, the expression for the propagation constant in Eq. (8) gives and the phase matching condition reads:In classical mechanics, the latter effect is the well known parametric resonance, which, in our case, is described by the Jacobi Eq. (4). In fact, with Eq. (13), Eq. (4) takes the form of the Mathieu equation:In the cylindrically symmetric case, , Eq. (16) shows that the geodesics adjacent to the helical trajectory s experience oscillations with period as a function of φ (or with period as a function of s). For , the classical parametric resonance [23] occurs when the period of harmonic perturbation, , is close to a half of period , i.e., , . Substitution of the values of periods into the latter equation yields the condition of the classical parametric resonance for a MFC:For , this equation coincides with the LPG phase matching condition, Eq. (15) that will be further considered in Section 6. Clearly, the rotational orders are excluded from this equation because . More detailed calculation [24] shows that the solutions of Eq. (16) are stable forand are unstable otherwise. In particular, for the microfiber radius close to the parametric resonance value, , , Eqs. (17) - (19) yield the following condition of stability:Usually, so that and Eq. (19) is reduced to .5. Ray dynamics: stability and the classical parametric resonance
It is interesting to investigate the evolution of optical rays, which propagate inside the MFC near the geodesic s and possess finite but small angle of reflection from the MFC surface. The problem of 3D evolution of an optical ray is solved as follows. Parameter θ is eliminated from Eq. (9) of the MFC surface and the latter is represented in the cylindrical coordinatesas. The equations of a straight optical ray are represented parametrically as a function of φ in the form , , , where and, , , , , are constants. Then, the point of intersection between the MFC surface and aray is found by numerical solution of equation . This approach allows us to perform much faster numerical calculations compared to a more general method, in which determination of the intersection between a straight line and a surface requests solution of two equations with two unknowns (in our case, θ and φ).
Here, the model of an MFC having the n-fold rotational symmetry with is considered. It is assumed that the parameters of an MFC are close to the parametric resonance. The radius of the central rod is set to μm (the length units considered in this Section are not important and can be rescaled). Then, from Eq. (17), we have μm. The MFC pitch is put equal to . The MFC coil with μm, μm, and μm is shown in Fig. 2 . The amplitude of the central rod radius perturbation in Eq. (14) is set to . The point plots in Figs. 3(a)-(g) show fragments of the Poincaré surfaces of sections, i.e., the sine of the reflection angle at reflection points, (left hand side figures) and the relative deviation of the reflection points (right hand side figures) versus the number of turns. The fragments are calculated for the ray (red line in the inset of Fig. 2) launched parallel to the geodesic s from the point at axis n, which is spaced by 0.07 from s. The spacing 0.07 is chosen to arrive at the unstable and stochastic behavior of the ray for μm shown in Fig. 3(d). In our numerical simulations, the evolution of the similarly defined rays was investigated for different microfiber radii . Comparison of this behavior for the MFC with microfiber radii μm, μm, μm, and μm plotted in Fig. 3(a)-(d) shows that there is a sharp lower threshold of instability, so that for μm the ray behavior is stable (for the case of Fig. 3(a), simulation beyond the plotted data up to 2000 turns did not show any instability). Alternatively, Figs. 3(d)-(g) show that the unstable behavior smoothly vanishes with increasing and the evolution of the ray becomes stable for μm. Figure 3(h) and (i) plot the sine of incidence angle and the relative deviation for a single reflection point, which immediately follows completion of the 200th turn, as a function of the microfiber radius. The latter plots are useful for comparison with the wave calculations of Section 6. The size of the region of instability found from Fig. 3 is ~0.45 μm in reasonable agreement with Eq. (19), which gives for the region of parametric resonance μm.
6. Propagation of WGMs: mode coupling and the wave parametric resonance
Let us investigate the effect of the introduced periodic perturbation on the propagation of WGMs along the MFC. As previously, it is assumed that the coiled microfiber has the uniform radius , while the central rod is cylindrically asymmetric and possesses the n-fold rotational symmetry defined by Eq. (14). Similar to Ref [14], calculation is performed with the beam propagation model (BPM) where the local microfiber radius of curvature along the geodesic s is introduced by replacing of the refractive index with the effective index . In this approximation, the pitch p is set to zero, i.e., the torsion effect is ignored. Therefore, the essentially vector effects (polarization dependence, Berry phase [1–3], etc.) are neglected. The local radius of curvature in the latter equation is expressed through as . For and defined from Eq. (14) with , we have:
Close to the MFC parameters considered section 5, here we set and μm. In addition, we set μm, which, as follows from Eq. (20), approximately corresponds to μm of Section 5. The length of the coiled microfiber, , is set to 10 turns, which is equal to 100 grating periods or mm.The fundamental WGM is launched along the defined MFC. The transmission power of this mode, , is found by calculation of the squared overlap integral between the field distribution at position z along the microfiber and the initial distribution of the WGM. Generally, the behavior of depends on MFC parameters , , , and radiation wavelength λ.
Figure 4(a) plots the fundamental WGM power at the end of the MFC as a function of microfiber radius for the wavelength μm. It is seen that the transmission power has a resonant behavior in the region between 11 and 13 μm. Shifting of the resonant region from μm defined in Section 5 and from the unstable region in Fig. 3(h) and (i) towards larger microfiber radius values is due to the smaller effective radius of the WGM trajectory compared to .
The WGM propagation in the resonant region exhibits remarkable features. In particular, for certain values of , the fundamental WGM experiences stable periodic exchange of power with a single WGM, the phenomenon well known from the LPG theory [16]. An example of such behavior is shown in Fig. 5(a) for μm. This figure shows the transmission power of the fundamental WGM, , as a function of propagation distance z(upper plots) and the corresponding distribution of the radiation along the microfiber (lower plots). Inset surface plots are the cross-sectional radiation distribution at the positions indicated by arrows. For other values of , the wave parametric resonance distributes the fundamental WGM power among multiple neighboring WGMs. This situation is illustrated in Fig. 5(b) for μm. It is seen that now the behavior of transmission power is aperiodic and coupling between multiple modes accompanies the WGM evolution. The wavelength dependence on the transmission power P for μm and μm considered in Fig. 5(a) and (b) is shown in Fig. 4 (b). Notice the sharp peak in Fig. 4(b) cooresponding to λ = 0.82 μm and r 0 = 11.53 μm, which is attributed to the resonance transmission to the mode E 01 (see Eq. (8)).
The analysis and explanation of the obtained numerical results can be performed with the coupled wave equations. Let us introduce the adiabatic cross-sectional modes of the MFC, , which parametrically depends on the coordinate s, and represent the propagating radiation field as:
The coupled wave equations for the coefficients areThe asymptotic expressions for the WGMs following from Eqs. (1), (2) and (8) are . Here we are interested in resonant transitions between the series of modes with quantum numbers . For this series, the harmonic dependence of local radius of microfiber curvature defined by Eq. (20) yieldsand equidistant spacing between propagation constants, . Consider, for example, the situation shown in Fig. 5(a) when only two modes, and , are coupling. Substituting Eq. (23) into Eq. (22) and neglecting the fast oscillating terms yield the following coupling wave equations for these modes:From Eq. (24), the condition of phase matching is , which coincides with the classical parametric resonance condition described in Section 4. At this condition, modes and experience periodic full power transfer depicted in Fig. 5(a). The period of the power transfer determined from Eq. (24) isFor the parameters of the MFC considered in this Section, mm. The value of this period in Fig. 5(a), determined numerically, is 13.2 mm, which is ~25% less than the estimate given by Eq. (24). The reason of the discrepancy is in the asymptotic nature of Eqs. (1) and (2). In fact, for the considered MFC parameters, the size of mode along the vertical direction in the insets of Fig. 5(a) is comparable with the microfiber diameter, so that the parabolic approximation for the circular microfiber profile in the vicinity of geodesic s, which is assumed in Eqs. (1) and (2), is not accurate. However, Eq. (25) has a reasonable accuracy and is useful for estimates of the characteristic length along which the intermode power transfer takes place. Similarly, for the situation depicted in Fig. 5(b), the propagation constants of modes are not accurately equidistant as follows from Eq. (2). Nevertheless, the considered relatively small harmonic perturbation defined by Eq. (20) causes the multi-mode wave parametric resonance. Further analysis of the behavior of WGMs based on the coupled wave equations, Eq. (22), is beyond the scope of this paper.7. Discussion and summary
An ideal uniform MFC possesses an infinite helical symmetry. Deformations of an ideal MFC can be accidentally introduced in the process of its fabrication or specially introduced in order to arrive at certain transmission characteristics of this device. Application and further development of the results of this paper can be useful for the analysis of transmission losses caused by MFC deformations as well as for engineering of a specially deformed MFC.
The effect of periodic perturbations on the MFC transmission loss is of special interest. In fact, a conventional drawing station allows to fabricate 20-40 μm diameter silica microfibers, which are uniform enough to ensure the adiabatic suppression of the intermode transitions. In addition, the translation symmetry of the central rod can be maintained along the distances of a few tens of mm very precisely. However, the accuracy of rotation symmetry of the rod can be hardly made better than a few tenth of a percent. As the result, the light propagating along a coiled microfiber is subject to periodic bending perturbation caused by the axial asymmetry of the central rod which can be expressed by its local radius variation:
where the terms, which determine the axial asymmetry, are relatively small, . For small asymmetry of the central rod, the perturbation theory allows to separate the contribution of harmonics in Eq. (26) and to reduce the problem to the case of a single harmonic with period considered in this paper. In a more general case, both the microfiber and central rod random nonuniformities should be taken into account. These nonuniformities contain a periodic component, which can be expressed similar to Eq. (26).For the MFC engineering applications, the partial harmonic amplitude and phase parameters in Eq. (26), and , which determine the central rod shape, can be used to arrive at the desired transmission characteristics of the fundamental WGM propagating along the MFC. Here, again, for small , the design procedure is reduced to the case of a single harmonic investigated above.
In summary, the classical whispering gallery ray and WGM evolution in an MFC with n-fold helical symmetry is an interesting new example of the behavior of wave systems that are chaotic in the classical limit [25]. A simple condition of instability is found. For this condition, the classical stochastization, the classical and wave parametric resonance, and the stable WGM coupling, which results in periodic exchange of power, are investigated. Propagation of WGMs along an MFC with more complex nonuniformities can be considered similarly and applied to the evaluation of losses in MFC-based optical delay lines [14] and to the analysis of experimental data for fabricated MFCs [13,15].
References and links
1. M. Berry, “Interpreting the anholonomy of coiled light,” Nature 326(6110), 277–278 (1987). [CrossRef]
2. R. Y. Chiao and Y. S. Wu, “Manifestations of Berry’s topological phase for the photon,” Phys. Rev. Lett. 57(8), 933–936 (1986). [CrossRef] [PubMed]
3. G. B. Malykin and V. I. Pozdnyakova, “Geometric phases in single mode fiber light guides and fiber ring interferometers,” Phys-Usp. 47(3), 289–308 (2004). [CrossRef]
4. M. Sumetsky, “Optical fiber microcoil resonators,” Opt. Express 12(10), 2303–2316 (2004), http://www.opticsinfobase.org/abstract.cfm?URI=OPEX-12-10-2303. [CrossRef] [PubMed]
5. D. J. Hill, “The evolution and exploitation of the fibre-optic hydrophone,” Proc. SPIE 6619, 661907 (2007). [CrossRef]
6. C. N. Alexeyev, B. A. Lapin, and M. A. Yavorsky, “Optical vortices and topological phase in strongly anisotropic coiled few-mode optical fibers,” J. Opt. Soc. Am. B 24(10), 2666–2675 (2007). [CrossRef]
7. V. Atanasov and R. Dandoloff, “Curvature-induced quantum behaviour on a helical nanotube,” Phys. Lett. A 372(40), 6141–6144 (2008). [CrossRef]
8. M. Sumetsky, “Basic elements for microfiber photonics: micro/nanofibers and microfiber coil resonators,” J. Lightwave Technol. 26(1), 21–27 (2008). [CrossRef]
9. N. G. Broderick, “Optical snakes and ladders: dispersion and nonlinearity in microcoil resonators,” Opt. Express 16(20), 16247–16254 (2008), http://www.opticsinfobase.org/abstract.cfm?&uri=oe-16-20-16247. [CrossRef] [PubMed]
10. J. Scheuer, “Fiber microcoil optical gyroscope,” Opt. Lett. 34(11), 1630–1632 (2009). [CrossRef] [PubMed]
11. M. Sumetsky, “Optical microfiber loop and coil resonators,” in Practical applications of microresonators in optics and photonics, Ed. A. Matsko (CRC Press, 2009) p. 355–384.
12. G. Brambilla, F. Xu, P. Horak, Y. Jung, F. Koizumi, N. P. Sessions, E. Koukharenko, X. Feng, G. S. Murugan, J. S. Wilkinson, and D. J. Richardson, “Optical fiber nanowires and microwires: fabrication and applications,” Adv. Opt. Photon. 1(1), 107–161 (2009). [CrossRef]
13. L. Yao, T. A. Birks, and J. C. Knight, “Low bend loss in tightly-bent fibers through adiabatic bend transitions,” Opt. Express 17(4), 2962–2967 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-17-4-2962. [CrossRef] [PubMed]
14. M. Sumetsky, “Optical microfiber coil delay line,” Opt. Express 17(9), 7196–7205 (2009), http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-17-9-7196. [CrossRef] [PubMed]
15. M. Sumetsky, Y. Dulashko, and S. Ghalmi, “Fabrication of miniature optical fiber and microfiber coils,” Opt. Lasers Eng. 48 , 272-275 (2010) [CrossRef]
16. T. Erdogan, “Fiber grating spectra,” J. Lightwave Technol. 15(8), 1277–1294 (1997). [CrossRef]
17. S. Ramachandran, “Dispersion-tailored few-mode fibers: a versatile platform for in-fiber photonic devices,” J. Lightwave Technol. 23(11), 3426–3443 (2005). [CrossRef]
18. H. Xuan, W. Jin, and M. Zhang, “CO2 laser induced long period gratings in optical microfibers,” Opt. Express 17(24), 21882–21890 (2009), http://www.opticsinfobase.org/abstract.cfm?URI=oe-17-24-21882. [CrossRef] [PubMed]
19. V. M. Babich, and V. S. Buldyrev, Short-Wavelength Diffraction Theory. Asymptotic Methods. (Springer, Berlin, 1991).
20. M. Berry and K. Mount, “Semiclassical approximation in wave mechanics,” Rep. Prog. Phys. 35(1), 315–397 (1972). [CrossRef]
21. The Gaussian curvature of a surface is the inversed product of its principal radii of curvature [22].
22. See, e.g., M. P. Do Carmo, Differential geometry of curves and surfaces (Prentice Hall, Inc., Englewood Cliffs, New Jersey, 1976).
23. L. D. Landau, and E. M. Lifshitz, Mechanics (Pergamon Press, 1976).
24. A. M. Perelomov, and Y. B. Zeldovich, Quantum mechanics: selected topics (World Scientific Publishing Co, Singapore, 1998).
25. Quantum chaos,http://en.wikipedia.org/wiki/Quantum_chaos.