Dense photonic integration promises to revolutionize optical computing and communications. However, efforts towards this goal face unacceptable attenuation of light caused by surface roughness in microscopic devices. Here we address this problem by introducing Surface Nanoscale Axial Photonics (SNAP). The SNAP platform is based on whispering gallery modes circulating around the optical fiber surface and undergoing slow axial propagation readily described by the one-dimensional Schrödinger equation. These modes can be steered with dramatically small nanoscale variation of the fiber radius, which is quite simple to introduce in practice. Extremely low loss of SNAP devices is achieved due to the low surface roughness inherent in a drawn fiber surface. In excellent agreement with the developed theory, we experimentally demonstrate localization of light in quantum wells, halting light by a point source, tunneling through potential barriers, dark states, etc. This demonstration has intriguing potential applications in filtering, switching, slowing light, and sensing.
© 2011 OSA
For the past few decades, research and development in photonic integrated circuits [1–9] has struggled to build a platform with the miniature dimensions, flexibility, and control needed to deliver breakthrough capabilities in optical computing, communications, and fundamental science. The platform must include waveguides for propagation and routing light, but must also provide means for trapping light. Building microscopic high-quality-factor traps for photons is non-trivial. One cannot simply construct a microscopic version of a traditional optical resonator, since naturally reflecting materials (metals) are highly absorbing. Also, one cannot construct photonic quantum wells using a straightforward analogy to electronics and quantum mechanics. Indeed, this would require total internal reflection of light everywhere at the boundary surrounding the trap, a challenge for the low-loss optical cavities with microscopic dimensions .
Two basic platforms for fabrication of photonic circuits with microscale elements have been developed. The first, a ring resonator platform, is built on thin transparent photonic wires that tightly confine light due to the high refractive index contrast between the wire and surrounding material. In this platform, wires form closed rings that are coupled together to provide higher level of functionality [1,2,8,9]. The second platform is based on photonic crystals, materials with periodic modulation of refractive index . Thin waveguides are created using a linear disruption of periodicity while microresonators are formed by localized defects [3–7]. With the advances in lithographic technology and design methods, these two platforms have steadily increased in complexity, achieved lowed loss, and broadened functional capability. The miniature integrated photonic circuits that have been proposed and demonstrated using the ring-resonator and photonic crystals platforms can potentially revolutionize optical communications and computing by implementation of superfast all-optical processors of light pulses on a chip and, in particular, miniature optical buffers [1,3] and all-optical switches [2,6].
Yet, despite the remarkable accomplishments, the existing platforms still face a severe limitation: It is necessary but very difficult to fabricate photonic circuits that have both microscopic dimensions and ultra-small losses (required, e.g., for creating the ultra-high quality-factor microresonators). Indeed, the microscopic dimensions of photonic elements can only be achieved with high-index-contrast material interfaces. However, precisely in this regime of high index contrast, the sensitivity to surface roughness and fabrication errors is daunting and potentially fundamentally limiting. To reduce the sensitivity to the interface roughness and uncontrolled attenuation, photonic circuits with low-index contrast can be exploited [9,10]. However, to avoid radiation losses at low-contrast waveguide bends, the size of photonic resonators has to have the millimeter scale [9,10]. These photonic circuits might be low-loss and easier to fabricate, but they are no longer microscopic.
In this paper, we present a new photonic platform with ultra-low loss, flexibility, and elements having microscopic dimensions . Our solution consists in Surface Nanoscale Axial Photonics (SNAP) employing whispering gallery modes (WGMs) [12–15] that circulate circumferentially around the surface of a thin optical fiber while also undergoing slow propagation along the fiber axis (Fig. 1 ). SNAP elegantly manifests a number of key properties:
- (a) Remarkably low surface roughness effortlessly results from the fiber draw process. The key benefit of SNAP is the orders of magnitude lower loss (~10−5 dB/cm [14,15]) achieved at a high-index-contrast silica surface compared to high-index-contrast lithography-based technologies (~0.1 dB/cm [1–4]).
- (b) Periodicity and slow light without the periodic modulation of the refractive index (as for photonic crystals): Periodicity is introduced automatically by each revolution around the fiber surface. Axial propagation naturally has features of slow light, since the propagation of slow WGMs is primarily azimuthal.
- (c) SNAP is readily described by the one dimensional Schrödinger equation: The axial propagation of slow WGMs exhibits turning points, localization in quantum wells, tunneling through barriers, etc.
- (d) Precise control of light achieved with nanoscale variation of the effective fiber radius, which define the “potential” that steers, localizes, and engineers WGMs.
- (e) Novel phenomenon of surface WGMs completely halted with a point source of light.
- (f) Microscopic length scale in each dimension: For slow WGMs, the characteristic axial wavelength is much larger than the wavelength of light and has the order of a few tens of microns. This dramatically simplifies fabrication of SNAP devices and, in particular, the ring microresonators. As opposed to the macroscopic increase of the diameter of ring resonators in low-contrast photonic circuits, we increase the third dimension along the fiber axis and keep it microscopic.
2. SNAP theory: prototyping one-dimensional quantum mechanics with light
Modes in an optical fiber are characterized by the propagation constant which depends on the radiation wavelength λ and variations of the fiber radius, and refractive index . In conventional optical fibers, light is directed along the interior fiber core and has the propagation constant close to . In contrast, SNAP employs transverse WGMs wrapped around the fiber surface by total internal reflection. The propagation constant of these modes is much smaller than and the speed of their axial propagation is much smaller than the speed of light in the fiber material, . In fact, the axial speed of a WGM and its propagation constant can be zero at the resonance wavelength defined by the condition of stopped axial propagation at a point z0, where the resonance width determines the propagation loss.
The central idea of SNAP is to exploit the sensitivity of WGMs to the extremely small variations of the fiber radius and index near the resonance . Generally, a variation in radius causes coupling between modes and intermodal transitions, a complex problem which is usually addressed with the system of coupled wave equations . In SNAP, this problem is absent: Variation and is so small and smooth that the coupled wave equations become decoupled and the dependence of the WGM field on the axial coordinate, A(z), is defined by a single differential equation. Near the resonance wavelength , this equation takes the form of the stationary Schrödinger equation (see Appendix A, Subsection A.1):Eq. (1) is proportional to the effective radius variation,
Near the resonance wavelength , dramatically small effective radius variations are enough to reflect and confine light. Indeed, it follows from Eq. (1) that for the typical experimentally determined value of μm at wavelength μm and fiber radius μm, light can be confined with an effective radius variation of the order of pm. This corresponds to feature perturbations of pm and The characteristic axial wavelength in this case is 1mm, i.e., 1000 times greater than the wavelength of light. In the experiments described below in Sections 3 and 4, the characteristic axial wavelength is ~50 μm.
Slow WGMs can be excited in an optical fiber using a microfiber (MF) , specifically, a micrometer diameter waist of a biconical fiber taper, which is attached normal to the SNAP device and connected to the light source and detector as illustrated in Fig. 1. The MF waveguide acts as a point source in Eq. (1) so that the amplitude of the WGM excited by the MF positioned at is expressed through the Green’s function of Eq. (1),
The spectrum of excited WGMs can be observed by measuring the transmission amplitude of the MF,
Parameter C in Eqs. (3) and (4) determines coupling to the MF. It is a weak function of wavelength and, with good accuracy, is constant in the vicinity of resonance . The derivation of Eqs. (3) and (4) based on the formalism of the Lippmann-Schwinger equation [18,19] is given in Subsection A.2 of Appendix A.
SNAP devices include the three basic building blocks illustrated in Fig. 2 . The first is the WGM bottle microresonator  coupled to a MF (Fig. 2(a)). This corresponds to a quantum well in Eq. (1), and leads to the formation of states localized between turning points and at a discrete series of wavelengths corresponding to resonances in the transmission amplitude . From Eq. (4), the resonance transmission amplitude is proportional to the WGM amplitude at the MF position . In particular, if coincides with a node of the WGM then the coupling to the MF vanishes and the WGM becomes dark, i.e., it practically does not show up in (for details, see Appendix A, Subsection A.3.3).
A second building block is a concave fiber waist (Fig. 2(b)). In this case, for wavelengths in the underbarrier region, , the WGM is localized due to the exponential decay of its amplitude away from the MF position . Alternatively, above the barrier, for , the excited WGM is delocalized. In the particular case of a uniform fiber, , Fig. 2(b) offers a simple explanation of the localization of light in a uniform cylindrical microresonator described with a less general semiclassical theory in Ref. . (see Appendix A, Subsection A.3.1).
In the third building block (Fig. 2(c)), the MF is positioned near a turning point where the SNAP device has monotonically increasing radius. In this case, the wave that is launched by the MF along the positive axial direction, interferes with the wave that is launched along the negative direction and reflects from the turning point . At discrete wavelengths when the condition of destructive interference of these two waves is fulfilled, the distribution of light is fully localized between the turning point and MF position . This effect is similar to the situation in the Mach-Zehnder interferometer when light propagating along two arms of the interferometer approaches the region of coupling with opposite phases and equal amplitudes so that the waves in two arms cancel each other. Figure 2(c) offers a simple explanation of the remarkable effect when a point contact with a MF source completely halts light propagating along the SNAP device, while the detailed theory of this effect is given in Subsection A.3.2 of Appendix A. In particular, Fig. 2(c) clarifies the appearance of localized conical modes discovered in Ref. .
These three basic building block elements can be used to create a wide range of densely integrated photonic devices by coupling them together. They can also be coupled to one or more transverse MFs to provide multiple inputs and outputs. In such integrated assemblies, the field distribution along the SNAP device as well as the transmission spectrum, group delay, and dispersion of light transmitted through MFs can be engineered using Eq. (1), (3), and (4). Recently we have shown that it is possible to reproducibly fabricate SNAP microresonators at the surface of an optical fiber with sub-angstrom accuracy . Further development and application of methods reported in  potentially allows creation of low-loss, densely packed devices for filtering, switching, slowing light, and sensing.
In the experiments below we consider SNAP devices which reproduce the nanoscale radius variation of Fig. 2. In doing this, we experimentally demonstrate all the described basic elements and find excellent agreement of their performance with the developed theory.
3. Experiment: Localization in quantum wells, tunneling, dark states
The described basic SNAP phenomena, familiar from elementary quantum mechanics , are confirmed experimentally in excellent agreement with the developed theory. In our experiments, to arrive at small dimensions of SNAP devices, we fabricated them from a regular optical fiber drawn down to The nanoscale variation of the effective fiber radius, was introduced using a simple technique of controlled local heating with CO2 laser and pulling. Thus we avoid the surface roughness that might result from other methods of fiber post-processing. It is known that post-processing of silica by annealing and melting causes variations of both the density and refractive index [25, 26]. In addition, post-processing of an optical fiber result in its refractive index and physical radius variation due to the relaxation of the frozen-in tension, which is introduced by the pulling force in the process of fiber drawing [27, 28]. From Eqs. (1) and (2), only the combined variation of the refractive index and physical fiber radius entering the expression for the effective fiber radius affects the performance of SNAP devices. Due to the drastic elongation along the fiber axis, the field distribution in SNAP devices is much easier to access, control, and engineer as compared to the WGM-based microdevices demonstrated previously including spherical, toroidal, and ring microresonators [14, 15] as well as microresonators directly fabricated by relatively small deformation of the optical fiber surface .
We fabricated samples of SNAP devices reproducing the characteristic nanoscale fiber radius variation illustrated in Fig. 2, the shape of an elongated bottle (quantum well) with a neck (barrier). These samples include bottle microresonators with multiple axial localized states (Figs. 3(b), (c) ), three axial localized states (Figs. 3(d), (e)) and a single axial localized state (Fig. 3(f), (g)). The samples were experimentally characterized as follows. First, an MF was translated along the test fiber in 20 μm steps where the transmission amplitude spectra (vertical plots in Figs. 3(b), (d), and (f)) were measured with an optical vector analyzer (1.3 pm wavelength resolution). In our experiments, the resonances corresponding to different polarizations of light did not overlap in transmission spectra and only one polarization was considered. The effective fiber radius variations (bold curves in Figs. 3(c), (e), and (g)) are found from these plots by applying the MF scanning method [30, 31] corrected with Eqs. (1) and (4) (see Appendix A, Subsection A.5).
In excellent agreement with theory, we observed full localization of light in dramatically shallow bottle microresonators. In Figs. 3(b), (c), the fiber shape features an elongated (300 μm in axial length) and extremely shallow (only 7 nm in radius variation) bottle microcavity. In agreement with the developed theory, the bottle resonances in Fig. 3(b) are localized between turning points zt1 and zt2. Depending on the value of coupling with the MF, the Q-factor of these resonances varies from a relatively low () or greater than 106 (with the resolution limited by the measurement device). The resonance state can be dark (i.e., practically uncoupled from the MF) if the MF position approaches a node of the localized state or is deep in the underbarrier region (see Appendix A, Subsection A.3.3). Approximating the effective radius variation in Figs. 3(b) and (d), with the quadratic dependence, , we find the axial radius-of-curvature m of the multi-level bottle microresonator (Fig. 3(c)) and a smaller m for the three-level bottle microresonator (Fig. 3(e)). The experimental free spectral range of the resonances in Figs. 3(b) and (d) is 0.082 nm and 0.38 nm, which is in good agreement with the theoretical values 0.08 nm and 0.39 nm found from Eq. (1) for the harmonic oscillator, (see Appendix A, Subsection A.4). At the field maxima, the resonance width is primarily determined by coupling to the MF (see Appendix A, Subsection A.3.3). For example, in Fig. 3(b) the maximum resonance width nm. Then, using Eq. (A.3.15) the value of coupling to MF is estimated, . Finally, Figs. 3(f) and (g) demonstrate a shallow bottle resonator with a single axial state. The length of this resonator is 70 μm and its height is only 5 angstrom. It is remarkable that such a minor deviation from uniformity is able to fully confine light.
Oscillations of spectra outside the quantum wells in Fig. 3 are easily explained. These oscillations result from the interference between light propagating directly through the MF and light, which couples into the test fiber, reflects from a turning point, and then couples back into the MF. Interestingly, on the right hand side of the barriers in Figs. 3(b), (d), and (f) there are high Q-factor resonances, which coincide with the quantum well resonances but are situated outside of the quantum wells. Appearance of these resonances is explained with Fig. 2: Light excited by a MF in the region of Fig. 2(c) propagates in direction of the quantum well, tunnels through the barrier, and resonantly attenuates. In the vicinity of the tip of the barrier, the tunneling rate may affect the Q-factor of resonances. Finally, if the MF is positioned deeply in the underbarrier region (bottom of Figs. 3(b), (d), and (f)), the outgoing WGMs strongly decay, interference is absent, and the transmission spectra are smooth.
The characteristic axial length of the low-lying WGMs of the SNAP bottle microresonators, μm, is surprisingly small for so large axial radii-of-curvature R. This fact is simply explained by the expression for the WGM localized in a quadratic quantum well, Eq. (A.4.6) of Appendix A, yielding . For example, the axial size Lz increases only by a factor of ten when radius R grows from the typical μm [14,15] to the gigantic m considered here. For the same reason, dramatic increase of the axial radius does not significantly affect the value of the WGM volume, an important parameter in cavity quantum electrodynamics [7,14]. The volume of the mode is usually defined as the ratio of the integral of the field power over the cavity volume and the maximum field power. From Eq. (A.4.6), we find where for the first axial states the dimensionless parameter , and .
4. Experiment: Halting light with a point source
We have experimentally demonstrated the remarkable effect of halting light with a point source, as predicted from our theory (see Fig. 2(c) and the description above). This effect takes place takes for the SNAP device illustrated in Fig. 4(a) when light excited at reflects from and destructively interferes at . The transmission spectra and radius variation of the device, measured as in the previous case, are shown in Figs. 4(b) and (c). The localization of light is confirmed experimentally with two microfibers, MF1 and MF2. MF2 is translated along the test fiber and probes the field excited by MF1 as illustrated in Fig. 4(a). The measurement results are shown in Fig. 4(d), where the vertical axis lines of all spectra correspond to zero transmission. It is seen that at the discrete values of wavelengths indicated by horizontal arrows, the field distribution is fully localized along the fiber segments with a length of less than 100 μm and a radius variation of less than 5 nm. To compare the measurement results with theory, the fiber radius variation is approximated by the quadratic dependence , and the Green’s function of Eq. (1) is found analytically  (see Appendix A, Subsection A.4). The axial fiber radius, m, is found from Fig. 4(c). With this value, the comparison of the transmission spectra measured experimentally and calculated with Eq. (4) and Eq. (A.4.2) of Appendix A (respectively, black and blue curves in Fig. 4(b)) shows excellent agreement. A minor deviation near the principal peak maxima is explained by the deviation of the actual radius variation from the quadratic dependence away from its minimum. In addition, in our calculations, we ignored the effect of the probe MF2, which can be made relatively small by decreasing the coupling to MF2 (experimentally, this can be achieved by translating MF2 along its length and increasing the MF2 diameter at the contact point). Finally, for comparison, the surface plot of the theoretical calculation of the WGM field amplitude distribution found from Eq. (A.4.2) for the parameters of the experiment, Fig. 4(f), is placed in the background of Fig. 4(d). Good agreement is found both for the positions of the localized states and for the field distribution. Note that coupling to the MF2 vanishes together with the field amplitude and, hence, does not change the condition of full localization. A deviation near the principal peaks is, again, due to the assumed quadratic approximation for the radius variation and coupling to the MF2.
Nanoscale variation of the effective optical fiber radius (including the variation of the fiber radius and/or refractive index) enables formation and integration of independent or coupled microdevices (e.g., a sequence of microresonators) on a SNAP platform. The axial radiation wavelength of SNAP microdevices is significantly larger than the wavelength of light, which simplifies their operation and broadens the field of potential applications. The further decrease of the axial size of SNAP elements can be achieved with larger variation of the effective fiber radius while retaining adiabatic behavior and adherence to the developed theory. Due to the very low attenuation of light propagation along the optical fiber surface, the SNAP microdevices can exhibit significantly improved performance in filtering, time delay, slowing light, switching, sensing, etc., compared to lithographically fabricated photonic circuits.
A.1. Derivation of Eq. (3) for the propagation constant of a slow WGM
The propagation constant of optical modes in a fiber with radius and refractive index , which is situated in the surrounding medium with index , is defined by the equation 
where and are the Bessel and Hankel function, respectively. In the zero approximation, the propagation constant is set to zero, , and Eq. (A.1.1) is split into two equations for TE and TM modes:
where is a root of the Airy function (, , ), , and signs + and – correspond to TE and TM polarizations. From Eq. (A.1.7), the free spectral range (FSR) is
Assume now that the wavelength , the fiber radius r, and the fiber refractive index are slightly shifted from and , respectively. For a WGM propagating close to the fiber surface, in the first approximation in and , iteration of Eq. (A.1.1) near the resonance yields
This equation leads to the expression for in Eq. (1).
A.2. Derivation of Eq. (3) for the WGM field distribution and Eq. (4) for the resonant transmission amplitude
Generally, the field excited by a MF in the SNAP device is localized in the vicinity of the MF/SNAP fiber contact point and does not have the rotational symmetry. However, near the resonance, , the beam launched by the MF constructively interferes in the process of circulation along the SNAP fiber surface. Then, after a large number of turns, the beam acquires axial symmetry. The axially symmetric component of the beam becomes much greater in amplitude compared to its original asymmetric part, so that the latter can be ignored.
For the nanoscale and adiabatic variation of the SNAP fiber radius and outside the region of coupling with the MF, i.e., in the absence of source, Eq. (1) is the known uncoupled wave equation  with the propagation constant defined by Eq. (A.1.9). For a MF with the diameter ~1 μm, the axial size of the coupling region is μm. Near the resonance , is much smaller than the axial wavelength . Indeed, from Eq. (A.1.9), for μm, μm, nm, and nm, we have μm. Thus, the effect of the MF positioned at can be described by adding a δ-function point source to the right hand side of Eq. (1). This immediately yields the expression for WGM amplitude defined by Eq. (3). More accurately, Eq. (3) for the WGM distribution and Eq. (4) for the transmission amplitude can be derived based on the formalism of the Lippmann-Schwinger equation. Following Ref. , the Hamiltonian describing the electromagnetic waves in the optical fiber coupled to a MF waveguide is approximated by , where is the Hamiltonian of the uncoupled states in the MF and SNAP fiber and defines the MF/SNAP fiber coupling. Here is an uncoupled waveguide mode of the single mode MF and is an uncoupled WGM. For the axially symmetric fiber with very small radius variation, in the adiabatic approximation, the WGM factors in the cylindrical coordinates into a product of axial , radial, , and azimuthal, , components:
m is the discrete azimuthal quantum number, p is the discrete radial quantum number, q is the discrete or continuous axial quantum number, and the amplitude satisfied Eq. (1). All the factors in Eq. (A.2.1) are assumed to be normalized (e.g., ). The coupling operator V is spatially localized in a small vicinity of the MF contact point with the coordinate . As mentioned, the axial length of this region is much smaller than the characteristic variation length of so that the coupling elements of this operator, , are recast as
The total wave function is determined by the Lippmann-Schwinger equation,
where ε is a positive infinitesimal number enforcing the outgoing boundary condition and T is the scattering T-matrix. After approximate diagonalization and renormalization, the elements of the T-matrix, which determine the transition from the incoming wave into the WGM are [18, 19]:
where is the full resonant width of the state , which takes into account both transmission losses in the SNAP fiber and leakage in the MF/SNAP fiber coupling region. The part of localized in the SNAP fiber which is excited by the income wave is
where is the Green’s function of Eq. (1) at . In the vicinity of the resonance all the terms in the second line of this equation except for those with can be ignored so that
The transmission amplitude through the MF is found as 
and in the neighborhood of a resonance, , to
which coincides with Eq. (4) for and .
A.3. The WGM distribution and transmission amplitude in the semiclassical approximation
A.3.1. No turning points. Localization in a uniform SNAP device
For a uniform SNAP fiber, , the propagation constant is independent of z,
Substitution of this equation into Eq. (A.3.1) shows that, as expected, the WGM field exponentially vanish at both sides of the MF due to the presence of loss defined by . In the absence of attenuation, , in accordance with illustration in Fig. 2(b), the WGM amplitude is delocalized and is uniform along the axis z for , i.e., above the potential barrier. It is localized and exponentially decaying for , i.e., in the underbarrier region. From Eqs. (A.3.1), (A.3.2), and (A.3.3) the following expressions for the field amplitude and transmission amplitude are found:
The tree-dimensional problem of resonant transmission through a uniform SNAP device coupled to a MF was solved in Ref.  by calculation of the sum over the turns of circulating beam excited in the SNAP fiber. Equation (A.3.2) coincides with the expression for the WGM field amplitude obtained in Ref. . However, Eq. (A.3.4) corrects the expression for the transmission amplitude of Ref. , which mistakenly contained the factor i in front of . With this correction, the shape of the transmission resonance appears to be asymmetric rather than symmetric. In particular, for the small MF/SNAP fiber coupling C, the transmission power is
A.3.2. One turning point. Localization enforced by a point contact
If the radius variation of an optical fiber is monotonic, the presence of a MF may still lead to full localization of a WGM as illustrated in Fig. 2(c). In this case, solution of Eq. (1) exponentially vanishes to the left hand side of the turning point, , and is an outgoing wave to the right hand side of the MF position, . The Green’s function of Eq. (1) found in the semiclassical approximation with these boundary conditions is
From this equation, light is confined between the turning point and the MF position if the following condition of destructive interference at is fulfilled:
where n is a large positive integer. Then, in the considered approximation, the transmission amplitude defined by Eq. (4) is equal to unity, i.e., light does not couple into the SNAP fiber and the localized state becomes dark.
A.3.3. Two turning points. A bottle microresonator. Dark states
If the MF is located inside the bottle microresonator, which contains multiple localized states as illustrated in Fig. 2(a), and is separated from the turning points by the distance greater than the axial wavelength , then Eq. (1) can be solved in the semiclassical approximation. In particular, the transmission amplitude, Eq. (4), is
with defined by Eq. (A.3.8). The resonant peaks of the transmission amplitude correspond to the zeros of denominator in this equation , which are determined by the quantization rule
Expansion of Eq. (A.3.10) near yields
If the numerator in Eq. (A.3.12) vanishes then the bottle microresonator becomes a dark state uncoupled from the MF.
It is instructive to compare the expression for the resonant transmission amplitude through a localized state, Eq. (A.3.12), with the known formula
As expected, the coupling coefficient achieves maxima at the field maxima when cos2(…) in Eq. (A.3.14) equals unity. For example, for the parabolic quantum well, when the radius dependence is approximated by we have
A.4. Exact solution for the quadratic variation of the effective fiber radius
The semiclassical approximation of Section 3 fails near turning points and also near the bottom of quantum wells and the top of potential barriers. Solution of Eq. (1) in these regions is critical because it corresponds to the edge peaks of the resonant transmission amplitude of the bottle microresonators (Fig. 3) and also to largest principal peaks featuring slopes and neck regions of the fiber (Fig. 4). The effective radius variation defined by Eq. (2) in these cases can be approximated by the quadratic dependence:
where the axial radius R can be positive for the neck-shaped SNAP fiber and negative for the bottle-shaped fiber. Then, the Green function of Eq. (1) can be expressed through the parabolic cylinder functions. For numerical simulations, it is more convenient to use the integral representation :
are the dimensionless axial coordinate and wavelength, and
are the characteristic lengths along and .
For negative R, the WGMs of the SNAP bottle microresonator are found as
where Hq(x) is the Hermite polynomial, Ai(x) is the Airy function and tp is the root of Airy function (see, e.g., ).
A.5. Experimental measurement of the effective radius variation
In the semiclassical approximation, the relation between the resonance amplitude and effective radius variation of a SNAP fiber can be clarified with Fig. 2. If the MF is situated in the region where the SNAP fiber has a bottle shape (Fig. 2(a)), the spectra possess sharp resonances corresponding to localized WGMs (see, e.g., Fig. 3). For each MF position , the resonant peaks of the amplitude are located in the interval between the turning point wavelength, and the wavelength corresponding to the top of the potential barrier and independent of . Thus, the effective radius variation is found as
In the case of extremely slow variation of the fiber radius, reflection of WGMs from the turning points can be ignored. The peak of this amplitude is shifting proportionally to so that
Equation (A.5.2) resembles the approach of Refs [30, 31]. Finally, for a SNAP fiber segment, which does not contain localized bottle states, (e.g., the regions near the fiber neck (Fig. 2(b)) and to the right hand side of the neck (Fig. 2(c)), the resonance amplitude has the Airy-type oscillations (see, e.g., Fig. 4 (b)) It can be shown numerically, that, in this case, the fiber radius variation is proportional to the shift of the largest main peak of the resonance amplitude, :
The authors are grateful to Y. Dulashko for assisting in the experiments, to D. J. DiGiovanni for useful discussions and suggestions, and to M. Fishteyn for useful consultations.
References and links
1. F. N. Xia, L. Sekaric, and Y. Vlasov, “Ultracompact optical buffers on a silicon chip,” Nat. Photonics 1(1), 65–71 (2007). [CrossRef]
2. Y. Vlasov, W. M. J. Green, and F. Xia, “High-throughput silicon nanophotonic wavelength-insensitive switch for on-chip optical networks,” Nat. Photonics 2(4), 242–246 (2008). [CrossRef]
3. M. Notomi, E. Kuramochi, and T. Tanabe, “Large-scale arrays of ultrahigh-Q coupled nanocavities,” Nat. Photonics 2(12), 741–747 (2008). [CrossRef]
4. T. Baba, “Slow light in photonic crystals,” Nat. Photonics 2(8), 465–473 (2008). [CrossRef]
5. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the Flow of Light (Princeton University Press, 2008).
6. M. Notomi, “Manipulating light with strongly modulated photonic crystals,” Rep. Prog. Phys. 73(9), 096501 (2010). [CrossRef]
7. U. K. Khankhoje, S.-H. Kim, B. C. Richards, J. Hendrickson, J. Sweet, J. D. Olitzky, G. Khitrova, H. M. Gibbs, and A. Scherer, “Modelling and fabrication of GaAs photonic-crystal cavities for cavity quantum electrodynamics,” Nanotechnology 21(6), 065202 (2010). [CrossRef] [PubMed]
8. W. Bogaerts, S. K. Selvaraja, P. Dumon, J. Brouckaert, K. De Vos, D. Van Thourhout, and R. Baets, “Silicon-on-insulator spectral filters fabricated with CMOS technology,” IEEE J. Sel. Top. Quantum Electron. 16(1), 33–44 (2010). [CrossRef]
9. A. Melloni, A. Canciamilla, C. Ferrari, F. Morichetti, L. O'Faolain, T. F. Krauss, R. De La Rue, A. Samarelli, and M. Sorel, “Tunable delay lines in silicon photonics: coupled resonators and photonic crystals, a comparison,” IEEE Photonics J. 2(2), 181–194 (2010). [CrossRef]
10. C. R. Doerr and K. Okamoto, “Advances in silica planar lightwave circuits,” J. Lightwave Technol. 24(12), 4763–4789 (2006). [CrossRef]
11. M. Sumetsky, “Localization of light in an optical fiber with nanoscale radius variation,” in CLEO/Europe and EQEC 2011 Conference Digest, OSA Technical Digest (CD) (Optical Society of America, 2011), Postdeadline Paper PDA_8.
12. L. Rayleigh, “The problem of the whispering gallery,” Philos. Mag. 20, 1001–1004 (1910).
13. A. N. Oraevsky, “Whispering-gallery waves,” Quantum Electron. 32(5), 377–400 (2002). [CrossRef]
15. A. B. Matsko and V. S. Ilchenko, “Optical resonators with whispering-gallery modes—part I: Basics,” IEEE J. Sel. Top. Quantum Electron. 12(1), 3–14 (2006). [CrossRef]
16. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman and Hall, London, 1983).
18. S. Fan, P. R. Villeneuve, J. D. Joannopoulos, M. Khan, C. Manolatou, and H. Haus, “Theoretical analysis of channel drop tunneling processes,” Phys. Rev. B 59(24), 15882–15892 (1999). [CrossRef]
19. Y. Xu, Y. Li, R. K. Lee, and A. Yariv, “Scattering-theory analysis of waveguide-resonator coupling,” Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics 62(55 Pt B), 7389–7404 (2000). [CrossRef] [PubMed]
23. M. Sumetsky, D. J. DiGiovanni, Y. Dulashko, J. M. Fini, X. Liu, E. M. Monberg, and T. F. Taunay, “Surface nanoscale axial photonics: robust fabrication of high quality factor microresonators,” Opt. Lett. (to be published).
24. L. D. Landau and E. M. Lifshitz, Quantum Mechanics, (Pergamon Press, 1977).
25. A. Q. Tool, L. W. Tilton, and J. B. Saunders, “Changes caused in the refractivity and density of glass by annealing,” J Res Natl Bur Stand (1934) 38(5), 519–526 (1947). [PubMed]
26. H. Bach and N. Neuroth eds, The Properties of Optical Glass (Springer Verlag, 1995).
27. A. D. Yablon, M. F. Yan, P. Wisk, F. V. DiMarcello, J. W. Fleming, W. A. Reed, E. M. Monberg, D. J. DiGiovanni, J. Jasapara, and M. E. Lines, “Refractive index perturbations in optical fibers resulting from frozen-in viscoelasticity,” Appl. Phys. Lett. 84(1), 19–21 (2004). [CrossRef]
28. J. W. Fleming, “Sub glass transition relaxation in optical fibers,” in Optical Fiber Communication Conference, Technical Digest (CD) (Optical Society of America, 2004), paper TuB2.
29. A. A. Savchenkov, I. S. Grudinin, A. B. Matsko, D. Strekalov, M. Mohageg, V. S. Ilchenko, and L. Maleki, “Morphology-dependent photonic circuit elements,” Opt. Lett. 31(9), 1313–1315 (2006). [CrossRef] [PubMed]
30. T. A. Birks, J. C. Knight, and T. E. Dimmick, “High-resolution measurement of the fiber diameter variations using whispering gallery modes and no optical alignment,” IEEE Photon. Technol. Lett. 12(2), 182–183 (2000). [CrossRef]
32. R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, New York, 1965).